Integrative analysis and risk model construction for super‑enhancer‑related immune genes in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common type of kidney cancer associated with poor prognosis, and accounts for the majority of RCC-related deaths. The lack of comprehensive diagnostic and prognostic biomarkers has limited further understanding of the pathophysiology of ccRCC. Super-enhancers (SEs) are congregated enhancer clusters that have a key role in tumor processes such as epithelial-mesenchymal transition, metabolic reprogramming, immune escape and resistance to apoptosis. RCC may also be immunogenic and sensitive to immunotherapy. In the present study, an Arraystar human SE-long non-coding RNA (lncRNA) microarray was first employed to profile the differentially expressed SE-lncRNAs and mRNAs in 5 paired ccRCC and peritumoral tissues and to identify SE-related genes. The overlap of these genes with immune genes was then determined to identify SE-related immune genes. A model for predicting clinical prognosis and response to immunotherapy was built following the comprehensive analysis of a ccRCC gene expression dataset from The Cancer Genome Atlas (TCGA) database. The patients from TCGA were divided into high- and low-risk groups based on the median score derived from the risk model, and the Kaplan-Meier survival analysis showed that the low-risk group had a higher survival probability. In addition, according to the receiver operating characteristic curve analysis, the risk model had more advantages than other clinical factors in predicting the overall survival (OS) rate of patients with ccRCC. Using this model, it was demonstrated that the high-risk group had a more robust immune response. Furthermore, 61 potential drugs with half-maximal inhibitory concentration values that differed significantly between the two patient groups were screened to investigate potential drug treatment of ccRCC. In summary, the present study provided a novel index for predicting the survival probability of patients with ccRCC and may provide some insights into the mechanisms through which SE-related immune genes influence the diagnosis, prognosis and potential treatment drugs of ccRCC.

Super-enhancers: drivers of cells' identities and cells' debacles

Precise spatiotemporal regulations of gene expression are essential for determining cells' fates and functions. Enhancers are cis-acting DNA elements that act as periodic transcriptional thrusters and their activities are cell type specific. Clusters of enhancers, called super-enhancers, are more densely occupied by transcriptional activators than enhancers, driving stronger expression of their target genes, which have prominent roles in establishing and maintaining cellular identities. Here we review the current knowledge on the composition and structure of super-enhancers to understand how they robustly stimulate the expression of cellular identity genes. We also review their involvement in the development of various cell types and both noncancerous and cancerous disorders, implying the therapeutic interest of targeting them to fight against various diseases.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

Transcription factor condensates and signaling driven transcription

Transcription Factor (TF) condensates are a heterogenous mix of RNA, DNA, and multiple co-factor proteins capable of modulating the transcriptional response of the cell. The dynamic nature and the spatial location of TF-condensates in the 3D nuclear space is believed to provide a fast response, which is on the same pace as the signaling cascade and yet ever-so-specific in the crowded environment of the nucleus. However, the current understanding of how TF-condensates can achieve these feet so quickly and efficiently is still unclear. In this review, we draw parallels with other protein condensates and share our speculations on how the nucleus uses these TF-condensates to achieve high transcriptional specificity and fidelity. We discuss the various constituents of TF-condensates, their properties, and the known and unknown functions of TF-condensates with a particular focus on steroid signaling-induced transcriptional programs.

A body map of super-enhancers and their function in pig

Introduction

Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs.

Methods

Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression.

Results

We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression.

Discussion

In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Chromatin Landscape Governing Murine Epidermal Differentiation

Chromatin landscape and regulatory networks are determinants in lineage specification and differentiation. To define the temporospatial differentiation axis in murine epidermal cells in vivo, we generated datasets profiling expression dynamics (RNA sequencing), chromatin accessibility (assay for transposase-accessible chromatin using sequencing), architecture (Hi-C), and histone modifications (chromatin immunoprecipitation followed by sequencing) in the epidermis. We show that many differentially regulated genes are suppressed during the differentiation process, with superenhancers controlling differentiation-specific epigenomic changes. Our data shows the relevance of the Dlx/Klf/Grhl combinatorial regulatory network in maintaining correct temporospatial gene expression during epidermal differentiation. We determined differential open compartments, topologically associating domain score, and looping in the basal cell and suprabasal cell epidermal fractions, with the evolutionarily conserved epidermal differentiation complex region showing distinct suprabasal cell-specific topologically associating domain and loop formation that coincided with superenhancer sites. Overall, our study provides a global genome-wide resource of chromatin dynamics that define unrecognized regulatory networks and the epigenetic control of Dlx3-bound superenhancer elements during epidermal differentiation.

HOXDeRNA activates a cancerous transcription program and super-enhancers genome-wide

Background

The origin and genesis of highly malignant and heterogenous glioblastoma brain tumors remain unknown. We previously identified an enhancer-associated long non-coding RNA, LINC01116 (named HOXDeRNA here), that is absent in the normal brain but is commonly expressed in malignant glioma. HOXDeRNA has a unique capacity to transform human astrocytes into glioma-like cells. This work aimed to investigate molecular events underlying the genome-wide function of this lncRNA in glial cell fate and transformation.

Results

Using a combination of RNA-Seq, ChIRP-Seq, and ChIP-Seq, we now demonstrate that HOXDeRNA binds in trans to the promoters of genes encoding 44 glioma-specific transcription factors distributed throughout the genome and derepresses them by removing the Polycomb repressive complex 2 (PRC2). Among the activated transcription factors are the core neurodevelopmental regulators SOX2, OLIG2, POU3F2, and SALL2. This process requires an RNA quadruplex structure of HOXDeRNA that interacts with EZH2. Moreover, HOXDeRNA-induced astrocyte transformation is accompanied by the activation of multiple oncogenes such as EGFR, PDGFR, BRAF, and miR-21, and glioma-specific super-enhancers enriched for binding sites of glioma master transcription factors SOX2 and OLIG2.

Conclusions

Our results demonstrate that HOXDeRNA overrides PRC2 repression of glioma core regulatory circuitry with RNA quadruplex structure. These findings help reconstruct the sequence of events underlying the process of astrocyte transformation and suggest a driving role for HOXDeRNA and a unifying RNA-dependent mechanism of gliomagenesis.

Complexity of enhancer networks predicts cell identity and disease genes revealed by single-cell multi-omics analysis

Many enhancers exist as clusters in the genome and control cell identity and disease genes; however, the underlying mechanism remains largely unknown. Here, we introduce an algorithm, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression profiles. The complexity of enhancer networks is assessed by two metrics: the number of enhancers and the frequency of predicted enhancer interactions (PEIs) based on chromatin co-accessibility. We apply eNet algorithm to a human blood dataset and find cell identity and disease genes tend to be regulated by complex enhancer networks. The network hub enhancers (enhancers with frequent PEIs) are the most functionally important. Compared with super-enhancers, enhancer networks show better performance in predicting cell identity and disease genes. eNet is robust and widely applicable in various human or mouse tissues datasets. Thus, we propose a model of enhancer networks containing three modes: Simple, Multiple and Complex, which are distinguished by their complexity in regulating gene expression. Taken together, our work provides an unsupervised approach to simultaneously identify key cell identity and disease genes and explore the underlying regulatory relationships among enhancers in single cells.

