BACKGROUND: Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.
METHODS: Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define \'shared\' (between both samples) and \'private\' (present in one sample) mutations.
RESULTS: Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.
CONCLUSIONS: This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

BACKGROUND: A reliable preclinical model of patient-derived organoids (PDOs) was developed in a case study of a 69-year-old woman diagnosed with breast cancer (BC) to investigate the tumour evolution before and after neoadjuvant chemotherapy and surgery. The results were achieved due to the development of PDOs from tissues collected before (O-PRE) and after (O-POST) treatment.
METHODS: PDO cultures were characterized by histology, immunohistochemistry (IHC), transmission electron microscopy (TEM), scanning electron microscopy (SEM), confocal microscopy, flow cytometry, real-time PCR, bulk RNA-seq, single-cell RNA sequencing (scRNA-seq) and drug screening.
RESULTS: Both PDO cultures recapitulated the histological and molecular profiles of the original tissues, and they showed typical mammary gland organization, confirming their reliability as a personalized in vitro model. Compared with O-PRE, O-POST had a greater proliferation rate with a significant increase in the Ki67 proliferation index. Moreover O-POST exhibited a more stem-like and aggressive phenotype, with increases in the CD24low/CD44low and EPCAMlow/CD49fhigh cell populations characterized by increased tumour initiation potential and multipotency and metastatic potential in invasive lobular carcinoma. Analysis of ErbB receptor expression indicated a decrease in HER-2 expression coupled with an increase in EGFR expression in O-POST. In this context, deregulation of the PI3K/Akt signalling pathway was assessed by transcriptomic analysis, confirming the altered transcriptional profile. Finally, transcriptomic single-cell analysis identified 11 cell type clusters, highlighting the selection of the luminal component and the decrease in the number of Epithelial-mesenchymal transition cell types in O-POST.
CONCLUSIONS: Neoadjuvant treatment contributed to the enrichment of cell populations with luminal phenotypes that were more resistant to chemotherapy in O-POST. PDOs represent an excellent 3D cell model for assessing disease evolution.

All tumour cells in a patient have shared and non-shared genetic alterations, and the diversity of mutations is described as intratumoural heterogeneity (ITH). Multiregion sequencing is a genome sequencing analytical technique used for multiple, spatially-separated biopsy tissues that may further our understanding of ITH and tumour evolution. Although genetic mutations in extramammary Paget\'s disease (EMPD) have recently been detected by next-generation sequencing analysis, there have been no reports of ITH based on multiregion sequencing in EMPD. Thus, we clarified the landscape of ITH and tumour evolution in EMPD. We performed whole-exome sequencing on 35 tissues (30 tumour tissues and five normal skin samples as a paired control), collected from five patients with EMPD. The rate of private mutations was significantly higher than that of ubiquitous and shared mutations. Ubiquitous mutations were not present in driver genes, and most driver genes exhibited private and shared mutations. The most frequent base substitution was C>T in almost all lesions, and most mutational signatures corresponded to signature 1, 2, 3, and 8. The types of proposed aetiology in most lesions were based on age and AID/APOBEC family and BRCA1/BRCA2 mutations. Evolutionary trees were characterized by short trunks and long branches due to the extremely high ratio of private mutations. In contrast, pathogenic factors, such as base substitutions, mutational signatures, and proposed aetiology, were shared. Tumour evolution in EMPD appears to be characterized by a high level of genetic ITH with shared background factors.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

Follicular lymphoma (FL) develops through a stepwise acquisition of cooperative genetic changes with t(14;18)(q32;q21)/IGH::BCL2 occurring early at the pre-B stage of B-cell development. Patients with FL typically show an indolent clinical course, remitting and relapsing with the eventual development of resistance to treatments. Interestingly, the majority of transformed FL do not progress directly from FL but originate from their clonally related lymphoma precursor (CLP) cells. To examine whether such divergent tumour evolution also underpins the relapses in patients with early-stage FL, we investigated by targeted next-generation sequencing 13 cases (stage I = 9, stage II = 4), who showed complete remission (mean: 5 years; range: 1-11.5 years) following local radiotherapy but subsequently relapsed (≥2 in 5). A clonal relationship between the diagnostic FL and relapses was confirmed in 11 cases. In six cases, common and distinct variants were seen between the paired diagnostic and relapsed lymphomas, indicating their divergent evolution from a CLP. In two cases, different B-cell clones were involved in the diagnostic and relapsed lymphomas, including one case involving two different BCL2 translocations. In the remaining five cases, the relapsed lymphoma developed via a linear progression (n = 4) or a mixed evolutionary path (n = 1). These findings may bear important implications in the routine diagnosis and management of relapsed FL. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Mutation accumulation in tumour evolution is one major cause of intra-tumour heterogeneity (ITH), which often leads to drug resistance during treatment. Previous studies with multi-region sequencing have shown that mutation divergence among samples within the patient is common, and the importance of spatial sampling to obtain a complete picture in tumour measurements. However, quantitative comparisons of the relationship between mutation heterogeneity and tumour expansion modes, sampling distances as well as the sampling methods are still few. Here, we investigate how mutations diverge over space by varying the sampling distance and tumour expansion modes using individual-based simulations. We measure ITH by the Jaccard index between samples and quantify how ITH increases with sampling distance, the pattern of which holds in various sampling methods and sizes. We also compare the inferred mutation rates based on the distributions of variant allele frequencies under different tumour expansion modes and sampling sizes. In exponentially fast expanding tumours, a mutation rate can always be inferred for any sampling size. However, the accuracy compared with the true value decreases when the sampling size decreases, where small sampling sizes result in a high estimate of the mutation rate. In addition, such an inference becomes unreliable when the tumour expansion is slow, such as in surface growth.

