Search

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. Agency for Science, Technology and Research (A*STAR) Ministry of Education (MOE) National Medical Research Council (NMRC) National Research Foundation (NRF) This study was supported by National Medical Research Council grants NMRC/STaR/0026/2015, MOH-000967-00 and OFLCG18May-0003, and A*ccelerate GAP fund ETPL/15-R15 GAP-0021 (to P. Tan) and MOE tier 2 grant (MOE2017-T2-1-105) and NMRC CS-IRG grant (NMRC/CIRG/1481/2017 to S. Li) and SCISSOR (A*STAR IAF-PP) H18/01/a0/020. Funding was also provided by Cancer Science Institute of Singapore, NUS, under the National Research Foundation Singapore and the Singapore Ministry of Education under its Research Centers of Excellence initiative, and block funding from Duke-NUS Medical School. A.M. Nargund is also supported by St. Baldrick's Foundation Research Award. R. Sundar is supported by a National Medical Research Council (NMRC) Fellowship, Singapore. M.J. Fullwood is supported by the RNA Biology Center at the Cancer Science Institute of Singapore, NUS, as part of funding under the Singapore Ministry of Education Academic Research Fund Tier 3 awarded to Daniel Tenen (MOE2014-T3-1-006) and the National Research Foundation Singapore and the Singapore Ministry of Education under its Research Centers of Excellence initiative.

ObjectiveGastric 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.DesignTranscriptomic 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.ResultsWe 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, withTEAD1as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlightedNUAK1kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.ConclusionOur 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.

Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. We analyzed 1256 gastric samples (1152 IMs) from 692 subjects through a prospective 10-year study. We identified 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9), largely occurring as small clonal events. Analysis of clonal dynamics between and within subjects, and also longitudinally across time, revealed that IM clones are likely transient but increase in size upon progression to dysplasia, with eventual transmission of somatic events to paired GCs. Single-cell and spatial profiling highlighted changes in tissue ecology and lineage heterogeneity in IM, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling revealed expression-based molecular subtypes of IM, including a body-resident “pseudoantralized” subtype associated with incomplete histology, antral/intestinal cell types,ARID1Amutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical- genomic models outperform clinical-only models in predicting IMs likely to progress. Our results raise opportunities for GC precision prevention and interception by highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression.

Supplementary text S1. Supplementary Figure 1: Chromatin Profiles of Primary GC. Supplementary Figure 2: Epithelial features of GC promoters. Supplementary Figure 3: GC Somatic Promoter Features. Supplementary Figure 4: Association of Somatic Promoters with Gene Expression in GC and Other Tumor Types. Supplementary Figure 5: Changes in DNA methylation at CpG island containing promoters. Supplementary Figure 6: Expression distribution of alternative and canonical isoforms. Supplementary Figure 7: Characterization of RASA3 Isoform. Supplementary Figure 8: Characterization of MET Isoforms. Supplementary Figure 9: Immunogenicity of N-terminal peptides. Supplementary Figure 10: Immunogenicity Assay and Nanostring Profiling. Supplementary Figure 11: Functional Assessment of Peptide Immunogenicity. Supplementary Figure 12: EZH2 Inhibition. Supplementary Figure 13: Unannotated somatic promoters.

Supplementary Figure 5

Legends for Supplementary Figures and Tables

Supplementary Figure 4

Supplementary Figure 1

Supplementary Figure 3

Supplementary Figure 2

Supplementary Tables S1-S5

