Dynamics of broad H3K4me3 domains uncover an epigenetic switch between cell identity and cancer-related genes

Chromatin profiling and state predictions reveal insights into epigenetic regulation during early porcine development

BACKGROUND

Given their physiological similarities to humans, pigs are increasingly used as model organisms in human-oriented biomedical studies. Additionally, their value to animal agriculture across the globe has led to the development of numerous studies to investigate how to improve livestock welfare and production efficiency. As such, pigs are uniquely poised as compelling models that can yield findings with potential implications in both human and animal contexts. Despite this, many gaps remain in our knowledge about the foundational mechanisms that govern gene expression in swine across different developmental stages, particularly in early development. To address some of these gaps, we profiled the histone marks H3K4me3, H3K27ac, and H3K27me3 and the SWI/SNF central ATPase BRG1 in two porcine cell lines representing discrete early developmental time points and used the resulting information to construct predicted chromatin state maps for these cells. We combined this approach with analysis of publicly available RNA-seq data to examine the relationship between epigenetic status and gene expression in these cell types.

RESULTS

In porcine fetal fibroblast (PFF) and trophectoderm cells (PTr2), we saw expected patterns of enrichment for each of the profiled epigenetic features relative to specific genomic regions. H3K4me3 was primarily enriched at and around global gene promoters, H3K27ac was enriched in promoter and intergenic regions, H3K27me3 had broad stretches of enrichment across the genome and narrower enrichment patterns in and around the promoter regions of some genes, and BRG1 primarily had detectable enrichment at and around promoter regions and in intergenic stretches, with many instances of H3K27ac co-enrichment. We used this information to perform genome-wide chromatin state predictions for 10 different states using ChromHMM. Using the predicted chromatin state maps, we identified a subset of genomic regions marked by broad H3K4me3 enrichment, and annotation of these regions revealed that they were highly associated with essential developmental processes and consisted largely of expressed genes. We then compared the identities of the genes marked by these regions to genes identified as cell-type-specific using transcriptome data and saw that a subset of broad H3K4me3-marked genes was also specifically expressed in either PFF or PTr2 cells.

CONCLUSIONS

These findings enhance our understanding of the epigenetic landscape present in early swine development and provide insight into how variabilities in chromatin state are linked to cell identity. Furthermore, this data captures foundational epigenetic details in two valuable porcine cell lines and contributes to the growing body of knowledge surrounding the epigenetic landscape in this species.

Histone H4 asymmetrically dimethylated at arginine 3 (H4R3me2a), a mark of super-enhancers

Histone H4 asymmetrically dimethylated at arginine 3 (H4R3me2a) is an active histone mark catalyzed by protein arginine methyltransferase 1 (PRMT1), a major arginine methyltransferase in vertebrates catalyzing asymmetric dimethylation of arginine. H4R3me2a stimulates the activity of lysine acetyltransferases such as CBP/p300, which catalyze the acetylation of H3K27, a mark of active enhancers, super-enhancers, and promoters. There are a few studies on the genomic location of H4R3me2a. In chicken polychromatic erythrocytes, H4R3me2a is found in introns and intergenic regions and binds to the globin locus control region (a super-enhancer) and globin regulatory regions. In this report, we analyzed chromatin immunoprecipitation sequencing data for the genomic location of H4R3me2a in the breast cancer cell line MCF7. As in avian cells, MCF7 H4R3me2a is present in intronic and intergenic regions. Nucleosomes with H4R3me2a and H3K27ac next to nucleosome-free regions are found at super-enhancers, enhancers, and promoter regions of expressed genes. Genes with critical roles in breast cancer cells have broad domains of nucleosomes with H4R3me2a, H3K27ac, and H3K4me3. Our results are consistent with PRMT1-mediated H4R3me2a playing a key role in the chromatin organization of regulatory regions of vertebrate genomes.

Dynamic and distinct histone modifications facilitate human trophoblast lineage differentiation

