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.

Extrachromosomal DNA in the cancerous transformation of Barrett's oesophagus.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

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.

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.

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.

PRMT1 Is Recruited via DNA-PK to Chromatin Where It Sustains the Senescence-Associated Secretory Phenotype in Response to Cisplatin.

Protein arginine methyltransferase 1 (PRMT1) is overexpressed in various human cancers and linked to poor response to chemotherapy. Various PRMT1 inhibitors are currently under development; yet, we do not fully understand the mechanisms underpinning PRMT1 involvement in tumorigenesis and chemoresistance. Using mass spectrometry-based proteomics, we identified PRMT1 as regulator of arginine methylation in ovarian cancer cells treated with cisplatin. We showed that DNA-dependent protein kinase (DNA-PK) binds to and phosphorylates PRMT1 in response to cisplatin, inducing its chromatin recruitment and redirecting its enzymatic activity toward Arg3 of histone H4 (H4R3). On chromatin, the DNA-PK/PRMT1 axis induces senescence-associated secretory phenotype through H4R3me2a deposition at pro-inflammatory gene promoters. Finally, PRMT1 inhibition reduces the clonogenic growth of cancer cells exposed to low doses of cisplatin, sensitizing them to apoptosis. While unravelling the role of PRMT1 in response to genotoxic agents, our findings indicate the possibility of targeting PRMT1 to overcome chemoresistance in cancer.

FBXW5 Promotes Tumorigenesis and Metastasis in Gastric Cancer via Activation of the FAK-Src Signaling Pathway.

: F-box/WD repeat-containing protein 5 (FBXW5) is a member of the FBXW subclass of F-box proteins. Despite its known function as a component of the Skp1-Cullin-F-box (SCF) ubiquitin ligase complex, the role of FBXW5 in gastric cancer tumorigenesis and metastasis has not been investigated. The present study investigates the role of FBXW5 in tumorigenesis and metastasis, as well as the regulation of key signaling pathways in gastric cancer; using in-vitro FBXW5 knockdown/overexpression cell line and in-vivo models. In-vitro knockdown of FBXW5 results in a decrease in cell proliferation and cell cycle progression, with a concomitant increase in cell apoptosis and caspase-3 activity. Furthermore, knockdown of FBXW5 also leads to a down regulation in cell migration and adhesion, characterized by a reduction in actin polymerization, focal adhesion turnover and traction forces. This study also delineates the mechanistic role of FBXW5 in oncogenic signaling as its inhibition down regulates RhoA-ROCK 1 (Rho-associated protein kinase 1) and focal adhesion kinase (FAK) signaling cascades. Overexpression of FBXW5 promotes in-vivo tumor growth, whereas its inhibition down regulates in-vivo tumor metastasis. When considered together, our study identifies the novel oncogenic role of FBXW5 in gastric cancer and draws further interest regarding its clinical utility as a potential therapeutic target.

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.

hmSEEKER: Identification of hmSILAC Doublets in MaxQuant Output Data.

