Genetic immune escape landscape in primary and metastatic cancer

Quantification and Functional Studies of Neutrophilic Cells Identifies Distinct Papilloma Phenotypes

OBJECTIVES

To characterize the distribution of immune cell subsets within laryngeal papillomas and to study the function of potentially immunosuppressive neutrophilic and regulatory T cells (Tregs).

METHODS

Fresh clinical papilloma specimens were collected at the time of surgery and studied with multiparameter flow cytometry. Papilloma infiltrating neutrophilic cells and Tregs were sorted and studied functionally with ex vivo T cell suppression assays.

RESULTS

Flow cytometric analysis of fresh laryngeal papillomas samples from 18 adult patients with recurrent respiratory papillomatosis revealed patterns in immune constituency between patients. Clearly divergent phenotypes based primarily on the degree of neutrophilic and T cell infiltration were identified. Relative neutrophilic cell enrichment and T cell depletion were observed in 50% of samples and neutrophilic cell depletion and T cell enrichment were observed in the others. Greater papilloma neutrophilic cell enrichment was positively associated with the number of clinically indicated interventions required in the 12 months prior to sample collection, linking papilloma neutrophil inflammation to disease severity. Functional assays revealed the ability of both papilloma infiltrating neutrophilic and Tregs to suppress T cell function at roughly equal magnitudes, but substantially increased infiltration of neutrophilic cells compared to Tregs across samples.

CONCLUSION

Neutrophilic cells are an important contributor to immunosuppression within the respiratory papilloma microenvironment. Given these data and the association between greater neutrophilic cell infiltration and lack of clinical response to therapeutic vaccination, additional study of strategies aimed at limiting neutrophilic cell infiltration or function within papillomas is warranted.

LEVEL OF EVIDENCE

4 Laryngoscope, 134:3238-3244, 2024.

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

Epigenetic regulation of breast cancer metastasis

Breast cancer is the most frequently diagnosed malignancy and the second leading cause of cancer-related mortality among women worldwide. Recurrent metastasis is associated with poor patient outcomes and poses a significant challenge in breast cancer therapies. Cancer cells adapting to a new tissue microenvironment is the key event in distant metastasis development, where the disseminating tumor cells are likely to acquire genetic and epigenetic alterations during the process of metastatic colonization. Despite several decades of research in this field, the exact mechanisms governing metastasis are not fully understood. However, emerging body of evidence indicates that in addition to genetic changes, epigenetic reprogramming of cancer cells and the metastatic niche are paramount toward successful metastasis. Here, we review and discuss the latest knowledge about the salient attributes of metastasis and epigenetic regulation in breast cancer and crucial research domains that need further investigation.

Germline-mediated immunoediting sculpts breast cancer subtypes and metastatic proclivity

Tumors with the same diagnosis can have different molecular profiles and response to treatment. It remains unclear when and why these differences arise. Somatic genomic aberrations occur within the context of a highly variable germline genome. Interrogating 5870 breast cancer lesions, we demonstrated that germline-derived epitopes in recurrently amplified genes influence somatic evolution by mediating immunoediting. Individuals with a high germline-epitope burden in human epidermal growth factor receptor 2 (HER2/ERBB2) are less likely to develop HER2-positive breast cancer compared with other subtypes. The same holds true for recurrent amplicons defining three aggressive estrogen receptor (ER)-positive subgroups. Tumors that overcome such immune-mediated negative selection are more aggressive and demonstrate an "immune cold" phenotype. These data show that the germline genome plays a role in dictating somatic evolution.

Genetic landscape of interval and screen detected breast cancer

Interval breast cancers (IBCs) are cancers diagnosed between screening episodes. Understanding the biological differences between IBCs and screen-detected breast-cancers (SDBCs) has the potential to improve mammographic screening and patient management. We analysed and compared the genomic landscape of 288 IBCs and 473 SDBCs by whole genome sequencing of paired tumour-normal patient samples collected as part of the UK 100,000 Genomes Project. Compared to SDBCs, IBCs were more likely to be lobular, higher grade, and triple negative. A more aggressive clinical phenotype was reflected in IBCs displaying features of genomic instability including a higher mutation rate and number of chromosomal structural abnormalities, defective homologous recombination and TP53 mutations. We did not however, find evidence to indicate that IBCs are associated with a significantly different immune response. While IBCs do not represent a unique molecular class of invasive breast cancer they exhibit a more aggressive phenotype, which is likely to be a consequence of the timing of tumour initiation. This information is relevant both with respect to treatment as well as informing the screening interval for mammography.

