Epigenomic profiling of primary gastric adenocarcinoma reveals super-enhancer heterogeneity

Super-enhancer-driven ZFP36L1 promotes PD-L1 expression in infiltrative gastric cancer

Gastric cancer (GC) is a major cause of cancer-related mortality worldwide. Despite the widespread recognition of tumor immunotherapy in treating unresectable GC, challenges, including ineffective immunotherapy and drug resistance, persist. Therefore, understanding the regulatory mechanisms of PD-L1, particularly in the context of super-enhancers (SEs) and zinc finger protein 36 ring finger protein-like 1 (ZFP36L1) RNA-binding protein, is crucial. In this study, we performed H3K27ac Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing, investigated the heterogeneity of SEs between two GC subtypes with differential growth patterns, and revealed the immune escape signatures driven by ZFP36L1-SE in infiltrative GC through SEs inhibitors treatment. The regulation of ZFP36L1 to PD-L1 was evaluated by quantitative PCR, western blot, flow cytometry, and immunohistochemistry. Furthermore, we explored its regulatory mechanisms using a combination of molecular biology techniques, including luciferase reporter assay, GST/RNA pull-down, chromatin immunoprecipitation (ChIP)/RIP experiments, and in vivo functional assays. We demonstrated that ZFP36L1, driven by an SE, enhances IFN-γ-induced PD-L1 expression, with SPI1 identified as the specific transcription factor binding to ZFP36L1-SE. Mechanistically, ZFP36L1 binds to the adenylate uridylate-rich element in the 3' untranslated region (3'UTR) of HDAC3 mRNA, exacerbating its mRNA decay, and thereby facilitating PD-L1 abnormal transcriptional activation. Collectively, our findings provide mechanistic insights into the role of the SPI1-ZFP36L1-HDAC3-PD-L1 signaling axis in orchestrating immune escape mechanisms in GC, thereby offering valuable insights into the potential targets for immune checkpoint therapy in GC management.

CASCADES, a novel SOX2 super-enhancer-associated long noncoding RNA, regulates cancer stem cell specification and differentiation in glioblastoma

Glioblastoma is the most common primary malignant brain tumor in adults, with a median survival of just over 1 year. The failure of available treatments to achieve remission in patients with glioblastoma (GBM) has been attributed to the presence of cancer stem cells (CSCs), which are thought to play a central role in tumor development and progression and serve as a treatment-resistant cell repository capable of driving tumor recurrence. In fact, the property of "stemness" itself may be responsible for treatment resistance. In this study, we identify a novel long noncoding RNA (lncRNA), cancer stem cell-associated distal enhancer of SOX2 (CASCADES), that functions as an epigenetic regulator in glioma CSCs (GSCs). CASCADES is expressed in isocitrate dehydrogenase (IDH)-wild-type GBM and is significantly enriched in GSCs. Knockdown of CASCADES in GSCs results in differentiation towards a neuronal lineage in a cell- and cancer-specific manner. Bioinformatics analysis reveals that CASCADES functions as a super-enhancer-associated lncRNA epigenetic regulator of SOX2. Our findings identify CASCADES as a critical regulator of stemness in GSCs that represents a novel epigenetic and therapeutic target for disrupting the CSC compartment in glioblastoma.

Gastric cancer-Epidemiology, modifiable and non-modifiable risk factors, challenges and opportunities: An updated review

Gastric cancer represents a significant global health challenge due to its high mortality and incidence rates, particularly in Eastern Asia, Eastern Europe, and South America. This comprehensive review synthesizes the latest epidemiological data and explores both modifiable and non-modifiable risk factors associated with gastric cancer, aiming to delineate the multifactorial etiology of this disease. Modifiable risk factors include Helicobacter pylori infection, obesity, dietary habits, smoking and alcohol consumption, whereas nonmodifiable factors comprise genetic predispositions, age, family history and male gender. The interplay of these factors significantly impacts the risk and progression of gastric cancer, suggesting potential preventive strategies. The challenges in treating gastric cancer are considerable, largely because of the late-stage diagnosis and the heterogeneity of the disease, which complicate effective treatment regimens. Current treatment strategies involve a combination of surgery, chemotherapy, radiotherapy, and targeted therapies. The FLOT regimen (5-FU, Leucovorin, Oxaliplatin and Docetaxel) is now a standard for resectable cases in Europe and the US, showing superior survival and response rates over ECF and ECX regimens. For HER2-positive gastric cancer, trastuzumab combined with chemotherapy improves overall survival, as demonstrated by the ToGA trial. Additionally, immune checkpoint inhibitors like pembrolizumab and nivolumab offer promising results. However, the five-year survival rate remains low, underscoring the urgency for improved therapeutic approaches. Recent advancements in molecular biology and cancer genomics have begun to pave the way for personalized medicine in gastric cancer care, focusing on molecular targeted therapies and immunotherapy. This review also highlights the critical need for better screening methods that could facilitate early detection and treatment, potentially improving the prognosis. By integrating epidemiological insights with new therapeutic strategies, this article aims to thoroughly understand of gastric cancer's dynamics and outline a framework for future research and clinical management, advocating for a multidisciplinary approach to tackle this formidable disease.

Single-cell RNA sequencing reveals the differentiation and regulation of endplate cells in human intervertebral disc degeneration

Low back pain (LBP) is largely attributed to intervertebral disc degeneration (IVDD), of which the endplate changes are an important component. However, the alterations in cell fate and properties within the endplates during degeneration remain unknown. Here, we firstly performed the single-cell RNA-sequencing analysis (scRNA-seq) of the cells focusing on degenerative human endplates. By unsupervised clustering of the 8,534 single-cell based on the gene expression, we identified nine distinct cell types. We employed Gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, and the single-cell regulatory network inference and clustering (SCENIC) to determine the enriched pathways and transcriptional activities across seven chondrocyte subpopulations. Furthermore, two cell fates of chondrocyte differentiation were found by trajectory analysis, one was enriched in inflammation-related genes, and the other was related to extracellular matrix (ECM). Additionally, the intercellular interactions of macrophages (MA) and chondrocytes, T cells/natural killer cells (T/NK) and chondrocytes were examined by ligand-receptor pairs analysis, showing the important regulative function of FN1 from MA and CD74 from T/NK during endplate degeneration. Overall, our findings provide novel perspectives on the endplate degeneration at the single-cell level and a whole-transcriptome size.

A new high-throughput screening methodology for the discovery of cancer-testis antigen using multi-omics data

