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.

Low-cost and clinically applicable copy number profiling using repeat DNA.

BACKGROUND: Somatic copy number alterations (SCNAs) are an important class of genomic alteration in cancer. They are frequently observed in cancer samples, with studies showing that, on average, SCNAs affect 34% of a cancer cell's genome. Furthermore, SCNAs have been shown to be major drivers of tumour development and have been associated with response to therapy and prognosis. Large-scale cancer genome studies suggest that tumours are driven by somatic copy number alterations (SCNAs) or single-nucleotide variants (SNVs). Despite the frequency of SCNAs and their clinical relevance, the use of genomics assays in the clinic is biased towards targeted gene panels, which identify SNVs but provide limited scope to detect SCNAs throughout the genome. There is a need for a comparably low-cost and simple method for high-resolution SCNA profiling. RESULTS: We present conliga, a fully probabilistic method that infers SCNA profiles from a low-cost, simple, and clinically-relevant assay (FAST-SeqS). When applied to 11 high-purity oesophageal adenocarcinoma samples, we obtain good agreement (Spearman's rank correlation coefficient, rs=0.94) between conliga's inferred SCNA profiles using FAST-SeqS data (approximately £14 per sample) and those inferred by ASCAT using high-coverage WGS (gold-standard). We find that conliga outperforms CNVkit (rs=0.89), also applied to FAST-SeqS data, and is comparable to QDNAseq (rs=0.96) applied to low-coverage WGS, which is approximately four-fold more expensive, more laborious and less clinically-relevant. By performing an in silico dilution series experiment, we find that conliga is particularly suited to detecting SCNAs in low tumour purity samples. At two million reads per sample, conliga is able to detect SCNAs in all nine samples at 3% tumour purity and as low as 0.5% purity in one sample. Crucially, we show that conliga's hidden state information can be used to decide when a sample is abnormal or normal, whereas CNVkit and QDNAseq cannot provide this critical information. CONCLUSIONS: We show that conliga provides high-resolution SCNA profiles using a convenient, low-cost assay. We believe conliga makes FAST-SeqS a more clinically valuable assay as well as a useful research tool, enabling inexpensive and fast copy number profiling of pre-malignant and cancer samples.

Multiparametric High-Content Cell Painting Identifies Copper Ionophores as Selective Modulators of Esophageal Cancer Phenotypes.

Esophageal adenocarcinoma is of increasing global concern due to increasing incidence, a lack of effective treatments, and poor prognosis. Therapeutic target discovery and clinical trials have been hindered by the heterogeneity of the disease, the lack of "druggable" driver mutations, and the dominance of large-scale genomic rearrangements. We have previously undertaken a comprehensive small-molecule phenotypic screen using the high-content Cell Painting assay to quantify the morphological response to a total of 19,555 small molecules across a panel of genetically distinct human esophageal cell lines to identify new therapeutic targets and small molecules for the treatment of esophageal adenocarcinoma. In this current study, we report for the first time the dose-response validation studies for the 72 screening hits from the target-annotated LOPAC and Prestwick FDA-approved compound libraries and the full list of 51 validated esophageal adenocarcinoma-selective small molecules (71% validation rate). We then focus on the most potent and selective hit molecules, elesclomol, disulfiram, and ammonium pyrrolidinedithiocarbamate. Using a multipronged, multitechnology approach, we uncover a unified mechanism of action and a vulnerability in esophageal adenocarcinoma toward copper-dependent cell death that could be targeted in the future.

A suspension technique for efficient large-scale cancer organoid culturing and perturbation screens.

