Large-scale manipulation of promoter DNA methylation reveals context-specific transcriptional responses and stability

BACKGROUND: Cytosine DNA methylation is widely described as a transcriptional repressive mark with the capacity to silence promoters. Epigenome engineering techniques enable direct testing of the effect of induced DNA methylation on endogenous promoters; however, the downstream effects have not yet been comprehensively assessed. RESULTS: Here, we simultaneously induce methylation at thousands of promoters in human cells using an engineered zinc finger-DNMT3A fusion protein, enabling us to test the effect of forced DNA methylation upon transcription, chromatin accessibility, histone modifications, and DNA methylation persistence after the removal of the fusion protein. We find that transcriptional responses to DNA methylation are highly context-specific, including lack of repression, as well as cases of increased gene expression, which appears to be driven by the eviction of methyl-sensitive transcriptional repressors. Furthermore, we find that some regulatory networks can override DNA methylation and that promoter methylation can cause alternative promoter usage. DNA methylation deposited at promoter and distal regulatory regions is rapidly erased after removal of the zinc finger-DNMT3A fusion protein, in a process combining passive and TET-mediated demethylation. Finally, we demonstrate that induced DNA methylation can exist simultaneously on promoter nucleosomes that possess the active histone modification H3K4me3, or DNA bound by the initiated form of RNA polymerase II. CONCLUSIONS: These findings have important implications for epigenome engineering and demonstrate that the response of promoters to DNA methylation is more complex than previously appreciated. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13059-022-02728-5

Genome-wide binding of the orphan nuclear receptor TR4 suggests its general role in fundamental biological processes

<p>Abstract</p> <p>Background</p> <p>The orphan nuclear receptor TR4 (human testicular receptor 4 or NR2C2) plays a pivotal role in a variety of biological and metabolic processes. With no known ligand and few known target genes, the mode of TR4 function was unclear.</p> <p>Results</p> <p>We report the first genome-wide identification and characterization of TR4 <it>in vivo </it>binding. Using chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq), we identified TR4 binding sites in 4 different human cell types and found that the majority of target genes were shared among different cells. TR4 target genes are involved in fundamental biological processes such as RNA metabolism and protein translation. In addition, we found that a subset of TR4 target genes exerts cell-type specific functions. Analysis of the TR4 binding sites revealed that less than 30% of the peaks from any of the cell types contained the DR1 motif previously derived from <it>in vitro </it>studies, suggesting that TR4 may be recruited to the genome via interaction with other proteins. A bioinformatics analysis of the TR4 binding sites predicted a <it>cis </it>regulatory module involving TR4 and ETS transcription factors. To test this prediction, we performed ChIP-seq for the ETS factor ELK4 and found that 30% of TR4 binding sites were also bound by ELK4. Motif analysis of the sites bound by both factors revealed a lack of the DR1 element, suggesting that TR4 binding at a subset of sites is facilitated through the ETS transcription factor ELK4. Further studies will be required to investigate the functional interdependence of these two factors.</p> <p>Conclusions</p> <p>Our data suggest that TR4 plays a pivotal role in fundamental biological processes across different cell types. In addition, the identification of cell type specific TR4 binding sites enables future studies of the pathways underlying TR4 action and its possible role in metabolic diseases.</p

Farnesoid X Receptor Induces Murine Scavenger Receptor Class B Type I via Intron Binding

Farnesoid X receptor (FXR) is a nuclear receptor and a key regulator of liver cholesterol and triglyceride homeostasis. Scavenger receptor class B type I (SR-BI) is critical for reverse cholesterol transport (RCT) by transporting high-density lipoprotein (HDL) into liver. FXR induces SR-BI, however, the underlying molecular mechanism of this induction is not known. The current study confirmed induction of SR-BI mRNA by activated FXR in mouse livers, a human hepatoma cell line, and primary human hepatocytes. Genome-wide FXR binding analysis in mouse livers identified 4 putative FXR response elements in the form of inverse repeat separated by one nucleotide (IR1) at the first intron and 1 IR1 at the downstream of the mouse Sr-bi gene. ChIP-qPCR analysis revealed FXR binding to only the intronic IR1s, but not the downstream one. Luciferase assays and site-directed mutagenesis further showed that 3 out of 4 IR1s were able to activate gene transcription. A 16-week high-fat diet (HFD) feeding in mice increased hepatic Sr-bi gene expression in a FXR-dependent manner. In addition, FXR bound to the 3 bona fide IR1s in vivo, which was increased following HFD feeding. Serum total and HDL cholesterol levels were increased in FXR knockout mice fed the HFD, compared to wild-type mice. In conclusion, the Sr-bi/SR-BI gene is confirmed as a FXR target gene in both mice and humans, and at least in mice, induction of Sr-bi by FXR is via binding to intronic IR1s. This study suggests that FXR may serve as a promising molecular target for increasing reverse cholesterol transport

EBP1 Is a Novel E2F Target Gene Regulated by Transforming Growth Factor-β

