BACKGROUND:The risk of developing advanced neoplasia (AN; colorectal cancer and/or high-grade dysplasia) in ulcerative colitis (UC) patients with a low-grade dysplasia (LGD) lesion is variable and difficult to predict. This is a major challenge for effective clinical management. OBJECTIVE:We aimed to provide accurate AN risk stratification in UC patients with LGD. We hypothesised that the pattern and burden of somatic genomic copy number alterations (CNAs) in LGD lesions could predict future AN risk. DESIGN:We performed a retrospective multicentre validated case-control study using n=270 LGD samples from n=122 patients with UC. Patients were designated progressors (n=40) if they had a diagnosis of AN in the ~5 years following LGD diagnosis or non-progressors (n=82) if they remained AN-free during follow-up. DNA was extracted from the baseline LGD lesion, low-coverage whole genome sequencing performed and data processed to detect CNAs. Survival analysis was used to evaluate CNAs as predictors of future AN risk. RESULTS:CNA burden was significantly higher in progressors than non-progressors (p=2×10-6 in discovery cohort) and was a very significant predictor of AN risk in univariate analysis (OR=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. Within-LGD lesion genetic heterogeneity did not confound risk prediction. CONCLUSION:Measurement of CNAs in LGD is an accurate predictor of AN risk in inflammatory bowel disease and is likely to support clinical management.

Digistain measures the elevated DNA levels in Cancer biopsies using mid-IR transmission. We show, in a 801 patient trial, that it provides accurate prognostic indicators, that can save lives by obviating unnecessary chemotherapy.

Purpose Digistain Index (DI), measured using an inexpensive mid-infrared spectrometer, reflects the level of aneuploidy in unstained tissue sections and correlates with tumor grade. We investigated whether incorporating DI with other clinicopathological variables could predict outcomes in patients with early breast cancer. Methods DI was calculated in 801 patients with hormone receptor-positive, HER2-negative primary breast cancer and ≤ 3 positive lymph nodes. All patients were treated with systemic endocrine therapy and no chemotherapy. Multivariable proportional hazards modeling was used to incorporate DI with clinicopathological variables to generate the Digistain Prognostic Score (DPS). DPS was assessed for prediction of 5- and 10-year outcomes (recurrence, recurrence-free survival [RFS] and overall survival [OS]) using receiver operating characteristics and Cox proportional hazards regression models. Kaplan–Meier analysis evaluated the ability of DPS to stratify risk. Results DPS was consistently highly accurate and had negative predictive values for all three outcomes, ranging from 0.96 to 0.99 at 5 years and 0.84 to 0.95 at 10 years. DPS demonstrated statistically significant prognostic ability with significant hazard ratios (95% CI) for low- versus high-risk classification for RFS, recurrence and OS (1.80 [CI 1.31–2.48], 1.83 [1.32–2.52] and 1.77 [1.28–2.43], respectively; all P < 0.001). Conclusion DPS showed high accuracy and predictive performance, was able to stratify patients into low or high-risk, and considering its cost and rapidity, has the potential to offer clinical utility.

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.

Supplementary Methods, Figures 1-4 from Breast Cancer Dormancy Can Be Maintained by Small Numbers of Micrometastases

Background & Aims: While normal human liver is thought to be generally quiescent, clonal hepatocyte expansions have been observed, though neither their cellular source nor their expansion dynamics have been determined. Knowing the hepatocyte cell of origin, and their subsequent dynamics and trajectory within the human liver will provide an important basis to understand disease associated dysregulation.Methods: Herein, we use in vivo lineage tracing and methylation sequence analysis to demonstrate normal human hepatocyte ancestry. We exploit next-generation mitochondrial sequencing to determine hepatocyte clonal expansion dynamics across spatially distinct areas of laser-captured, microdissected, clones, in tandem with computational modelling in morphologically normal human liver.Results: Hepatocyte clones and rare SOX9+ hepatocyte progenitors commonly associate with portal tracts and we present evidence that clones can lineage-trace with cholangiocytes, indicating the presence of a bipotential common ancestor at this niche. Within clones, we demonstrate methylation CpG sequence diversity patterns indicative of periportal not pericentral ancestral origins, indicating a portal to central vein expansion trajectory. Using spatial analysis of mitochondrial DNA variants by next-generation sequencing coupled with mathematical modelling and Bayesian inference across the portal-central axis, we demonstrate that patterns of mitochondrial DNA variants reveal large numbers of spatially restricted mutations in conjunction with limited numbers of clonal mutations.Conclusions: These datasets support the existence of a periportal progenitor niche and indicate that clonal patches exhibit punctuated but slow growth, then quiesce, likely due to acute environmental stimuli. These findings crucially contribute to our understanding of hepatocyte dynamics in the normal human liver.& COPY; 2023 The Author(s). Published by Elsevier B.V. on behalf of European Association for the Study of the Liver. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Supplementary Figure 1 from Ectopic Expression of P-Cadherin Correlates with Promoter Hypomethylation Early in Colorectal Carcinogenesis and Enhanced Intestinal Crypt Fission In vivo

BACKGROUND & AIMS: Barrett's esophagus (BE) is a risk factor for esophageal adenocarcinoma but our understanding of how it evolves is poorly understood. We investigated BE gland phenotype distribution, the clonal nature of phenotypic change, and how phenotypic diversity plays a role in progression. METHODS: Using immunohistochemistry and histology, we analyzed the distribution and the diversity of gland phenotype between and within biopsy specimens from patients with nondysplastic BE and those who had progressed to dysplasia or had developed postesophagectomy BE. Clonal relationships were determined by the presence of shared mutations between distinct gland types using laser capture microdissection sequencing of the mitochondrial genome. RESULTS: We identified 5 different gland phenotypes in a cohort of 51 nondysplastic patients where biopsy specimens were taken at the same anatomic site (1.0-2.0 cm superior to the gastroesophageal junction. Here, we observed the same number of glands with 1 and 2 phenotypes, but 3 phenotypes were rare. We showed a common ancestor between parietal cell containing, mature gastric (oxyntocardiac) and goblet cell-containing, intestinal (specialized) gland phenotypes. Similarly, we have shown a clonal relationship between cardiac-type glands and specialized and mature intestinal glands. Using the Shannon diversity index as a marker of gland diversity, we observed significantly increased phenotypic diversity in patients with BE adjacent to dysplasia and predysplasia compared to nondysplastic BE and postesophagectomy BE, suggesting that diversity develops over time. CONCLUSIONS: We showed that the range of BE phenotypes represents an evolutionary process and that changes in gland diversity may play a role in progression. Furthermore, we showed a common ancestry between gastric and intestinal-type glands in BE.

BACKGROUND:Lung squamous cell carcinoma (LUSC) accounts for a significant proportion of cancer deaths worldwide, and is preceded by the appearance of progressively disorganised pre-invasive lesions in the airway epithelium. Yet the biological mechanisms underlying progression of pre-invasive lesions into invasive LUSC are not fully understood. LRIG1 (leucine-rich repeats and immunoglobulin-like domains 1) is downregulated in pre-invasive airway lesions and invasive LUSC tumours and this correlates with decreased lung cancer patient survival.METHODS AND RESULTS:Using an Lrig1 knock-in reporter mouse and human airway epithelial cells collected at bronchoscopy, we show that during homeostasis LRIG1 is heterogeneously expressed in the airway epithelium. In basal airway epithelial cells, the suspected cell of origin of LUSC, LRIG1 identifies a subpopulation of progenitor cells with higher in vitro proliferative and self-renewal potential in both the mouse and human. Using the N-nitroso-tris-chloroethylurea (NTCU)-induced murine model of LUSC, we find that Lrig1 loss-of-function leads to abnormally high cell proliferation during the earliest stages of pre-invasive disease and to the formation of significantly larger invasive tumours, suggesting accelerated disease progression.CONCLUSION:Together, our findings identify LRIG1 as a marker of basal airway progenitor cells with high proliferative potential and as a regulator of pre-invasive lung cancer progression. This work highlights the clinical relevance of LRIG1 and the potential of the NTCU-induced LUSC model for functional assessment of candidate tumour suppressors and oncogenes.

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.

The majority of human liver research is disease-focused such that far less is known of cellular dynamics within normal human liver. We have leveraged cytochrome c oxidase deficiency as a marker of clonal hepatocyte populations in such tissues. We demonstrate these populations commonly associate with portal tracts and lineage-trace hepatocytes with cholangiocytes, indicating the presence of a bipotential common ancestor at this niche. We also observe rare periportal SOX9+ hepatocytes progenitor candidates in our human tissues. To understand clonal expansion dynamics, we measured methylation diversity and identified mtDNA variants by next-generation sequencing within spatially-defined clonal hepatocyte patches. We coupled our sequencing with mathematical modelling and Bayesian inference to compare spatial patterns of mtDNA variants under assumptions with or without faster expansion from a portal-associated niche. These datasets support the existence of a periportal progenitor niche and indicate that clonal patches slowly expand, perhaps due to acute environmental stimuli, then quiesce. These findings crucially contribute to our understanding of hepatocyte dynamics in normal human liver and provide a baseline for understanding how such dynamics may be modulated in diseased liver.

Aneuploidy, defined as the loss and gain of whole and part chromosomes, is a near-ubiquitous feature of cancer genomes, is prognostic, and likely an important determinant of cancer cell biology. In colorectal cancer (CRC), aneuploidy is found in virtually all tumours, including precursor adenomas. However, the temporal evolutionary dynamics that select for aneuploidy remain broadly uncharacterised. Here we perform genomic analysis of 755 samples from a total of 167 patients with colorectal-derived neoplastic lesions that cross-sectionally represent the distinct stages of tumour evolution, and longitudinally track individual tumours through metastasis and treatment. Precancer lesions (adenomas) exhibited low levels of aneuploidy but high intra-tumour heterogeneity, whereas cancers had high aneuploidy but low heterogeneity, indicating that progression is through a genetic bottleneck that suppresses diversity. Individual CRC glands from the same tumour have similar karyotypes, despite prior evidence of ongoing instability at the cell level. Pseudo-stable aneuploid genomes were observed in metastatic lesions sampled from liver and other organs, after chemo- or targeted therapies, and late recurrences detected many years after the diagnosis of a primary tumour. Modelling indicates that these data are consistent with the action of stabilising selection that ‘traps’ cancer cell genomes on a fitness peak defined by the specific pattern of aneuploidy. These data show that the initial progression of CRC requires the traversal of a rugged fitness landscape and subsequent genomic evolution, including metastatic dissemination and therapeutic resistance, is constrained by stabilising selection.

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.

Objective I BD confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multiregion sequencing. Design E xome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 patients with CA-CRC (n= 55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using single nucleotide polymorphism arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n= 9), low-grade dysplasia (LGD; n= 30), high-grade dysplasia (HGD; n= 13), mixed LGD/HGD (n= 7) and CA-CRC (n= 19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. Results 10/12 tumours were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNA) per Mb, 20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated ageingassociated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighbouring CACRC, indicating a precancer mutational field. CA-CRCs were often near triploid (40%) or near tetraploid (20%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated ' catastrophic' CNA increase. Conclusions E volutionary genomic analysis revealed precancer clones bearing extensive SNAs and CNAs, with progression to cancer involving a dramatic accrual of CNAs at HGD. Detection of the cancerised field is an encouraging prospect for surveillance, but punctuated evolution may limit the window for early detection.