The super-enhancer repertoire in porcine liver

The transcriptional initiation of genes is inextricably bound with the functions of cis-regulatory sequences. The pig is one of the most important livestock species and an ideal animal model for biomedical studies. At the same time, the liver is a critical organ with diverse and complex metabolic functions. Here, we performed Cleavage Under Targets and Tagmentation (CUT&Tag) coupled with high-throughput sequencing to profile the chromatin landscape of histone H3 lysine 27 acetylation (H3K27ac), histone H3 lysine 4 monomethylation (H3K4me1), and CCAAT enhancer-binding protein β (C-EBPβ) in the 70-d-old porcine liver, compared the different profiles among the three markers and their associated stitched-enhancers by stitching and sorting the peaks within 12.5 kb (Pott and Lieb, 2015) and generated the porcine liver-specific super-enhancers (SEs) by the combination of three markers. Compared to typical enhancers (TEs) and other stitched-enhancers, liver-specific SEs showed a higher density of cis-motifs and SNPs, which may recruit more tissue-specific vital TFs. The expression profiles in fetal and 70-d-old pigs proved that a large proportion of SE-associated genes were up-regulated and were more related to hepatic metabolisms and detoxification pathways. Our results illustrated the difference and connection among promoter and enhancer markers, identified the features of liver SEs and their associated genes, and provided novel insight into cis-element identification, function, and liver transcriptional regulation.

RANKL-responsive epigenetic mechanism reprograms macrophages into bone-resorbing osteoclasts

Monocyte/macrophage lineage cells are highly plastic and can differentiate into various cells under different environmental stimuli. Bone-resorbing osteoclasts are derived from the monocyte/macrophage lineage in response to receptor activator of NF-κB ligand (RANKL). However, the epigenetic signature contributing to the fate commitment of monocyte/macrophage lineage differentiation into human osteoclasts is largely unknown. In this study, we identified RANKL-responsive human osteoclast-specific superenhancers (SEs) and SE-associated enhancer RNAs (SE-eRNAs) by integrating data obtained from ChIP-seq, ATAC-seq, nuclear RNA-seq and PRO-seq analyses. RANKL induced the formation of 200 SEs, which are large clusters of enhancers, while suppressing 148 SEs in macrophages. RANKL-responsive SEs were strongly correlated with genes in the osteoclastogenic program and were selectively increased in human osteoclasts but marginally presented in osteoblasts, CD4+ T cells, and CD34+ cells. In addition to the major transcription factors identified in osteoclasts, we found that BATF binding motifs were highly enriched in RANKL-responsive SEs. The depletion of BATF1/3 inhibited RANKL-induced osteoclast differentiation. Furthermore, we found increased chromatin accessibility in SE regions, where RNA polymerase II was significantly recruited to induce the extragenic transcription of SE-eRNAs, in human osteoclasts. Knocking down SE-eRNAs in the vicinity of the NFATc1 gene diminished the expression of NFATc1, a major regulator of osteoclasts, and osteoclast differentiation. Inhibiting BET proteins suppressed the formation of some RANKL-responsive SEs and NFATc1-associated SEs, and the expression of SE-eRNA:NFATc1. Moreover, SE-eRNA:NFATc1 was highly expressed in the synovial macrophages of rheumatoid arthritis patients exhibiting high-osteoclastogenic potential. Our genome-wide analysis revealed RANKL-inducible SEs and SE-eRNAs as osteoclast-specific signatures, which may contribute to the development of osteoclast-specific therapeutic interventions.

Enhancer and super-enhancer landscape in polycystic kidney disease

Widespread aberrant gene expression is a pathological hallmark of polycystic kidney disease (PKD). Numerous pathogenic signaling cascades, including c-Myc, Fos, and Jun, are transactivated. However, the underlying epigenetic regulators are poorly defined. Here we show that H3K27ac, an acetylated modification of DNA packing protein histone H3 that marks active enhancers, is elevated in mouse and human samples of autosomal dominant PKD. Using comparative H3K27ac ChIP-Seq analysis, we mapped over 16000 active intronic and intergenic enhancer elements in Pkd1-mutant mouse kidneys. We found that the cystic kidney epigenetic landscape resembles that of a developing kidney, and over 90% of upregulated genes in Pkd1-mutant kidneys are co-housed with activated enhancers in the same topologically associated domains. Furthermore, we identified an evolutionarily conserved enhancer cluster downstream of the c-Myc gene and super-enhancers flanking both Jun and Fos loci in mouse and human models of autosomal dominant PKD. Deleting these regulatory elements reduced c-Myc, Jun, or Fos abundance and suppressed proliferation and 3D cyst growth of Pkd1-mutant cells. Finally, inhibiting glycolysis and glutaminolysis or activating Ppara in Pkd1-mutant cells lowerd global H3K27ac levels and its abundance on c-Myc enhancers. Thus, our work suggests that epigenetic rewiring mediates the transcriptomic dysregulation in PKD, and the regulatory elements can be targeted to slow cyst growth.

Expression of risk genes linked to vitamin D receptor super-enhancer regions and their association with phenotype severity in multiple sclerosis

Multiple sclerosis (MS) is a chronic debilitating neurological condition with a wide range of phenotype variability. A complex interplay of genetic and environmental factors contributes to disease onset and progression in MS patients. Vitamin D deficiency is a known susceptibility factor for MS, however the underlying mechanism of vitamin D-gene interactions in MS etiology is still poorly understood. Vitamin D receptor super-enhancers (VSEs) are enriched in MS risk variants and may modulate these environment-gene interactions. mRNA expression in total of 64 patients with contrasting MS severity was quantified in select genes. First, RNA-seq was performed on a discovery cohort (10 mild, 10 severe MS phenotype) and ten genes regulated by VSEs that have been linked to MS risk were analyzed. Four candidates showed a significant positive association (GRINA, PLEC, PARP10, and LRG1) in the discovery cohort and were then quantified using digital droplet PCR (ddPCR) in a validation cohort (33 mild, 11 severe MS phenotype). A significant differential expression persisted in the validation cohort for three of the VSE-MS genes: GRINA (p = 0.0138), LRG1 (p = 0.0157), and PLEC (p = 0.0391). In summary, genes regulated by VSE regions that contain known MS risk variants were shown to have differential expression based on disease severity (p<0.05). The findings implicate a role for vitamin D super-enhancers in modulating disease activity. In addition, expression levels may have some utility as prognostic biomarkers in the future.

MCIBox: a toolkit for single-molecule multi-way chromatin interaction visualization and micro-domains identification

