Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Chemokine CCL19 promotes type 2 T-cell differentiation and allergic airway inflammation

BACKGROUND

Allergic asthma is driven largely by allergen-specific TH2 cells, which develop in regional lymph nodes on the interaction of naive CD4+ T cells with allergen-bearing dendritic cells that migrate from the lung. This migration event is dependent on CCR7 and its chemokine ligand, CCL21. However, is has been unclear whether the other CCR7 ligand, CCL19, has a role in allergic airway disease.

OBJECTIVE

This study sought to define the role of CCL19 in TH2 differentiation and allergic airway disease.

METHODS

Ccl19-deficient mice were studied in an animal model of allergic asthma. Dendritic cells or fibroblastic reticular cells from wild-type and Ccl19-deficient mice were cultured with naive CD4+ T cells, and cytokine production was measured by ELISA. Recombinant CCL19 was added to CD4+ T-cell cultures, and gene expression was assessed by RNA-sequencing and quantitative PCR. Transcription factor activation was assessed by flow cytometry.

RESULTS

Lungs of Ccl19-deficient mice had less allergic airway inflammation, reduced airway hyperresponsiveness, and less IL-4 and IL-13 production compared with lungs of Ccl19-sufficient animals. Naive CD4+ T cells cocultured with Ccl19-deficient dendritic cells or fibroblastic reticular cells produced lower amounts of type 2 cytokines than did T cells cocultured with their wild-type counterparts. Recombinant CCL19 increased phosphorylation of STAT5 and induced expression of genes associated with TH2 cell and IL-2 signaling pathways.

CONCLUSIONS

These results reveal a novel, TH2 cell-inducing function of CCL19 in allergic airway disease and suggest that strategies to block this pathway might help to reduce the incidence or severity of allergic asthma.

GLIS3 expression in the thyroid gland in relation to TSH signaling and regulation of gene expression

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time at which GLIS3 target genes, such as Slc5a5 (Nis), become expressed. This, together with observations showing that ubiquitous Glis3KO mice do not display major changes in prenatal thyroid gland morphology, indicated that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of GLIS3 in postnatal thyroid suggested a link between GLIS3 protein expression and blood TSH levels. This was supported by data showing that treatment with TSH, cAMP, or adenylyl cyclase activators or expression of constitutively active PKA enhanced GLIS3 protein stability and transcriptional activity, indicating that GLIS3 activity is regulated at least in part by TSH/TSHR-mediated activation of PKA. The TSH-dependent increase in GLIS3 transcriptional activity would be critical for the induction of GLIS3 target gene expression, including several thyroid hormone (TH) biosynthetic genes, in thyroid follicular cells of mice fed a low iodine diet (LID) when blood TSH levels are highly elevated. Like TH biosynthetic genes, the expression of cell cycle genes is suppressed in ubiquitous Glis3KO mice fed a LID; however, in thyroid-specific Glis3 knockout mice, the expression of cell cycle genes was not repressed, in contrast to TH biosynthetic genes. This indicated that the inhibition of cell cycle genes in ubiquitous Glis3KO mice is dependent on changes in gene expression in GLIS3 target tissues other than the thyroid.

Evaluating Stacked Methylation Markers for Blood-Based Multicancer Detection

The ability to detect several types of cancer using a non-invasive, blood-based test holds the potential to revolutionize oncology screening. We mined tumor methylation array data from the Cancer Genome Atlas (TCGA) covering 14 cancer types and identified two novel, broadly-occurring methylation markers at TLX1 and GALR1. To evaluate their performance as a generalized blood-based screening approach, along with our previously reported methylation biomarker, ZNF154, we rigorously assessed each marker individually or combined. Utilizing TCGA methylation data and applying logistic regression models within each individual cancer type, we found that the three-marker combination significantly increased the average area under the ROC curve (AUC) across the 14 tumor types compared to single markers (p = 1.158 × 10-10; Friedman test). Furthermore, we simulated dilutions of tumor DNA into healthy blood cell DNA and demonstrated increased AUC of combined markers across all dilution levels. Finally, we evaluated assay performance in bisulfite sequenced DNA from patient tumors and plasma, including early-stage samples. When combining all three markers, the assay correctly identified nine out of nine lung cancer plasma samples. In patient plasma from hepatocellular carcinoma, ZNF154 alone yielded the highest combined sensitivity and specificity values averaging 68% and 72%, whereas multiple markers could achieve higher sensitivity or specificity, but not both. Altogether, this study presents a comprehensive pipeline for the identification, testing, and validation of multi-cancer methylation biomarkers with a considerable potential for detecting a broad range of cancer types in patient blood samples.

DNASE1L3 enhances antitumor immunity and suppresses tumor progression in colon cancer

DNASE1L3, an enzyme highly expressed in DCs, is functionally important for regulating autoimmune responses to self-DNA and chromatin. Deficiency of DNASE1L3 leads to development of autoimmune diseases in both humans and mice. However, despite the well-established causal relationship between DNASE1L3 and immunity, little is known about the involvement of DNASE1L3 in regulation of antitumor immunity, the foundation of modern antitumor immunotherapy. In this study, we identify DNASE1L3 as a potentially new regulator of antitumor immunity and a tumor suppressor in colon cancer. In humans, DNASE1L3 is downregulated in tumor-infiltrating DCs, and this downregulation is associated with poor patient prognosis and reduced tumor immune cell infiltration in many cancer types. In mice, Dnase1l3 deficiency in the tumor microenvironment enhances tumor formation and growth in several colon cancer models. Notably, the increased tumor formation and growth in Dnase1l3-deficient mice are associated with impaired antitumor immunity, as evidenced by a substantial reduction of cytotoxic T cells and a unique subset of DCs. Consistently, Dnase1l3-deficient DCs directly modulate cytotoxic T cells in vitro. To our knowledge, our study unveils a previously unknown link between DNASE1L3 and antitumor immunity and further suggests that restoration of DNASE1L3 activity may represent a potential therapeutic approach for anticancer therapy.

Role of GLIS3 in thyroid development and in the regulation of gene expression in thyroid specific Glis3KO mice

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time when GLIS3 target genes, such as Slc5a5 (Nis), become also expressed. We further show that Glis3KO mice do not display any major changes in prenatal thyroid gland morphology indicating that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of thyroid-specific Glis3 knockout (Glis3-Pax8Cre) mice fed either a normal or low-iodine diet (ND or LID) revealed that, in contrast to ubiquitous Glis3KO mice, thyroid follicular cell proliferation and the expression of cell cycle genes were not repressed suggesting that the inhibition of thyroid follicular cell proliferation in ubiquitous Glis3KO mice is related to loss of GLIS3 function in other cell types. However, the expression of several thyroid hormone biosynthesis-, extracellular matrix (ECM)-, and inflammation-related genes was still suppressed in Glis3-Pax8Cre mice particularly under conditions of high blood levels of thyroid stimulating hormone (TSH). We further demonstrate that treatment with TSH, protein kinase A (PKA) or adenylyl cyclase activators or expression of constitutively active PKA enhances GLIS3 protein and activity, suggesting that GLIS3 transcriptional activity is regulated in part by TSH/TSHR-mediated activation of the PKA pathway. This mechanism of regulation provides an explanation for the dramatic increase in GLIS3 protein expression and the subsequent induction of GLIS3 target genes, including several thyroid hormone biosynthetic genes, in thyroid follicular cells of mice fed a LID.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here we determine the roles of CHD4 to enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility at transcription factor binding sites; its depletion leads to altered motif scanning and redistribution of transcription factors to sites not previously occupied. During GATA3-mediated cellular reprogramming, CHD4 activity is necessary to prevent inappropriate chromatin opening and enhancer licensing. Mechanistically, CHD4 competes with transcription factor-DNA interaction by promoting nucleosome positioning over binding motifs. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents inappropriate gene expression by editing binding site selection by transcription factors.

GLIS3 regulates transcription of thyroid hormone biosynthetic genes in coordination with other thyroid transcription factors

BACKGROUND

Loss of the transcription factor GLI-Similar 3 (GLIS3) function causes congenital hypothyroidism (CH) in both humans and mice due to decreased expression of several thyroid hormone (TH) biosynthetic genes in thyroid follicular cells. Whether and to what extent, GLIS3 regulates thyroid gene transcription in coordination with other thyroid transcriptional factors (TFs), such as PAX8, NKX2.1 and FOXE1, is poorly understood.