Organoid cell culture methodologies are enabling the generation of cell models from healthy and diseased tissue. Patient-derived cancer organoids that recapitulate the genetic and histopathological diversity of patient tumours are being systematically generated, providing an opportunity to investigate new cancer biology and therapeutic approaches. The use of organoid cultures for many applications, including genetic and chemical perturbation screens, is limited due to the technical demands and cost associated with their handling and propagation. Here we report and benchmark a suspension culture technique for cancer organoids which allows for the expansion of models to tens of millions of cells with increased efficiency in comparison to standard organoid culturing protocols. Using whole-genome DNA and RNA sequencing analyses, as well as medium-throughput drug sensitivity testing and genome-wide CRISPR-Cas9 screening, we demonstrate that cancer organoids grown as a suspension culture are genetically and phenotypically similar to their counterparts grown in standard conditions. This culture technique simplifies organoid cell culture and extends the range of organoid applications, including for routine use in large-scale perturbation screens.

Computational pathology aids derivation of microRNA biomarker signals from Cytosponge samples.

BACKGROUND: Non-endoscopic cell collection devices combined with biomarkers can detect Barrett's intestinal metaplasia and early oesophageal cancer. However, assays performed on multi-cellular samples lose information about the cell source of the biomarker signal. This cross-sectional study examines whether a bespoke artificial intelligence-based computational pathology tool could ascertain the cellular origin of microRNA biomarkers, to inform interpretation of the disease pathology, and confirm biomarker validity. METHODS: The microRNA expression profiles of 110 targets were assessed with a custom multiplexed panel in a cohort of 117 individuals with reflux that took a Cytosponge test. A computational pathology tool quantified the amount of columnar epithelium present in pathology slides, and results were correlated with microRNA signals. An independent cohort of 139 Cytosponges, each from an individual patient, was used to validate the findings via qPCR. FINDINGS: Seventeen microRNAs are upregulated in BE compared to healthy squamous epithelia, of which 13 remain upregulated in dysplasia. A pathway enrichment analysis confirmed association to neoplastic and cell cycle regulation processes. Ten microRNAs positively correlated with columnar epithelium content, with miRNA-192-5p and -194-5p accurately detecting the presence of gastric cells (AUC 0.97 and 0.95). In contrast, miR-196a-5p is confirmed as a specific BE marker. INTERPRETATION: Computational pathology tools aid accurate cellular attribution of molecular signals. This innovative design with multiplex microRNA coupled with artificial intelligence has led to discovery of a quality control metric suitable for large scale application of the Cytosponge. Similar approaches could aid optimal interpretation of biomarkers for clinical use. FUNDING: Funded by the NIHR Cambridge Biomedical Research Centre, the Medical Research Council, the Rosetrees and Stoneygate Trusts, and CRUK core grants.

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.

A thiol-bound drug reservoir enhances APR-246-induced mutant p53 tumor cell death.

The tumor suppressor gene TP53 is the most frequently mutated gene in cancer. The compound APR-246 (PRIMA-1Met/Eprenetapopt) is converted to methylene quinuclidinone (MQ) that targets mutant p53 protein and perturbs cellular antioxidant balance. APR-246 is currently tested in a phase III clinical trial in myelodysplastic syndrome (MDS). By in vitro, ex vivo, and in vivo models, we show that combined treatment with APR-246 and inhibitors of efflux pump MRP1/ABCC1 results in synergistic tumor cell death, which is more pronounced in TP53 mutant cells. This is associated with altered cellular thiol status and increased intracellular glutathione-conjugated MQ (GS-MQ). Due to the reversibility of MQ conjugation, GS-MQ forms an intracellular drug reservoir that increases availability of MQ for targeting mutant p53. Our study shows that redox homeostasis is a critical determinant of the response to mutant p53-targeted cancer therapy.

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.

Genomic evidence supports a clonal diaspora model for metastases of esophageal adenocarcinoma.

The poor outcomes in esophageal adenocarcinoma (EAC) prompted us to interrogate the pattern and timing of metastatic spread. Whole-genome sequencing and phylogenetic analysis of 388 samples across 18 individuals with EAC showed, in 90% of patients, that multiple subclones from the primary tumor spread very rapidly from the primary site to form multiple metastases, including lymph nodes and distant tissues-a mode of dissemination that we term 'clonal diaspora'. Metastatic subclones at autopsy were present in tissue and blood samples from earlier time points. These findings have implications for our understanding and clinical evaluation of EAC.

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.

