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

The somatic mutation landscape of normal gastric epithelium

The landscapes of somatic mutation in normal cells inform us about the processes of mutation and selection operative throughout life, providing insight into normal ageing and the earliest stages of cancer development1. Here, by whole-genome sequencing of 238 microdissections2 from 30 individuals, including 18 with gastric cancer, we elucidate the developmental trajectories of normal and malignant gastric epithelium. We find that gastric glands are units of monoclonal cell populations that accrue roughly 28 somatic single-nucleotide variants per year, predominantly attributable to endogenous mutational processes. In individuals with gastric cancer, metaplastic glands often show elevated mutation burdens due to acceleration of mutational processes linked to proliferation and oxidative damage. Unusually for normal cells, gastric epithelial cells often carry recurrent trisomies of specific chromosomes, which are highly enriched in a subset of individuals. Surveying 829 polyclonal gastric microbiopsies by targeted sequencing, we find somatic 'driver' mutations in a distinctive repertoire of known cancer genes, including ARID1A, ARID1B, ARID2, CTNNB1 and KDM6A. The prevalence of mutant clones increases with age to occupy roughly 8% of the gastric epithelial lining by age 60 years and is significantly increased by the presence of severe chronic inflammation. Our findings provide insights into intrinsic and extrinsic influences on somatic evolution in the gastric epithelium in healthy, precancerous and malignant states.

p53 mutation biases squamocolumnar junction progenitor cells towards dysplasia rather than metaplasia in Barrett's oesophagus

BACKGROUND

While p53 mutations occur early in Barrett's oesophagus (BE) progression to oesophageal adenocarcinoma (EAC), their role in gastric cardia stem cells remains unclear.

OBJECTIVE

This study investigates the impact of p53 mutation on the fate and function of cardia progenitor cells in BE to EAC progression, particularly under the duress of chronic injury.

DESIGN

We used a BE mouse model (L2-IL1β) harbouring a Trp53 mutation (R172H) to study the effects of p53 on Cck2r+ cardia progenitor cells. We employed lineage tracing, pathological analysis, organoid cultures, single-cell RNA sequencing (scRNA-seq) and computational analyses to investigate changes in progenitor cell behaviour, differentiation patterns and tumour progression. Additionally, we performed orthotopic transplantation of sorted metaplastic and mutant progenitor cells to assess their tumourigenic potential in vivo.

RESULTS

The p53 mutation acts as a switch to expand progenitor cells and inhibit their differentiation towards metaplasia, but only amidst chronic injury. In L2-IL1β mice, p53 mutation increased progenitors expansion and lineage-tracing with a shift from metaplasia to dysplasia. scRNA-seq revealed dysplastic cells arise directly from mutant progenitors rather than progressing through metaplasia. In vitro, p53 mutation enhanced BE progenitors' organoid-forming efficiency, growth, DNA damage resistance and progression to aneuploidy. Sorted metaplastic cells grew poorly with no progression to dysplasia, while mutant progenitors gave rise to dysplasia in orthotopic transplantation. Computational analyses indicated that p53 mutation inhibited stem cell differentiation through Notch activation.

CONCLUSIONS

p53 mutation contributes to BE progression by increasing expansion and fitness of undifferentiated cardia progenitors and preventing their differentiation towards metaplasia.

Community assembly modeling of microbial evolution within Barrett's esophagus and esophageal adenocarcinoma

Mathematical modeling of somatic evolution, a process impacting both host cells and microbial communities in the human body, can capture important dynamics driving carcinogenesis. Here we considered models for esophageal adenocarcinoma (EAC), a cancer that has dramatically increased in incidence over the past few decades in Western populations, with high case fatality rates due to late-stage diagnoses. Despite advancements in genomic analyses of the precursor Barrett's esophagus (BE), prevention of late-stage EAC remains a significant clinical challenge. Previous microbiome studies in BE and EAC have focused on quantifying static microbial abundance differences rather than evolutionary dynamics. Using whole genome sequencing data from esophageal tissues, we first applied a robust bioinformatics pipeline to extract non-host DNA reads, mapped these putative reads to microbial taxa, and retained those taxa with high genomic coverage. When applying mathematical models of microbial evolution to sequential stages of progression to EAC, we observed evidence of neutral dynamics in community assembly within normal esophageal tissue and BE, but not EAC. In a case-control study of BE patients who progressed to EAC cancer outcomes (CO) versus those who had non-cancer outcomes (NCO) during follow-up (mean=10.5 years), we found that Helicobacter pylori deviated significantly from the neutral expectation in BE NCO, suggesting that factors related to H. pylori or H. pylori infection itself may influence EAC risk. Additionally, simulations incorporating selection recapitulated non-neutral behaviors observed in the datasets. Formally modeling dynamics during progression holds promise in clinical applications by offering a deeper understanding of microbial involvement in cancer development.

