Homopolymer switches mediate adaptive mutability in mismatch repair-deficient colorectal cancer.

Mismatch repair (MMR)-deficient cancer evolves through the stepwise erosion of coding homopolymers in target genes. Curiously, the MMR genes MutS homolog 6 (MSH6) and MutS homolog 3 (MSH3) also contain coding homopolymers, and these are frequent mutational targets in MMR-deficient cancers. The impact of incremental MMR mutations on MMR-deficient cancer evolution is unknown. Here we show that microsatellite instability modulates DNA repair by toggling hypermutable mononucleotide homopolymer runs in MSH6 and MSH3 through stochastic frameshift switching. Spontaneous mutation and reversion modulate subclonal mutation rate, mutation bias and HLA and neoantigen diversity. Patient-derived organoids corroborate these observations and show that MMR homopolymer sequences drift back into reading frame in the absence of immune selection, suggesting a fitness cost of elevated mutation rates. Combined experimental and simulation studies demonstrate that subclonal immune selection favors incremental MMR mutations. Overall, our data demonstrate that MMR-deficient colorectal cancers fuel intratumor heterogeneity by adapting subclonal mutation rate and diversity to immune selection.

Low coverage whole genome sequencing of low-grade dysplasia strongly predicts colorectal cancer risk in ulcerative colitis.

Patients with inflammatory bowel disease (IBD) are at increased risk of colorectal cancer (CRC), and this risk increases dramatically in those who develop low-grade dysplasia (LGD). However, there is currently no accurate way to risk-stratify patients with LGD, leading to both over- and under-treatment of cancer risk. Here we show that the burden of somatic copy number alterations (CNAs) within resected LGD lesions strongly predicts future cancer development. We performed a retrospective multi-centre validated case-control study of n=122 patients (40 progressors, 82 non-progressors, 270 LGD regions). Low coverage whole genome sequencing revealed CNA burden was significantly higher in progressors than non-progressors (p=2×10-6 in discovery cohort) and was a very significant predictor of CRC risk in univariate analysis (odds ratio = 36; p=9×10-7), outperforming existing clinical risk factors such as lesion size, shape and focality. Optimal risk prediction was achieved with a multivariate model combining CNA burden with the known clinical risk factor of incomplete LGD resection. The measurement of CNAs in LGD lesions is a robust, low-cost and rapidly translatable predictor of CRC risk in IBD that can be used to direct management and so prevent CRC in high-risk individuals whilst sparing those at low-risk from unnecessary intervention.

Profiling of Copy Number Alterations Using Low-Coverage Whole-Genome Sequencing Informs Differential Diagnosis and Prognosis in Primary Cutaneous Follicle Center Lymphoma.

Primary cutaneous follicle center lymphoma (PCFCL) has an excellent prognosis using local treatment, whereas nodal follicular lymphoma (nFL), occasionally presenting with cutaneous spread, often requires systemic therapy. Distinction of the 2 diseases based on histopathology alone might be challenging. Copy number alterations (CNAs) have scarcely been explored on a genome-wide scale in PCFCL; however, they might serve as potential biomarkers during differential diagnosis and risk stratification. Low-coverage whole-genome sequencing is a robust, high-throughput method for genome-wide copy number profiling. In this study, we analyzed 28 PCFCL samples from 20 patients and compared the copy number profiles with a cohort of diagnostic samples of 64 nFL patients. Although the copy number profile of PCFCL was similar to that of nFL, PCFCL lacked amplifications of 18q, with the frequency peaking at 18q21.33 in nFL cases involving the BCL2 locus (PCFCL: 5.0% vs nFL: 31.3%, P = .018, Fisher exact test). Development of distant cutaneous spread was significantly associated with higher genomic instability including the proportion of genome altered (0.02 vs 0.13, P = .033) and number of CNAs (2 vs 9 P = .017), as well as the enrichment of 2p22.2-p15 amplification involving REL and XPO1 (6.3% vs 60.0%, P = .005), 3q23-q24 amplification (0.0% vs 50.0%, P = .004), 6q16.1-q23.3 deletion (6.3% vs 50.0%, P = .018), and 9p21.3 deletion covering CDKN2A and CDKN2B loci (0.0% vs 40.0%, P = .014, all Fisher exact test) in PCFCL. Analysis of sequential tumor samples in 2 cases harboring an unfavorable clinical course pointed to the acquisition of 2p amplification in the earliest common progenitor underlining its pivotal role in malignant transformation. By performing genome-wide copy number profiling on the largest patient cohort to date, we identified distinctive CNA alterations conceivably facilitating the differential diagnosis of PCFCL and secondary cutaneous involvement of nFL and potentially aiding the risk stratification of patients with PCFCL in the future.

