Substitution mutational signatures in whole-genome-sequenced cancers in the UK population

Increased Synthetic Cytotoxicity of Combinatorial Chemoradiation Therapy in Homologous Recombination Deficient Tumors

PURPOSE

Homologous recombination deficient (HRD) tumors are exquisitely sensitive to platinum-based chemotherapy and when combined with radiation therapy (RT), leads to improved overall survival in multiple cancer types. Whether a subset of tumors with distinct molecular characteristics demonstrate increased benefit from cisplatin and RT (c-RT) is unclear. We hypothesized that HRD tumors, whether associated with BRCA mutations or genomic scars of HRD, exhibit exquisite sensitivity to c-RT, and that HRD may be a significant driver of c-RT benefit.

METHODS AND MATERIALS

Sensitivity to c-RT was examined using isogenic and sporadic breast cancer cell lines. HRD was assessed using 4 assays: RT-induced Rad51 foci, a DR-GFP reporter assay, a genomic scar score (large-scale state transitions [LST]), and clonogenic survival assays. Whole-genome sequencing of 4 breast tumors from a phase 2 clinical trial of neoadjuvant c-RT in triple-negative breast cancer was performed and HRD was defined using HRDetect.

RESULTS

BRCA1/2 deficient cell lines displayed functional HRD based on the Rad51 functional assay, with c-RT to RT or cisplatin interaction ratios (IR) of 1.11 and 26.84 for the BRCA1 isogenic pair at 2 μM cisplatin and 6 Gy, respectively. The highest LST lines demonstrated HRD and synthetic cytotoxicity to c-RT with IR at 2 Gy and cisplatin 20 μM of 7.50, and the lowest LST line with IR of 0.65. Of 4 evaluable patients in the phase 2 trial, one achieved a pathologic complete response with corresponding HRD based on multiple genomic scar scores including HRDetect and LST scores, compared with patients without a pathologic complete response.

CONCLUSIONS

HRD breast cancers, whether identified by BRCA1/2 mutation status, functional tests, or mutational signatures, appear to be significantly more sensitive to c-RT compared with isogenic controls or tumors without HRD mutational signatures. HRD tumors may be exquisitely sensitive to c-RT which warrants further clinical investigation to guide a precision oncology approach.

Deciphering the Foundations of Mitochondrial Mutational Spectra: Replication-Driven and Damage-Induced Signatures Across Chordate Classes

Mitochondrial DNA (mtDNA) mutagenesis remains poorly understood despite its crucial role in disease, aging, and evolutionary tracing. In this study, we reconstructed a comprehensive 192-component mtDNA mutational spectrum for chordates by analyzing 118,397 synonymous mutations in the CytB gene across 1,697 species and five classes. This analysis revealed three primary forces shaping mtDNA mutagenesis: (i) symmetrical, replication-driven errors by mitochondrial polymerase (POLG), resulting in C > T and A > G mutations that are highly conserved across classes; (ii) asymmetrical, damage-driven C > T mutations on the single-stranded heavy strand with clock-like dynamics; and (iii) asymmetrical A > G mutations on the heavy strand, with dynamics suggesting sensitivity to oxidative damage. The third component, sensitive to oxidative damage, positions mtDNA mutagenesis as a promising marker for metabolic and physiological processes across various classes, species, organisms, tissues, and cells. The deconvolution of the mutational spectra into mutational signatures uncovered deficiencies in both base excision repair (BER) and mismatch repair (MMR) pathways. Further analysis of mutation hotspots, abasic sites, and mutational asymmetries underscores the critical role of single-stranded DNA damage (components ii and iii), which, uncorrected due to BER and MMR deficiencies, contributes roughly as many mutations as POLG-induced errors (component i).

IARC Workshop on the Key Characteristics of Carcinogens: Assessment of End Points for Evaluating Mechanistic Evidence of Carcinogenic Hazards

BACKGROUND

The 10 key characteristics (KCs) of carcinogens form the basis of a framework to identify, organize, and evaluate mechanistic evidence relevant to carcinogenic hazard identification. The 10 KCs are related to mechanisms by which carcinogens cause cancer. The International Agency for Research on Cancer (IARC) Monographs programme has successfully applied the KCs framework for the mechanistic evaluation of different types of exposures, including chemicals, metals, and complex exposures, such as environmental, occupational, or dietary exposures. The use of this framework has significantly enhanced the identification and organization of relevant mechanistic data, minimized bias in evaluations, and enriched the knowledge base regarding the mechanisms of known and suspected carcinogens.

OBJECTIVES

We sought to report the main outcomes of an IARC Scientific Workshop convened by the IARC to establish appropriate, transparent, and uniform application of the KCs in future IARC Monographs evaluations.

METHODS

A group of experts from different disciplines reviewed the IARC Monographs experience with the KCs of carcinogens, discussing three main themes: a) the interpretation of end points forming the evidence base for the KCs, b) the incorporation of data from novel assays on the KCs, and c) the integration of the mechanistic evidence as part of cancer hazard identification. The workshop participants assessed the relevance and the informativeness of multiple KCs-associated end points for the evaluation of mechanistic evidence in studies of exposed humans and experimental systems.

DISCUSSION

Consensus was reached on how to enhance the use of in silico, molecular, and cellular high-output and high-throughput data. In addition, approaches to integrate evidence across the KCs and opportunities to improve methodologies of mechanistic evaluation of cancer hazards were explored. The findings described herein and in a forthcoming IARC technical report will support future working groups of experts in reporting and interpreting results under the KCs framework within the IARC Monographs or in other contexts. https://doi.org/10.1289/EHP15389.

Ongoing chromothripsis underpins osteosarcoma genome complexity and clonal evolution

Osteosarcoma is the most common primary cancer of the bone, with a peak incidence in children and young adults. Using multi-region whole-genome sequencing, we find that chromothripsis is an ongoing mutational process, occurring subclonally in 74% of osteosarcomas. Chromothripsis generates highly unstable derivative chromosomes, the ongoing evolution of which drives the acquisition of oncogenic mutations, clonal diversification, and intra-tumor heterogeneity across diverse sarcomas and carcinomas. In addition, we characterize a new mechanism, termed loss-translocation-amplification (LTA) chromothripsis, which mediates punctuated evolution in about half of pediatric and adult high-grade osteosarcomas. LTA chromothripsis occurs when a single double-strand break triggers concomitant TP53 inactivation and oncogene amplification through breakage-fusion-bridge cycles. It is particularly prevalent in osteosarcoma and is not detected in other cancers driven by TP53 mutation. Finally, we identify the level of genome-wide loss of heterozygosity as a strong prognostic indicator for high-grade osteosarcoma.

Leveraging mRNA technology for antigen based immuno-oncology therapies

