Strand-resolved mutagenicity of DNA damage and repair.

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Identification of Clonal Hematopoiesis Driver Mutations through In Silico Saturation Mutagenesis.

Clonal hematopoiesis (CH) is a phenomenon of clonal expansion of hematopoietic stem cells driven by somatic mutations affecting certain genes. Recently, CH has been linked to the development of hematologic malignancies, cardiovascular diseases, and other conditions. Although the most frequently mutated CH driver genes have been identified, a systematic landscape of the mutations capable of initiating this phenomenon is still lacking. Here, we trained machine-learning models for 12 of the most recurrent CH genes to identify their driver mutations. These models outperform expert-curated rules based on prior knowledge of the function of these genes. Moreover, their application to identify CH driver mutations across almost half a million donors of the UK Biobank reproduces known associations between CH driver mutations and age, and the prevalence of several diseases and conditions. We thus propose that these models support the accurate identification of CH across healthy individuals.

Genomic deletions explain the generation of alternative BRAF isoforms conferring resistance to MAPK inhibitors in melanoma.

Resistance to MAPK inhibitors (MAPKi), the main cause of relapse in BRAF-mutant melanoma, is associated with the production of alternative BRAF mRNA isoforms (altBRAFs) in up to 30% of patients receiving BRAF inhibitor monotherapy. These altBRAFs have been described as being generated by alternative pre-mRNA splicing, and splicing modulation has been proposed as a therapeutic strategy to overcome resistance. In contrast, we report that altBRAFs are generated through genomic deletions. Using different in vitro models of altBRAF-mediated melanoma resistance, we demonstrate the production of altBRAFs exclusively from the BRAF V600E allele, correlating with corresponding genomic deletions. Genomic deletions are also detected in tumor samples from melanoma and breast cancer patients expressing altBRAFs. Along with the identification of altBRAFs in BRAF wild-type and in MAPKi-naive melanoma samples, our results represent a major shift in our understanding of mechanisms leading to the generation of BRAF transcripts variants associated with resistance in melanoma.

Origins of Second Malignancies in Children and Mutational Footprint of Chemotherapy in Normal Tissues.

Pediatric cancers are rare diseases, and children without known germline predisposing conditions who develop a second malignancy during developmental ages are extremely rare. We present four such clinical cases and, through whole-genome and error-correcting ultra-deep duplex sequencing of tumor and normal samples, we explored the origin of the second malignancy in four children, uncovering different routes of development. The exposure to cytotoxic therapies was linked to the emergence of a secondary acute myeloid leukemia. A common somatic mutation acquired early during embryonic development was the driver of two solid malignancies in another child. In two cases, the two tumors developed from completely independent clones diverging during embryogenesis. Importantly, we demonstrate that platinum-based therapies contributed at least one order of magnitude more mutations per day of exposure than aging to normal tissues in these children. SIGNIFICANCE: Using whole-genome and error-correcting ultra-deep duplex sequencing, we uncover different origins for second neoplasms in four children. We also uncover the presence of platinum-related mutations across 10 normal tissues of exposed individuals, highlighting the impact that the use of cytotoxic therapies may have on cancer survivors. See related commentary by Pacyna and Nangalia, p. 900. This article is featured in Selected Articles from This Issue, p. 897.

Hotspot propensity across mutational processes.

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Cancer Evolution: A Multifaceted Affair.

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery. SIGNIFICANCE: Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Mutational signature dynamics shaping the evolution of oesophageal adenocarcinoma.

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma.

Copy number footprints of platinum-based anticancer therapies.

Recently, distinct mutational footprints observed in metastatic tumors, secondary malignancies and normal human tissues have been demonstrated to be caused by the exposure to several chemotherapeutic drugs. These characteristic mutations originate from specific lesions caused by these chemicals to the DNA of exposed cells. However, it is unknown whether the exposure to these chemotherapies leads to a specific footprint of larger chromosomal aberrations. Here, we address this question exploiting whole genome sequencing data of metastatic tumors obtained from patients exposed to different chemotherapeutic drugs. As a result, we discovered a specific copy number footprint across tumors from patients previously exposed to platinum-based therapies. This footprint is characterized by a significant increase in the number of chromosomal fragments of copy number 1-4 and size smaller than 10 Mb in exposed tumors with respect to their unexposed counterparts (median 14-387% greater across tumor types). The number of chromosomal fragments characteristic of the platinum-associated CN footprint increases significantly with the activity of the well known platinum-related footprint of single nucleotide variants across exposed tumors.