The emerging ligation-free three-dimensional (3D) genome mapping technologies can identify multiplex chromatin interactions with single-molecule precision. These technologies not only offer new insight into high-dimensional chromatin organization and gene regulation, but also introduce new challenges in data visualization and analysis. To overcome these challenges, we developed MCIBox, a toolkit for multi-way chromatin interaction (MCI) analysis, including a visualization tool and a platform for identifying micro-domains with clustered single-molecule chromatin complexes. MCIBox is based on various clustering algorithms integrated with dimensionality reduction methods that can display multiplex chromatin interactions at single-molecule level, allowing users to explore chromatin extrusion patterns and super-enhancers regulation modes in transcription, and to identify single-molecule chromatin complexes that are clustered into micro-domains. Furthermore, MCIBox incorporates a two-dimensional kernel density estimation algorithm to identify micro-domains boundaries automatically. These micro-domains were stratified with distinctive signatures of transcription activity and contained different cell-cycle-associated genes. Taken together, MCIBox represents an invaluable tool for the study of multiple chromatin interactions and inaugurates a previously unappreciated view of 3D genome structure.

Effective Quality Breeding Directions-Comparison and Conservative Analysis of Hepatic Super-Enhancers between Chinese and Western Pig Breeds

The transcriptional initiation of genes is closely bound to the functions of cis-regulatory elements, including promoters, typical enhancers (TEs), and recently-identified super-enhancers (SEs). In this study, we identified these cis-regulatory elements in the livers of two Chinese (Meishan and Enshi Black) and two Western (Duroc and Large White) pig breeds using ChIP-seq data, then explored their similarities and differences. In addition, we analyzed the conservation of SEs among different tissues and species (pig, human, and mouse). We observed that SEs were more significantly enriched by transcriptional initiation regions, TF binding sites, and SNPs than other cis-elements. Western breeds included fewer SEs in number, while more growth-related QTLs were associated with these SEs. Additionally, the SEs were highly tissue-specific, and were conserved in the liver among humans, pigs, and mice. We concluded that intense selection could concentrate functional SEs; thus, SEs could be applied as effective detection regions in genomic selection breeding.

CDK12 orchestrates super-enhancer-associated CCDC137 transcription to direct hepatic metastasis in colorectal cancer

BACKGROUND

Hepatic metastasis is the primary and direct cause of death in individuals with colorectal cancer (CRC) attribute to lack of effective therapeutic targets. The present study aimed to identify potential druggable candidate targets for patients with liver metastatic CRC.

METHODS

The transcriptional profiles of super-enhancers (SEs) in primary and liver metastatic CRC were evaluated in publicly accessible CRC datasets. Immunohistochemistry of human CRC tissues was conducted to determine the expression level of CDK12. Cellular proliferation, survival and stemness were examined upon CDK12 inhibition by shCDK12 or a selective CDK12 inhibitor named SR-4835 with multiple in vitro and in vivo assays. RNA sequencing and bioinformatics analyses were carried out to investigate the mechanisms of CDK12 inhibition in CRC cells.

RESULTS

We identified CDK12 as a driver gene for direct hepatic metastasis in CRC. Suppression of CDK12 led to robust inhibition of proliferation, survival and stemness. Mechanistically, CDK12 intervention preferentially repressed the transcription of SE-associated genes. Integration of the SE landscape and RNA sequencing, BCL2L1 and CCDC137 were identified as SE-associated oncogenic genes to strengthen the abilities of cellular survival, proliferation and stemness, eventually increasing liver metastasis of CRC.

CONCLUSIONS

Our data highlight the potential of CDK12 and SE-associated oncogenic transcripts as therapeutic targets for patients with liver metastatic CRC.

Acetylated HOXB13 Regulated Super Enhancer Genes Define Therapeutic Vulnerabilities of Castration-Resistant Prostate Cancer

PURPOSE

Androgen receptor (AR) antagonism is exacerbated by HOXB13 in castration-resistant prostate cancers (CRPC). However, it is unclear when and how HOXB13 primes CRPCs for AR antagonism. By mass-spectrometry analysis of CRPC extract, we uncovered a novel lysine 13 (K13) acetylation in HOXB13 mediated by CBP/p300. To determine whether acetylated K13-HOXB13 is a clinical biomarker of CRPC development, we characterized its role in prostate cancer biology.

EXPERIMENTAL DESIGN

We identified tumor-specific acK13-HOXB13 signal enriched super enhancer (SE)-regulated targets. We analyzed the effect of loss of HOXB13K13-acetylation on chromatin binding, SE proximal target gene expression, self-renewal, enzalutamide sensitivity, and CRPC tumor growth by employing isogenic parental and HOXB13K13A mutants. Finally, using primary human prostate organoids, we evaluated whether inhibiting an acK13-HOXB13 target, ACK1, with a selective inhibitor (R)-9b is superior to AR antagonists in inhibiting CRPC growth.

RESULTS

acK13-HOXB13 promotes increased expression of lineage (AR, HOXB13), prostate cancer diagnostic (FOLH1), CRPC-promoting (ACK1), and angiogenesis (VEGFA, Angiopoietins) genes early in prostate cancer development by establishing tumor-specific SEs. acK13-HOXB13 recruitment to key SE-regulated targets is insensitive to enzalutamide. ACK1 expression is significantly reduced in the loss of function HOXB13K13A mutant CRPCs. Consequently, HOXB13K13A mutants display reduced self-renewal, increased sensitivity to enzalutamide, and impaired xenograft tumor growth. Primary human prostate tumor organoids expressing HOXB13 are significantly resistant to AR antagonists but sensitive to (R)-9b.

CONCLUSIONS

In summary, acetylated HOXB13 is a biomarker of clinically significant prostate cancer. Importantly, PSMA-targeting agents and (R)-9b could be new therapeutic modalities to target HOXB13-ACK1 axis regulated prostate cancers.

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.

Functional characterization and clinical significance of super-enhancers in lung adenocarcinoma

Super-enhancers (SEs) are important transcriptional regulators in tumorigenesis; however, the functional characterization and clinical significance of SEs in lung adenocarcinoma (LUAD) remain unclear. By using H3K27ac ChIP-seq data of two LUAD cell lines and eight lung tissues, we detected 1045 cancer-specific and 5032 normal-specific SEs. Compared to normal-specific SEs, cancer-specific SEs have different regulatory mechanisms where associated target genes were enriched in critical tumor-related pathways and tended to be regulated by transcription factors of Fos Proto-Oncogene, AP-1 Transcription Factor Subunit and Jun Proto-Oncogene, AP-1 Transcription Factor Subunit families. By using expression data of 513 LUAD and 57 adjacent samples from The Cancer Genome Atlas and 80 tumor-normal paired LUAD samples from the Nanjing Lung Cancer Cohort study, we performed differential expression analysis of target genes for SEs and defined 243 crucial SEs. Unsupervised clustering of crucial SEs revealed two subtypes with different levels of genomic aberrations (i.e., mutation and copy number alteration) and clinical outcomes (progression-free interval: p = 0.030; disease-free interval: p = 0.047). In addition, patients with adverse clinical outcomes were more sensitive to three small molecule inhibitors (bortezomib, doxorubicin, and etoposide), and their targets (PSMB5 and TOP2A) also have elevated expression levels among these patients. Taken together, our findings provided a comprehensive characterization of SEs in LUAD and emphasized their clinical significance in LUAD therapy.

Super-Enhancers Dysregulations in Hematological Malignancies