The placenta serves as an essential organ for fetal growth throughout pregnancy. Histone modification is a crucial regulatory mechanism involved in numerous biological processes and development. Nevertheless, there remains a significant gap in our understanding regarding the epigenetic regulations that influence trophoblast lineage differentiation, a fundamental aspect of placental development. Here, through comprehensive mapping of H3K4me3, H3K27me3, H3K9me3, and H3K27ac loci during the differentiation of trophoblast stem cells (TSCs) into syncytiotrophoblasts (STs) and extravillous trophoblasts (EVTs), we reveal dynamic reconfiguration in H3K4me3 and H3K27ac patterns that establish an epigenetic landscape conducive to proper trophoblast lineage differentiation. We observe that broad H3K4me3 domains are associated with trophoblast lineage-specific gene expression. Unlike embryonic stem cells, TSCs lack robust bivalent domains. Notably, the repression of ST- and EVT-active genes in TSCs is primarily attributed to the weak H3K4me3 signal rather than bivalent domains. We also unveil the inactivation of TSC enhancers precedes the activation of ST enhancers during ST formation. Our results provide a comprehensive global map of diverse histone modifications, elucidating the dynamic histone modifications during trophoblast lineage differentiation.

Roles and Regulation of H3K4 Methylation During Mammalian Early Embryogenesis and Embryonic Stem Cell Differentiation

From generation of germ cells, fertilization, and throughout early mammalian embryonic development, the chromatin undergoes significant alterations to enable precise regulation of gene expression and genome use. Methylation of histone 3 lysine 4 (H3K4) correlates with active regions of the genome, and it has emerged as a dynamic mark throughout this timeline. The pattern and the level of H3K4 methylation are regulated by methyltransferases and demethylases. These enzymes, as well as their protein partners, play important roles in early embryonic development and show phenotypes in embryonic stem cell self-renewal and differentiation. The various roles of H3K4 methylation are interpreted by dedicated chromatin reader proteins, linking this modification to broader molecular and cellular phenotypes. In this review, we discuss the regulation of different levels of H3K4 methylation, their distinct accumulation pattern, and downstream molecular roles with an early embryogenesis perspective.

Broad H3K4me3 domains: Maintaining cellular identity and their implication in super-enhancer hijacking

The human and mouse genomes are complex from a genomic standpoint. Each cell has the same genomic sequence, yet a wide array of cell types exists due to the presence of a plethora of regulatory elements in the non-coding genome. Recent advances in epigenomic profiling have uncovered non-coding gene proximal promoters and distal enhancers of transcription genome-wide. Extension of promoter-associated H3K4me3 histone mark across the gene body, known as a broad H3K4me3 domain (H3K4me3-BD), is a signature of constitutive expression of cell-type-specific regulation and of tumour suppressor genes in healthy cells. Recently, it has been discovered that the presence of H3K4me3-BDs over oncogenes is a cancer-specific feature associated with their dysregulated gene expression and tumourigenesis. Moreover, it has been shown that the hijacking of clusters of enhancers, known as super-enhancers (SE), by proto-oncogenes results in the presence of H3K4me3-BDs over the gene body. Therefore, H3K4me3-BDs and SE crosstalk in healthy and cancer cells therefore represents an important mechanism to identify future treatments for patients with SE driven cancers.

Distinct chromatin signatures in the Arabidopsis male gametophyte

Epigenetic reprogramming in the germline contributes to the erasure of epigenetic inheritance across generations in mammals but remains poorly characterized in plants. Here we profiled histone modifications throughout Arabidopsis male germline development. We find that the sperm cell has widespread apparent chromatin bivalency, which is established by the acquisition of H3K27me3 or H3K4me3 at pre-existing H3K4me3 or H3K27me3 regions, respectively. These bivalent domains are associated with a distinct transcriptional status. Somatic H3K27me3 is generally reduced in sperm, while dramatic loss of H3K27me3 is observed at only ~700 developmental genes. The incorporation of the histone variant H3.10 facilitates the establishment of sperm chromatin identity without a strong impact on resetting of somatic H3K27me3. Vegetative nuclei harbor thousands of specific H3K27me3 domains at repressed genes, while pollination-related genes are highly expressed and marked by gene body H3K4me3. Our work highlights putative chromatin bivalency and restricted resetting of H3K27me3 at developmental regulators as key features in plant pluripotent sperm.

