Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

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.

Significance of Crypt Atypia in Barrett's Esophagus: A Clinical, Molecular, and Outcome Study.

BACKGROUND & AIMS: The aim of this study was to characterize baseline morphologic features of crypts in nondysplastic Barrett's esophagus and correlate them with DNA content abnormalities and risk of progression to high-grade dysplasia (HGD) or esophageal adenocarcinoma (EAC). METHODS: The morphologic features of nondysplastic crypts in baseline biopsy specimens from 212 BE patients (2956 biopsy specimens) were graded histologically using a 4-point scale (crypt atypia levels, 0-3). DNA content abnormalities were detected using flow cytometry. RESULTS: In patients who had dysplasia in their baseline biopsy specimens, dysplasia was associated significantly with increasing grades of crypt atypia in the background nondysplastic Barrett's esophagus (P < .001). In a subset of patients without dysplasia at baseline (N = 149), a higher grade of crypt atypia was associated with longer Barrett's esophagus segment length (5.5 vs 3.3 cm; P = .0095), and a higher percentage of cells with 4N DNA content (3.67 ± 1.27 vs 2.93 ± 1.22; P = .018). Crypt atypia was associated with the development of any neoplasia (low-grade dysplasia and HGD/EAC). Although no significant association was noted between the grade of crypt atypia and increased 4N, aneuploidy, or progression to HGD/EAC, only patients with grade 2 or 3 crypt atypia showed increased 4N, aneuploidy, or progression to HGD/EAC. CONCLUSIONS: Patients with Barrett's esophagus likely develop dysplasia via a progressive increase in the level of crypt atypia before the onset of dysplasia, and these changes may reflect some alteration of DNA content.

Genetic Insights in Barrett's Esophagus and Esophageal Adenocarcinoma.

Beginning in the 1980s, an alarming rise in the incidence of esophageal adenocarcinoma (EA) led to screening of patients with reflux to detect Barrett's esophagus (BE) and surveillance of BE to detect early EA. This strategy, based on linear progression disease models, resulted in selective detection of BE that does not progress to EA over a lifetime (overdiagnosis) and missed BE that rapidly progresses to EA (underdiagnosis). Here we review the historical thought processes that resulted in this undesired outcome and the transformation in our understanding of genetic and evolutionary principles governing neoplastic progression that has come from application of modern genomic technologies to cancers and their precursors. This new synthesis provides improved strategies for prevention and early detection of EA by addressing the environmental and mutational processes that can determine "windows of opportunity" in time to detect rapidly progressing BE and distinguish it from slowly or nonprogressing BE.

Studying Cancer Evolution in Barrett's Esophagus and Esophageal Adenocarcinoma.

Technological advances in genome sequencing and copy number analysis have allowed researchers to catalog the wide variety of genomic alterations that occur across diverse cancer types. For most cancer types, the lack of high-frequency alterations and the heterogeneity observed both within and between tumors suggest neoplastic progression proceeds through a branched evolutionary pathway as proposed by Nowell in 1976, as opposed to the linear pathway that has dominated medical science for the last century. To understand how cancer evolves over time and space in the body, new study designs are needed that can distinguish between alterations that develop in patients who progress to cancer from to those who don't. Here we present approaches developed in the study of Barrett's esophagus, a premalignant precursor of esophageal adenocarcinoma, and discuss strategies for applying the results from these analyses to address the critical clinical problems of overdiagnosis of benign disease, early detection of life-threatening cancer, and effective risk stratification.

NSAIDs modulate clonal evolution in Barrett's esophagus.