METHODS

PAX8, NKX2.1, and FOXE1 ChIP-Seq analysis with mouse thyroid glands and rat thyrocyte PCCl3 cells was performed and compared to that of GLIS3 to analyze the co-regulation of gene transcription in thyroid follicular cells by these TFs.

RESULTS

Analysis of the PAX8, NKX2.1, and FOXE1 cistromes identified extensive overlaps between these TF binding loci and those of GLIS3 indicating that GLIS3 shares many of the same regulatory regions with PAX8, NKX2.1, and FOXE1, particularly in genes associated with TH biosynthesis, induced by thyroid stimulating hormone (TSH), and suppressed in Glis3KO thyroid glands, including Slc5a5 (Nis), Slc26a4, Cdh16, and Adm2. ChIP-QPCR analysis showed that loss of GLIS3 did not significantly affect PAX8 or NKX2.1 binding and did not cause major alterations in H3K4me3 and H3K27me3 epigenetic signals.

CONCLUSIONS

Our study indicates that GLIS3 regulates transcription of TH biosynthetic and TSH-inducible genes in thyroid follicular cells in coordination with PAX8, NKX2.1, and FOXE1 by binding within the same regulatory hub. GLIS3 does not cause major changes in chromatin structure at these common regulatory regions. GLIS3 may induce transcriptional activation by enhancing the interaction of these regulatory regions with other enhancers and/or RNA Polymerase II (Pol II) complexes.

Beadchip technology to detect DNA methylation in mouse faithfully recapitulates whole-genome bisulfite sequencing

Aim: To facilitate wide-scale implementation of Illumina Mouse Methylation BeadChip (MMB) technology, array-based measurement of cytosine methylation was compared with the gold-standard assessment of DNA methylation by whole-genome bisulfite sequencing (WGBS). Methods: DNA methylation across two mouse strains (C57B6 and C3H) and both sexes was assessed using the MMB and compared with previously existing deep-coverage WGBS of mice of the same strain and sex. Results & conclusion: The findings demonstrated that 93.3-99.2% of sites had similar measurements of methylation across technologies and that differentially methylated cytosines and regions identified by each technology overlap and enrich for similar biological functions, suggesting that the MMB faithfully recapitulates the findings of WGBS.

Contribution of the Wolffian duct mesenchyme to the formation of the female reproductive tract

The female reproductive tract develops from its embryonic precursor, the Müllerian duct. In close proximity to the Müllerian duct lies the precursor for the male reproductive tract, the Wolffian duct, which is eliminated in the female embryo during sexual differentiation. We discovered that a component of the Wolffian duct, its mesenchyme, is not eliminated after sexual differentiation. Instead, the Wolffian duct mesenchyme underwent changes in transcriptome and chromatin accessibility from male tract to female tract identity, and became a unique mesenchymal population in the female reproductive tract with localization and transcriptome distinct from the mesenchyme derived from the Müllerian duct. Partial ablation of the Wolffian duct mesenchyme stunted the growth of the fetal female reproductive tract in ex vivo organ culture. These findings reveal a new fetal origin of mesenchymal tissues for female reproductive tract formation and reshape our understanding of sexual differentiation of reproductive tracts.

CUT&Tag-BS for simultaneous profiling of histone modification and DNA methylation with high efficiency and low cost

It remains a challenge to decipher the complex relationship between DNA methylation, histone modification, and the underlying DNA sequence with limited input material. Here, we developed an efficient, low-input, and low-cost method for the simultaneous profiling of genomic localization of histone modification and methylation status of the underlying DNA at single-base resolution from the same cells in a single experiment by integrating cleavage under targets and tagmentation (CUT&Tag) with tagmentation-based bisulfite sequencing (CUT&Tag-BS). We demonstrated the validity of our method using representative histone modifications of euchromatin and constitutive and facultative heterochromatin (H3K4me1, H3K9me3, and H3K27me3, respectively). Similar histone modification enrichment patterns were observed in CUT&Tag-BS compared with non-bisulfite-treated control, and H3K4me1-marked regions were found to mostly be CpG poor, lack methylation concordance, and exhibit prevalent DNA methylation heterogeneity among mouse embryonic stem cells (mESCs). We anticipate that CUT&Tag-BS will be widely applied to directly address the genomic relationship between DNA methylation and histone modification, especially in low-input scenarios with precious biological samples.

Low-input ATAC&mRNA-seq protocol for simultaneous profiling of chromatin accessibility and gene expression

We present a simple, fast, and robust protocol (low-input ATAC&mRNA-seq) to simultaneously generate ATAC-seq and mRNA-seq libraries from the same cells in limited cell numbers by coupling a simplified ATAC procedure using whole cells with a novel mRNA-seq approach that features a seamless on-bead process including direct mRNA isolation from the cell lysate, solid-phase cDNA synthesis, and direct tagmentation of mRNA/cDNA hybrids for library preparation. It enables dual-omics profiling from limited material when joint epigenome and transcriptome analyses are needed.

For complete details on the use and execution of this protocol, please refer to Li et al. (2021).

•

A simple and fast dual-omics profiling protocol for low-input samples

•

Simultaneous measurement of open chromatin and gene expression from the same cells

•

Suitable for limited numbers of cells ranging from 5,000 to 20,000 cells

GLIS1 regulates trabecular meshwork function and intraocular pressure and is associated with glaucoma in humans

Chronically elevated intraocular pressure (IOP) is the major risk factor of primary open-angle glaucoma, a leading cause of blindness. Dysfunction of the trabecular meshwork (TM), which controls the outflow of aqueous humor (AqH) from the anterior chamber, is the major cause of elevated IOP. Here, we demonstrate that mice deficient in the Krüppel-like zinc finger transcriptional factor GLI-similar-1 (GLIS1) develop chronically elevated IOP. Magnetic resonance imaging and histopathological analysis reveal that deficiency in GLIS1 expression induces progressive degeneration of the TM, leading to inefficient AqH drainage from the anterior chamber and elevated IOP. Transcriptome and cistrome analyses identified several glaucoma- and extracellular matrix-associated genes as direct transcriptional targets of GLIS1. We also identified a significant association between GLIS1 variant rs941125 and glaucoma in humans (P = 4.73 × 10-6), further supporting a role for GLIS1 into glaucoma etiology. Our study identifies GLIS1 as a critical regulator of TM function and maintenance, AqH dynamics, and IOP.

A simple and robust method for simultaneous dual-omics profiling with limited numbers of cells

Deciphering epigenetic regulation of gene expression requires measuring the epigenome and transcriptome jointly. Single-cell multi-omics technologies have been developed for concurrent profiling of chromatin accessibility and gene expression. However, multi-omics profiling of low-input bulk samples remains challenging. Therefore, we developed low-input ATAC&mRNA-seq, a simple and robust method for studying the role of chromatin structure in gene regulation in a single experiment with thousands of cells, to maximize insights from limited input material by obtaining ATAC-seq and mRNA-seq data simultaneously from the same cells with data quality comparable to that of conventional mono-omics assays. Integrative data analysis revealed similar strong association between promoter accessibility and gene expression when using the data of low-input ATAC&mRNA-seq as when using single-assay data, underscoring the accuracy and reliability of our dual-omics assay to generate both datum types simultaneously with just thousands of cells. We envision our method to be widely applied in many biological disciplines with limited materials.

Single Nucleotide Resolution Analysis Reveals Pervasive, Long-Lasting DNA Methylation Changes by Developmental Exposure to a Mitochondrial Toxicant

Mitochondrial-driven alterations of the epigenome have been reported, but whether they are relevant at the organismal level remains unknown. The viable yellow agouti mouse (A^vy) is a powerful epigenetic biosensor model that reports on the DNA methylation status of the A^vy locus, which is established prior to the three-germ-layer separation, through the coat color of the animals. Here we show that maternal exposure to rotenone, a potent mitochondrial complex I inhibitor, not only changes the DNA methylation status of the A^vy locus in the skin but broadly affects the liver DNA methylome of the offspring. These effects are accompanied by altered gene expression programs that persist throughout life, and which associate with impairment of antioxidant activity and mitochondrial function in aged animals. These pervasive and lasting genomic effects suggest a putative role for mitochondria in regulating life-long gene expression programs through developmental nuclear epigenetic remodeling.

