Super-enhancers: critical roles and therapeutic targets in hematologic malignancies

Super-enhancer profiling reveals ThPOK/ZBTB7B, a CD4 + cell lineage commitment factor, as a master regulator that restricts breast cancer cells to a luminal non-migratory phenotype

Despite efforts to understand breast cancer biology, metastatic disease remains a clinical challenge. Identifying suppressors of breast cancer progression and mechanisms of transition to more invasive phenotypes could provide game changing therapeutic opportunities. Transcriptional deregulation is central to all malignancies, highlighted by the extensive reprogramming of regulatory elements that underlie oncogenic programs. Among these, super-enhancers (SEs) stand out due to their enrichment in genes controlling cancer hallmarks. To reveal novel breast cancer dependencies, we integrated the analysis of the SE landscape with master regulator activity inference for a series of breast cancer cell lines. As a result, we identified T - h elper-inducing Poxviruses and Zinc-finger ( PO Z)/ K rüppel-like factor (ThPOK, ZBTB7B ), a CD4 + cell lineage commitment factor, as a breast cancer master regulator that is recurrently associated with a SE. ThPOK expression is highest in luminal breast cancer but is significantly reduced in the basal subtype. Manipulation of ThPOK levels in cell lines shows that its repressive function restricts breast cancer cells to an epithelial phenotype by suppressing the expression of genes involved in the epithelial-mesenchymal transition (EMT), WNT/β-catenin target genes, and the pro-metastatic TGFβ pathway. Our study reveals ThPOK as a master transcription factor that restricts the acquisition of metastatic features in breast cancer cells.

Integration of chromosome locations and functional aspects of enhancers and topologically associating domains in knowledge graphs enables versatile queries about gene regulation

Knowledge about transcription factor binding and regulation, target genes, cis-regulatory modules and topologically associating domains is not only defined by functional associations like biological processes or diseases but also has a determinative genome location aspect. Here, we exploit these location and functional aspects together to develop new strategies to enable advanced data querying. Many databases have been developed to provide information about enhancers, but a schema that allows the standardized representation of data, securing interoperability between resources, has been lacking. In this work, we use knowledge graphs for the standardized representation of enhancers and topologically associating domains, together with data about their target genes, transcription factors, location on the human genome, and functional data about diseases and gene ontology annotations. We used this schema to integrate twenty-five enhancer datasets and two domain datasets, creating the most powerful integrative resource in this field to date. The knowledge graphs have been implemented using the Resource Description Framework and integrated within the open-access BioGateway knowledge network, generating a resource that contains an interoperable set of knowledge graphs (enhancers, TADs, genes, proteins, diseases, GO terms, and interactions between domains). We show how advanced queries, which combine functional and location restrictions, can be used to develop new hypotheses about functional aspects of gene expression regulation.

A deep learning model for DNA enhancer prediction based on nucleotide position aware feature encoding

Enhancers, genomic DNA elements, regulate neighboring gene expression crucial for biological processes like cell differentiation and stress response. However, current machine learning methods for predicting DNA enhancers often underutilize hidden features in gene sequences, limiting model accuracy. Hence, this article proposes the PDCNN model, a deep learning-based enhancer prediction method. PDCNN extracts statistical nucleotide representations from gene sequences, discerning positional distribution information of nucleotides in modifier-like DNA sequences. With a convolutional neural network structure, PDCNN employs dual convolutional and fully connected layers. The cross-entropy loss function iteratively updates using a gradient descent algorithm, enhancing prediction accuracy. Model parameters are fine-tuned to select optimal combinations for training, achieving over 95% accuracy. Comparative analysis with traditional methods and existing models demonstrates PDCNN's robust feature extraction capability. It outperforms advanced machine learning methods in identifying DNA enhancers, presenting an effective method with broad implications for genomics, biology, and medical research.

Increased H19/miR-675 Expression in Adult T-Cell Leukemia Is Associated with a Unique Notch Signature Pathway

The Notch pathway is a key cancer driver and is important in tumor progression. Early research suggested that Notch activity was highly dependent on the expression of the intracellular cleaved domain of Notch-1 (NICD). However, recent insights into Notch signaling reveal the presence of Notch pathway signatures, which may vary depending on different cancer types and tumor microenvironments. Herein, we perform a comprehensive investigation of the Notch signaling pathway in adult T-cell leukemia (ATL) primary patient samples. Using gene arrays, we demonstrate that the Notch pathway is constitutively activated in ATL patient samples. Furthermore, the activation of Notch in ATL cells remains elevated irrespective of the presence of activating mutations in Notch itself or its repressor, FBXW7, and that ATL cells are dependent upon Notch-1 expression for proliferation and survival. We demonstrate that ATL cells exhibit the expression of pivotal Notch-related genes, including notch-1, hes1, c-myc, H19, and hes4, thereby defining a critical Notch signature associated with ATL disease. Finally, we demonstrate that lncRNA H19 is highly expressed in ATL patient samples and ATL cells and contributes to Notch signaling activation. Collectively, our results shed further light on the Notch pathway in ATL leukemia and reveal new therapeutic approaches to inhibit Notch activation in ATL cells.

Super-enhancer mediated upregulation of MYEOV suppresses ferroptosis in lung adenocarcinoma

Super-enhancers (SEs) exerted a crucial role in regulating the transcription of oncogenes across various malignancies while the roles of SEs driven genes and the core regulatory elements remain elusive in LUAD. In this study, cancer-specific-SE-genes of lung adenocarcinoma (LUAD) were profiled through H3K27ac ChIP-seq data of cancer cell lines and normal lung tissues, which enriched in in biological processes and pathways integral to the pathophysiology of LUAD. Based on this study, LUAD cells were susceptible to SEs inhibitors, with a reduction of cell proliferation as well as an elevation of apoptosis upon JQ1 or THZ1 intervention. Moreover, the integration of SEs landscapes, CRISPRi, ChIP-PCR, Hi-C data analysis and dual-luciferase reporter assays revealed that myeloma overexpressed gene (MYEOV) was aberrantly overexpressed in LUAD via transcriptional activation by the core SE elements. Functionally, the knockdown of MYEOV undermined cell proliferation in vitro and tumor growth in vivo. In addition, the knockdown of MYEOV generated a prominent ferroptotic phenotype, characterized by elevation of intracellular ferrous iron, reactive oxygen species and lipid peroxidation, together with alteration in marker proteins (SLC7A11, GPX4, FTH1, and ACSL4). Instead, the overexpression of MYEOV accelerated cell proliferation and abrogated ferroptosis. Clinically, the overexpression of MYEOV was observed in LUAD tissues indicating a poor prognosis in patients with LUAD. Mechanistically, SMPD1-induced autophagic degradation of GPX4 assumed a crucial role in the process of ferroptosis triggered by MYEOV knockdown. Serving as an oncogene repressing ferroptosis, promoting proliferation as well as shortening survival in LUAD, SEs-mediated activation of MYEOV might distinguish as a promising therapeutic target.

Multiorgan locked-state model of chronic diseases and systems pharmacology opportunities

With increasing human life expectancy, the global medical burden of chronic diseases is growing. Hence, chronic diseases are a pressing health concern and will continue to be in decades to come. Chronic diseases often involve multiple malfunctioning organs in the body. An imminent question is how interorgan crosstalk contributes to the etiology of chronic diseases. We conceived the locked-state model (LoSM), which illustrates how interorgan communication can give rise to body-wide memory-like properties that 'lock' healthy or pathological conditions. Next, we propose cutting-edge systems biology and artificial intelligence strategies to decipher chronic multiorgan locked states. Finally, we discuss the clinical implications of the LoSM and assess the power of systems-based therapies to dismantle pathological multiorgan locked states while improving treatments for chronic diseases.

MYB: A Key Transcription Factor in the Hematopoietic System Subject to Many Levels of Control