Bioinformatic Analyses of Broad H3K79me2 Domains in Different Leukemia Cell Line Data Sets

A subset of expressed genes is associated with a broad H3K4me3 (histone H3 trimethylated at lysine 4) domain that extends throughout the gene body. Genes marked in this way in normal cells are involved in cell-identity and tumor-suppressor activities, whereas in cancer cells, genes driving the cancer phenotype (oncogenes) have this feature. Other histone modifications associated with expressed genes that display a broad domain have been less studied. Here, we identified genes with the broadest H3K79me2 (histone H3 dimethylated at lysine 79) domain in human leukemic cell lines representing different forms of leukemia. Taking a bioinformatic approach, we provide evidence that genes with the broadest H3K79me2 domain have known roles in leukemia (e.g., JMJD1C). In the mixed-lineage leukemia cell line MOLM-13, the HOXA9 gene is in a 100 kb broad H3K79me2 domain with other HOXA protein-coding and oncogenic long non-coding RNA genes. The genes in this domain contribute to leukemia. This broad H3K79me2 domain has an unstable chromatin structure, as was evident by enhanced chromatin accessibility throughout. Together, we provide evidence that identification of genes with the broadest H3K79me2 domain will aid in generating a panel of genes in the diagnosis and therapeutic treatment of leukemia in the future.

Functional genomics uncovers the transcription factor BNC2 as required for myofibroblastic activation in fibrosis

Tissue injury triggers activation of mesenchymal lineage cells into wound-repairing myofibroblasts, whose unrestrained activity leads to fibrosis. Although this process is largely controlled at the transcriptional level, whether the main transcription factors involved have all been identified has remained elusive. Here, we report multi-omics analyses unraveling Basonuclin 2 (BNC2) as a myofibroblast identity transcription factor. Using liver fibrosis as a model for in-depth investigations, we first show that BNC2 expression is induced in both mouse and human fibrotic livers from different etiologies and decreases upon human liver fibrosis regression. Importantly, we found that BNC2 transcriptional induction is a specific feature of myofibroblastic activation in fibrotic tissues. Mechanistically, BNC2 expression and activities allow to integrate pro-fibrotic stimuli, including TGFβ and Hippo/YAP1 signaling, towards induction of matrisome genes such as those encoding type I collagen. As a consequence, Bnc2 deficiency blunts collagen deposition in livers of mice fed a fibrogenic diet. Additionally, our work establishes BNC2 as potentially druggable since we identified the thalidomide derivative CC-885 as a BNC2 inhibitor. Altogether, we propose that BNC2 is a transcription factor involved in canonical pathways driving myofibroblastic activation in fibrosis.

Myofibroblasts contribute to the development of liver fibrosis. Here, the authors report that the transcription factor Basonuclin 2 (BNC2) integrates fibrogenic signals and drives myofibroblastic transcriptional activation in liver fibrosis.

Histone H3K4me3 breadth in hypoxia reveals endometrial core functions and stress adaptation linked to endometriosis

Trimethylation of histone H3 at lysine 4 (H3K4me3) is a marker of active promoters. Broad H3K4me3 promoter domains have been associated with cell type identity, but H3K4me3 dynamics upon cellular stress have not been well characterized. We assessed this by exposing endometrial stromal cells to hypoxia, which is a major cellular stress condition. We observed that hypoxia modifies the existing H3K4me3 marks and that promoter H3K4me3 breadth rather than height correlates with transcription. Broad H3K4me3 domains mark genes for endometrial core functions and are maintained or selectively extended upon hypoxia. Hypoxic extension of H3K4me3 breadth associates with stress adaptation genes relevant for the survival of endometrial cells including transcription factor KLF4, for which we found increased protein expression in the stroma of endometriosis lesions. These results substantiate the view on broad H3K4me3 as a marker of cell identity genes and reveal participation of H3K4me3 extension in cellular stress adaptation.