Selection and Application of Tissue microRNAs for Nonendoscopic Diagnosis of Barrett's Esophagus.

BACKGROUND & AIMS: MicroRNA (miRNA) is highly stable in biospecimens and provides tissue-specific profiles, making it a useful biomarker of carcinogenesis. We aimed to discover a set of miRNAs that could accurately discriminate Barrett's esophagus (BE) from normal esophageal tissue and to test its diagnostic accuracy when applied to samples collected by a noninvasive esophageal cell sampling device. METHODS: We analyzed miRNA expression profiles of 2 independent sets of esophageal biopsy tissues collected during endoscopy from 38 patients with BE and 26 patients with normal esophagus (controls) using Agilent microarray and Nanostring nCounter assays. Consistently up-regulated miRNAs were quantified by real-time polymerase chain reaction in esophageal tissues collected by Cytosponge from patients with BE vs without BE. miRNAs were expressed from plasmids and antisense oligonucleotides were expressed in normal esophageal squamous cells; effects on proliferation and gene expression patterns were analyzed. RESULTS: We identified 15 miRNAs that were significantly up-regulated in BE vs control tissues. Of these, 11 (MIR215, MIR194, MIR 192, MIR196a, MIR199b, MIR10a, MIR145, MIR181a, MIR30a, MIR7, and MIR199a) were validated in Cytosponge samples. The miRNAs with the greatest increases in BE tissues (7.9-fold increase in expression or more, P < .0001: MIR196a, MIR192, MIR194, and MIR215) each identified BE vs control tissues with area under the curve (AUC) values of 0.82 or more. We developed an optimized multivariable logistic regression model, based on expression levels of 6 miRNAs (MIR7, MIR30a, MIR181a, MIR192, MIR196a, and MIR199a), that identified patients with BE with an AUC value of 0.89, 86.2% sensitivity, and 91.6% specificity. Expression level of MIR192, MIR196a, MIR199a, combined that of trefoil factor 3, identified patients with BE with an AUC of 0.93, 93.1% sensitivity, and 93.7% specificity. Hypomethylation was observed in the promoter region of the highly up-regulated cluster MIR192-MIR194. Overexpression of these miRNAs in normal esophageal squamous cells increased their proliferation, via GRHL3 and PTEN signaling. CONCLUSIONS: In analyses of miRNA expression patterns of BE vs non-BE tissues, we identified a profile that can identify Cytosponge samples from patients with BE with an AUC of 0.93. Expression of MIR194 is increased in BE samples via epigenetic mechanisms that might be involved in BE pathogenesis.

Open chromatin profiling identifies AP1 as a transcriptional regulator in oesophageal adenocarcinoma.

Oesophageal adenocarcinoma (OAC) is one of the ten most prevalent forms of cancer and is showing a rapid increase in incidence and yet exhibits poor survival rates. Compared to many other common cancers, the molecular changes that occur in this disease are relatively poorly understood. However, genes encoding chromatin remodeling enzymes are frequently mutated in OAC. This is consistent with the emerging concept that cancer cells exhibit reprogramming of their chromatin environment which leads to subsequent changes in their transcriptional profile. Here, we have used ATAC-seq to interrogate the chromatin changes that occur in OAC using both cell lines and patient-derived material. We demonstrate that there are substantial changes in the regulatory chromatin environment in the cancer cells and using this data we have uncovered an important role for ETS and AP1 transcription factors in driving the changes in gene expression found in OAC cells.

Impact of mutations in Toll-like receptor pathway genes on esophageal carcinogenesis.