BACKGROUND

Cancer/testis antigens (CTAs), also known as tumor-specific antigens (TSAs) are specifically expressed in cancer cells and exhibit high immunogenicity, making them promising targets for immunotherapy and cancer vaccines.

METHODS

A new integrated high-throughput screening methodology for CTAs was proposed in this study through combining DNA methylation and RNA sequencing data. Briefly, the genes with increased transcript level and decreased DNA methylation were identified by multi-omics analysis. RNA sequencing studies in cell lines exposed to DNA methyltransferase (DNMT) inhibitors were performed to validate the inherent causal relationship between DNA hypomethylation and gene expression upregulation.

RESULTS

We proposed a new integrated high-throughput screening methodology for identification of CTAs using multi-omics analysis. In addition, we tested the feasibility of this method using gastric cancer (GC) as an example. In GC, we identified over 2000 primary candidate CTAs and ultimately identified 20 CTAs with significant tissue-specificity, including a testis-specific serine protease TESSP1/PRSS41. Integrated analysis confirmed that PRSS41 expression was reactivated in gastrointestinal cancers by promoter DNA hypomethylation at the CpG site (cg08104780). Additionally, DNA hypomethylation of PRSS41 predicted a poor prognosis in GC.

CONCLUSION

We propose a new high-throughput screening method for the identification of CTAs in cancer and validate its effectiveness. Our work emphasizes that serine protease PRSS41 is a novel TSA that is reactivated in GC due to promoter DNA hypomethylation.

Integrated enhancer regulatory network by enhancer-promoter looping in gastric cancer

Enhancer cis-regulatory elements play critical roles in gene regulation at many stages of cell growth. Enhancers in cancer cells also regulate the transcription of oncogenes. In this study, we performed a comprehensive analysis of long-range chromatin interactions, histone modifications, chromatin accessibility and expression in two gastric cancer (GC) cell lines compared to normal gastric epithelial cells. We found that GC-specific enhancers marked by histone modifications can activate a population of genes, including some oncogenes, by interacting with their proximal promoters. In addition, motif analysis of enhancer-promoter interacting enhancers showed that GC-specific transcription factors are enriched. Among them, we found that MYB is crucial for GC cell growth and activated by the enhancer with an enhancer-promoter loop and TCF7 upregulation. Clinical GC samples showed epigenetic activation of enhancers at the MYB locus and significant upregulation of TCF7 and MYB, regardless of molecular GC subtype and clinicopathological factors. Single-cell RNA sequencing of gastric mucosa with intestinal metaplasia showed high expression of TCF7 and MYB in intestinal stem cells. When we inactivated the loop-forming enhancer at the MYB locus using CRISPR interference (dCas9-KRAB), GC cell growth was significantly inhibited. In conclusion, we identified MYB as an oncogene activated by a loop-forming enhancer and contributing to GC cell growth.

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.

Histone 3 lysine 9 acetylation-specific reprogramming regulates esophageal squamous cell carcinoma progression and metastasis

Dysregulation of histone acetylation is widely implicated in tumorigenesis, yet its specific roles in the progression and metastasis of esophageal squamous cell carcinoma (ESCC) remain unclear. Here, we profiled the genome-wide landscapes of H3K9ac for paired adjacent normal (Nor), primary ESCC (EC) and metastatic lymph node (LNC) esophageal tissues from three ESCC patients. Compared to H3K27ac, we identified a distinct epigenetic reprogramming specific to H3K9ac in EC and LNC samples relative to Nor samples. This H3K9ac-related reprogramming contributed to the transcriptomic aberration of targeting genes, which were functionally associated with tumorigenesis and metastasis. Notably, genes with gained H3K9ac signals in both primary and metastatic lymph node samples (common-gained gene) were significantly enriched in oncogenes. Single-cell RNA-seq analysis further revealed that the corresponding top 15 common-gained genes preferred to be enriched in mesenchymal cells with high metastatic potential. Additionally, in vitro experiment demonstrated that the removal of H3K9ac from the common-gained gene MSI1 significantly downregulated its transcription, resulting in deficiencies in ESCC cell proliferation and migration. Together, our findings revealed the distinct characteristics of H3K9ac in esophageal squamous cell carcinogenesis and metastasis, and highlighted the potential therapeutic avenue for intervening ESCC through epigenetic modulation via H3K9ac.

Decoding gene regulatory circuitry underlying TNBC chemoresistance reveals biomarkers for therapy response and therapeutic targets

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype characterised by extensive intratumoral heterogeneity, high rates of metastasis and chemoresistance, leading to poor clinical outcomes. Despite progress, the mechanistic basis of chemotherapy resistance in TNBC patients remains poorly understood. Here, leveraging single-cell transcriptome datasets of matched longitudinal TNBC chemoresponsive and chemoresistant patient cohorts, we unravel distinct cell subpopulations intricately associated with chemoresistance and the signature genes defining these populations. Notably, using genome-wide mapping of the H3K27ac mark, we show that the expression of these chemoresistance genes is driven via a set of TNBC super-enhancers and associated transcription factor networks across TNBC subtypes. Furthermore, genetic screens reveal that a subset of these transcription factors is essential for the survival of TNBC cells, and their loss increases sensitivity to chemotherapeutic agents. Overall, our study has revealed epigenetic and transcription factor networks underlying chemoresistance and suggests novel avenues to stratify and improve the treatment of patients with a high risk of developing resistance.

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.

Hepatocyte nuclear factor 4 a (HNF4α): A perspective in cancer

HNF4α, a transcription factor, plays a vital role in regulating functional genes and biological processes. Its alternative splicing leads to various transcript variants encoding different isoforms. The spotlight has shifted towards the extensive discussion on tumors interplayed withHNF4α abnormalities. Aberrant HNF4α expression has emerged as sentinel markers of epigenetic shifts, casting reverberations upon downstream target genes and intricate signaling pathways, most notably with cancer. This review provides a comprehensive overview of HNF4α's involvement in tumor progression and metastasis, elucidating its role and underlying mechanisms.

Super Enhancer-Regulated LncRNA LINC01089 Induces Alternative Splicing of DIAPH3 to Drive Hepatocellular Carcinoma Metastasis

