Repressive Histone Modifications Research Articles

Rewiring of chromatin state and gene expression by transposable elements.

Transposable elements form a major fraction of the genome in various eukaryotic species. Although deleterious effects of transpositions within the genome have been reported, recent findings suggest that transposable elements can function as novel regulatory elements to fine-tune gene expression. Transposable elements can impact the chromatin state through processes such as heterochromatin formation, enhancer-promoter interactions, and chromatin boundary formation, mainly because of the functions of chromatin-based pathways that regulate the expression of transposable elements via DNA methylation and repressive histone modification. Therefore, transposable elements can rewire the chromatin state and gene expression depending on their insertions. Here, we review the findings that reveal the role of transposable elements as modifiers of the chromatin state and gene expression as well as the molecular mechanisms capable of inducing these changes.

Replication timing organizes epigenome

Epigenome The temporal order of DNA replication is conserved from yeast to humans, but its biological significance remains unclear. Klein et al. eliminated the protein RIF1, a master regulator of replication timing, in several human cell lines. RIF1 loss during the G1 phase of the cell cycle resulted in a heterogeneous, nearly random replication timing program from the first S phase that persisted even in stable RIF1-null clones. Altered replication timing was followed by replication-dependent redistribution of active and repressive histone modifications and alterations in genome architecture. These results support a model in which replication timing orchestrates the epigenetic state of newly replicated chromatin. Science , this issue p. [371][1] [1]: /lookup/doi/10.1126/science.aba5545

TRIM28 Expression on Dendritic Cells Prevents Excessive T Cell Priming by Silencing Endogenous Retrovirus.

Acquired immune reaction is initiated by dendritic cells (DCs), which present Ags to a few naive Ag-specific T cells. Deregulation of gene expression in DCs may alter the outcome of the immune response toward immunodeficiency and/or autoimmune diseases. Expression of TRIM28, a nuclear protein that mediates gene silencing through heterochromatin, decreased in DCs from old mice, suggesting alteration of gene regulation. Mice specifically lacking TRIM28 in DCs show increased DC population in the spleen and enhanced T cell priming toward inflammatory effector T cells, leading to acceleration and exacerbation in experimental autoimmune encephalomyelitis. TRIM28-deficient DCs were found to ectopically transcribe endogenous retrovirus (ERV) elements. Combined genome-wide analysis revealed a strong colocalization among the decreased repressive histone mark H3K9me3-transcribed ERV elements and the derepressed host genes that were related to inflammation in TRIM28-deficient DCs. This suggests that TRIM28 occupancy of ERV elements critically represses expression of proximal inflammatory genes on the genome. We propose that gene silencing through repressive histone modification by TRIM28 plays a role in maintaining the integrity of precise gene regulation in DCs, which prevents aberrant T cell priming to inflammatory effector T cells.

Genome-wide programmable transcriptional memory by CRISPR-based epigenome editing

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

H3K36 methylation reprograms gene expression to drive early gametocyte development in Plasmodium falciparum

BackgroundThe Plasmodium sexual gametocyte stages are the only transmissible form of the malaria parasite and are thus responsible for the continued transmission of the disease. Gametocytes undergo extensive functional and morphological changes from commitment to maturity, directed by an equally extensive control program. However, the processes that drive the differentiation and development of the gametocyte post-commitment, remain largely unexplored. A previous study reported enrichment of H3K36 di- and tri-methylated (H3K36me2&3) histones in early-stage gametocytes. Using chromatin immunoprecipitation followed by high-throughput sequencing, we identify a stage-specific association between these repressive histone modifications and transcriptional reprogramming that define a stage II gametocyte transition point.ResultsHere, we show that H3K36me2 and H3K36me3 from stage II gametocytes are associated with repression of genes involved in asexual proliferation and sexual commitment, indicating that H3K36me2&3-mediated repression of such genes is essential to the transition from early gametocyte differentiation to intermediate development. Importantly, we show that the gene encoding the transcription factor AP2-G as commitment master regulator is enriched with H3K36me2&3 and actively repressed in stage II gametocytes, providing the first evidence of ap2-g gene repression in post-commitment gametocytes. Lastly, we associate the enhanced potency of the pan-selective Jumonji inhibitor JIB-04 in gametocytes with the inhibition of histone demethylation including H3K36me2&3 and a disruption of normal transcriptional programs.ConclusionsTaken together, our results provide the first description of an association between global gene expression reprogramming and histone post-translational modifications during P. falciparum early sexual development. The stage II gametocyte-specific abundance of H3K36me2&3 manifests predominantly as an independent regulatory mechanism targeted towards genes that are repressed post-commitment. H3K36me2&3-associated repression of genes is therefore involved in key transcriptional shifts that accompany the transition from early gametocyte differentiation to intermediate development.

