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.

Comprehensive molecular phenotyping of ARID1A-deficient gastric cancer reveals pervasive epigenomic reprogramming and therapeutic opportunities.

OBJECTIVE: Gastric cancer (GC) is a leading cause of cancer mortality, with ARID1A being the second most frequently mutated driver gene in GC. We sought to decipher ARID1A-specific GC regulatory networks and examine therapeutic vulnerabilities arising from ARID1A loss. DESIGN: Genomic profiling of GC patients including a Singapore cohort (>200 patients) was performed to derive mutational signatures of ARID1A inactivation across molecular subtypes. Single-cell transcriptomic profiles of ARID1A-mutated GCs were analysed to examine tumour microenvironmental changes arising from ARID1A loss. 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 of ARID1A-mutated GCs were converged at the genomic, single-cell and epigenomic level, and targeted by pharmacological inhibition. RESULTS: We observed prevalent ARID1A inactivation 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 at ARID1A-occupied distal enhancers, but unexpectedly gain of H3K27ac at ARID1A-occupied promoters in genes such as NFKB1 and NFKB2. Promoter activation in ARID1A-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 to ARID1A-genomic status. CONCLUSION: Our results suggest a therapeutic strategy for ARID1A-mutated GCs targeting both tumour-intrinsic (BRD4-assocatiated promoter activation) and extrinsic (NFKB immunomodulation) cancer phenotypes.

Epigenetic promoter alterations in GI tumour immune-editing and resistance to immune checkpoint inhibition.

OBJECTIVES: Epigenomic 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. DESIGN: Alternate 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. RESULTS: APBhigh 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). CONCLUSION: These findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

DNA damage signalling as an anti-cancer barrier in gastric intestinal metaplasia.

OBJECTIVE: Intestinal metaplasia (IM) is a premalignant stage that poses a greater risk for subsequent gastric cancer (GC). However, factors regulating IM to GC progression remain unclear. Previously, activated DNA damage response (DDR) signalling factors were shown to engage tumour-suppressive networks in premalignant lesions. Here, we interrogate the relationship of DDR signalling to mutational accumulation in IM lesions. DESIGN: IM biopsies were procured from the gastric cancer epidemiology programme, an endoscopic surveillance programme where biopsies have been subjected to (epi)genomic characterisation. IM samples were classified as genome-stable or genome-unstable based on their mutational burden/somatic copy-number alteration (CNA) profiles. Samples were probed for DDR signalling and cell proliferation, using the markers γH2AX and MCM2, respectively. The expression of the gastric stem cell marker, CD44v9, was also assessed. Tissue microarrays representing the GC progression spectrum were included. RESULTS: MCM2-positivity increased during GC progression, while γH2AX-positivity showed modest increase from normal to gastritis and IM stages, with further increase in GC. γH2AX levels correlated with the extent of chronic inflammation. Interestingly, genome-stable IM lesions had higher γH2AX levels underscoring a protective anti-cancer role for DDR signalling. In contrast, genome-unstable IM lesions with higher mutational burden/CNAs had lower γH2AX levels, elevated CD44v9 expression and modest promoter hypermethylation of DNA repair genes WRN, MLH1 and RAD52. CONCLUSIONS: Our data suggest that IM lesions with active DDR will likely experience a longer latency at the premalignant state until additional hits that override DDR signalling clonally expand and promote progression. These observations provide insights on the factors governing IM progression.

HNF4α pathway mapping identifies wild-type IDH1 as a targetable metabolic node in gastric cancer.

OBJECTIVE: Gastric 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. DESIGN: We 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. RESULTS: Gene 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. CONCLUSIONS: Our 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.

An LCM-based genomic analysis of SPEM, Gastric Cancer and Pyloric Gland Adenoma in an Asian cohort.