Hematological malignancies affecting either the lymphoid or the myeloid lineages involve epigenetic mutations or dysregulation in the majority of cases. These epigenetic abnormalities can affect regulatory elements in the genome and, particularly, enhancers. Recently, large regulatory elements known as super-enhancers, initially identified for their critical roles in cell-type specific expression regulation of genes controlling cell identity, have been shown to also be involved in tumorigenesis in many cancer types and hematological malignancies via the regulation of numerous oncogenes, including MYC. In this review, we highlight the existing links between super-enhancers and hematological malignancies, with a particular focus on acute myeloid leukemia, a clonal hematopoietic neoplasm with dismal outcomes, resulting in an uncontrolled proliferation of myeloblasts, abnormally blocked during differentiation and accumulating within the patient's bone marrow. We report recent works, performed during the last few years, treating this subject and consider the possibility of targeting oncogenic regulatory elements, as well as the effectiveness and limitations reported so far for such strategies.

Super-enhancer associated core regulatory circuits mediate susceptibility to retinoic acid in neuroblastoma cells

Neuroblastoma is a pediatric tumour that accounts for more than 15% of cancer-related deaths in children. High-risk tumours are often difficult to treat, and patients' survival chances are less than 50%. Retinoic acid treatment is part of the maintenance therapy given to neuroblastoma patients; however, not all tumours differentiate in response to retinoic acid. Within neuroblastoma tumors, two phenotypically distinct cell types have been identified based on their super-enhancer landscape and transcriptional core regulatory circuitries: adrenergic (ADRN) and mesenchymal (MES). We hypothesized that the distinct super-enhancers in these different tumour cells mediate differential response to retinoic acid. To this end, three different neuroblastoma cell lines, ADRN (MYCN amplified and non-amplified) and MES cells, were treated with retinoic acid, and changes in the super-enhancer landscape upon treatment and after subsequent removal of retinoic acid was studied. Using ChIP-seq for the active histone mark H3K27ac, paired with RNA-seq, we compared the super-enhancer landscape in cells that undergo neuronal differentiation in response to retinoic acid versus those that fail to differentiate and identified unique super-enhancers associated with neuronal differentiation. Among the ADRN cells that respond to treatment, MYCN-amplified cells remain differentiated upon removal of retinoic acid, whereas MYCN non-amplified cells revert to an undifferentiated state, allowing for the identification of super-enhancers responsible for maintaining differentiation. This study identifies key super-enhancers that are crucial for retinoic acid-mediated differentiation.

Super-Enhancer-Associated Long Non-Coding RNA LINC01485 Promotes Osteogenic Differentiation of Human Bone Marrow Mesenchymal Stem Cells by Regulating MiR-619-5p/RUNX2 Axis

OBJECTIVE

To investigate the mechanisms of super-enhancer-associated LINC01485/miR-619-5p/RUNX2 signaling axis involvement in osteogenic differentiation of human bone marrow mesenchymal stem cells (hBMSCs).

METHODS

Osteogenic differentiation of hBMSCs was induced in vitro. The expression levels of LINC01485 and miR-619-5p during osteogenesis were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Osteogenic differentiation was examined by qRT-PCR, western blot, alkaline phosphatase (ALP) staining, ALP activity measurement, and Alizarin Red S (ARS) staining assays. Thereafter, the effects of LINC01485 and miR-619-5p on osteogenic differentiation of hBMSCs were evaluated by performing loss- and gain-of-function experiments. Subsequently, a fluorescence in situ hybridization (FISH) assay was employed to determine the cellular localization of LINC01485. Bioinformatics analysis, RNA antisense purification (RAP) assay, and dual-luciferase reporter assays were conducted to analyze the interactions of LINC01485, miR-619-5p, and RUNX2. Rescue experiments were performed to further delineate the role of the competitive endogenous RNA (ceRNA) signaling axis consisting of LINC01485/miR-619-5p/RUNX2 in osteogenic differentiation of hBMSCs.

RESULTS

The expression of LINC01485 was up-regulated during osteogenic differentiation of hBMSCs. The overexpression of LINC01485 promoted osteogenic differentiation of hBMSCs by up-regulating the expression of osteogenesis-related genes [e.g., runt-related transcription factor 2 (RUNX2), osterix (OSX), collagen type 1 alpha 1 (COL1A1), osteocalcin (OCN), and osteopontin (OPN)], and increasing the activity of ALP. ALP staining and ARS staining were also found to be increased upon overexpression of LINC01485. The opposing results were obtained upon LINC01485 interference in hBMSCs. miR-619-5p was found to inhibit osteogenic differentiation. FISH assay displayed that LINC01485 was mainly localized in the cytoplasm. RAP assay results showed that LINC01485 bound to miR-619-5p, and dual-luciferase reporter assay verified that LINC01485 bound to miR-619-5p, while miR-619-5p and RUNX2 bound to each other. Rescue experiments illustrated that LINC01485 could promote osteogenesis by increasing RUNX2 expression by sponging miR-619-5p.

CONCLUSION

LINC01485 could influence RUNX2 expression by acting as a ceRNA of miR-619-5p, thereby promoting osteogenic differentiation of hBMSCs. The LINC01485/miR-619-5p/RUNX2 axis might comprise a novel target in the bone tissue engineering field.

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

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

Primordial super-enhancers: heat shock-induced chromatin organization in yeast

Specialized mechanisms ensure proper expression of critically important genes such as those specifying cell identity or conferring protection from environmental stress. Investigations of the heat shock response have been critical in elucidating basic concepts of transcriptional control. Recent studies demonstrate that in response to thermal stress, heat shock-responsive genes associate with high levels of transcriptional activators and coactivators and those in yeast intensely interact across and between chromosomes, coalescing into condensates. In mammalian cells, cell identity genes that are regulated by super-enhancers (SEs) are also densely occupied by transcriptional machinery that form phase-separated condensates. We suggest that the stress-remodeled yeast nucleome bears functional and structural resemblance to mammalian SEs, and will reveal fundamental mechanisms of gene control by transcriptional condensates.

The Role of rDNA Clusters in Global Epigenetic Gene Regulation

The regulation of gene expression has been studied for decades, but the underlying mechanisms are still not fully understood. As well as local and distant regulation, there are specific mechanisms of regulation during development and physiological modulation of gene activity in differentiated cells. Current research strongly supports a role for the 3D chromosomal structure in the regulation of gene expression. However, it is not known whether the genome structure reflects the formation of active or repressed chromosomal domains or if these structures play a primary role in the regulation of gene expression. During early development, heterochromatinization of ribosomal DNA (rDNA) is coupled with silencing or activation of the expression of different sets of genes. Although the mechanisms behind this type of regulation are not known, rDNA clusters shape frequent inter-chromosomal contacts with a large group of genes controlling development. This review aims to shed light on the involvement of clusters of ribosomal genes in the global regulation of gene expression. We also discuss the possible role of RNA-mediated and phase-separation mechanisms in the global regulation of gene expression by nucleoli.

Epigenomic features revealed by ATAC-seq impact transgene expression in CHO cells