Gastric cancer cases are often diagnosed at an advanced stage with poor prognosis. Platinum-based chemotherapy has been internationally accepted as first-line therapy for inoperable or metastatic gastric cancer. To achieve greater benefits, selection of patients eligible for this treatment is critical. Although gene expression profiling has been widely used as a genomic classifier to identify molecular subtypes of gastric cancer and to stratify patients for different chemotherapy regimens, its prediction accuracy can be improved. Adenosine-to-inosine (A-to-I) RNA editing has emerged as a new player contributing to gastric cancer development and progression, offering potential clinical utility for diagnosis and treatment. Using a systematic computational approach followed by both in vitro validations and in silico validations in The Cancer Genome Atlas (TCGA), we conducted a transcriptome-wide RNA editing analysis of a cohort of 104 patients with advanced gastric cancer and identified an RNA editing (GCRE) signature to guide gastric cancer chemotherapy. RNA editing events stood as a prognostic and predictive biomarker in advanced gastric cancer. A GCRE score based on the GCRE signature consisted of 50 editing sites associated with 29 genes, predicting response to chemotherapy with a high accuracy (84%). Of note, patients demonstrating higher editing levels of this panel of sites presented a better overall response. Consistently, gastric cancer cell lines with higher editing levels showed higher chemosensitivity. Applying the GCRE score on TCGA dataset confirmed that responders had significantly higher levels of editing in advanced gastric cancer. Overall, this newly defined GCRE signature reliably stratifies patients with advanced gastric cancer and predicts response from chemotherapy.Significance:This study describes a novel A-to-I RNA editing signature as a prognostic and predictive biomarker in advanced gastric cancer, providing a new tool to improve patient stratification and response to therapy.

Supplementary Figure 6

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.

ObjectivesEpigenomic alterations in cancer interact with the immune microenvironment to dictate tumour evolution and therapeutic response. We aimed to study the regulation of the tumour immune microenvironment through epigenetic alternate promoter use in gastric cancer and to expand our findings to other gastrointestinal tumours.DesignAlternate promoter burden (APB) was quantified using a novel bioinformatic algorithm (proActiv) to infer promoter activity from short-read RNA sequencing and samples categorised into APBhigh, APBint and APBlow. Single-cell RNA sequencing was performed to analyse the intratumour immune microenvironment. A humanised mouse cancer in vivo model was used to explore dynamic temporal interactions between tumour kinetics, alternate promoter usage and the human immune system. Multiple cohorts of gastrointestinal tumours treated with immunotherapy were assessed for correlation between APB and treatment outcomes.ResultsAPBhigh gastric cancer tumours expressed decreased levels of T-cell cytolytic activity and exhibited signatures of immune depletion. Single-cell RNAsequencing analysis confirmed distinct immunological populations and lower T-cell proportions in APBhigh tumours. Functional in vivo studies using ‘humanised mice’ harbouring an active human immune system revealed distinct temporal relationships between APB and tumour growth, with APBhigh tumours having almost no human T-cell infiltration. Analysis of immunotherapy-treated patients with GI cancer confirmed resistance of APBhigh tumours to immune checkpoint inhibition. APBhigh gastric cancer exhibited significantly poorer progression-free survival compared with APBlow (median 55 days vs 121 days, HR 0.40, 95% CI 0.18 to 0.93, p=0.032).ConclusionThese findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

ObjectiveGastric cancer (GC) is a leading cause of cancer mortality, withARID1Abeing the second most frequently mutated driver gene in GC. We sought to decipherARID1A-specific GC regulatory networks and examine therapeutic vulnerabilities arising fromARID1Aloss.DesignGenomic profiling of GC patients including a Singapore cohort (>200 patients) was performed to derive mutational signatures ofARID1Ainactivation across molecular subtypes. Single-cell transcriptomic profiles ofARID1A-mutated GCs were analysed to examine tumour microenvironmental changes arising fromARID1Aloss. Genome-wide ARID1A binding and chromatin profiles (H3K27ac, H3K4me3, H3K4me1, ATAC-seq) were generated to identify gastric-specific epigenetic landscapes regulated by ARID1A. Distinct cancer hallmarks ofARID1A-mutated GCs were converged at the genomic, single-cell and epigenomic level, and targeted by pharmacological inhibition.ResultsWe observed prevalentARID1Ainactivation across GC molecular subtypes, with distinct mutational signatures and linked to a NFKB-driven proinflammatory tumour microenvironment.ARID1A-depletion caused loss of H3K27ac activation signals atARID1A-occupied distal enhancers, but unexpectedly gain of H3K27ac at ARID1A-occupied promoters in genes such asNFKB1andNFKB2. Promoter activation inARID1A-mutated GCs was associated with enhanced gene expression, increased BRD4 binding, and reduced HDAC1 and CTCF occupancy. Combined targeting of promoter activation and tumour inflammation via bromodomain and NFKB inhibitors confirmed therapeutic synergy specific toARID1A-genomic status.ConclusionOur results suggest a therapeutic strategy forARID1A-mutated GCs targeting both tumour-intrinsic (BRD4-assocatiated promoter activation) and extrinsic (NFKB immunomodulation) cancer phenotypes.