Spasmolytic polypeptide-expressing metaplasia (SPEM) and pyloric gland adenoma (PGA) in the stomach are metaplastic and neoplastic lesions, respectively, in which gastric body glands are replaced by pyloric glands. The aim of this study was to evaluate the genomic profile of SPEM and compare it with intestinal-type gastric cancer (GC) and PGA. Thirteen gastrectomies showing PGA with or without dysplasia, GC and SPEM were retrospectively selected. MUC5AC, MUC6, gastrin, and TFF2 IHC were performed. Lesions were subjected to laser capture microdissection followed by DNA extraction. Forty-three DNA samples were extracted from PGA without cytological dysplasia, PGA with low-grade and high-grade dysplasia and pyloric gland adenocarcinoma, GC, SPEM, and adjacent normal tissue from the body of the stomach and were subjected to exome sequencing for 49 genes that are commonly dysregulated in GC. Sanger sequencing was performed for confirmation. Twenty nonsynonymous mutations were identified in SPEM, and none of these were frameshifts or indels. PGA with or without cytological dysplasia showed a significantly higher number of mutations compared with SPEM. As cytological dysplasia increased from no dysplasia to dysplasia in PGA, the percentage of frameshift mutations, indels, and missense variations increased. Further missense or frameshift mutations were observed in the KRAS, APC, TP53, and CTNNB1 genes in the PGA group. In GC, mutations were observed in the TP53 gene (p.Arg248Gln). Missense mutations in the MUC5AC, KRAS, BRAF, and EZH2 genes were common between SPEM and GC. SPEM showed fewer genomic variations than GC and PGA, and was genomically distinct from the pyloric epithelium in PGA. Stepwise progression of PGA from PGA without dysplasia to PGA with dysplasia/adenocarcinoma was associated an increase in mutations. SPEM appears to be more genomically similar to GC than PGA.

Activation of Transforming Growth Factor Beta 1 Signaling in Gastric Cancer-associated Fibroblasts Increases Their Motility, via Expression of Rhomboid 5 Homolog 2, and Ability to Induce Invasiveness of Gastric Cancer Cells.

BACKGROUND & AIMS: Fibroblasts that interact with cancer cells are called cancer-associated fibroblasts (CAFs), which promote progression of different tumor types. We investigated the characteristics and functions of CAFs in diffuse-type gastric cancers (DGCs) by analyzing features of their genome and gene expression patterns. METHODS: We isolated CAFs and adjacent non-cancer fibroblasts (NFs) from 110 gastric cancer (GC) tissues from patients who underwent gastrectomy in Japan from 2008 through 2016. Cells were identified using specific markers of various cell types by immunoblot and flow cytometry. We selected pairs of CAFs and NFs for whole-exome and RNA sequencing analyses, and compared expression of specific genes using quantitative reverse transcription PCR. Protein levels and phosphorylation were compared by immunoblot and immunofluorescence analyses. Rhomboid 5 homolog 2 (RHBDF2) was overexpressed from a transgene in fibroblasts or knocked down using small interfering RNAs. Motility and invasiveness of isolated fibroblasts and GC cell lines (AGS, KATOIII, MKN45, NUGC3, NUGC4, OCUM-2MD3 and OCUM-12 cell lines) were quantified by real-time imaging analyses. We analyzed 7 independent sets of DNA microarray data from patients with GC and associated expression levels of specific genes with patient survival times. Nude mice were given injections of OCUM-2MD3 in the stomach wall; tumors and metastases were collected and analyzed by immunohistochemistry. RESULTS: Many of the genes with increased expression in CAFs compared with NFs were associated with transforming growth factor beta 1 (TGFB1) activity. When CAFs were cultured in extracellular matrix, they became more motile than NFs; DGC cells incubated with CAFs were also more motile and invasive in vitro than DGC cells not incubated with CAFs. When injected into nude mice, CAF-incubated DGC cells invaded a greater number of lymphatic vessels than NF-incubated DGC cells. We identified RHBDF2 as a gene overexpressed in CAFs compared with NFs. Knockdown of RHBDF2 in CAFs reduced their elongation and motility in response to TGFB1, whereas overexpression of RHBDF2 in NFs increased their motility in extracellular matrix. RHBDF2 appeared to regulate oncogenic and non-canonical TGFB1 signaling. Knockdown of RHBDF2 in CAFs reduced cleavage of the TGFB receptor 1 (TGFBR1) by ADAM metallopeptidase domain 17 (ADAM17 or TACE) and reduced expression of genes that regulate motility. Incubation of NFs with in interleukin 1 alpha (IL1A), IL1B or tumor necrosis factor, secreted by DGCs, increased fibroblast expression of RHBDF2. Simultaneous high expression of these cytokines in GC samples was associated with shorter survival times of patients. CONCLUSIONS: In CAFs isolated from human DGCs, we observed increased expression of RHBDF2, which regulates TGFB1 signaling. Expression of RHBDF2 in fibroblasts is induced by inflammatory cytokines (such as IL1A, IL1B, and tumor necrosis factor) secreted by DGCs. RHBDF2 promotes cleavage of TGFBR1 by activating TACE and motility of CAFs in response to TGFB1. These highly motile CAFs induce DGCs to invade extracellular matrix and lymphatic vessels in nude mice.