Hepatocellular carcinoma (HCC) is one of the most lethal neoplasms and has a 5-year survival rate of only 18% in patients with metastatic diseases. Epigenetic modifiers and alterations, including histone modifications, long noncoding RNAs (lncRNA), RNA alternative splicing, and N6-methyladenosine (m6A) modification, are key regulators of HCC development, highlighting the importance of understanding the cross-talk between these biological processes. In the current study, we identified LINC01089 as a super enhancer (SE)-driven lncRNA that promotes epithelial-mesenchymal transition (EMT), migration, invasion, and metastasis of HCC cells in vivo and in vitro. The transcription factor E2F1 bound to a LINC01089 SE, promoting LINC01089 transcription and overexpression. LINC01089 interacted with heterogeneous nuclear ribonucleoprotein M (hnRNPM) and led to hnRNPM-mediated skipping of DIAPH3 exon 3. Knockdown of LINC01089 increased the inclusion of DIAPH3 exon 3, which contains an important m6A-modification site that is recognized by IGF2BP3 to increase DIAPH3 mRNA stability. Thus, LINC01089 loss increased DIAPH3 protein levels, which suppressed the ERK/Elk1/Snail axis and inhibited EMT of HCC cells. In conclusion, this study revealed cross-talk between different epigenetics modifiers and alterations that drives HCC progression and identified LINC01089 as a potential prognostic marker and therapeutic target for HCC.

SIGNIFICANCE

LINC01089 is a super enhancer-driven long noncoding RNA that induces ERK signaling and epithelial-mesenchymal transition by regulating DIAPH3 alternative splicing that blocks N6-methyladenosine-mediated mRNA stabilization, establishing an epigenetic network that promotes hepatocellular carcinoma metastasis.

Super-enhancer driven SOX2 promotes tumor formation by chromatin re-organization in nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) is a malignant head and neck cancer with a high incidence in Southern China and Southeast Asia. Patients with remote metastasis and recurrent NPC have poor prognosis. Thus, a better understanding of NPC pathogenesis may identify novel therapies to address the unmet clinical needs.

METHODS

H3K27ac ChIP-seq and HiChIP was applied to understand the enhancer landscapes and the chromosome interactions. Whole genome sequencing was conducted to analyze the relationship between genomic variations and epigenetic dysregulation. CRISPRi and JQ1 treatment were used to evaluate the transcriptional regulation of SOX2 SEs. Colony formation assay, survival analysis and in vivo subcutaneous patient-derived xenograft assays were applied to explore the function and clinical relevance of SOX2 in NPC.

FINDINGS

We globally mapped the enhancer landscapes and generated NPC enhancer connectomes, linking NPC specific enhancers and SEs. We found five overlapped genes, including SOX2, among super-enhancer regulated genes, survival related genes and NPC essential genes. The mRNA expression of SOX2 was repressed when applying CRISPRi targeting different SOX2 SEs or JQ1 treatment. Next, we identified a genetic variation (Chr3:181422197, G > A) in SOX2 SE which is correlated with higher expression of SOX2 and poor survival. In addition, SOX2 was highly expressed in NPC and is correlated with short survival in patients with NPC. Knock-down of SOX2 suppressed tumor growth in vitro and in vivo.

INTERPRETATION

Our study demonstrated the super-enhancer landscape with chromosome interactions and identified super-enhancer driven SOX2 promotes tumorigenesis, suggesting that SOX2 is a potential therapeutic target for patients with NPC.

FUNDING

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

LINE-1 repression in Epstein-Barr virus-associated gastric cancer through viral-host genome interaction

Long INterspersed Element 1 (LINE-1 or L1) acts as a major remodeling force in genome regulation and evolution. Accumulating evidence shows that virus infection impacts L1 expression, potentially impacting host antiviral response and diseases. The underlying regulation mechanism is unclear. Epstein-Barr virus (EBV), a double-stranded DNA virus linked to B-cell and epithelial malignancies, is known to have viral-host genome interaction, resulting in transcriptional rewiring in EBV-associated gastric cancer (EBVaGC). By analyzing publicly available datasets from the Gene Expression Omnibus (GEO), we found that EBVaGC has L1 transcriptional repression compared with EBV-negative gastric cancer (EBVnGC). More specifically, retrotransposition-associated young and full-length L1s (FL-L1s) were among the most repressed L1s. Epigenetic alterations, especially increased H3K9me3, were observed on FL-L1s. H3K9me3 deposition was potentially attributed to increased TASOR expression, a key component of the human silencing hub (HUSH) complex for H3K9 trimethylation. The 4C- and HiC-seq data indicated that the viral DNA interacted in the proximity of the TASOR enhancer, strengthening the loop formation between the TASOR enhancer and its promoter. These results indicated that EBV infection is associated with increased H3K9me3 deposition, leading to L1 repression. This study uncovers a regulation mechanism of L1 expression by chromatin topology remodeling associated with viral-host genome interaction in EBVaGC.

In silico enhancer mining reveals SNS-032 and EHMT2 inhibitors as therapeutic candidates in high-grade serous ovarian cancer

BACKGROUND

Epigenomic dysregulation has been linked to solid tumour malignancies, including ovarian cancers. Profiling of re-programmed enhancer locations associated with disease has the potential to improve stratification and thus therapeutic choices. Ovarian cancers are subdivided into histological subtypes that have significant molecular and clinical differences, with high-grade serous carcinoma representing the most common and aggressive subtype.

METHODS

We interrogated the enhancer landscape(s) of normal ovary and subtype-specific ovarian cancer states using publicly available data. With an initial focus on H3K27ac histone mark, we developed a computational pipeline to predict drug compound activity based on epigenomic stratification. Lastly, we substantiated our predictions in vitro using patient-derived clinical samples and cell lines.

RESULTS

Using our in silico approach, we highlighted recurrent and privative enhancer landscapes and identified the differential enrichment of a total of 164 transcription factors involved in 201 protein complexes across the subtypes. We pinpointed SNS-032 and EHMT2 inhibitors BIX-01294 and UNC0646 as therapeutic candidates in high-grade serous carcinoma, as well as probed the efficacy of specific inhibitors in vitro.

CONCLUSION

Here, we report the first attempt to exploit ovarian cancer epigenomic landscapes for drug discovery. This computational pipeline holds enormous potential for translating epigenomic profiling into therapeutic leads.

Super-enhancer landscape rewiring in cancer: The epigenetic control at distal sites

Super-enhancers evolve as elements at the top of the hierarchical control of gene expression. They are important end-gatherers of signaling pathways that control stemness, differentiation or adaptive responses. Many epigenetic regulations focus on these regions, and not surprisingly, during the process of tumorigenesis, various alterations can account for their dysfunction. Super-enhancers are emerging as key drivers of the aberrant gene expression landscape that sustain the aggressiveness of cancer cells. In this review, we will describe and discuss about the structure of super-enhancers, their epigenetic regulation, and the major changes affecting their functionality in cancer.