MYB is a master regulator and pioneer factor highly expressed in hematopoietic progenitor cells (HPCs) where it contributes to the reprogramming processes operating during hematopoietic development. MYB plays a complex role being involved in several lineages of the hematopoietic system. At the molecular level, the MYB gene is subject to intricate regulation at many levels through several enhancer and promoter elements, through transcriptional elongation control, as well as post-transcriptional regulation. The protein is modulated by post-translational modifications (PTMs) such as SUMOylation restricting the expression of its downstream targets. Together with a range of interaction partners, cooperating transcription factors (TFs) and epigenetic regulators, MYB orchestrates a fine-tuned symphony of genes expressed during various stages of haematopoiesis. At the same time, the complex MYB system is vulnerable, being a target for unbalanced control and cancer development.

Super-enhancers: Implications in gastric cancer

Gastric cancer (GC) is the fifth most prevalent malignancy and the third leading cause of cancer-related mortality globally. Despite intensive efforts to enhance the efficiencies of various therapeutics (chemotherapy, surgical interventions, molecular-targeted therapies, immunotherapies), the prognosis for patients with GC remains poor. This might be predominantly due to the limited understanding of the complicated etiology of GC. Importantly, epigenetic modifications and alterations are crucial during GC development. Super-enhancers (SEs) are a large cluster of adjacent enhancers that greatly activate transcription. SEs sustain cell-specific identity by enhancing the transcription of specific oncogenes. In this review, we systematically summarize how SEs are involved in GC development, including the SE landscape in GC, the SE target genes in GC, and the interventions related to SE functions for treating GC.

The role of enhancers in psoriasis and atopic dermatitis

Regulatory elements, particularly enhancers, play a crucial role in disease susceptibility and progression. Enhancers are DNA sequences that activate gene expression and can be affected by epigenetic modifications, interactions with transcription factors (TFs) or changes to the enhancer DNA sequence itself. Altered enhancer activity impacts gene expression and contributes to disease. In this review, we define enhancers and the experimental techniques used to identify and characterize them. We also discuss recent studies that examine how enhancers contribute to atopic dermatitis (AD) and psoriasis. Articles in the PubMed database were identified (from 1 January 2010 to 28 February 2023) that were relevant to enhancer variants, enhancer-associated TFs and enhancer histone modifications in psoriasis or AD. Most enhancers associated with these conditions regulate genes affecting epidermal homeostasis or immune function. These discoveries present potential therapeutic targets to complement existing treatment options for AD and psoriasis.

Molecular basis for SOX2-dependent regulation of super-enhancer activity

Pioneer transcription factors (TFs) like SOX2 are vital for stemness and cancer through enhancing gene expression within transcriptional condensates formed with coactivators, RNAs and mediators on super-enhancers (SEs). Despite their importance, how these factors work together for transcriptional condensation and activation remains unclear. SOX2, a pioneer TF found in SEs of pluripotent and cancer stem cells, initiates SE-mediated transcription by binding to nucleosomes, though the mechanism isn't fully understood. To address SOX2's role in SEs, we identified mSE078 as a model SOX2-enriched SE and p300 as a coactivator through bioinformatic analysis. In vitro and cell assays showed SOX2 forms condensates with p300 and SOX2-binding motifs in mSE078. We further proved that SOX2 condensation is highly correlated with mSE078's enhancer activity in cells. Moreover, we successfully demonstrated that p300 not only elevated transcriptional activity but also triggered chromatin acetylation via its direct interaction with SOX2 within these transcriptional condensates. Finally, our validation of SOX2-enriched SEs showcased their contribution to target gene expression in both stem cells and cancer cells. In its entirety, this study imparts valuable mechanistic insights into the collaborative interplay of SOX2 and its coactivator p300, shedding light on the regulation of transcriptional condensation and activation within SOX2-enriched SEs.

MYC function and regulation in physiological perspective

MYC, a key member of the Myc-proto-oncogene family, is a universal transcription amplifier that regulates almost every physiological process in a cell including cell cycle, proliferation, metabolism, differentiation, and apoptosis. MYC interacts with several cofactors, chromatin modifiers, and regulators to direct gene expression. MYC levels are tightly regulated, and deregulation of MYC has been associated with numerous diseases including cancer. Understanding the comprehensive biology of MYC under physiological conditions is an utmost necessity to demark biological functions of MYC from its pathological functions. Here we review the recent advances in biological mechanisms, functions, and regulation of MYC. We also emphasize the role of MYC as a global transcription amplifier.

Chromatin dynamics and subnuclear gene positioning for transcriptional regulation

Plants have been found to exhibit diverse characteristics and functions of chromatin organization, showing both similarities and differences to animals. It is becoming clear how chromatin organization is linked to transcriptional regulation in response to environmental stresses. Regulation of specific chromatin positions in the nuclear space is important for transcription, and the mechanisms that enable such chromatin dynamics are gradually being unveiled. Genes move between subdomains responsible for transcriptional activation or suppression in the subnuclear space in a gene repositioning cycle. We propose a model of localized chromatin interaction in nuclear subdomains, in which the dynamics of local chromatin interactions have a more important impact on the regulation of gene expression than large-scale chromatin organization. In this mini-review, we highlight recent findings on chromatin dynamics, particularly involving transcriptional regulation, and discuss future directions in the study of chromatin organization in plants.

Novel protein products encoded by upstream open reading frames of the MYCN gene in pediatric embryonal tumors

The MYCC and MYCN loci are each associated with two upstream open reading frames (uORFs) potentially encoding small proteins (9-21 kDa). We previously demonstrated that uORFs mrtl and MYCHEX1 of MYCC are translated, and their protein products may function to regulate the expression of the "parent" oncogene. We hypothesized that a similar relationship might exist between MYCN and its two uORFs: MYCNOT and MNOP, and investigated the uORF-encoded proteins associated with MYCN to confirm their expression and intracellular location in neuroblastoma and medulloblastoma cells and tissues. MNOP, MYCNOT, mrtl, and MYCHEX1 were readily detected via reverse transcription polymerase chain reaction and Western blot analysis in tumor cell lines. In tumor tissue, MNOP protein expression was confirmed; however, MCYNOT generated from alternative splicing MYCNΔ1b mRNA was not detected. Immunofluorescence staining of MYCNOT displayed multiple bright foci in the nucleus and diffuse staining in the cytoplasm, suggesting that this small protein may function in both the nucleus and cytoplasm. Upon JQ1 treatment, MYCN, MYCNOT, and mrtl decreased substantially or disappeared completely in three different tumor cell lines. Significant levels of apoptosis were observed in each pediatric embryonal tumor cell line but not T47D breast carcinoma cells, suggesting that response to JQ1 transcriptional inhibition is greatest in tumor cells, which depend on MYC to maintain an undifferentiated phenotype. In conclusion, both MYCN uORF-encoded proteins MNOP and MYCNOT, together with the two MYCC uORF-encoded proteins mrtl and MYCHEX1 were detected simultaneously in tumor cell lines and tumor tissues. These four distinct proteins are translated from the "5'-untranslated region" of MYCN or MYCC mRNA and display consistent distribution patterns within the cell. Additional studies to further elucidate the physiological and pathological roles of these uORF-encoded proteins are warranted, as insights gained could inform new strategies for modulating MYC-family oncogenes by targeting their uORFs.

Super Enhancer Regulatory Gene FYB1 Promotes the Progression of T Cell Acute Lymphoblastic Leukemia by Activating IGLL1

Background

Arising from T progenitor cells, T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignant tumor, accounting for 15% of childhood ALL and 25% of adult ALL. Composing of putative enhancers in close genomic proximity, super enhancer (SE) is critical for cell identity and the pathogenesis of multiple cancers. Belonging to the cytosolute linker protein group, FYB1 is essential for TCR signaling and extensively studied in terms of tumor pathogenesis and metastasis. Dissecting the role of FYN binding protein 1 (FYB1) in T-ALL holds the potential to improve the treatment outcome and prognosis of T-ALL.

Methods

In this study, SEs were explored using public H3K27ac ChIP-seq data derived from T-ALL cell lines, AML cell lines and hematopoietic stem and progenitor cells (HSPCs). Downstream target of FYB1 gene was identified by RNA-seq. Effects of shRNA-mediated downregulation of FYB1 and immunoglobulin lambda-like polypeptide 1 (IGLL1) on self-renewal of T-ALL cells were evaluated in vitro and/or in vivo.

Results

