Broad H3K4me3 is associated with increased transcription elongation and enhancer activity at tumor-suppressor genes

Decoding the dynamic DNA methylation and hydroxymethylation landscapes in endodermal lineage intermediates during pancreatic differentiation of hESC

Dynamic changes in DNA methylation and demethylation reprogram transcriptional outputs to instruct lineage specification during development. Here, we applied an integrative epigenomic approach to unveil DNA (hydroxy)methylation dynamics representing major endodermal lineage intermediates during pancreatic differentiation of human embryonic stem cells (hESCs). We found that 5-hydroxymethylcytosine (5hmC) marks genomic regions to be demethylated in the descendent lineage, thus reshaping the DNA methylation landscapes during pancreatic lineage progression. DNA hydroxymethylation is positively correlated with enhancer activities and chromatin accessibility, as well as the selective binding of lineage-specific pioneer transcription factors, during pancreatic differentiation. We further discovered enrichment of hydroxymethylated regions (termed '5hmC-rim') at the boundaries of large hypomethylated functional genomic regions, including super-enhancer, DNA methylation canyon and broad-H3K4me3 peaks. We speculate that '5hmC-rim' might safeguard low levels of cytosine methylation at these regions. Our comprehensive analysis highlights the importance of dynamic changes of epigenetic landscapes in driving pancreatic differentiation of hESC.

EpiAlign: an alignment-based bioinformatic tool for comparing chromatin state sequences

The availability of genome-wide epigenomic datasets enables in-depth studies of epigenetic modifications and their relationships with chromatin structures and gene expression. Various alignment tools have been developed to align nucleotide or protein sequences in order to identify structurally similar regions. However, there are currently no alignment methods specifically designed for comparing multi-track epigenomic signals and detecting common patterns that may explain functional or evolutionary similarities. We propose a new local alignment algorithm, EpiAlign, designed to compare chromatin state sequences learned from multi-track epigenomic signals and to identify locally aligned chromatin regions. EpiAlign is a dynamic programming algorithm that novelly incorporates varying lengths and frequencies of chromatin states. We demonstrate the efficacy of EpiAlign through extensive simulations and studies on the real data from the NIH Roadmap Epigenomics project. EpiAlign is able to extract recurrent chromatin state patterns along a single epigenome, and many of these patterns carry cell-type-specific characteristics. EpiAlign can also detect common chromatin state patterns across multiple epigenomes, and it will serve as a useful tool to group and distinguish epigenomic samples based on genome-wide or local chromatin state patterns.

The national blueprint for pregnancy/birth longitudinal cohorts to study factor VIII immunogenicity: NHLBI State of the Science (SOS) Workshop on factor VIII inhibitors

INTRODUCTION

Patients with haemophilia can develop inhibitors to exogenous coagulation factors. Some patients are tolerant to factor, while those who develop inhibitors do so early in life. Genetics and environmental factors are known to contribute to inhibitor risk. However, it is not yet possible to predict inhibitor formation or treatment responsiveness in individuals. We hypothesize that factors in the antenatal/neonatal period inform inhibitor risk development.

AIM

To consider the design of longitudinal studies beginning in the antenatal/neonatal period and the use of new technologies to better understand haemophilia inhibitors.

METHODS

A working group was formed for the NHLBI State of the Science Workshop: Factor VIII Inhibitors: Generating a National Blueprint for Future Research to solicit input from the US haemophilia community and international collaborators to consider design of pregnancy/birth longitudinal cohorts that leverage -omics, existing phenotypic data, and in silico modelling to study inhibitors.

RESULTS

An antenatal/neonatal longitudinal cohort should begin with enrolment of pregnant genetic carriers of haemophilia and span the at-risk period for inhibitor development in the child. Data and samples from the mother, placenta, neonate and young child can be obtained that are amenable to existing assays, genomics and other -omics studies. Data can inform in silico prediction and mathematical models.

CONCLUSION

A longitudinal study beginning before birth offers the unique opportunity to study factors that influence inhibitor development prior to exposure. Advances in -omics and computational biology can study complex phenotypes in this rare disease. This study could be accomplished through interdisciplinary efforts and patient community engagement.

Crosstalk between chromatin structure, cohesin activity and transcription

BACKGROUND

A complex interplay between chromatin and topological machineries is critical for genome architecture and function. However, little is known about these reciprocal interactions, even for cohesin, despite its multiple roles in DNA metabolism.

RESULTS

We have used genome-wide analyses to address how cohesins and chromatin structure impact each other in yeast. Cohesin inactivation in scc1-73 mutants during the S and G2 phases causes specific changes in chromatin structure that preferentially take place at promoters; these changes include a significant increase in the occupancy of the - 1 and + 1 nucleosomes. In addition, cohesins play a major role in transcription regulation that is associated with specific promoter chromatin architecture. In scc1-73 cells, downregulated genes are enriched in promoters with short or no nucleosome-free region (NFR) and a fragile "nucleosome - 1/RSC complex" particle. These results, together with a preferential increase in the occupancy of nucleosome - 1 of these genes, suggest that cohesins promote transcription activation by helping RSC to form the NFR. In sharp contrast, the scc1-73 upregulated genes are enriched in promoters with an "open" chromatin structure and are mostly at cohesin-enriched regions, suggesting that a local accumulation of cohesins might help to inhibit transcription. On the other hand, a dramatic loss of chromatin integrity by histone depletion during DNA replication has a moderate effect on the accumulation and distribution of cohesin peaks along the genome.

CONCLUSIONS

Our analyses of the interplay between chromatin integrity and cohesin activity suggest that cohesins play a major role in transcription regulation, which is associated with specific chromatin architecture and cohesin-mediated nucleosome alterations of the regulated promoters. In contrast, chromatin integrity plays only a minor role in the binding and distribution of cohesins.

NUCLIZE for quantifying epigenome: generating histone modification data at single-nucleosome resolution using genuine nucleosome positions

Background

Defining histone modification at single-nucleosome resolution provides accurate epigenomic information in individual nucleosomes. However, most of histone modification data deposited in current databases, such as ENCODE and Roadmap, have low resolution with peaks of several kilo-base pairs (kb), which due to the technical defects of regular ChIP-Seq technology.

Results

To generate histone modification data at single-nucleosome resolution, we developed a novel approach, NUCLIZE, using synergistic analyses of histone modification data from ChIP-Seq and high-resolution nucleosome mapping data from native MNase-Seq. With this approach, we generated quantitative epigenomics data of single and multivalent histone modification marks in each nucleosome. We found that the dominant trivalent histone mark (H3K4me3/H3K9ac/H3K27ac) and others showed defined and specific patterns near each TSS, indicating potential epigenetic codes regulating gene transcription.

Conclusions

Single-nucleosome histone modification data render epigenomic data become quantitative, which is essential for investigating dynamic changes of epigenetic regulation in the biological process or for functional epigenomics studies. Thus, NUCLIZE turns current epigenomic mapping studies into genuine functional epigenomics studies with quantitative epigenomic data.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5932-6) contains supplementary material, which is available to authorized users.

Mediator complex interaction partners organize the transcriptional network that defines neural stem cells

The Mediator complex regulates transcription by connecting enhancers to promoters. High Mediator binding density defines super enhancers, which regulate cell-identity genes and oncogenes. Protein interactions of Mediator may explain its role in these processes but have not been identified comprehensively. Here, we purify Mediator from neural stem cells (NSCs) and identify 75 protein-protein interaction partners. We identify super enhancers in NSCs and show that Mediator-interacting chromatin modifiers colocalize with Mediator at enhancers and super enhancers. Transcription factor families with high affinity for Mediator dominate enhancers and super enhancers and can explain genome-wide Mediator localization. We identify E-box transcription factor Tcf4 as a key regulator of NSCs. Tcf4 interacts with Mediator, colocalizes with Mediator at super enhancers and regulates neurogenic transcription factor genes with super enhancers and broad H3K4me3 domains. Our data suggest that high binding-affinity for Mediator is an important organizing feature in the transcriptional network that determines NSC identity.