Lozoya et al. provide in vivo evidence of the epigenetic effects of mitochondrial dysfunction. Developmental-only exposure to rotenone through the mother’s diet inhibits mitochondrial complex I in the dams and results in lifelong nuclear DNA methylation and gene expression changes in the offspring. Aged offspring also show functional outcomes.

UDP-glucose and P2Y14 receptor amplify allergen-induced airway eosinophilia

Airway eosinophilia is a hallmark of allergic asthma and is associated with mucus production, airway hyperresponsiveness, and shortness of breath. Although glucocorticoids are widely used to treat asthma, their prolonged use is associated with several side effects. Furthermore, many individuals with eosinophilic asthma are resistant to glucocorticoid treatment, and they have an unmet need for novel therapies. Here, we show that UDP-glucose (UDP-G), a nucleotide sugar, is selectively released into the airways of allergen-sensitized mice upon their subsequent challenge with that same allergen. Mice lacking P2Y14R, the receptor for UDP-G, had decreased airway eosinophilia and airway hyperresponsiveness compared with wild-type mice in a protease-mediated model of asthma. P2Y14R was dispensable for allergic sensitization and for the production of type 2 cytokines in the lung after challenge. However, UDP-G increased chemokinesis in eosinophils and enhanced their response to the eosinophil chemoattractant, CCL24. In turn, eosinophils triggered the release of UDP-G into the airway, thereby amplifying eosinophilic recruitment. This positive feedback loop was sensitive to therapeutic intervention, as a small molecule antagonist of P2Y14R inhibited airway eosinophilia. These findings thus reveal a pathway that can be therapeutically targeted to treat asthma exacerbations and glucocorticoid-resistant forms of this disease.

Molecular Actions Underlying Wolffian Duct Regression in Sexual Differentiation of Murine Reproductive Tracts

Sexually dimorphic establishment of the reproductive tract system requires sex-specific regression of the Wolffian duct and Müllerian duct in the mesonephros. In an XX embryo, the Wolffian duct regresses under the control of the mesenchymal transcription factor COUP-TFII. To understand cellular and molecular actions underlying Wolffian duct regression, we performed transcriptomic analyses of XX mesonephroi with or without Coup-tfII and genome-wide analysis of COUP-TFII chromatin occupancy in XX mesonephroi. The integrative analysis of COUP-TFII genome-wide binding and transcriptomic analysis revealed the suppression of muscle differentiation and extracellular matrix genes by COUP-TFII and identified a group of potential transcriptional partners of COUP-TFII in the mesenchyme that potentially facilitate Wolffian duct regression. These findings provide insights into the molecular action of COUP-TFII in the Wolffian duct mesenchyme and identify a list of biologically relevant candidate genes and pathways for future functional analyses in sexual differentiation of reproductive tracts.

The genomic regulatory elements for estrogen receptor alpha transactivation-function-1 regulated genes

Estrogen receptor alpha (ERα) is a ligand-dependent transcription regulator, containing two transactivation functional domains, AF-1 and AF-2. The selective estrogen receptor modulators (SERMs), including 4-hydroxytamoxifen (4OHT), activate AF-1 preferentially rather than AF-2. However, it is unclear whether this specific function is related to the tissue-selective functionality of SERMs. Moreover, there is no information determining AF-1-dependent estrogenic-genes existing in tissues. We sought to identify AF-1-dependent estrogenic-genes using the AF-2 mutated knock-in (KI) mouse model, AF2ERKI. AF2ER is an AF-2 disrupted estradiol (E2)-insensitive mutant ERα, but AF-1-dependent transcription can be activated by the estrogen-antagonists, fulvestrant (ICI) and 4OHT. Gene profiling and ChIP-Seq analysis identified Klk1b21 as an ICI-inducible gene in AF2ERKI uterus. The regulatory activity was analyzed further using a cell-based reporter assay. The 5'-flanking 0.4k bp region of Klk1b21 gene responded as an ERα AF-1-dependent estrogen-responsive promoter. The 150 bp minimum ERα binding element (EBE) consists of three direct repeats. These three half-site sequences were essential for the ERα-dependent transactivation and were differentially recognized by E2 and 4OHT for the gene activation. This response was impaired when the minimum EBE was fused with a thymidine-kinase promoter but could be restored by fusion with the 100 bp minimum transcription initiation element (TIE) of Klk1b21, suggesting that the cooperative function of EBE and TIE is essential for mediating AF-1-dependent transactivation. These findings provide the first in vivo evidence that endogenous ERα AF-1 dominant estrogenic-genes exist in estrogen-responsive organs. Such findings will aid in understanding the mechanism of ERα-dependent tissue-selective activity of SERMs.

Cancer-specific mutation of GATA3 disrupts the transcriptional regulatory network governed by Estrogen Receptor alpha, FOXA1 and GATA3

Estrogen receptors (ER) are activated by the steroid hormone 17β-estradiol. Estrogen receptor alpha (ER-α) forms a regulatory network in mammary epithelial cells and in breast cancer with the transcription factors FOXA1 and GATA3. GATA3 is one of the most frequently mutated genes in breast cancer and is capable of specifying chromatin localization of FOXA1 and ER-α. How GATA3 mutations found in breast cancer impact genomic localization of ER-α and the transcriptional network downstream of ER-α and FOXA1 remains unclear. Here, we investigate the function of a recurrent patient-derived GATA3 mutation (R330fs) on this regulatory network. Genomic analysis indicates that the R330fs mutant can disrupt localization of ER-α and FOXA1. Loci co-bound by all three factors are enriched for genes integral to mammary gland development as well as epithelial cell biology. This gene set is differentially regulated in GATA3 mutant cells in culture and in tumors bearing similar mutations in vivo. The altered distribution of ER-α and FOXA1 in GATA3-mutant cells is associated with altered chromatin architecture, which leads to differential gene expression. These results suggest an active role for GATA3 zinc finger 2 mutants in ER-α positive breast tumors.

Interaction of the pioneer transcription factor GATA3 with nucleosomes

During cellular reprogramming, the pioneer transcription factor GATA3 binds chromatin, and in a context-dependent manner directs local chromatin remodeling and enhancer formation. Here, we use high-resolution nucleosome mapping in human cells to explore the impact of the position of GATA motifs on the surface of nucleosomes on productive enhancer formation, finding productivity correlates with binding sites located near the nucleosomal dyad axis. Biochemical experiments with model nucleosomes demonstrate sufficiently stable transcription factor-nucleosome interaction to empower cryo-electron microscopy structure determination of the complex at 3.15 Å resolution. The GATA3 zinc fingers efficiently bind their target 5′-GAT-3′ sequences in the nucleosome when they are located in solvent accessible, consecutive major grooves without significant changes in nucleosome structure. Analysis of genomic loci bound by GATA3 during reprogramming suggests a correlation of recognition motif sequence and spacing that may distinguish productivity of new enhancer formation.

GATA 3 functions as a pioneer factor during cellular reprogramming. Here the authors delineate nucleosome positioning relative to GATA3 binding motifs and describe the structure of a GATA3–nucleosome complex; providing insight into how a pioneer factor interacts with nucleosomes and catalyze their local remodelling to produce an accessible enhancer.

Comprehensive structure-function characterization of DNMT3B and DNMT3A reveals distinctive de novo DNA methylation mechanisms

Mammalian DNA methylation patterns are established by two de novo DNA methyltransferases, DNMT3A and DNMT3B, which exhibit both redundant and distinctive methylation activities. However, the related molecular basis remains undetermined. Through comprehensive structural, enzymology and cellular characterization of DNMT3A and DNMT3B, we here report a multi-layered substrate-recognition mechanism underpinning their divergent genomic methylation activities. A hydrogen bond in the catalytic loop of DNMT3B causes a lower CpG specificity than DNMT3A, while the interplay of target recognition domain and homodimeric interface fine-tunes the distinct target selection between the two enzymes, with Lysine 777 of DNMT3B acting as a unique sensor of the +1 flanking base. The divergent substrate preference between DNMT3A and DNMT3B provides an explanation for site-specific epigenomic alterations seen in ICF syndrome with DNMT3B mutations. Together, this study reveals distinctive substrate-readout mechanisms of the two DNMT3 enzymes, implicative of their differential roles during development and pathogenesis.