As an SE-driven gene, overexpression of FYB1 was observed in T-ALL, according to the Cancer Cell Line Encyclopedia database. In vitro, knocking down FYB1 led to comprised growth and enhanced apoptosis of T-ALL cells. In vivo, downregulation of FYB1 significantly decreased the disease burden by suppressing tumor growth and improved survival rate. Knocking down FYB1 resulted in significantly decreased expression of IGLL1 that was also an SE-driven gene in T-ALL. As a downstream target of FYB1, IGLL1 exerted similar role as FYB1 in inhibiting growth of T-ALL cells.

Conclusion

Our results suggested that FYB1 gene played important role in regulating self-renewal of T-ALL cells by activating IGLL1, representing a promising therapeutic target for T-ALL patients.

Cytoprotective, Cytotoxic and Cytostatic Roles of Autophagy in Response to BET Inhibitors

The bromodomain and extra-terminal domain (BET) family inhibitors are small molecules that target the dysregulated epigenetic readers, BRD2, BRD3, BRD4 and BRDT, at various transcription-related sites, including super-enhancers. BET inhibitors are currently under investigation both in pre-clinical cell culture and tumor-bearing animal models, as well as in clinical trials. However, as is the case with other chemotherapeutic modalities, the development of resistance is likely to constrain the therapeutic benefits of this strategy. One tumor cell survival mechanism that has been studied for decades is autophagy. Although four different functions of autophagy have been identified in the literature (cytoprotective, cytotoxic, cytostatic and non-protective), primarily the cytoprotective and cytotoxic forms appear to function in different experimental models exposed to BET inhibitors (with some evidence for the cytostatic form). This review provides an overview of the cytoprotective, cytotoxic and cytostatic functions of autophagy in response to BET inhibitors in various tumor models. Our aim is to determine whether autophagy targeting or modulation could represent an effective therapeutic strategy to enhance the response to these modalities and also potentially overcome resistance to BET inhibition.

BRD4 Inhibition as a Strategy to Prolong the Response to Standard of Care in Estrogen Receptor-Positive Breast Cancer

Breast cancer is the most commonly occurring malignancy in women and the second most common cause of cancer-related deaths. ER+ breast cancer constitutes approximately 70% of all breast cancer cases. The standard of care for ER+ breast cancer involves estrogen antagonists such as tamoxifen or fulvestrant in combination with CDK4/6 inhibitors such as palbociclib. However, these treatments are often not curative, with disease recurrence and metastasis being responsible for patient mortality. Overexpression of the epigenetic regulator, BRD4, has been shown to be a negative prognostic indicator in breast cancer, and BET family inhibitors such as ARV-825 and ABBV-744 have garnered interest for their potential to improve and prolong the response to current therapeutic strategies. The current work examined the potential of utilizing ARV-825 and ABBV-744 to increase the effectiveness of tamoxifen or fulvestrant plus palbociclib. ARV-825 was effective in both p53 wild-type (WT) breast tumor cells and in cells lacking functional p53 either alone or in combination with tamoxifen, while the effectiveness of ABBV-744 was limited to fulvestrant plus palbociclib in p53 WT cells. These differential effects may be related to the capacity to suppress c-Myc, a downstream target of BRD4.

Super-enhancers complexes zoom in transcription in cancer

Super-enhancers (SEs) consist of multiple typical enhancers enriched at high density with transcription factors, histone-modifying enzymes and cofactors. Oncogenic SEs promote tumorigenesis and malignancy by altering protein-coding gene expression and noncoding regulatory element function. Therefore, they play central roles in the treatment of cancer. Here, we review the structural characteristics, organization, identification, and functions of SEs and the underlying molecular mechanism by which SEs drive oncogenic transcription in tumor cells. We then summarize abnormal SE complexes, SE-driven coding genes, and noncoding RNAs involved in tumor development. In summary, we believe that SEs show great potential as biomarkers and therapeutic targets.

Super-enhancer-driven lncRNA LIMD1-AS1 activated by CDK7 promotes glioma progression

Long non-coding RNAs (lncRNAs) are tissue-specific expression patterns and dysregulated in cancer. How they are regulated still needs to be determined. We aimed to investigate the functions of glioma-specific lncRNA LIMD1-AS1 activated by super-enhancer (SE) and identify the potential mechanisms. In this paper, we identified a SE-driven lncRNA, LIMD1-AS1, which is expressed at significantly higher levels in glioma than in normal brain tissue. High LIMD1-AS1 levels were significantly associated with a shorter survival time of glioma patients. LIMD1-AS1 overexpression significantly enhanced glioma cells proliferation, colony formation, migration, and invasion, whereas LIMD1-AS1 knockdown inhibited their proliferation, colony formation, migration, and invasion, and the xenograft tumor growth of glioma cells in vivo. Mechanically, inhibition of CDK7 significantly attenuates MED1 recruitment to the super-enhancer of LIMD1-AS1 and then decreases the expression of LIMD1-AS1. Most importantly, LIMD1-AS1 could directly bind to HSPA5, leading to the activation of interferon signaling. Our findings support the idea that CDK7 mediated-epigenetically activation of LIMD1-AS1 plays a crucial role in glioma progression and provides a promising therapeutic approach for patients with glioma.

SENet: A deep learning framework for discriminating super- and typical enhancers by sequence information

Super-enhancers are large domains on the genome where multiple short typical enhancers within a specific genomic distance are stitched together. Typically, they are cell type-specific and responsible for defining cell identity and regulating gene transcription. Numerous studies have demonstrated that super-enhancers are enriched for trait-associated variants, and mutations in super-enhancers are possibly related to known diseases. Recently, several machine learning-based methods have been used to distinguish super-enhancers from typical enhancers by using high-throughput data from various experimental methods. The acquisition of such experimental data is usually costly and time-consuming. In this paper, we innovatively proposed SENet, a groundbreaking method based on a deep neural network model, for discriminating between the two categories solely utilizing sequence information. SENet employs dna2vec feature embedding, convolution for local feature extraction, attention pooling for refined feature retention, and Transformer for contextual information extraction. Experiments demonstrate that SENet outperforms all current state-of-the-art computational methods and shows satisfactory performance in cross-species validation. Our method pioneers the distinction between super-enhancers and typical ones using only sequence information. The source code and datasets are stored in https://github.com/lhy0322/SENet.

Super-enhancer-associated gene CAPG promotes AML progression

Acute myeloid leukemia is the most common acute leukemia in adults, the barrier of refractory and drug resistance has yet to be conquered in the clinical. Abnormal gene expression and epigenetic changes play an important role in pathogenesis and treatment. A super-enhancer is an epigenetic modifier that promotes pro-tumor genes and drug resistance by activating oncogene transcription. Multi-omics integrative analysis identifies the super-enhancer-associated gene CAPG and its high expression level was correlated with poor prognosis in AML. CAPG is a cytoskeleton protein but has an unclear function in AML. Here we show the molecular function of CAPG in regulating NF-κB signaling pathway by proteomic and epigenomic analysis. Knockdown of Capg in the AML murine model resulted in exhausted AML cells and prolonged survival of AML mice. In conclusion, SEs-associated gene CAPG can contributes to AML progression through NF-κB.

Anchoring super-enhancer-driven oncogenic lncRNAs for anti-tumor therapy in hepatocellular carcinoma

Super-enhancer (SE) plays a vital role in the determination of cell identity and fate. Up-regulated expression of coding genes is frequently associated with SE. However, the transcription dysregulation driven by SE, from the viewpoint of long non-coding RNA (lncRNA), remains unclear. Here, SE-associated lncRNAs in HCC are comprehensively outlined for the first time. This study integrally screens and identifies several novel SE-associated lncRNAs that are highly abundant and sensitive to JQ1. Especially, HSAL3 is identified as an uncharacterized SE-driven oncogenic lncRNA, which is activated by transcription factors HCFC1 and HSF1 via its super-enhancer. HSAL3 interference negatively regulates NOTCH signaling, implying the potential mechanism of its tumor-promoting role. The expression of HSAL3 is increased in HCC samples, and higher HSAL3 expression indicates an inferior overall survival of HCC patients. Furthermore, siHSAL3 loaded nanoparticles exert anti-tumor effect on HCC in vitro and in vivo. In conclusion, this is the first comprehensive survey of SE-associated lncRNAs in HCC. HSAL3 is a novel SE-driven oncogenic lncRNA, and siHSAL3 loaded nanoparticles are therapeutic candidates for HCC. This work sheds lights on the merit of anchoring SE-driven oncogenic lncRNAs for HCC treatment.