The Mediator complex regulates transcription by connecting enhancers to promoters. Here, the authors purify Mediator from neural stem cells (NSCs), identify 75 novel protein-protein interaction partners and characterize the Mediator-interacting network that regulates transcription and establishes NSC identity.

LncRNA-p21 alters the antiandrogen enzalutamide-induced prostate cancer neuroendocrine differentiation via modulating the EZH2/STAT3 signaling

While the antiandrogen enzalutamide (Enz) extends the castration resistant prostate cancer (CRPC) patients' survival an extra 4.8 months, it might also result in some adverse effects via inducing the neuroendocrine differentiation (NED). Here we found that lncRNA-p21 is highly expressed in the NEPC patients derived xenograft tissues (NEPC-PDX). Results from cell lines and human clinical sample surveys also revealed that lncRNA-p21 expression is up-regulated in NEPC and Enz treatment could increase the lncRNA-p21 to induce the NED. Mechanism dissection revealed that Enz could promote the lncRNA-p21 transcription via altering the androgen receptor (AR) binding to different androgen-response-elements, which switch the EZH2 function from histone-methyltransferase to non-histone methyltransferase, consequently methylating the STAT3 to promote the NED. Preclinical studies using the PDX mouse model proved that EZH2 inhibitor could block the Enz-induced NED. Together, these results suggest targeting the Enz/AR/lncRNA-p21/EZH2/STAT3 signaling may help urologists to develop a treatment for better suppression of the human CRPC progression.

Disruption of Broad Epigenetic Domains in PDAC Cells by HAT Inhibitors

The spreading of epigenetic domains has emerged as a distinguishing epigenomic phenotype for diverse cell types. In particular, clusters of H3K27ac- and H3K4me3-marked elements, referred to as super-enhancers, and broad H3K4me3 domains, respectively, have been linked to cell identity and disease states. Here, we characterized the broad domains from different pancreatic ductal adenocarcinoma (PDAC) cell lines that represent distinct histological grades. Our integrative genomic analysis found that human derived cell line models for distinct PDAC grades exhibit characteristic broad epigenetic features associated with gene expression patterns that are predictive of patient prognosis and provide insight into pancreatic cancer cell identity. In particular, we find that genes marked by overlapping Low-Grade broad domains correspond to an epithelial phenotype and hold potential as markers for patient stratification. We further utilize ChIP-seq to compare the effects of histone acetyltransferase (HAT) inhibitors to detect global changes in histone acetylation and methylation levels. We found that HAT inhibitors impact certain broad domains of pancreatic cancer cells. Overall, our results reveal potential roles for broad domains in cells from distinct PDAC grades and demonstrate the plasticity of particular broad epigenomic domains to epigenetic inhibitors.

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

Aberrant activation of CYR61 enhancers in colorectal cancer development

BACKGROUND

High expression of secreted matricellular protein cysteine-rich 61 (CYR61) correlates with poor prognosis in colorectal cancer (CRC). Aberrant enhancer activation has been shown to correlate with expression of key genes involved in cancer progression. However, such mechanisms in CYR61 transcription regulation remain unexplored.

METHODS

Expression of CYR61 was determined by immunohistochemistry (IHC), quantitative real-time PCR (qRT-PCR) and western blotting (WB) in CRC patients paraffin specimens and colon cell lines. ChIP-seq data of enhancer-characteristic histone modifications, in CRC tissues from the Gene Expression Omnibus (GEO) database, were reanalyzed to search for putative enhancers of CYR61. Dual-luciferase reporter assay was used to detected enhancer activity. Physical interactions between putative enhancers and CYR61 promoter were detected by chromosome conformation capture (3C) assay. Histone modification and transcription factors (TFs) enrichment were detected by ChIP-qPCR. Additionally, biological function of enhancers was investigated by transwell migration assays.

RESULTS

CRC tissues and cell lines expressed higher level of CYR61 than normal colon mucosa. Three putative enhancers located downstream of CYR61 were found in CRC tissues by ChIP-seq data reanalysis. Consistent with the ChIP-seq analysis results in the GEO database, the normal colon mucosal epithelial cell line NCM460 possessed no active CYR61 enhancers, whereas colon cancer cells exhibited different patterns of active CYR61 enhancers. HCT116 cells had an active Enhancer3, whereas RKO cells had both Enhancer1 and Enhancer3 active. Pioneer factor FOXA1 promoted CYR61 expression by recruiting CBP histone acetyltransferase binding and increasing promoter-enhancer looping frequencies and enhancer activity. CBP knockdown attenuated H3K27ac enrichment, promoter-enhancer looping frequencies, and enhancer activity. Small molecule compound 12-O-tetradecanoyl phorbol-13-acetate (TPA) treatment, which stimulated CYR61 expression, and verteporfin (VP) treatment, which inhibited CYR61 expression, confirmed that the enhancers regulated CYR61 expression. Knockdown and ectopic expression of CYR61 rescued cell migration changes induced by over-expressing and knockdown of FOXA1, respectively.

CONCLUSIONS

CYR61 enhancer activation, mediated by FOXA1 and CBP, occurs during CRC progression to up-regulate CYR61 expression and promote cell migration in CRC, suggesting inhibition of recruitment of FOXA1 and/or CBP to CYR61 enhancers may have therapeutic implications.

A Role for Chromatin Remodeling in Cohesin Loading onto Chromosomes

Cohesin is a conserved, ring-shaped protein complex that topologically embraces DNA. Its central role in genome organization includes functions in sister chromatid cohesion, DNA repair, and transcriptional regulation. Cohesin loading onto chromosomes requires the Scc2-Scc4 cohesin loader, whose presence on chromatin in budding yeast depends on the RSC chromatin remodeling complex. Here we reveal a dual role of RSC in cohesin loading. RSC acts as a chromatin receptor that recruits Scc2-Scc4 by a direct protein interaction independent of chromatin remodeling. In addition, chromatin remodeling is required to generate a nucleosome-free region that is the substrate for cohesin loading. An engineered cohesin loading module can be created by fusing the Scc2 C terminus to RSC or to other chromatin remodelers, but not to unrelated DNA binding proteins. These observations demonstrate the importance of nucleosome-free DNA for cohesin loading and provide insight into how cohesin accesses DNA during its varied chromosomal activities.

Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

The liver plays a key role in metabolism and affects pig production. However, the functional annotation of noncoding regions of the pig liver remains poorly understood. We revealed the landscape of cis-regulatory elements and their functional characterization in pig liver. We identified 102,373 cis-regulatory elements in the pig liver, including enhancers, promoters, super-enhancers, and broad H3K4me3 domains, and highlighted 26 core transcription regulatory factors in the pig liver as well. We found similarity of cis-regulatory elements among those of pigs, humans, and cattle. Despite the low proportion of functionally conserved enhancers (~30%) between pig and human liver tissue, ~78% of the pig liver enhancer orthologues sequence could play an enhancer role in other human tissues. Additionally, we observed that the ratio of consistent super-enhancer-associated genes was significantly higher than the ratio of functionally conserved super-enhancers. Approximately 54% of the core regulation factors driven by super-enhancers were consistent across the liver from these three species. Our pig liver annotation and functional characterization studies provide a system and resource for noncoding annotation for future gene regulatory studies in pigs. Furthermore, our study also showed the high level functional conservation of cis-regulatory elements in mammals; it also improved our understanding of regulation function of mammal cis-regulatory elements.

Broad domains of histone 3 lysine 4 trimethylation are associated with transcriptional activation in CA1 neurons of the hippocampus during memory formation