Molecular characterization of the evolution of premalignant lesions in the upper aerodigestive tract

Introduction

Early relapse and development of metastatic disease are some of the primary reasons for the poor prognosis of patients with head and neck squamous cell carcinoma (HNSCC). HNSCC is a heterogeneous disease which may develop in large premalignant fields of genetically altered cells. Yet knowing which individuals will progress and develop clinically significant cancers during their lifetimes remains one of the most important challenges of reducing HNSCC morbidity and mortality. To further elucidate the molecular mechanisms, we performed a focused analysis of the genome and immune microenvironment from multiple, matched normal squamous tissue, premalignant lesions, as well as primary and recurrent tumors from seven patients with p16-negative HNSCC.

Methods

We performed targeted panel Next Generation Sequencing (161 genes) to analyze somatic variants from sequentially collected, matched formalin-fixed paraffin-embedded tissue (normal, premalignant, HNSCC) from two patients. These samples plus samples from five additional patients were analyzed with the Nanostring PanCancer Immune Panel. In addition, we performed shallow whole genome sequencing (0.5x coverage on average) on samples from three of these patients. Patients were, apart from one case, primarily treated with curative-intent surgery, and received subsequent adjuvant treatment, if indicated.

Results

The most frequently mutated genes were TP53 and NOTCH1. Other mutated genes included NOTCH3 and CDKN2A, among others. A significant number of mutations were private to dysplasia and invasive carcinoma, respectively, however, almost 20% were shared between them. Increasing genomic instability was observed when comparing histologically normal squamous mucosa with higher levels of dysplasia. High-grade dysplasia showed similarly rearranged genomes as invasive carcinoma. Pathways related to interferon alpha and gamma response were upregulated even in moderate dysplastic lesions with increasing expression in higher grades of dysplasia and carcinoma. SPINK5, a known tumor suppressor gene in HNSCC, was already downregulated in low-grade dysplastic lesions, indicating an early deactivation in the evolution of the disease.

Conclusion

Genomic alterations as well as aberrant immune gene expression can be observed early in the evolution of tumors of the upper aerodigestive tract, highlighting the potential for targeting early mechanisms of disease progression.

Insights into the mechanisms and structure of breakage-fusion-bridge cycles in cervical cancer using long-read sequencing

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and our study applies long-read sequencing to this important chromosomal rearrangement type. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes, and only one BFB breakpoint showed chromothripsis. Five cell lines have a chromosomal region 11q BFB event, with YAP1-BIRC3-BIRC2 amplification. Indeed, YAP1 amplification is associated with a 10-year-earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that individuals with cervical cancer and YAP1-BIRC3-BIRC2 amplification, especially those of African ancestry, might benefit from targeted therapy. In summary, we uncovered valuable insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

Optimized whole-genome sequencing workflow for tumor diagnostics in routine pathology practice

Two decades after the genomics revolution, oncology is rapidly transforming into a genome-driven discipline, yet routine cancer diagnostics is still mainly microscopy based, except for tumor type-specific predictive molecular tests. Pathology laboratories struggle to quickly validate and adopt biomarkers identified by genomics studies of new targeted therapies. Consequently, clinical implementation of newly approved biomarkers suffers substantial delays, leading to unequal patient access to these therapies. Whole-genome sequencing (WGS) can successfully address these challenges by providing a stable molecular diagnostic platform that allows detection of a multitude of genomic alterations in a single cost-efficient assay and facilitating rapid implementation, as well as by the development of new genomic biomarkers. Recently, the Whole-genome sequencing Implementation in standard Diagnostics for Every cancer patient (WIDE) study demonstrated that WGS is a feasible and clinically valid technique in routine clinical practice with a turnaround time of 11 workdays. As a result, WGS was successfully implemented at the Netherlands Cancer Institute as part of routine diagnostics in January 2021. The success of implementing WGS has relied on adhering to a comprehensive protocol including recording patient information, sample collection, shipment and storage logistics, sequencing data interpretation and reporting, integration into clinical decision-making and data usage. This protocol describes the use of fresh-frozen samples that are necessary for WGS but can be challenging to implement in pathology laboratories accustomed to using formalin-fixed paraffin-embedded samples. In addition, the protocol outlines key considerations to guide uptake of WGS in routine clinical care in hospitals worldwide.

