Disentangling oncogenic amplicons in esophageal adenocarcinoma.

Esophageal adenocarcinoma is a prominent example of cancer characterized by frequent amplifications in oncogenes. However, the mechanisms leading to amplicons that involve breakage-fusion-bridge cycles and extrachromosomal DNA are poorly understood. Here, we use 710 esophageal adenocarcinoma cases with matched samples and patient-derived organoids to disentangle complex amplicons and their associated mechanisms. Short-read sequencing identifies ERBB2, MYC, MDM2, and HMGA2 as the most frequent oncogenes amplified in extrachromosomal DNAs. We resolve complex extrachromosomal DNA and breakage-fusion-bridge cycles amplicons by integrating of de-novo assemblies and DNA methylation in nine long-read sequenced cases. Complex amplicons shared between precancerous biopsy and late-stage tumor, an enrichment of putative enhancer elements and mobile element insertions are potential drivers of complex amplicons' origin. We find that patient-derived organoids recapitulate extrachromosomal DNA observed in the primary tumors and single-cell DNA sequencing capture extrachromosomal DNA-driven clonal dynamics across passages. Prospectively, long-read and single-cell DNA sequencing technologies can lead to better prediction of clonal evolution in esophageal adenocarcinoma.

Understanding the malignant potential of gastric metaplasia of the oesophagus and its relevance to Barrett's oesophagus surveillance: individual-level data analysis.

OBJECTIVE: Whether gastric metaplasia (GM) of the oesophagus should be considered as Barrett's oesophagus (BO) is controversial. Given concern intestinal metaplasia (IM) may be missed due to sampling, the UK guidelines include GM as a type of BO. Here, we investigated whether the risk of misdiagnosis and the malignant potential of GM warrant its place in the UK surveillance. DESIGN: We performed a thorough pathology and endoscopy review to follow clinical outcomes in a novel UK cohort of 244 patients, covering 1854 person years of follow-up. We complemented this with a comparative genomic analysis of 160 GM and IM specimens, focused on early molecular hallmarks of BO and oesophageal adenocarcinoma (OAC). RESULTS: We found that 58 of 77 short-segment (<3 cm) GM (SS-GM) cases (75%) continued to be observed as GM-only across a median of 4.4 years of follow-up. We observed that disease progression in GM-only cases and GM+IM cases (cases with reported GM on some occasions, IM on others) was significantly lower than in the IM-only cases (Kaplan-Meier, p=0.03). Genomic analysis revealed that the mutation burden in GM is significantly lower than in IM (p<0.01). Moreover, GM does not bear the mutational hallmarks of OAC, with an absence of associated signatures and driver gene mutations. Finally, we established that GM found adjacent to OAC is evolutionarily distant from cancer. CONCLUSION: SS-GM is a distinct entity from SS-IM and the malignant potential of GM is lower than IM. It is questionable whether SS-GM warrants inclusion in BO surveillance.

SMAD4 and KCNQ3 alterations are associated with lymph node metastases in oesophageal adenocarcinoma.

Metastasis in oesophageal adenocarcinoma (OAC) is an important predictor of survival. Radiological staging is used to stage metastases in patients, and guide treatment selection, but is limited by the accuracy of the approach. Improvements in staging will lead to improved clinical decision making and patient outcomes. Sequencing studies on primary tumours and pre-cancerous tissue have revealed the mutational landscape of OAC, and increasingly cheap and widespread sequencing approaches offer the potential to improve staging assessment. In this work we present an analysis of lymph node metastases found by radiological and pathological sampling, identifying new roles of the genes SMAD4 and KCNQ3 in metastasis. Through transcriptomic analysis we find that both genes are associated with canonical Wnt pathway activity, but KCNQ3 is uniquely associated with changes in planar cell polaritiy associated with non-canonical Wnt signalling. We go on to validate our observations in KCNQ3 in cell line and xenograph systems, showing that overexpression of KCNQ3 reduces wound closure and the number of metastases observed. Our results suggest both genes as novel biomarkers of metastatic risk and offer new potential routes to drug targeting.

KCNQ potassium channels modulate Wnt activity in gastro-oesophageal adenocarcinomas.