The application of messenger RNA (mRNA) technology in antigen-based immuno-oncology therapies represents a significant advancement in cancer treatment. Cancer vaccines are an effective combinatorial partner to sensitize the host immune system to the tumor and boost the efficacy of immune therapies. Selecting suitable tumor antigens is the key step to devising effective vaccinations and amplifying the immune response. Tumor neoantigens are de novo epitopes derived from somatic mutations, avoiding T-cell central tolerance of self-epitopes and inducing immune responses to tumors. The identification and prioritization of patient-specific tumor neoantigens are based on advanced computational algorithms taking advantage of the profiling with next-generation sequencing considering factors involved in human leukocyte antigen (HLA)-peptide-T-cell receptor (TCR) complex formation, including peptide presentation, HLA-peptide affinity, and TCR recognition. This review discusses the development and clinical application of mRNA vaccines in oncology, with a particular focus on recent clinical trials and the computational workflows and methodologies for identifying both shared and individual antigens. While this review centers on therapeutic mRNA vaccines targeting existing tumors, it does not cover preventative vaccines. Preclinical experimental validations are crucial in cancer vaccine development, but we emphasize the computational approaches that facilitate neoantigen selection and design, highlighting their role in advancing mRNA vaccine development. The versatility and rapid development potential of mRNA make it an ideal platform for personalized neoantigen immunotherapy. We explore various strategies for antigen target identification, including tumor-associated and tumor-specific antigens and the computational tools used to predict epitopes capable of eliciting strong immune responses. We address key design considerations for enhancing the immunogenicity and stability of mRNA vaccines, as well as emerging trends and challenges in the field. This comprehensive overview highlights the therapeutic potential of mRNA-based cancer vaccines and underscores ongoing research efforts aimed at optimizing these therapies for improved clinical outcomes.

Prolonged persistence of mutagenic DNA lesions in somatic cells

DNA is subject to continual damage, leaving each cell with thousands of individual DNA lesions at any given moment1-3. The efficiency of DNA repair means that most known classes of lesion have a half-life of minutes to hours3,4, but the extent to which DNA damage can persist for longer durations remains unknown. Here, using high-resolution phylogenetic trees from 89 donors, we identified mutations arising from 818 DNA lesions that persisted across multiple cell cycles in normal human stem cells from blood, liver and bronchial epithelium5-12. Persistent DNA lesions occurred at increased rates, with distinctive mutational signatures, in donors exposed to tobacco or chemotherapy, suggesting that they can arise from exogenous mutagens. In haematopoietic stem cells, persistent DNA lesions, probably from endogenous sources, generated the characteristic mutational signature SBS1913; occurred steadily throughout life, including in utero; and endured for 2.2 years on average, with 15-25% of lesions lasting at least 3 years. We estimate that on average, a haematopoietic stem cell has approximately eight such lesions at any moment in time, half of which will generate a mutation with each cell cycle. Overall, 16% of mutations in blood cells are attributable to SBS19, and similar proportions of driver mutations in blood cancers exhibit this signature. These data indicate the existence of a family of DNA lesions that arise from endogenous and exogenous mutagens, are present in low numbers per genome, persist for months to years, and can generate a substantial fraction of the mutation burden of somatic cells.

Assessment of candidate high-grade serous ovarian carcinoma predisposition genes through integrated germline and tumour sequencing

High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, only half of which is explained. Previously, we performed germline exome sequencing on BRCA1 and BRCA2-negative HGSOC patients, revealing three proposed and 43 novel candidate genes enriched with rare loss-of-function variants. For validation, we undertook case-control analyses using genomic data from disease-free controls. This confirms enrichment for nearly all previously identified genes. Additionally, one-hundred-and-eleven HGSOC tumours from variant carriers were sequenced alongside other complementary studies, seeking evidence of biallelic inactivation as supportive evidence. PALB2 and ATM validate as HGSOC predisposition genes, with 6/8 germline carrier tumours exhibiting biallelic inactivation accompanied by characteristic mutational signatures. Among candidate genes, only LLGL2 consistently shows biallelic inactivation and protein expression loss, supporting it as a novel HGSOC susceptibility gene. The remaining candidate genes fail to validate. Integrating case-control analyses with tumour sequencing is thus crucial for accurate gene discovery in familial cancer studies.

Investigating the origins of the mutational signatures in cancer

Most of the risk factors associated with chronic and complex diseases, such as cancer, stem from exogenous and endogenous exposures experienced throughout an individual's life, collectively known as the exposome. These exposures can modify DNA, which can subsequently lead to the somatic mutations found in all normal and tumor tissues. Understanding the precise origins of specific somatic mutations has been challenging due to multitude of DNA adducts (i.e. the DNA adductome) and their diverse positions within the genome. Thus far, this limitation has prevented researchers from precisely linking exposures to DNA adducts and DNA adducts to subsequent mutational outcomes. Indeed, many common mutations observed in human cancers appear to originate from error-prone endogenous processes. Consequently, it remains unclear whether these mutations result from exposure-induced DNA adducts, or arise indirectly from endogenous processes or are a combination of both. In this review, we summarize approaches that aim to bridge our understanding of the mechanism by which exposure leads to DNA damage and then to mutation and highlight some of the remaining challenges and shortcomings to fully supporting this paradigm. We emphasize the need to integrate cellular DNA adductomics, long read-based mapping, single-molecule duplex sequencing of native DNA molecules and advanced computational analysis. This proposed holistic approach aims to unveil the causal connections between key DNA modifications and the mutational landscape, whether they originate from external exposures, internal processes or a combination of both, thereby addressing key questions in cancer biology.

Establishment and characterization of a sigmoid colon cancer organoid with spinal metastasis

Background

Sigmoid colon cancer with spinal metastases is rare in distant metastasis. In addition, the prognosis of colon cancer patients with spinal metastases is extremely poor. In order to find effective therapeutic agents, we need to know the biological characteristics of such patients from related models.

Methods

We collected sigmoid colon cancer tissue from a young female subject who was diagnosed with sigmoid colon cancer with multiple spinal metastases. We successfully established a sigmoid colon cancer organoid using this tissue and investigated drug screening in the patient. HE staining, immunohistochemistry, and DNA sequencing were utilized to compare the biological characteristics between the original tumor and the organoid. Furthermore, we investigated the drug screening of the sigmoid colon cancer organoid in vitro.

Results

A colon cancer organoid from sigmoid colon cancer with spinal metastases was successfully established. The organoid culture maintained the morphological features, histological features, and genomic landscape of the corresponding sigmoid colon cancer cells. Moreover, we performed drug screening tests to evaluate the effects of chemotherapeutic drugs and targeted drugs.

Conclusion

The sigmoid colon cancer organoid with spinal metastases was a favorable preclinical model to explore the clinicopathologic characteristics of colon cancer patients with spinal metastases.

Tumor Mutation Signature Reveals the Risk Factors of Lung Adenocarcinoma with EGFR or KRAS Mutation

INTRODUCTION

EGFR and KRAS mutations are frequently detected in lung adenocarcinoma (LUAD). Tumor mutational signature (TMS) determination is an approach to identify somatic mutational patterns associated with pathogenic factors. In this study, through the analysis of TMS, the underlying pathogenic factors of LUAD with EGFR and KRAS mutations were traced.

METHODS

This was a retrospective study. TMS of LUAD with KRAS and EGFR mutations from the TCGA, OncoSG, and MSK datasets was determined by two bioinformatics tools, namely the "MutationalPatterns" and "FitMS" packages. Elevated microsatellite alterations at selected tetranucleotide repeats (EMAST) of LUAD clinical specimens was analyzed using capillary electrophoresis.

RESULTS

In LUAD with KRAS mutations, TMS analysis indicated that the smoking-related SBS4 signature was enriched. For LUAD with EGFR L858R mutation, the smoking-related SBS4 signature was enriched in the Western population from the TCGA database; however, the smoking-related SBS4 signature was not obvious in Asian LUAD patients. LUAD with EGFR exon19 deletion (19Del) exhibited stronger SBS15 signature, which was related to defective DNA mismatch repair. Capillary electrophoresis analysis showed that an EMAST locus was frequently instable in LUAD with EGFR 19Del. Different from the Western population, Asian LUAD patients with EGFR mutations exhibited the enrichment of SBS1, SBS2, and SBS13 signatures, which were associated with the endogenous mutation process of cytidine deamination.

CONCLUSIONS

