Papillary renal cell carcinoma (pRCC) is an important subtype of kidney cancer with a problematic pathological classification and highly variable clinical behaviour. Here we sequence the genomes or exomes of 31 pRCCs, and in four tumours, multi-region sequencing is undertaken. We identify BAP1, SETD2, ARID2 and Nrf2 pathway genes (KEAP1, NHE2L2 and CUL3) as probable drivers, together with at least eight other possible drivers. However, only ~10% of tumours harbour detectable pathogenic changes in any one driver gene, and where present, the mutations are often predicted to be present within cancer sub-clones. We specifically detect parallel evolution of multiple SETD2 mutations within different sub-regions of the same tumour. By contrast, large copy number gains of chromosomes 7, 12, 16 and 17 are usually early, monoclonal changes in pRCC evolution. The predominance of large copy number variants as the major drivers for pRCC highlights an unusual mode of tumorigenesis that may challenge precision medicine approaches.

Genomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution.

Mitotic recombination is important for inactivating tumour suppressor genes by copy-neutral loss of heterozygosity (LOH). Although meiotic recombination maps are plentiful, little is known about mitotic recombination. The APC gene (chr5q21) is mutated in most colorectal tumours and its usual mode of LOH is mitotic recombination.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. However, evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease.

Checkpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

In order to identify new genes with differential expression in early intestinal tumours, we performed mRNA (messenger ribonucleic acid) expression profiling of 16 human and 63 mouse adenomas. All individuals had germline APC mutations to ensure that tumorigenesis was driven by 'second hits' at APC. Using stringent filtering to identify changes consistent between humans and mice, we identified 60 genes up-regulated and 151 down-regulated in tumours. For 22 selected genes--including known Wnt targets--expression differences were confirmed by qRT-PCR (quantitative reverse transcription polymerase chain reaction). Most, but not all, differences were also present in colorectal carcinomas. In situ analysis showed a complex picture. Expression of up-regulated genes in adenomas was usually uniform/diffuse (e.g. ITGA6) or prominent in the tumour core (e.g. LGR5); in normal tissue, these genes were expressed at crypt bases or the transit amplifying zone. Down-regulated genes were often undetectable in adenomas, but in normal tissue were expressed in mesenchyme (e.g. GREM1/2) or differentiated cells towards crypt tops (e.g. SGK1). In silico analysis of TCF4-binding motifs showed that some of our genes were probably direct Wnt targets. Previous studies, mostly focused on human tumours, showed partial overlap with our 'expression signature', but 37 genes were unique to our study, including TACSTD2, SEMA3F, HOXA9 and IER3 (up-regulated), and TAGLN, GREM1, GREM2, MAB21L2 and RARRES2 (down-regulated). Combined analysis of our and published human data identified additional genes differentially expressed in adenomas, including decreased BMPs (bone morphogenetic proteins) and increased BUB1/BUB1B. Several of the newly identified, differentially expressed genes represent potential diagnostic or therapeutic targets for intestinal tumours.

Although thousands of solid tumors have been sequenced to date, a fundamental under-sampling bias is inherent in current methodologies. This is caused by a tissue sample input of fixed dimensions (e.g., 6 mm biopsy), which becomes grossly under-powered as tumor volume scales. Here, we demonstrate representative sequencing (Rep-Seq) as a new method to achieve unbiased tumor tissue sampling. Rep-Seq uses fixed residual tumor material, which is homogenized and subjected to next-generation sequencing. Analysis of intratumor tumor mutation burden (TMB) variability shows a high level of misclassification using current single-biopsy methods, with 20% of lung and 52% of bladder tumors having at least one biopsy with high TMB but low clonal TMB overall. Misclassification rates by contrast are reduced to 2% (lung) and 4% (bladder) when a more representative sampling methodology is used. Rep-Seq offers an improved sampling protocol for tumor profiling, with significant potential for improved clinical utility and more accurate deconvolution of clonal structure.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Genetic intra-tumour heterogeneity fuels clonal evolution, but our understanding of clinically relevant clonal dynamics remain limited. We investigated spatial and temporal features of clonal diversification in clear cell renal cell carcinoma through a combination of modelling and real tumour analysis. We observe that the mode of tumour growth, surface or volume, impacts the extent of subclonal diversification, enabling interpretation of clonal diversity in patient tumours. Specific patterns of proliferation and necrosis explain clonal expansion and emergence of parallel evolution and microdiversity in tumours. In silico time-course studies reveal the appearance of budding structures before detectable subclonal diversification. Intriguingly, we observe radiological evidence of budding structures in early-stage clear cell renal cell carcinoma, indicating that future clonal evolution may be predictable from imaging. Our findings offer a window into the temporal and spatial features of clinically relevant clonal evolution.

ADAPTeR is a prospective, phase II study of nivolumab (anti-PD-1) in 15 treatment-naive patients (115 multiregion tumor samples) with metastatic clear cell renal cell carcinoma (ccRCC) aiming to understand the mechanism underpinning therapeutic response. Genomic analyses show no correlation between tumor molecular features and response, whereas ccRCC-specific human endogenous retrovirus expression indirectly correlates with clinical response. T cell receptor (TCR) analysis reveals a significantly higher number of expanded TCR clones pre-treatment in responders suggesting pre-existing immunity. Maintenance of highly similar clusters of TCRs post-treatment predict response, suggesting ongoing antigen engagement and survival of families of T cells likely recognizing the same antigens. In responders, nivolumab-bound CD8+ T cells are expanded and express GZMK/B. Our data suggest nivolumab drives both maintenance and replacement of previously expanded T cell clones, but only maintenance correlates with response. We hypothesize that maintenance and boosting of a pre-existing response is a key element of anti-PD-1 mode of action.

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