Context-dependent effects of CDKN2A and other 9p21 gene losses during the evolution of esophageal cancer

CDKN2A is a tumor suppressor located in chromosome 9p21 and frequently lost in Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC). How CDKN2A and other 9p21 gene co-deletions affect EAC evolution remains understudied. We explored the effects of 9p21 loss in EACs and cancer progressor and non-progressor BEs with matched genomic, transcriptomic and clinical data. Despite its cancer driver role, CDKN2A loss in BE prevents EAC initiation by counterselecting subsequent TP53 alterations. 9p21 gene co-deletions predict poor patient survival in EAC but not BE through context-dependent effects on cell cycle, oxidative phosphorylation and interferon response. Immune quantifications using bulk transcriptome, RNAscope and high-dimensional tissue imaging showed that IFNE loss reduces immune infiltration in BE, but not EAC. Mechanistically, CDKN2A loss suppresses the maintenance of squamous epithelium, contributing to a more aggressive phenotype. Our study demonstrates context-dependent roles of cancer genes during disease evolution, with consequences for cancer detection and patient management.

Tumor initiation and early tumorigenesis: molecular mechanisms and interventional targets

Tumorigenesis is a multistep process, with oncogenic mutations in a normal cell conferring clonal advantage as the initial event. However, despite pervasive somatic mutations and clonal expansion in normal tissues, their transformation into cancer remains a rare event, indicating the presence of additional driver events for progression to an irreversible, highly heterogeneous, and invasive lesion. Recently, researchers are emphasizing the mechanisms of environmental tumor risk factors and epigenetic alterations that are profoundly influencing early clonal expansion and malignant evolution, independently of inducing mutations. Additionally, clonal evolution in tumorigenesis reflects a multifaceted interplay between cell-intrinsic identities and various cell-extrinsic factors that exert selective pressures to either restrain uncontrolled proliferation or allow specific clones to progress into tumors. However, the mechanisms by which driver events induce both intrinsic cellular competency and remodel environmental stress to facilitate malignant transformation are not fully understood. In this review, we summarize the genetic, epigenetic, and external driver events, and their effects on the co-evolution of the transformed cells and their ecosystem during tumor initiation and early malignant evolution. A deeper understanding of the earliest molecular events holds promise for translational applications, predicting individuals at high-risk of tumor and developing strategies to intercept malignant transformation.

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.

Decoding the basis of histological variation in human cancer

Molecular abnormalities that shape human neoplasms dissociate their phenotypic landscape from that of the healthy counterpart. Through the lens of a microscope, tumour pathology optically captures such aberrations projected onto a tissue slide and has categorized human epithelial neoplasms into distinct histological subtypes based on the diverse morphogenetic and molecular programmes that they manifest. Tumour histology often reflects tumour aggressiveness, patient prognosis and therapeutic vulnerability, and thus has been used as a de facto diagnostic tool and for making clinical decisions. However, it remains elusive how the diverse histological subtypes arise and translate into pleiotropic biological phenotypes. Molecular analysis of clinical tumour tissues and their culture, including patient-derived organoids, and add-back genetic reconstruction of tumorigenic pathways using gene engineering in culture models and rodents further elucidated molecular mechanisms that underlie morphological variations. Such mechanisms include genetic mutations and epigenetic alterations in cellular identity codes that erode hard-wired morphological programmes and histologically digress tumours from the native tissues. Interestingly, tumours acquire the ability to grow independently of the niche-driven stem cell ecosystem along with these morphological alterations, providing a biological rationale for histological diversification during tumorigenesis. This Review comprehensively summarizes our current understanding of such plasticity in the histological and lineage commitment fostered cooperatively by molecular alterations and the tumour environment, and describes basic and clinical implications for future cancer therapy.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Most large structural variants in cancer genomes can be detected without long reads