Estrogen receptor α (ERα)-binding super-enhancers drive key mediators that control uterine estrogen responses in mice

Estrogen receptor α (ERα) modulates gene expression by interacting with chromatin regions that are frequently distal from the promoters of estrogen-regulated genes. Active chromatin-enriched "super-enhancer" (SE) regions, mainly observed in in vitro culture systems, often control production of key cell type-determining transcription factors. Here, we defined super-enhancers that bind to ERα in vivo within hormone-responsive uterine tissue in mice. We found that SEs are already formed prior to estrogen exposure at the onset of puberty. The genes at SEs encoded critical developmental factors, including retinoic acid receptor α (RARA) and homeobox D (HOXD). Using high-throughput chromosome conformation capture (Hi-C) along with DNA sequence analysis, we demonstrate that most SEs are located at a chromatin loop end and that most uterine genes in loop ends associated with these SEs are regulated by estrogen. Although the SEs were formed before puberty, SE-associated genes acquired optimal ERα-dependent expression after reproductive maturity, indicating that pubertal processes that occur after SE assembly and ERα binding are needed for gene responses. Genes associated with these SEs affected key estrogen-mediated uterine functions, including transforming growth factor β (TGFβ) and LIF interleukin-6 family cytokine (LIF) signaling pathways. To the best of our knowledge, this is the first identification of SE interactions that underlie hormonal regulation of genes in uterine tissue and optimal development of estrogen responses in this tissue.

Alterations in promoter interaction landscape and transcriptional network underlying metabolic adaptation to diet

Metabolic adaptation to nutritional state requires alterations in gene expression in key tissues. Here, we investigated chromatin interaction dynamics, as well as alterations in cis-regulatory loci and transcriptional network in a mouse model system. Chronic consumption of a diet high in saturated fat, when compared to a diet high in carbohydrate, led to dramatic reprogramming of the liver transcriptional network. Long-range interaction of promoters with distal regulatory loci, monitored by promoter capture Hi-C, was regulated by metabolic status in distinct fashion depending on diet. Adaptation to a lipid-rich diet, mediated largely by nuclear receptors including Hnf4α, relied on activation of preformed enhancer/promoter loops. Adaptation to carbohydrate-rich diet led to activation of preformed loops and to de novo formation of new promoter/enhancer interactions. These results suggest that adaptation to nutritional changes and metabolic stress occurs through both de novo and pre-existing chromatin interactions which respond differently to metabolic signals.

Metabolic adaptation to different diets results in changes to gene expression. Here, the authors characterise the chromatin landscape and transcriptional network in mice on a diet of high saturated fat, compared to a diet high in carbohydrate, finding a dramatic reprogramming of the liver transcriptional network.

Transcriptome and DNA Methylome Analysis in a Mouse Model of Diet-Induced Obesity Predicts Increased Risk of Colorectal Cancer

Colorectal cancer (CRC) tends to occur at older age; however, CRC incidence rates have been rising sharply among young age groups. The increasing prevalence of obesity is recognized as a major risk, yet the mechanistic underpinnings remain poorly understood. Using a diet-induced obesity mouse model, we identified obesity-associated molecular changes in the colonic epithelium of young and aged mice, and we further investigated whether the changes were reversed after weight loss. Transcriptome analysis indicated that obesity-related colonic cellular metabolic switch favoring long-chain fatty acid oxidation happened in young mice, while obesity-associated downregulation of negative feedback regulators of pro-proliferative signaling pathways occurred in older mice. Strikingly, colonic DNA methylome was pre-programmed by obesity at young age, priming for a tumor-prone gene signature after aging. Furthermore, obesity-related changes were substantially preserved after short-term weight loss, but they were largely reversed after long-term weight loss. We provided mechanistic insights into increased CRC risk in obesity.

Revealing a human p53 universe

p53 transcriptional networks are well-characterized in many organisms. However, a global understanding of requirements for in vivo p53 interactions with DNA and relationships with transcription across human biological systems in response to various p53 activating situations remains limited. Using a common analysis pipeline, we analyzed 41 data sets from genome-wide ChIP-seq studies of which 16 have associated gene expression data, including our recent primary data with normal human lymphocytes. The resulting extensive analysis, accessible at p53 BAER hub via the UCSC browser, provides a robust platform to characterize p53 binding throughout the human genome including direct influence on gene expression and underlying mechanisms. We establish the impact of spacers and mismatches from consensus on p53 binding in vivo and propose that once bound, neither significantly influences the likelihood of expression. Our rigorous approach revealed a large p53 genome-wide cistrome composed of >900 genes directly targeted by p53. Importantly, we identify a core cistrome signature composed of genes appearing in over half the data sets, and we identify signatures that are treatment- or cell-specific, demonstrating new functions for p53 in cell biology. Our analysis reveals a broad homeostatic role for human p53 that is relevant to both basic and translational studies.

A distal super enhancer mediates estrogen-dependent mouse uterine-specific gene transcription of Igf1 (insulin-like growth factor 1)

Insulin-like growth factor 1 (IGF1) is primarily synthesized in and secreted from the liver; however, estrogen (E2), through E2 receptor α (ERα), increases uterine Igf1 mRNA levels. Previous ChIP-seq analyses of the murine uterus have revealed a potential enhancer region distal from the Igf1 transcription start site (TSS) with multiple E2-dependent ERα-binding regions. Here, we show E2-dependent super enhancer-associated characteristics and suggest contact between the distal enhancer and the Igf1 TSS. We hypothesized that this distal super-enhancer region controls E2-responsive induction of uterine Igf1 transcripts. We deleted 430 bp, encompassing one of the ERα-binding sites, thereby disrupting interactions of the enhancer with gene-regulatory factors. As a result, E2-mediated induction of mouse uterine Igf1 mRNA is completely eliminated, whereas hepatic Igf1 expression remains unaffected. This highlights the central role of a distal enhancer in the assembly of the factors necessary for E2-dependent interaction with the Igf1 TSS and induction of uterus-specific Igf1 transcription. Of note, loss of the enhancer did not affect fertility or uterine growth responses. Deletion of uterine Igf1 in a PgrCre;Igf1f/f model decreased female fertility but did not impact the E2-induced uterine growth response. Moreover, E2-dependent activation of uterine IGF1 signaling was not impaired by disrupting the distal enhancer or by deleting the coding transcript. This indicated a role for systemic IGF1, suggested that other growth mediators drive uterine response to E2, and suggested that uterine-derived IGF1 is essential for reproductive success. Our findings elucidate the role of a super enhancer in Igf1 regulation and uterine growth.

Genome-wide identification of FOXL2 binding and characterization of FOXL2 feminizing action in the fetal gonads

The identity of the gonads is determined by which fate, ovarian granulosa cell or testicular Sertoli cell, the bipotential somatic cell precursors choose to follow. In most vertebrates, the conserved transcription factor FOXL2 contributes to the fate of granulosa cells. To understand FOXL2 functions during gonad differentiation, we performed genome-wide analysis of FOXL2 chromatin occupancy in fetal ovaries and established a genetic mouse model that forces Foxl2 expression in the fetal testis. When FOXL2 was ectopically expressed in the somatic cell precursors in the fetal testis, FOXL2 was sufficient to repress Sertoli cell differentiation, ultimately resulting in partial testis-to-ovary sex-reversal. Combining genome-wide analysis of FOXL2 binding in the fetal ovary with transcriptomic analyses of our Foxl2 gain-of-function and previously published Foxl2 loss-of-function models, we identified potential pathways responsible for the feminizing action of FOXL2. Finally, comparison of FOXL2 genome-wide occupancy in the fetal ovary with testis-determining factor SOX9 genome-wide occupancy in the fetal testis revealed extensive overlaps, implying that antagonistic signals between FOXL2 and SOX9 occur at the chromatin level.

Ly-6C+CD11b+conventional dendritic cells accumulate in inflamed lung and differentiate into type 2 dendritic cells