Transcriptional changes in the hippocampus are required for memory formation, and these changes are regulated by numerous post-translational modifications of chromatin-associated proteins. One of the epigenetic marks that has been implicated in memory formation is histone 3 lysine 4 trimethylation (H3K4me3), and this modification is found at the promoters of actively transcribed genes. The total levels of H3K4me3 are increased in the CA1 region of the hippocampus during memory formation, and genetic perturbation of the K4 methyltransferases and demethylases interferes with forming memories. Previous chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) analyses failed to detect changes in H3K4me3 levels at the promoters of memory-linked genes. Since the breadth of H3K4me3 marks was recently reported to be associated with the transcriptional outcome of a gene, we re-analyzed H3K4me3 ChIP-seq data sets to identify the role of H3K4me3 broad domains in CA1 neurons, as well as identify differences in breadth that occur during contextual fear conditioning. We found that, under baseline conditions, broad H3K4me3 peaks mark important learning and memory genes and are often regulated by super-enhancers. The peaks at many learning-associated genes become broader during novel environment exposure and memory formation. Furthermore, the important learning- and memory-associated lysine methyltransferases, Kmt2a and Kmt2b, are involved in maintaining H3K4me3 peak width. Our findings highlight the importance of analyzing H3K4me3 peak shape, and demonstrate that breadth of H3K4me3 marks in neurons of the hippocampus is regulated during memory formation.

Targeting Super-Enhancers as a Therapeutic Strategy for Cancer Treatment

Super-enhancers (SEs) refer to large clusters of enhancers that drive gene expressions. Recent data has provided novel insights in elucidating the roles of SEs in many diseases, including cancer. Many mechanisms involved in tumorigenesis and progression, ranging from internal gene mutation and rearrangement to external damage and inducement, have been demonstrated to be highly associated with SEs. Moreover, translocation, formation, deletion, or duplication of SEs themselves could lead to tumor development. It has been reported that various oncogenic molecules and pathways are tightly regulated by SEs. Moreover, several clinical trials on novel SEs blockers, such as BET inhibitor and CDK7i, have indicated the potential roles of SEs in cancer therapy. In this review, we highlighted the underlying mechanism of action of SEs in cancer development and the corresponding novel potential therapeutic strategies. It is speculated that targeting SEs could complement the traditional approaches and lead to more effective treatment for cancer patients.

Genome-Wide Identification of Histone Modifications Involved in Placental Development in Pigs

Development of placental folds is a critical event affecting placental function in pigs because it can increase surface area for improvement in capillary density as gestation advances. However, the molecular mechanisms of the event are not well defined. Histone modifications have important roles in gene regulation. To investigate their effects on regulation of genes controlling porcine placental development, RNA-seq and ChIP-seq of porcine placental tissues from gestational days 50 (establishment stage of placental folds) and 95 (expanding stage of placental folds) were carried out in this study. The differentially expressed genes were identified and of which the down- and up-regulated genes are related to endoplasmic reticulum (ER) stress and angiogenesis, respectively. In addition, we mapped the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), which are associated with transcriptional activation. A number of differential modification regions between the 2 gestational stages were identified and majority of them are those with increased signals of H3K4me3 (14,576 out of 16,931). Furthermore, we observed that the increase of H3K4me3 is significantly correlated with the elevated expression levels of the neighboring genes, and notably, these genes were enriched in pathways related to blood vessel formation and microvascular permeability. Taken together, the findings suggest important roles of histone modifications on placental remolding in response to developmental changes.

Drosophila Histone Demethylase KDM5 Regulates Social Behavior through Immune Control and Gut Microbiota Maintenance

Loss-of-function mutations in the histone demethylases KDM5A, KDM5B, or KDM5C are found in intellectual disability (ID) and autism spectrum disorders (ASD) patients. Here, we use the model organism Drosophila melanogaster to delineate how KDM5 contributes to ID and ASD. We show that reducing KDM5 causes intestinal barrier dysfunction and changes in social behavior that correlates with compositional changes in the gut microbiota. Therapeutic alteration of the dysbiotic microbiota through antibiotic administration or feeding with a probiotic Lactobacillus strain partially rescues the behavioral, lifespan, and cellular phenotypes observed in kdm5-deficient flies. Mechanistically, KDM5 was found to transcriptionally regulate component genes of the immune deficiency (IMD) signaling pathway and subsequent maintenance of host-commensal bacteria homeostasis in a demethylase-dependent manner. Together, our study uses a genetic approach to dissect the role of KDM5 in the gut-microbiome-brain axis and suggests that modifying the gut microbiome may provide therapeutic benefits for ID and ASD patients.

Remodeling of epigenome and transcriptome landscapes with aging in mice reveals widespread induction of inflammatory responses

Aging is accompanied by the functional decline of tissues. However, a systematic study of epigenomic and transcriptomic changes across tissues during aging is missing. Here, we generated chromatin maps and transcriptomes from four tissues and one cell type from young, middle-aged, and old mice—yielding 143 high-quality data sets. We focused on chromatin marks linked to gene expression regulation and cell identity: histone H3 trimethylation at lysine 4 (H3K4me3), a mark enriched at promoters, and histone H3 acetylation at lysine 27 (H3K27ac), a mark enriched at active enhancers. Epigenomic and transcriptomic landscapes could easily distinguish between ages, and machine-learning analysis showed that specific epigenomic states could predict transcriptional changes during aging. Analysis of data sets from all tissues identified recurrent age-related chromatin and transcriptional changes in key processes, including the up-regulation of immune system response pathways such as the interferon response. The up-regulation of the interferon response pathway with age was accompanied by increased transcription and chromatin remodeling at specific endogenous retroviral sequences. Pathways misregulated during mouse aging across tissues, notably innate immune pathways, were also misregulated with aging in other vertebrate species—African turquoise killifish, rat, and humans—indicating common signatures of age across species. To date, our data set represents the largest multitissue epigenomic and transcriptomic data set for vertebrate aging. This resource identifies chromatin and transcriptional states that are characteristic of young tissues, which could be leveraged to restore aspects of youthful functionality to old tissues.

The chromatin remodeling protein Lsh alters nucleosome occupancy at putative enhancers and modulates binding of lineage specific transcription factors

Dynamic regulation of chromatin accessibility is a key feature of cellular differentiation during embryogenesis, but the precise factors that control access to chromatin remain largely unknown. Lsh/HELLS is critical for normal development and mutations of Lsh in human cause the ICF (Immune deficiency, Centromeric instability, Facial anomalies) syndrome, a severe immune disorder with multiple organ deficiencies. We report here that Lsh, previously known to regulate DNA methylation level, has a genome wide chromatin remodeling function. Using micrococcal nuclease (MNase)-seq analysis, we demonstrate that Lsh protects MNase accessibility at transcriptional regulatory regions characterized by DNase I hypersensitivity and certain histone 3 (H3) tail modifications associated with enhancers. Using an auxin-inducible degron system, allowing proteolytical degradation of Lsh, we show that Lsh mediated changes in nucleosome occupancy are independent of DNA methylation level and are characterized by reduced H3 occupancy. While Lsh mediated nucleosome occupancy prevents binding sites for transcription factors in wild type cells, depletion of Lsh leads to an increase in binding of ectopically expressed tissue specific transcription factors to their respective binding sites. Our data suggests that Lsh mediated chromatin remodeling can modulate nucleosome positioning at a subset of putative enhancers contributing to the preservation of cellular identity through regulation of accessibility.

Characterizing the nuclease accessibility of DNA in human cells to map higher order structures of chromatin