Regulation of gene expression requires transcription factor binding to specific DNA elements, and a large body of work has focused on the identification of such sequences. However, it is becoming increasingly clear that eukaryotic transcription factors can exhibit widespread, nonfunctional binding to genomic DNA sites. Conversely, some of these proteins, such as E2F, can also modulate gene expression by binding to non-consensus elements. E2F comprises a family of transcription factors that play key roles in a wide variety of cellular functions, including survival, differentiation, activation during tissue regeneration, metabolism, and proliferation. E2F factors bind to the Erb3-binding protein 1 (EBP1) promoter in live cells. We now show that E2F binding to the EBP1 promoter occurs through two tandem DNA elements that do not conform to typical consensus E2F motifs. Exogenously expressed E2F1 activates EBP1 reporters lacking one, but not both sites, suggesting a degree of redundancy under certain conditions. E2F1 increases the levels of endogenous EBP1 mRNA in breast carcinoma and other transformed cell lines. In contrast, in non-transformed primary epidermal keratinocytes, E2F, together with the retinoblastoma family of proteins, appears to be involved in decreasing EBP1 mRNA abundance in response to growth inhibition by transforming growth factor-β1. Thus, E2F is likely a central coordinator of multiple responses that culminate in regulation of EBP1 gene expression, and which may vary depending on cell type and context

Search for Gravitational Waves from Primordial Black Hole Binary Coalescences in the Galactic Halo

We use data from the second science run of the LIGO gravitational-wave detectors to search for the gravitational waves from primordial black hole (PBH) binary coalescence with component masses in the range 0.2--$1.0 M_\odot$. The analysis requires a signal to be found in the data from both LIGO observatories, according to a set of coincidence criteria. No inspiral signals were found. Assuming a spherical halo with core radius 5 kpc extending to 50 kpc containing non-spinning black holes with masses in the range 0.2--$1.0 M_\odot$, we place an observational upper limit on the rate of PBH coalescence of 63 per year per Milky Way halo (MWH) with 90% confidence.Comment: 7 pages, 4 figures, to be submitted to Phys. Rev.

Clustered ChIP-Seq-defined transcription factor binding sites and histone modifications map distinct classes of regulatory elements

<p>Abstract</p> <p>Background</p> <p>Transcription factor binding to DNA requires both an appropriate binding element and suitably open chromatin, which together help to define regulatory elements within the genome. Current methods of identifying regulatory elements, such as promoters or enhancers, typically rely on sequence conservation, existing gene annotations or specific marks, such as histone modifications and p300 binding methods, each of which has its own biases.</p> <p>Results</p> <p>Herein we show that an approach based on clustering of transcription factor peaks from high-throughput sequencing coupled with chromatin immunoprecipitation (Chip-Seq) can be used to evaluate markers for regulatory elements. We used 67 data sets for 54 unique transcription factors distributed over two cell lines to create regulatory element clusters. By integrating the clusters from our approach with histone modifications and data for open chromatin, we identified general methylation of lysine 4 on histone H3 (H3K4me) as the most specific marker for transcription factor clusters. Clusters mapping to annotated genes showed distinct patterns in cluster composition related to gene expression and histone modifications. Clusters mapping to intergenic regions fall into two groups either directly involved in transcription, including miRNAs and long noncoding RNAs, or facilitating transcription by long-range interactions. The latter clusters were specifically enriched with H3K4me1, but less with acetylation of lysine 27 on histone 3 or p300 binding.</p> <p>Conclusion</p> <p>By integrating genomewide data of transcription factor binding and chromatin structure and using our data-driven approach, we pinpointed the chromatin marks that best explain transcription factor association with different regulatory elements. Our results also indicate that a modest selection of transcription factors may be sufficient to map most regulatory elements in the human genome.</p

Multiple Wnt/ß-Catenin Responsive Enhancers Align with the MYC Promoter through Long-Range Chromatin Loops

Inappropriate activation of c-Myc (MYC) gene expression by the Wnt/ß-catenin signaling pathway is required for colorectal carcinogenesis. The elevated MYC levels in colon cancer cells are attributed in part to ß-catenin/TCF4 transcription complexes that are assembled at proximal Wnt/ß-catenin responsive enhancers (WREs). Recent studies suggest that additional WREs that control MYC expression reside far upstream of the MYC transcription start site. Here, I report the characterization of five novel WREs that localize to a region over 400 kb upstream from MYC. These WREs harbor nucleosomes with post-translational histone modifications that demarcate enhancer and gene promoter regions. Using quantitative chromatin conformation capture, I show that the distal WREs are aligned with the MYC promoter through large chromatin loops. The chromatin loops are not restricted to colon cancer cells, but are also found in kidney epithelial and lung fibroblast cell lines that lack de-regulated Wnt signaling and nuclear ß-catenin/TCF4 complexes. While each chromatin loop is detected in quiescent cells, the positioning of three of the five distal enhancers with the MYC promoter is induced by serum mitogens. These findings suggest that the architecture of the MYC promoter is comprised of distal elements that are juxtaposed through large chromatin loops and that ß-catenin/TCF4 complexes utilize this conformation to activate MYC expression in colon cancer cells

Joint Binding of OTX2 and MYC in Promotor Regions Is Associated with High Gene Expression in Medulloblastoma