FUME-TCRseq Enables Sensitive and Accurate Sequencing of the T-cell Receptor from Limited Input of Degraded RNA.

Genomic analysis of the T-cell receptor (TCR) reveals the strength, breadth, and clonal dynamics of the adaptive immune response to pathogens or cancer. The diversity of the TCR repertoire, however, means that sequencing is technically challenging, particularly for samples with low-quality, degraded nucleic acids. Here, we developed and validated FUME-TCRseq, a robust and sensitive RNA-based TCR sequencing methodology that is suitable for formalin-fixed paraffin-embedded samples and low amounts of input material. FUME-TCRseq incorporates unique molecular identifiers into each molecule of cDNA, allowing correction for sequencing errors and PCR bias. Using RNA extracted from colorectal and head and neck cancers to benchmark the accuracy and sensitivity of FUME-TCRseq against existing methods demonstrated excellent concordance between the datasets. Furthermore, FUME-TCRseq detected more clonotypes than a commercial RNA-based alternative, with shorter library preparation time and significantly lower cost. The high sensitivity and the ability to sequence RNA of poor quality and limited amount enabled quantitative analysis of small numbers of cells from archival tissue sections, which is not possible with other methods. Spatially resolved FUME-TCRseq analysis of colorectal cancers using macrodissected archival samples revealed the shifting T-cell landscapes at the transition to an invasive phenotype and between tumor subclones containing distinct driver alterations. In summary, FUME-TCRseq represents an accurate, sensitive, and low-cost tool for the characterization of T-cell repertoires, particularly in samples with low-quality RNA that have not been accessible using existing methodology. SIGNIFICANCE: FUME-TCRseq is a TCR sequencing methodology that supports sensitive and spatially resolved detection of TCR clones in archival clinical specimens, which can facilitate longitudinal tracking of immune responses through disease course and treatment.

Epigenome and early selection determine the tumour-immune evolutionary trajectory of colorectal cancer.

Immune system control is a major hurdle that cancer evolution must circumvent. The relative timing and evolutionary dynamics of subclones that have escaped immune control remain incompletely characterized, and how immune-mediated selection shapes the epigenome has received little attention. Here, we infer the genome- and epigenome-driven evolutionary dynamics of tumour-immune coevolution within primary colorectal cancers (CRCs). We utilise our existing CRC multi-region multi-omic dataset that we supplement with high-resolution spatially-resolved neoantigen sequencing data and highly multiplexed imaging of the tumour microenvironment (TME). Analysis of somatic chromatin accessibility alterations (SCAAs) reveals frequent somatic loss of accessibility at antigen presenting genes, and that SCAAs contribute to silencing of neoantigens. We observe that strong immune escape and exclusion occur at the outset of CRC formation, and that within tumours, including at the microscopic level of individual tumour glands, additional immune escape alterations have negligible consequences for the immunophenotype of cancer cells. Further minor immuno-editing occurs during local invasion and is associated with TME reorganisation, but that evolutionary bottleneck is relatively weak. Collectively, we show that immune evasion in CRC follows a "Big Bang" evolutionary pattern, whereby genetic, epigenetic and TME-driven immune evasion acquired by the time of transformation defines subsequent cancer-immune evolution.

Molecular characteristics of early-onset pancreatic ductal adenocarcinoma.