Packaging of DNA into chromatin regulates DNA accessibility and consequently all DNA-dependent processes. The nucleosome is the basic packaging unit of DNA forming arrays that are suggested, by biochemical studies, to fold hierarchically into ordered higher-order structures of chromatin. This organization has been recently questioned using microscopy techniques, proposing an irregular structure. To address the principles of chromatin organization, we applied an in situ differential MNase-seq strategy and analyzed in silico the results of complete and partial digestions of human chromatin. We investigated whether different levels of chromatin packaging exist in the cell. We assessed the accessibility of chromatin within distinct domains of kb to Mb genomic regions, performed statistical analyses and computer modelling. We found no difference in MNase accessibility, suggesting no difference in fiber folding between domains of euchromatin and heterochromatin or between other sequence and epigenomic features of chromatin. Thus, our data suggests the absence of differentially organized domains of higher-order structures of chromatin. Moreover, we identified only local structural changes, with individual hyper-accessible nucleosomes surrounding regulatory elements, such as enhancers and transcription start sites. The regulatory sites per se are occupied with structurally altered nucleosomes, exhibiting increased MNase sensitivity. Our findings provide biochemical evidence that supports an irregular model of large-scale chromatin organization.

Cell Specificity of Human Regulatory Annotations and Their Genetic Effects on Gene Expression

Epigenomic signatures from histone marks and transcription factor (TF)-binding sites have been used to annotate putative gene regulatory regions. However, a direct comparison of these diverse annotations is missing, and it is unclear how genetic variation within these annotations affects gene expression. Here, we compare five widely used annotations of active regulatory elements that represent high densities of one or more relevant epigenomic marks-"super" and "typical" (nonsuper) enhancers, stretch enhancers, high-occupancy target (HOT) regions, and broad domains-across the four matched human cell types for which they are available. We observe that stretch and super enhancers cover cell type-specific enhancer "chromatin states," whereas HOT regions and broad domains comprise more ubiquitous promoter states. Expression quantitative trait loci (eQTL) in stretch enhancers have significantly smaller effect sizes compared to those in HOT regions. Strikingly, chromatin accessibility QTL in stretch enhancers have significantly larger effect sizes compared to those in HOT regions. These observations suggest that stretch enhancers could harbor genetically primed chromatin to enable changes in TF binding, possibly to drive cell type-specific responses to environmental stimuli. Our results suggest that current eQTL studies are relatively underpowered or could lack the appropriate environmental context to detect genetic effects in the most cell type-specific "regulatory annotations," which likely contributes to infrequent colocalization of eQTL with genome-wide association study signals.

Guide Positioning Sequencing identifies aberrant DNA methylation patterns that alter cell identity and tumor-immune surveillance networks

Aberrant DNA methylation is a distinguishing feature of cancer. Yet, how methylation affects immune surveillance and tumor metastasis remains ambiguous. We introduce a novel method, Guide Positioning Sequencing (GPS), for precisely detecting whole-genome DNA methylation with cytosine coverage as high as 96% and unbiased coverage of GC-rich and repetitive regions. Systematic comparisons of GPS with whole-genome bisulfite sequencing (WGBS) found that methylation difference between gene body and promoter is an effective predictor of gene expression with a correlation coefficient of 0.67 (GPS) versus 0.33 (WGBS). Moreover, Methylation Boundary Shift (MBS) in promoters or enhancers is capable of modulating expression of genes associated with immunity and tumor metabolism. Furthermore, aberrant DNA methylation results in tissue-specific enhancer switching, which is responsible for altering cell identity during liver cancer development. Altogether, we demonstrate that GPS is a powerful tool with improved accuracy and efficiency over WGBS in simultaneously detecting genome-wide DNA methylation and genomic variation. Using GPS, we show that aberrant DNA methylation is associated with altering cell identity and immune surveillance networks, which may contribute to tumorigenesis and metastasis.

Histone H3 lysine K4 methylation and its role in learning and memory

Epigenetic modifications such as histone methylation permit change in chromatin structure without accompanying change in the underlying genomic sequence. A number of studies in animal models have shown that dysregulation of various components of the epigenetic machinery causes cognitive deficits at the behavioral level, suggesting that proper epigenetic control is necessary for the fundamental processes of learning and memory. Histone H3 lysine K4 (H3K4) methylation comprises one component of such epigenetic control, and global levels of this mark are increased in the hippocampus during memory formation. Modifiers of H3K4 methylation are needed for memory formation, shown through animal studies, and many of the same modifiers are mutated in human cognitive diseases. Indeed, all of the known H3K4 methyltransferases and four of the known six H3K4 demethylases have been associated with impaired cognition in a neurologic or psychiatric disorder. Cognitive impairment in such patients often manifests as intellectual disability, consistent with a role for H3K4 methylation in learning and memory. As a modification quintessentially, but not exclusively, associated with transcriptional activity, H3K4 methylation provides unique insights into the regulatory complexity of writing, reading, and erasing chromatin marks within an activated neuron. The following review will discuss H3K4 methylation and connect it to transcriptional events required for learning and memory within the developed nervous system. This will include an initial discussion of the most recent advances in the developing methodology to analyze H3K4 methylation, namely mass spectrometry and deep sequencing, as well as how these methods can be applied to more deeply understand the biology of this mark in the brain. We will then introduce the core enzymatic machinery mediating addition and removal of H3K4 methylation marks and the resulting epigenetic signatures of these marks throughout the neuronal genome. We next foray into the brain, discussing changes in H3K4 methylation marks within the hippocampus during memory formation and retrieval, as well as the behavioral correlates of H3K4 methyltransferase deficiency in this region. Finally, we discuss the human cognitive diseases connected to each H3K4 methylation modulator and summarize advances in developing drugs to target them.

Chorioamnionitis exposure remodels the unique histone modification landscape of neonatal monocytes and alters the expression of immune pathway genes

Chorioamnionitis is an intrauterine infection involving inflammation of the chorion, amnion, and placenta. It leads to a fetal systemic inflammatory response that can alter the transcription of neonatal immune genes. We have previously shown that neonatal monocytes gain the activating histone tail modification H3K4me3 at promoter sites of immunologically important genes as development progresses from preterm neonate to adult. In this study, we applied ChIP-seq and RNA-seq to evaluate the impact of chorioamnionitis on the neonatal monocyte H3K4me3 histone modification landscape over the course of fetal and neonatal immune system development. Chorioamnionitis exposure in neonatal monocytes resulted in a net increase in total monocyte H3K4me3, primarily in introns and intergenic regions. Immune gene expression was decreased in chorioamnionitis-exposed monocytes, with the majority of enriched transcripts falling into pathways that are not linked to the immune system. Over half of all neonatal monocyte H3K4me3 peaks, independent of their location, were associated with active gene transcription. Overall, chorioamnionitis exposure resulted in the global remodeling of the neonatal monocyte H3K4me3 landscape and changes in the expression of known immune genes. These changes resulted in a less robust inflammatory response upon exposure to a secondary challenge, which may explain why chorioamnionitis-exposed neonates have an increased risk of sepsis. DATABASE: ChIP-seq data for U30/O30/Term: GEO GSE81957 ChIP-seq data for U30C/O30C/TermC: GEO GSE111873 RNA-seq data for U/L/CU/CL: GEO GSE111927.

Drug targeting of NR4A nuclear receptors for treatment of acute myeloid leukemia

NR4As are AML tumor suppressors that are frequently silenced in human acute myeloid leukemia (AML). Despite their potential as novel targets for therapeutic intervention, mechanisms of NR4A silencing and strategies for their reactivation remain poorly defined. Here we show that NR4A silencing in AML occurs through blockade of transcriptional elongation rather than epigenetic promoter silencing. By intersection of NR4A-regulated gene signatures captured upon acute, exogenous expression of NR4As in human AML cells with in silico chemical genomics screening, we identify several FDA-approved drugs including dihydroergotamine (DHE) that reactivate NR4A expression and regulate NR4A-dependent gene signatures. We show that DHE induces NR4A expression via recruitment of the super elongation complex to enable elongation of NR4A promoter paused RNA polymerase II. Finally, DHE exhibits AML selective NR4A-dependent anti-leukemic activity in cytogenetically distinct human AML cells in vitro and delays AML progression in mice revealing its potential as a novel therapeutic agent in AML.