Epigenomic mapping identifies an enhancer repertoire that regulates cell identity in bladder cancer through distinct transcription factor networks

Muscle-invasive bladder cancer (BLCA) is an aggressive disease. Consensus BLCA transcriptomic subtypes have been proposed, with two major Luminal and Basal subgroups, presenting distinct molecular and clinical characteristics. However, how these distinct subtypes are regulated remains unclear. We hypothesized that epigenetic activation of distinct super-enhancers could drive the transcriptional programs of BLCA subtypes. Through integrated RNA-sequencing and epigenomic profiling of histone marks in primary tumours, cancer cell lines, and normal human urothelia, we established the first integrated epigenetic map of BLCA and demonstrated the link between subtype and epigenetic control. We identified the repertoire of activated super-enhancers and highlighted Basal, Luminal and Normal-associated SEs. We revealed super-enhancer-regulated networks of candidate master transcription factors for Luminal and Basal subgroups including FOXA1 and ZBED2, respectively. FOXA1 CRISPR-Cas9 mutation triggered a shift from Luminal to Basal phenotype, confirming its role in Luminal identity regulation and induced ZBED2 overexpression. In parallel, we showed that both FOXA1 and ZBED2 play concordant roles in preventing inflammatory response in cancer cells through STAT2 inhibition. Our study furthers the understanding of epigenetic regulation of muscle-invasive BLCA and identifies a co-regulated network of super-enhancers and associated transcription factors providing potential targets for the treatment of this aggressive disease.

eRNA profiling uncovers the enhancer landscape of oesophageal adenocarcinoma and reveals new deregulated pathways

Cancer is driven by both genetic and epigenetic changes that impact on gene expression profiles and the resulting tumourigenic phenotype. Enhancers are transcriptional regulatory elements that are key to our understanding of how this rewiring of gene expression is achieved in cancer cells. Here, we have harnessed the power of RNA-seq data from hundreds of patients with oesophageal adenocarcinoma (OAC) or its precursor state Barrett's oesophagus coupled with open chromatin maps to identify potential enhancer RNAs and their associated enhancer regions in this cancer. We identify ~1000 OAC-specific enhancers and use these data to uncover new cellular pathways that are operational in OAC. Among these are enhancers for JUP, MYBL2, and CCNE1, and we show that their activity is required for cancer cell viability. We also demonstrate the clinical utility of our dataset for identifying disease stage and patient prognosis. Our data therefore identify an important set of regulatory elements that enhance our molecular understanding of OAC and point to potential new therapeutic directions.

Repression of DERL3 via DNA methylation by Epstein-Barr virus latent membrane protein 1 in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is Epstein-Barr virus (EBV)-associated invasive malignancy. Increasing evidence indicates that epigenetic abnormalities, including DNA methylation, play important roles in the development of NPC. In particular, the EBV principal oncogene, latent membrane protein 1 (LMP1), is considered a key factor in inducing aberrant DNA methylation of several tumour suppressor genes in NPC, although the mechanism remains unclear. Herein, we comprehensively analysed the methylome data of Infinium BeadArray from 51 NPC and 52 normal nasopharyngeal tissues to identify LMP1-inducible methylation genes. Using hierarchical clustering analysis, we classified NPC into the high-methylation, low-methylation, and normal-like subgroups. We defined high-methylation genes as those that were methylated in the high-methylation subgroup only and common methylation genes as those that were methylated in both high- and low-methylation subgroups. Subsequently, we identified 715 LMP1-inducible methylation genes by observing the methylome data of the nasopharyngeal epithelial cell line with or without LMP1 expression. Because high-methylation genes were enriched with LMP1-inducible methylation genes, we extracted 95 high-methylation genes that overlapped with the LMP1-inducible methylation genes. Among them, we identified DERL3 as the most significantly methylated gene affected by LMP1 expression. DERL3 knockdown in cell lines resulted in significantly increased cell proliferation, migration, and invasion. Lower DERL3 expression was more frequently detected in the advanced T-stage NPC than in early T-stage NPC. These results indicate that DERL3 repression by DNA methylation contributes to NPC tumour progression.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities

OBJECTIVE

Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities.

DESIGN

Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition.

RESULTS

We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.

CONCLUSION

Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

Epigenetic Regulation of MAP3K8 in EBV-Associated Gastric Carcinoma

Super-enhancers (SEs) regulate gene expressions, which are critical for cell type-identity and tumorigenesis. Although genome wide H3K27ac profiling have revealed the presence of SE-associated genes in gastric cancer (GC), their roles remain unclear. In this study, ChIP-seq and HiChIP-seq experiments revealed mitogen-activated protein kinase 8 (MAP3K8) to be an SE-associated gene with chromosome interactions in Epstein-Barr virus-associated gastric carcinoma (EBVaGC) cells. CRISPRi mediated repression of the MAP3K8 SEs attenuated MAP3K8 expression and EBVaGC cell proliferation. The results were validated by treating EBVaGC cells with bromodomain and the extra-terminal motif (BET) inhibitor, OTX015. Further, functional analysis of MAP3K8 in EBVaGC revealed that silencing MAP3K8 could inhibit the cell proliferation, colony formation, and migration of EBVaGC cells. RNA-seq and pathway analysis indicated that knocking down MAP3K8 obstructed the notch signaling pathway and epithelial-mesenchymal transition (EMT) in EBVaGC cells. Further, analysis of the cancer genome atlas (TCGA) and GSE51575 databases exhibited augmented MAP3K8 expression in gastric cancer and it was found to be inversely correlated with the disease-free progression of GC. Moreover, Spearman's correlation revealed that MAP3K8 expression was positively correlated with the expressions of notch pathway and EMT related genes, such as, Notch1, Notch2, C-terminal binding protein 2 (CTBP2), alpha smooth muscle actin isotype 2 (ACTA2), transforming growth factor beta receptor 1 (TGFβR1), and snail family transcriptional repressors 1/2 (SNAI1/SNAI2) in GC. Taken together, we are the first to functionally interrogate the mechanism of SE-mediated regulation of MAP3K8 in EBVaGC cell lines.

Chromatin and noncoding RNA-mediated mechanisms of gastric tumorigenesis