Voltage-sensitive potassium channels play an important role in controlling membrane potential and ionic homeostasis in the gut and have been implicated in gastrointestinal (GI) cancers. Through large-scale analysis of 897 patients with gastro-oesophageal adenocarcinomas (GOAs) coupled with in vitro models, we find KCNQ family genes are mutated in ∼30% of patients, and play therapeutically targetable roles in GOA cancer growth. KCNQ1 and KCNQ3 mediate the WNT pathway and MYC to increase proliferation through resultant effects on cadherin junctions. This also highlights novel roles of KCNQ3 in non-excitable tissues. We also discover that activity of KCNQ3 sensitises cancer cells to existing potassium channel inhibitors and that inhibition of KCNQ activity reduces proliferation of GOA cancer cells. These findings reveal a novel and exploitable role of potassium channels in the advancement of human cancer, and highlight that supplemental treatments for GOAs may exist through KCNQ inhibitors.

The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes.

Most cancer types exhibit aberrant transcriptional activity, including derepression of retrotransposable elements (RTEs). However, the degree, specificity and potential consequences of RTE transcriptional activation may differ substantially among cancer types and subtypes. Representing one extreme of the spectrum, we characterize the transcriptional activity of RTEs in cohorts of esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) from the OCCAMS (Oesophageal Cancer Clinical and Molecular Stratification) consortium, and from TCGA (The Cancer Genome Atlas). We found exceptionally high RTE inclusion in the EAC transcriptome, driven primarily by transcription of genes incorporating intronic or adjacent RTEs, rather than by autonomous RTE transcription. Nevertheless, numerous chimeric transcripts straddling RTEs and genes, and transcripts from stand-alone RTEs, particularly KLF5- and SOX9-controlled HERVH proviruses, were overexpressed specifically in EAC. Notably, incomplete mRNA splicing and EAC-characteristic intronic RTE inclusion was mirrored by relative loss of the respective fully-spliced, functional mRNA isoforms, consistent with compromised cellular fitness. Defective RNA splicing was linked with strong transcriptional activation of a HERVH provirus on Chr Xp22.32 and defined EAC subtypes with distinct molecular features and prognosis. Our study defines distinguishable RTE transcriptional profiles of EAC, reflecting distinct underlying processes and prognosis, thus providing a framework for targeted studies.

Mutational signature dynamics shaping the evolution of oesophageal adenocarcinoma.

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma.

Single-Cell RNA Sequencing Unifies Developmental Programs of Esophageal and Gastric Intestinal Metaplasia.