CRlncRC: a machine learning-based method for cancer-related long noncoding RNA identification using integrated features

BACKGROUND

Long noncoding RNAs (lncRNAs) are widely involved in the initiation and development of cancer. Although some computational methods have been proposed to identify cancer-related lncRNAs, there is still a demanding to improve the prediction accuracy and efficiency. In addition, the quick-update data of cancer, as well as the discovery of new mechanism, also underlay the possibility of improvement of cancer-related lncRNA prediction algorithm. In this study, we introduced CRlncRC, a novel Cancer-Related lncRNA Classifier by integrating manifold features with five machine-learning techniques.

RESULTS

CRlncRC was built on the integration of genomic, expression, epigenetic and network, totally in four categories of features. Five learning techniques were exploited to develop the effective classification model including Random Forest (RF), Naïve bayes (NB), Support Vector Machine (SVM), Logistic Regression (LR) and K-Nearest Neighbors (KNN). Using ten-fold cross-validation, we showed that RF is the best model for classifying cancer-related lncRNAs (AUC = 0.82). The feature importance analysis indicated that epigenetic and network features play key roles in the classification. In addition, compared with other existing classifiers, CRlncRC exhibited a better performance both in sensitivity and specificity. We further applied CRlncRC to lncRNAs from the TANRIC (The Atlas of non-coding RNA in Cancer) dataset, and identified 121 cancer-related lncRNA candidates. These potential cancer-related lncRNAs showed a certain kind of cancer-related indications, and many of them could find convincing literature supports.

CONCLUSIONS

Our results indicate that CRlncRC is a powerful method for identifying cancer-related lncRNAs. Machine-learning-based integration of multiple features, especially epigenetic and network features, had a great contribution to the cancer-related lncRNA prediction. RF outperforms other learning techniques on measurement of model sensitivity and specificity. In addition, using CRlncRC method, we predicted a set of cancer-related lncRNAs, all of which displayed a strong relevance to cancer as a valuable conception for the further cancer-related lncRNA function studies.

Single-molecule nascent RNA sequencing identifies regulatory domain architecture at promoters and enhancers

Eukaryotic RNA polymerase II (Pol II) has been found at both promoters and distal enhancers, suggesting additional functions beyond mRNA production. To understand this role, we sequenced nascent RNAs at single-molecule resolution to unravel the interplay between Pol II initiation, capping and pausing genome-wide. Our analyses identify two pause classes that are associated with different RNA capping profiles. More proximal pausing is associated with less complete capping, less elongation and a more enhancer-like complement of transcription factors than later pausing. Unexpectedly, transcription start sites (TSSs) are predominantly found in constellations composed of multiple divergent pairs. TSS clusters are intimately associated with precise arrays of nucleosomes and correspond with boundaries of transcription factor binding and chromatin modification at promoters and enhancers. TSS architecture is largely unchanged during the dramatic transcriptional changes induced by heat shock. Together, our results suggest that promoter- and enhancer-associated Pol II is a regulatory nexus for integrating information across TSS ensembles.

Super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers

Super-enhancers and stretch enhancers represent classes of transcriptional enhancers that have been shown to control the expression of cell identity genes and carry disease- and trait-associated variants. Specifically, super-enhancers are clusters of enhancers defined based on the binding occupancy of master transcription factors, chromatin regulators, or chromatin marks, while stretch enhancers are large chromatin-defined regulatory regions of at least 3,000 base pairs. Several studies have characterized these regulatory regions in numerous cell types and tissues to decipher their functional importance. However, the differences and similarities between these regulatory regions have not been fully assessed. We integrated genomic, epigenomic, and transcriptomic data from ten human cell types to perform a comparative analysis of super and stretch enhancers with respect to their chromatin profiles, cell type-specificity, and ability to control gene expression. We found that stretch enhancers are more abundant, more distal to transcription start sites, cover twice as much the genome, and are significantly less conserved than super-enhancers. In contrast, super-enhancers are significantly more enriched for active chromatin marks and cohesin complex, and more transcriptionally active than stretch enhancers. Importantly, a vast majority of super-enhancers (85%) overlap with only a small subset of stretch enhancers (13%), which are enriched for cell type-specific biological functions, and control cell identity genes. These results suggest that super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers, and importantly, most of the stretch enhancers that are distinct from super-enhancers do not show an association with cell identity genes, are less active, and more likely to be poised enhancers.

Examining the Roles of H3K4 Methylation States with Systematically Characterized Antibodies

Histone post-translational modifications (PTMs) are important genomic regulators often studied by chromatin immunoprecipitation (ChIP), whereby their locations and relative abundance are inferred by antibody capture of nucleosomes and associated DNA. However, the specificity of antibodies within these experiments has not been systematically studied. Here, we use histone peptide arrays and internally calibrated ChIP (ICeChIP) to characterize 52 commercial antibodies purported to distinguish the H3K4 methylforms (me1, me2, and me3, with each ascribed distinct biological functions). We find that many widely used antibodies poorly distinguish the methylforms and that high- and low-specificity reagents can yield dramatically different biological interpretations, resulting in substantial divergence from the literature for numerous H3K4 methylform paradigms. Using ICeChIP, we also discern quantitative relationships between enhancer H3K4 methylation and promoter transcriptional output and can measure global PTM abundance changes. Our results illustrate how poor antibody specificity contributes to the "reproducibility crisis," demonstrating the need for rigorous, platform-appropriate validation.

Mutational analysis of uterine cervical cancer that survived multiple rounds of radiotherapy

Radiotherapy is an essential component of cancer therapy. Despite advances in cancer genomics, the mutation signatures of radioresistant tumors have not yet been fully elucidated. To address this issue, we analyzed a unique set of clinical specimens from a uterine cervical cancer that repeatedly locally recurred after multiple rounds of radiotherapy. Exon sequencing of 409 cancer-related genes in the treatment-naïve tumor and the tumors that recurred after initial and secondary radiotherapy identified (i) activating mutations in PIK3CA and KRAS, and putative inactivating mutations in SMAD4, as trunk mutation signatures that persisted over the clinical course; and (ii) mutations in KMT2A, TET1, and NLRP1 as acquired mutation signatures observed only in recurrent tumors after radiotherapy. Comprehensive mining of published in vitro genomics data pertaining to radiosensitivity revealed that simultaneous mutations in KRAS and SMAD4, which have not been described previously in uterine cervical cancer, are associated with cancer cell radioresistance. The association between this mutation signature and radioresistance was validated by isogenic cell-based experiments. These results provide proof-of-principle for the analytical pipeline employed in this study, which explores clinically relevant mutation signatures for radioresistance, and demonstrate that this approach is worth pursuing with larger cohorts in the future.

ChromTime: modeling spatio-temporal dynamics of chromatin marks

To model spatial changes of chromatin mark peaks over time we develop and apply ChromTime, a computational method that predicts peaks to be either expanding, contracting, or holding steady between time points. Predicted expanding and contracting peaks can mark regulatory regions associated with transcription factor binding and gene expression changes. Spatial dynamics of peaks provide information about gene expression changes beyond localized signal density changes. ChromTime detects asymmetric expansions and contractions, which for some marks associate with the direction of transcription. ChromTime facilitates the analysis of time course chromatin data in a range of biological systems.

Germline Duplication of SNORA18L5 Increases Risk for HBV-related Hepatocellular Carcinoma by Altering Localization of Ribosomal Proteins and Decreasing Levels of p53

BACKGROUND & AIMS