Aneuploidy is defined as the cellular state of having a number of chromosomes that deviates from a multiple of the normal haploid chromosome number of a given organism. Aneuploidy can be present in a static state: Down syndrome individuals stably maintain an extra copy of chromosome 21 in their cells. In cancer cells, however, aneuploidy is usually present in combination with chromosomal instability (CIN) which leads to a continual generation of new chromosomal alterations and the development of intratumour heterogeneity (ITH). The prevalence of cells with specific chromosomal alterations is further shaped by evolutionary selection, for example, during the administration of cancer therapies. Aneuploidy, CIN and ITH have each been individually associated with poor prognosis in cancer, and a wealth of evidence suggests they contribute, either alone or in combination, to cancer therapy resistance by providing a reservoir of potential resistant states, or the ability to rapidly evolve resistance. A full understanding of the contribution and interplay between aneuploidy, CIN and ITH is required to tackle therapy resistance in cancer patients. However, these characteristics often co-occur and are intrinsically linked, presenting a major challenge to defining their individual contributions. Moreover, their accurate measurement in both experimental and clinical settings is a technical hurdle. Here, we attempt to deconstruct the contribution of the individual and combined roles of aneuploidy, CIN and ITH to therapy resistance in cancer, and outline emerging approaches to measure and disentangle their roles as a step towards integrating these principles into cancer therapeutic strategy.

Liquid biopsies and the dynamic tracking of somatic mutations within circulating tumour DNA (ctDNA) can provide insight into the dynamics of cancer evolution and the intra-tumour heterogeneity that fuels treatment resistance. However, identifying and tracking dynamic changes in somatic copy number alterations (SCNAs), which have been associated with poor outcome and metastasis, using ctDNA is challenging. Pancreatic adenocarcinoma is a disease which has been considered to harbour early punctuated events in its evolution, leading to an early fitness peak, with minimal further subclonal evolution.
To interrogate the role of SCNAs in pancreatic adenocarcinoma cancer evolution, we applied whole-exome sequencing of 55 longitudinal cell-free DNA (cfDNA) samples taken from 24 patients (including 8 from whom a patient-derived xenograft (PDX) was derived) with metastatic disease prospectively recruited into a clinical trial. We developed a method, Aneuploidy in Circulating Tumour DNA (ACT-Discover), that leverages haplotype phasing of paired tumour biopsies or PDXs to identify SCNAs in cfDNA with greater sensitivity.
SCNAs were observed within 28 of 47 evaluable cfDNA samples. Of these events, 30% could only be identified by harnessing the haplotype-aware approach leveraged in ACT-Discover. The exceptional purity of PDX tumours enabled near-complete phasing of genomic regions in allelic imbalance, highlighting an important auxiliary function of PDXs. Finally, although the classical model of pancreatic cancer evolution emphasises the importance of early, homogenous somatic events as a key requirement for cancer development, ACT-Discover identified substantial heterogeneity of SCNAs, including parallel focal and arm-level events, affecting different parental alleles within individual tumours. Indeed, ongoing acquisition of SCNAs was identified within tumours throughout the disease course, including within an untreated metastatic tumour.
This work demonstrates the power of haplotype phasing to study genomic variation in cfDNA samples and reveals undiscovered intra-tumour heterogeneity with important scientific and clinical implications. Implementation of ACT-Discover could lead to important insights from existing cohorts or underpin future prospective studies seeking to characterise the landscape of tumour evolution through liquid biopsy.

Metastasis is a complex process and the leading cause of cancer-related death globally. Recent studies have demonstrated that genomic sequencing data from paired primary and metastatic tumours can be used to trace the evolutionary origins of cells responsible for metastasis. This approach has yielded new insights into the genomic alterations that engender metastatic potential, and the mechanisms by which cancer spreads. Given that the reliability of these approaches is contingent upon how representative the samples are of primary and metastatic tumour heterogeneity, we review insights from studies that have reconstructed the evolution of metastasis within the context of their cohorts and designs. We discuss the role of research autopsies in achieving the comprehensive sampling necessary to advance the current understanding of metastasis.

Cancer development is an evolutionary process. A key selection pressure is exerted by therapy, one of the few players in cancer evolution that can be controlled. As such, an understanding of how treatment acts to sculpt the tumour and its microenvironment and how this influences a tumour\'s subsequent evolutionary trajectory is critical. In this review, we examine cancer evolution and intra-tumour heterogeneity in the context of therapy. We focus on how radiotherapy, chemotherapy and immunotherapy shape both tumour development and the environment in which tumours evolve and how resistance can develop or be selected for during treatment.