Genomics improves risk stratification of adults with T-cell acute lymphoblastic leukemia enrolled in measurable residual disease-oriented trials.

Genetic information has been crucial to understand the pathogenesis of T-cell acute lymphoblastic leukemia (T-ALL) at diagnosis and at relapse, but still nowadays has a limited value in a clinical context. Few genetic markers are associated with the outcome of T-ALL patients, independently of measurable residual disease (MRD) status after therapy. In addition, the prognostic relevance of genetic features may be modulated by the specific treatment used. We analyzed the genetic profile of 145 T-ALL patients by targeted deep sequencing. Genomic information was integrated with the clinicalbiological and survival data of a subset of 116 adult patients enrolled in two consecutive MRD-oriented trials of the Spanish PETHEMA (Programa Español de Tratamientos en Hematología) group. Genetic analysis revealed a mutational profile defined by DNMT3A/ N/KRAS/ MSH2/ U2AF1 gene mutations that identified refractory/resistant patients. Mutations in the DMNT3A gene were also found in the non-leukemic cell fraction of patients with T-ALL, revealing a possible mutational-driven clonal hematopoiesis event to prime T-ALL in elderly. The prognostic impact of this adverse genetic profile was independent of MRD status on day +35 of induction therapy. The combined worse-outcome genetic signature and MRD on day +35 allowed risk stratification of T-ALL into standard or high-risk groups with significantly different 5- year overall survival (OS) of 52% (95% confidence interval: 37-67) and 17% (95% confidence interval: 1-33), respectively. These results confirm the relevance of the tumor genetic profile in predicting patient outcome in adult T-ALL and highlight the need for novel gene-targeted chemotherapeutic schedules to improve the OS of poor-prognosis T-ALL patients.

Identification of Clonal Hematopoiesis Driver Mutations through In Silico Saturation Mutagenesis.

Clonal hematopoiesis (CH) is a phenomenon of clonal expansion of hematopoietic stem cells driven by somatic mutations affecting certain genes. Recently, CH has been linked to the development of a number of hematologic malignancies, cardiovascular diseases and other conditions. Although the most frequently mutated CH driver genes have been identified, a systematic landscape of the mutations capable of initiating this phenomenon is still lacking. Here, we train high-quality machine-learning models for 12 of the most recurrent CH driver genes to identify their driver mutations. These models outperform an experimental base-editing approach and expert-curated rules based on prior knowledge of the function of these genes. Moreover, their application to identify CH driver mutations across almost half a million donors of the UK Biobank reproduces known associations between CH driver mutations and age, and the prevalence of several diseases and conditions. We thus propose that these models support the accurate identification of CH across healthy individuals.

High p16 expression and heterozygous RB1 loss are biomarkers for CDK4/6 inhibitor resistance in ER+ breast cancer.

CDK4/6 inhibitors combined with endocrine therapy have demonstrated higher antitumor activity than endocrine therapy alone for the treatment of advanced estrogen receptor-positive breast cancer. Some of these tumors are de novo resistant to CDK4/6 inhibitors and others develop acquired resistance. Here, we show that p16 overexpression is associated with reduced antitumor activity of CDK4/6 inhibitors in patient-derived xenografts (n = 37) and estrogen receptor-positive breast cancer cell lines, as well as reduced response of early and advanced breast cancer patients to CDK4/6 inhibitors (n = 89). We also identified heterozygous RB1 loss as biomarker of acquired resistance and poor clinical outcome. Combination of the CDK4/6 inhibitor ribociclib with the PI3K inhibitor alpelisib showed antitumor activity in estrogen receptor-positive non-basal-like breast cancer patient-derived xenografts, independently of PIK3CA, ESR1 or RB1 mutation, also in drug de-escalation experiments or omitting endocrine therapy. Our results offer insights into predicting primary/acquired resistance to CDK4/6 inhibitors and post-progression therapeutic strategies.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Discovering the drivers of clonal hematopoiesis.