DNA methylation differences in noncoding regions in ER negative breast tumors between Black and White women

Introduction

Incidence of estrogen receptor (ER)-negative breast cancer, an aggressive tumor subtype associated with worse prognosis, is higher among African American/Black women than other US racial and ethnic groups. The reasons for this disparity remain poorly understood but may be partially explained by differences in the epigenetic landscape.

Methods

We previously conducted genome-wide DNA methylation profiling of ER- breast tumors from Black and White women and identified a large number of differentially methylated loci (DML) by race. Our initial analysis focused on DML mapping to protein-coding genes. In this study, motivated by increasing appreciation for the biological importance of the non-protein coding genome, we focused on 96 DMLs mapping to intergenic and noncoding RNA regions, using paired Illumina Infinium Human Methylation 450K array and RNA-seq data to assess the relationship between CpG methylation and RNA expression of genes located up to 1Mb away from the CpG site.

Results

Twenty-three (23) DMLs were significantly correlated with the expression of 36 genes (FDR<0.05), with some DMLs associated with the expression of single gene and others associated with more than one gene. One DML (cg20401567), hypermethylated in ER- tumors from Black versus White women, mapped to a putative enhancer/super-enhancer element located 1.3 Kb downstream of HOXB2. Increased methylation at this CpG correlated with decreased expression of HOXB2 (Rho=-0.74, FDR<0.001) and other HOXB/HOXB-AS genes. Analysis of an independent set of 207 ER- breast cancers from TCGA similarly confirmed hypermethylation at cg20401567 and reduced HOXB2 expression in tumors from Black versus White women (Rho=-0.75, FDR<0.001).

Discussion

Our findings indicate that epigenetic differences in ER- tumors between Black and White women are linked to altered gene expression and may hold functional significance in breast cancer pathogenesis.

New genetic and epigenetic insights into the chemokine system: the latest discoveries aiding progression toward precision medicine

Over the past thirty years, the importance of chemokines and their seven-transmembrane G protein-coupled receptors (GPCRs) has been increasingly recognized. Chemokine interactions with receptors trigger signaling pathway activity to form a network fundamental to diverse immune processes, including host homeostasis and responses to disease. Genetic and nongenetic regulation of both the expression and structure of chemokines and receptors conveys chemokine functional heterogeneity. Imbalances and defects in the system contribute to the pathogenesis of a variety of diseases, including cancer, immune and inflammatory diseases, and metabolic and neurological disorders, which render the system a focus of studies aiming to discover therapies and important biomarkers. The integrated view of chemokine biology underpinning divergence and plasticity has provided insights into immune dysfunction in disease states, including, among others, coronavirus disease 2019 (COVID-19). In this review, by reporting the latest advances in chemokine biology and results from analyses of a plethora of sequencing-based datasets, we outline recent advances in the understanding of the genetic variations and nongenetic heterogeneity of chemokines and receptors and provide an updated view of their contribution to the pathophysiological network, focusing on chemokine-mediated inflammation and cancer. Clarification of the molecular basis of dynamic chemokine-receptor interactions will help advance the understanding of chemokine biology to achieve precision medicine application in the clinic.

Genomic targets of the IRE1-XBP1s pathway in mediating metabolic adaptation in epithelial plasticity

Epithelial mesenchymal plasticity (EMP) is a complex cellular reprogramming event that plays a major role in tissue homeostasis. Recently we observed the unfolded protein response (UPR) triggers EMP through the inositol-requiring protein 1 (IRE1α)-X-box-binding protein 1 spliced (XBP1s) axis, enhancing glucose shunting to protein N glycosylation. To better understand the genomic targets of XBP1s, we identified its genomic targets using Cleavage Under Targets and Release Using Nuclease (CUT&RUN) of a FLAG-epitope tagged XBP1s in RSV infection. CUT&RUN identified 7086 binding sites in chromatin that were enriched in AP-1 motifs and GC-sequences. Of these binding sites, XBP1s peaks mapped to 4827 genes controlling Rho-GTPase signaling, N-linked glycosylation and ER-Golgi transport. Strikingly, XBP1s peaks were within 1 kb of transcription start sites of 2119 promoters. In addition to binding core mesenchymal transcription factors SNAI1 and ZEB1, we observed that hexosamine biosynthetic pathway (HBP) enzymes were induced and contained proximal XBP1s peaks. We demonstrate that IRE1α -XBP1s signaling is necessary and sufficient to activate core enzymes by recruiting elongation-competent phospho-Ser2 CTD modified RNA Pol II. We conclude that the IRE1α-XBP1s pathway coordinately regulates mesenchymal transcription factors and hexosamine biosynthesis in EMP by a mechanism involving recruitment of activated pSer2-Pol II to GC-rich promoters.

Tumorigenic activation around HPV integrated sites in head and neck squamous cell carcinoma

Human papillomavirus (HPV) is causally involved in the development of head and neck squamous cell carcinoma (HNSCC). The integration of HPV drives tumorigenesis through expression of oncogenic viral genes as well as genomic alterations in surrounding regions. To elucidate involvement of epigenetic dysregulation in tumorigenesis, we here performed integrated analyses of the epigenome, transcriptome and interactome using ChIP-seq, RNA-seq and Hi-C and 4C-seq for HPV(+) HNSCCs. We analyzed clinical HNSCC using The Cancer Genome Atlas data and found that genes neighboring HPV integration sites were significantly upregulated and were correlated with oncogenic phenotypes in HPV(+) HNSCCs. While we found four HPV integration sites in HPV(+) HNSCC cell line UPCI-SCC-090 through target enrichment sequencing, 4C-seq revealed 0.5 to 40 Mb of HPV-interacting regions (HPVIRs) where host genomic regions interacted with integrated HPV genomes. While 9% of the HPVIRs were amplified and activated epigenetically forming super-enhancers, the remaining non-amplified regions were found to show a significant increase in H3K27ac levels and an upregulation of genes associated with GO terms, for example, Signaling by WNT and Cell Cycle. Among those genes, ITPR3 was significantly upregulated, involving UPCI-SCC-090-specific super-enhancer formation around the ITPR3 promoter and in the 80-kb-downstream region. The knockdown of ITPR3 by siRNA or CRISPR deletions of the distant enhancer region led to a significant suppression of cell proliferation. The epigenetic activation of HPVIRs was also confirmed in other cell lines, UM-SCC-47 and UM-SCC-104. These data indicate that epigenetic activation in HPVIRs contributes, at least partially, to genesis of HPV(+) HNSCC.

Super-enhancer-driven TOX2 mediates oncogenesis in Natural Killer/T Cell Lymphoma

BACKGROUND

Extranodal natural killer/T-cell lymphoma (NKTL) is an aggressive type of non-Hodgkin lymphoma with dismal outcome. A better understanding of disease biology and key oncogenic process is necessary for the development of targeted therapy. Super-enhancers (SEs) have been shown to drive pivotal oncogenes in various malignancies. However, the landscape of SEs and SE-associated oncogenes remain elusive in NKTL.

METHODS

We used Nano-ChIP-seq of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs NKTL primary tumor samples. Integrative analysis of RNA-seq and survival data further pinned down high value, novel SE oncogenes. We utilized shRNA knockdown, CRISPR-dCas9, luciferase reporter assay, ChIP-PCR to investigate the regulation of transcription factor (TF) on SE oncogenes. Multi-color immunofluorescence (mIF) staining was performed on an independent cohort of clinical samples. Various function experiments were performed to evaluate the effects of TOX2 on the malignancy of NKTL in vitro and in vivo.

RESULTS