Gastric cancer (GC) is one of the most common and deadly cancers in the world. It is a multifactorial disease highly influenced by environmental factors, which include radiation, smoking, diet, and infectious pathogens. Accumulating evidence suggests that epigenetic regulators are frequently altered in GC, playing critical roles in gastric tumorigenesis. Epigenetic regulation involves DNA methylation, histone modification, and noncoding RNAs. While it is known that environmental factors cause widespread alterations in DNA methylation, promoting carcinogenesis, the chromatin- and noncoding RNA-mediated mechanisms of gastric tumorigenesis are still poorly understood. In this review, we focus on discussing recent discoveries addressing the roles of histone modifiers and noncoding RNAs and the mechanisms of their interactions in gastric tumorigenesis. A better understanding of epigenetic regulation would likely facilitate the development of novel therapeutic approaches targeting specific epigenetic regulators in GC.

The Emerging Role of Super-enhancers as Therapeutic Targets in The Digestive System Tumors

Digestive system tumors include malignancies of the stomach, pancreas, colon, rectum, and the esophagus, and are associated with high morbidity and mortality. Aberrant epigenetic modifications play a vital role in the progression of digestive system tumors. The aberrant transcription of key oncogenes is driven by super-enhancers (SEs), which are characterized by large clusters of enhancers with significantly high density of transcription factors, cofactors, and epigenetic modulatory proteins. The SEs consist of critical epigenetic regulatory elements, which modulate the biological characteristics of digestive system tumors including tumor cell identity and differentiation, tumorigenesis, environmental response, immune response, and chemotherapeutic resistance. The core transcription regulatory loop of the digestive system tumors is complex and a high density of transcription regulatory complexes in the SEs and the crosstalk between SEs and the noncoding RNAs. In this review, we summarized the known characteristics and functions of the SEs in the digestive system tumors. Furthermore, we discuss the oncogenic roles and regulatory mechanisms of SEs in the digestive system tumors. We highlight the role of SE-driven genes, enhancer RNAs (eRNAs), lncRNAs, and miRNAs in the digestive system tumor growth and progression. Finally, we discuss clinical significance of the CRISPR-Cas9 gene editing system and inhibitors of SE-related proteins such as BET and CDK7 as potential cancer therapeutics.

Oncogenic LINC00857 recruits TFAP2C to elevate FAT1 expression in gastric cancer

FAT atypical cadherin 1 (FAT1) is a mutant gene frequently found in human cancers and mainly accumulates at the plasma membrane of cancer cells. Emerging evidence has implicated FAT1 in the progression of gastric cancer (GC). This study intended to identify a regulatory network related to FAT1 in GC development. Upregulated expression of FAT1 was confirmed in GC tissues, and silencing FAT1 was observed to result in suppression of GC cell oncogenic phenotypes. Mechanistic investigation results demonstrated that FAT1 upregulated AP-1 expression by phosphorylating c-JUN and c-FOS, whereas LINC00857 elevated the expression of FAT1 by recruiting a transcription factor TFAP2C. Functional experiments further suggested that LINC00857 enhanced the malignant biological characteristics of GC cells through TFAP2C-mediated promotion of FAT1. More importantly, LINC00857 silencing delayed the tumor growth and blocked epithelial-mesenchymal transition in tumor-bearing mice, which was associated with downregulated expression of TFAP2C/FAT1. To conclude, LINC00857 plays an oncogenic role in GC through regulating the TFAP2C/FAT1/AP-1 axis. Therefore, this study contributes to extended the understanding of gastric carcinogenesis and LINC00857 may serve as a therapeutic target for GC.

A local tumor microenvironment acquired super-enhancer induces an oncogenic driver in colorectal carcinoma

Tumors exhibit enhancer reprogramming compared to normal tissue. The etiology is largely attributed to cell-intrinsic genomic alterations. Here, using freshly resected primary CRC tumors and patient-matched adjacent normal colon, we find divergent epigenetic landscapes between CRC tumors and cell lines. Intriguingly, this phenomenon extends to highly recurrent aberrant super-enhancers gained in CRC over normal. We find one such super-enhancer activated in epithelial cancer cells due to surrounding inflammation in the tumor microenvironment. We restore this super-enhancer and its expressed gene, PDZK1IP1, following treatment with cytokines or xenotransplantation into nude mice, thus demonstrating cell-extrinsic etiology. We demonstrate mechanistically that PDZK1IP1 enhances the reductive capacity CRC cancer cells via the pentose phosphate pathway. We show this activation enables efficient growth under oxidative conditions, challenging the previous notion that PDZK1IP1 acts as a tumor suppressor in CRC. Collectively, these observations highlight the significance of epigenomic profiling on primary specimens.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

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.

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.

Identifying common transcriptome signatures of cancer by interpreting deep learning models

Background

Cancer is a set of diseases characterized by unchecked cell proliferation and invasion of surrounding tissues. The many genes that have been genetically associated with cancer or shown to directly contribute to oncogenesis vary widely between tumor types, but common gene signatures that relate to core cancer pathways have also been identified. It is not clear, however, whether there exist additional sets of genes or transcriptomic features that are less well known in cancer biology but that are also commonly deregulated across several cancer types.

Results

Here, we agnostically identify transcriptomic features that are commonly shared between cancer types using 13,461 RNA-seq samples from 19 normal tissue types and 18 solid tumor types to train three feed-forward neural networks, based either on protein-coding gene expression, lncRNA expression, or splice junction use, to distinguish between normal and tumor samples. All three models recognize transcriptome signatures that are consistent across tumors. Analysis of attribution values extracted from our models reveals that genes that are commonly altered in cancer by expression or splicing variations are under strong evolutionary and selective constraints. Importantly, we find that genes composing our cancer transcriptome signatures are not frequently affected by mutations or genomic alterations and that their functions differ widely from the genes genetically associated with cancer.

Conclusions

Our results highlighted that deregulation of RNA-processing genes and aberrant splicing are pervasive features on which core cancer pathways might converge across a large array of solid tumor types.

Supplementary Information

The online version contains supplementary material available at (10.1186/s13059-022-02681-3).

Immune Infiltration and Clinical Outcome of Super-Enhancer-Associated lncRNAs in Stomach Adenocarcinoma