TMS analysis reveals that smoking is associated with LUAD with KRAS mutations. Defective DNA mismatch repair and endogenous cytidine deamination are associated with LUAD with EGFR mutations, especially for the EGFR 19Del. The endogenous mutational process is stronger in Asian LUAD patients than Western LUAD patients.

Enhancing the detection of clinically relevant biomarkers in advanced uterine and tubo-ovarian carcinomas through genome-wide analysis

BACKGROUND

Advanced-stage tube-ovarian cancers (TOC) and uterine cancers (UC) significantly contribute to cancer mortality. While surgery achieves clinical remission in most cases, recurrence often necessitates systemic therapy. Recent molecular phenotype studies have advanced targeted therapies. We employed whole genome sequencing (WGS) to investigate biomarkers in gynecologic malignancies.

DESIGN

Ninety-one tumor samples (45 TOC, 46 UC) were analyzed for germline mutations, somatic driver mutations, rearrangements, genome-wide signatures, and molecular phenotypes. A WGS-based high-confidence classifier for homologous recombination deficiency (HRD) was applied. Genomic profiles were correlated with clinical outcomes.

RESULTS

The HRD phenotype was identified in serous carcinoma components, with 50 % of HRD cases showing BRCA1/2 wildtype (33 %) or variants of unknown significance (17 %). HRD correlated with better overall survival in tubo-ovarian serous carcinoma and was linked to responses to platinum therapy and PARP inhibitors. Endometrioid carcinomas showed no HRD phenotype despite BRCA1/2 variants. CDK12-type genomic instability (10 %) occurred in cases with CDK12 rearrangements, and CCNE1 gain (8 %) was observed in CCNE1-amplified cases, independent of copy number. In endometrioid carcinoma, mismatch repair (MMR) deficiency single base substitution signatures, particularly SBS44, contributed to high tumor mutation burden (TMB). Ultra-high TMB was found in two cases with pathogenic POLE and POLQ mutations, exhibiting multiple SBS signatures, including MMR deficiency.

CONCLUSION

Comprehensive genomic profiling of advanced-stage TOC and UC via WGS reveals key biomarkers and therapeutic targets, enhancing diagnostic accuracy and advancing personalized medicine in gynecological cancers.

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

SNMF: Integrated Learning of Mutational Signatures and Prediction of DNA Repair Deficiencies

Motivation

Many tumours show deficiencies in DNA damage response (DDR), which influence tumorigenesis and progression, but also expose vulnerabilities with therapeutic potential. Assessing which patients might benefit from DDR-targeting therapy requires knowledge of tumour DDR deficiency status, with mutational signatures reportedly better predictors than loss of function mutations in select genes. However, signatures are identified independently using unsupervised learning, and therefore not optimised to distinguish between different pathway or gene deficiencies.

Results

We propose SNMF, a supervised non-negative matrix factorisation that jointly optimises the learning of signatures: (1) shared across samples, and (2) predictive of DDR deficiency. We applied SNMF to mutation profiles of human induced pluripotent cell lines carrying gene knockouts linked to three DDR pathways. The SNMF model achieved high accuracy (0.971) and learned more complete signatures of the DDR status of a sample, further discerning distinct mechanisms within a pathway. Cell line SNMF signatures recapitulated tumour-derived COSMIC signatures and predicted DDR pathway deficiency of TCGA tumours with high recall, suggesting that SNMF-like models can leverage libraries of induced DDR deficiencies to decipher intricate DDR signatures underlying patient tumours.

Availability

https://github.com/joanagoncalveslab/SNMF .

Genomic landscape of cancer in racially and ethnically diverse populations

Cancer incidence and mortality rates can vary widely among different racial and ethnic groups, attributed to a complex interplay of genetic, environmental and social factors. Recently, substantial progress has been made in investigating hereditary genetic risk factors and in characterizing tumour genomes. However, most research has been conducted in individuals of European ancestries and, increasingly, in individuals of Asian ancestries. The study of germline and somatic genetics in cancer across racial and ethnic groups using omics technologies offers opportunities to identify similarities and differences in both heritable traits and the molecular features of cancer genomes. An improved understanding of population-specific cancer genomics, as well as translation of those findings across populations, will help reduce cancer disparities and ensure that personalized medicine and public health approaches are equitable across racial and ethnic groups.

Benchmarking 13 tools for mutational signature attribution, including a new and improved algorithm

Mutational signatures are characteristic patterns of mutations caused by endogenous mutational processes or by exogenous mutational exposures. There has been little benchmarking of approaches for determining which signatures are present in a sample and estimating the number of mutations due to each signature. This problem is referred to as "signature attribution." We show that there are often many combinations of signatures that can reconstruct the patterns of mutations in a sample reasonably well, even after encouraging sparse solutions. We benchmarked 13 approaches to signature attribution, including a new approach called Presence Attribute Signature Activity (PASA), on large synthetic data sets (2700 synthetic samples in total). These data sets recapitulated the single-base, insertion-deletion, and doublet-base mutational signature repertoires of nine cancer types. For single-base substitution mutations, PASA and MuSiCal outperformed other approaches on all the cancer types combined. However, the ranking of approaches varied by cancer type. For doublet-base substitutions and small insertions and deletions, while PASA outperformed the other approaches in most of the nine cancer types, the ranking of approaches again varied by cancer type. We believe that this variation reflects inherent difficulties in signature attribution. These difficulties stem from the fact that there are often many attributions that can reasonably explain the pattern of mutations in a sample and from the combinatorial search space due to the need to impose sparsity. Tables herein can provide guidance on the selection of mutational signature attribution approaches that are best suited to particular cancer types and study objectives.

Mutagenic impact and evolutionary influence of radiotherapy in hematologic malignancies

Ionizing radiotherapy (RT) is a widely used palliative and curative treatment strategy for malignancies. In solid tumors, RT-induced double strand breaks lead to the accumulation of indels, and their repair by non-homologous end-joining has been linked to the ID8 mutational signature in resistant cells. However, the extent of RT-induced DNA damage in hematologic malignancies and its impact on their evolution and interplay with commonly used chemotherapies has not yet been explored. Here, we interrogated 580 whole genome sequencing (WGS) from patients with large B-cell lymphoma, multiple myeloma, and myeloid neoplasms and identified ID8 only in relapsed disease. Yet, it was detected after exposure to both RT and mutagenic chemotherapy (i.e., platinum). Using WGS of single-cell colonies derived from treated lymphoma cells, we revealed a dose-response relationship between RT and platinum and ID8. Finally, using ID8 as a genomic barcode we demonstrate that a single RT-resistant cell may seed systemic relapse.

Base-excision repair pathway shapes 5-methylcytosine deamination signatures in pan-cancer genomes

Transition of cytosine to thymine in CpG dinucleotides is the most frequent type of mutation in cancer. This increased mutability is commonly attributed to the spontaneous deamination of 5-methylcytosine (5mC), which is normally repaired by the base-excision repair (BER) pathway. However, the contribution of 5mC deamination in the increasing diversity of cancer mutational signatures remains poorly explored. We integrate mutational signatures analysis in a large series of tumor whole genomes with lineage-specific epigenomic data to draw a detailed view of 5mC deamination in cancer. We uncover tumor type-specific patterns of 5mC deamination signatures in CpG and non-CpG contexts. We demonstrate that the BER glycosylase MBD4 preferentially binds to active chromatin and early replicating DNA, which correlates with lower mutational burden in these domains. We validate our findings by modeling BER deficiencies in isogenic cell models. Here, we establish MBD4 as the main actor responsible for 5mC deamination repair in humans.