Cancer is considered an outcome of decades-long clonal evolution fueled by acquisition of somatic genomic abnormalities (SGAs). Non-steroidal anti-inflammatory drugs (NSAIDs) have been shown to reduce cancer risk, including risk of progression from Barrett's esophagus (BE) to esophageal adenocarcinoma (EA). However, the cancer chemopreventive mechanisms of NSAIDs are not fully understood. We hypothesized that NSAIDs modulate clonal evolution by reducing SGA acquisition rate. We evaluated thirteen individuals with BE. Eleven had not used NSAIDs for 6.2±3.5 (mean±standard deviation) years and then began using NSAIDs for 5.6±2.7 years, whereas two had used NSAIDs for 3.3±1.4 years and then discontinued use for 7.9±0.7 years. 161 BE biopsies, collected at 5-8 time points over 6.4-19 years, were analyzed using 1Million-SNP arrays to detect SGAs. Even in the earliest biopsies there were many SGAs (284±246 in 10/13 and 1442±560 in 3/13 individuals) and in most individuals the number of SGAs changed little over time, with both increases and decreases in SGAs detected. The estimated SGA rate was 7.8 per genome per year (95% support interval [SI], 7.1-8.6) off-NSAIDs and 0.6 (95% SI 0.3-1.5) on-NSAIDs. Twelve individuals did not progress to EA. In ten we detected 279±86 SGAs affecting 53±30 Mb of the genome per biopsy per time point and in two we detected 1,463±375 SGAs affecting 180±100 Mb. In one individual who progressed to EA we detected a clone having 2,291±78 SGAs affecting 588±18 Mb of the genome at three time points in the last three of 11.4 years of follow-up. NSAIDs were associated with reduced rate of acquisition of SGAs in eleven of thirteen individuals. Barrett's cells maintained relative equilibrium level of SGAs over time with occasional punctuations by expansion of clones having massive amount of SGAs.

Genetic clonal diversity predicts progression to esophageal adenocarcinoma.

Neoplasms are thought to progress to cancer through genetic instability generating cellular diversity and clonal expansions driven by selection for mutations in cancer genes. Despite advances in the study of molecular biology of cancer genes, relatively little is known about evolutionary mechanisms that drive neoplastic progression. It is unknown, for example, which may be more predictive of future progression of a neoplasm: genetic homogenization of the neoplasm, possibly caused by a clonal expansion, or the accumulation of clonal diversity. Here, in a prospective study, we show that clonal diversity measures adapted from ecology and evolution can predict progression to adenocarcinoma in the premalignant condition known as Barrett's esophagus, even when controlling for established genetic risk factors, including lesions in TP53 (p53; ref. 6) and ploidy abnormalities. Progression to cancer through accumulation of clonal diversity, on which natural selection acts, may be a fundamental principle of neoplasia with important clinical implications.

Breakage fusion bridge cycles drive high oncogene copy number, but not intratumoral genetic heterogeneity or rapid cancer genome change.

Oncogene amplification is a major driver of cancer pathogenesis. Breakage fusion bridge (BFB) cycles, like extrachromosomal DNA (ecDNA), can lead to high copy numbers of oncogenes, but their impact on intratumoral heterogeneity, treatment response, and patient survival are not well understood due to difficulty in detecting them by DNA sequencing. We describe a novel algorithm that detects and reconstructs BFB amplifications using optical genome maps (OGMs), called OM2BFB. OM2BFB showed high precision (>93%) and recall (92%) in detecting BFB amplifications in cancer cell lines, PDX models and primary tumors. OM-based comparisons demonstrated that short-read BFB detection using our AmpliconSuite (AS) toolkit also achieved high precision, albeit with reduced sensitivity. We detected 371 BFB events using whole genome sequences from 2,557 primary tumors and cancer lines. BFB amplifications were preferentially found in cervical, head and neck, lung, and esophageal cancers, but rarely in brain cancers. BFB amplified genes show lower variance of gene expression, with fewer options for regulatory rewiring relative to ecDNA amplified genes. BFB positive (BFB (+)) tumors showed reduced heterogeneity of amplicon structures, and delayed onset of resistance, relative to ecDNA(+) tumors. EcDNA and BFB amplifications represent contrasting mechanisms to increase the copy numbers of oncogene with markedly different characteristics that suggest different routes for intervention.

Somatic whole genome dynamics of precancer in Barrett's esophagus reveals features associated with disease progression.

While the genomes of normal tissues undergo dynamic changes over time, little is understood about the temporal-spatial dynamics of genomes in premalignant tissues that progress to cancer compared to those that remain cancer-free. Here we use whole genome sequencing to contrast genomic alterations in 427 longitudinal samples from 40 patients with stable Barrett's esophagus compared to 40 Barrett's patients who progressed to esophageal adenocarcinoma (ESAD). We show the same somatic mutational processes are active in Barrett's tissue regardless of outcome, with high levels of mutation, ESAD gene and focal chromosomal alterations, and similar mutational signatures. The critical distinction between stable Barrett's versus those who progress to cancer is acquisition and expansion of TP53-/- cell populations having complex structural variants and high-level amplifications, which are detectable up to six years prior to a cancer diagnosis. These findings reveal the timing of common somatic genome dynamics in stable Barrett's esophagus and define key genomic features specific to progression to esophageal adenocarcinoma, both of which are critical for cancer prevention and early detection strategies.