An immunogenetic basis for lung cancer risk

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

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.

Tumor-localized interleukin-2 and interleukin-12 combine with radiation therapy to safely potentiate regression of advanced malignant melanoma in pet dogs

The clinical use of interleukin-2 and -12 cytokines against cancer is limited by their narrow therapeutic windows due to on-target, off-tumor activation of immune cells when delivered systemically. Engineering IL-2 and IL-12 to bind to extracellular matrix collagen allows these cytokines to be retained within tumors after intralesional injection, overcoming these clinical safety challenges. While this approach has potentiated responses in syngeneic mouse tumors without toxicity, the complex tumor-immune interactions in human cancers are difficult to recapitulate in mouse models of cancer. This has driven an increased role for comparative oncology clinical trials in companion (pet) dogs with spontaneous cancers that feature analogous tumor and immune biology to human cancers. Here, we report the results from a dose-escalation clinical trial of intratumoral collagen-binding IL-2 and IL-12 cytokines in pet dogs with malignant melanoma, observing encouraging local and regional responses to therapy that may suggest human clinical benefit with this approach.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Navigating the Immune Maze: Pioneering Strategies for Unshackling Cancer Immunotherapy Resistance

Cancer immunotherapy has ushered in a transformative era in oncology, offering unprecedented promise and opportunities. Despite its remarkable breakthroughs, the field continues to grapple with the persistent challenge of treatment resistance. This resistance not only undermines the widespread efficacy of these pioneering treatments, but also underscores the pressing need for further research. Our exploration into the intricate realm of cancer immunotherapy resistance reveals various mechanisms at play, from primary and secondary resistance to the significant impact of genetic and epigenetic factors, as well as the crucial role of the tumor microenvironment (TME). Furthermore, we stress the importance of devising innovative strategies to counteract this resistance, such as employing combination therapies, tailoring immune checkpoints, and implementing real-time monitoring. By championing these state-of-the-art methods, we anticipate a paradigm that blends personalized healthcare with improved treatment options and is firmly committed to patient welfare. Through a comprehensive and multifaceted approach, we strive to tackle the challenges of resistance, aspiring to elevate cancer immunotherapy as a beacon of hope for patients around the world.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing the most detailed somatic mutational landscape to date. We identify new driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for drug repurposing. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The twin observations that higher T-cell infiltration is associated with better outcome and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

Advances in tumor immunomodulation based on nanodrug delivery systems

Immunotherapy is a therapeutic approach that employs immunological principles and techniques to enhance and amplify the body's immune response, thereby eradicating tumor cells. Immunotherapy has demonstrated effective antitumor effects on a variety of malignant tumors. However, when applied to humans, many immunotherapy drugs fail to target lesions with precision, leading to an array of adverse immune-related reactions that profoundly limit the clinical application of immunotherapy. Nanodrug delivery systems enable the precise delivery of immunotherapeutic drugs to targeted tissues or specific immune cells, enhancing the immune antitumor effect while reducing the number of adverse reactions. A nanodrug delivery system provides a feasible strategy for activating the antitumor immune response by the following mechanisms: 1) increased targeting and uptake of vaccines by DCs, which enhances the efficacy of the immune response; 2) increased tumor cell immunogenicity; 3) regulation of TAMs and other cells by, for example, regulating the polarization of TAMs and interfering with TAN formation, and ECM remodeling by CAFs; and 4) interference with tumor immune escape signaling pathways, namely, the PD-1/PD-L1, FGL1/LAG-3 and IDO signaling pathways. This paper reviews the progress of nanodrug delivery system research with respect to tumor immunotherapy based on tumor immunomodulation over the last few years, discussing the promising future of these delivery systems under this domain.

Immunoediting of KEAP1-NRF2 mutant tumours is required to circumvent NRF2-mediated immune surveillance

In human cancer, activating mutations in the KEAP1-NRF2 pathway are frequently observed, and positively selected for, as they confer the cytoprotective functions of the transcription factor NRF2 on the cancer cells. This results in the development of aggressive tumours which are resistant to treatment with chemotherapeutic compounds. Recent clinical developments have also revealed that NRF2-activated cancers are similarly resistant to immune checkpoint inhibitor drugs. As the mechanism of action of these immune modulating therapies is tangential to the classical cytoprotective function of NRF2, it is unclear how aberrant NRF2 activity could impact the anti-cancer functionality of the immune system. In this context, we found that in human cancer, NRF2-activated cells are highly immunoedited, which allows the cancer cells to escape immune surveillance and develop into malignant tumours. This immunoediting takes the form of reduced antigen presentation by the MHC-I complex, coupled with reduced expression of activating ligands for NK cells. Together, these modifications to the immunogenicity of NRF2-activated cancers inhibit immune effector cell infiltration and engagement, and contribute to the formation of the immunologically cold tumour microenvironment which is a characteristic feature of NRF2-activated cancers.