Single nucleotide polymorphisms could affect risk for hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). We performed a germline copy number variation (CNV)-based genome-wide association study (GWAS) in populations of Chinese ancestry to search for germline CNVs that increase risk of HCC.

METHODS

We conducted a CNV-based GWAS of 1583 HCC cases (persons with chronic HBV infection and HCC) and 1540 controls (persons with chronic HBV infection without HCC) in Chinese populations. Identified candidates were expressed in L-02, HepG2, or TP53^-/- or wild-type HCT116 cells, and knocked down with short hairpin RNAs in HepG2, Bel-7402, and SMMC-7721 cells; proliferation, colony formation, and apoptosis were measured. Formation of xenograft tumors from cell lines was monitored in nude mice. Subcellular localization of ribosome proteins and levels or activity of p53 were investigated by co-immunoprecipitation, immunofluorescence, and immunoblot analyses. Levels of small nucleolar RNA H/ACA box 18-like 5 (SNORA18L5) were quantified by quantitative reverse transcription polymerase chain reaction.

RESULTS

We identified a low-frequency duplication at chromosome 15q13.3 strongly associated with risk of HBV-related HCC (overall P = 3.17 × 10^-8; odds ratio, 12.02). Copy numbers of the 15q13.3 duplication correlated with the expression of SNORA18L5 in liver tissues. Overexpression of SNORA18L5 increased HCC cell proliferation and growth of xenograft tumors in mice; knockdown reduced HCC proliferation and tumor growth. SNORA18L5 overexpression in HepG2 and SMMC-7721 cells inhibited p53-dependent cell cycle arrest and apoptosis. Overexpression of SNORA18L5 led to hyperactive ribosome biogenesis, increasing levels of mature 18S and 28S ribosomal RNAs and causing the ribosomal proteins RPL5 and RPL11 to stay in the nucleolus, which kept them from binding to MDM2. This resulted in increased MDM2-mediated ubiquitination and degradation of p53. Levels of SNORA18L5 were increased in HCC tissues compared with nontumor liver tissues and associated with shorter survival times of patients.

CONCLUSIONS

In a CNV-based GWAS, we associated duplication at 15q13.3 with increased risk of HBV-related HCC. We found SNORA18L5 at this location to promote HCC cell proliferation and tumor growth in mice. SNORA18L5 increases ribosome biogenesis, facilitates ribosomal RNA maturation, and alters localization of RPL5 and RPL11, allowing for increased MDM2-mediated proteolysis of p53 and cell cycle arrest.

Regulation of H3K4me3 at Transcriptional Enhancers Characterizes Acquisition of Virus-Specific CD8+ T Cell-Lineage-Specific Function

Infection triggers large-scale changes in the phenotype and function of T cells that are critical for immune clearance, yet the gene regulatory mechanisms that control these changes are largely unknown. Using ChIP-seq for specific histone post-translational modifications (PTMs), we mapped the dynamics of ∼25,000 putative CD8⁺ T cell transcriptional enhancers (TEs) differentially utilized during virus-specific T cell differentiation. Interestingly, we identified a subset of dynamically regulated TEs that exhibited acquisition of a non-canonical (H3K4me3⁺) chromatin signature upon differentiation. This unique TE subset exhibited characteristics of poised enhancers in the naive CD8⁺ T cell subset and demonstrated enrichment for transcription factor binding motifs known to be important for virus-specific CD8⁺ T cell differentiation. These data provide insights into the establishment and maintenance of the gene transcription profiles that define each stage of virus-specific T cell differentiation.

Repression of human and mouse brain inflammaging transcriptome by broad gene-body histone hyperacetylation

Brain "inflammaging," a low-grade and chronic inflammation, is a major hallmark for aging-related neurodegenerative diseases. Here, by profiling H3K27ac and gene expression patterns in human and mouse brains, we found that age-related up-regulated (Age-Up) and down-regulated (Age-Down) genes have distinct H3K27ac patterns. Although both groups show promoter H3K27ac, the Age-Up genes, enriched for inflammation-related functions, are additionally marked by broad H3K27ac distribution over their gene bodies, which is progressively reduced during aging. Age-related gene expression changes can be predicted by gene-body H3K27ac level. Contrary to the presumed transcription activation function of promoter H3K27ac, we found that broad gene-body hyper H3K27ac suppresses overexpression of inflammaging genes. Altogether, our findings revealed opposite regulations by H3K27ac of Age-Up and Age-Down genes and a mode of broad gene-body H3K27ac in repressing transcription.

Long genes linked to autism spectrum disorders harbor broad enhancer-like chromatin domains

Genetic variants associated with autism spectrum disorders (ASDs) are enriched in genes encoding synaptic proteins and chromatin regulators. Although the role of synaptic proteins in ASDs is widely studied, the mechanism by which chromatin regulators contribute to ASD risk remains poorly understood. Upon profiling and analyzing the transcriptional and epigenomic features of genes expressed in the cortex, we uncovered a unique set of long genes that contain broad enhancer-like chromatin domains (BELDs) spanning across their entire gene bodies. Analyses of these BELD genes show that they are highly transcribed with frequent RNA polymerase II (Pol II) initiation and low Pol II pausing, and they exhibit frequent chromatin-chromatin interactions within their gene bodies. These BELD features are conserved from rodents to humans, are enriched in genes involved in synaptic function, and appear post-natally concomitant with synapse development. Importantly, we find that BELD genes are highly implicated in neurodevelopmental disorders, particularly ASDs, and that their expression is preferentially down-regulated in individuals with idiopathic autism. Finally, we find that the transcription of BELD genes is particularly sensitive to alternations in ASD-associated chromatin regulators. These findings suggest that the epigenomic regulation of BELD genes is important for post-natal cortical development and lend support to a model by which mutations in chromatin regulators causally contribute to ASDs by preferentially impairing BELD gene transcription.

WDR5 supports colon cancer cells by promoting methylation of H3K4 and suppressing DNA damage

BACKGROUND

KMT2/MLL proteins are commonly overexpressed or mutated in cancer and have been shown to support cancer maintenance. These proteins are responsible for methylating histone 3 at lysine 4 and promoting transcription and DNA synthesis; however, they are inactive outside of a multi-protein complex that requires WDR5. WDR5 has been implicated in cancer for its role in the COMPASS complex and its interaction with Myc; however, the role of WDR5 in colon cancer has not yet been elucidated.

METHODS

WDR5 expression was evaluated using RT-qPCR and western blot analysis. Cell viability and colony forming assays were utilized to evaluate the effects of WDR5 depletion or inhibition in colon cancer cells. Downstream effects of WDR5 depletion and inhibition were observed by western blot.

RESULTS

WDR5 is overexpressed in colon tumors and colon cancer cell lines at the mRNA and protein level. WDR5 depletion reduces cell viability in HCT116, LoVo, RKO, HCT15, SW480, SW620, and T84 colon cancer cells. Inhibition of the WDR5:KMT2/MLL interaction using OICR-9429 reduces cell viability in the same panel of cell lines albeit not to the same extent as RNAi-mediated WDR5 depletion. WDR5 depletion reduced H3K4Me3 and increased phosphorylation of H2AX in HCT116, SW620, and RKO colon cancer cells; however, OICR-9429 treatment did not recapitulate these effects in all cell lines potentially explaining the reduced toxicity of OICR-9429 treatment as compared to WDR5 depletion. WDR5 depletion also sensitized colon cancer cells to radiation-induced DNA damage.

CONCLUSIONS

These data demonstrate a clear role for WDR5 in colon cancer and future studies should examine its potential to serve as a therapeutic target in cancer. Additional studies are needed to fully elucidate if the requirement for WDR5 is independent of or consistent with its role within the COMPASS complex. OICR-9429 treatment was particularly toxic to SW620 and T84 colon cancer cells, two cell lines without mutations in WDR5 and KMT2/MLL proteins suggesting COMPASS complex inhibition may be particularly effective in tumors lacking KMT2 mutations. Additionally, the ability of WDR5 depletion to amplify the toxic effects of radiation presents the possibility of targeting WDR5 to sensitize cells to DNA-damaging therapies.