Within-patient phylogenetic reconstruction reveals early events in Barrett's Esophagus.

Barrett's Esophagus is a neoplastic condition which progresses to esophageal adenocarcinoma in 5% of cases. Key events affecting the outcome likely occur before diagnosis of Barrett's and cannot be directly observed; we use phylogenetic analysis to infer such past events. We performed whole-genome sequencing on 4-6 samples from 40 cancer outcome and 40 noncancer outcome patients with Barrett's Esophagus, and inferred within-patient phylogenies of deconvoluted clonal lineages. Spatially proximate lineages clustered in the phylogenies, but temporally proximate ones did not. Lineages with inferred loss-of-function mutations in both copies of TP53 and CDKN2A showed enhanced spatial spread, whereas lineages with loss-of-function mutations in other frequently mutated loci did not. We propose a two-phase model with expansions of TP53 and CKDN2A mutant lineages during initial growth of the segment, followed by relative stasis. Subsequent to initial expansion, mutations in these loci as well as ARID1A and SMARCA4 may show a local selective advantage but do not expand far: The spatial structure of the Barrett's segment remains stable during surveillance even in patients who go on to cancer. We conclude that the cancer/noncancer outcome is strongly affected by early steps in formation of the Barrett's segment.

Chromosomal instability and copy number alterations in Barrett's esophagus and esophageal adenocarcinoma.

PURPOSE: Chromosomal instability, as assessed by many techniques, including DNA content aneuploidy, loss of heterozygosity, and comparative genomic hybridization, has consistently been reported to be common in cancer and rare in normal tissues. Recently, a panel of chromosome instability biomarkers, including loss of heterozygosity and DNA content, has been reported to identify patients at high and low risk of progression from Barrett's esophagus (BE) to esophageal adenocarcinoma (EA), but required multiple platforms for implementation. Although chromosomal instability involving amplifications and deletions of chromosome regions have been observed in nearly all cancers, copy number alterations (CNA) in premalignant tissues have not been well characterized or evaluated in cohort studies as biomarkers of cancer risk. EXPERIMENTAL DESIGN: We examined CNAs in 98 patients having either BE or EA using Bacterial Artificial Chromosome (BAC) array comparative genomic hybridization to characterize CNAs at different stages of progression ranging from early BE to advanced EA. RESULTS: CNAs were rare in early stages (less than high-grade dysplasia) but were progressively more frequent and larger in later stages (high-grade dysplasia and EA), including high-level amplifications. The number of CNAs correlated highly with DNA content aneuploidy. Patients whose biopsies contained CNAs involving >70 Mbp were at increased risk of progression to DNA content abnormalities or EA (hazards ratio, 4.9; 95% confidence interval, 1.6-14.8; P = 0.0047), and the risk increased as more of the genome was affected. CONCLUSIONS: Genome-wide analysis of CNAs provides a common platform for the evaluation of chromosome instability for cancer risk assessment as well as for the identification of common regions of alteration that can be further studied for biomarker discovery.

Cell proliferation, cell cycle abnormalities, and cancer outcome in patients with Barrett's esophagus: a long-term prospective study.

PURPOSE: Elevated cellular proliferation and cell cycle abnormalities, which have been associated with premalignant lesions, may be caused by inactivation of tumor suppressor genes. We measured proliferative and cell cycle fractions of biopsies from a cohort of patients with Barrett's esophagus to better understand the role of proliferation in early neoplastic progression and the association between cell cycle dysregulation and tumor suppressor gene inactivation. EXPERIMENTAL DESIGN: Cell proliferative fractions (determined by Ki67/DNA multiparameter flow cytometry) and cell cycle fractions (DNA content flow cytometry) were measured in 853 diploid biopsies from 362 patients with Barrett's esophagus. The inactivation status of CDKN2A and TP53 was assessed in a subset of these biopsies in a cross-sectional study. A prospective study followed 276 of the patients without detectable aneuploidy for an average of 6.3 years with esophageal adenocarcinoma as an end point. RESULTS: Diploid S and 4N (G(2)/tetraploid) fractions were significantly higher in biopsies with TP53 mutation and loss of heterozygosity. CDKN2A inactivation was not associated with higher Ki67-positive, diploid S, G(1), or 4N fractions. High Ki67-positive and G(1)-phase fractions were not associated with the future development of esophageal adenocarcinoma (P=0.13 and P=0.15, respectively), whereas high diploid S-phase and 4N fractions were (P=0.03 and P<0.0001, respectively). CONCLUSIONS: High Ki67-positive proliferative fractions were not associated with inactivation of CDKN2A and TP53 or future development of cancer in our cohort of patients with Barrett's esophagus. Biallelic inactivation of TP53 was associated with elevated 4N fractions, which have been associated with the future development of esophageal adenocarcinoma.