Allergic asthma is a complex inflammatory disease characterized by airway obstruction and airway hyperresponsiveness. Pulmonary dendritic cells are critical for allergic sensitization as they take up inhaled allergens and direct differentiation of allergen-specific helper T cell lineages in lung-draining lymph nodes. In this study, we used mass cytometry to study lung DC subsets that increase in frequency after allergen inhalation. This approach revealed a novel cell population displaying Ly-6C, in addition to surface markers routinely used to identify type 2 conventional DCs (cDC2), including CD11b and CD26. These Ly-6C+ cells had dendrites, lacked monocyte/macrophage markers such as F4/80 and CD115, and were absent in the mice lacking FMS-like tyrosine kinase 3ligand (FLT3L), a growth factor essential for cDC differentiation. Based on these observations, we conclude that this novel population represents bona fide cDCs. Their numbers in the lung were low at steady state, but increased dramatically after allergic sensitization. Cultured DC precursors (preDCs) displayed Ly-6C, but lost that marker by time. Ly-6C+CD11b+cDCs freshly isolated from mouse lung also lost Ly-6C expression after ex vivo culture. Together, these results suggest that Ly-6C+CD11b+cDCs are the precursor of cDC2. Given that the number of these DCs increase dramatically in the lung after allergen inhalation, they might play important roles in allergic sensitization and the induction of allergic asthma.

Crohn's disease IRGM risk alleles are associated with altered gene expression in human tissues

Crohn's disease (CD) is a chronic inflammatory gastrointestinal disorder. Genetic association studies have implicated dysregulated autophagy in CD. Among risk loci identified are a promoter single nucleotide polymorphism (SNP)( rs13361189 ) and two intragenic SNPs ( rs9637876 , rs10065172 ) in immunity-related GTPase family M ( IRGM) a gene that encodes a protein of the autophagy initiation complex. All three SNPs have been proposed to modify IRGM expression, but reports have been divergent and largely derived from cell lines. Here, analyzing RNA-Sequencing data of human tissues from the Genotype-Tissue Expression Project, we found that rs13361189 minor allele carriers had reduced IRGM expression in whole blood and terminal ileum, and upregulation in ileum of ZNF300P1, a locus adjacent to IRGM on chromosome 5q33.1 that encodes a long noncoding RNA. Whole blood and ileum from minor allele carriers had altered expression of multiple additional genes that have previously been linked to colitis and/or autophagy. Notable among these was an increase in ileum of LTF (lactoferrin), an established fecal inflammatory biomarker of CD, and in whole blood of TNF, a key cytokine in CD pathogenesis. Last, we confirmed that risk alleles at all three loci associated with increased risk for CD but not ulcerative colitis in a case-control study. Taken together, our findings suggest that genetically encoded IRGM deficiency may predispose to CD through dysregulation of inflammatory gene networks. Gene expression profiling of disease target tissues in genetically susceptible populations is a promising strategy for revealing new leads for the study of molecular pathogenesis and, potentially, for precision medicine. NEW & NOTEWORTHY Single nucleotide polymorphisms in immunity-related GTPase family M ( IRGM), a gene that encodes an autophagy initiation protein, have been linked epidemiologically to increased risk for Crohn's disease (CD). Here, we show for the first time that subjects with risk alleles at two such loci, rs13361189 and rs10065172 , have reduced IRGM expression in whole blood and terminal ileum, as well as dysregulated expression of a wide array of additional genes that regulate inflammation and autophagy.

Abstract 5418: GATA3 zinc-finger mutation induces transcriptional reprogramming in breast cancer through gain and loss of function

Recent large-scale breast cancer genomic profiling identified frequent mutations in a pioneer transcription factor, GATA3. The expression of GATA3 is a prominent marker of luminal breast tumors, and loss of GATA3 expression is associated with aggressive tumor phenotypes. In addition, GATA3 mutations are frequently found in these luminal tumors. However, the clinical and molecular outcomes of GATA3 mutations are poorly understood. Our novel classification scheme of GATA3 mutations defines distinct clinical features of breast tumors with mutations in the second GATA3 zinc-finger (ZnFn2): frequent observation in luminal B subtype and worse prognosis than other GATA3 mutant cases. To dissect the molecular impacts of these ZnFn2 mutations, we generated a GATA3 mutant breast cancer cell clone by CRISPR-Cas9 gene editing. An engineered ZnFn2 mutant cell line manifested increased tumor growth as xenografts and more aggressive phenotypes in vitro. The ZnFn2 mutation led to loss of GATA3 binding and decreased expression at a subset of genes, including progesterone receptor (PR). The mutant cells also exhibited gain of GATA3 binding at other loci, correlated with increased expression of mesenchymal marker genes. Decreased expression of PR and impaired response to progesterone were crucial for cancer-promoting functions in GATA3 ZnFn2 mutant cells. Downregulation of PR and its downstream genes was also observed in the clinical gene expression data. These results illuminate tumor-promoting functions of GATA3 ZnFn2 mutations in breast cancer.

Citation Format: Motoki Takaku, Sara A. Grimm, John D. Roberts, Kaliopi Chrysovergis, Page Myers, Charles M. Perou, Paul A. Wade. GATA3 zinc-finger mutation induces transcriptional reprogramming in breast cancer through gain and loss of function [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5418.

Loss of Glis3 causes dysregulation of retrotransposon silencing and germ cell demise in fetal mouse testis

Fetal germ cell development is regulated by an elaborate combination of cell-extrinsic and cell-intrinsic signals. Here we identify a novel role for the Krüppel-like transcription factor Gli-Similar 3 (Glis3) in male germ cell development in the mouse embryos. Glis3 is expressed in male germ cells during the brief window of time prior to initiation of piRNA-dependent retrotransposon surveillance. Disruption of Glis3 function led to a widespread reduction in retrotransposon silencing factors, aberrant retrotransposon expression and pronounced germ cell loss. Experimental induction of precocious Glis3 expression in vivo before its normal expression resulted in premature expression of several piRNA pathway members, suggesting that GLIS3 is necessary for the activation of the retrotransposon silencing programs. Our findings reveal an unexpected role for GLIS3 in the development of male germ cells and point to a central role for GLIS3 in the control of retrotransposon silencing in the fetal germline.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

GATA3 zinc finger 2 mutations reprogram the breast cancer transcriptional network

GATA3 is frequently mutated in breast cancer; these mutations are widely presumed to be loss-of function despite a dearth of information regarding their effect on disease course or their mechanistic impact on the breast cancer transcriptional network. Here, we address molecular and clinical features associated with GATA3 mutations. A novel classification scheme defines distinct clinical features for patients bearing breast tumors with mutations in the second GATA3 zinc-finger (ZnFn2). An engineered ZnFn2 mutant cell line by CRISPR–Cas9 reveals that mutation of one allele of the GATA3 second zinc finger (ZnFn2) leads to loss of binding and decreased expression at a subset of genes, including Progesterone Receptor. At other loci, associated with epithelial to mesenchymal transition, gain of binding correlates with increased gene expression. These results demonstrate that not all GATA3 mutations are equivalent and that ZnFn2 mutations impact breast cancer through gain and loss-of function.

In breast cancer GATA3 is known to be frequently mutated, but the function of these mutations is unclear. Here, the authors utilise CRISPR-Cas9 to model frame-shift mutations in zinc finger 2 of GATA3, highlighting that GATA3 mutation can have gain- or loss-of function effects in breast cancer.

An obesity-associated gut microbiome reprograms the intestinal epigenome and leads to altered colonic gene expression

Background

The gut microbiome, a key constituent of the colonic environment, has been implicated as an important modulator of human health. The eukaryotic epigenome is postulated to respond to environmental stimuli through alterations in chromatin features and, ultimately, gene expression. How the host mediates epigenomic responses to gut microbiota is an emerging area of interest. Here, we profile the gut microbiome and chromatin characteristics in colon epithelium from mice fed either an obesogenic or control diet, followed by an analysis of the resultant changes in gene expression.

Results

The obesogenic diet shapes the microbiome prior to the development of obesity, leading to altered bacterial metabolite production which predisposes the host to obesity. This microbiota–diet interaction leads to changes in histone modification at active enhancers that are enriched for binding sites for signal responsive transcription factors. These alterations of histone methylation and acetylation are associated with signaling pathways integral to the development of colon cancer. The transplantation of obesogenic diet-conditioned microbiota into germ free mice, combined with an obesogenic diet, recapitulates the features of the long-term diet regimen. The diet/microbiome-dependent changes are reflected in both the composition of the recipient animals’ microbiome as well as in the set of transcription factor motifs identified at diet-influenced enhancers.

Conclusions

These findings suggest that the gut microbiome, under specific dietary exposures, stimulates a reprogramming of the enhancer landscape in the colon, with downstream effects on transcription factors. These chromatin changes may be associated with those seen during colon cancer development.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1389-1) contains supplementary material, which is available to authorized users.