Unique patterns of trimethylation of histone H3 lysine 4 are prone to changes during aging in Caenorhabditis elegans somatic cells

Tri-methylation on histone H3 lysine 4 (H3K4me3) is associated with active gene expression but its regulatory role in transcriptional activation is unclear. Here we used Caenorhabditis elegans to investigate the connection between H3K4me3 and gene expression regulation during aging. We uncovered around 30% of H3K4me3 enriched regions to show significant and reproducible changes with age. We further showed that these age-dynamic H3K4me3 regions largely mark gene-bodies and are acquired during adult stages. We found that these adult-specific age-dynamic H3K4me3 regions are correlated with gene expression changes with age. In contrast, H3K4me3 marking established during developmental stages remained largely stable with age, even when the H3K4me3 associated genes exhibited RNA expression changes during aging. Importantly, the genes associated with changes in H3K4me3 and RNA levels with age are enriched for functional groups commonly implicated in aging biology. Therefore, our findings suggested divergent roles of H3K4me3 in gene expression regulation during aging, with important implications on aging-dependent pathophysiologies.

Histone modifications, the specific chemical modifications on histone proteins, are key for regulating the packing of DNA, and thus have important influence on diverse biological processes. An intensely studied function of histone modifications is their contribution to regulating gene expression. Recent studies in diverse model organisms demonstrated that the global alterations of particular histone modifications, for instance H3K4me3, extend the lifespan of the organism. However, the underlying molecular mechanisms remain largely unclear. In this study, we monitored whether and how the genome-wide pattern of the histone modification H3K4me3 changes during aging in the somatic cells of the model organism C. elegans. We identified interesting and non-conventional patterns of H3K4me3, which span gene-bodies and are acquired during adulthood, that are particularly prone to changes with aging. This is contrasted to the well-studied H3K4me3 patterns that span transcriptional start sites and 5’ promoter regions and are established early during development, which remain stable with age. Consistent with the close association between H3K4me3 marking and active transcription, we observed that the age-dynamic H3K4me3 markings are highly correlated with corresponding RNA expression changes. Importantly, the genes that are associated with both H3K4me3 and RNA expression changes with age are over-represented for functional groups commonly implicated in aging biology. In summary, our findings revealed a lesser known pattern of H3K4me3 modification that can have important biological roles in aging.

MLL4 Is Required to Maintain Broad H3K4me3 Peaks and Super-Enhancers at Tumor Suppressor Genes

Super-enhancers are large clusters of enhancers that activate gene expression. Broad trimethyl histone H3 lysine 4 (H3K4me3) often defines active tumor suppressor genes. However, how these epigenomic signatures are regulated for tumor suppression is little understood. Here we show that brain-specific knockout of the H3K4 methyltransferase MLL4 (a COMPASS-like enzyme, also known as KMT2D) in mice spontaneously induces medulloblastoma. Mll4 loss upregulates oncogenic Ras and Notch pathways while downregulating neuronal gene expression programs. MLL4 enhances DNMT3A-catalyzed DNA methylation and SIRT1/BCL6-mediated H4K16 deacetylation, which antagonize expression of Ras activators and Notch pathway components, respectively. Notably, Mll4 loss downregulates tumor suppressor genes (e.g., Dnmt3a and Bcl6) by diminishing broad H3K4me3 and super-enhancers and also causes widespread impairment of these epigenomic signatures during medulloblastoma genesis. These findings suggest an anti-tumor role for super-enhancers and provide a unique tumor-suppressive mechanism in which MLL4 is necessary to maintain broad H3K4me3 and super-enhancers at tumor suppressor genes.

Distinct epigenomic patterns are associated with haploinsufficiency and predict risk genes of developmental disorders

Haploinsufficiency is a major mechanism of genetic risk in developmental disorders. Accurate prediction of haploinsufficient genes is essential for prioritizing and interpreting deleterious variants in genetic studies. Current methods based on mutation intolerance in population data suffer from inadequate power for genes with short transcripts. Here we show haploinsufficiency is strongly associated with epigenomic patterns, and develop a computational method (Episcore) to predict haploinsufficiency leveraging epigenomic data from a broad range of tissue and cell types by machine learning methods. Based on data from recent exome sequencing studies on developmental disorders, Episcore achieves better performance in prioritizing likely-gene-disrupting (LGD) de novo variants than current methods. We further show that Episcore is less-biased by gene size, and complementary to mutation intolerance metrics for prioritizing LGD variants. Our approach enables new applications of epigenomic data and facilitates discovery and interpretation of novel risk variants implicated in developmental disorders.

Predicting haploinsufficient genes helps to understand the genetic risk underlying developmental disorders. Here, the authors develop a Random Forest-based method that uses epigenomic data to predict haploinsufficiency, Episcore, which is complementary to methods based on mutation intolerance scores.

Methionine metabolism influences genomic architecture and gene expression through H3K4me3 peak width

Nutrition and metabolism are known to influence chromatin biology and epigenetics through post-translational modifications, yet how this interaction influences genomic architecture and connects to gene expression is unknown. Here we consider, as a model, the metabolically-driven dynamics of H3K4me3, a histone methylation mark that is known to encode information about active transcription, cell identity, and tumor suppression. We analyze the genome-wide changes in H3K4me3 and gene expression in response to alterations in methionine availability in both normal mouse physiology and human cancer cells. Surprisingly, we find that the location of H3K4me3 peaks is largely preserved under methionine restriction, while the response of H3K4me3 peak width encodes almost all aspects of H3K4me3 biology including changes in expression levels, and the presence of cell identity and cancer-associated genes. These findings may reveal general principles for how nutrient availability modulates specific aspects of chromatin dynamics to mediate biological function.

OCEAN-C: mapping hubs of open chromatin interactions across the genome reveals gene regulatory networks

We develop a method called open chromatin enrichment and network Hi-C (OCEAN-C) for antibody-independent mapping of global open chromatin interactions. By integrating FAIRE-seq and Hi-C, OCEAN-C detects open chromatin interactions enriched by active cis-regulatory elements. We identify more than 10,000 hubs of open chromatin interactions (HOCIs) in human cells, which are mainly active promoters and enhancers bound by many DNA-binding proteins and form interaction networks crucial for gene transcription. In addition to identifying large-scale topological structures, including topologically associated domains and A/B compartments, OCEAN-C can detect HOCI-mediated chromatin interactions that are strongly associated with gene expression, super-enhancers, and broad H3K4me3 domains.

Human TFIIH Kinase CDK7 Regulates Transcription-Associated Chromatin Modifications

CDK7 phosphorylates the RNA polymerase II (pol II) C-terminal domain CTD and activates the P-TEFb-associated kinase CDK9, but its regulatory roles remain obscure. Here, using human CDK7 analog-sensitive (CDK7as) cells, we observed reduced capping enzyme recruitment, increased pol II promoter-proximal pausing, and defective termination at gene 3' ends upon CDK7 inhibition. We also noted that CDK7 regulates chromatin modifications downstream of transcription start sites. H3K4me3 spreading was restricted at gene 5' ends and H3K36me3 was displaced toward gene 3' ends in CDK7as cells. Mass spectrometry identified factors that bound TFIIH-phosphorylated versus P-TEFb-phosphorylated CTD (versus unmodified); capping enzymes and H3K4 methyltransferase complexes, SETD1A/B, selectively bound phosphorylated CTD, and the H3K36 methyltransferase SETD2 specifically bound P-TEFb-phosphorylated CTD. Moreover, TFIIH-phosphorylated CTD stimulated SETD1A/B activity toward nucleosomes, revealing a mechanistic basis for CDK7 regulation of H3K4me3 spreading. Collectively, these results implicate a CDK7-dependent "CTD code" that regulates chromatin marks in addition to RNA processing and pol II pausing.