Correction to: NSAID use and somatic exomic mutations in Barrett's esophagus.

It was highlighted that in the original article [1] the Availability of data and materials section was incorrect.

Direct inference of SNP heterozygosity rates and resolution of LOH detection.

Single nucleotide polymorphisms (SNPs) have been increasingly utilized to investigate somatic genetic abnormalities in premalignancy and cancer. LOH is a common alteration observed during cancer development, and SNP assays have been used to identify LOH at specific chromosomal regions. The design of such studies requires consideration of the resolution for detecting LOH throughout the genome and identification of the number and location of SNPs required to detect genetic alterations in specific genomic regions. Our study evaluated SNP distribution patterns and used probability models, Monte Carlo simulation, and real human subject genotype data to investigate the relationships between the number of SNPs, SNP HET rates, and the sensitivity (resolution) for detecting LOH. We report that variances of SNP heterozygosity rate in dbSNP are high for a large proportion of SNPs. Two statistical methods proposed for directly inferring SNP heterozygosity rates require much smaller sample sizes (intermediate sizes) and are feasible for practical use in SNP selection or verification. Using HapMap data, we showed that a region of LOH greater than 200 kb can be reliably detected, with losses smaller than 50 kb having a substantially lower detection probability when using all SNPs currently in the HapMap database. Higher densities of SNPs may exist in certain local chromosomal regions that provide some opportunities for reliably detecting LOH of segment sizes smaller than 50 kb. These results suggest that the interpretation of the results from genome-wide scans for LOH using commercial arrays need to consider the relationships among inter-SNP distance, detection probability, and sample size for a specific study. New experimental designs for LOH studies would also benefit from considering the power of detection and sample sizes required to accomplish the proposed aims.

Evolution of Barrett's esophagus through space and time at single-crypt and whole-biopsy levels.

The low risk of progression of Barrett's esophagus (BE) to esophageal adenocarcinoma can lead to over-diagnosis and over-treatment of BE patients. This may be addressed through a better understanding of the dynamics surrounding BE malignant progression. Although genetic diversity has been characterized as a marker of malignant development, it is still unclear how BE arises and develops. Here we uncover the evolutionary dynamics of BE at crypt and biopsy levels in eight individuals, including four patients that experienced malignant progression. We assay eight individual crypts and the remaining epithelium by SNP array for each of 6-11 biopsies over 2 time points per patient (358 samples in total). Our results indicate that most Barrett's segments are clonal, with similar number and inferred rates of alterations observed for crypts and biopsies. Divergence correlates with geographical location, being higher near the gastro-esophageal junction. Relaxed clock analyses show that genomic instability precedes and is enhanced by genome doubling. These results shed light on the clinically relevant evolutionary dynamics of BE.

NSAID use and somatic exomic mutations in Barrett's esophagus.

BACKGROUND: Use of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) has been shown to protect against tetraploidy, aneuploidy, and chromosomal alterations in the metaplastic condition Barrett's esophagus (BE) and to lower the incidence and mortality of esophageal adenocarcinoma (EA). The esophagus is exposed to both intrinsic and extrinsic mutagens resulting from gastric reflux, chronic inflammation, and exposure to environmental carcinogens such as those found in cigarettes. Here we test the hypothesis that NSAID use inhibits accumulation of point mutations/indels during somatic genomic evolution in BE. METHODS: Whole exome sequences were generated from 82 purified epithelial biopsies and paired blood samples from a cross-sectional study of 41 NSAID users and 41 non-users matched by sex, age, smoking, and continuous time using or not using NSAIDs. RESULTS: NSAID use reduced overall frequency of point mutations across the spectrum of mutation types, lowered the frequency of mutations even when adjusted for both TP53 mutation and smoking status, and decreased the prevalence of clones with high variant allele frequency. Never smokers who consistently used NSAIDs had fewer point mutations in signature 17, which is commonly found in EA. NSAID users had, on average, a 50% reduction in functional gene mutations in nine cancer-associated pathways and also had less diversity in pathway mutational burden compared to non-users. CONCLUSIONS: These results indicate NSAID use functions to limit overall mutations on which selection can act and supports a model in which specific mutant cell populations survive or expand better in the absence of NSAIDs.