A comprehensive comparison of tools for fitting mutational signatures

Mutational signatures connect characteristic mutational patterns in the genome with biological or chemical processes that take place in cancers. Analysis of mutational signatures can help elucidate tumor evolution, prognosis, and therapeutic strategies. Although tools for extracting mutational signatures de novo have been extensively benchmarked, a similar effort is lacking for tools that fit known mutational signatures to a given catalog of mutations. We fill this gap by comprehensively evaluating twelve signature fitting tools on synthetic mutational catalogs with empirically driven signature weights corresponding to eight cancer types. On average, SigProfilerSingleSample and SigProfilerAssignment/MuSiCal perform best for small and large numbers of mutations per sample, respectively. We further show that ad hoc constraining the list of reference signatures is likely to produce inferior results. Evaluation of real mutational catalogs suggests that the activity of signatures that are absent in the reference catalog poses considerable problems to all evaluated tools.

Origins and impact of extrachromosomal DNA

Extrachromosomal DNA (ecDNA) is a major contributor to treatment resistance and poor outcome for patients with cancer1,2. Here we examine the diversity of ecDNA elements across cancer, revealing the associated tissue, genetic and mutational contexts. By analysing data from 14,778 patients with 39 tumour types from the 100,000 Genomes Project, we demonstrate that 17.1% of tumour samples contain ecDNA. We reveal a pattern highly indicative of tissue-context-based selection for ecDNAs, linking their genomic content to their tissue of origin. We show that not only is ecDNA a mechanism for amplification of driver oncogenes, but it also a mechanism that frequently amplifies immunomodulatory and inflammatory genes, such as those that modulate lymphocyte-mediated immunity and immune effector processes. Moreover, ecDNAs carrying immunomodulatory genes are associated with reduced tumour T cell infiltration. We identify ecDNAs bearing only enhancers, promoters and lncRNA elements, suggesting the combinatorial power of interactions between ecDNAs in trans. We also identify intrinsic and environmental mutational processes linked to ecDNA, including those linked to its formation, such as tobacco exposure, and progression, such as homologous recombination repair deficiency. Clinically, ecDNA detection was associated with tumour stage, more prevalent after targeted therapy and cytotoxic treatments, and associated with metastases and shorter overall survival. These results shed light on why ecDNA is a substantial clinical problem that can cooperatively drive tumour growth signals, alter transcriptional landscapes and suppress the immune system.

TP53 mutation status and consensus molecular subtypes of colorectal cancer in patients from Rwanda

BACKGROUND

Mutations in the TP53 tumor suppressor gene are well-established drivers of colorectal cancer (CRC) development. However, data on the prevalence of TP53 variants and their association with consensus molecular subtype (CMS) classification in patients with CRC from Rwanda are currently lacking. This study addressed this knowledge gap by investigating TP53 mutation status concerning CMS classification in a CRC cohort from Rwanda.

METHODS

Formalin-fixed paraffin-embedded (FFPE) tissue blocks were obtained from 51 patients with CRC at the University Teaching Hospital of Kigali, Rwanda. Exons 4 to 11 and their flanking intron-exon boundaries in the TP53 gene were sequenced using Sanger sequencing to identify potential variants. The recently established immunohistochemistry-based classifier was employed to determine the CMS of each tumor.

RESULTS

Sequencing analysis of cancerous tissue DNA revealed TP53 pathogenic variants in 23 of 51 (45.1%) patients from Rwanda. These variants were predominantly missense types (18/23, 78.3%). The most frequent were c.455dup (p.P153Afs*28), c.524G > A (p.R175H), and c.733G > A (p.G245S), each identified in three tumors. Trinucleotide sequence context analysis of the 23 mutations (20 of which were single-base substitutions) revealed a predominance of the [C > N] pattern among single-base substitutions (SBSs) (18/20; 90.0%), with C[C > T]G being the most frequent mutation (5/18, 27.8%). Furthermore, pyrimidine bases (C and T) were preferentially found at the 5' flanking position of the mutated cytosine (13/18; 72.2%). Analysis of CMS subtypes revealed the following distribution: CMS1 (microsatellite instability-immune) (6/51, 11.8%), CMS2 (canonical) (28/51, 54.9%), CMS3 (metabolic) (9/51, 17.6%), and CMS4 (mesenchymal) (8/51, 15.7%). Interestingly, the majority of TP53 variants were in the CMS2 subgroup (14/23; 60.1%).

CONCLUSION

Our findings indicate a high frequency of TP53 variants in CRC patients from Rwanda. Importantly, these variants are enriched in the CMS2 subtype. This study, representing the second investigation into molecular alterations in patients with CRC from Rwanda and the first to explore TP53 mutations and CMS classification, provides valuable insights into the molecular landscape of CRC in this understudied population.

Next-Generation Integrated Sequencing Identifies Poor Prognostic Factors in Patients with MYD88-Mutated Chronic Lymphocytic Leukemia in Taiwan

INTRODUCTION

Chronic lymphocytic leukemia (CLL) is the most common type of leukemia in the Western countries and is very rare in Asia.

METHODS

Peripheral blood or bone marrow mononuclear cells obtained at initial diagnosis from 215 patients with CLL were analyzed by using next-generation sequencing to investigate the ethnic differences in genetic abnormalities.

RESULTS

Whole-genome sequencing and whole-exome sequencing analyses on 30 cases showed that 9 genes, including IGLL5, MYD88, TCHH, DSCAM, AXDND1, BICRA, KMT2D, MYT1L, and RBM43, were more frequently mutated in our Taiwanese cohort compared with those of the Western cohorts. IGLL5, MYD88, and KMT2D genes were further analyzed by targeted sequencing in another 185 CLL patients, unraveling frequencies of 29.3%, 20.9%, and 15.0%, respectively. The most frequent positional mutation of MYD88 was V217F (26/45, 57.8%), followed by L265P (9/45, 20.0%). MYD88 mutations were significantly associated with IGLL5 mutations (p = 0.0004), mutated IGHV (p < 0.0001) and 13q deletion (p = 0.0164). CLL patients with co-occurrence of MYD88 mutations with KMT2D or/and IGLL5 mutations were associated with a significantly inferior survival compared to those with MYD88 mutation alone (not reached vs. 131.8 months, p = 0.007). In multivariate analysis, MYD88 mutation without KMT2D or IGLL5 mutations was an independently favorable predictor.

CONCLUSIONS

IGLL5, MYD88, and KMT2D mutations were enriched in Taiwanese CLL, and co-occurrence of MYD88 mutations with KMT2D or/and IGLL5 mutations was associated with a poorer prognosis.

Tumour mutational burden: clinical utility, challenges and emerging improvements

Tumour mutational burden (TMB), defined as the total number of somatic non-synonymous mutations present within the cancer genome, varies across and within cancer types. A first wave of retrospective and prospective research identified TMB as a predictive biomarker of response to immune-checkpoint inhibitors and culminated in the disease-agnostic approval of pembrolizumab for patients with TMB-high tumours based on data from the Keynote-158 trial. Although the applicability of outcomes from this trial to all cancer types and the optimal thresholds for TMB are yet to be ascertained, research into TMB is advancing along three principal avenues: enhancement of TMB assessments through rigorous quality control measures within the laboratory process, including the mitigation of confounding factors such as limited panel scope and low tumour purity; refinement of the traditional TMB framework through the incorporation of innovative concepts such as clonal, persistent or HLA-corrected TMB, tumour neoantigen load and mutational signatures; and integration of TMB with established and emerging biomarkers such as PD-L1 expression, microsatellite instability, immune gene expression profiles and the tumour immune contexture. Given its pivotal functions in both the pathogenesis of cancer and the ability of the immune system to recognize tumours, a profound comprehension of the foundational principles and the continued evolution of TMB are of paramount relevance for the field of oncology.