Role of ERα in Mediating Female Uterine Transcriptional Responses to IGF1

Estrogen (E2) signaling through its nuclear receptor, E2 receptor α (ERα) increases insulinlike growth factor 1 (IGF1) in the rodent uterus, which then initiates further signals via the IGF1 receptor. Directly administering IGF1 results in similar biological and transcriptional uterine responses. Our studies using global ERα-null mice demonstrated a loss of uterine biological responses of the uterus to E2 or IGF1 treatment, while maintaining transcriptional responses to IGF1. To address this discrepancy in the need for uterine ERα in mediating the IGF1 transcriptional vs growth responses, we assessed the IGF1 transcriptional responses in PgrCre+Esr1f/f (called ERαUtcKO) mice, which selectively lack ERα in progesterone receptor (PGR) expressing cells, including all uterine cells, while maintaining ERα expression in other tissues and cells that do not express Pgr. Additionally, we profiled IGF1-induced ERα binding sites in uterine chromatin using chromatin immunoprecipitation sequencing. Herein, we explore the transcriptional and molecular signaling that underlies our findings to refine our understanding of uterine IGF1 signaling and identify ERα-mediated and ERα-independent uterine transcriptional responses. Defining these mechanisms in vivo in whole tissue and animal contexts provides details of nuclear receptor mediated mechanisms that impact biological systems and have potential applicability to reproductive processes of humans, livestock and wildlife.

The novel p53 target TNFAIP8 variant 2 is increased in cancer and offsets p53-dependent tumor suppression

Tumor necrosis factor-α-induced protein 8 (TNFAIP8) is a stress-response gene that has been associated with cancer, but no studies have differentiated among or defined the regulation or function of any of its several recently described expression variants. We found that TNFAIP8 variant 2 (v2) is overexpressed in multiple human cancers, whereas other variants are commonly downregulated in cancer (v1) or minimally expressed in cancer or normal tissue (v3-v6). Silencing v2 in cancer cells induces p53-independent inhibition of DNA synthesis, widespread binding of p53, and induction of target genes and p53-dependent cell cycle arrest and DNA damage sensitization. Cell cycle arrest induced by v2 silencing requires p53-dependent induction of p21. In response to the chemotherapeutic agent doxorubicin, p53 regulates v2 through binding to an intragenic enhancer, together indicating that p53 and v2 engage in complex reciprocal regulation. We propose that TNFAIP8 v2 promotes human cancer by broadly repressing p53 function, in essence offsetting p53-dependent tumor suppression.

The Impact of Environmental and Endogenous Damage on Somatic Mutation Load in Human Skin Fibroblasts

Accumulation of somatic changes, due to environmental and endogenous lesions, in the human genome is associated with aging and cancer. Understanding the impacts of these processes on mutagenesis is fundamental to understanding the etiology, and improving the prognosis and prevention of cancers and other genetic diseases. Previous methods relying on either the generation of induced pluripotent stem cells, or sequencing of single-cell genomes were inherently error-prone and did not allow independent validation of the mutations. In the current study we eliminated these potential sources of error by high coverage genome sequencing of single-cell derived clonal fibroblast lineages, obtained after minimal propagation in culture, prepared from skin biopsies of two healthy adult humans. We report here accurate measurement of genome-wide magnitude and spectra of mutations accrued in skin fibroblasts of healthy adult humans. We found that every cell contains at least one chromosomal rearrangement and 600–13,000 base substitutions. The spectra and correlation of base substitutions with epigenomic features resemble many cancers. Moreover, because biopsies were taken from body parts differing by sun exposure, we can delineate the precise contributions of environmental and endogenous factors to the accrual of genetic changes within the same individual. We show here that UV-induced and endogenous DNA damage can have a comparable impact on the somatic mutation loads in skin fibroblasts.

Somatic genomes are constantly accumulating changes caused by endogenous lesions, errors in DNA replication and repair, as well as environmental insults. Despite the importance of somatic genome instability in aging and age-related pathologies, including cancers, accurate measurements of mutation loads in healthy cells is still missing. In this study, we developed an experimental approach to accurately determine the somatic genome changes accrued in cell lineages over the lifetime of healthy humans. We show that the amounts and types of mutations in skin cells resemble many cancers, thus indicating that the mechanisms that lead to carcinogenesis are also functional in healthy cells. Moreover, sun-exposed skin cells have a higher mutation load attributable to ultraviolet radiation (UV) unlike cells from hips that were protected by clothing. Our work provides precise measurements of the mutation loads in single cells in human skin. Furthermore our data allowed defining the mutagenic impacts of environmental and endogenous processes within the same individual and led to conclusion that these processes have a comparable impact on the somatic mutation load.

GATA3-dependent cellular reprogramming requires activation-domain dependent recruitment of a chromatin remodeler

Background

Transcription factor-dependent cellular reprogramming is integral to normal development and is central to production of induced pluripotent stem cells. This process typically requires pioneer transcription factors (TFs) to induce de novo formation of enhancers at previously closed chromatin. Mechanistic information on this process is currently sparse.

Results

Here we explore the mechanistic basis by which GATA3 functions as a pioneer TF in a cellular reprogramming event relevant to breast cancer, the mesenchymal to epithelial transition (MET). In some instances, GATA3 binds previously inaccessible chromatin, characterized by stable, positioned nucleosomes where it induces nucleosome eviction, alters local histone modifications, and remodels local chromatin architecture. At other loci, GATA3 binding induces nucleosome sliding without concomitant generation of accessible chromatin. Deletion of the transactivation domain retains the chromatin binding ability of GATA3 but cripples chromatin reprogramming ability, resulting in failure to induce MET.

Conclusions

These data provide mechanistic insights into GATA3-mediated chromatin reprogramming during MET, and suggest unexpected complexity to TF pioneering. Successful reprogramming requires stable binding to a nucleosomal site; activation domain-dependent recruitment of co-factors including BRG1, the ATPase subunit of the SWI/SNF chromatin remodeling complex; and appropriate genomic context. The resulting model provides a new conceptual framework for de novo enhancer establishment by a pioneer TF.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0897-0) contains supplementary material, which is available to authorized users.

Abstract 2862: MBD3 regulates chromatin accessibility at active promoters

Chromatin structure is tightly regulated in cells and its dysregulation is associated with diseases such as cancer. The Mi-2/NuRD (Nucleosome Remodeling Deacetylase) complex is postulated to organize chromatin structure using its nucleosome remodeling and histone deacetylase activities. MBD3 is an integral component of NuRD, which potentially targets the complex to the specific sites of the genome. Recently, we have demonstrated that MBD3/NuRD targets regulatory elements of active genes, a finding diametrically opposed to historical models depicting NuRD as a co-repressor localized to transcriptionally non-permissive chromatin. Moreover, MBD3 co-localized with RNA polymerase II (pol II) at loci where promoter-enhancer loops are formed. Although our previous findings shed light on the role of NuRD in the regulation of active genes, the detailed mechanism of how NuRD is involved in transcriptional regulation is yet not fully understood. To investigate this question, we have developed a next generation, tagmentation based chromatin immunoprecipitation method followed by massively parallel sequencing (ChIP-nexo), which substantially increased the resolution of the MBD3/NuRD mapping. ChIP-nexo defined MBD3/NuRD localization at high resolution, showing accumulation of MBD3 at active promoters with a dip at transcription start site, suggesting an active role of NuRD in modulation of the cell-type specific transcriptional network. To interrogate the detailed functions of NuRD at active promoters, we conducted MNase-seq using conditions which measures both nucleosome positioning and sensitivity, in MBD3 depleted cells. Depletion of MBD3 changed the accessibility of nucleosomes at promoters, while maintaining overall nucleosome positioning. These data indicate an active regulatory function of NuRD in fine tuning the transcriptional network by modulating transcription rates of pol II. We propose that NuRD may be involved in establishing cell type specific gene expression patterns in diverse cell types by organizing promoter nucleoprotein architecture.