Neosquamous epithelium does not typically arise from Barrett's epithelium.

PURPOSE: Neosquamous epithelium (NSE) can arise within Barrett's esophagus as a consequence of medical or surgical acid reduction therapy, as well as after endoscopic ablation. Morphologic studies have suggested that NSE can develop from adjacent squamous epithelium, submucosal gland ducts, or multipotent progenitor cell(s) that can give rise to either squamous or Barrett's epithelium, depending on the luminal environment. The cells responsible for Barrett's epithelium self-renewal are frequently mutated during neoplastic progression. If NSE arises from the same cells that self-renew the Barrett's epithelium, the two tissues should be clonally related and share genetic alterations; if NSE does not originate in the self-renewing Barrett's, NSE and Barrett's esophagus should be genetically independent. EXPERIMENTAL DESIGN: We isolated islands of NSE and the surrounding Barrett's epithelium from 20 patients by microdissection and evaluated each tissue for genetic alterations in exon 2 of CDKN2A or exons 5 to 9 of the TP53 gene. Nine patients had p16 mutations and 11 had TP53 mutations within the Barrett's epithelium. RESULTS: In 1 of 20 patients, a focus of NSE had a 146 bp deletion in p16 identical to that found in surrounding Barrett's epithelium. The NSE in the remaining 19 patients was wild-type for p16 or TP53. CONCLUSION: Our mutational data support the hypothesis that, in most circumstances, NSE originates in cells different from those responsible for self-renewal of Barrett's epithelium. However, in one case, NSE and Barrett's epithelium seem to have arisen from a progenitor cell that was capable of differentiating into either intestinal metaplasia or NSE.

Quantification of Multiple Tumor Clones Using Gene Array and Sequencing Data.

Cancer development is driven by genomic alterations, including copy number aberrations. The detection of copy number aberrations in tumor cells is often complicated by possible contamination of normal stromal cells in tumor samples and intratumor heterogeneity, namely the presence of multiple clones of tumor cells. In order to correctly quantify copy number aberrations, it is critical to successfully de-convolute the complex structure of the genetic information from tumor samples. In this article, we propose a general Bayesian method for estimating copy number aberrations when there are normal cells and potentially more than one tumor clones. Our method provides posterior probabilities for the proportions of tumor clones and normal cells. We incorporate prior information on the distribution of the copy numbers to prioritize biologically more plausible solutions and alleviate possible identifiability issues that have been observed by many researchers. Our model is flexible and can work for both SNP array and next-generation sequencing data. We compare our method to existing ones and illustrate the advantage of our approach in multiple datasets.

Mutagen sensitivity and neoplastic progression in patients with Barrett's esophagus: a prospective analysis.

BACKGROUND: Defects in DNA damage recognition and repair have been associated with a wide variety of cancers. We conducted a prospective study to determine whether mutagen sensitivity, as determined by an in vitro assay, was associated with the future development of cancer in patients with Barrett's esophagus, which is associated with increased risk of progression to esophageal adenocarcinoma. METHODS: We measured sensitivity to bleomycin in peripheral blood lymphocytes in a cohort of 220 patients with Barrett's esophagus. We followed these patients for 1,230 person-years (range, 3 months to 10.1 years; median, 6.4 years), using development of cancer and aneuploidy as end points. A subset of these patients was evaluated for inactivation of tumor-suppressor genes CDKN2A/p16 and TP53 [by mutation and loss of heterozygosity (LOH)] in their Barrett's segments at the time of, or before, the bleomycin test, and the patients were stratified by CDKN2A/p16 and TP53 status in an analysis of mutagen sensitivity and progression. RESULTS: Bleomycin-sensitive patients were found to be at significantly greater risk of developing aneuploidy (adjusted hazard ratio, 3.71; 95% confidence interval, 1.44-9.53) and nonsignificantly greater risk of cancer (adjusted hazard ratio, 1.63; 95% confidence interval, 0.71-3.75). Among patients with detectable LOH at the TP53 locus (on chromosome 17p), increasing bleomycin sensitivity was associated with increased risk of developing cancer (P(trend) < 0.001) and aneuploidy (P(trend) = 0.005). CONCLUSIONS: This study supports the hypothesis that sensitivity to mutagens increases the risk of neoplastic progression in persons with Barrett's esophagus, particularly those with 17p LOH including TP53.