Heavy methyl Stable Isotope Labeling with Amino acids in Cell culture (hmSILAC) is a metabolic labeling strategy employed in proteomics to increase the confidence of global identification of methylated peptides by MS. However, to this day, the automatic and robust identification of heavy and light peak doublets from MS-raw data of hmSILAC experiments is a challenging task, for which the choice of computational methods is very limited. Here, hmSEEKER, a software designed to work downstream of a MaxQuant analysis for in-depth search of MS peak pairs that correspond to light and heavy methyl-peptide within MaxQuant-generated tables is described with good sensitivity and specificity. The software is written in Perl, and its code and user manual are freely available at Bitbucket (https://bit.ly/2scCT9u).

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.

Fam60a defines a variant Sin3a-Hdac complex in embryonic stem cells required for self-renewal.

Sin3a is the central scaffold protein of the prototypical Hdac1/2 chromatin repressor complex, crucially required during early embryonic development for the growth of pluripotent cells of the inner cell mass. Here, we compare the composition of the Sin3a-Hdac complex between pluripotent embryonic stem (ES) and differentiated cells by establishing a method that couples two independent endogenous immunoprecipitations with quantitative mass spectrometry. We define the precise composition of the Sin3a complex in multiple cell types and identify the Fam60a subunit as a key defining feature of a variant Sin3a complex present in ES cells, which also contains Ogt and Tet1. Fam60a binds on H3K4me3-positive promoters in ES cells, together with Ogt, Tet1 and Sin3a, and is essential to maintain the complex on chromatin. Finally, we show that depletion of Fam60a phenocopies the loss of Sin3a, leading to reduced proliferation, an extended G1-phase and the deregulation of lineage genes. Taken together, Fam60a is an essential core subunit of a variant Sin3a complex in ES cells that is required to promote rapid proliferation and prevent unscheduled differentiation.

EpiMINE, a computational program for mining epigenomic data.

BACKGROUND: In epigenetic research, both the increasing ease of high-throughput sequencing and a greater interest in genome-wide studies have resulted in an exponential flooding of epigenetic-related data in public domain. This creates an opportunity for exploring data outside the limits of any specific query-centred study. Such data have to undergo standard primary analyses that are accessible with multiple well-stabilized programs. Further downstream analyses, such as genome-wide comparative, correlative and quantitative analyses, are critical in deciphering key biological features. However, these analyses are only accessible for computational researchers and completely lack platforms capable of handling, analysing and linking multiple interdisciplinary datasets with efficient analytical methods. RESULTS: Here, we present EpiMINE, a program for mining epigenomic data. It is a user-friendly, stand-alone computational program designed to support multiple datasets, for performing genome-wide correlative and quantitative analysis of ChIP-seq and RNA-seq data. Using data available from the ENCODE project, we illustrated several features of EpiMINE through different biological scenarios to show how easy some known observations can be verified. These results highlight how these approaches can be helpful in identifying novel biological features. CONCLUSIONS: EpiMINE performs different kinds of genome-wide quantitative and correlative analyses, using ChIP-seq- and RNA-seq-related datasets. Its framework enables it to be used by both experimental and computational researchers. EpiMINE can be downloaded from https://sourceforge.net/projects/epimine/.

Maintenance of leukemic cell identity by the activity of the Polycomb complex PRC1 in mice.

Leukemia is a complex heterogeneous disease often driven by the expression of oncogenic fusion proteins with different molecular and biochemical properties. Whereas several fusion proteins induce leukemogenesis by activating Hox gene expression (Hox-activating fusions), others impinge on different pathways that do not involve the activation of Hox genes (non-Hox-activating fusions). It has been postulated that one of the main oncogenic properties of the HOXA9 transcription factor is its ability to control the expression of the p16/p19 tumor suppressor locus (Cdkn2a), thereby compensating Polycomb-mediated repression, which is dispensable for leukemias induced by Hox-activating fusions. We show, by genetically depleting the H2A ubiquitin ligase subunits of the Polycomb repressive complex 1 (PRC1), Ring1a and Ring1b, that Hoxa9 activation cannot repress Cdkn2a expression in the absence of PRC1 and its dependent deposition of H2AK119 monoubiquitination (H2AK119Ub). This demonstrates the essential role of PRC1 activity in supporting the oncogenic potential of Hox-activating fusion proteins. By combining genetic tools with genome-wide location and transcription analyses, we further show that PRC1 activity is required for the leukemogenic potential of both Hox-activating and non-Hox-activating fusions, thus preventing the differentiation of leukemic cells independently of the expression of the Cdkn2a locus. Overall, our results genetically demonstrate that PRC1 activity and the deposition of H2AK119Ub are critical factors that maintain the undifferentiated identity of cancer cells, positively sustaining the progression of different types of leukemia.

PRC2 preserves intestinal progenitors and restricts secretory lineage commitment.

Chromatin modifications shape cell heterogeneity by activating and repressing defined sets of genes involved in cell proliferation, differentiation and development. Polycomb-repressive complexes (PRCs) act synergistically during development and differentiation by maintaining transcriptional repression of common genes. PRC2 exerts this activity by catalysing H3K27 trimethylation. Here, we show that in the intestinal epithelium PRC2 is required to sustain progenitor cell proliferation and the correct balance between secretory and absorptive lineage differentiation programs. Using genetic models, we show that PRC2 activity is largely dispensable for intestinal stem cell maintenance but is strictly required for radiation-induced regeneration by preventing Cdkn2a transcription. Combining these models with genomewide molecular analysis, we further demonstrate that preferential accumulation of secretory cells does not result from impaired proliferation of progenitor cells induced by Cdkn2a activation but rather from direct regulation of transcription factors responsible for secretory lineage commitment. Overall, our data uncover a dual role of PRC2 in intestinal homeostasis highlighting the importance of this repressive layer in controlling cell plasticity and lineage choices in adult tissues.

Polycomb Complex PRC1 Preserves Intestinal Stem Cell Identity by Sustaining Wnt/ß-Catenin Transcriptional Activity.

Polycomb repressive complexes (PRCs) are among the most important gatekeepers of establishing and maintaining cell identity in metazoans. PRC1, which plays a dominant role in this context, executes its functions via multiple subcomplexes, which all contribute to H2AK119 mono-ubiquitination (H2Aubq). Despite our comprehensive knowledge of PRC1-dependent H2Aubq in embryonic stem cells and during early development, its role in adult stem cells still remains poorly characterized. Here we show that PRC1 activity is required for the integrity of the intestinal epithelium, regulating stem cell self-renewal via a cell-autonomous mechanism that is independent from Cdkn2a expression. By dissecting the PRC1-dependent transcription program in intestinal stem cells, we demonstrate that PRC1 represses a large number of non-lineage-specific transcription factors that directly affect ß-catenin/Tcf transcriptional activity. Our data reveal that PRC1 preserves Wnt/ß-catenin activity in adult stem cells to maintain intestinal homeostasis and supports tumor formation induced by the constitutive activation of this pathway.

Polycomb proteins control proliferation and transformation independently of cell cycle checkpoints by regulating DNA replication.

The ability of PRC1 and PRC2 to promote proliferation is a main feature that links polycomb (PcG) activity to cancer. PcGs silence the expression of the tumour suppressor locus Ink4a/Arf, whose products positively regulate pRb and p53 functions. Enhanced PcG activity is a frequent feature of human tumours, and PcG inhibition has been proposed as a strategy for cancer treatment. However, the recurrent inactivation of pRb/p53 responses in human cancers raises a question regarding the ability of PcG proteins to affect cellular proliferation independently from this checkpoint. Here we demonstrate that PRCs regulate cellular proliferation and transformation independently of the Ink4a/Arf-pRb-p53 pathway. We provide evidence that PRCs localize at replication forks, and that loss of their function directly affects the progression and symmetry of DNA replication forks. Thus, we have identified a novel activity by which PcGs can regulate cell proliferation independently of major cell cycle restriction checkpoints.

Polycomb-dependent H3K27me1 and H3K27me2 regulate active transcription and enhancer fidelity.

H3K27me3 is deposited at promoters by the preferential association of Polycomb repressive complex 2 (PRC2) with CpG-rich DNA elements regulating development by repressing gene transcription. H3K27 is also present in mono- and dimethylated states; however, the functional roles of H3K27me1 and H3K27me2 deposition remain poorly characterized. Here, we show that PRC2 activity is not only associated with H3K27me3 but also regulates all forms of H3K27 methylation in a spatially defined manner, contributing to different genomic functions in mouse embryonic stem cells. H3K27me1 accumulates within transcribed genes, promotes transcription, and is regulated by Setd2-dependent H3K36me3 deposition. Contrarily, H3K27me2 is present on approximately 70% of total histone H3 and is distributed in large chromatin domains, exerting protective functions by preventing firing of non-cell-type-specific enhancers. Considering that only 5%-10% of deregulated genes in PRC2-deficient cells are direct H3K27me3 targets, our data support an active role for all H3K27 methylated forms in regulating transcription and determining cell identity.

Tet proteins connect the O-linked N-acetylglucosamine transferase Ogt to chromatin in embryonic stem cells.

O-linked N-acetylglucosamine (O-GlcNAc) transferase (Ogt) activity is essential for embryonic stem cell (ESC) viability and mouse development. Ogt is present both in the cytoplasm and the nucleus of different cell types and catalyzes serine and threonine glycosylation. We have characterized the biochemical features of nuclear Ogt and identified the ten-eleven translocation (TET) proteins Tet1 and Tet2 as stable partners of Ogt in the nucleus of ESCs. We show at a genome-wide level that Ogt preferentially associates with Tet1 to genes promoters in close proximity of CpG-rich transcription start sites. These regions are characterized by low levels of DNA modification, suggesting a link between Tet1 and Ogt activities in regulating CpG island methylation. Finally, we show that Tet1 is required for binding of Ogt to chromatin affecting Tet1 activity. Taken together, our data characterize how O-GlcNAcylation is recruited to chromatin and interacts with the activity of 5-methylcytosine hydroxylases.