SETD2 histone modifier loss in aggressive GI stromal tumours.

BACKGROUND: GI stromal tumours (GISTs) are clinically heterogenous exhibiting varying degrees of disease aggressiveness in individual patients. OBJECTIVES: We sought to identify genetic alterations associated with high-risk GIST, explore their molecular consequences, and test their utility as prognostic markers. DESIGNS: Exome sequencing of 18 GISTs was performed (9 patients with high-risk/metastatic and 5 patients with low/intermediate-risk), corresponding to 11 primary and 7 metastatic tumours. Candidate alterations were validated by prevalence screening in an independent patient cohort (n=120). Functional consequences of SETD2 mutations were investigated in primary tissues and cell lines. Transcriptomic profiles for 8 GISTs (4 SETD2 mutated, 4 SETD2 wild type) and DNA methylation profiles for 22 GISTs (10 SETD2 mutated, 12 SETD2 wild type) were analysed. Statistical associations between molecular, clinicopathological factors, and relapse-free survival were determined. RESULTS: High-risk GISTs harboured increased numbers of somatic mutations compared with low-risk GISTs (25.2 mutations/high-risk cases vs 6.8 mutations/low-risk cases; two sample t test p=3.1×10-5). Somatic alterations in the SETD2 histone modifier gene occurred in 3 out of 9 high-risk/metastatic cases but no low/intermediate-risk cases. Prevalence screening identified additional SETD2 mutations in 7 out of 80 high-risk/metastatic cases but no low/intermediate-risk cases (n=29). Combined, the frequency of SETD2 mutations was 11.2% (10/89) and 0% (0/34) in high-risk and low-risk GISTs respectively. SETD2 mutant GISTs exhibited decreased H3K36me3 expression while SETD2 silencing promoted DNA damage in GIST-T1 cells. In gastric GISTs, SETD2 mutations were associated with overexpression of HOXC cluster genes and a DNA methylation signature of hypomethylated heterochromatin. Gastric GISTs with SETD2 mutations, or GISTs with hypomethylated heterochromatin, showed significantly shorter relapse-free survival on univariate analysis (log rank p=4.1×10-5). CONCLUSIONS: Our data suggest that SETD2 is a novel GIST tumour suppressor gene associated with disease progression. Assessing SETD2 genetic status and SETD2-associated epigenomic phenotypes may guide risk stratification and provide insights into mechanisms of GIST clinical aggressiveness.

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.

Integrated environmental, lifestyle, and epigenetic risk prediction of primary gastric neoplasia using the longitudinally monitored cohorts.