Super-enhancers (SEs) comprise large clusters of enhancers that highly enhance gene expression. Long non-coding RNAs (lncRNAs) tend to be dysregulated in cases of stomach adenocarcinoma (STAD) and are vital for balancing tumor immunity. However, whether SE-associated lncRNAs play a role in the immune infiltration of STAD remains unknown. In the present study, we identified SE-associated lncRNAs in the H3K27ac ChIP-seq datasets from 11 tumor tissues and two cell lines. We found that the significantly dysregulated SE-associated lncRNAs were strongly correlated with immune cell infiltration through the application of six algorithms (ImmuncellAI, CIBERSORT, EPIC, quantiSeq, TIMER, and xCELL), as well as immunomodulators and chemokines. We found that the expression of SE-associated lncRNA TM4SF1-AS1 was negatively correlated with the proportion of CD8+ T cells present in STAD. TM4SF1-AS1 suppresses T cell-mediated immune killing function and predicts immune response to anti-PD1 therapy. ChIP-seq, Hi-C and luciferase assay results verified that TM4SF1-AS1 was regulated by its super-enhancer. RNA-seq data showed that TM4SF1-AS1 is involved in immune and cancer-related processes or pathways. In conclusion, SE-associated lncRNAs are involved in the tumor immune microenvironment and act as indicators of clinical outcomes in STAD. This study highlights the importance of SE-associated lncRNAs in the immune regulation of STAD.

Single-Cell Atlas of Lineage States, Tumor Microenvironment, and Subtype-Specific Expression Programs in Gastric Cancer

Gastric cancer heterogeneity represents a barrier to disease management. We generated a comprehensive single-cell atlas of gastric cancer (>200,000 cells) comprising 48 samples from 31 patients across clinical stages and histologic subtypes. We identified 34 distinct cell-lineage states including novel rare cell populations. Many lineage states exhibited distinct cancer-associated expression profiles, individually contributing to a combined tumor-wide molecular collage. We observed increased plasma cell proportions in diffuse-type tumors associated with epithelial-resident KLF2 and stage-wise accrual of cancer-associated fibroblast subpopulations marked by high INHBA and FAP coexpression. Single-cell comparisons between patient-derived organoids (PDO) and primary tumors highlighted inter- and intralineage similarities and differences, demarcating molecular boundaries of PDOs as experimental models. We complemented these findings by spatial transcriptomics, orthogonal validation in independent bulk RNA-sequencing cohorts, and functional demonstration using in vitro and in vivo models. Our results provide a high-resolution molecular resource of intra- and interpatient lineage states across distinct gastric cancer subtypes.

SIGNIFICANCE

We profiled gastric malignancies at single-cell resolution and identified increased plasma cell proportions as a novel feature of diffuse-type tumors. We also uncovered distinct cancer-associated fibroblast subtypes with INHBA-FAP-high cell populations as predictors of poor clinical prognosis. Our findings highlight potential origins of deregulated cell states in the gastric tumor ecosystem. This article is highlighted in the In This Issue feature, p. 587.

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.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Mapping the genomic diaspora of gastric cancer

Gastric cancer (GC) is a leading contributor to global cancer incidence and mortality. Pioneering genomic studies, focusing largely on primary GCs, revealed driver alterations in genes such as ERBB2, FGFR2, TP53 and ARID1A as well as multiple molecular subtypes. However, clinical efforts targeting these alterations have produced variable results, hampered by complex co-alteration patterns in molecular profiles and intra-patient genomic heterogeneity. In this Review, we highlight foundational and translational advances in dissecting the genomic cartography of GC, including non-coding variants, epigenomic aberrations and transcriptomic alterations, and describe how these alterations interplay with environmental influences, germline factors and the tumour microenvironment. Mapping of these alterations over the GC life cycle in normal gastric tissues, metaplasia, primary carcinoma and distant metastasis will improve our understanding of biological mechanisms driving GC development and promoting cancer hallmarks. On the translational front, integrative genomic approaches are identifying diverse mechanisms of GC therapy resistance and emerging preclinical targets, enabled by technologies such as single-cell sequencing and liquid biopsies. Validating these insights will require specifically designed GC cohorts, converging multi-modal genomic data with longitudinal data on therapeutic challenges and patient outcomes. Genomic findings from these studies will facilitate 'next-generation' clinical initiatives in GC precision oncology and prevention.

Gastric cancer and genomics: review of literature

Gastric cancer (GC) is a major health concern in many countries. GC is a heterogeneous disease stratified by histopathological differences. However, these variations are not used to determine GC management. Next-generation sequencing (NGS) technologies have become widely used, and cancer genomic analysis has recently revealed the relationships between various malignant tumors and genomic information. In 2014, studies using whole-exome sequencing (WES) and whole-genome sequencing (WGS) for GC revealed the entire structure of GC genomics. Genomics with NGS has been used to identify new therapeutic targets for GC. Moreover, personalized medicine to provide specific therapy for targets based on multiplex gene panel testing of tumor tissues has become of clinical use. Recently, immune checkpoint inhibitors (ICIs) have been used for GC treatment; however, their response rates are limited. To predict the anti-tumor effects of ICIs for GC and to select patients suitable for ICI treatment, genomics also provides informative data not only of tumors but also of tumor microenvironments, such as tumor-infiltrating lymphocytes. In therapeutic strategies for unresectable or recurrent malignant tumors, the target is not only the primary lesion but also metastatic lesions, and metastatic lesions are often resistant to chemotherapy. Unlike colorectal carcinoma, there is a heterogeneous status of genetic variants between the primary and metastatic lesions in GC. Liquid biopsy analysis is also helpful for predicting the genomic status of both primary and metastatic lesions. Genomics has become an indispensable tool for GC treatment and is expected to be further developed in the future.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Genome-wide profiling in colorectal cancer identifies PHF19 and TBC1D16 as oncogenic super enhancers

Colorectal cancer is one of the most common cancers in the world. Although genomic mutations and single nucleotide polymorphisms have been extensively studied, the epigenomic status in colorectal cancer patient tissues remains elusive. Here, together with genomic and transcriptomic analysis, we use ChIP-Seq to profile active enhancers at the genome wide level in colorectal cancer paired patient tissues (tumor and adjacent tissues from the same patients). In total, we sequence 73 pairs of colorectal cancer tissues and generate 147 H3K27ac ChIP-Seq, 144 RNA-Seq, 147 whole genome sequencing and 86 H3K4me3 ChIP-Seq samples. Our analysis identifies 5590 gain and 1100 lost variant enhancer loci in colorectal cancer, and 334 gain and 121 lost variant super enhancer loci. Multiple key transcription factors in colorectal cancer are predicted with motif analysis and core regulatory circuitry analysis. Further experiments verify the function of the super enhancers governing PHF19 and TBC1D16 in regulating colorectal cancer tumorigenesis, and KLF3 is identified as an oncogenic transcription factor in colorectal cancer. Taken together, our work provides an important epigenomic resource and functional factors for epigenetic studies in colorectal cancer.