Prevalence of germline variants in Brazilian pancreatic carcinoma patients

We evaluated the prevalence of pathogenic/likely pathogenic germline variants (PGV) in Brazilian pancreatic adenocarcinoma (PC) patients, that represent a multiethnic population, in a cross-sectional study. We included 192 PC patients unselected for family history of cancer. We evaluated a panel of 113 cancer genes, through genomic DNA sequencing and 46 ancestry-informative markers, through multiplex PCR. The median age was 61 years; 63.5% of the patients presented disease clinical stages III or IV; 8.3% reported personal history of cancer; 4.7% and 16.1% reported first-degree relatives with PC or breast and/or prostate cancer, respectively. Although the main ancestry was European, there was considerable genetic composition admixture. Twelve patients (6.25%) were PGV carriers in PC predisposition genes (ATM, BRCA1, BRCA2, CDKN2A, MSH2, PALB2) and another 25 (13.0%) were PGV carriers in genes with a limited association or not previously associated with PC (ACD, BLM, BRIP1, CHEK2, ERCC4, FANCA, FANCE, FANCM, GALNT12, MITF, MRE11, MUTYH, POLE, RAD51B, RAD51C, RECQL4, SDHA, TERF2IP). The most frequently affected genes were CHEK2, ATM and FANC. In tumor samples from PGV carriers in ACD, BRIP1, MRE11, POLE, SDHA, TERF2IP, which were examined through exome sequencing, the main single base substitutions (SBS) mutational signature was SBS1+5+18, probably associated with age, tobacco smoking and reactive oxygen species. SBS3 associated with homologous repair deficiency was also represented, but on a lower scale. There was no difference in the frequency of PGV carriers between: (a) patients with or without first-degree relatives with cancer; and (b) patients with admixed ancestry versus those with predominantly European ancestry. Furthermore, there was no difference in overall survival between PGV carriers and non-carriers. Therefore, genetic testing should be offered to all Brazilian pancreatic cancer patients, regardless of their ancestry. Genes with limited or previously unrecognized associations with pancreatic cancer should be further investigated to clarify their role in cancer risk.

Mutational Features and Tumor Microenvironment Alterations in High-Grade Appendiceal Cancers Treated With Iterative Hyperthermic Intraperitoneal Chemotherapy

PURPOSE

High-grade appendiceal adenocarcinomas (HGAA) with peritoneal metastases (PMs) are associated with poor survival. Hyperthermic intraperitoneal chemotherapy (HIPEC) is a novel treatment approach for unresectable HGAA-PM. However, its influence on immunogenomic profiles has not yet been fully explored.

MATERIALS AND METHODS

We obtained 79 samples of metastatic peritoneal tumor deposits from patients diagnosed with HGAA and performed whole-exome sequencing, RNA sequencing, and immunoprofiling before and after HIPEC. Tumor biopsies were subjected to immunogenomic profiling to detect mutational signatures and immune populations associated with oncologic outcomes.

RESULTS

Fifteen patients with HGAA-PMs were included in the study. The median progression-free survival (PFS) was 6.7 months (2.7-25.3) and the median overall survival was 11.4 months (4.7-42). Mucin-associated genes (MUC16, MUC3A, and MUC5AC) and titin (TTN) had the highest mutation frequencies. Mutational signatures such as single-base substitution 29 and doublet-base substitution 11 were present in >50% of single-base and double-base mutations. Higher PD-L1 coexpression on CD8+ T cells demonstrated a higher PFS both intratumorally (P = .019) and at the margin (P = .025).

CONCLUSION

HIPEC-associated mutational signatures were identified in HGAA-PMs. Elevated PD-L1+ cytotoxic T-cell populations after HIPEC had better PFS, offering valuable insights for prognostication in the context of HIPEC treatment.

Prognostic genome and transcriptome signatures in colorectal cancers

Colorectal cancer is caused by a sequence of somatic genomic alterations affecting driver genes in core cancer pathways1. Here, to understand the functional and prognostic impact of cancer-causing somatic mutations, we analysed the whole genomes and transcriptomes of 1,063 primary colorectal cancers in a population-based cohort with long-term follow-up. From the 96 mutated driver genes, 9 were not previously implicated in colorectal cancer and 24 had not been linked to any cancer. Two distinct patterns of pathway co-mutations were observed, timing analyses identified nine early and three late driver gene mutations, and several signatures of colorectal-cancer-specific mutational processes were identified. Mutations in WNT, EGFR and TGFβ pathway genes, the mitochondrial CYB gene and 3 regulatory elements along with 21 copy-number variations and the COSMIC SBS44 signature correlated with survival. Gene expression classification yielded five prognostic subtypes with distinct molecular features, in part explained by underlying genomic alterations. Microsatellite-instable tumours divided into two classes with different levels of hypoxia and infiltration of immune and stromal cells. To our knowledge, this study constitutes the largest integrated genome and transcriptome analysis of colorectal cancer, and interlinks mutations, gene expression and patient outcomes. The identification of prognostic mutations and expression subtypes can guide future efforts to individualize colorectal cancer therapy.

Mutational Signatures in Colorectal Cancer: Translational Insights, Clinical Applications, and Limitations

A multitude of exogenous and endogenous processes have the potential to result in DNA damage. While the repair mechanisms are typically capable of correcting this damage, errors in the repair process can result in mutations. The findings of research conducted in 2012 indicate that mutations do not occur randomly but rather follow specific patterns that can be attributed to known or inferred mutational processes. The process of mutational signature analysis allows for the inference of the predominant mutational process for a given cancer sample, with significant potential for clinical applications. A deeper comprehension of these mutational signatures in CRC could facilitate enhanced prevention strategies, facilitate the comprehension of genotoxic drug activity, predict responses to personalized treatments, and, in the future, inform the development of targeted therapies in the context of precision oncology. The efforts of numerous researchers have led to the identification of several mutational signatures, which can be categorized into different mutational signature references. In CRC, distinct mutational signatures are identified as correlating with mismatch repair deficiency, polymerase mutations, and chemotherapy treatment. In this context, a mutational signature analysis offers considerable potential for enhancing minimal residual disease (MRD) tests in stage II (high-risk) and stage III CRC post-surgery, stratifying CRC based on the impacts of genetic and epigenetic alterations for precision oncology, identifying potential therapeutic vulnerabilities, and evaluating drug efficacy and guiding therapy, as illustrated in a proof-of-concept clinical trial.

The cytidine deaminase APOBEC3C has unique sequence and genome feature preferences

APOBEC proteins are cytidine deaminases that restrict the replication of viruses and transposable elements. Several members of the APOBEC3 family, APOBEC3A, APOBEC3B, and APOBEC3H-I, can access the nucleus and cause what is thought to be indiscriminate deamination of the genome, resulting in mutagenesis and genome instability. Although APOBEC3C is also present in the nucleus, the full scope of its deamination target preferences is unknown. By expressing human APOBEC3C in a yeast model system, I have defined the APOBEC3C mutation signature, as well as the preferred genome features of APOBEC3C targets. The APOBEC3C mutation signature is distinct from those of the known cancer genome mutators APOBEC3A and APOBEC3B. APOBEC3C produces DNA strand-coordinated mutation clusters, and APOBEC3C mutations are enriched near the transcription start sites of active genes. Surprisingly, APOBEC3C lacks the bias for the lagging strand of DNA replication that is seen for APOBEC3A and APOBEC3B. The unique preferences of APOBEC3C constitute a mutation profile that will be useful in defining sites of APOBEC3C mutagenesis in human genomes.