Association of HLA-A, B, and C alleles and cancer susceptibility in 179 solid malignancies

BACKGROUND

The major histocompatibility complex (MHC) genes are known to be capable of influencing the susceptibility of many cancers. All mammalian cells, including cancer cells, express MHC class I molecules consisting of human leukocyte antigens (HLA) A, B, and C. The tumor susceptibility of HLA-A, B, and C alleles has not been studied extensively in solid tumors.

METHODS

HLA-A, B, and C genotypes of 179 solid tumors were collected from Caris Comprehensive Tumor Profiling reports, including 45 GU, 44 GI, 28 pancreaticobiliary, 21 thoracic, 15 breast, 13 Gyn, among others. The tumors were mainly from Caucasians (82%). The HLA allele frequencies in the tumors were compared to those of respective ethnic populations in the US National Marrow Donor Program (NMDP) database. Fisher's exact tests were performed, adjusted P values were calculated using Benjamini-Hochberg's method for false discovery rate (FDR), and Prevalence ratios (PRs) were calculated to quantify associations.

RESULTS

Twenty-one alleles were not listed in the NMDP. Among them, A*11:303 alone was present in 11 carcinomas, and B*08:222 was seen in 4 tumors. Among the alleles listed in the NMDP, C*08:02, B*14:02, A*03:02, and B*44:06 were significantly associated with tumors in Caucasian Americans (PR: 2.50-170), while B*44:02 appeared protective (PR: 0.36). Alleles with less significant associations were listed.

CONCLUSIONS

From the HLA-A, B, and C data of the 179 tumors, we identified several susceptible alleles and one protective allele. Of interest, 21 alleles were not listed in the NMDP. The limited cases prevented our analysis from identifying cancer-susceptible alleles in other races.

Insights into the Mechanisms and Structure of Breakage-Fusion-Bridge Cycles in Cervical Cancer using Long-Read Sequencing

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and this is one of the first analyses of these events using long-read sequencing. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes and only one BFB breakpoint showed chromothripsis. Five cell lines have a Chr11q BFB event, with YAP1/BIRC2/BIRC3 gene amplification. Indeed, YAP1 amplification is associated with a 10-year earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that cervical cancer patients with YAP1/BIRC2/BIRC3-amplification, especially those of African American ancestry, might benefit from targeted therapy. In summary, we uncovered new insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

Pan-cancer whole-genome comparison of primary and metastatic solid tumours

Metastatic cancer remains an almost inevitably lethal disease1-3. A better understanding of disease progression and response to therapies therefore remains of utmost importance. Here we characterize the genomic differences between early-stage untreated primary tumours and late-stage treated metastatic tumours using a harmonized pan-cancer analysis (or reanalysis) of two unpaired primary4 and metastatic5 cohorts of 7,108 whole-genome-sequenced tumours. Metastatic tumours in general have a lower intratumour heterogeneity and a conserved karyotype, displaying only a modest increase in mutations, although frequencies of structural variants are elevated overall. Furthermore, highly variable tumour-specific contributions of mutational footprints of endogenous (for example, SBS1 and APOBEC) and exogenous mutational processes (for example, platinum treatment) are present. The majority of cancer types had either moderate genomic differences (for example, lung adenocarcinoma) or highly consistent genomic portraits (for example, ovarian serous carcinoma) when comparing early-stage and late-stage disease. Breast, prostate, thyroid and kidney renal clear cell carcinomas and pancreatic neuroendocrine tumours are clear exceptions to the rule, displaying an extensive transformation of their genomic landscape in advanced stages. Exposure to treatment further scars the tumour genome and introduces an evolutionary bottleneck that selects for known therapy-resistant drivers in approximately half of treated patients. Our data showcase the potential of pan-cancer whole-genome analysis to identify distinctive features of late-stage tumours and provide a valuable resource to further investigate the biological basis of cancer and resistance to therapies.