Histological and mutational profile of diffuse gastric cancer: current knowledge and future challenges

Gastric cancer (GC) pathogenesis is complex and heterogeneous, reflecting morphological, molecular and genetic diversity. Diffuse gastric cancer (DGC) and intestinal gastric cancer (IGC) are the major histological types. GC may be sporadic or hereditary; sporadic GC is related to environmental and genetic low-risk factors and hereditary GC is caused by inherited high-risk mutations, so far identified only for the diffuse histotype. DGC phenotypic heterogeneity challenges the current understanding of molecular mechanisms underlying carcinogenesis. The definition of a DGC-specific mutational profile remains controversial, possibly reflecting the heterogeneity of DGC-related histological subtypes [signet-ring cell carcinoma (SRCC) and poorly cohesive carcinoma not otherwise specified (PCC-NOS)]. Indeed, DGC and DGC-related subtypes may present specific mutational profiles underlying the particularly aggressive behaviour and dismal prognosis of DGC vs IGC and PCC-NOS vs SRCC. In this systematic review, we revised the histological presentations, molecular classifications and approved therapies for gastric cancer, with a focus on DGC. We then analysed results from the most relevant studies, reporting mutational analysis data specifying mutational frequencies, and their relationship with DGC and IGC histological types, and with specific DGC subtypes (SRCC and PCC-NOS). We aimed at identifying histology-associated mutational profiles with an emphasis in DGC and its subtypes (DGC vs IGC; sporadic vs hereditary DGC; and SRCC vs PCC-NOS). We further used these mutational profiles to identify the most commonly affected molecular pathways and biological functions, and explored the clinical trials directed specifically to patients with DGC. This systematic analysis is expected to expose a DGC-specific molecular profile and shed light into potential targets for therapeutic intervention, which are currently missing.

Deciphering the molecular mechanism of the cancer formation by chromosome structural dynamics

Cancer reflects the dysregulation of the underlying gene network, which is strongly related to the 3D genome organization. Numerous efforts have been spent on experimental characterizations of the structural alterations in cancer genomes. However, there is still a lack of genomic structural-level understanding of the temporal dynamics for cancer initiation and progression. Here, we use a landscape-switching model to investigate the chromosome structural transition during the cancerization and reversion processes. We find that the chromosome undergoes a non-monotonic structural shape-changing pathway with initial expansion followed by compaction during both of these processes. Furthermore, our analysis reveals that the chromosome with a more expanding structure than those at both the normal and cancer cell during cancerization exhibits a sparse contact pattern, which shows significant structural similarity to the one at the embryonic stem cell in many aspects, including the trend of contact probability declining with the genomic distance, the global structural shape geometry and the spatial distribution of loci on the chromosome. In light of the intimate structure-function relationship at the chromosomal level, we further describe the cell state transition processes by the chromosome structural changes, suggesting an elevated cell stemness during the formation of the cancer cells. We show that cell cancerization and reversion are highly irreversible processes in terms of the chromosome structural transition pathways, spatial repositioning of chromosomal loci and hysteresis loop of contact evolution analysis. Our model draws a molecular-scale picture of cell cancerization from the chromosome structural perspective. The process contains initial reprogramming towards the stem cell followed by the differentiation towards the cancer cell, accompanied by an initial increase and subsequent decrease of the cell stemness.

Hepatocyte nuclear factor HNF1A is a potential regulator in shaping the super-enhancer landscape in colorectal cancer liver metastasis

Super-enhancers (SEs) play essential roles in colorectal cancer (CRC) progression. However, how the SE landscape is orchestrated by transcriptional regulators and evolves is not clear. Using de novo motif analysis, we show that the hepatocyte nuclear factor 1 (HNF1)-binding motif is enriched in SEs in cell lines derived from liver metastases, but not in those from primary tumors. This finding was further validated by extending the method to pancreatic cancer and a pair of isogenic CRC lines. Next, we revealed HNF1-alpha (HNF1A) was majorly expressed and upregulated in CRC liver metastatic cell lines. Clinically, HNF1A was remarkably upregulated in synchronous liver metastases as compared to localized tumors. Collectively, our study implicates HNF1A as a key regulator in shaping the SE landscape in CRC liver metastasis.

No Need to Stick Together to Be Connected: Multiple Types of Enhancers' Networking

Simple Summary

Transcription regulation programs require the functional interaction of distal and proximal regulatory regions, interacting by specific 3D chromatin configurations. Enhancers are cis-acting regulatory elements able to promote gene expression regardless their orientation and distance from the transcription starting site. Their systematic mapping by genome-wide chromatin profiling and chromosome conformation analysis, combined with the development of gene-editing approaches to modulate their function, revealed that many enhancers work together to fine-tune the expression of their target genes. This review aim to describe the functions of different types of enhancers and the modalities of enhancers’ interaction, focusing on their role in the regulation of complex biological processes like cancer development.

Abstract

The control of gene expression at a transcriptional level requires a widespread landscape of regulatory elements. Central to these regulatory circuits are enhancers (ENHs), which are defined as cis-acting DNA elements able to increase the transcription of a target gene in a distance- and orientation-independent manner. ENHs are not independent functional elements but work in a complex and dynamic cooperative network, constituting the building blocks of multimodular domains of gene expression regulation. The information from each of these elements converges on the target promoter, contributing to improving the precision and sharpness of gene modulation. ENHs’ interplay varies in its nature and extent, ranging from an additive to redundant effect depending on contexts. Moving from super-enhancers that drive the high expression levels of identity genes, to shadow-enhancers, whose redundant functions contribute to buffering the variation in gene expression, this review aims to describe the different modalities of ENHs’ interaction and their role in the regulation of complex biological processes like cancer development.

Integrative epigenomic and high-throughput functional enhancer profiling reveals determinants of enhancer heterogeneity in gastric cancer

BACKGROUND

Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity-however, most predicted enhancer regions remain to be functionally tested.

METHODS

We analyzed 132 epigenomic histone modification profiles of 18 primary gastric cancer (GC) samples, 18 normal gastric tissues, and 28 GC cell lines using Nano-ChIP-seq technology. We applied Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to assess functional enhancer activity. An Activity-by-contact (ABC) model was employed to explore the effects of histone acetylation and CapSTARR-seq levels on enhancer-promoter interactions.