Different regions of a mammalian genome have different accessibilities to transcriptional machinery. The integration site of a transgene affects how actively it is transcribed. Highly accessible genomic regions called super-enhancers have been recently described as strong regulatory elements that shape cell identity. Super-enhancers have been identified in Chinese hamster ovary (CHO) cells using the Assay for Transposase-Accessible Chromatin Sequencing (ATAC-seq). Genes near super-enhancer regions had high transcript levels and were enriched for oncogenic signaling and proliferation functions, consistent with an immortalized phenotype. Inaccessible regions in the genome with low ATAC signal also had low transcriptional activity. Genes in inaccessible regions were enriched for remote tissue functions such as taste, smell, and neuronal activation. A lentiviral reporter integration assay showed integration into super-enhancer regions conferred higher reporter expression than insertion into inaccessible regions. Targeted integration of an IgG vector into the Plec super-enhancer region yielded clones that expressed the immunoglobulin light chain gene mostly in the top 20% of all transcripts with the majority in the top 5%. The results suggest the epigenomic landscape of CHO cells can guide the selection of integration sites in the development of cell lines for therapeutic protein production.

Epigenetics and expression of key genes associated with cardiac fibrosis: NLRP3, MMP2, MMP9, CCN2/CTGF and AGT

Aims:

Excessive inflammatory signaling and pathological remodeling of the extracellular matrix drive cardiac fibrosis and require changes in gene expression.

Materials and methods:

Using bioinformatics, both tissue-specific expression profiles and epigenomic profiles of some genes critical for cardiac fibrosis were examined, namely, NLRP3, MMP2, MMP9, CCN2/CTGF, AGT (encodes angiotensin II precursors) and hsa-mir-223 (post-transcriptionally regulates NLRP3).

Results:

In monocytes, neutrophils, fibroblasts, venous cells, liver and brain, enhancers or super-enhancers were found that correlate with high expression of these genes. One enhancer extended into a silent gene neighbor. These enhancers harbored tissue-specific foci of DNA hypomethylation, open chromatin and transcription factor binding.

Conclusions:

This study identified previously undescribed enhancers containing hypomethylated transcription factor binding subregions that are predicted to regulate expression of these cardiac fibrosis-inducing genes.

Epigenome screening highlights that JMJD6 confers an epigenetic vulnerability and mediates sunitinib sensitivity in renal cell carcinoma

Aberrant epigenetic reprogramming represents a hallmark of renal cell carcinoma (RCC) tumorigenesis and progression. Whether there existed other epigenetic vulnerabilities that could serve as therapeutic targets remained unclear and promising. Here, we combined the clustered regularly interspaced short palindromic repeats functional screening results and multiple RCC datasets to identify JMJD6 as the potent target in RCC. JMJD6 expression correlated with poor survival outcomes of RCC patients and promoted RCC progression in vitro and in vivo. Mechanistically, aberrant p300 led to high JMJD6 expression, which activated a series of oncogenic crosstalk. Particularly, high-throughput sequencing data revealed that JMJD6 could assemble super-enhancers to drive a list of identity genes in kidney cancer, including VEGFA, β-catenin, and SRC. Moreover, this JMJD6-mediated oncogenic effect could be suppressed by a novel JMJD6 inhibitor (SKLB325), which was further demonstrated in RCC cells, patient-derived organoid models, and in vivo. Given the probable overlapped crosstalk between JMJD6 signature and tyrosine kinase inhibitors downstream targets, targeting JMJD6 sensitized RCC to sunitinib and was synergistic when they were combined together. Collectively, this study indicated that targeting JMJD6 was an effective approach to treat RCC patients.

Downregulated METTL14 accumulates BPTF that reinforces super-enhancers and distal lung metastasis via glycolytic reprogramming in renal cell carcinoma

Background: Methyltransferase-like 14 (METTL14) participates in tumorigenesis in several malignancies, but how METTL14 mediates the metastasis of renal cell carcinoma (RCC) has never been reported. Methods: Western blotting, quantitative real-time PCR, and immunohistochemistry were used to determine the mRNA and protein levels of relevant genes. Methylated RNA immunoprecipitation sequencing and RNA sequencing were utilized to screen potential targets of METTL14. Chromatin immunoprecipitation sequencing and assay for transposase-accessible chromatin sequencing were performed to investigate epigenetic alterations. The biological roles and mechanisms of METTL14/BPTF in promoting lung metastasis were confirmed in vitro and in vivo using cell lines, patient samples, xenograft models, and organoids. Results: Utilizing the TCGA-KIRC and Ruijin-RCC datasets, we found low expression of METTL14 in mRCC samples, which predicted poor prognosis. METTL14 deficiency promoted RCC metastasis in vitro and in vivo. Mechanistically, METTL14-mediated m⁶A modification negatively regulated the mRNA stability of bromodomain PHD finger transcription factor (BPTF) and depended on BPTF to drive lung metastasis. Accumulated BPTF in METTL14-deficient cells remodeled the enhancer landscape to reinforce several oncogenic crosstalk. Particularly, BPTF constituted super-enhancers that activate downstream targets like enolase 2 and SRC proto-oncogene nonreceptor tyrosine kinase, leading to glycolytic reprogramming of METTL14^-/- cells. Finally, we determined the efficacy of the BPTF inhibitor AU1 in suppressing mRCC of patient-derived cells, mRCC-derived organoids (MDOs), and orthotopic xenograft models. Conclusions: Our study is the first to investigate the essential role of m⁶A modification and the METTL14/BPTF axis in the epigenetic and metabolic remodeling of mRCC, highlighting AU1 as a vital therapeutic candidate.

Super-Enhancers in the Regulation of Gene Transcription: General Aspects and Antitumor Targets

Super-enhancers (genome elements that activate gene transcription) are DNA regions with an elevated concentration of transcriptional complexes. These multiprotein structures contain, among other components, the cyclin-dependent kinases 8 and 19. These and other transcriptional protein kinases are regarded as novel targets for pharmacological inhibition by antitumor drug candidates.

Genome editing demonstrates that the -5 kb Nanog enhancer regulates Nanog expression by modulating RNAPII initiation and/or recruitment

Transcriptional enhancers have been defined by their ability to operate independent of distance and orientation in plasmid-based reporter assays of gene expression. At present, histone marks are used to identify and define enhancers but do not consider the endogenous role of an enhancer in the context of native chromatin. We employed a combination of genomic editing, single cell analyses, and sequencing approaches to investigate a Nanog-associated cis-regulatory element, which has been reported by others to be either an alternative promoter or a super-enhancer. We first demonstrate both distance and orientation independence in native chromatin, eliminating the issues raised with plasmid-based approaches. We next demonstrate that the dominant super-enhancer modulates Nanog globally and operates by recruiting and/or initiating RNA Polymerase II. Our studies have important implications to how transcriptional enhancers are defined and how they regulate gene expression.

Mediator's Kinase Module: A Modular Regulator of Cell Fate