Role of DNA methylation-related chromatin remodeling in aryl hydrocarbon receptor-dependent regulation of T-2 toxin highly inducible Cytochrome P450 1A4 gene.

Mycotoxins are toxic secondary metabolites produced by food-contaminating fungi, which lead to global epigenetic changes and cause toxicity to both farm animals and humans. However, whether mycotoxins induce gene-specific epigenetic alterations associated with inducible downstream gene expression is unclear as are the underlying regulatory mechanisms. Here, we found that T-2 toxin and its deacetylated metabolites but not deoxynivalenol (DON) or other representative mycotoxins highly induced the expression of cytochrome P450 1A4 (CYP1A4) in both Leghorn male hepatoma (LMH) cells and chicken primary hepatocytes, and this effect was related to the regulation of both aryl hydrocarbon receptor (AhR) and DNA methylation. We used methylation-sensitive restriction enzyme digestion-qPCR (MSRE-qPCR) and chromatin immunoprecipitation (ChIP) assays and found that the binding of DNA methyltransferase 1 (DNMT1) and histone deacetylase 2 (HDAC2) to highly methylated CpG island 3-2 at the enhancer of CYP1A4 was accompanied by the recruitment of the repressive histone modification marker H3K27me3, inducing a silent state. In turn, T-2 toxin stimulation enriched the binding of AhR to demethylated CpG island 3-2, which facilitated p300 and H3K9ac recruitment and ultimately generated an activated chromatin structure at the enhancer by increasing the active histone modification markers, including H3K4me3, H3K27ac, and H3K14ac. Interestingly, T-2 toxin-induced AhR activation also facilitated RNA polymerase II binding to CpG island 2, which may form a transcriptionally active chromatin structure at the promoter and ultimately transactivate CYP1A4. Our findings provide novel insights into the epigenetic regulation of T-2 toxin-induced gene expression.

Integrative pan cancer analysis reveals epigenomic variation in cancer type and cell specific chromatin domains

Epigenetic mechanisms contribute to the initiation and development of cancer, and epigenetic variation promotes dynamic gene expression patterns that facilitate tumor evolution and adaptation. While the NCI-60 panel represents a diverse set of human cancer cell lines that has been used to screen chemical compounds, a comprehensive epigenomic atlas of these cells has been lacking. Here, we report an integrative analysis of 60 human cancer epigenomes, representing a catalog of activating and repressive histone modifications. We identify genome-wide maps of canonical sharp and broad H3K4me3 domains at promoter regions of tumor suppressors, H3K27ac-marked conventional enhancers and super enhancers, and widespread inter-cancer and intra-cancer specific variability in H3K9me3 and H4K20me3-marked heterochromatin domains. Furthermore, we identify features of chromatin states, including chromatin state switching along chromosomes, correlation of histone modification density with genetic mutations, DNA methylation, enrichment of DNA binding motifs in regulatory regions, and gene activity and inactivity. These findings underscore the importance of integrating epigenomic maps with gene expression and genetic variation data to understand the molecular basis of human cancer. Our findings provide a resource for mining epigenomic maps of human cancer cells and for identifying epigenetic therapeutic targets.

Bisphenol A Exposure Disrupts Enamel Formation via EZH2-Mediated H3K27me3

Enamel formation is a serial and complex biological process, during which related genes are expressed progressively in a spatiotemporal manner. This process is vulnerable to environmental cues, resulting in developmental defects of enamel (DDE). However, how environmental factors are biologically integrated during enamel formation is still poorly understood. Here, we investigated the mechanism of DDE elicited by a model endocrine-disrupting chemical, bisphenol A (BPA), in mouse incisors. We show that BPA exposure leads to DDE in mouse incisors, as well as excessive proliferation in dental epithelial stem/progenitor cells. Western blotting, chromatin immunoprecipitation sequencing, and immunofluorescence staining revealed that this effect was accompanied by upregulation of a repressive mark, H3K27me3, in the labial cervical loop of mouse incisors. Perturbation of H3K27me3 methyltransferase EZH2 repressed the level of H3K27me3 and partially attenuated the excessive proliferation in dental epithelial stem/progenitor cells and DDE induced by BPA exposure. Overall, our results demonstrate the essential role of repressive histone modification H3K27me3 in DDE elicited by exposure to an endocrine-disrupting chemical.