Mutations in genes that confer a selective advantage to hematopoietic stem cells (HSCs) drive clonal hematopoiesis (CH). While some CH drivers have been identified, the compendium of all genes able to drive CH upon mutations in HSCs remains incomplete. Exploiting signals of positive selection in blood somatic mutations may be an effective way to identify CH driver genes, analogously to cancer. Using the tumor sample in blood/tumor pairs as reference, we identify blood somatic mutations across more than 12,000 donors from two large cancer genomics cohorts. The application of IntOGen, a driver discovery pipeline, to both cohorts, and more than 24,000 targeted sequenced samples yields a list of close to 70 genes with signals of positive selection in CH, available at http://www.intogen.org/ch . This approach recovers known CH genes, and discovers other candidates.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Combined presentation and immunogenicity analysis reveals a recurrent RAS.Q61K neoantigen in melanoma.

Neoantigens are now recognized drivers of the antitumor immune response. Recurrent neoantigens, shared among groups of patients, have thus become increasingly coveted therapeutic targets. Here, we report on the data-driven identification of a robustly presented, immunogenic neoantigen that is derived from the combination of HLA-A*01:01 and RAS.Q61K. Analysis of large patient cohorts indicated that this combination applies to 3% of patients with melanoma. Using HLA peptidomics, we were able to demonstrate robust endogenous presentation of the neoantigen in 10 tumor samples. We detected specific reactivity to the mutated peptide within tumor-infiltrating lymphocytes (TILs) from 2 unrelated patients, thus confirming its natural immunogenicity. We further investigated the neoantigen-specific clones and their T cell receptors (TCRs) via a combination of TCR sequencing, TCR overexpression, functional assays, and single-cell transcriptomics. Our analysis revealed a diverse repertoire of neoantigen-specific clones with both intra- and interpatient TCR similarities. Moreover, 1 dominant clone proved to cross-react with the highly prevalent RAS.Q61R variant. Transcriptome analysis revealed a high association of TCR clones with specific T cell phenotypes in response to cognate melanoma, with neoantigen-specific cells showing an activated and dysfunctional phenotype. Identification of recurrent neoantigens and their reactive TCRs can promote "off-the-shelf" precision immunotherapies, alleviating limitations of personalized treatments.

The evolution of hematopoietic cells under cancer therapy.

Chemotherapies may increase mutagenesis of healthy cells and change the selective pressures in tissues, thus influencing their evolution. However, their contributions to the mutation burden and clonal expansions of healthy somatic tissues are not clear. Here, exploiting the mutational footprint of some chemotherapies, we explore their influence on the evolution of hematopoietic cells. Cells of Acute Myeloid Leukemia (AML) secondary to treatment with platinum-based drugs show the mutational footprint of these drugs, indicating that non-malignant blood cells receive chemotherapy mutations. No trace of the 5-fluorouracil (5FU) mutational signature is found in AMLs secondary to exposure to 5FU, suggesting that cells establishing the leukemia could be quiescent during treatment. Using the platinum-based mutational signature as a barcode, we determine that the clonal expansion originating the secondary AMLs begins after the start of the cytotoxic treatment. Its absence in clonal hematopoiesis cases is consistent with the start of the clonal expansion predating the exposure to platinum-based drugs.

In silico saturation mutagenesis of cancer genes.

Despite the existence of good catalogues of cancer genes1,2, identifying the specific mutations of those genes that drive tumorigenesis across tumour types is still a largely unsolved problem. As a result, most mutations identified in cancer genes across tumours are of unknown significance to tumorigenesis3. We propose that the mutations observed in thousands of tumours-natural experiments testing their oncogenic potential replicated across individuals and tissues-can be exploited to solve this problem. From these mutations, features that describe the mechanism of tumorigenesis of each cancer gene and tissue may be computed and used to build machine learning models that encapsulate these mechanisms. Here we demonstrate the feasibility of this solution by building and validating 185 gene-tissue-specific machine learning models that outperform experimental saturation mutagenesis in the identification of  driver and passenger mutations. The models and their assessment of each mutation are designed to be interpretable, thus avoiding a black-box prediction device. Using these models, we outline the blueprints of potential driver mutations in cancer genes, and demonstrate the role of mutation probability in shaping the landscape of observed driver mutations. These blueprints will support the interpretation of newly sequenced tumours in patients and the study of the mechanisms of tumorigenesis of cancer genes across tissues.