Gastric cancer cases are often diagnosed at an advanced stage with poor prognosis. Platinum-based chemotherapy has been internationally accepted as first-line therapy for inoperable or metastatic gastric cancer. To achieve greater benefits, selection of patients eligible for this treatment is critical. Although gene expression profiling has been widely used as a genomic classifier to identify molecular subtypes of gastric cancer and to stratify patients for different chemotherapy regimens, its prediction accuracy can be improved. Adenosine-to-inosine (A-to-I) RNA editing has emerged as a new player contributing to gastric cancer development and progression, offering potential clinical utility for diagnosis and treatment. Using a systematic computational approach followed by both in vitro validations and in silico validations in The Cancer Genome Atlas (TCGA), we conducted a transcriptome-wide RNA editing analysis of a cohort of 104 patients with advanced gastric cancer and identified an RNA editing (GCRE) signature to guide gastric cancer chemotherapy. RNA editing events stood as a prognostic and predictive biomarker in advanced gastric cancer. A GCRE score based on the GCRE signature consisted of 50 editing sites associated with 29 genes, predicting response to chemotherapy with a high accuracy (84%). Of note, patients demonstrating higher editing levels of this panel of sites presented a better overall response. Consistently, gastric cancer cell lines with higher editing levels showed higher chemosensitivity. Applying the GCRE score on TCGA dataset confirmed that responders had significantly higher levels of editing in advanced gastric cancer. Overall, this newly defined GCRE signature reliably stratifies patients with advanced gastric cancer and predicts response from chemotherapy. Significance: This study describes a novel A-to-I RNA editing signature as a prognostic and predictive biomarker in advanced gastric cancer, providing a new tool to improve patient stratification and response to therapy.

Epstein-Barr virus (EBV) is associated with several human malignancies including 8-10% of gastric cancers (GCs). Genome-wide analysis of 3D chromatin topologies across GC lines, primary tissue and normal gastric samples revealed chromatin domains specific to EBV-positive GC, exhibiting heterochromatin-to-euchromatin transitions and long-range human-viral interactions with non-integrated EBV episomes. EBV infection in vitro suffices to remodel chromatin topology and function at EBV-interacting host genomic loci, converting H3K9me3+ heterochromatin to H3K4me1+/H3K27ac+ bivalency and unleashing latent enhancers to engage and activate nearby GC-related genes (for example TGFBR2 and MZT1). Higher-order epigenotypes of EBV-positive GC thus signify a novel oncogenic paradigm whereby non-integrative viral genomes can directly alter host epigenetic landscapes ('enhancer infestation'), facilitating proto-oncogene activation and tumorigenesis.