Genome Triplication Leads to Transcriptional Divergence of FLOWERING LOCUS C Genes During Vernalization in the Genus Brassica.

The genus Brassica includes oil crops, vegetables, condiments, fodder crops, and ornamental plants. Brassica species underwent a whole genome triplication event after speciation between ancestral species of Brassica and closely related genera including Arabidopsis thaliana. Diploid species such as Brassica rapa and Brassica oleracea have three copies of genes orthologous to each A. thaliana gene, although deletion in one or two of the three homologs has occurred in some genes. The floral transition is one of the crucial events in a plant’s life history, and time of flowering is an important agricultural trait. There is a variation in flowering time within species of the genus Brassica, and this variation is largely dependent on a difference in vernalization requirements. In Brassica, like in A. thaliana, the key gene of vernalization is FLOWERING LOCUS C (FLC). In Brassica species, the vernalization response including the repression of FLC expression by cold treatment and the enrichment of the repressive histone modification tri-methylated histone H3 lysine 27 (H3K27me3) at the FLC locus is similar to A. thaliana. B. rapa and B. oleracea each have four paralogs of FLC, and the allotetraploid species, Brassica napus, has nine paralogs. The increased number of paralogs makes the role of FLC in vernalization more complicated; in a single plant, paralogs vary in the expression level of FLC before and after vernalization. There is also variation in FLC expression levels between accessions. In this review, we focus on the regulatory circuits of the vernalization response of FLC expression in the genus Brassica.

Histone modification dynamics at H3K27 are associated with altered transcription of in planta induced genes in Magnaporthe oryzae.

Transcriptional dynamic in response to environmental and developmental cues are fundamental to biology, yet many mechanistic aspects are poorly understood. One such example is fungal plant pathogens, which use secreted proteins and small molecules, termed effectors, to suppress host immunity and promote colonization. Effectors are highly expressed in planta but remain transcriptionally repressed ex planta, but our mechanistic understanding of these transcriptional dynamics remains limited. We tested the hypothesis that repressive histone modification at H3-Lys27 underlies transcriptional silencing ex planta, and that exchange for an active chemical modification contributes to transcription of in planta induced genes. Using genetics, chromatin immunoprecipitation and sequencing and RNA-sequencing, we determined that H3K27me3 provides significant local transcriptional repression. We detail how regions that lose H3K27me3 gain H3K27ac, and these changes are associated with increased transcription. Importantly, we observed that many in planta induced genes were marked by H3K27me3 during axenic growth, and detail how altered H3K27 modification influences transcription. ChIP-qPCR during in planta growth suggests that H3K27 modifications are generally stable, but can undergo dynamics at specific genomic locations. Our results support the hypothesis that dynamic histone modifications at H3K27 contributes to fungal genome regulation and specifically contributes to regulation of genes important during host infection.

Chromatin Regulation in Development: Current Understanding and Approaches.

The regulation of mammalian stem cell fate during differentiation is complex and can be delineated across many levels. At the chromatin level, the replacement of histone variants by chromatin-modifying proteins, enrichment of specific active and repressive histone modifications, long-range gene interactions, and topological changes all play crucial roles in the determination of cell fate. These processes control regulatory elements of critical transcriptional factors, thereby establishing the networks unique to different cell fates and initiate waves of distinctive transcription events. Due to the technical challenges posed by previous methods, it was difficult to decipher the mechanism of cell fate determination at early embryogenesis through chromatin regulation. Recently, single-cell approaches have revolutionised the field of developmental biology, allowing unprecedented insights into chromatin structure and interactions in early lineage segregation events during differentiation. Here, we review the recent technological advancements and how they have furthered our understanding of chromatin regulation during early differentiation events.

Dicer-independent snRNA/snoRNA-derived nuclear RNA 3 regulates tumor-associated macrophage function by epigenetically repressing inducible nitric oxide synthase transcription.