The median age of patients with pancreatic ductal adenocarcinoma (PDAC) at diagnosis is 71 years; however, around 10% present with early-onset pancreatic cancer (EOPC), i.e., before age 50. The molecular mechanisms underlying such an early onset are unknown. We assessed the role of common PDAC drivers (KRAS, TP53, CDKN2A and SMAD4) and determined their mutational status and protein expression in 90 formalin-fixed, paraffin-embedded tissues, including multiple primary and matched metastases, from 37 EOPC patients. KRAS was mutated in 88% of patients; p53 was altered in 94%, and p16 and SMAD4 were lost in 86% and 71% of patients, respectively. Meta-synthesis showed a higher rate of p53 alterations in EOPC than in late-onset PDAC (94% vs. 69%, P = 0.0009) and significantly higher loss of SMAD4 (71% vs. 44%, P = 0.0025). The majority of EOPC patients accumulated aberrations in all four drivers; in addition, high tumour heterogeneity was observed across all tissues. The cumulative effect of an exceptionally high rate of alterations in all common PDAC driver genes combined with high tumour heterogeneity suggests an important mechanism underlying the early onset of PDAC.

Phenotypic noise and plasticity in cancer evolution.

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

Adaptive therapy achieves long-term control of chemotherapy resistance in high grade ovarian cancer.

Drug resistance results in poor outcomes for most patients with metastatic cancer. Adaptive Therapy (AT) proposes to address this by exploiting presumed fitness costs incurred by drug-resistant cells when drug is absent, and prescribing dose reductions to allow fitter, sensitive cells to re-grow and re-sensitise the tumour. However, empirical evidence for treatment-induced fitness change is lacking. We show that fitness costs in chemotherapy-resistant ovarian cancer cause selective decline and apoptosis of resistant populations in low-resource conditions. Moreover, carboplatin AT caused fluctuations in sensitive/resistant tumour population size in vitro and significantly extended survival of tumour-bearing mice. In sequential blood-derived cell-free DNA and tumour samples obtained longitudinally from ovarian cancer patients during treatment, we inferred resistant cancer cell population size through therapy and observed it correlated strongly with disease burden. These data have enabled us to launch a multicentre, phase 2 randomised controlled trial (ACTOv) to evaluate AT in ovarian cancer.

Evolutionary and immune microenvironment dynamics during neoadjuvant treatment of oesophagael adenocarcinoma.

Locally advanced oesophageal adenocarcinoma (EAC) remains difficult to treat because of common resistance to neoadjuvant therapy and high recurrence rates. The ecological and evolutionary dynamics responsible for treatment failure are incompletely understood. Here, we performed a comprehensive multi-omic analysis of samples collected from EAC patients in the MEMORI clinical trial, revealing major changes in gene expression profiles and immune microenvironment composition that did not appear to be driven by changes in clonal composition. Multi-region multi-timepoint whole exome (300x depth) and paired transcriptome sequencing was performed on 27 patients pre-, during and after neoadjuvant treatment. EAC showed major transcriptomic changes during treatment with upregulation of immune and stromal pathways and oncogenic pathways such as KRAS, Hedgehog and WNT. However, genetic data revealed that clonal sweeps were rare, suggesting that gene expression changes were not clonally driven. Additional longitudinal image mass cytometry was performed in a subset of 15 patients and T-cell receptor sequencing in 10 patients, revealing remodelling of the T-cell compartment during treatment and other shifts in microenvironment composition. The presence of immune escape mechanisms and a lack of clonal T-cell expansions were linked to poor clinical treatment response. This study identifies profound transcriptional changes during treatment with limited evidence that clonal replacement is the cause, suggesting phenotypic plasticity and immune dynamics as mechanisms for therapy resistance with pharmacological relevance.

First passage time analysis of spatial mutation patterns reveals sub-clonal evolutionary dynamics in colorectal cancer.