Citation Format: Takashi Shimbo, Christopher Lavender, Sara A. Grimm, Makiko I. Doi, Telmo Henriques, Kimberly R. Cannady, Kevin J. Murphy, Daniel A. Gilchrist, Adam Burkholder, Jeffrey J. Hayes, Karen Adelman, Trevor K. Archer, Kenneth S. Zaret, Paul A. Wade. MBD3 regulates chromatin accessibility at active promoters. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2862. doi:10.1158/1538-7445.AM2015-2862

GATA3 in Breast Cancer: Tumor Suppressor or Oncogene?

GATA3 is a highly conserved, essential transcription factor expressed in a number of tissues, including the mammary gland. GATA3 expression is required for normal development of the mammary gland where it is estimated to be the most abundant transcription factor in luminal epithelial cells. In breast cancer, GATA3 expression is highly correlated with the luminal transcriptional program. Recent genomic analysis of human breast cancers has revealed high-frequency mutation in GATA3 in luminal tumors, suggesting "driver" function(s). Here we discuss mutation of GATA3 in breast cancer and the potential mechanism(s) by which mutation may lead to a growth advantage in cancer.

Abstract 405: Mbd3 localizes at promoters, gene bodies and enhancers of active genes

The structure of chromatin is tightly regulated; its disregulation is causal for developmental abnormalities and for cancer. The nucleosome remodeling and deacetylase (NuRD) complex regulates chromatin structure and gene expression via its nucleosome remodeling and histone deacetylase activities. MBD3, a component of NuRD, has pleiotropic roles in embryonic stem (ES) cells and tumor biology. However, the results of recent studies are conflicting, resulting in an unclear model for NuRD function. Here, we investigated the localization of MBD3 in two human breast cancer cell lines (MCF-7 and MDA-MB-231) with two independent genomic technologies (DamID and ChIP-seq). DamID and ChIP-seq data were in an excellent agreement and led consistent conclusions in both cell lines. MBD3 localized not only at promoters and enhancers of active genes but also accumulated at gene bodies; this previously undescribed localization suggests an active role of MBD3 in transcription. MBD3 preferred hypomethylated, CpG rich promoters marked by H3K4me3. We also noted a subset of MBD3 that localized at silent promoters marked by H3K9me3 or H3K27me3, consistent with biochemical studies. Another subset of MBD3 co-localized with H3K27ac and was in physiological proximity with promoters, suggesting a role of MBD3 in the regulation of high-order structure of the genome. Functional analysis of chromatin demonstrated that MBD3 regulates nucleosome occupancy near promoters and in gene bodies. These data suggest multiple roles of MBD3/NuRD in the regulation of chromatin structure and gene expression for both active and silent genes. The findings serve as an important initial step in defining regulatory mechanism of MBD3/NuRD in ES cells and in tumor biology.

Citation Format: Takashi Shimbo, Ying Du, Sara A. Grimm, Archana Dhasarathy, Deepak Mav, Ruchir R. Shah, Huidong Shi, Paul A. Wade. Mbd3 localizes at promoters, gene bodies and enhancers of active genes. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 405. doi:10.1158/1538-7445.AM2014-405

Abstract 475: GATA3 mutations in breast cancer

Breast cancer is one of the most common cancers in women worldwide, resulting in over 400,000 deaths per year. Approximately 90% of breast cancers occur in women who have no family history of breast cancer, indicating that most breast cancers occur following spontaneous mutations resulting from the aging process or the action of environmental agents. Breast cancers can be divided into several subtypes base on gene expression profiles. Luminal breast cancer, the most common subtype, is characterized by expression of estrogen receptor alpha (ER-alpha). The zinc finger transcription factor GATA3, which is known as a regulator of mammary gland development, is frequently expressed in luminal breast cancer cells, and its expression is highly correlated with ER-alpha. In breast cancers, the expression levels of GATA3 directly correlate with favorable prognosis, and its ectopic expression leads to the suppression of tumor metastasis. Comprehensive genomic analysis of breast tumors has revealed frequent somatic mutations of GATA3 in luminal breast cancers. However it is unclear how these GATA3 mutations affect the biological properties of breast cancer. These mutations are located at exclusively in exons 5 and 6 that encode the carboxyl terminus of the protein. Exon 5 contains the second of two zinc finger domains (ZnF2), which is important for the DNA binding activity of GATA3. Exon 6 has no previously described biochemical function.

In this study, we focus on two mutations found in cancer patients - a frame shift at tyrosine 345 (Y345fs) and a frame shift at arginine 330 (R330fs). These two mutants are of particular interest as they produce truncated proteins, and both mutants lack exon 6. Y345fs preserves ZnF2, while R330fs interrupts ZnF2. These two mutants permit dissection of the molecular outcomes flowing from loss of exon 6 (Y345fs) versus the combination of loss of ZnF2 and exon 6 (R330fs). To assess the impact of these mutations, we established breast cancer cells expressing Y345fs and R330fs respectively, and analyzed their phenotype at the cellular and molecular levels. These studies begin to address the biochemical and functional outcomes of commonly occurring mutations in breast cancer.

Citation Format: Motoki Takaku, Aleksandra Adomas, Sara A. Grimm, Shimbo Takashi, Paul A. Wade. GATA3 mutations in breast cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 475. doi:10.1158/1538-7445.AM2014-475

Breast tumor specific mutation in GATA3 affects physiological mechanisms regulating transcription factor turnover

BACKGROUND

The transcription factor GATA3 is a favorable prognostic indicator in estrogen receptor-α (ERα)-positive breast tumors in which it participates with ERα and FOXA1 in a complex transcriptional regulatory program driving tumor growth. GATA3 mutations are frequent in breast cancer and have been classified as driver mutations. To elucidate the contribution(s) of GATA3 alterations to cancer, we studied two breast cancer cell lines, MCF7, which carries a heterozygous frameshift mutation in the second zinc finger of GATA3, and T47D, wild-type at this locus.

METHODS

Immunofluorescence staining and subcellular fractionation were employed to verify cellular localization of GATA3 in T47D and MCF7 cells. To test protein stability, cells were treated with translation inhibitor, cycloheximide or proteasome inhibitor, MG132, and GATA3 abundance was measured over time using immunoblot. GATA3 turn-over in response to hormone was determined by treating the cells with estradiol or ERα agonist, ICI 182,780. DNA binding ability of recombinant GATA3 was evaluated using electrophoretic mobility shift assay and heparin chromatography. Genomic location of GATA3 in MCF7 and T47D cells was assessed by chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq).

RESULTS

GATA3 localized in the nucleus in T47D and MCF7 cells, regardless of the mutation status. The truncated protein in MCF7 had impaired interaction with chromatin and was easily released from the nucleus. Recombinant mutant GATA3 was able to bind DNA to a lesser degree than the wild-type protein. Heterozygosity for the truncating mutation conferred protection from regulated turnover of GATA3, ERα and FOXA1 following estrogen stimulation in MCF7 cells. Thus, mutant GATA3 uncoupled protein-level regulation of master regulatory transcription factors from hormone action. Consistent with increased protein stability, ChIP-seq profiling identified greater genome-wide accumulation of GATA3 in MCF7 cells bearing the mutation, albeit with a similar distribution across the genome, comparing to T47D cells.

CONCLUSIONS

We propose that this specific, cancer-derived mutation in GATA3 deregulates physiologic protein turnover, stabilizes GATA3 binding across the genome and modulates the response of breast cancer cells to estrogen signaling.

Novel DNA motif binding activity observed in vivo with an estrogen receptor α mutant mouse

Estrogen receptor α (ERα) interacts with DNA directly or indirectly via other transcription factors, referred to as "tethering." Evidence for tethering is based on in vitro studies and a widely used "KIKO" mouse model containing mutations that prevent direct estrogen response element DNA- binding. KIKO mice are infertile, due in part to the inability of estradiol (E2) to induce uterine epithelial proliferation. To elucidate the molecular events that prevent KIKO uterine growth, regulation of the pro-proliferative E2 target gene Klf4 and of Klf15, a progesterone (P4) target gene that opposes the pro-proliferative activity of KLF4, was evaluated. Klf4 induction was impaired in KIKO uteri; however, Klf15 was induced by E2 rather than by P4. Whole uterine chromatin immunoprecipitation-sequencing revealed enrichment of KIKO ERα binding to hormone response elements (HREs) motifs. KIKO binding to HRE motifs was verified using reporter gene and DNA-binding assays. Because the KIKO ERα has HRE DNA-binding activity, we evaluated the "EAAE" ERα, which has more severe DNA-binding domain mutations, and demonstrated a lack of estrogen response element or HRE reporter gene induction or DNA-binding. The EAAE mouse has an ERα null-like phenotype, with impaired uterine growth and transcriptional activity. Our findings demonstrate that the KIKO mouse model, which has been used by numerous investigators, cannot be used to establish biological functions for ERα tethering, because KIKO ERα effectively stimulates transcription using HRE motifs. The EAAE-ERα DNA-binding domain mutant mouse demonstrates that ERα DNA-binding is crucial for biological and transcriptional processes in reproductive tissues and that ERα tethering may not contribute to estrogen responsiveness in vivo.