BackgroundSmall RNAs (sRNAs) extensively mediate gene‐specific chromatin regulation in lower organisms. As a dominant type of functional sRNAs in mature mammals, microRNAs mainly regulate gene expression at post‐transcription level in the cytoplasm. Currently, whether there exists a type of nuclear‐localized sRNAs mediating gene‐specific epigenetic regulation in mature mammalian cells remains largely unclear. Here, we profiled sRNAs enriched in the nucleus and investigated their function in mediating gene‐specific epigenetic regulation in anti‐tumor immunity.MethodsWe established cytoplasmic and nuclear transcriptomes of sRNAs of dendritic cells (DCs) using high‐throughput sequencing. Transcription abundances of sRNAs and mRNAs were analyzed by reverse transcription‐quantitative polymerase chain reaction (RT‐qPCR) assay. The associations between sRNAs and Argonaute (AGO) proteins were detected by RNA immunoprecipitation analysis. Synthesized sRNAs and locked nucleic acid (LNA) ‐modified sRNA inhibitors were used to screen the function of sRNAs in innate immune cells. The effect of sRNA on the enrichment of either chromatin remodeler or histone modification at the gene promoter was analyzed by chromatin immunoprecipitation (ChIP)‐qPCR assay. Chromatin accessibility qPCR assay was used to detect the accessibility of gene promoters. A B16 melanoma‐bearing mouse model was established to determine the function of sRNAs in tumor‐associated macrophages (TAMs) and their effect on tumor growth.ResultsWe identified a new class of nucleus‐localized sRNAs, named snRNA/snoRNA‐derived nuclear RNAs (sdnRNAs). Some sdnRNAs were Dicer‐independent and had no association with Argonaute proteins. sdnRNA‐3, the most abundant Dicer‐independent sdnRNAs identified in our analysis, was selected as a representative to examine the biological function of sdnRNAs. sdnRNA‐3 selectively inhibited the transcription of Nos2 in macrophages during innate immune response by repressing the chromatin accessibility at Nos2 gene promoter. sdnRNA‐3 promoted the enrichments of repressive chromatin‐remodeling regulator Mi‐2β and the repressive histone modification H3K27me3 at Nos2 gene promoter. In the B16 melanoma mouse model, we found higher expression of sdnRNA‐3 in M2 TAMs than M1 TAMs and DCs. Transfer of sdnRNA‐3‐silenced macrophages inhibited tumor growth with increased expression of inducible nitric oxide synthase (iNOS) in TAMs.ConclusionsOur results demonstrated that the sdnRNA‐3 repressed the transcription of Nos2 by repressing chromatin accessibility at the promoter, providing new insights into the regulation of macrophage function in tumor immunity.

Effect of Initial Seeding Density on Cell Behavior-Driven Epigenetic Memory and Preferential Lineage Differentiation of Human iPSCs

Understanding the cellular behavioral mechanisms underlying memory formation and maintenance in human induced pluripotent stem cell (hiPSC) culture provides key strategies for achieving stability and robustness of cell differentiation. Here, we show that changes in cell behavior-driven epigenetic memory of hiPSC cultures alter their pluripotent state and subsequent differentiation. Interestingly, pluripotency-associated genes were activated during the entire cell growth phase along with increased active modifications and decreased repressive modifications. This memory effect can last several days in the long-term stationary phase and is sustained in the aspect of cell behavioral changes after subculture. Further, changes in growth-related cell behavior were found to induce nucleoskeletal reorganization, thereby enabling hiPSCs to change their ratio of active and repressive histone modifications and their differentiation potential. Overall, we discuss the cell behavior-driven epigenetic memory induced by the culture environment, and the effect of previous memory on cell lineage specification in the process of hiPSC differentiation.

Mathematical Modeling of Histone Modifications Reveals the Formation Mechanism and Function of Bivalent Chromatin