SE landscape was substantially different in NKTL samples in comparison with normal tonsils. Several SEs at key transcriptional factor (TF) genes, including TOX2, TBX21(T-bet), EOMES, RUNX2, and ID2, were identified. We confirmed that TOX2 was aberrantly overexpressed in NKTL relative to normal NK cells and high expression of TOX2 was associated with worse survival. Modulation of TOX2 expression by shRNA, CRISPR-dCas9 interference of SE function impacted on cell proliferation, survival and colony formation ability of NKTL cells. Mechanistically, we found that RUNX3 regulates TOX2 transcription by binding to the active elements of its SE. Silencing TOX2 also impaired tumor formation of NKTL cells in vivo. Metastasis-associated phosphatase PRL-3 has been identified and validated as a key downstream effector of TOX2-mediated oncogenesis.

CONCLUSIONS

Our integrative SE profiling strategy revealed the landscape of SEs, novel targets and insights into molecular pathogenesis of NKTL. The RUNX3-TOX2-SE-TOX2-PRL-3 regulatory pathway may represent a hallmark of NKTL biology. Targeting TOX2 could be a valuable therapeutic intervene for NKTL patients and warrants further study in clinic.

Pan-cancer analysis of super enhancer-induced PRR7-AS1 as a potential prognostic and immunological biomarker

Introduction: Systematic pan-cancer analysis of the roles and regulatory mechanisms for PRR7-AS1 is currently not available. Methods: In the present study, a comprehensive bioinformatic approach was used to mine the underlying oncogenic effects of PRR7-AS1, including expression status, prognostic value and immune characteristics. Results: We discovered that PRR7-AS1 expression was remarkably upregulated in most cancer types and exhibited a negative correlation with the prognosis. Furthermore, PRR7-AS1 expression was inversely connected with the majority of tumor-infiltrating immune cells, immune scores and immune checkpoint gene expression in pancancer. There was also a significant correlation between PRR7-AS1 expression status and tumor mutational burden, microsatellite instability, and neoantigens in certain tumors. PRR7-AS1 had the best predictive power for immune checkpoint blockade efficacy compared to other well-recognized biomarkers. PRR7-AS1 overexpression could affect cytotoxic T cells-mediated antitumor responses. Functional enrichment analysis revealed that PRR7-AS1 might be involved in the metabolic pathways. Super enhancer activity might have participated in the regulation of PRR7-AS1 expression. And we constructed the competitive endogenous RNA networks for PRR7-AS1. Discussion: In general, PRR7-AS1 had the potential to be a diagnostic, prognostic and immune biomarker for pan cancer. PRR7-AS1 was correlated with an immunosuppressive microenvironment and was a new potential target for immunotherapy. Epigenetic factors were the driving forces for PRR7-AS1 overexpression in tumors.

CDK9 inhibition induces epigenetic reprogramming revealing strategies to circumvent resistance in lymphoma

Diffuse large B-cell lymphoma (DLBCL) exhibits significant genetic heterogeneity which contributes to drug resistance, necessitating development of novel therapeutic approaches. Pharmacological inhibitors of cyclin-dependent kinases (CDK) demonstrated pre-clinical activity in DLBCL, however many stalled in clinical development. Here we show that AZD4573, a selective inhibitor of CDK9, restricted growth of DLBCL cells. CDK9 inhibition (CDK9i) resulted in rapid changes in the transcriptome and proteome, with downmodulation of multiple oncoproteins (eg, MYC, Mcl-1, JunB, PIM3) and deregulation of phosphoinotiside-3 kinase (PI3K) and senescence pathways. Following initial transcriptional repression due to RNAPII pausing, we observed transcriptional recovery of several oncogenes, including MYC and PIM3. ATAC-Seq and ChIP-Seq experiments revealed that CDK9i induced epigenetic remodeling with bi-directional changes in chromatin accessibility, suppressed promoter activation and led to sustained reprograming of the super-enhancer landscape. A CRISPR library screen suggested that SE-associated genes in the Mediator complex, as well as AKT1, confer resistance to CDK9i. Consistent with this, sgRNA-mediated knockout of MED12 sensitized cells to CDK9i. Informed by our mechanistic findings, we combined AZD4573 with either PIM kinase or PI3K inhibitors. Both combinations decreased proliferation and induced apoptosis in DLBCL and primary lymphoma cells in vitro as well as resulted in delayed tumor progression and extended survival of mice xenografted with DLBCL in vivo. Thus, CDK9i induces reprogramming of the epigenetic landscape, and super-enhancer driven recovery of select oncogenes may contribute to resistance to CDK9i. PIM and PI3K represent potential targets to circumvent resistance to CDK9i in the heterogeneous landscape of DLBCL.

Somatic Super-Enhancer Epigenetic Signature for Overall Survival Prediction in Patients with Breast Invasive Carcinoma

To analyze genome-wide super-enhancers (SEs) methylation signature of breast invasive carcinoma (BRCA) and its clinical value. Differential methylation sites (DMS) between BRCA and adjacent tissues from The Cancer Genome Atlas (TCGA) database were identified by using ChAMP package in R software. Super-enhancers were identified sing ROSE software. Overlap analysis was used to assess the potential DMS in SEs region. Feature selection was performed by Cox regression and least absolute shrinkage and selection operator (LASSO) algorithm based on TCGA training cohort. Prognosis model validation was performed in TCGA training cohort, TCGA validation cohort, and gene expression omnibus (GEO) test cohort. The gene ontology and KEGG analysis revealed that SEs target genes were significantly enriched in cell-migration-associated processes and pathways. A total of 83 654 DMS were identified between BRCA and adjacent tissues. Around 2397 DMS in SEs region were identified by overlap study and used to feature selection. By using Cox regression and LASSO algorithm, 42 features were selected to develop a clinical prediction model (CPM). Both training (TCGA) and validation cohorts (TCGA and GEO) show that the CPM has ideal discrimination and calibration. The CPM based on DMS at SE regions has ideal discrimination and calibration, which combined with tumor node metastasis (TNM) stage could improve prognostication, and thus contribute to individualized medicine.

Bioinformatics Pipelines for Identification of Super-Enhancers and 3D Chromatin Contacts

Precise regulation of gene expression is integral in development. Emerging studies have highlighted that super-enhancers (SEs), which are clusters of multiple enhancers, play critical roles in regulating cell type-specific gene expression via 3D chromatin, thereby defining the cellular identities of given cells. Here we provide optimized bioinformatics pipelines to identify SEs and 3D chromatin contacts. Our pipelines encompass the processing of chromatin immunoprecipitation sequencing (ChIP-seq) data to identify SEs and the processing of genome-wide chromosome conformation capture (Hi-C) data. We can then infer long-range chromatin contacts between SEs and other genomic regions. This integrative computational approach, which can be applied to CUT&RUN and CUT&Tag, alternative technologies to ChIP-seq, allows us to identify genomic locations of SEs and their 3D genome configuration, whereby multiple SEs act in concert. We show an analysis of mouse spermatogenesis as an example of this application.

Targeting cereblon in hematologic malignancies

The protein cereblon (CRBN) is a substrate receptor of the cullin 4-really interesting new gene (RING) E3 ubiquitin ligase complex CRL4^CRBN. Targeting CRBN mediates selective protein ubiquitination and subsequent degradation via the proteasome. This review describes novel thalidomide analogs, immunomodulatory drugs, also known as CRBN E3 ubiquitin ligase modulators or molecular glues (avadomide, iberdomide, CC-885, CC-90009, BTX-1188, CC-92480, CC-99282, CFT7455, and CC-91633), and CRBN-based proteolysis targeting chimeras (PROTACs) with increased efficacy and potent activity for application in hematologic malignancies. Both types of CRBN-binding drugs, molecular glues, and PROTACs stimulate the interaction between CRBN and its neosubstrates, recruiting target disease-promoting proteins and the E3 ubiquitin ligase CRL4^CRBN. Proteins that are traditionally difficult to target (transcription factors and oncoproteins) can be polyubiquitinated and degraded in this way. The competition of CRBN neosubstrates with endogenous CRBN-interacting proteins and the pharmacology and rational combination therapies of and mechanisms of resistance to CRL4^CRBN modulators or CRBN-based PROTACs are described.

Reversible Myc hypomorphism identifies a key Myc-dependency in early cancer evolution