BACKGROUND: DNA methylation accumulates in non-malignant gastric mucosa after exposure to pathogens. To elucidate how environmental, methylation, and lifestyle factors interplay to influence primary gastric neoplasia (GN) risk, we analyzed longitudinally monitored cohorts in Japan and Singapore. METHODS: Asymptomatic subjects who underwent a gastric mucosal biopsy on the health check-up were enrolled. We analyzed the association between clinical factors and GN development using Cox hazard models. We further conducted comprehensive methylation analysis on selected tissues, including (i) mucosae from subjects developing GN later, (ii) mucosae from subjects not developing GN later, and (iii) GN tissues and surrounding mucosae. We also use the methylation data of mucosa collected in Singapore. The association between methylation and GN risk, as well as lifestyle and methylation, were analyzed. FINDINGS: Among 4234 subjects, GN was developed in 77 subjects. GN incidence was correlated with age, drinking, smoking, and Helicobacter pylori (HP) status. Accumulation of methylation in biopsied gastric mucosae was predictive of higher future GN risk and shorter duration to GN incidence. Whereas methylation levels were associated with HP positivity, lifestyle, and morphological alterations, DNA methylation remained an independent GN risk factor through multivariable analyses. Pro-carcinogenic epigenetic alterations initiated by HP exposure were amplified by unfavorable but modifiable lifestyle choices. Adding DNA methylation to the model with clinical factors improved the predictive ability for the GN risk. INTERPRETATION: The integration of environmental, lifestyle, and epigenetic information can provide increased resolution in the stratification of primary GN risk. FUNDING: The funds are listed in Acknowledgements section.

Spatiotemporal genomic profiling of intestinal metaplasia reveals clonal dynamics of gastric cancer progression.

Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. Analyzing 1,256 gastric samples (1,152 IMs) across 692 subjects from a prospective 10-year study, we identify 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9). Single-cell and spatial profiles highlight changes in tissue ecology and IM lineage heterogeneity, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling reveals expression-based molecular subtypes of IM associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, 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 transform to GC. By highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression, our results raise opportunities for GC precision prevention and interception.

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.

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.

Long-read transcriptome sequencing reveals abundant promoter diversity in distinct molecular subtypes of gastric cancer.

BACKGROUND: Deregulated gene expression is a hallmark of cancer; however, most studies to date have analyzed short-read RNA sequencing data with inherent limitations. Here, we combine PacBio long-read isoform sequencing (Iso-Seq) and Illumina paired-end short-read RNA sequencing to comprehensively survey the transcriptome of gastric cancer (GC), a leading cause of global cancer mortality. RESULTS: We performed full-length transcriptome analysis across 10 GC cell lines covering four major GC molecular subtypes (chromosomal unstable, Epstein-Barr positive, genome stable and microsatellite unstable). We identify 60,239 non-redundant full-length transcripts, of which > 66% are novel compared to current transcriptome databases. Novel isoforms are more likely to be cell line and subtype specific, expressed at lower levels with larger number of exons, with longer isoform/coding sequence lengths. Most novel isoforms utilize an alternate first exon, and compared to other alternative splicing categories, are expressed at higher levels and exhibit higher variability. Collectively, we observe alternate promoter usage in 25% of detected genes, with the majority (84.2%) of known/novel promoter pairs exhibiting potential changes in their coding sequences. Mapping these alternate promoters to TCGA GC samples, we identify several cancer-associated isoforms, including novel variants of oncogenes. Tumor-specific transcript isoforms tend to alter protein coding sequences to a larger extent than other isoforms. Analysis of outcome data suggests that novel isoforms may impart additional prognostic information. CONCLUSIONS: Our results provide a rich resource of full-length transcriptome data for deeper studies of GC and other gastrointestinal malignancies.

CRISPRi enables isoform-specific loss-of-function screens and identification of gastric cancer-specific isoform dependencies.