Bivalent chromatin is characterized by occupation of both activating histone modifications and repressive histone modifications. While bivalent chromatin is known to link with many biological processes, the mechanisms responsible for its multiple functions remain unclear. Here, we develop a mathematical model that involves antagonistic histone modifications H3K4me3 and H3K27me3 to capture the key features of bivalent chromatin. Three necessary conditions for the emergence of bivalent chromatin are identified, including advantageous methylating activity over demethylating activity, frequent noise conversions of modifications, and sufficient nonlinearity. The first condition is further confirmed by analyzing the experimental data from a recent study. Investigation of the composition of bivalent chromatin reveals that bivalent nucleosomes carrying both H3K4me3 and H3K27me3 account for no more than half of nucleosomes at the bivalent chromatin domain. We identify that bivalent chromatin not only allows transitions to multiple states but also serves as a stepping stone to facilitate a step-wise transition between repressive chromatin state and activating chromatin state, and thus elucidate crucial roles of bivalent chromatin in mediating phenotypical plasticity during cell fate determination.

Pathogenic Impacts of Dysregulated Polycomb Repressive Complex Function in Hematological Malignancies.

Polycomb repressive complexes (PRCs) are epigenetic regulators that mediate repressive histone modifications. PRCs play a pivotal role in the maintenance of hematopoietic stem cells through repression of target genes involved in cell proliferation and differentiation. Next-generation sequencing technologies have revealed that various hematologic malignancies harbor mutations in PRC2 genes, such as EZH2, EED, and SUZ12, and PRC1.1 genes, such as BCOR and BCORL1. Except for the activating EZH2 mutations detected in lymphoma, most of these mutations compromise PRC function and are frequently associated with resistance to chemotherapeutic agents and poor prognosis. Recent studies have shown that mutations in PRC genes are druggable targets. Several PRC2 inhibitors, including EZH2-specific inhibitors and EZH1 and EZH2 dual inhibitors have shown therapeutic efficacy for tumors with and without activating EZH2 mutations. Moreover, EZH2 loss-of-function mutations appear to be attractive therapeutic targets for implementing the concept of synthetic lethality. Further understanding of the epigenetic dysregulation associated with PRCs in hematological malignancies should improve treatment outcomes.

Identification of a unique epigenetic profile in women with diminished ovarian reserve

To investigate whether epigenetic profiles of mural granulosa cells (MGC) and leukocytes from women with diminished ovarian reserve (DOR) differ from those of women with normal or high ovarian reserve. Prospectively collected material from a multicenter cohort of women undergoing fertility treatment. Private and university-based facilities for clinical services and research. One hundred and nineteen women of various ages and ovarian reserve status (antimüllerian hormone level) who provided blood samples and MGC. None. Measures of epigenetic aging rates from whole-genome methylation array data: DNA methylation variability, age acceleration, DNA methylation telomere length estimator (DNAmTL), and accumulation of epimutations. Comparison of DOR or high ovarian reserve samples to controls (normal ovarian reserve) showed differential methylation variability between DOR and normal samples at 4,199 CpGs in MGC, and 447 between high and normal (false-discovery rate < 0.05). Variable sites in MGC from DOR were enriched in regions marked with the repressive histone modification H3K27me3, and also included genes involved in folliculogenesis, such as insulin growth factor 2 (IGF2) and antimüllerian hormone (AMH). Regardless of ovarian reserve, very few signals were detected in leukocytes, and no overlaps with those in MGC were found. Furthermore, we found a higher number of epimutations in MGC from women with DOR (Kruskal-Wallis test, difference in mean = 3,485). The somatic cells of human ovarian follicles have a distinctive epigenetic profile in women with DOR. A high frequency of epimutations suggests premature aging. Ovarian reserve status was not reflected in the leukocyte epigenetic profile.

Epigenetic Mechanisms Mediating Cell State Transitions in Chondrocytes.

Epigenetic modifications play critical roles in regulating cell lineage differentiation, but the epigenetic mechanisms guiding specific differentiation steps within a cell lineage have rarely been investigated. To decipher such mechanisms, we used the defined transition from proliferating (PC) into hypertrophic chondrocytes (HC) during endochondral ossification as a model. We established a map of activating and repressive histone modifications for each cell type. ChromHMM state transition analysis and Pareto-based integration of differential levels of mRNA and epigenetic marks revealed that differentiation-associated gene repression is initiated by the addition of H3K27me3 to promoters still carrying substantial levels of activating marks. Moreover, the integrative analysis identified genes specifically expressed in cells undergoing the transition into hypertrophy. Investigation of enhancer profiles detected surprising differences in enhancer number, location, and transcription factor binding sites between the two closely related cell types. Furthermore, cell type-specific upregulation of gene expression was associated with increased numbers of H3K27ac peaks. Pathway analysis identified PC-specific enhancers associated with chondrogenic genes, whereas HC-specific enhancers mainly control metabolic pathways linking epigenetic signature to biological functions. Since HC-specific enhancers show a higher conservation in postnatal tissues, the switch to metabolic pathways seems to be a hallmark of differentiated tissues. Surprisingly, the analysis of H3K27ac levels at super-enhancers revealed a rapid adaption of H3K27ac occupancy to changes in gene expression, supporting the importance of enhancer modulation for acute alterations in gene expression. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