Epigenomic translocation of H3K4me3 broad domains over oncogenes following hijacking of super-enhancers

Chromosomal translocations are important drivers of haematological malignancies whereby proto-oncogenes are activated by juxtaposition with enhancers, often called enhancer hijacking. We analyzed the epigenomic consequences of rearrangements between the super-enhancers of the immunoglobulin heavy locus (IGH) and proto-oncogene CCND1 that are common in B cell malignancies. By integrating BLUEPRINT epigenomic data with DNA breakpoint detection, we characterized the normal chromatin landscape of the human IGH locus and its dynamics after pathological genomic rearrangement. We detected an H3K4me3 broad domain (BD) within the IGH locus of healthy B cells that was absent in samples with IGH-CCND1 translocations. The appearance of H3K4me3-BD over CCND1 in the latter was associated with overexpression and extensive chromatin accessibility of its gene body. We observed similar cancer-specific H3K4me3-BDs associated with hijacking of super-enhancers of other common oncogenes in B cell (MAF, MYC, and FGFR3/NSD2) and T cell malignancies (LMO2, TLX3, and TAL1). Our analysis suggests that H3K4me3-BDs can be created by super-enhancers and supports the new concept of epigenomic translocation, in which the relocation of H3K4me3-BDs from cell identity genes to oncogenes accompanies the translocation of super-enhancers.

A Comparative Analysis of Super-Enhancers and Broad H3K4me3 Domains in Pig, Human, and Mouse Tissues

Super-enhancers (SEs) and broad H3K4me3 domains (BDs) are crucial regulators in the control of tissue identity in human and mouse. However, their features in pig remain largely unknown. In this study, by integrative computational analyses of epigenomic and transcriptomic data, we have characterized SEs and BDs in six pig tissues and analyzed their conservation in comparison with human and mouse tissues. Similar to human and mouse, pig SEs and BDs display higher tissue specificity than their typical counterparts. Genes proximal to SEs and BDs are associated with tissue identity in most tissues. About 55-182 SEs (5-17% in total) and 99-309 BDs (8-16% in total) across pig tissues are considered as functionally conserved elements because they have orthologous SEs and BDs in human and mouse. However, these elements do not necessarily exhibit sequence conservation. The functionally conserved SEs are correlated to tissue identity in majority of pig tissues, while those conserved BDs are linked to tissue identity in a few tissues. Our study provides resources for future gene regulatory studies in pig. It highlights that SEs are more effective in defining tissue identity than BDs, which is contrasting to a previous study. It also provides novel insights on understanding the sequence features of functionally conserved elements.

The dynamic broad epigenetic (H3K4me3, H3K27ac) domain as a mark of essential genes

Transcriptionally active chromatin is marked by tri-methylation of histone H3 at lysine 4 (H3K4me3) located after first exons and around transcription start sites. This epigenetic mark is typically restricted to narrow regions at the 5`end of the gene body, though a small subset of genes have a broad H3K4me3 domain which extensively covers the coding region. Although most studies focus on the H3K4me3 mark, the broad H3K4me3 domain is associated with a plethora of histone modifications (e.g., H3 acetylated at K27) and is therein termed broad epigenetic domain. Genes marked with the broad epigenetic domain are involved in cell identity and essential cell functions and have clinical potential as biomarkers for patient stratification. Reducing expression of genes with the broad epigenetic domain may increase the metastatic potential of cancer cells. Enhancers and super-enhancers interact with the broad epigenetic domain marked genes forming a hub of interactions involving nucleosome-depleted regions. Together, the regulatory elements coalesce with transcription factors, chromatin modifying/remodeling enzymes, coactivators, and the Mediator and/or Integrator complex into a transcription factory which may be analogous to a liquid–liquid phase-separated condensate. The broad epigenetic domain has a dynamic chromatin structure which supports frequent transcription bursts. In this review, we present the current knowledge of broad epigenetic domains.