Short-read sequencing is the workhorse of cancer genomics yet is thought to miss many structural variants (SVs), particularly large chromosomal alterations. To characterize missing SVs in short-read whole genomes, we analyzed 'loose ends'-local violations of mass balance between adjacent DNA segments. In the landscape of loose ends across 1,330 high-purity cancer whole genomes, most large (>10-kb) clonal SVs were fully resolved by short reads in the 87% of the human genome where copy number could be reliably measured. Some loose ends represent neotelomeres, which we propose as a hallmark of the alternative lengthening of telomeres phenotype. These pan-cancer findings were confirmed by long-molecule profiles of 38 breast cancer and melanoma cases. Our results indicate that aberrant homologous recombination is unlikely to drive the majority of large cancer SVs. Furthermore, analysis of mass balance in short-read whole genome data provides a surprisingly complete picture of cancer chromosomal structure.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

Deterministic evolution and stringent selection during preneoplasia

The earliest events during human tumour initiation, although poorly characterized, may hold clues to malignancy detection and prevention1. Here we model occult preneoplasia by biallelic inactivation of TP53, a common early event in gastric cancer, in human gastric organoids. Causal relationships between this initiating genetic lesion and resulting phenotypes were established using experimental evolution in multiple clonally derived cultures over 2 years. TP53 loss elicited progressive aneuploidy, including copy number alterations and structural variants prevalent in gastric cancers, with evident preferred orders. Longitudinal single-cell sequencing of TP53-deficient gastric organoids similarly indicates progression towards malignant transcriptional programmes. Moreover, high-throughput lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programmes repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a marked degree of phenotypic convergence in premalignant epithelial organoids. These data imply predictability in the earliest stages of tumorigenesis and show evolutionary constraints and barriers to malignant transformation, with implications for earlier detection and interception of aggressive, genome-instable tumours.

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.

AI-assisted discovery of an ethnicity-influenced driver of cell transformation in esophageal and gastroesophageal junction adenocarcinomas

Although Barrett's metaplasia of the esophagus (BE) is the only known precursor lesion to esophageal adenocarcinomas (EACs), drivers of cellular transformation in BE remain incompletely understood. We use an artificial intelligence-guided network approach to study EAC initiation and progression. Key predictions are subsequently validated in a human organoid model, in patient-derived biopsy specimens of BE, a case-control study of genomics of BE progression, and in a cross-sectional study of 113 patients with BE and EACs. Our model classified healthy esophagus from BE and BE from EACs in several publicly available gene expression data sets (n = 932 samples). The model confirmed that all EACs must originate from BE and pinpointed a CXCL8/IL8↔neutrophil immune microenvironment as a driver of cellular transformation in EACs and gastroesophageal junction adenocarcinomas. This driver is prominent in White individuals but is notably absent in African Americans (AAs). Network-derived gene signatures, independent signatures of neutrophil processes, CXCL8/IL8 expression, and an absolute neutrophil count (ANC) are associated with risk of progression. SNPs associated with changes in ANC by ethnicity (e.g., benign ethnic neutropenia [BEN]) modify that risk. Findings define a racially influenced immunological basis for cell transformation and suggest that BEN in AAs may be a deterrent to BE→EAC progression.

Characterizing isoform switching events in esophageal adenocarcinoma

Isoform switching events with predicted functional consequences are common in many cancers, but characterization of switching events in esophageal adenocarcinoma (EAC) is lacking. Next-generation sequencing was used to detect levels of RNA transcripts and identify specific isoforms in treatment-naïve esophageal tissues ranging from premalignant Barrett’s esophagus (BE), BE with low- or high-grade dysplasia (BE.LGD, BE.HGD), and EAC. Samples were stratified by histopathology and TP53 mutation status, identifying significant isoform switching events with predicted functional consequences. Comparing BE.LGD with BE.HGD, a histopathology linked to cancer progression, isoform switching events were identified in 75 genes including KRAS, RNF128, and WRAP53. Stratification based on TP53 status increased the number of significant isoform switches to 135, suggesting switching events affect cellular functions based on TP53 mutation and tissue histopathology. Analysis of isoforms agnostic, exclusive, and shared with mutant TP53 revealed unique signatures including demethylation, lipid and retinoic acid metabolism, and glucuronidation, respectively. Nearly half of isoform switching events were identified without significant gene-level expression changes. Importantly, two TP53-interacting isoforms, RNF128 and WRAP53, were significantly linked to patient survival. Thus, analysis of isoform switching events may provide new insight for the identification of prognostic markers and inform new potential therapeutic targets for EAC.