KRAB zinc finger protein diversification drives mammalian interindividual methylation variability

Most transposable elements (TEs) in the mouse genome are heavily modified by DNA methylation and repressive histone modifications. However, a subset of TEs exhibit variable methylation levels in genetically identical individuals, and this is associated with epigenetically conferred phenotypic differences, environmental adaptability, and transgenerational epigenetic inheritance. The evolutionary origins and molecular mechanisms underlying interindividual epigenetic variability remain unknown. Using a repertoire of murine variably methylated intracisternal A-particle (VM-IAP) epialleles as a model, we demonstrate that variable DNA methylation states at TEs are highly susceptible to genetic background effects. Taking a classical genetics approach coupled with genome-wide analysis, we harness these effects and identify a cluster of KRAB zinc finger protein (KZFP) genes that modifies VM-IAPs in trans in a sequence-specific manner. Deletion of the cluster results in decreased DNA methylation levels and altered histone modifications at the targeted VM-IAPs. In some cases, these effects are accompanied by dysregulation of neighboring genes. We find that VM-IAPs cluster together phylogenetically and that this is linked to differential KZFP binding, suggestive of an ongoing evolutionary arms race between TEs and this large family of epigenetic regulators. These findings indicate that KZFP divergence and concomitant evolution of DNA binding capabilities are mechanistically linked to methylation variability in mammals, with implications for phenotypic variation and putative paradigms of mammalian epigenetic inheritance.

Interactive analysis of single-cell epigenomic landscapes with ChromSCape

Chromatin modifications orchestrate the dynamic regulation of gene expression during development and in disease. Bulk approaches have characterized the wide repertoire of histone modifications across cell types, detailing their role in shaping cell identity. However, these population-based methods do not capture cell-to-cell heterogeneity of chromatin landscapes, limiting our appreciation of the role of chromatin in dynamic biological processes. Recent technological developments enable the mapping of histone marks at single-cell resolution, opening up perspectives to characterize the heterogeneity of chromatin marks in complex biological systems over time. Yet, existing tools used to analyze bulk histone modifications profiles are not fit for the low coverage and sparsity of single-cell epigenomic datasets. Here, we present ChromSCape, a user-friendly interactive Shiny/R application distributed as a Bioconductor package, that processes single-cell epigenomic data to assist the biological interpretation of chromatin landscapes within cell populations. ChromSCape analyses the distribution of repressive and active histone modifications as well as chromatin accessibility landscapes from single-cell datasets. Using ChromSCape, we deconvolve chromatin landscapes within the tumor micro-environment, identifying distinct H3K27me3 landscapes associated with cell identity and breast tumor subtype.

Broad genic repression domains signify enhanced silencing of oncogenes

Cancers result from a set of genetic and epigenetic alterations. Most known oncogenes were identified by gain-of-function mutations in cancer, yet little is known about their epigenetic features. Through integrative analysis of 11,596 epigenomic profiles and mutations from >8200 tumor-normal pairs, we discover broad genic repression domains (BGRD) on chromatin as an epigenetic signature for oncogenes. A BGRD is a widespread enrichment domain of the repressive histone modification H3K27me3 and is further enriched with multiple other repressive marks including H3K9me3, H3K9me2, and H3K27me2. Further, BGRD displays widespread enrichment of repressed cis-regulatory elements. Shortening of BGRDs is linked to derepression of transcription. BGRDs at oncogenes tend to be conserved across normal cell types. Putative tumor-promoting genes and lncRNAs defined using BGRDs are experimentally verified as required for cancer phenotypes. Therefore, BGRDs play key roles in epigenetic regulation of cancer and provide a direction for mutation-independent discovery of oncogenes.