Germ-line hypomorphism of the pleiotropic transcription factor Myc in mice, either through Myc gene haploinsufficiency or deletion of Myc enhancers, delays onset of various cancers while mice remain viable and exhibit only relatively mild pathologies. Using a genetically engineered mouse model in which Myc expression may be systemically and reversibly hypomorphed at will, we asked whether this resistance to tumour progression is also emplaced when Myc hypomorphism is acutely imposed in adult mice. Indeed, adult Myc hypomorphism profoundly blocked KRasG12D-driven lung and pancreatic cancers, arresting their evolution at the early transition from indolent pre-tumour to invasive cancer. We show that such arrest is due to the incapacity of hypomorphic levels of Myc to drive release of signals that instruct the microenvironmental remodelling necessary to support invasive cancer. The cancer protection afforded by long-term adult imposition of Myc hypomorphism is accompanied by only mild collateral side effects, principally in haematopoiesis, but even these are circumvented if Myc hypomorphism is imposed metronomically whereas potent cancer protection is retained.

Non-coding genome in small cell lung cancer between theoretical view and clinical applications

Small cell lung cancer (SCLC) is a highly aggressive cancer of the neuroendocrine system, characterized by poor differentiation, rapid growth, and poor overall survival (OS) of patients. Despite the recent advances in the treatment of SCLC recently, the 2-year survival rate of patients with the cancer is only 14-15%, occasioned by the acquired resistance to drugs and serious off-target effects. In humans, the coding region is only 2% of the total genome, and 20% of that is associated with human diseases. Beyond the coding genome are RNAs, promoters, enhancers, and other intricate elements. The non-coding regulatory regions, mainly the non-coding RNAs (ncRNAs), regulate numerous biological activities including cell proliferation, metastasis, and drug resistance. As such, they are potential diagnostic or prognostic biomarkers, and also potential therapeutic targets for SCLC. Therefore, understanding how non-coding elements regulate SCLC development and progression holds significant clinical implications. Herein, we summarized the recent discoveries on the relationship between the non-coding elements including long non-coding RNAs (lncRNA), microRNAs (miRNAs), circular RNA (circRNA), enhancers as well as promotors, and the pathogenesis of SCLC and their potential clinical applications.

Super enhancers: Pathogenic roles and potential therapeutic targets for acute myeloid leukemia (AML)

Acute myeloid leukemia (AML) is a malignant hematological tumor with disordered oncogenes/tumor suppressor genes and limited treatments. The potent anti-cancer effects of bromodomain and extra-terminal domain (BET) inhibitors, targeting the key component of super enhancers, in early clinical trials on AML patients, implies the critical role of super enhancers in AML. Here, we review the concept and characteristic of super enhancer, and then summarize the current researches about super enhancers in AML pathogenesis, diagnosis and classification, followed by illustrate the potential super enhancer-related targets and drugs, and propose the future directions of super enhancers in AML. This information provides integrated insight into the roles of super enhancers in this disease.

BET-Independent Murine Leukemia Virus Integration Is Retargeted In Vivo and Selects Distinct Genomic Elements for Lymphomagenesis

Moloney murine leukemia virus (MLV) infects BALB/c mice and induces T-cell lymphoma in mice. Retroviral integration is mediated by the interaction of the MLV integrase (IN) with members of the bromodomain and extraterminal motif (BET) protein family (BRD2, BRD3, and BRD4). The introduction of the W390A mutation into MLV IN abolishes the BET interaction. Here, we compared the replication of W390A MLV to that of wild-type (WT) MLV in adult BALB/c mice to study the role of BET proteins in replication, integration, and tumorigenesis in vivo. Comparing WT and W390A MLV infections revealed similar viral loads in the blood, thymus, and spleen cells. Interestingly, W390A MLV integration was retargeted away from GC-enriched genomic regions. However, both WT MLV- and W390A MLV-infected mice developed T-cell lymphoma after similar latencies represented by an enlarged thymus and spleen and multiorgan tumor infiltration. Integration site sequencing from splenic tumor cells revealed clonal expansion in all WT MLV- and W390A MLV-infected mice. However, the integration profiles of W390A MLV and WT MLV differed significantly. Integrations were enriched in enhancers and promoters, but compared to the WT, W390A MLV integrated less frequently into enhancers and more frequently into oncogene bodies such as Notch1 and Ppp1r16b. We conclude that host factors direct MLV in vivo integration site selection. Although BET proteins target WT MLV integration preferentially toward enhancers and promoters, insertional lymphomagenesis can occur independently from BET, likely due to the intrinsically strong enhancer/promoter of the MLV long terminal repeat (LTR).

IMPORTANCE In this study, we have shown that the in vivo replication of murine leukemia virus happens independently of BET proteins, which are key host determinants involved in retroviral integration site selection. This finding opens a new research line in the discovery of alternative viral or host factors that may complement the dominant host factor. In addition, our results show that BET-independent murine leukemia virus uncouples insertional mutagenesis from gene enhancers, although lymphomagenesis still occurs despite the lack of an interaction with BET proteins. Our findings also have implications for the engineering of BET-independent MLV-based vectors for gene therapy, which may not be a safe alternative.

Chromatin Rewiring by Mismatch Repair Protein MSH2 Alters Cell Adhesion Pathways and Sensitivity to BET Inhibition in Gastric Cancer

Mutations in the DNA mismatch repair gene MSH2 are causative of microsatellite instability (MSI) in multiple cancers. Here, we discovered that besides its well-established role in DNA repair, MSH2 exerts a novel epigenomic function in gastric cancer. Unbiased CRISPR-based mass spectrometry combined with genome-wide CRISPR functional screening revealed that in early-stage gastric cancer MSH2 genomic binding is not randomly distributed but rather is associated specifically with tumor-associated super-enhancers controlling the expression of cell adhesion genes. At these loci, MSH2 genomic binding was required for chromatin rewiring, de novo enhancer-promoter interactions, maintenance of histone acetylation levels, and regulation of cell adhesion pathway expression. The chromatin function of MSH2 was independent of its DNA repair catalytic activity but required MSH6, another DNA repair gene, and recruitment to gene loci by the SWI/SNF chromatin remodeler SMARCA4/BRG1. Loss of MSH2 in advanced gastric cancers was accompanied by deficient cell adhesion pathway expression, epithelial-mesenchymal transition, and enhanced tumorigenesis in vitro and in vivo. However, MSH2-deficient gastric cancers also displayed addiction to BAZ1B, a bromodomain-containing family member, and consequent synthetic lethality to bromodomain and extraterminal motif (BET) inhibition. Our results reveal a role for MSH2 in gastric cancer epigenomic regulation and identify BET inhibition as a potential therapy in MSH2-deficient gastric malignancies.

SIGNIFICANCE

DNA repair protein MSH2 binds and regulates cell adhesion genes by enabling enhancer-promoter interactions, and loss of MSH2 causes deficient cell adhesion and bromodomain and extraterminal motif inhibitor synthetic lethality in gastric cancer.

Super-Enhancers, Phase-Separated Condensates, and 3D Genome Organization in Cancer

3D chromatin organization plays an important role in transcription regulation and gene expression. The 3D genome is highly maintained by several architectural proteins, such as CTCF, Yin Yang 1, and cohesin complex. This structural organization brings regulatory DNA elements in close proximity to their target promoters. In this review, we discuss the 3D chromatin organization of super-enhancers and their relationship to phase-separated condensates. Super-enhancers are large clusters of DNA elements. They can physically contact with their target promoters by chromatin looping during transcription. Multiple transcription factors can bind to enhancer and promoter sequences and recruit a complex array of transcriptional co-activators and RNA polymerase II to effect transcriptional activation. Phase-separated condensates of transcription factors and transcriptional co-activators have been implicated in assembling the transcription machinery at particular enhancers. Cancer cells can hijack super-enhancers to drive oncogenic transcription to promote cell survival and proliferation. These dysregulated transcriptional programs can cause cancer cells to become highly dependent on transcriptional regulators, such as Mediator and BRD4. Moreover, the expression of oncogenes that are driven by super-enhancers is sensitive to transcriptional perturbation and often occurs in phase-separated condensates, supporting therapeutic rationales of targeting SE components, 3D genome organization, or dysregulated condensates in cancer.

Manipulation and Mechanical Deformation of Leukemia Cells by High-Frequency Ultrasound Single Beam

