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

Super enhancer acquisition drives expression of oncogenic PPP1R15B that regulates protein homeostasis in multiple myeloma

Multiple myeloma is a hematological malignancy arising from immunoglobulin-secreting plasma cells. It remains poorly understood how chromatin rewiring of regulatory elements contributes to tumorigenesis and therapy resistance in myeloma. Here we generate a high-resolution contact map of myeloma-associated super-enhancers by integrating H3K27ac ChIP-seq and HiChIP from myeloma cell lines, patient-derived myeloma cells and normal plasma cells. Our comprehensive transcriptomic and phenomic analyses prioritize candidate genes with biological and clinical implications in myeloma. We show that myeloma cells frequently acquire SE that transcriptionally activate an oncogene PPP1R15B, which encodes a regulatory subunit of the holophosphatase complex that dephosphorylates translation initiation factor eIF2α. Epigenetic silencing or knockdown of PPP1R15B activates pro-apoptotic eIF2α-ATF4-CHOP pathway, while inhibiting protein synthesis and immunoglobulin production. Pharmacological inhibition of PPP1R15B using Raphin1 potentiates the anti-myeloma effect of bortezomib. Our study reveals that myeloma cells are vulnerable to perturbation of PPP1R15B-dependent protein homeostasis, highlighting a promising therapeutic strategy.

Evolutionarily conserved enhancer-associated features within the MYEOV locus suggest a regulatory role for this non-coding DNA region in cancer

The myeloma overexpressed gene (MYEOV) has been proposed to be a proto-oncogene due to high RNA transcript levels found in multiple cancers, including myeloma, breast, lung, pancreas and esophageal cancer. The presence of an open reading frame (ORF) in humans and other primates suggests protein-coding potential. Yet, we still lack evidence of a functional MYEOV protein. It remains undetermined how MYEOV overexpression affects cancerous tissues. In this work, we show that MYEOV has likely originated and may still function as an enhancer, regulating CCND1 and LTO1. Firstly, MYEOV 3' enhancer activity was confirmed in humans using publicly available ATAC-STARR-seq data, performed on B-cell-derived GM12878 cells. We detected enhancer histone marks H3K4me1 and H3K27ac overlapping MYEOV in multiple healthy human tissues, which include B cells, liver and lung tissue. The analysis of 3D genome datasets revealed chromatin interactions between a MYEOV-3'-putative enhancer and the proto-oncogene CCND1. BLAST searches and multi-sequence alignment results showed that DNA sequence from this human enhancer element is conserved from the amphibians/amniotes divergence, with a 273 bp conserved region also found in all mammals, and even in chickens, where it is consistently located near the corresponding CCND1 orthologues. Furthermore, we observed conservation of an active enhancer state in the MYEOV orthologues of four non-human primates, dogs, rats, and mice. When studying this homologous region in mice, where the ORF of MYEOV is absent, we not only observed an enhancer chromatin state but also found interactions between the mouse enhancer homolog and Ccnd1 using 3D-genome interaction data. This is similar to the interaction observed in humans and, interestingly, coincides with CTCF binding sites in both species. Taken together, this suggests that MYEOV is a primate-specific gene with a de novo ORF that originated at an evolutionarily older enhancer region. This deeply conserved putative enhancer element could regulate CCND1 in both humans and mice, opening the possibility of studying MYEOV regulatory functions in cancer using non-primate animal models.

1q amplification and PHF19 expressing high-risk cells are associated with relapsed/refractory multiple myeloma

Multiple Myeloma is an incurable plasma cell malignancy with a poor survival rate that is usually treated with immunomodulatory drugs (iMiDs) and proteosome inhibitors (PIs). The malignant plasma cells quickly become resistant to these agents causing relapse and uncontrolled growth of resistant clones. From whole genome sequencing (WGS) and RNA sequencing (RNA-seq) studies, different high-risk translocation, copy number, mutational, and transcriptional markers can be identified. One of these markers, PHF19, epigenetically regulates cell cycle and other processes and is already studied using RNA-seq. In this study, we generate a large (325,025 cells and 49 patients) single cell multi-omic dataset and jointly quantify ATAC- and RNA-seq for each cell and matched genomic profiles for each patient. We identify an association between one plasma cell subtype with myeloma progression that we call relapsed/refractory plasma cells (RRPCs). These cells are associated with chromosome 1q alterations, TP53 mutations, and higher expression of PHF19. We also identify downstream regulation of cell cycle inhibitors in these cells, possible regulation by the transcription factor (TF) PBX1 on chromosome 1q, and determine that PHF19 may be acting primarily through this subset of cells.