A complete understanding of how exposure to environmental substances promotes cancer formation is lacking. More than 70 years ago, tumorigenesis was proposed to occur in a two-step process: an initiating step that induces mutations in healthy cells, followed by a promoter step that triggers cancer development1. Here we propose that environmental particulate matter measuring ≤2.5 μm (PM2.5), known to be associated with lung cancer risk, promotes lung cancer by acting on cells that harbour pre-existing oncogenic mutations in healthy lung tissue. Focusing on EGFR-driven lung cancer, which is more common in never-smokers or light smokers, we found a significant association between PM2.5 levels and the incidence of lung cancer for 32,957 EGFR-driven lung cancer cases in four within-country cohorts. Functional mouse models revealed that air pollutants cause an influx of macrophages into the lung and release of interleukin-1β. This process results in a progenitor-like cell state within EGFR mutant lung alveolar type II epithelial cells that fuels tumorigenesis. Ultradeep mutational profiling of histologically normal lung tissue from 295 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 53% of healthy tissue samples, respectively. These findings collectively support a tumour-promoting role for  PM2.5 air pollutants  and provide impetus for public health policy initiatives to address air pollution to reduce disease burden.

Understanding the evolutionary pathways to metastasis and resistance to immune-checkpoint inhibitors (ICI) in melanoma is critical for improving outcomes. Here, we present the most comprehensive intrapatient metastatic melanoma dataset assembled to date as part of the Posthumous Evaluation of Advanced Cancer Environment (PEACE) research autopsy program, including 222 exome sequencing, 493 panel-sequenced, 161 RNA sequencing, and 22 single-cell whole-genome sequencing samples from 14 ICI-treated patients. We observed frequent whole-genome doubling and widespread loss of heterozygosity, often involving antigen-presentation machinery. We found KIT extrachromosomal DNA may have contributed to the lack of response to KIT inhibitors of a KIT-driven melanoma. At the lesion-level, MYC amplifications were enriched in ICI nonresponders. Single-cell sequencing revealed polyclonal seeding of metastases originating from clones with different ploidy in one patient. Finally, we observed that brain metastases that diverged early in molecular evolution emerge late in disease. Overall, our study illustrates the diverse evolutionary landscape of advanced melanoma.

Adenomatous polyposis coli (APC) mutation is the hallmark of colorectal cancer (CRC), resulting in constitutive WNT activation. Despite decades of research, targeting WNT signaling in cancer remains challenging due to its on-target toxicity. We have previously shown that the deubiquitinating enzyme USP7 is a tumor-specific WNT activator in APC-truncated cells by deubiquitinating and stabilizing β-catenin, but its role in gut tumorigenesis is unknown. Here, we show in vivo that deletion of Usp7 in Apc-truncated mice inhibits crypt hyperproliferation and intestinal tumor development. Loss of Usp7 prolongs the survival of the sporadic intestinal tumor model. Genetic deletion, but not pharmacological inhibition, of Usp7 in Apc+/- intestine induces colitis and enteritis. USP7 inhibitor treatment suppresses growth of patient-derived cancer organoids carrying APC truncations in vitro and in xenografts. Our findings provide direct evidence that USP7 inhibition may offer a safe and efficacious tumor-specific therapy for both sporadic and germline APC-mutated CRC.

Intra-tumor heterogeneity concerns the existence of genetically different subclones within the same tumor. Single sample quantification of heterogeneity relies on precise determination of chromosomal copy numbers throughout the genome, and an assessment of whether identified mutation variant allele fractions match clonal or subclonal copy numbers. We discuss these issues using data from SNP arrays, whole exome sequencing and pathologist purity estimates on several breast cancers characterized by ERBB2 amplification. We show that chromosomal copy numbers can only be estimated from SNP array signals or sequencing depths for subclonal tumor samples with simple subclonal architectures under certain assumptions.

Dogs are ubiquitous in human society and attempts to manage their populations are common to most countries. Managing dog populations is achieved through a range of interventions to suit the dog population dynamics and dog ownership characteristics of the location, with a number of potential impacts or goals in mind. Impact assessment provides the opportunity for interventions to identify areas of inefficiencies for improvement and build evidence of positive change.

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Frameshift insertion/deletions (fs-indels) are an infrequent but highly immunogenic mutation subtype. Although fs-indels are degraded through the nonsense-mediated decay (NMD) pathway, we hypothesise that some fs-indels escape degradation and elicit anti-tumor immune responses. Using allele-specific expression analysis, expressed fs-indels are enriched in genomic positions predicted to escape NMD, and associated with higher protein expression, consistent with degradation escape (NMD-escape). Across four independent melanoma cohorts, NMD-escape mutations are significantly associated with clinical-benefit to checkpoint inhibitor (CPI) therapy (Pmeta = 0.0039). NMD-escape mutations are additionally found to associate with clinical-benefit in the low-TMB setting. Furthermore, in an adoptive cell therapy treated melanoma cohort, NMD-escape mutation count is the most significant biomarker associated with clinical-benefit. Analysis of functional T cell reactivity screens from personalized vaccine studies shows direct evidence of fs-indel derived neoantigens eliciting immune response, particularly those with highly elongated neo open reading frames. NMD-escape fs-indels represent an attractive target for biomarker optimisation and immunotherapy design.

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.