Both OTX2 and MYC are important oncogenes in medulloblastoma, the most common malignant brain tumor in childhood. Much is known about MYC binding to promoter regions, but OTX2 binding is hardly investigated. We used ChIP-on-chip data to analyze the binding patterns of both transcription factors in D425 medulloblastoma cells. When combining the data for all promoter regions in the genome, OTX2 binding showed a remarkable bi-modal distribution pattern with peaks around −250 bp upstream and +650 bp downstream of the transcription start sites (TSSs). Indeed, 40.2% of all OTX2-bound TSSs had more than one significant OTX2-binding peak. This OTX2-binding pattern was very different from the TSS-centered single peak binding pattern observed for MYC and other known transcription factors. However, in individual promoter regions, OTX2 and MYC have a strong tendency to bind in proximity of each other. OTX2-binding sequences are depleted near TSSs in the genome, providing an explanation for the observed bi-modal distribution of OTX2 binding. This contrasts to the enrichment of E-box sequences at TSSs. Both OTX2 and MYC binding independently correlated with higher gene expression. Interestingly, genes of promoter regions with multiple OTX2 binding as well as MYC binding showed the highest expression levels in D425 cells and in primary medulloblastomas. Genes within this class of promoter regions were enriched for medulloblastoma and stem cell specific genes. Our data suggest an important functional interaction between OTX2 and MYC in regulating gene expression in medulloblastoma

Genome-Wide Identification of Calcium-Response Factor (CaRF) Binding Sites Predicts a Role in Regulation of Neuronal Signaling Pathways

Calcium-Response Factor (CaRF) was first identified as a transcription factor based on its affinity for a neuronal-selective calcium-response element (CaRE1) in the gene encoding Brain-Derived Neurotrophic Factor (BDNF). However, because CaRF shares no homology with other transcription factors, its properties and gene targets have remained unknown. Here we show that the DNA binding domain of CaRF has been highly conserved across evolution and that CaRF binds DNA directly in a sequence-specific manner in the absence of other eukaryotic cofactors. Using a binding site selection screen we identify a high-affinity consensus CaRF response element (cCaRE) that shares significant homology with the CaRE1 element of Bdnf. In a genome-wide chromatin immunoprecipitation analysis (ChIP-Seq), we identified 176 sites of CaRF-specific binding (peaks) in neuronal genomic DNA. 128 of these peaks are within 10kB of an annotated gene, and 60 are within 1kB of an annotated transcriptional start site. At least 138 of the CaRF peaks contain a common 10-bp motif with strong statistical similarity to the cCaRE, and we provide evidence predicting that CaRF can bind independently to at least 64.5% of these motifs in vitro. Analysis of this set of putative CaRF targets suggests the enrichment of genes that regulate intracellular signaling cascades. Finally we demonstrate that expression of a subset of these target genes is altered in the cortex of Carf knockout (KO) mice. Together these data strongly support the characterization of CaRF as a unique transcription factor and provide the first insight into the program of CaRF-regulated transcription in neurons

Analysis of Xq27-28 linkage in the international consortium for prostate cancer genetics (ICPCG) families.

BACKGROUND: Genetic variants are likely to contribute to a portion of prostate cancer risk. Full elucidation of the genetic etiology of prostate cancer is difficult because of incomplete penetrance and genetic and phenotypic heterogeneity. Current evidence suggests that genetic linkage to prostate cancer has been found on several chromosomes including the X; however, identification of causative genes has been elusive. METHODS: Parametric and non-parametric linkage analyses were performed using 26 microsatellite markers in each of 11 groups of multiple-case prostate cancer families from the International Consortium for Prostate Cancer Genetics (ICPCG). Meta-analyses of the resultant family-specific linkage statistics across the entire 1,323 families and in several predefined subsets were then performed. RESULTS: Meta-analyses of linkage statistics resulted in a maximum parametric heterogeneity lod score (HLOD) of 1.28, and an allele-sharing lod score (LOD) of 2.0 in favor of linkage to Xq27-q28 at 138 cM. In subset analyses, families with average age at onset less than 65 years exhibited a maximum HLOD of 1.8 (at 138 cM) versus a maximum regional HLOD of only 0.32 in families with average age at onset of 65 years or older. Surprisingly, the subset of families with only 2-3 affected men and some evidence of male-to-male transmission of prostate cancer gave the strongest evidence of linkage to the region (HLOD = 3.24, 134 cM). For this subset, the HLOD was slightly increased (HLOD = 3.47 at 134 cM) when families used in the original published report of linkage to Xq27-28 were excluded. CONCLUSIONS: Although there was not strong support for linkage to the Xq27-28 region in the complete set of families, the subset of families with earlier age at onset exhibited more evidence of linkage than families with later onset of disease. A subset of families with 2-3 affected individuals and with some evidence of male to male disease transmission showed stronger linkage signals. Our results suggest that the genetic basis for prostate cancer in our families is much more complex than a single susceptibility locus on the X chromosome, and that future explorations of the Xq27-28 region should focus on the subset of families identified here with the strongest evidence of linkage to this region.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are