Esophageal adenocarcinoma (EAC) develops in an inflammatory microenvironment with reduced microbial diversity, but mechanisms for these influences remain poorly characterized. We hypothesized that mutations targeting the Toll-like receptor (TLR) pathway could disrupt innate immune signaling and promote a microenvironment that favors tumorigenesis. Through interrogating whole genome sequencing data from 171 EAC patients, we showed that non-synonymous mutations collectively affect the TLR pathway in 25/171 (14.6%, PathScan p = 8.7x10-5) tumors. TLR mutant cases were associated with more proximal tumors and metastatic disease, indicating possible clinical significance of these mutations. Only rare mutations were identified in adjacent Barrett's esophagus samples. We validated our findings in an external EAC dataset with non-synonymous TLR pathway mutations in 33/149 (22.1%, PathScan p = 0.05) tumors, and in other solid tumor types exposed to microbiomes in the COSMIC database (10,318 samples), including uterine endometrioid carcinoma (188/320, 58.8%), cutaneous melanoma (377/988, 38.2%), colorectal adenocarcinoma (402/1519, 26.5%), and stomach adenocarcinoma (151/579, 26.1%). TLR4 was the most frequently mutated gene with eleven mutations in 10/171 (5.8%) of EAC tumors. The TLR4 mutants E439G, S570I, F703C and R787H were confirmed to have impaired reactivity to bacterial lipopolysaccharide with marked reductions in signaling by luciferase reporter assays. Overall, our findings show that TLR pathway genes are recurrently mutated in EAC, and TLR4 mutations have decreased responsiveness to bacterial lipopolysaccharide and may play a role in disease pathogenesis in a subset of patients.

Mutational signatures in esophageal adenocarcinoma define etiologically distinct subgroups with therapeutic relevance.

Esophageal adenocarcinoma (EAC) has a poor outcome, and targeted therapy trials have thus far been disappointing owing to a lack of robust stratification methods. Whole-genome sequencing (WGS) analysis of 129 cases demonstrated that this is a heterogeneous cancer dominated by copy number alterations with frequent large-scale rearrangements. Co-amplification of receptor tyrosine kinases (RTKs) and/or downstream mitogenic activation is almost ubiquitous; thus tailored combination RTK inhibitor (RTKi) therapy might be required, as we demonstrate in vitro. However, mutational signatures showed three distinct molecular subtypes with potential therapeutic relevance, which we verified in an independent cohort (n = 87): (i) enrichment for BRCA signature with prevalent defects in the homologous recombination pathway; (ii) dominant T>G mutational pattern associated with a high mutational load and neoantigen burden; and (iii) C>A/T mutational pattern with evidence of an aging imprint. These subtypes could be ascertained using a clinically applicable sequencing strategy (low coverage) as a basis for therapy selection.

Inhibition of mTOR-kinase destabilizes MYCN and is a potential therapy for MYCN-dependent tumors.

MYC oncoproteins deliver a potent oncogenic stimulus in several human cancers, making them major targets for drug development, but efforts to deliver clinically practical therapeutics have not yet been realized. In childhood cancer, aberrant expression of MYC and MYCN genes delineates a group of aggressive tumours responsible for a major proportion of pediatric cancer deaths. We designed a chemical-genetic screen that identifies compounds capable of enhancing proteasomal elimination of MYCN oncoprotein. We isolated several classes of compound that selectively kill MYCN expressing cells and we focus on inhibitors of PI3K/mTOR pathway in this study. We show that PI3K/mTOR inhibitors selectively killed MYCN-expressing neuroblastoma tumor cells, and induced significant apoptosis of transgenic MYCN-driven neuroblastoma tumors concomitant with elimination of MYCN protein in vivo. Mechanistically, the ability of these compounds to degrade MYCN requires complete blockade of mTOR but not PI3 kinase activity and we highlight NVP-BEZ235 as a PI3K/mTOR inhibitor with an ideal activity profile. These data establish that MYCN expression is a marker indicative of likely clinical sensitivity to mTOR inhibition, and provide a rationale for the selection of clinical candidate MYCN-destabilizers likely to be useful for the treatment of MYCN-driven cancers.