Selective gene expression is crucial in maintaining the self-renewing and multipotent properties of stem cells. Mediator is a large, evolutionarily conserved, multi-subunit protein complex that modulates gene expression by relaying signals from cell type-specific transcription factors to RNA polymerase II. In humans, this complex consists of 30 subunits arranged in four modules. One critical module of the Mediator complex is the kinase module consisting of four subunits: MED12, MED13, CDK8, and CCNC. The kinase module exists in variable association with the 26-subunit Mediator core and affects transcription through phosphorylation of transcription factors and by controlling Mediator structure and function. Many studies have shown the kinase module to be a key player in the maintenance of stem cells that is distinct from a general role in transcription. Genetic studies have revealed that dysregulation of this kinase subunit contributes to the development of many human diseases. In this review, we discuss the importance of the Mediator kinase module by examining how this module functions with the more recently identified transcriptional super-enhancers, how changes in the kinase module and its activity can lead to the development of human disease, and the role of this unique module in directing and maintaining cell state. As we look to use stem cells to understand human development and treat human disease through both cell-based therapies and tissue engineering, we need to remain aware of the on-going research and address critical gaps in knowledge related to the molecular mechanisms that control cell fate.

Epigenetics of Muscle- and Brain-Specific Expression of KLHL Family Genes

KLHL and the related KBTBD genes encode components of the Cullin-E3 ubiquitin ligase complex and typically target tissue-specific proteins for degradation, thereby affecting differentiation, homeostasis, metabolism, cell signaling, and the oxidative stress response. Despite their importance in cell function and disease (especially, KLHL40, KLHL41, KBTBD13, KEAP1, and ENC1), previous studies of epigenetic factors that affect transcription were predominantly limited to promoter DNA methylation. Using diverse tissue and cell culture whole-genome profiles, we examined 17 KLHL or KBTBD genes preferentially expressed in skeletal muscle or brain to identify tissue-specific enhancer and promoter chromatin, open chromatin (DNaseI hypersensitivity), and DNA hypomethylation. Sixteen of the 17 genes displayed muscle- or brain-specific enhancer chromatin in their gene bodies, and most exhibited specific intergenic enhancer chromatin as well. Seven genes were embedded in super-enhancers (particularly strong, tissue-specific clusters of enhancers). The enhancer chromatin regions typically displayed foci of DNA hypomethylation at peaks of open chromatin. In addition, we found evidence for an intragenic enhancer in one gene upregulating expression of its neighboring gene, specifically for KLHL40/HHATL and KLHL38/FBXO32 gene pairs. Many KLHL/KBTBD genes had tissue-specific promoter chromatin at their 5' ends, but surprisingly, two (KBTBD11 and KLHL31) had constitutively unmethylated promoter chromatin in their 3' exons that overlaps a retrotransposed KLHL gene. Our findings demonstrate the importance of expanding epigenetic analyses beyond the 5' ends of genes in studies of normal and abnormal gene regulation.

Rapid Enhancer Remodeling and Transcription Factor Repurposing Enable High Magnitude Gene Induction upon Acute Activation of NK Cells

Innate immune responses rely on rapid and precise gene regulation mediated by accessibility of regulatory regions to transcription factors (TFs). In natural killer (NK) cells and other innate lymphoid cells, competent enhancers are primed during lineage acquisition, and formation of de novo enhancers characterizes the acquisition of innate memory in activated NK cells and macrophages. Here, we investigated how primed and de novo enhancers coordinate to facilitate high-magnitude gene induction during acute activation. Epigenomic and transcriptomic analyses of regions near highly induced genes (HIGs) in NK cells both in vitro and in a model of Toxoplasma gondii infection revealed de novo chromatin accessibility and enhancer remodeling controlled by signal-regulated TFs STATs. Acute NK cell activation redeployed the lineage-determining TF T-bet to de novo enhancers, independent of DNA-sequence-specific motif recognition. Thus, acute stimulation reshapes enhancer function through the combinatorial usage and repurposing of both lineage-determining and signal-regulated TFs to ensure an effective response.

Truncated BRPF1 Cooperates with Smoothened to Promote Adult Shh Medulloblastoma

The transition of neural progenitors to differentiated postmitotic neurons is mainly considered irreversible in physiological conditions. In the present work, we show that Shh pathway activation through SmoM2 expression promotes postmitotic neurons dedifferentiation, re-entering in the cell cycle and originating medulloblastoma in vivo. Notably, human adult patients present inactivating mutations of the chromatin reader BRPF1 that are associated with SMO mutations and absent in pediatric and adolescent patients. Here, we found that truncated BRPF1 protein, as found in human adult patients, is able to induce medulloblastoma in adult mice upon SmoM2 activation. Indeed, postmitotic neurons re-entered the cell cycle and proliferated as a result of chromatin remodeling of neurons by BRPF1. Our model of brain cancer explains the onset of a subset of human medulloblastoma in adult individuals where granule neuron progenitors are no longer present.

A Tumor-Specific Super-Enhancer Drives Immune Evasion by Guiding Synchronous Expression of PD-L1 and PD-L2

PD-L1 and PD-L2 are important targets for immune checkpoint blockade, but how tumor cells achieve their expression remains to be addressed. Here, we find that PD-L1 and PD-L2 are co-expressed in cancer cell lines and tissues across different cancer types. In breast cancer, MDA-MB-231 and SUM-159 cells show high expression of both PD-L1 and PD-L2. The expression of both PD-L1 and PD-L2 is greatly reduced upon treatment of inhibitors of super-enhancers. Bioinformatic analysis identifies a potential super-enhancer (PD-L1L2-SE) that is located between the CD274 and CD273 genes. Genetic deletion of PD-L1L2-SE profoundly reduces the expression of PD-L1 and PD-L2. PD-L1L2-SE-deficient cancer cells fail to generate immune evasion and are sensitive to T cell-mediated killing. Notably, epigenetic activation of such a region (PD-L1L2-SE) is correlated with PD-L1 and PD-L2. Taken together, we identify a super-enhancer (PD-L1L2-SE) that is responsible for the overexpression of PD-L1 and PD-L2 as well as immune evasion in cancer.

Chromatin profiling reveals relocalization of lysine-specific demethylase 1 by an oncogenic fusion protein

Paediatric cancers commonly harbour quiet mutational landscapes and are instead characterized by single driver events such as the mutation of critical chromatin regulators, expression of oncohistones, or expression of oncogenic fusion proteins. These events ultimately promote malignancy through disruption of normal gene regulation and development. The driver protein in Ewing sarcoma, EWS/FLI, is an oncogenic fusion and transcription factor that reshapes the enhancer landscape, resulting in widespread transcriptional dysregulation. Lysine-specific demethylase 1 (LSD1) is a critical functional partner for EWS/FLI as inhibition of LSD1 reverses the transcriptional activity of EWS/FLI. However, how LSD1 participates in fusion-directed epigenomic regulation and aberrant gene activation is unknown. We now show EWS/FLI causes dynamic rearrangement of LSD1 and we uncover a role for LSD1 in gene activation through colocalization at EWS/FLI binding sites throughout the genome. LSD1 is integral to the establishment of Ewing sarcoma super-enhancers at GGAA-microsatellites, which ubiquitously overlap non-microsatellite loci bound by EWS/FLI. Together, we show that EWS/FLI induces widespread changes to LSD1 distribution in a process that impacts the enhancer landscape throughout the genome.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