Ultrasound single-beam acoustic tweezer system has attracted increasing attention in the field of biomechanics. Cell biomechanics play a pivotal role in leukemia cell functions. To better understand and compare the cell mechanics of the leukemia cells, herein, we fabricated an acoustic tweezer system in-house connected with a 50-MHz high-frequency cylinder ultrasound transducer. Selected leukemia cells (Jurkat, K562, and MV-411 cells) were cultured, trapped, and manipulated by high-frequency ultrasound single beam, which was transmitted from the ultrasound transducer without contacting any cells. The relative deformability of each leukemia cell was measured, characterized, and compared, and the leukemia cell (Jurkat cell) gaining the highest deformability was highlighted. Our results demonstrate that the high-frequency ultrasound single beam can be utilized to manipulate and characterize leukemia cells, which can be applied to study potential mechanisms in the immune system and cell biomechanics in other cell types.

The molecular understanding of super-enhancer dysregulation in cancer

Abnormalities in the regulation of gene expression are associated with various pathological conditions. Among the distal regulatory elements in the genome, the activation of target genes by enhancers plays a central role in the formation of cell type-specific gene expression patterns. Super-enhancers are a subclass of enhancers that frequently contain multiple enhancer-like elements and are characterized by dense binding of master transcription factors and Mediator complexes and high signals of active histone marks. Super-enhancers have been studied in detail as important regulatory regions that control cell identity and contribute to the pathogenesis of diverse diseases. In cancer, super-enhancers have multifaceted roles by activating various oncogenes and other cancer-related genes and shaping characteristic gene expression patterns in cancer cells. Alterations in super-enhancer activities in cancer involve multiple mechanisms, including the dysregulation of transcription factors and the super-enhancer-associated genomic abnormalities. The study of super-enhancers could contribute to the identification of effective biomarkers and the development of cancer therapeutics targeting transcriptional addiction. In this review, we summarize the roles of super-enhancers in cancer biology, with a particular focus on hematopoietic malignancies, in which multiple super-enhancer alteration mechanisms have been reported.

Super enhancers as master gene regulators in the pathogenesis of hematologic malignancies

Transcriptional deregulation of multiple oncogenes, tumor suppressors and survival pathways is a cancer cell hallmark. Super enhancers (SE) are long stretches of active enhancers in close linear proximity that ensure extraordinarily high expression levels of key genes associated with cell lineage, function and survival. SE landscape is intrinsically prone to changes and reorganization during the course of normal cell differentiation. This functional plasticity is typically utilized by cancer cells, which remodel their SE landscapes to ensure oncogenic transcriptional reprogramming. Multiple recent studies highlighted structural genetic mechanisms in non-coding regions that create new SE or hijack already existing ones. In addition, alterations in abundance/activity of certain SE-associated proteins or certain viral infections can elicit new super enhancers and trigger SE-driven transcriptional changes. For these reasons, SE profiling emerged as a powerful tool for discovering the core transcriptional regulatory circuits in tumor cells. This, in turn, provides new insights into cancer cell biology, and identifies main nodes of key cellular pathways to be potentially targeted. Since SEs are susceptible to inhibition, their disruption results in exponentially amassing 'butterfly' effect on gene expression and cell function. Moreover, many of SE elements are druggable, opening new therapeutic opportunities. Indeed, SE targeting drugs have been studied preclinically in various hematologic malignancies with promising effects. Herein, we review the unique features of SEs, present different cis- and trans-acting mechanisms through which hematologic tumor cells acquire SEs, and finally, discuss the potential of SE targeting in the therapy of hematologic malignancies.

Super Enhancer-Mediated Upregulation of HJURP Promotes Growth and Survival of t(4;14)-Positive Multiple Myeloma

Multiple myeloma is an incurable malignancy with marked clinical and genetic heterogeneity. The cytogenetic abnormality t(4;14) (p16.3;q32.3) confers aggressive behavior in multiple myeloma. Recently, essential oncogenic drivers in a wide range of cancers have been shown to be controlled by super-enhancers (SE). We used chromatin immunoprecipitation sequencing of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs in t(4;14)-translocated multiple myeloma. The histone chaperone HJURP was aberrantly overexpressed in t(4;14)-positive multiple myeloma due to transcriptional activation by a distal SE induced by the histone lysine methyltransferase NSD2. Silencing of HJURP with short hairpin RNA or CRISPR interference of SE function impaired cell viability and led to apoptosis. Conversely, HJURP overexpression promoted cell proliferation and abrogated apoptosis. Mechanistically, the NSD2/BRD4 complex positively coregulated HJURP transcription by binding the promoter and active elements of its SE. In summary, this study introduces SE profiling as an efficient approach to identify new targets and understand molecular pathogenesis in specific subtypes of cancer. Moreover, HJURP could be a valuable therapeutic target in patients with t(4;14)-positive myeloma. SIGNIFICANCE: A super-enhancer screen in t(4;14) multiple myeloma serves to identify genes that promote growth and survival of myeloma cells, which may be evaluated in future studies as therapeutic targets.

Molecular Mechanisms of lncRNAs in the Dependent Regulation of Cancer and Their Potential Therapeutic Use

Deep whole genome and transcriptome sequencing have highlighted the importance of an emerging class of non-coding RNA longer than 200 nucleotides (i.e., long non-coding RNAs (lncRNAs)) that are involved in multiple cellular processes such as cell differentiation, embryonic development, and tissue homeostasis. Cancer is a prime example derived from a loss of homeostasis, primarily caused by genetic alterations both in the genomic and epigenetic landscape, which results in deregulation of the gene networks. Deregulation of the expression of many lncRNAs in samples, tissues or patients has been pointed out as a molecular regulator in carcinogenesis, with them acting as oncogenes or tumor suppressor genes. Herein, we summarize the distinct molecular regulatory mechanisms described in literature in which lncRNAs modulate carcinogenesis, emphasizing epigenetic and genetic alterations in particular. Furthermore, we also reviewed the current strategies used to block lncRNA oncogenic functions and their usefulness as potential therapeutic targets in several carcinomas.

OUP accepted manuscript

Super enhancers (SEs) are broad enhancer domains usually containing multiple constituent enhancers that hold elevated activities in gene regulation. Disruption in one or more constituent enhancers causes aberrant SE activities that lead to gene dysregulation in diseases. To quantify SE aberrations, differential analysis is performed to compare SE activities between cell conditions. The state-of-art strategy in estimating differential SEs relies on overall activities and neglect the changes in length and structure of SEs. Here, we propose a novel computational method to identify differential SEs by weighting the combinatorial effects of constituent-enhancer activities and locations (i.e. internal dynamics). In addition to overall activity changes, our method identified four novel classes of differential SEs with distinct enhancer structural alterations. We demonstrate that these structure alterations hold distinct regulatory impact, such as regulating different number of genes and modulating gene expression with different strengths, highlighting the differentiated regulatory roles of these unexplored SE features. When compared to the existing method, our method showed improved identification of differential SEs that were linked to better discernment of cell-type-specific SE activity and functional interpretation.

Runt-Related Transcription Factor 3 Promotes Acute Myeloid Leukemia Progression

Acute myeloid leukemia (AML) is an aggressive hematological malignancy with high relapse/refractory rate. Genetic and epigenetic abnormalities are driving factors for leukemogenesis. RUNX1 and RUNX2 from the Runt-related transcription factor (RUNX) family played important roles in AML pathogenesis. However, the relationship between RUNX3 and AML remains unclear. Here, we found that RUNX3 was a super-enhancer-associated gene and highly expressed in AML cells. The Cancer Genome Atlas (TCGA) database showed high expression of RUNX3 correlated with poor prognosis of AML patients. We observed that Runx3 knockdown significantly inhibited leukemia progression by inducing DNA damage to enhance apoptosis in murine AML cells. By chromatin immunoprecipitation sequencing (ChIP-seq) analysis, we discovered that RUNX3 in AML cells mainly bound more genes involved in DNA-damage repair and antiapoptosis pathways compared to that in normal bone marrow cells. Runx3 knockdown obviously inhibited the expression of these genes in AML cells. Overall, we identified RUNX3 as an oncogene overexpressed in AML cells, and Runx3 knockdown suppressed AML progression by inducing DNA damage and apoptosis.