Genomic and epigenetic characterization of the arsenic-induced oncogenic microRNA-21

As one of the first identified oncogenic microRNAs, the precise details concerning the transcriptional regulation and function of microRNA-21 (miR-21) are still not completely established. The miR-21 gene is situated on chromosome 17q23.2, positioned at the 3'-UTR of the gene that encodes vacuole membrane protein-1 (VMP1). In this current study, we presented evidence indicating that miR-21 possesses its own gene promoter, which can be found in the intron 10 of the VMP1 gene. Chromatin immunoprecipitation followed by global DNA sequencing (ChIP-seq) revealed the presence of a broad H3K4me3 peak spanning the entire gene body of the primary miR-21 and the existence of super-enhancer clusters in the close proximity to both the miR-21 gene promoter and the transcription termination site in arsenic (As3+)-induced cancer stem-like cells (CSCs) and human induced pluripotent stem cells (hiPSCs). In non-transformed human bronchial epithelial cells (BEAS-2B), As3+ treatment enhanced Nrf2 binding to both the host gene VMP1 of miR-21 and the miR-21 gene. Knockout of Nrf2 inhibited both the basal and As3+-induced expressions of miR-21. Furthermore, the As3+-enhanced Nrf2 peaks in ChIP-seq fully overlap with these super-enhancers enriched with H3K4me1 and H3K27ac in the miR-21 gene, suggesting that Nrf2 may coordinate with other transcription factors through the super-enhancers to regulate the expression of miR-21 in cellular response to As3+. These findings demonstrate the unique genetic and epigenetic characteristics of miR-21 and may provide insights into understanding the novel mechanisms linking environmental As3+ exposure and human cancers.

Targeting super-enhancer activity for colorectal cancer therapy

In addition to genetic variants and copy number alterations, epigenetic deregulation of oncogenes and tumor suppressors is a major contributor in cancer development and propagation. Regulatory elements for gene transcription regulation can be found in promoters which are located in the vicinity of transcription start sites but also at a distance, in enhancer sites, brought to interact with proximal sites when occupied by enhancer protein complexes. These sites provide most of the specific regulatory sequences recognized by transcription factors. A sub-set of enhancers characterized by a longer structure and stronger activity, called super-enhancers, are critical for the expression of specific genes, usually associated with individual cell type identity and function. Super-enhancers show deregulation in cancer, which may have profound repercussions for cancer cell survival and response to therapy. Dysfunction of super-enhancers may result from multiple mechanisms that include changes in their sequence, alterations in the topological neighborhoods where they belong, and alterations in the proteins that mediate their function, such as transcription factors and epigenetic modifiers. These can become potential targets for therapeutic interventions. Genes that are targets of super-enhancers are cell and cancer type specific and could also be of interest for therapeutic targeting. In colorectal cancer, a super-enhancer regulated and over-expressed oncogene is MYC, under the influence of the WNT/β-catenin pathway. Identification and targeting of additional oncogenes regulated by super-enhancers in colorectal cancer may pave the way for combination therapies targeting the super-enhancer machinery and signal transduction pathways that regulate the specific transcription factors operative on them.

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.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

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

Broad domains of H3K4 methylation have been associated with consistent expression of tissue-specific, cell identity, and tumor suppressor genes. Here, we identified broad domain-associated genes in healthy human thymic T cell populations and a collection of T cell acute lymphoblastic leukemia (T-ALL) primary samples and cell lines. We found that broad domains are highly dynamic throughout T cell differentiation, and their varying breadth allows the distinction between normal and neoplastic cells. Although broad domains preferentially associate with cell identity and tumor suppressor genes in normal thymocytes, they flag key oncogenes in T-ALL samples. Moreover, the expression of broad domain-associated genes, both coding and noncoding, is frequently deregulated in T-ALL. Using two distinct leukemic models, we showed that the ectopic expression of T-ALL oncogenic transcription factor preferentially impacts the expression of broad domain-associated genes in preleukemic cells. Finally, an H3K4me3 demethylase inhibitor differentially targets T-ALL cell lines depending on the extent and number of broad domains. Our results show that the regulation of broad H3K4me3 domains is associated with leukemogenesis, and suggest that the presence of these structures might be used for epigenetic prioritization of cancer-relevant genes, including long noncoding RNAs.