Obesity, rather than diet, drives epigenomic alterations in colonic epithelium resembling cancer progression

While obesity represents one of several risk factors for colorectal cancer in humans, the mechanistic underpinnings of this association remain unresolved. Environmental stimuli, including diet, can alter the epigenetic landscape of DNA cis-regulatory elements affecting gene expression and phenotype. Here, we explored the impact of diet and obesity on gene expression and the enhancer landscape in murine colonic epithelium. Obesity led to the accumulation of histone modifications associated with active enhancers at genomic loci downstream of signaling pathways integral to the initiation and progression of colon cancer. Meanwhile, colon-specific enhancers lost the same histone mark, poising cells for loss of differentiation. These alterations reflect a transcriptional program with many features shared with the program driving colon cancer progression. The interrogation of enhancer alterations by diet in colonic epithelium provides insights into the biology underlying high-fat diet and obesity as risk factors for colon cancer.

Fine-tuning of epigenetic regulation with respect to promoter CpG content in a cell type-specific manner

Epigenetic regulation of gene expression is fundamental for cell type-specific gene expression. However, integrated comparative transcriptomic and epigenomic analyses in various adult primary differentiated cells remain underrepresented. We generated promoter landscapes of DNA methylation and three important histone methylation marks (H3K4me3, H3K9me2, and H3K27me3) in two primary cell types (B lymphocytes and liver) from adult mice. In line with previous studies, we also observed distinct H3K4me3 patterns at promoters dictated by CpG content in differentiated primary cells. We further explored the distribution of initiating RNA polymerase II and elongating RNA polymerase II across genes within different promoter classes, suggesting different rate-limiting steps at CpG-rich vs. CpG-poor genes. Examination of differentially expressed genes revealed that regulation of tissue-specific genes is closely related to gene function regardless of promoter type. Although repressive chromatin marks displayed differential preference to promoters based on CpG content, we observed fine-tuning of the pattern of association of these marks with specific promoter types in a cell type-specific manner. The distribution of H3K9me2 and H3K27me3, relative to CpG content, differed substantially between the two cell types. Cell-type specific accumulation of repressive chromatin marks was also observed at silent genes in both cell types, suggesting that differentiated primary cells may exhibit cell-type specificity in the distribution of repressive chromatin marks. Epigenetic regulation of gene expression and the association of specific histone marks with promoter sequence classes are fine-tuned in a cell type-specific manner. This unexpected finding underscores the value of extensive study of epigenetic marks across cell and tissue types.

MBD3 localizes at promoters, gene bodies and enhancers of active genes

The Mi-2/nucleosome remodeling and histone deacetylase (NuRD) complex is a multiprotein machine proposed to regulate chromatin structure by nucleosome remodeling and histone deacetylation activities. Recent reports describing localization of NuRD provide new insights that question previous models on NuRD action, but are not in complete agreement. Here, we provide location analysis of endogenous MBD3, a component of NuRD complex, in two human breast cancer cell lines (MCF-7 and MDA-MB-231) using two independent genomic techniques: DNA adenine methyltransferase identification (DamID) and ChIP-seq. We observed concordance of the resulting genomic localization, suggesting that these studies are converging on a robust map for NuRD in the cancer cell genome. MBD3 preferentially associated with CpG rich promoters marked by H3K4me3 and showed cell-type specific localization across gene bodies, peaking around the transcription start site. A subset of sites bound by MBD3 was enriched in H3K27ac and was in physical proximity to promoters in three-dimensional space, suggesting function as enhancers. MBD3 enrichment was also noted at promoters modified by H3K27me3. Functional analysis of chromatin indicated that MBD3 regulates nucleosome occupancy near promoters and in gene bodies. These data suggest that MBD3, and by extension the NuRD complex, may have multiple roles in fine tuning expression for both active and silent genes, representing an important step in defining regulatory mechanisms by which NuRD complex controls chromatin structure and modification status.

Chromatin structure is tightly regulated by multiple mechanisms; its dysregulation is associated with developmental abnormalities and disease. The Mi-2/nucleosome remodeling and histone deacetylase (NuRD) complex is proposed to regulate chromatin structure by changing the location and/or the chemical properties of the fundamental building block of chromatin, the nucleosome. NuRD has been shown by genetics to be important for normal development, yet the detailed mechanism of how NuRD regulates chromatin structure is still unclear. Here, we study the localization and function of MBD3, a component of NuRD, in two human breast cancer cell lines using two independent genomic technologies. Our data demonstrate that existing models, which associate NuRD with transcriptional repression, are not completely correct. Rather, MBD3 showed cell-type specific localization at active genes. Moreover, we found a previously unidentified localization of MBD3 across gene bodies and identified a regulatory role for MBD3 in nucleosome organization. Our data provide a reliable starting point from which to address mechanisms by which NuRD controls chromatin structure and nuclear biology.

DNA methylation profiling in human B cells reveals immune regulatory elements and epigenetic plasticity at Alu elements during B-cell activation

Memory is a hallmark of adaptive immunity, wherein lymphocytes mount a superior response to a previously encountered antigen. It has been speculated that epigenetic alterations in memory lymphocytes contribute to their functional distinction from their naive counterparts. However, the nature and extent of epigenetic alterations in memory compartments remain poorly characterized. Here we profile the DNA methylome and the transcriptome of B-lymphocyte subsets representing stages of the humoral immune response before and after antigen exposure in vivo from multiple humans. A significant percentage of activation-induced losses of DNA methylation mapped to transcription factor binding sites. An additional class of demethylated loci mapped to Alu elements across the genome and accompanied repression of DNA methyltransferase 3A. The activation-dependent DNA methylation changes were largely retained in the progeny of activated B cells, generating a similar epigenetic signature in downstream memory B cells and plasma cells with distinct transcriptional programs. These findings provide insights into the methylation dynamics of the genome during cellular differentiation in an immune response.

Research resource: whole-genome estrogen receptor α binding in mouse uterine tissue revealed by ChIP-seq

To advance understanding of mechanisms leading to biological and transcriptional endpoints related to estrogen action in the mouse uterus, we have mapped ERα and RNA polymerase II (PolII) binding sites using chromatin immunoprecipitation followed by sequencing of enriched chromatin fragments. In the absence of hormone, 5184 ERα-binding sites were apparent in the vehicle-treated ovariectomized uterine chromatin, whereas 17,240 were seen 1 h after estradiol (E₂) treatment, indicating that some sites are occupied by unliganded ERα, and that ERα binding is increased by E₂. Approximately 15% of the uterine ERα-binding sites were adjacent to (<10 kb) annotated transcription start sites, and many sites are found within genes or are found more than 100 kb distal from mapped genes; however, the density (sites per base pair) of ERα-binding sites is significantly greater adjacent to promoters. An increase in quantity of sites but no significant positional differences were seen between vehicle and E₂-treated samples in the overall locations of ERα-binding sites either distal from, adjacent to, or within genes. Analysis of the PolII data revealed the presence of poised promoter-proximal PolII on some highly up-regulated genes. Additionally, corecruitment of PolII and ERα to some distal enhancer regions was observed. A de novo motif analysis of sequences in the ERα-bound chromatin confirmed that estrogen response elements were significantly enriched. Interestingly, in areas of ERα binding without predicted estrogen response element motifs, homeodomain transcription factor-binding motifs were significantly enriched. The integration of the ERα- and PolII-binding sites from our uterine sequencing of enriched chromatin fragments data with transcriptional responses revealed in our uterine microarrays has the potential to greatly enhance our understanding of mechanisms governing estrogen response in uterine and other estrogen target tissues.