BackgroundEnhancers 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.ResultsAnalyzing 128 epigenomic histone modification profiles of primary GC samples, normal gastric tissues, and GC cell lines, we report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are tumor-associated in vivo (>50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying Capture-based Self-Transcribing Active Regulatory Region sequencing (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. Specifically, we identified cancer-relevant genes (e.g. 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.ConclusionsOur study indicates that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, and provides insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.

ObjectiveGenomic structural variations (SVs) causing rewiring of cis-regulatory elements remain largely unexplored in gastric cancer (GC). To identify SVs affecting enhancer elements in GC (enhancer-based SVs), we integrated epigenomic enhancer profiles revealed by paired-end H3K27ac ChIP-sequencing from primary GCs with tumour whole-genome sequencing (WGS) data (PeNChIP-seq/WGS).DesignWe applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of >200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs.ResultsPeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targeting CCNE1, a key driver of therapy resistance, occurred in 8% of patients frequently juxtaposing diverse distal enhancers to CCNE1 proximal regions. CCNE1-rearranged GCs were associated with high CCNE1 expression, disrupted CCNE1 topologically associating domain (TAD) boundaries, and novel TAD interactions in CCNE1-rearranged primary tumours. We also observed IGF2 enhancer-based SVs, previously noted in colorectal cancer, highlighting a common non-coding genetic driver alteration in gastric and colorectal malignancies.ConclusionIntegrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.

ObjectiveGastric cancer (GC) is a leading cause of cancer mortality. Previous studies have shown that hepatocyte nuclear factor-4α (HNF4α) is specifically overexpressed in GC and functionally required for GC development. In this study, we investigated, on a genome-wide scale, target genes of HNF4α and oncogenic pathways driven by HNF4α and HNF4α target genes.DesignWe performed HNF4α chromatin immunoprecipitation followed by sequencing across multiple GC cell lines, integrating HNF4α occupancy data with (epi)genomic and transcriptome data of primary GCs to define HNF4α target genes of in vitro and in vivo relevance. To investigate mechanistic roles of HNF4α and HNF4α targets, we performed cancer metabolic measurements, drug treatments and functional assays including murine xenograft experiments.ResultsGene expression analysis across 19 tumour types revealed HNF4α to be specifically upregulated in GCs. Unbiased pathway analysis revealed organic acid metabolism as the top HNF4α-regulated pathway, orthogonally supported by metabolomic analysis. Isocitrate dehydrogenase 1 (IDH1) emerged as a convergent HNF4α direct target gene regulating GC metabolism. We show that wild-type IDH1 is essential for GC cell survival, and that certain GC cells can be targeted by IDH1 inhibitors.ConclusionsOur results highlight a role for HNF4α in sustaining GC oncogenic metabolism, through the regulation of IDH1. Drugs targeting wild-type IDH1 may thus have clinical utility in GCs exhibiting HNF4α overexpression, expanding the role of IDH1 in cancer beyond IDH1/2 mutated malignancies.

The human genome contains more than 200,000 gene isoforms. However, different isoforms can be highly similar, and with an average length of 1.5kb remain difficult to study with short read sequencing. To systematically evaluate the ability to study the transcriptome at a resolution of individual isoforms we profiled 5 human cell lines with short read cDNA sequencing and Nanopore long read direct RNA, amplification-free direct cDNA, PCR-cDNA sequencing. The long read protocols showed a high level of consistency, with amplification-free RNA and cDNA sequencing being most similar. While short and long reads generated comparable gene expression estimates, they differed substantially for individual isoforms. We find that increased read length improves read-to-transcript assignment, identifies interactions between alternative promoters and splicing, enables the discovery of novel transcripts from repetitive regions, facilitates the quantification of full-length fusion isoforms and enables the simultaneous profiling of m6A RNA modifications when RNA is sequenced directly. Our study demonstrates the advantage of long read RNA sequencing and provides a comprehensive resource that will enable the development and benchmarking of computational methods for profiling complex transcriptional events at isoform-level resolution.

Additional file 25. Review history.

Additional file 1. Supplementary figures.

Transcriptional reactivation of telomerase catalytic subunit (TERT) is a frequent hallmark of cancer, occurring in 90% of human malignancies. However, specific mechanisms driving TERT reactivation remain obscure for many tumor types and in particular gastric cancer (GC), a leading cause of global cancer mortality. Here, through comprehensive genomic and epigenomic analysis of primary GCs and GC cell lines, we identified the transcription factor early B cell factor 1 (EBF1) as a TERT transcriptional repressor and inactivation of EBF1 function as a major cause of TERT upregulation. Abolishment of EBF1 function occurs through 3 distinct (epi)genomic mechanisms. First, EBF1 is epigenetically silenced via DNA methyltransferase, polycomb-repressive complex 2 (PRC2), and histone deacetylase activity in GCs. Second, recurrent, somatic, and heterozygous EBF1 DNA-binding domain mutations result in the production of dominant-negative EBF1 isoforms. Third, more rarely, genomic deletions and rearrangements proximal to the TERT promoter remobilize or abolish EBF1-binding sites, derepressing TERT and leading to high TERT expression. EBF1 is also functionally required for various malignant phenotypes in vitro and in vivo, highlighting its importance for GC development. These results indicate that multimodal genomic and epigenomic alterations underpin TERT reactivation in GC, converging on transcriptional repressors such as EBF1.

Genomic structural variations (SVs) causing rewiring ofWe applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs.PeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targetingIntegrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.