H3K36me3-mediated mismatch repair preferentially protects actively transcribed genes from mutation

Histone H3 trimethylation at lysine 36 (H3K36me3) is an important histone mark involved in both transcription elongation and DNA mismatch repair (MMR). It is known that H3K36me3 recruits the mismatch-recognition protein MutSα to replicating chromatin via its physical interaction with MutSα's PWWP domain, but the exact role of H3K36me3 in transcription is undefined. Using ChIP combined with whole-genome DNA sequencing analysis, we demonstrate here that H3K36me3, together with MutSα, is involved in protecting against mutation, preferentially in actively transcribed genomic regions. We found that H3K36me3 and MutSα are much more co-enriched in exons and actively transcribed regions than in introns and nontranscribed regions. The H3K36me3-MutSα co-enrichment correlated with a much lower mutation frequency in exons and actively transcribed regions than in introns and nontranscribed regions. Correspondingly, depleting H3K36me3 or disrupting the H3K36me3-MutSα interaction elevated the spontaneous mutation frequency in actively transcribed genes, but it had little influence on the mutation frequency in nontranscribed or transcriptionally inactive regions. Similarly, H₂O₂-induced mutations, which mainly cause base oxidations, preferentially occurred in actively transcribed genes in MMR-deficient cells. The data presented here suggest that H3K36me3-mediated MMR preferentially safeguards actively transcribed genes not only during replication by efficiently correcting mispairs in early replicating chromatin but also during transcription by directly or indirectly removing DNA lesions associated with a persistently open chromatin structure.

Genome-wide analysis reveals a role for TDG in estrogen receptor-mediated enhancer RNA transcription and 3-dimensional reorganization

Background

The estrogen receptor (ER) is a ligand-dependant transcription factor expressed in many breast cancers and is the target of many endocrine-based cancer therapies. Genome-wide studies have shown that the ER binds to gene-specific enhancer regions in response to β-estradiol (E2) which undergo transcription producing noncoding enhancer RNA (eRNA). While eRNAs are important for transcriptional activation of neighboring genes, the mechanism remains poorly understood.

Results

Using ChIP-Seq we generate a global profile of thymine DNA glycosylase (TDG), an ER coactivator that plays an essential role in DNA demethylation, in response to E2 in the MCF7 breast cancer cell line. Remarkably, we found that in response to E2 TDG localized to enhancers which also recruit ERα, RNA Pol II and other coregulators and which are marked by histone modifications indicative of active enhancers. Importantly, depletion of TDG inhibits E2-mediated transcription of eRNAs and transcription of ER-target genes. Functionally, we find that TDG both sensitizes MCF7 cells to tamoxifen-mediated cytostasis and increases migration and invasion of MCF7 cells.

Conclusions

Taken together we find that TDG plays a central role in mediating transcription at a subset of enhancers and governs how MCF7 cells respond to both estrogenic and anti-estrogenic compounds and may be an effective therapeutic target.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0176-2) contains supplementary material, which is available to authorized users.

Epigenome in Early Mammalian Development: Inheritance, Reprogramming and Establishment

Drastic epigenetic reprogramming takes place during preimplantation development, leading to the conversion of terminally differentiated gametes to a totipotent embryo. Deficiencies in remodeling of the epigenomes can cause severe developmental defects, including embryonic lethality. However, how chromatin modifications and chromatin organization are reprogrammed upon fertilization in mammals has long remained elusive. Here, we review recent progress in understanding how the epigenome is dynamically regulated during early mammalian development. The latest studies, including many from genome-wide perspectives, have revealed unusual principles of reprogramming for histone modifications, chromatin accessibility, and 3D chromatin architecture. These advances have shed light on the regulatory network controlling the earliest development and maternal-zygotic transition.

H3K4 demethylase KDM5B regulates global dynamics of transcription elongation and alternative splicing in embryonic stem cells

Epigenetic regulation of chromatin plays a critical role in controlling embryonic stem (ES) cell self-renewal and pluripotency. However, the roles of histone demethylases and activating histone modifications such as trimethylated histone 3 lysine 4 (H3K4me3) in transcriptional events such as RNA polymerase II (RNAPII) elongation and alternative splicing are largely unknown. In this study, we show that KDM5B, which demethylates H3K4me3, plays an integral role in regulating RNAPII occupancy, transcriptional initiation and elongation, and alternative splicing events in ES cells. Depletion of KDM5B leads to altered RNAPII promoter occupancy, and decreased RNAPII initiation and elongation rates at active genes and at genes marked with broad H3K4me3 domains. Moreover, our results demonstrate that spreading of H3K4me3 from promoter to gene body regions, which is mediated by depletion of KDM5B, modulates RNAPII elongation rates and RNA splicing in ES cells. We further show that KDM5B is enriched nearby alternatively spliced exons, and depletion of KDM5B leads to altered levels of H3K4 methylation in alternatively spliced exon regions, which is accompanied by differential expression of these alternatively splice exons. Altogether, our data indicate an epigenetic role for KDM5B in regulating RNAPII elongation and alternative splicing, which may support the diverse mRNA repertoire in ES cells.

A Functional Role for the Epigenetic Regulator ING1 in Activity-induced Gene Expression in Primary Cortical Neurons

Epigenetic regulation of activity-induced gene expression involves multiple levels of molecular interaction, including histone and DNA modifications, as well as mechanisms of DNA repair. Here we demonstrate that the genome-wide deposition of inhibitor of growth family member 1 (ING1), which is a central epigenetic regulatory protein, is dynamically regulated in response to activity in primary cortical neurons. ING1 knockdown leads to decreased expression of genes related to synaptic plasticity, including the regulatory subunit of calcineurin, Ppp3r1. In addition, ING1 binding at a site upstream of the transcription start site (TSS) of Ppp3r1 depends on yet another group of neuroepigenetic regulatory proteins, the Piwi-like family, which are also involved in DNA repair. These findings provide new insight into a novel mode of activity-induced gene expression, which involves the interaction between different epigenetic regulatory mechanisms traditionally associated with gene repression and DNA repair.

Determinants of Histone H3K4 Methylation Patterns

Various factors differentially recognize trimethylated histone H3 lysine 4 (H3K4me3) near promoters, H3K4me2 just downstream, and promoter-distal H3K4me1 to modulate gene expression. This methylation "gradient" is thought to result from preferential binding of the H3K4 methyltransferase Set1/complex associated with Set1 (COMPASS) to promoter-proximal RNA polymerase II. However, other studies have suggested that location-specific cues allosterically activate Set1. Chromatin immunoprecipitation sequencing (ChIP-seq) experiments show that H3K4 methylation patterns on active genes are not universal or fixed and change in response to both transcription elongation rate and frequency as well as reduced COMPASS activity. Fusing Set1 to RNA polymerase II results in H3K4me2 throughout transcribed regions and similarly extended H3K4me3 on highly transcribed genes. Tethered Set1 still requires histone H2B ubiquitylation for activity. These results show that higher-level methylations reflect not only Set1/COMPASS recruitment but also multiple rounds of transcription. This model provides a simple explanation for non-canonical methylation patterns at some loci or in certain COMPASS mutants.

The impact of cellular metabolism on chromatin dynamics and epigenetics

The substrates used to modify nucleic acids and chromatin are affected by nutrient availability and the activity of metabolic pathways. Thus, cellular metabolism constitutes a fundamental component of chromatin status and thereby of genome regulation. Here we describe the biochemical and genetic principles of how metabolism can influence chromatin biology and epigenetics, discuss the functional roles of this interplay in developmental and cancer biology, and present future directions in this rapidly emerging area.