The signature of early cancer dynamics on the spatial arrangement of tumour cells is poorly understood, and yet could encode information about how sub-clones grew within the expanding tumour. Novel methods of quantifying spatial tumour data at the cellular scale are required to link evolutionary dynamics to the resulting spatial architecture of the tumour. Here, we propose a framework using first passage times of random walks to quantify the complex spatial patterns of tumour cell population mixing. First, using a simple model of cell mixing we demonstrate how first passage time statistics can distinguish between different pattern structures. We then apply our method to simulated patterns of mutated and non-mutated tumour cell population mixing, generated using an agent-based model of expanding tumours, to explore how first passage times reflect mutant cell replicative advantage, time of emergence and strength of cell pushing. Finally, we explore applications to experimentally measured human colorectal cancer, and estimate parameters of early sub-clonal dynamics using our spatial computational model. We infer a wide range of sub-clonal dynamics, with mutant cell division rates varying between 1 and 4 times the rate of non-mutated cells across our sample set. Some mutated sub-clones emerged after as few as 100 non-mutant cell divisions, and others only after 50,000 divisions. The majority were consistent with boundary driven growth or short-range cell pushing. By analysing multiple sub-sampled regions in a small number of samples, we explore how the distribution of inferred dynamics could inform about the initial mutational event. Our results demonstrate the efficacy of first passage time analysis as a new methodology in spatial analysis of solid tumour tissue, and suggest that patterns of sub-clonal mixing can provide insights into early cancer dynamics.

The co-evolution of the genome and epigenome in colorectal cancer.

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

Immunosuppressive niche engineering at the onset of human colorectal cancer.

The evolutionary dynamics of tumor initiation remain undetermined, and the interplay between neoplastic cells and the immune system is hypothesized to be critical in transformation. Colorectal cancer (CRC) presents a unique opportunity to study the transition to malignancy as pre-cancers (adenomas) and early-stage cancers are frequently resected. Here, we examine tumor-immune eco-evolutionary dynamics from pre-cancer to carcinoma using a computational model, ecological analysis of digital pathology data, and neoantigen prediction in 62 patient samples. Modeling predicted recruitment of immunosuppressive cells would be the most common driver of transformation. As predicted, ecological analysis reveals that progressed adenomas co-localized with immunosuppressive cells and cytokines, while benign adenomas co-localized with a mixed immune response. Carcinomas converge to a common immune "cold" ecology, relaxing selection against immunogenicity and high neoantigen burdens, with little evidence for PD-L1 overexpression driving tumor initiation. These findings suggest re-engineering the immunosuppressive niche may prove an effective immunotherapy in CRC.

Lineage tracing in human tissues.

The dynamical process of cell division that underpins homeostasis in the human body cannot be directly observed in vivo, but instead is measurable from the pattern of somatic genetic or epigenetic mutations that accrue in tissues over an individual's lifetime. Because somatic mutations are heritable, they serve as natural lineage tracing markers that delineate clonal expansions. Mathematical analysis of the distribution of somatic clone sizes gives a quantitative readout of the rates of cell birth, death, and replacement. In this review we explore the broad range of somatic mutation types that have been used for lineage tracing in human tissues, introduce the mathematical concepts used to infer dynamical information from these clone size data, and discuss the insights of this lineage tracing approach for our understanding of homeostasis and cancer development. We use the human colon as a particularly instructive exemplar tissue. There is a rich history of human somatic cell dynamics surreptitiously written into the cell genomes that is being uncovered by advances in sequencing and careful mathematical analysis lineage of tracing data. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Evolutionary dynamics in Barrett oesophagus: implications for surveillance, risk stratification and therapy.

Cancer development is a dynamic evolutionary process characterized by marked intratumoural heterogeneity at the genetic, epigenetic and phenotypic levels. Barrett oesophagus, the pre-malignant condition to oesophageal adenocarcinoma (EAC), is an exemplary system to longitudinally study the evolution of malignancy. Evidence has emerged of Barrett oesophagus lesions pre-programmed for progression to EAC many years before clinical detection, indicating a considerable window for therapeutic intervention. In this Review, we explore the mechanisms underlying clonal expansion and contraction that establish the Barrett oesophagus clonal mosaicism over time and space and discuss intrinsic genotypic and extrinsic environmental drivers that direct the evolutionary trajectory of Barrett oesophagus towards a malignant phenotype. We propose that understanding and exploiting the evolutionary dynamics of Barrett oesophagus will identify novel therapeutic targets, improve prognostic tools and offer the opportunity for personalized surveillance programmes geared to prevent progression to EAC.