The non-coding genome in genetic brain disorders: new targets for therapy?

The non-coding genome, consisting of more than 98% of all genetic information in humans and once judged as ‘Junk DNA’, is increasingly moving into the spotlight in the field of human genetics. Non-coding regulatory elements (NCREs) are crucial to ensure correct spatio-temporal gene expression. Technological advancements have allowed to identify NCREs on a large scale, and mechanistic studies have helped to understand the biological mechanisms underlying their function. It is increasingly becoming clear that genetic alterations of NCREs can cause genetic disorders, including brain diseases. In this review, we concisely discuss mechanisms of gene regulation and how to investigate them, and give examples of non-coding alterations of NCREs that give rise to human brain disorders. The cross-talk between basic and clinical studies enhances the understanding of normal and pathological function of NCREs, allowing better interpretation of already existing and novel data. Improved functional annotation of NCREs will not only benefit diagnostics for patients, but might also lead to novel areas of investigations for targeted therapies, applicable to a wide panel of genetic disorders. The intrinsic complexity and precision of the gene regulation process can be turned to the advantage of highly specific treatments. We further discuss this exciting new field of ‘enhancer therapy’ based on recent examples.

Super-enhancer-associated TMEM44-AS1 aggravated glioma progression by forming a positive feedback loop with Myc

Background

Long non-coding RNAs (lncRNAs) have been considered as one type of gene expression regulator for cancer development, but it is not clear how these are regulated. This study aimed to identify a specific lncRNA that promotes glioma progression.

Methods

RNA sequencing (RNA-seq) and quantitative real-time PCR were performed to screen differentially expressed genes. CCK-8, transwell migration, invasion assays, and a mouse xenograft model were performed to determine the functions of TMEM44-AS1. Co-IP, ChIP, Dual-luciferase reporter assays, RNA pulldown, and RNA immunoprecipitation assays were performed to study the molecular mechanism of TMEM44-AS1 and the downstream target.

Results

We identified a novel lncRNA TMEM44-AS1, which was aberrantly expressed in glioma tissues, and that increased TMEM44-AS1 expression was correlated with malignant progression and poor survival for patients with glioma. Expression of TMEM44-AS1 increased the proliferation, colony formation, migration, and invasion of glioma cells. Knockdown of TMEM44-AS1 in glioma cells reduced cell proliferation, colony formation, migration and invasion, and tumor growth in a nude mouse xenograft model. Mechanistically, TMEM44-AS1 is directly bound to the SerpinB3, and sequentially activated Myc and EGR1/IL-6 signaling; Myc transcriptionally induced TMEM44-AS1 and directly bound to the promoter and super-enhancer of TMEM44-AS1, thus forming a positive feedback loop with TMEM44-AS. Further studies demonstrated that Myc interacts with MED1 regulates the super-enhancer of TMEM44-AS1. More importantly, a novel small-molecule Myc inhibitor, Myci975, alleviated TMEM44-AS1-promoted the growth of glioma cells.

Conclusions

Our study implicates a crucial role of the TMEM44-AS1-Myc axis in glioma progression and provides a possible anti-glioma therapeutic agent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02129-9.

A Comprehensive Analysis of SE-lncRNA/mRNA Differential Expression Profiles During Chondrogenic Differentiation of Human Bone Marrow Mesenchymal Stem Cells

Objective: Articular cartilage injury is common and difficult to treat clinically because of the characteristics of the cartilage. Bone marrow-derived mesenchymal stem cell (BMSC)-mediated cartilage regeneration is a promising therapy for treating articular cartilage injury. BMSC differentiation is controlled by numerous molecules and signaling pathways in the microenvironment at both the transcriptional and post-transcriptional levels. However, the possible function of super enhancer long non-coding RNAs (SE-lncRNAs) in the chondrogenic differentiation of BMSCs is still unclear. Our intention was to explore the expression profile of SE-lncRNAs and potential target genes regulated by SE-lncRNAs during chondrogenic differentiation in BMSCs.

Materials and Methods: In this study, we conducted a human Super-Enhancer LncRNA Microarray to investigate the differential expression profile of SE-lncRNAs and mRNAs during chondrogenic differentiation of BMSCs. Subsequent bioinformatic analysis was performed to clarify the important signaling pathways, SE-lncRNAs, and mRNAs associated with SE-lncRNAs regulating the chondrogenic differentiation of BMSCs.

Results: A total of 77 SE-lncRNAs were identified, of which 47 were upregulated and 30 were downregulated during chondrogenic differentiation. A total of 308 mRNAs were identified, of which 245 were upregulated and 63 were downregulated. Some pathways, such as focal adhesion, extracellular matrix (ECM)–receptor interaction, transforming growth factor-β (TGF-β) signaling pathway, and PI3K–Akt signaling pathway, were identified as the key pathways that may be implicated in the chondrogenic differentiation of BMSCs. Moreover, five potentially core regulatory mRNAs (PMEPA1, ENC1, TES, CDK6, and ADIRF) and 37 SE-lncRNAs in chondrogenic differentiation were identified by bioinformatic analysis.

Conclusion: We assessed the differential expression levels of SE-lncRNAs and mRNAs, along with the chondrogenic differentiation of BMSCs. By analyzing the interactions and co-expression, we identified the core SE-lncRNAs and mRNAs acting as regulators of the chondrogenic differentiation potential of BMSCs. Our study also provided novel insights into the mechanism of BMSC chondrogenic and cartilage regeneration.

A Comprehensive Toolbox to Analyze Enhancer-Promoter Functions

Knowledge in gene transcription and chromatin regulation has been intensely studied for decades, but thanks to next-generation sequencing (NGS) techniques there has been a major leap forward in the last few years. Historically, identification of specific enhancer elements has led to the identification of master transcription factors (TFs) in the 1990s. Genetic and biochemical experiments have identified the key regulators controlling RNA polymerase II (RNAPII) transcription and structurally analyses have elucidated detailed mechanisms. NGS and the development of chromatin immunoprecipitation (ChIP) have accelerated the gain of knowledge in the recent years. By now, we have a dazzling wealth of techniques that are currently used to put gene expression into a genome-wide context. This book is an attempt to assemble useful protocols for many researchers within and nearby research areas. In general, these innovative techniques focus on enhancer and promoter studies. The techniques should also be of interest for related fields such as DNA repair and replication.

Integrative analysis of genomic and epigenomic data reveals underlying super-enhancer-mediated microRNA regulatory network for human bone mineral density

Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined super-enhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range 'mSE-miRNA' approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).

Integrative machine learning framework for the identification of cell-specific enhancers from the human genome

Enhancers are deoxyribonucleic acid (DNA) fragments which when bound by transcription factors enhance the transcription of related genes. Due to its sporadic distribution and similar fractions, identification of enhancers from the human genome seems a daunting task. Compared to the traditional experimental approaches, computational methods with easy-to-use platforms could be efficiently applied to annotate enhancers' functions and physiological roles. In this aspect, several bioinformatics tools have been developed to identify enhancers. Despite their spectacular performances, existing methods have certain drawbacks and limitations, including fixed length of sequences being utilized for model development and cell-specificity negligence. A novel predictor would be beneficial in the context of genome-wide enhancer prediction by addressing the above-mentioned issues. In this study, we constructed new datasets for eight different cell types. Utilizing these data, we proposed an integrative machine learning (ML)-based framework called Enhancer-IF for identifying cell-specific enhancers. Enhancer-IF comprehensively explores a wide range of heterogeneous features with five commonly used ML methods (random forest, extremely randomized tree, multilayer perceptron, support vector machine and extreme gradient boosting). Specifically, these five classifiers were trained with seven encodings and obtained 35 baseline models. The output of these baseline models was integrated and again inputted to five classifiers for the construction of five meta-models. Finally, the integration of five meta-models through ensemble learning improved the model robustness. Our proposed approach showed an excellent prediction performance compared to the baseline models on both training and independent datasets in different cell types, thus highlighting the superiority of our approach in the identification of the enhancers. We assume that Enhancer-IF will be a valuable tool for screening and identifying potential enhancers from the human DNA sequences.