INTRODUCTION: Genes contain multiple promoters that can drive the expression of various transcript isoforms. Although transcript isoforms from the same gene could have diverse and non-overlapping functions, current loss-of-function methodologies are not able to differentiate between isoform-specific phenotypes. RESULTS: Here, we show that CRISPR interference (CRISPRi) can be adopted for targeting specific promoters within a gene, enabling isoform-specific loss-of-function genetic screens. We use this strategy to test functional dependencies of 820 transcript isoforms that are gained in gastric cancer (GC). We identify a subset of GC-gained transcript isoform dependencies, and of these, we validate CIT kinase as a novel GC dependency. We further show that some genes express isoforms with opposite functions. Specifically, we find that the tumour suppressor ZFHX3 expresses an isoform that has a paradoxical oncogenic role that correlates with poor patient outcome. CONCLUSIONS: Our work finds isoform-specific phenotypes that would not be identified using current loss-of-function approaches that are not designed to target specific transcript isoforms.

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.

Cross-species chromatin interactions drive transcriptional rewiring in Epstein-Barr virus-positive gastric adenocarcinoma.

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.

Genomic and epigenomic EBF1 alterations modulate TERT expression in gastric cancer.

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.

Acquired Resistance to FGFR Inhibitor in Diffuse-Type Gastric Cancer through an AKT-Independent PKC-Mediated Phosphorylation of GSK3ß.

Preclinical models of diffuse-type gastric cancer (DGC) that reliably predict clinical activity of novel compounds are lacking. To overcome the problem of poor tumor cellularity in DGC, we used cells from malignant ascites to establish DGC patient-derived xenograft (PDX) models that recapitulate the primary cancer. Cells in PDX model GAGA6 with FGFR2 amplification were sensitive to AZD4547, a potent FGFR inhibitor that is being clinically evaluated for FGFR-aberrant cancer types. Intermittent in vivo treatment of GAGA6 tumors with AZD4547 gave rise to PDX tumors with acquired resistance to AZD4547, GAGA6-R. Surprisingly, there were no mutations in the FGFR2 gene in GAGA6-R, negating gatekeeper mutations as a mechanism of drug resistance. Phosphorylation of FGFR2 and downstream signaling molecules AKT/PKB and MAPK/ERK remained inhibited by AZD4547. Further analysis of signaling pathways identified AKT-independent phosphorylation and inhibition of GSK3ß as a mechanism of drug resistance in GAGA6-R cells. Treatment of GAGA6-R cells with protein kinase C (PKC) inhibitor H7 in combination with AZD4547 led to dephosphorylation and activation of GSK3ß with concomitant downregulation of MCL-1 and BCL-XL. Combined treatment with AZD4547 and H7 in vitro synergistically enhanced cell death in GAGA6-R but not GAGA6 cells. Furthermore, midostaurin, a multikinase inhibitor with PKC-inhibiting activity, in part reversed resistance of GAGA6-R tumor to AZD4547 in vivo Our results suggest that upon challenge with FGFR inhibitors, FGFR2-amplified tumors that are highly dependent on FGFR2 signaling for survival rapidly develop resistance by switching to a PKC-mediated inhibition of GSK3ß to gain a survival advantage. Mol Cancer Ther; 17(1); 232-42. ©2017 AACR.

The Glaciozyma antarctica genome reveals an array of systems that provide sustained responses towards temperature variations in a persistently cold habitat.

Extremely low temperatures present various challenges to life that include ice formation and effects on metabolic capacity. Psyhcrophilic microorganisms typically have an array of mechanisms to enable survival in cold temperatures. In this study, we sequenced and analysed the genome of a psychrophilic yeast isolated in the Antarctic region, Glaciozyma antarctica. The genome annotation identified 7857 protein coding sequences. From the genome sequence analysis we were able to identify genes that encoded for proteins known to be associated with cold survival, in addition to annotating genes that are unique to G. antarctica. For genes that are known to be involved in cold adaptation such as anti-freeze proteins (AFPs), our gene expression analysis revealed that they were differentially transcribed over time and in response to different temperatures. This indicated the presence of an array of adaptation systems that can respond to a changing but persistent cold environment. We were also able to validate the activity of all the AFPs annotated where the recombinant AFPs demonstrated anti-freeze capacity. This work is an important foundation for further collective exploration into psychrophilic microbiology where among other potential, the genes unique to this species may represent a pool of novel mechanisms for cold survival.