Predicting Colorectal Cancer Occurrence in IBD.

Patients with colonic inflammatory bowel disease (IBD) are at an increased risk of developing colorectal cancer (CRC), and are therefore enrolled into a surveillance programme aimed at detecting dysplasia or early cancer. Current surveillance programmes are guided by clinical, endoscopic or histological predictors of colitis-associated CRC (CA-CRC). We have seen great progress in our understanding of these predictors of disease progression, and advances in endoscopic technique and management, along with improved medical care, has been mirrored by the falling incidence of CA-CRC over the last 50 years. However, more could be done to improve our molecular understanding of CA-CRC progression and enable better risk stratification for patients with IBD. This review summarises the known risk factors associated with CA-CRC and explores the molecular landscape that has the potential to complement and optimise the existing IBD surveillance programme.

The role of single-cell sequencing in studying tumour evolution.

Tumour evolution is a complex interplay between the acquisition of somatic (epi)genomic changes in tumour cells and the phenotypic consequences they cause, all in the context of a changing microenvironment. Single-cell sequencing offers a window into this dynamic process at the ultimate resolution of individual cells. In this review, we discuss the transformative insight offered by single-cell sequencing technologies for understanding tumour evolution.

Histological 3D reconstruction and in vivo lineage tracing of the human endometrium.

Regular menstrual shedding and repair of the endometrial functionalis is unique to humans and higher-order primates. The current consensus postulates endometrial glands to have a single-tubular architecture, where multi-potential stem cells reside in the blind-ending glandular-bases. Utilising fixed samples from patients, we have studied the three-dimensional (3D) micro-architecture of the human endometrium. We demonstrate that some non-branching, single, vertical functionalis glands originate from a complex horizontally interconnecting network of basalis glands. The existence of a multipotent endometrial epithelial stem cell capable of regenerating the entire complement of glandular lineages was demonstrated by in vivo lineage tracing, using naturally occurring somatic mitochondrial DNA mutations as clonal markers. Vertical tracking of mutated clones showed that at least one stem-cell population resides in the basalis glands. These novel findings provide insight into the efficient and scar-less regenerative potential of the human endometrium. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

In Situ Point Mutation Detection in FFPE Colorectal Cancers Using the BaseScope Assay.

In situ mutation detection (ISMD) is a powerful tool for the characterization of tumor heterogeneity at cellular resolution while preserving tissue morphology and spatial context. The BaseScope assay is a novel approach to ISMD, offering excellent specificity and sensitivity, with little requirement for assay optimization or technical expertise. Here we describe the validation and application of BaseScope ISMD probe sets to human formalin-fixed paraffin-embedded (FFPE) samples, firstly by testing the probes in well-characterized cell lines of known mutational status, and then by applying the assay to archival FFPE colorectal cancers.

Cancer associated fibroblast FAK regulates malignant cell metabolism.

Emerging evidence suggests that cancer cell metabolism can be regulated by cancer-associated fibroblasts (CAFs), but the mechanisms are poorly defined. Here we show that CAFs regulate malignant cell metabolism through pathways under the control of FAK. In breast and pancreatic cancer patients we find that low FAK expression, specifically in the stromal compartment, predicts reduced overall survival. In mice, depletion of FAK in a subpopulation of CAFs regulates paracrine signals that increase malignant cell glycolysis and tumour growth. Proteomic and phosphoproteomic analysis in our mouse model identifies metabolic alterations which are reflected at the transcriptomic level in patients with low stromal FAK. Mechanistically we demonstrate that FAK-depletion in CAFs increases chemokine production, which via CCR1/CCR2 on cancer cells, activate protein kinase A, leading to enhanced malignant cell glycolysis. Our data uncover mechanisms whereby stromal fibroblasts regulate cancer cell metabolism independent of genetic mutations in cancer cells.