How do maternal transcription factors interact with chromatin regions to coordinate the endodermal gene regulatory program? Paraiso et al. demonstrate that combinatorial binding of maternal Otx1, Vegt, and Foxh1 to select enhancers and super-enhancers in the genome controls endodermal cell fate specification during zygotic gene activation.

Functions and Clinical Significance of Super-Enhancers in Bone-Related Diseases

Super-enhancers (SEs) are a large cluster of cis-regulatory DNA elements that contain many binding motifs, which master transcription factors and cofactors bind to with high density. SEs usually regulate the expression of genes that can control the cell identity and fate, and SEs can be used to explain the patterns of the expression of cell-specific genes. Hence, it shows great potential for application in the treatment of diseases like cancer. At present, the clinical treatments for osteosarcoma, Ewing sarcoma, and other bone-related diseases remain challenging. The poor prognosis and difficult treatment of these diseases imposes heavy economic burden on patients and society. In recent years, research on SEs with respect to bone-related diseases has attracted increasing attention. In this paper, we first review the identification and functional mechanisms of SEs. Then, we integrate the findings of the emerging studies on SEs in bone-related diseases. Finally, we summarize recent strategies for targeting SEs for the treatment of bone-related diseases. This review aims to provide comprehensive insights into the roles of SEs in bone-related diseases.

Core regulatory circuitries in defining cancer cell identity across the malignant spectrum

Gene expression programmes driving cell identity are established by tightly regulated transcription factors that auto- and cross-regulate in a feed-forward manner, forming core regulatory circuitries (CRCs). CRC transcription factors create and engage super-enhancers by recruiting acetylation writers depositing permissive H3K27ac chromatin marks. These super-enhancers are largely associated with BET proteins, including BRD4, that influence higher-order chromatin structure. The orchestration of these events triggers accessibility of RNA polymerase machinery and the imposition of lineage-specific gene expression. In cancers, CRCs drive cell identity by superimposing developmental programmes on a background of genetic alterations. Further, the establishment and maintenance of oncogenic states are reliant on CRCs that drive factors involved in tumour development. Hence, the molecular dissection of CRC components driving cell identity and cancer state can contribute to elucidating mechanisms of diversion from pre-determined developmental programmes and highlight cancer dependencies. These insights can provide valuable opportunities for identifying and re-purposing drug targets. In this article, we review the current understanding of CRCs across solid and liquid malignancies and avenues of investigation for drug development efforts. We also review techniques used to understand CRCs and elaborate the indication of discussed CRC transcription factors in the wider context of cancer CRC models.

Perspectives on the use of super-enhancers as a defining feature of cell/tissue-identity genes

Super-enhancers (SE) have become a popular concept and are widely used as a feature defining key identity genes. Here, we provide perspectives on the use of SE to define and identify cell/tissue-identity genes. By mining SE and their associated genes using murine functional genomics data, we highlight and discuss current limitations and open questions regarding both the sensitivity and specificity of identity genes/transcription factors predicted by SE. In this context, we point to cell/tissue-specific promoters as an important additional level of information, which we propose to combine with SE when aiming to define potential identity genes.

Chromatin structure regulates cancer-specific alternative splicing events in primary HPV-related oropharyngeal squamous cell carcinoma

Human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV+ OPSCC) represents a unique disease entity within head and neck cancer with rising incidence. Previous work has shown that alternative splicing events (ASEs) are prevalent in HPV+ OPSCC, but further validation is needed to understand the regulation of this process and its role in these tumours. In this study, eleven ASEs (GIT2, CTNNB1, MKNK2, MRPL33, SIPA1L3, SNHG6, SYCP2, TPRG1, ZHX2, ZNF331, and ELOVL1) were selected for validation from 109 previously published candidate ASEs to elucidate the post-transcriptional mechanisms of oncogenesis in HPV+ disease. In vitro qRT-PCR confirmed differential expression of 9 of 11 ASE candidates, and in silico analysis within the TCGA cohort confirmed 8 of 11 candidates. Six ASEs (MRPL33, SIPA1L3, SNHG6, TPRG1, ZHX2, and ELOVL1) showed significant differential expression across both methods. Further evaluation of chromatin modification revealed that ASEs strongly correlated with cancer-specific distribution of acetylated lysine 27 of histone 3 (H3K27ac). Subsequent epigenetic treatment of HPV+ HNSCC cell lines (UM-SCC-047 and UPCI-SCC-090) with JQ1 not only induced downregulation of cancer-specific ASE isoforms, but also growth inhibition in both cell lines. The UPCI-SCC-090 cell line, with greater ASE expression, also showed more significant growth inhibition after JQ1 treatment. This study confirms several novel cancer-specific ASEs in HPV+OPSCC and provides evidence for the role of chromatin modifications in regulation of alternative splicing in HPV+OPSCC. This highlights the role of epigenetic changes in the oncogenesis of HPV+OPSCC, which represents a unique, unexplored target for therapeutics that can alter the global post-transcriptional landscape.

Epigenetics of Skeletal Muscle-Associated Genes in the ASB, LRRC, TMEM, and OSBPL Gene Families

Much remains to be discovered about the intersection of tissue-specific transcription control and the epigenetics of skeletal muscle (SkM), a very complex and dynamic organ. From four gene families, Leucine-Rich Repeat Containing (LRRC), Oxysterol Binding Protein Like (OSBPL), Ankyrin Repeat and Socs Box (ASB), and Transmembrane Protein (TMEM), we chose 21 genes that are preferentially expressed in human SkM relative to 52 other tissue types and analyzed relationships between their tissue-specific epigenetics and expression. We also compared their genetics, proteomics, and descriptions in the literature. For this study, we identified genes with little or no previous descriptions of SkM functionality (ASB4, ASB8, ASB10, ASB12, ASB16, LRRC14B, LRRC20, LRRC30, TMEM52, TMEM233, OSBPL6/ORP6, and OSBPL11/ORP11) and included genes whose SkM functions had been previously addressed (ASB2, ASB5, ASB11, ASB15, LRRC2, LRRC38, LRRC39, TMEM38A/TRIC-A, and TMEM38B/TRIC-B). Some of these genes have associations with SkM or heart disease, cancer, bone disease, or other diseases. Among the transcription-related SkM epigenetic features that we identified were: super-enhancers, promoter DNA hypomethylation, lengthening of constitutive low-methylated promoter regions, and SkM-related enhancers for one gene embedded in a neighboring gene (e.g., ASB8-PFKM, LRRC39-DBT, and LRRC14B-PLEKHG4B gene-pairs). In addition, highly or lowly co-expressed long non-coding RNA (lncRNA) genes probably regulate several of these genes. Our findings give insights into tissue-specific epigenetic patterns and functionality of related genes in a gene family and can elucidate normal and disease-related regulation of gene expression in SkM.

Non-overlapping Control of Transcriptome by Promoter- and Super-Enhancer-Associated Dependencies in Multiple Myeloma