Comprehensive investigation of the mutagenic potential of six pesticides classified by IARC as probably carcinogenic to humans

Pesticides are significant environmental pollutants, and many of them possess mutagenic potential, which is closely linked to carcinogenesis. Here we tested the mutagenicity of all six pesticides classified probably carcinogenic (Group 2A) by the International Agency of Research on Cancer: 4,4'-DDT, captafol, dieldrin, diazinon, glyphosate and malathion. Whole genome sequencing of TK6 human lymphoblastoid cell clones following 30-day exposure at subtoxic concentrations revealed a clear mutagenic effect of treatment with captafol or malathion when added at 200 nM or 100 μM initial concentrations, respectively. Each pesticide induced a specific base substitution mutational signature: captafol increased C to A mutations primarily, while malathion induced mostly C to T mutations. 4,4'-DDT, dieldrin, diazinon and glyphosate were not mutagenic. Whereas captafol induced chromosomal instability, H2A.X phosphorylation and cell cycle arrest in G2/M phase, all indicating DNA damage, malathion did not induce DNA damage markers or cell cycle alterations despite its mutagenic effect. Hypersensitivity of REV1 and XPA mutant DT40 chicken cell lines suggests that captafol induces DNA adducts that are bypassed by translesion DNA synthesis and are targets for nucleotide excision repair. The experimentally identified mutational signatures of captafol and malathion could shed light on the mechanism of action of these compounds. The signatures are potentially suitable for detecting past exposure in tumour samples, but the reanalysis of large cancer genome databases did not reveal any evidence of captafol or malathion exposure.

Pan-cancer mutational signature analysis of 111,711 targeted sequenced tumors using SATS

Tumor mutational signatures are informative for cancer diagnosis and treatment. However, targeted sequencing, commonly used in clinical settings, lacks specialized analytical tools and a dedicated catalogue of mutational signatures. Here, we introduce SATS, a scalable mutational signature analyzer for targeted sequencing data. SATS leverages tumor mutational burdens to identify and quantify signatures in individual tumors, overcoming the challenges of sparse mutations and variable gene panels. Validations across simulated data, pseudo-targeted sequencing data, and matched whole-genome and targeted sequencing samples show that SATS can accurately detect common mutational signatures and estimate their burdens. Applying SATS to 111,711 tumors from the AACR Project GENIE, we created a pan-cancer mutational signature catalogue specific to targeted sequencing. We further validated signatures in lung, breast and colorectal cancers using an additional 16,774 independent samples. This signature catalogue is a valuable resource for estimating signature burdens in individual targeted sequenced tumors, facilitating the integration of mutational signatures with clinical data.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Neotelomeres and telomere-spanning chromosomal arm fusions in cancer genomes revealed by long-read sequencing

Alterations in the structure and location of telomeres are pivotal in cancer genome evolution. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeats, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. These results provide a framework for the systematic study of telomeric repeats in cancer genomes, which could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Passive Smoking-Induced Mutagenesis as a Promoter of Lung Carcinogenesis

INTRODUCTION

The International Agency for Research on Cancer has classified passive smoking (PS) or secondhand smoke exposure as a group 1 carcinogen linked to lung cancer. However, in contrast to active smoking, the mutagenic properties of PS remain unclear.

METHODS

A consecutive cohort of 564 lung adenocarcinoma samples from female never-smokers, who provided detailed information about their exposure to PS during adolescence and in their thirties through a questionnaire, was prepared. Of these, all 291 cases for whom frozen tumor tissues were available were subjected to whole exome sequencing to estimate tumor mutational burden, and the top 84 cases who were exposed daily, or not, to PS during adolescence, in their thirties or in both periods, were further subjected to whole genome sequencing.

RESULTS

A modest yet statistically significant increase in tumor mutational burden was observed in the group exposed to PS compared with the group not exposed to PS (median values = 1.44 versus 1.29 per megabase, respectively; p = 0.020). Instead of inducing driver oncogene mutations, PS-induced substantial subclonal mutations exhibiting APOBEC-type signatures, including SMAD4 and ADGRG6 hotspot mutations. A polymorphic APOBEC3A/3B allele-specific to the Asian population that leads to up-regulated expression of APOBEC3A accentuated the mutational load in individuals exposed daily to PS during adolescence.

CONCLUSIONS

This study reveals that PS-induced mutagenesis can promote lung carcinogenesis. The APOBEC3A/3B polymorphism may serve as a biomarker for identifying passive nonsmoking individuals at high risk of developing lung cancer.

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Impact of genomic and epigenomic alterations of multigene on a multicancer pedigree

BACKGROUND

Germline mutations have been identified in a small number of hereditary cancers, but the genetic predisposition for many familial cancers remains to be elucidated.

METHODS

This study identified a Chinese pedigree that presented different cancers (breast cancer, BRCA; adenocarcinoma of the esophagogastric junction, AEG; and B-cell acute lymphoblastic leukemia, B-ALL) in each of the three generations. Whole-genome sequencing and whole-exome sequencing were performed on peripheral blood or bone marrow and cancer biopsy samples. Whole-genome bisulfite sequencing was conducted on the monozygotic twin brothers, one of whom developed B-ALL.

RESULTS

According to the ACMG guidelines, bioinformatic analysis of the genome sequencing revealed 20 germline mutations, particularly mutations in the DNAH11 (c.9463G > A) and CFH (c.2314G > A) genes that were documented in the COSMIC database and validated by Sanger sequencing. Forty-one common somatic mutated genes were identified in the cancer samples, displaying the same type of single nucleotide substitution Signature 5. Meanwhile, hypomethylation of PLEK2, MRAS, and RXRA as well as hypermethylation of CpG island associated with WT1 was shown in the twin with B-ALL.

CONCLUSIONS

These findings reveal genomic alterations in a pedigree with multiple cancers. Mutations found in the DNAH11, CFH genes, and other genes predispose to malignancies in this family. Dysregulated methylation of WT1, PLEK2, MRAS, and RXRA in the twin with B-ALL increases cancer susceptibility. The similarity of the somatic genetic changes among the three cancers indicates a hereditary impact on the pedigree. These familial cancers with germline and somatic mutations, as well as epigenomic alterations, represent a common molecular basis for many multiple cancer pedigrees.

Human cancer genomes harbor the mutational signature of tobacco-specific nitrosamines NNN and NNK

Tobacco usage is linked to multiple cancer types and accounts for a quarter of all cancer-related deaths. Tobacco smoke contains various carcinogenic compounds, including polycyclic aromatic hydrocarbons (PAH), though the mutagenic potential of many tobacco-related chemicals remains largely unexplored. In particular, the highly carcinogenic tobacco-specific nitrosamines NNN and NNK form pre-mutagenic pyridyloxobutyl (POB) DNA adducts. In the study presented here, we identified genome-scale POB-induced mutational signatures in cell lines and rat tumors, while also investigating their role in human cancer. These signatures are characterized by T>N and C>T mutations forming from specific POB adducts damaging dT and dC residues. Analysis of 2,780 cancer genomes uncovered POB signatures in ∼180 tumors; from cancer types distinct from the ones linked to smoking-related signatures SBS4 and SBS92. This suggests that, unlike PAH compounds, the POB pathway may contribute uniquely to the mutational landscapes of certain hematological malignancies and cancers of the kidney, breast, prostate and pancreas.