RESULTS

We report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are GC-associated in vivo (> 50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying CapSTARR-seq to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers, even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. We identified cancer-relevant genes (ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity.

CONCLUSIONS

Our results indicate that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, providing insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

Oncogenic enhancers drive esophageal squamous cell carcinogenesis and metastasis

The role of cis-elements and their aberrations remains unclear in esophageal squamous cell carcinoma (ESCC, further abbreviated EC). Here we survey 28 H3K27ac-marked active enhancer profiles and 50 transcriptomes in primary EC, metastatic lymph node cancer (LNC), and adjacent normal (Nor) esophageal tissues. Thousands of gained or lost enhancers and hundreds of altered putative super-enhancers are identified in EC and LNC samples respectively relative to Nor, with a large number of common gained or lost enhancers. Moreover, these differential enhancers contribute to the transcriptomic aberrations in ECs and LNCs. We also reveal putative driver onco-transcription factors, depletion of which diminishes cell proliferation and migration. The administration of chemical inhibitors to suppress the predicted targets of gained super-enhances reveals HSP90AA1 and PDE4B as potential therapeutic targets for ESCC. Thus, our epigenomic profiling reveals a compendium of reprogrammed cis-regulatory elements during ESCC carcinogenesis and metastasis for uncovering promising targets for cancer treatment.

ANLN Enhances Triple-Negative Breast Cancer Stemness Through TWIST1 and BMP2 and Promotes its Spheroid Growth

ANLN is frequently upregulated in triple-negative breast cancer (TNBC) and its high expression in tumors are significantly associated with poor survival and recurrence, thereby it has been proposed to function as a prognostic marker for breast cancer. However, the specific function and molecular mechanisms by which ANLN promotes TNBC tumorigenesis remain elusive. Using multiomic profiling, we recently uncovered ANLN as a TNBC-specific gene driven by super-enhancer. Here, by Crispr/Cas9 editing, we showed that knockout of ANLN inhibits spheroid growth of TNBC. Interestingly, its effect on cell proliferation in 2D cultures is minimal. ANLN depletion inhibits mammosphere formation and clonogenicity potently, suggesting its important function in regulating cancer stem cells (CSCs). We screened a panel of stem cell-related genes and uncovered several CSC genes regulated by ANLN. We further identify TWIST1 and BMP2 as essential genes that mediate ANLN’s function in stemness but not spheroid growth. These findings may contribute to search for effective targeted therapies to treat TNBC.

Epigenetic alteration contributes to the transcriptional reprogramming in T-cell prolymphocytic leukemia

T cell prolymphocytic leukemia (T-PLL) is a rare disease with aggressive clinical course. Cytogenetic analysis, whole-exome and whole-genome sequencing have identified primary structural alterations in T-PLL, including inversion, translocation and copy number variation. Recurrent somatic mutations were also identified in genes encoding chromatin regulators and those in the JAK-STAT signaling pathway. Epigenetic alterations are the hallmark of many cancers. However, genome-wide epigenomic profiles have not been reported in T-PLL, limiting the mechanistic study of its carcinogenesis. We hypothesize epigenetic mechanisms also play a key role in T-PLL pathogenesis. To systematically test this hypothesis, we generated genome-wide maps of regulatory regions using H3K4me3 and H3K27ac ChIP-seq, as well as RNA-seq data in both T-PLL patients and healthy individuals. We found that genes down-regulated in T-PLL are mainly associated with defense response, immune system or adaptive immune response, while up-regulated genes are enriched in developmental process, as well as WNT signaling pathway with crucial roles in cell fate decision. In particular, our analysis revealed a global alteration of regulatory landscape in T-PLL, with differential peaks highly enriched for binding motifs of immune related transcription factors, supporting the epigenetic regulation of oncogenes and genes involved in DNA damage response and T-cell activation. Together, our work reveals a causal role of epigenetic dysregulation in T-PLL.

Defining super-enhancer landscape in triple-negative breast cancer by multiomic profiling

Breast cancer is a heterogeneous disease, affecting over 3.5 million women worldwide, yet the functional role of cis-regulatory elements including super-enhancers in different breast cancer subtypes remains poorly characterized. Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer with a poor prognosis. Here we apply integrated epigenomic and transcriptomic profiling to uncover super-enhancer heterogeneity between breast cancer subtypes, and provide clinically relevant biological insights towards TNBC. Using CRISPR/Cas9-mediated gene editing, we identify genes that are specifically regulated by TNBC-specific super-enhancers, including FOXC1 and MET, thereby unveiling a mechanism for specific overexpression of the key oncogenes in TNBC. We also identify ANLN as a TNBC-specific gene regulated by super-enhancer. Our studies reveal a TNBC-specific epigenomic landscape, contributing to the dysregulated oncogene expression in breast tumorigenesis.

MLL3 suppresses tumorigenesis through regulating TNS3 enhancer activity

MLL3 is a histone H3K4 methyltransferase that is frequently mutated in cancer, but the underlying molecular mechanisms remain elusive. Here, we found that MLL3 depletion by CRISPR/sgRNA significantly enhanced cell migration, but did not elevate the proliferation rate of cancer cells. Through RNA-Seq and ChIP-Seq approaches, we identified TNS3 as the potential target gene for MLL3. MLL3 depletion caused downregulation of H3K4me1 and H3K27ac on an enhancer ~ 7 kb ahead of TNS3. 3C assay indicated the identified enhancer interacts with TNS3 promoter and repression of enhancer activity by dCas9-KRAB system impaired TNS3 expression. Exogenous expression of TNS3 in MLL3 deficient cells completely blocked the enhanced cell migration phenotype. Taken together, our study revealed a novel mechanism for MLL3 in suppressing cancer, which may provide novel targets for diagnosis or drug development.

MethylationToActivity: a deep-learning framework that reveals promoter activity landscapes from DNA methylomes in individual tumors

Although genome-wide DNA methylomes have demonstrated their clinical value as reliable biomarkers for tumor detection, subtyping, and classification, their direct biological impacts at the individual gene level remain elusive. Here we present MethylationToActivity (M2A), a machine learning framework that uses convolutional neural networks to infer promoter activities based on H3K4me3 and H3K27ac enrichment, from DNA methylation patterns for individual genes. Using publicly available datasets in real-world test scenarios, we demonstrate that M2A is highly accurate and robust in revealing promoter activity landscapes in various pediatric and adult cancers, including both solid and hematologic malignant neoplasms.