The relationship between promoter proximal transcription factor-associated gene expression and super-enhancer-driven transcriptional programs are not well defined. However, their distinct genomic occupancy suggests a mechanism for specific and separable gene control. We explored the transcriptional and functional interrelationship between E2F transcription factors and BET transcriptional co-activators in multiple myeloma. We found that the transcription factor E2F1 and its heterodimerization partner DP1 represent a dependency in multiple myeloma cells. Global chromatin analysis reveals distinct regulatory axes for E2F and BETs, with E2F predominantly localized to active gene promoters of growth and/or proliferation genes and BETs disproportionately at enhancer-regulated tissue-specific genes. These two separate gene regulatory axes can be simultaneously targeted to impair the myeloma proliferative program, providing an important molecular mechanism for combination therapy. This study therefore suggests a sequestered cellular functional control that may be perturbed in cancer with potential for development of a promising therapeutic strategy.

Uncontrolled proliferation is a hallmark of tumorigenesis and is associated with perturbed transcriptomic profile. Fulciniti et al. explored the interrelationship between E2F transcription factors and BET transcriptional co-activators in multiple myeloma, reporting the existence of two distinct regulatory axes that can be synergistically targeted to impact myeloma growth and survival.

Reconstruction of the Global Neural Crest Gene Regulatory Network In Vivo

Precise control of developmental processes is encoded in the genome in the form of gene regulatory networks (GRNs). Such multi-factorial systems are difficult to decode in vertebrates owing to their complex gene hierarchies and dynamic molecular interactions. Here we present a genome-wide in vivo reconstruction of the GRN underlying development of the multipotent neural crest (NC) embryonic cell population. By coupling NC-specific epigenomic and transcriptional profiling at population and single-cell levels with genome/epigenome engineering in vivo, we identify multiple regulatory layers governing NC ontogeny, including NC-specific enhancers and super-enhancers, novel trans-factors, and cis-signatures allowing reverse engineering of the NC-GRN at unprecedented resolution. Furthermore, identification and dissection of divergent upstream combinatorial regulatory codes has afforded new insights into opposing gene circuits that define canonical and neural NC fates early during NC ontogeny. Our integrated approach, allowing dissection of cell-type-specific regulatory circuits in vivo, has broad implications for GRN discovery and investigation.

Transcriptional Dysregulation of MYC Reveals Common Enhancer-Docking Mechanism

Transcriptional dysregulation of the MYC oncogene is among the most frequent events in aggressive tumor cells, and this is generally accomplished by acquisition of a super-enhancer somewhere within the 2.8 Mb TAD where MYC resides. We find that these diverse cancer-specific super-enhancers, differing in size and location, interact with the MYC gene through a common and conserved CTCF binding site located 2 kb upstream of the MYC promoter. Genetic perturbation of this enhancer-docking site in tumor cells reduces CTCF binding, super-enhancer interaction, MYC gene expression, and cell proliferation. CTCF binding is highly sensitive to DNA methylation, and this enhancer-docking site, which is hypomethylated in diverse cancers, can be inactivated through epigenetic editing with dCas9-DNMT. Similar enhancer-docking sites occur at other genes, including genes with prominent roles in multiple cancers, suggesting a mechanism by which tumor cell oncogenes can generally hijack enhancers. These results provide insights into mechanisms that allow a single target gene to be regulated by diverse enhancer elements in different cell types.

Dynamic Enhancer DNA Methylation as Basis for Transcriptional and Cellular Heterogeneity of ESCs

Variable levels of DNA methylation have been reported at tissue-specific differential methylation regions (DMRs) overlapping enhancers, including super-enhancers (SEs) associated with key cell identity genes, but the mechanisms responsible for this intriguing behavior are not well understood. We used allele-specific reporters at the endogenous Sox2 and Mir290 SEs in embryonic stem cells and found that the allelic DNA methylation state is dynamically switching, resulting in cell-to-cell heterogeneity. Dynamic DNA methylation is driven by the balance between DNA methyltransferases and transcription factor binding on one side and co-regulated with the Mediator complex recruitment and H3K27ac level changes at regulatory elements on the other side. DNA methylation at the Sox2 and the Mir290 SEs is independently regulated and has distinct consequences on the cellular differentiation state. Dynamic allele-specific DNA methylation at the two SEs was also seen at different stages in preimplantation embryos, revealing that methylation heterogeneity occurs in vivo.

Aberrant Transcriptional Regulation of Super-enhancers by RET Finger Protein-histone Deacetylase 1 Complex in Glioblastoma: Chemoresistance to Temozolomide

Glioblastoma (GBM), the most common primary brain tumor, is the most aggressive human cancers, with a median survival rate of only 14.6 months. Temozolomide (TMZ) is the frontline chemotherapeutic drug in GBM. Drug resistance is the predominant obstacle in TMZ therapy. Drug resistance occurs via multiple pathways such as DNA mismatch repair and base excision repair systems, by which glioma cells acquire chemoresistance to some extent (5% and 95%, respectively). Histone3 Lysin27 residue-acetylation (H3K27ac) status regulates cis-regulatory elements, which increases the likelihood of gene transcription. Histone deacetylase (HDAC) complex deacetylate lysine residues on core histones, leading to a decrease in gene transcription. In cis-regulatory element regions, complexes with HDAC repress histones by H3K27ac deacetylation. The cis-regulating and three-dimensional transcriptional mechanism is called "super-enhancer". RET finger protein (RFP) is a protein that is expressed in many kinds of cancer. RFP forms a protein complex with HDAC1. The disruption of the RFP-HDAC1 complex has resulted in increased drug sensitivity in other cancers. We conclude that the downregulation of RFP or the disruption of the RFP/HDAC1 complex leads to an increase in TMZ efficacy in glioblastoma by changing histone modifications which lead to changes in cell division, cell cycle and apoptosis.

Disruption of Broad Epigenetic Domains in PDAC Cells by HAT Inhibitors

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

Super-enhancers: novel target for pancreatic ductal adenocarcinoma

Super-enhancers (SEs) are unique areas of the genome which drive high-level of transcription and play a pivotal role in the cell physiology. Previous studies have established several important genes in cancer as SE-driven oncogenes. It is likely that oncogenes may hack the resident tissue regenerative program and interfere with SE-driven repair networks, leading to the specific pancreatic ductal adenocarcinoma (PDAC) phenotype. Here, we used ChIP-Seq to identify the presence of SE in PDAC cell lines. Differential H3K27AC marks were identified at enhancer regions of genes including c-MYC, MED1, OCT-4, NANOG, and SOX2 that can act as SE in non-cancerous, cancerous and metastatic PDAC cell lines. GZ17-6.02 affects acetylation of the genes, reduces transcription of major transcription factors, sonic hedgehog pathway proteins, and stem cell markers. In accordance with the decrease in Oct-4 expression, ChIP-Seq revealed a significant decrease in the occupancy of OCT-4 in the entire genome after GZ17-6.02 treatment suggesting the possible inhibitory effect of GZ17-6.02 on PDAC. Hence, SE genes are associated with PDAC and targeting their regulation with GZ17-6.02 offers a novel approach for treatment.