Tumor initiation and early tumorigenesis: molecular mechanisms and interventional targets

Tumorigenesis is a multistep process, with oncogenic mutations in a normal cell conferring clonal advantage as the initial event. However, despite pervasive somatic mutations and clonal expansion in normal tissues, their transformation into cancer remains a rare event, indicating the presence of additional driver events for progression to an irreversible, highly heterogeneous, and invasive lesion. Recently, researchers are emphasizing the mechanisms of environmental tumor risk factors and epigenetic alterations that are profoundly influencing early clonal expansion and malignant evolution, independently of inducing mutations. Additionally, clonal evolution in tumorigenesis reflects a multifaceted interplay between cell-intrinsic identities and various cell-extrinsic factors that exert selective pressures to either restrain uncontrolled proliferation or allow specific clones to progress into tumors. However, the mechanisms by which driver events induce both intrinsic cellular competency and remodel environmental stress to facilitate malignant transformation are not fully understood. In this review, we summarize the genetic, epigenetic, and external driver events, and their effects on the co-evolution of the transformed cells and their ecosystem during tumor initiation and early malignant evolution. A deeper understanding of the earliest molecular events holds promise for translational applications, predicting individuals at high-risk of tumor and developing strategies to intercept malignant transformation.

Cancer mutational signatures identification in clinical assays using neural embedding-based representations

While mutational signatures provide a plethora of prognostic and therapeutic insights, their application in clinical-setting, targeted gene panels is extremely limited. We develop a mutational representation model (which learns and embeds specific mutation signature connections) that enables prediction of dominant signatures with only a few mutations. We predict the dominant signatures across more than 60,000 tumors with gene panels, delineating their landscape across different cancers. Dominant signature predictions in gene panels are of clinical importance. These included UV, tobacco, and apolipoprotein B mRNA editing enzyme, catalytic polypeptide (APOBEC) signatures that are associated with better survival, independently from mutational burden. Further analyses reveal gene and mutation associations with signatures, such as SBS5 with TP53 and APOBEC with FGFR3S249C. In a clinical use case, APOBEC signature is a robust and specific predictor for resistance to epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs). Our model provides an easy-to-use way to detect signatures in clinical setting assays with many possible clinical implications for an unprecedented number of cancer patients.

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.

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

The mutagenic forces shaping the genomic landscape of lung cancer in never smokers

Lung cancer in never smokers (LCINS) accounts for up to 25% of all lung cancers and has been associated with exposure to secondhand tobacco smoke and air pollution in observational studies. Here, we evaluate the mutagenic exposures in LCINS by examining deep whole-genome sequencing data from a large international cohort of 871 treatment-naïve LCINS recruited from 28 geographical locations within the Sherlock-Lung study. KRAS mutations were 3.8-fold more common in adenocarcinomas of never smokers from North America and Europe, while a 1.6-fold higher prevalence of EGFR and TP53 mutations was observed in adenocarcinomas from East Asia. Signature SBS40a, with unknown cause, was found in most samples and accounted for the largest proportion of single base substitutions in adenocarcinomas, being enriched in EGFR-mutated cases. Conversely, the aristolochic acid signature SBS22a was almost exclusively observed in patients from Taipei. Even though LCINS exposed to secondhand smoke had an 8.3% higher mutational burden and 5.4% shorter telomeres, passive smoking was not associated with driver mutations in cancer driver genes or the activities of individual mutational signatures. In contrast, patients from regions with high levels of air pollution were more likely to have TP53 mutations while exhibiting shorter telomeres and an increase in most types of somatic mutations, including a 3.9-fold elevation of signature SBS4 (q-value=3.1 × 10-5), previously linked mainly to tobacco smoking, and a 76% increase of clock-like signature SBS5 (q-value=5.0 × 10-5). A positive dose-response effect was observed with air pollution levels, which correlated with both a decrease in telomere length and an elevation in somatic mutations, notably attributed to signatures SBS4 and SBS5. Our results elucidate the diversity of mutational processes shaping the genomic landscape of lung cancer in never smokers.

G-Quadruplex Forming DNA Sequence Context Is Enriched around Points of Somatic Mutations in a Subset of Multiple Myeloma Patients

Multiple myeloma (MM) is the second most common hematological malignancy, which remains incurable despite recent advances in treatment strategies. Like other forms of cancer, MM is characterized by genomic instability, caused by defects in DNA repair. Along with mutations in DNA repair genes and genotoxic drugs used to treat MM, non-canonical secondary DNA structures (four-stranded G-quadruplex structures) can affect accumulation of somatic mutations and chromosomal abnormalities in the tumor cells of MM patients. Here, we tested the hypothesis that G-quadruplex structures may influence the distribution of somatic mutations in the tumor cells of MM patients. We sequenced exomes of normal and tumor cells of 11 MM patients and analyzed the data for the presence of G4 context around points of somatic mutations. To identify molecular mechanisms that could affect mutational profile of tumors, we also analyzed mutational signatures in tumor cells as well as germline mutations for the presence of specific SNPs in DNA repair genes or in genes regulating G-quadruplex unwinding. In several patients, we found that sites of somatic mutations are frequently located in regions with G4 context. This pattern correlated with specific germline variants found in these patients. We discuss the possible implications of these variants for mutation accumulation and specificity in MM and propose that the extent of G4 context enrichment around somatic mutation sites may be a novel metric characterizing mutational processes in tumors.

MUSE-XAE: MUtational Signature Extraction with eXplainable AutoEncoder enhances tumour types classification

MOTIVATION

Mutational signatures are a critical component in deciphering the genetic alterations that underlie cancer development and have become a valuable resource to understand the genomic changes during tumorigenesis. Therefore, it is essential to employ precise and accurate methods for their extraction to ensure that the underlying patterns are reliably identified and can be effectively utilized in new strategies for diagnosis, prognosis, and treatment of cancer patients.

RESULTS

We present MUSE-XAE, a novel method for mutational signature extraction from cancer genomes using an explainable autoencoder. Our approach employs a hybrid architecture consisting of a nonlinear encoder that can capture nonlinear interactions among features, and a linear decoder which ensures the interpretability of the active signatures. We evaluated and compared MUSE-XAE with other available tools on both synthetic and real cancer datasets and demonstrated that it achieves superior performance in terms of precision and sensitivity in recovering mutational signature profiles. MUSE-XAE extracts highly discriminative mutational signature profiles by enhancing the classification of primary tumour types and subtypes in real world settings. This approach could facilitate further research in this area, with neural networks playing a critical role in advancing our understanding of cancer genomics.

AVAILABILITY AND IMPLEMENTATION

MUSE-XAE software is freely available at https://github.com/compbiomed-unito/MUSE-XAE.

A glycolytic metabolite bypasses "two-hit" tumor suppression by BRCA2

Knudson's "two-hit" paradigm posits that carcinogenesis requires inactivation of both copies of an autosomal tumor suppressor gene. Here, we report that the glycolytic metabolite methylglyoxal (MGO) transiently bypasses Knudson's paradigm by inactivating the breast cancer suppressor protein BRCA2 to elicit a cancer-associated, mutational single-base substitution (SBS) signature in nonmalignant mammary cells or patient-derived organoids. Germline monoallelic BRCA2 mutations predispose to these changes. An analogous SBS signature, again without biallelic BRCA2 inactivation, accompanies MGO accumulation and DNA damage in Kras-driven, Brca2-mutant murine pancreatic cancers and human breast cancers. MGO triggers BRCA2 proteolysis, temporarily disabling BRCA2's tumor suppressive functions in DNA repair and replication, causing functional haploinsufficiency. Intermittent MGO exposure incites episodic SBS mutations without permanent BRCA2 inactivation. Thus, a metabolic mechanism wherein MGO-induced BRCA2 haploinsufficiency transiently bypasses Knudson's two-hit requirement could link glycolysis activation by oncogenes, metabolic disorders, or dietary challenges to mutational signatures implicated in cancer evolution.