Intestinal metaplasia in the esophagus (Barrett's esophagus IM, or BE-IM) and stomach (GIM) are considered precursors for esophageal and gastric adenocarcinoma, respectively. We hypothesize that BE-IM and GIM follow parallel developmental trajectories in response to differing inflammatory insults. Here, we construct a single-cell RNA-sequencing atlas, supported by protein expression studies, of the entire gastrointestinal tract spanning physiologically normal and pathologic states including gastric metaplasia in the esophagus (E-GM), BE-IM, atrophic gastritis, and GIM. We demonstrate that BE-IM and GIM share molecular features, and individual cells simultaneously possess transcriptional properties of gastric and intestinal epithelia, suggesting phenotypic mosaicism. Transcriptionally E-GM resembles atrophic gastritis; genetically, it is clonal and has a lower mutational burden than BE-IM. Finally, we show that GIM and BE-IM acquire a protumorigenic, activated fibroblast microenvironment. These findings suggest that BE-IM and GIM can be considered molecularly similar entities in adjacent organs, opening the path for shared detection and treatment strategies. SIGNIFICANCE: Our data capture the gradual molecular and phenotypic transition from a gastric to intestinal phenotype (IM) in the esophagus and stomach. Because BE-IM and GIM can predispose to cancer, this new understanding of a common developmental trajectory could pave the way for a more unified approach to detection and treatment. See related commentary by Stachler, p. 1291. This article is highlighted in the In This Issue feature, p. 1275.

Rearrangement processes and structural variations show evidence of selection in oesophageal adenocarcinomas.

Oesophageal adenocarcinoma (OAC) provides an ideal case study to characterize large-scale rearrangements. Using whole genome short-read sequencing of 383 cases, for which 214 had matched whole transcriptomes, we observed structural variations (SV) with a predominance of deletions, tandem duplications and inter-chromosome junctions that could be identified as LINE-1 mobile element (ME) insertions. Complex clusters of rearrangements resembling breakage-fusion-bridge cycles or extrachromosomal circular DNA accounted for 22% of complex SVs affecting known oncogenes. Counting SV events affecting known driver genes substantially increased the recurrence rates of these drivers. After excluding fragile sites, we identified 51 candidate new drivers in genomic regions disrupted by SVs, including ETV5, KAT6B and CLTC. RUNX1 was the most recurrently altered gene (24%), with many deletions inactivating the RUNT domain but preserved the reading frame, suggesting an altered protein product. These findings underscore the importance of identification of SV events in OAC with implications for targeted therapies.

Molecular phenotyping reveals the identity of Barrett's esophagus and its malignant transition.

The origin of human metaplastic states and their propensity for cancer is poorly understood. Barrett's esophagus is a common metaplastic condition that increases the risk for esophageal adenocarcinoma, and its cellular origin is enigmatic. To address this, we harvested tissues spanning the gastroesophageal junction from healthy and diseased donors, including isolation of esophageal submucosal glands. A combination of single-cell transcriptomic profiling, in silico lineage tracing from methylation, open chromatin and somatic mutation analyses, and functional studies in organoid models showed that Barrett's esophagus originates from gastric cardia through c-MYC and HNF4A-driven transcriptional programs. Furthermore, our data indicate that esophageal adenocarcinoma likely arises from undifferentiated Barrett's esophagus cell types even in the absence of a pathologically identifiable metaplastic precursor, illuminating early detection strategies.

Genomic Analysis of Response to Neoadjuvant Chemotherapy in Esophageal Adenocarcinoma.

Neoadjuvant therapy followed by surgery is the standard of care for locally advanced esophageal adenocarcinoma (EAC). Unfortunately, response to neoadjuvant chemotherapy (NAC) is poor (20-37%), as is the overall survival benefit at five years (9%). The EAC genome is complex and heterogeneous between patients, and it is not yet understood whether specific mutational patterns may result in chemotherapy sensitivity or resistance. To identify associations between genomic events and response to NAC in EAC, a comparative genomic analysis was performed in 65 patients with extensive clinical and pathological annotation using whole-genome sequencing (WGS). We defined response using Mandard Tumor Regression Grade (TRG), with responders classified as TRG1-2 (n = 27) and non-responders classified as TRG4-5 (n =38). We report a higher non-synonymous mutation burden in responders (median 2.08/Mb vs. 1.70/Mb, p = 0.036) and elevated copy number variation in non-responders (282 vs. 136/patient, p < 0.001). We identified copy number variants unique to each group in our cohort, with cell cycle (CDKN2A, CCND1), c-Myc (MYC), RTK/PIK3 (KRAS, EGFR) and gastrointestinal differentiation (GATA6) pathway genes being specifically altered in non-responders. Of note, NAV3 mutations were exclusively present in the non-responder group with a frequency of 22%. Thus, lower mutation burden, higher chromosomal instability and specific copy number alterations are associated with resistance to NAC.

Evidence that polyploidy in esophageal adenocarcinoma originates from mitotic slippage caused by defective chromosome attachments.

Polyploidy is present in many cancer types and is increasingly recognized as an important factor in promoting chromosomal instability, genome evolution, and heterogeneity in cancer cells. However, the mechanisms that trigger polyploidy in cancer cells are largely unknown. In this study, we investigated the origin of polyploidy in esophageal adenocarcinoma (EAC), a highly heterogenous cancer, using a combination of genomics and cell biology approaches in EAC cell lines, organoids, and tumors. We found the EAC cells and organoids present specific mitotic defects consistent with problems in the attachment of chromosomes to the microtubules of the mitotic spindle. Time-lapse analyses confirmed that EAC cells have problems in congressing and aligning their chromosomes, which can ultimately culminate in mitotic slippage and polyploidy. Furthermore, whole-genome sequencing, RNA-seq, and quantitative immunofluorescence analyses revealed alterations in the copy number, expression, and cellular distribution of several proteins known to be involved in the mechanics and regulation of chromosome dynamics during mitosis. Together, these results provide evidence that an imbalance in the amount of proteins implicated in the attachment of chromosomes to spindle microtubules is the molecular mechanism underlying mitotic slippage in EAC. Our findings that the likely origin of polyploidy in EAC is mitotic failure caused by problems in chromosomal attachments not only improves our understanding of cancer evolution and diversification, but may also aid in the classification and treatment of EAC and possibly other highly heterogeneous cancers.

Alterations in the Ca2+ toolkit in oesophageal adenocarcinoma.

Aim: To investigate alterations in transcription of genes, encoding Ca2+ toolkit proteins, in oesophageal adenocarcinoma (OAC) and to assess associations between gene expression, tumor grade, nodal-metastatic stage, and patient survival. Methods: The expression of 275 transcripts, encoding components of the Ca2+ toolkit, was analyzed in two OAC datasets: the Cancer Genome Atlas [via the University of Alabama Cancer (UALCAN) portal] and the oesophageal-cancer, clinical, and molecular stratification [Oesophageal Cancer Clinical and Molecular Stratification (OCCAMS)] dataset. Effects of differential expression of these genes on patient survival were determined using Kaplan-Meier log-rank tests. OAC grade- and metastatic-stage status was investigated for a subset of genes. Adjustment for the multiplicity of testing was made throughout. Results: Of the 275 Ca2+-toolkit genes analyzed, 75 displayed consistent changes in expression between OAC and normal tissue in both datasets. The channel-encoding genes, N-methyl-D-aspartate receptor 2D (GRIN2D), transient receptor potential (TRP) ion channel classical or canonical 4 (TRPC4), and TRP ion channel melastatin 2 (TRPM2) demonstrated the greatest increase in expression in OAC in both datasets. Nine genes were consistently upregulated in both datasets and were also associated with improved survival outcomes. The 6 top-ranking genes for the weighted significance of altered expression and survival outcomes were selected for further analysis: voltage-gated Ca2+ channel subunit α 1D (CACNA1D), voltage-gated Ca2+ channel auxiliary subunit α2 δ4 (CACNA2D4), junctophilin 1 (JPH1), acid-sensing ion channel 4 (ACCN4), TRPM5, and secretory pathway Ca2+ ATPase 2 (ATP2C2). CACNA1D, JPH1, and ATP2C2 were also upregulated in advanced OAC tumor grades and nodal-metastatic stages in both datasets. Conclusions: This study has unveiled alterations of the Ca2+ toolkit in OAC, compared to normal tissue. Such Ca2+ signalling findings are consistent with those from studies on other cancers. Genes that were consistently upregulated in both datasets might represent useful markers for patient diagnosis. Genes that were consistently upregulated, and which were associated with improved survival, might be useful markers for patient outcome. These survival-associated genes may also represent targets for the development of novel chemotherapeutic agents.

The mutREAD method detects mutational signatures from low quantities of cancer DNA.

Mutational processes acting on cancer genomes can be traced by investigating mutational signatures. Because high sequencing costs limit current studies to small numbers of good-quality samples, we propose a robust, cost- and time-effective method, called mutREAD, to detect mutational signatures from small quantities of DNA, including degraded samples. We show that mutREAD recapitulates mutational signatures identified by whole genome sequencing, and will ultimately allow the study of mutational signatures in larger cohorts and, by compatibility with formalin-fixed paraffin-embedded samples, in clinical settings.

Identification of Subtypes of Barrett's Esophagus and Esophageal Adenocarcinoma Based on DNA Methylation Profiles and Integration of Transcriptome and Genome Data.

BACKGROUND & AIMS: Esophageal adenocarcinomas (EACs) are heterogeneous and often preceded by Barrett's esophagus (BE). Many genomic changes have been associated with development of BE and EAC, but little is known about epigenetic alterations. We performed epigenetic analyses of BE and EAC tissues and combined these data with transcriptome and genomic data to identify mechanisms that control gene expression and genome integrity. METHODS: In a retrospective cohort study, we collected tissue samples and clinical data from 150 BE and 285 EAC cases from the Oesophageal Cancer Classification and Molecular Stratification consortium in the United Kingdom. We analyzed methylation profiles of all BE and EAC tissues and assigned them to subgroups using non-negative matrix factorization with k-means clustering. Data from whole-genome sequencing and transcriptome studies were then incorporated; we performed integrative methylation and RNA-sequencing analyses to identify genes that were suppressed with increased methylation in promoter regions. Levels of different immune cell types were computed using single-sample gene set enrichment methods. We derived 8 organoids from 8 EAC tissues and tested their sensitivity to different drugs. RESULTS: BE and EAC samples shared genome-wide methylation features, compared with normal tissues (esophageal, gastric, and duodenum; controls) from the same patients and grouped into 4 subtypes. Subtype 1 was characterized by DNA hypermethylation with a high mutation burden and multiple mutations in genes in cell cycle and receptor tyrosine signaling pathways. Subtype 2 was characterized by a gene expression pattern associated with metabolic processes (ATP synthesis and fatty acid oxidation) and lack methylation at specific binding sites for transcription factors; 83% of samples of this subtype were BE and 17% were EAC. The third subtype did not have changes in methylation pattern, compared with control tissue, but had a gene expression pattern that indicated immune cell infiltration; this tumor type was associated with the shortest time of patient survival. The fourth subtype was characterized by DNA hypomethylation associated with structure rearrangements, copy number alterations, with preferential amplification of CCNE1 (cells with this gene amplification have been reported to be sensitive to CDK2 inhibitors). Organoids with reduced levels of MGMT and CHFR expression were sensitive to temozolomide and taxane drugs. CONCLUSIONS: In a comprehensive integrated analysis of methylation, transcriptome, and genome profiles of more than 400 BE and EAC tissues, along with clinical data, we identified 4 subtypes that were associated with patient outcomes and potential responses to therapy.

Transcriptomic profiling reveals three molecular phenotypes of adenocarcinoma at the gastroesophageal junction.

Cancers occurring at the gastroesophageal junction (GEJ) are classified as predominantly esophageal or gastric, which is often difficult to decipher. We hypothesized that the transcriptomic profile might reveal molecular subgroups which could help to define the tumor origin and behavior beyond anatomical location. The gene expression profiles of 107 treatment-naïve, intestinal type, gastroesophageal adenocarcinomas were assessed by the Illumina-HTv4.0 beadchip. Differential gene expression (limma), unsupervised subgroup assignment (mclust) and pathway analysis (gage) were undertaken in R statistical computing and results were related to demographic and clinical parameters. Unsupervised assignment of the gene expression profiles revealed three distinct molecular subgroups, which were not associated with anatomical location, tumor stage or grade (p > 0.05). Group 1 was enriched for pathways involved in cell turnover, Group 2 was enriched for metabolic processes and Group 3 for immune-response pathways. Patients in group 1 showed the worst overall survival (p = 0.019). Key genes for the three subtypes were confirmed by immunohistochemistry. The newly defined intrinsic subtypes were analyzed in four independent datasets of gastric and esophageal adenocarcinomas with transcriptomic data available (RNAseq data: OCCAMS cohort, n = 158; gene expression arrays: Belfast, n = 63; Singapore, n = 191; Asian Cancer Research Group, n = 300). The subgroups were represented in the independent cohorts and pooled analysis confirmed the prognostic effect of the new subtypes. In conclusion, adenocarcinomas at the GEJ comprise three distinct molecular phenotypes which do not reflect anatomical location but rather inform our understanding of the key pathways expressed.

The landscape of selection in 551 esophageal adenocarcinomas defines genomic biomarkers for the clinic.

Esophageal adenocarcinoma (EAC) is a poor-prognosis cancer type with rapidly rising incidence. Understanding of the genetic events driving EAC development is limited, and there are few molecular biomarkers for prognostication or therapeutics. Using a cohort of 551 genomically characterized EACs with matched RNA sequencing data, we discovered 77 EAC driver genes and 21 noncoding driver elements. We identified a mean of 4.4 driver events per tumor, which were derived more commonly from mutations than copy number alterations, and compared the prevelence of these mutations to the exome-wide mutational excess calculated using non-synonymous to synonymous mutation ratios (dN/dS). We observed mutual exclusivity or co-occurrence of events within and between several dysregulated EAC pathways, a result suggestive of strong functional relationships. Indicators of poor prognosis (SMAD4 and GATA4) were verified in independent cohorts with significant predictive value. Over 50% of EACs contained sensitizing events for CDK4 and CDK6 inhibitors, which were highly correlated with clinically relevant sensitivity in a panel of EAC cell lines and organoids.

Organoid cultures recapitulate esophageal adenocarcinoma heterogeneity providing a model for clonality studies and precision therapeutics.

Esophageal adenocarcinoma (EAC) incidence is increasing while 5-year survival rates remain less than 15%. A lack of experimental models has hampered progress. We have generated clinically annotated EAC organoid cultures that recapitulate the morphology, genomic, and transcriptomic landscape of the primary tumor including point mutations, copy number alterations, and mutational signatures. Karyotyping of organoid cultures has confirmed polyclonality reflecting the clonal architecture of the primary tumor. Furthermore, subclones underwent clonal selection associated with driver gene status. Medium throughput drug sensitivity testing demonstrates the potential of targeting receptor tyrosine kinases and downstream mediators. EAC organoid cultures provide a pre-clinical tool for studies of clonal evolution and precision therapeutics.