Novel insights into tumorigenesis revealed by molecular analysis of Lynch syndrome cases with multiple colorectal tumors

Background

Lynch syndrome (LS) is an autosomal dominant multi-organ cancer syndrome with a high lifetime risk of cancer. The number of cumulative colorectal adenomas in LS does not generally exceed ten, and removal of adenomas via routine screening minimizes the cancer burden. However, abnormal phenotypes may mislead initial diagnosis and subsequently cause suboptimal treatment.

Aim

Currently, there is no standard guide for the care of multiple colorectal adenomas in LS individuals. We aimed to shed insight into the molecular features and reasons for multiplicity of adenomas in LS patients.

Methods

We applied whole exome sequencing on nine adenomas (ten samples) and three assumed primary carcinomas (five samples) of an LS patient developing the tumors during a 21-year follow-up period. We compared the findings to the tumor profiles of two additional LS cases ascertained through colorectal tumor multiplicity, as well as to ten adenomas and 15 carcinomas from 23 unrelated LS patients with no elevated adenoma burden from the same population. As LS associated cancers can arise via several molecular pathways, we also profiled the tumors for CpG Island Methylator Phenotype (CIMP), and LINE-1 methylation.

Results

All tumors were microsatellite unstable (MSI), and MSI was present in several samples derived from normal mucosa as well. Interestingly, frequent frameshift variants in RNF43 were shared among substantial number of the tumors of our primary case and the tumors of LS cases with multiple tumors but almost absent in our control LS cases. The RNF43 variants were completely absent in the normal tissue, indicating tumor-associated mutational hotspots. The RNF43 status correlated with the mutational signature SBS96. Contrary to LS tumors from the reference set with no elevated colorectal tumor burden, the somatic variants occurred significantly more frequently at C>T in the CpG context, irrespective of CIMP or LINE-1 status, potentially indicating other, yet unknown methylation-related mechanisms. There were no signs of somatic mosaicism affecting the MMR genes. Somatic variants in APC and CTNNB1 were unique to each tumor.

Conclusion

Frequent somatic RNF43 hot spot variants combined with SBS96 signature and increased tendency to DNA methylation may contribute to tumor multiplicity in LS.

Molecular analysis of cancer genomes in children with Lynch syndrome: Exploring causal associations

Lynch syndrome (LS) predisposes to cancer in adulthood and is caused by heterozygous germline variants in a mismatch repair (MMR) gene. Recent studies show an increased prevalence of LS among children with cancer, suggesting a causal relationship. For LS-spectrum (LSS) cancers, including high-grade gliomas and colorectal cancer, causality has been supported by typical MMR-related tumor characteristics, but for non-LSS cancers, causality is unclear. We characterized 20 malignant tumors of 18 children with LS, including 16 non-LSS tumors. We investigated second hits, tumor mutational load, mutational signatures and MMR protein expression. In all LSS tumors and three non-LSS tumors, we detected MMR deficiency caused by second hit somatic alterations. Furthermore, these MMR-deficient tumors carried driver variants that likely originated as a consequence of MMR deficiency. However, in 13 non-LSS tumors (81%), a second hit and MMR deficiency were absent, thus a causal link between LS and cancer development in these children is lacking. These findings demonstrate that causality of LS in children with cancer, which can be determined by molecular tumor characterization, seems to be restricted to specific tumor types. Large molecular and epidemiological studies are needed to further refine the tumor spectrum in children with LS.

Genetic interactions reveal distinct biological and therapeutic implications in breast cancer

Co-occurrence and mutual exclusivity of genomic alterations may reflect the existence of genetic interactions, potentially shaping distinct biological phenotypes and impacting therapeutic response in breast cancer. However, our understanding of them remains limited. Herein, we investigate a large-scale multi-omics cohort (n = 873) and a real-world clinical sequencing cohort (n = 4,405) including several clinical trials with detailed treatment outcomes and perform functional validation in patient-derived organoids, tumor fragments, and in vivo models. Through this comprehensive approach, we construct a network comprising co-alterations and mutually exclusive events and characterize their therapeutic potential and underlying biological basis. Notably, we identify associations between TP53mut-AURKAamp and endocrine therapy resistance, germline BRCA1mut-MYCamp and improved sensitivity to PARP inhibitors, and TP53mut-MYBamp and immunotherapy resistance. Furthermore, we reveal that precision treatment strategies informed by co-alterations hold promise to improve patient outcomes. Our study highlights the significance of genetic interactions in guiding genome-informed treatment decisions beyond single driver alterations.

Familial uveal melanoma and other tumors in 25 families with monoallelic germline MBD4 variants

BACKGROUND

Monoallelic germline MBD4 pathogenic variants were recently reported to cause a predisposition to uveal melanoma, associated with a specific tumor mutational signature and good response to immunotherapy. Monoallelic tumor pathogenic variants have also been described in brain tumors, breast cancers, and myxofibrosarcomas, whereas biallelic germline MBD4 pathogenic variants have been involved in a recessive hereditary adenomatous polyposis and a specific type of acute myeloid leukemia.

METHODS

We analyzed MBD4 for all patients with a diagnosis of uveal melanoma at Institut Curie since July 2021 and in the 3240 consecutive female probands explored at the Institut Curie for suspicion of predisposition to breast cancer between July 2021 and February 2023.

RESULTS

We describe 25 families whose probands carry a monoallelic germline pathogenic variant in MBD4. Eighteen of these families presented with uveal melanoma (including a case patient with multiple uveal melanoma), and 7 families presented with breast cancer. Family histories showed the first familial case of uveal melanoma in monoallelic MBD4 pathogenic variant carriers and other various types of cancers in relatives, especially breast, renal, and colorectal tumors.

CONCLUSIONS

Monoallelic MBD4 pathogenic variant may explain some cases of familial and multiple uveal melanoma as well as various cancer types, expanding the tumor spectrum of this predisposition. Further genetic testing in relatives combined with molecular tumor analyses will help define the tumor spectrum and estimate each tumor's risk.

Extrachromosomal DNA in cancer

Extrachromosomal DNA (ecDNA) has recently been recognized as a major contributor to cancer pathogenesis that is identified in most cancer types and is associated with poor outcomes. When it was discovered over 60 years ago, ecDNA was considered to be rare, and its impact on tumour biology was not well understood. The application of modern imaging and computational techniques has yielded powerful new insights into the importance of ecDNA in cancer. The non-chromosomal inheritance of ecDNA during cell division results in high oncogene copy number, intra-tumoural genetic heterogeneity and rapid tumour evolution that contributes to treatment resistance and shorter patient survival. In addition, the circular architecture of ecDNA results in altered patterns of gene regulation that drive elevated oncogene expression, potentially enabling the remodelling of tumour genomes. The generation of clusters of ecDNAs, termed ecDNA hubs, results in interactions between enhancers and promoters in trans, yielding a new paradigm in oncogenic transcription. In this Review, we highlight the rapid advancements in ecDNA research, providing new insights into ecDNA biogenesis, maintenance and transcription and its role in promoting tumour heterogeneity. To conclude, we delve into a set of unanswered questions whose answers will pave the way for the development of ecDNA targeted therapeutic approaches.