Mutational signature analysis identifies MUTYH deficiency in colorectal cancers and adrenocortical carcinomas

Spatial intra-tumour heterogeneity and treatment-induced genomic evolution in oesophageal adenocarcinoma: implications for prognosis and therapy

BACKGROUND

Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.

METHODS

Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define 'shared' (between both samples) and 'private' (present in one sample) mutations.

RESULTS

Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.

CONCLUSIONS

This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

The spectrum of TP53 mutations in Rwandan patients with gastric cancer

BACKGROUND

Gastric cancer is the sixth most frequently diagnosed cancer and third in causing cancer-related death globally. The most frequently mutated gene in human cancers is TP53, which plays a pivotal role in cancer initiation and progression. In Africa, particularly in Rwanda, data on TP53 mutations are lacking. Therefore, this study intended to obtain TP53 mutation status in Rwandan patients with gastric cancer.

RESULTS

Formalin-fixed paraffin-embedded tissue blocks of 95 Rwandan patients with histopathologically proven gastric carcinoma were obtained from the University Teaching Hospital of Kigali. After DNA extraction, all coding regions of the TP53 gene and the exon-intron boundary region of TP53 were sequenced using the Sanger sequencing. Mutated TP53 were observed in 24 (25.3%) of the 95 cases, and a total of 29 mutations were identified. These TP53 mutations were distributed between exon 4 and 8 and most of them were missense mutations (19/29; 65.5%). Immunohistochemical analysis for TP53 revealed that most of the TP53 missense mutations were associated with TP53 protein accumulation. Among the 29 mutations, one was novel (c.459_477delCGGCACCCGCGTCCGCGCC). This 19-bp deletion mutation in exon 5 caused the production of truncated TP53 protein (p.G154Wfs*10). Regarding the spectrum of TP53 mutations, G:C > A:T at CpG sites was the most prevalent (10/29; 34.5%) and G:C > T:A was the second most prevalent (7/29; 24.1%). Interestingly, when the mutation spectrum of TP53 was compared to three previous TP53 mutational studies on non-Rwandan patients with gastric cancer, G:C > T:A mutations were significantly more frequent in this study than in our previous study (p = 0.013), the TCGA database (p = 0.017), and a previous study on patients from Hong Kong (p = 0.006). Even after correcting for false discovery, statistical significance was observed.

CONCLUSIONS

Our results suggested that TP53 G:C > T:A transversion mutation in Rwandan patients with gastric cancer is more frequent than in non-Rwandan patients with gastric cancer, indicating at an alternative etiological and carcinogenic progression of gastric cancer in Rwanda.

Cellular Repair of Synthetic Analogs of Oxidative DNA Damage Reveals a Key Structure-Activity Relationship of the Cancer-Associated MUTYH DNA Repair Glycosylase

The base excision repair glycosylase MUTYH prevents mutations associated with the oxidatively damaged base, 8-oxo-7,8-dihydroguanine (OG), by removing undamaged misincorporated adenines from OG:A mispairs. Defects in OG:A repair in individuals with inherited MUTYH variants are correlated with the colorectal cancer predisposition syndrome known as MUTYH-associated polyposis (MAP). Herein, we reveal key structural features of OG required for efficient repair by human MUTYH using structure-activity relationships (SAR). We developed a GFP-based plasmid reporter assay to define SAR with synthetically generated OG analogs in human cell lines. Cellular repair results were compared with kinetic parameters measured by adenine glycosylase assays in vitro. Our results show substrates lacking the 2-amino group of OG, such as 8OI:A (8OI = 8-oxoinosine), are not repaired in cells, despite being excellent substrates in in vitro adenine glycosylase assays, new evidence that the search and detection steps are critical factors in cellular MUTYH repair functionality. Surprisingly, modification of the O8/N7H of OG, which is the distinguishing feature of OG relative to G, was tolerated in both MUTYH-mediated cellular repair and in vitro adenine glycosylase activity. The lack of sensitivity to alterations at the O8/N7H in the SAR of MUTYH substrates is distinct from previous work with bacterial MutY, indicating that the human enzyme is much less stringent in its lesion verification. Our results imply that the human protein relies almost exclusively on detection of the unique major groove position of the 2-amino group of OG within OGsyn:Aanti mispairs to select contextually incorrect adenines for excision and thereby thwart mutagenesis. These results predict that MUTYH variants that exhibit deficiencies in OG:A detection will be severely compromised in a cellular setting. Moreover, the reliance of MUTYH on the interaction with the OG 2-amino group suggests that disrupting this interaction with small molecules may provide a strategy to develop potent and selective MUTYH inhibitors.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Molecular-guided therapy for the treatment of patients with relapsed and refractory childhood cancers: a Beat Childhood Cancer Research Consortium trial

BACKGROUND

Children with relapsed central nervous system (CNS tumors), neuroblastoma, sarcomas, and other rare solid tumors face poor outcomes. This prospective clinical trial examined the feasibility of combining genomic and transcriptomic profiling of tumor samples with a molecular tumor board (MTB) approach to make real‑time treatment decisions for children with relapsed/refractory solid tumors.

METHODS

Subjects were divided into three strata: stratum 1-relapsed/refractory neuroblastoma; stratum 2-relapsed/refractory CNS tumors; and stratum 3-relapsed/refractory rare solid tumors. Tumor samples were sent for tumor/normal whole-exome (WES) and tumor whole-transcriptome (WTS) sequencing, and the genomic data were used in a multi-institutional MTB to make real‑time treatment decisions. The MTB recommended plan allowed for a combination of up to 4 agents. Feasibility was measured by time to completion of genomic sequencing, MTB review and initiation of treatment. Response was assessed after every two cycles using Response Evaluation Criteria in Solid Tumors (RECIST). Patient clinical benefit was calculated by the sum of the CR, PR, SD, and NED subjects divided by the sum of complete response (CR), partial response (PR), stable disease (SD), no evidence of disease (NED), and progressive disease (PD) subjects. Grade 3 and higher related and unexpected adverse events (AEs) were tabulated for safety evaluation.

RESULTS

A total of 186 eligible patients were enrolled with 144 evaluable for safety and 124 evaluable for response. The average number of days from biopsy to initiation of the MTB-recommended combination therapy was 38 days. Patient benefit was exhibited in 65% of all subjects, 67% of neuroblastoma subjects, 73% of CNS tumor subjects, and 60% of rare tumor subjects. There was little associated toxicity above that expected for the MGT drugs used during this trial, suggestive of the safety of utilizing this method of selecting combination targeted therapy.

CONCLUSIONS

This trial demonstrated the feasibility, safety, and efficacy of a comprehensive sequencing model to guide personalized therapy for patients with any relapsed/refractory solid malignancy. Personalized therapy was well tolerated, and the clinical benefit rate of 65% in these heavily pretreated populations suggests that this treatment strategy could be an effective option for relapsed and refractory pediatric cancers.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT02162732. Prospectively registered on June 11, 2014.

Oxidative stress accelerates intestinal tumorigenesis by enhancing 8-oxoguanine-mediated mutagenesis in MUTYH-deficient mice

Oxidative stress-induced DNA damage and its repair systems are related to cancer etiology; however, the molecular basis triggering tumorigenesis is not well understood. Here, we aimed to explore the causal relationship between oxidative stress, somatic mutations in pre-tumor-initiated normal tissues, and tumor incidence in the small intestines of MUTYH-proficient and MUTYH-deficient mice. MUTYH is a base excision repair enzyme associated with human colorectal cancer. Mice were administered different concentrations of potassium bromate (KBrO3; an oxidizing agent)-containing water for 4 wk for mutagenesis studies or 16 wk for tumorigenesis studies. All Mutyh -/- mice treated with >0.1% KBrO3 developed multiple tumors, and the average tumor number increased dose dependently. Somatic mutation analysis of Mutyh -/-/rpsL transgenic mice revealed that G:C  > T:A transversion was the only mutation type correlated positively with KBrO3 dose and tumor incidence. These mutations preferentially occurred at 5'G in GG and GAA sequences in rpsL This characteristic mutation pattern was also observed in the genomic region of Mutyh -/- tumors using whole-exome sequencing. It closely corresponded to signature 18 and SBS36, typically caused by 8-oxo-guanine (8-oxoG). 8-oxoG-induced mutations were sequence context dependent, yielding a biased amino acid change leading to missense and stop-gain mutations. These mutations frequently occurred in critical amino acid codons of known cancer drivers, Apc or Ctnnb1, known for activating Wnt signal pathway. Our results indicate that oxidative stress contributes to increased tumor incidence by elevating the likelihood of gaining driver mutations by increasing 8-oxoG-mediated mutagenesis, particularly under MUTYH-deficient conditions.

Biochemical analysis of H2O2-induced mutation spectra revealed that multiple damages were involved in the mutational process

Reactive oxygen species (ROS) are a major threat to genomic integrity and believed to be one of the etiologies of cancers. Here we developed a cell-free system to analyze ROS-induced mutagenesis, in which DNA was exposed to H2O2 and then subjected to translesion DNA synthesis by various DNA polymerases. Then, frequencies of mutations on the DNA products were determined by using next-generation sequencing technology. The majority of observed mutations were either C>A or G>A, caused by dAMP insertion at G and C residues, respectively. These mutations showed similar spectra to COSMIC cancer mutational signature 18 and 36, which are proposed to be caused by ROS. The in vitro mutations can be produced by replicative DNA polymerases (yeast DNA polymerase δ and ε), suggesting that ordinary DNA replication is sufficient to produce them. Very little G>A mutation was observed immediately after exposure to H2O2, but the frequency was increased during the 24 h after the ROS was removed, indicating that the initial oxidation product of cytosine needs to be maturated into a mutagenic lesion. Glycosylase-sensitivities of these mutations suggest that the C>A were made on 8-oxoguanine or Fapy-guanine, and that G>A were most likely made on 5-hydroxycytosine modification.

Direct Measurement of 8OG syn-anti Flips in Mutagenic 8OG•A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter the propensities to form lowly-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn•Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly-populated (pop. ~ 5%) and short-lived (lifetime ~ 50 ms) non-mutagenic 8OGanti•Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on a labeled partner adenine. The Watson-Crick-like 8OGsyn•Aanti mismatch also rescued the kinetics of Hoogsteen motions at distance A-T base pairs, which the G•A mismatch had slowed down. The results lend further support for 8OGanti•Aanti as a minor conformational state of 8OG•A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

A maternal germline mutator phenotype in a family affected by heritable colorectal cancer

Variation in DNA repair genes can increase cancer risk by elevating the rate of oncogenic mutation. Defects in one such gene, MUTYH, are known to elevate the incidence of colorectal cancer in a recessive Mendelian manner, and some evidence has also linked MUTYH to elevated incidence of other cancers as well as elevated mutation rates in normal somatic and germline cells. Here, we use whole genome sequencing to measure germline de novo mutation rates in a large extended family affected by pathogenic MUTYH variation and a history of colorectal cancer. Although this family's genotype, p.Y179C/V234M (c.536A>G/700G>A on transcript NM_001128425), contains a variant with conflicting functional interpretations, we use an in vitro cell line assay to determine that it partially attenuates MUTYH's function. In the children of mothers affected by the Y179C/V234M genotype, we identify an elevation of the C>A mutation rate that is weaker than mutator effects previously reported to be caused by other pathogenic MUTYH genotypes, suggesting that mutation rates in normal tissues may be useful for classifying cancer-associated variation along a continuum of severity. Surprisingly, we detect no significant elevation of the C>A mutation rate in children born to a father with the same biallelic MUTYH genotype, despite calculating that we should have adequate power to detect such a mutator effect. This suggests that the oxidative stress repaired by MUTYH may contribute more to female reproductive aging than male reproductive aging in the general population.

Epistasis between mutator alleles contributes to germline mutation spectra variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair [1], mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs [2,3]. In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh [4]. Its effect depended on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci had greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Germline de novo mutations in families with Mendelian cancer syndromes caused by defects in DNA repair

DNA repair defects underlie many cancer syndromes. We tested whether de novo germline mutations (DNMs) are increased in families with germline defects in polymerase proofreading or base excision repair. A parent with a single germline POLE or POLD1 mutation, or biallelic MUTYH mutations, had 3-4 fold increased DNMs over sex-matched controls. POLE had the largest effect. The DNMs carried mutational signatures of the appropriate DNA repair deficiency. No DNM increase occurred in offspring of MUTYH heterozygous parents. Parental DNA repair defects caused about 20-150 DNMs per child, additional to the ~60 found in controls, but almost all extra DNMs occurred in non-coding regions. No increase in post-zygotic mutations was detected, excepting a child with bi-allelic MUTYH mutations who was excluded from the main analysis; she had received chemotherapy and may have undergone oligoclonal haematopoiesis. Inherited DNA repair defects associated with base pair-level mutations increase DNMs, but phenotypic consequences appear unlikely.

DNA Glycosylases Define the Outcome of Endogenous Base Modifications

Chemically modified nucleic acid bases are sources of genomic instability and mutations but may also regulate gene expression as epigenetic or epitranscriptomic modifications. Depending on the cellular context, they can have vastly diverse impacts on cells, from mutagenesis or cytotoxicity to changing cell fate by regulating chromatin organisation and gene expression. Identical chemical modifications exerting different functions pose a challenge for the cell's DNA repair machinery, as it needs to accurately distinguish between epigenetic marks and DNA damage to ensure proper repair and maintenance of (epi)genomic integrity. The specificity and selectivity of the recognition of these modified bases relies on DNA glycosylases, which acts as DNA damage, or more correctly, as modified bases sensors for the base excision repair (BER) pathway. Here, we will illustrate this duality by summarizing the role of uracil-DNA glycosylases, with particular attention to SMUG1, in the regulation of the epigenetic landscape as active regulators of gene expression and chromatin remodelling. We will also describe how epigenetic marks, with a special focus on 5-hydroxymethyluracil, can affect the damage susceptibility of nucleic acids and conversely how DNA damage can induce changes in the epigenetic landscape by altering the pattern of DNA methylation and chromatin structure.

Neuroblastoma arises in early fetal development and its evolutionary duration predicts outcome

Neuroblastoma, the most frequent solid tumor in infants, shows very diverse outcomes from spontaneous regression to fatal disease. When these different tumors originate and how they evolve are not known. Here we quantify the somatic evolution of neuroblastoma by deep whole-genome sequencing, molecular clock analysis and population-genetic modeling in a comprehensive cohort covering all subtypes. We find that tumors across the entire clinical spectrum begin to develop via aberrant mitoses as early as the first trimester of pregnancy. Neuroblastomas with favorable prognosis expand clonally after short evolution, whereas aggressive neuroblastomas show prolonged evolution during which they acquire telomere maintenance mechanisms. The initial aneuploidization events condition subsequent evolution, with aggressive neuroblastoma exhibiting early genomic instability. We find in the discovery cohort (n = 100), and validate in an independent cohort (n = 86), that the duration of evolution is an accurate predictor of outcome. Thus, insight into neuroblastoma evolution may prospectively guide treatment decisions.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

Genetic Predisposition to Colorectal Cancer: How Many and Which Genes to Test?

Colorectal cancer is one of the most common tumors, and genetic predisposition is one of the key risk factors in the development of this malignancy. Lynch syndrome and familial adenomatous polyposis are the best-known genetic diseases associated with hereditary colorectal cancer. However, some other genetic disorders confer an increased risk of colorectal cancer, such as Li-Fraumeni syndrome (TP53 gene), MUTYH-associated polyposis (MUTYH gene), Peutz-Jeghers syndrome (STK11 gene), Cowden syndrome (PTEN gene), and juvenile polyposis syndrome (BMPR1A and SMAD4 genes). Moreover, the recent advances in molecular techniques, in particular Next-Generation Sequencing, have led to the identification of many new genes involved in the predisposition to colorectal cancers, such as RPS20, POLE, POLD1, AXIN2, NTHL1, MSH3, RNF43 and GREM1. In this review, we summarized the past and more recent findings in the field of cancer predisposition genes, with insights into the role of the encoded proteins and into the associated genetic disorders. Furthermore, we discussed the possible clinical utility of genetic testing in terms of prevention protocols and therapeutic approaches.

DNA damage and somatic mutations in mammalian cells after irradiation with a nail polish dryer

Ultraviolet A light is commonly emitted by UV-nail polish dryers with recent reports suggesting that long-term use may increase the risk for developing skin cancer. However, no experimental evaluation has been conducted to reveal the effect of radiation emitted by UV-nail polish dryers on mammalian cells. Here, we show that irradiation by a UV-nail polish dryer causes high levels of reactive oxygen species, consistent with 8-oxo-7,8-dihydroguanine damage and mitochondrial dysfunction. Analysis of somatic mutations reveals a dose-dependent increase of C:G>A:T substitutions in irradiated samples with mutagenic patterns similar to mutational signatures previously attributed to reactive oxygen species. In summary, this study demonstrates that radiation emitted by UV-nail polish dryers can both damage DNA and permanently engrave mutations on the genomes of primary mouse embryonic fibroblasts, human foreskin fibroblasts, and human epidermal keratinocytes.

Unravelling the instability of mutational signatures extraction via archetypal analysis

The high cosine similarity between some single-base substitution mutational signatures and their characteristic flat profiles could suggest the presence of overfitting and mathematical artefacts. The newest version (v3.3) of the signature database available in the Catalogue Of Somatic Mutations In Cancer (COSMIC) provides a collection of 79 mutational signatures, which has more than doubled with respect to previous version (30 profiles available in COSMIC signatures v2), making more critical the associations between signatures and specific mutagenic processes. This study both provides a systematic assessment of the de novo extraction task through simulation scenarios based on the latest version of the COSMIC signatures and highlights, through a novel approach using archetypal analysis, which COSMIC signatures are redundant and more likely to be considered as mathematical artefacts. 29 archetypes were able to reconstruct the profile of all the COSMIC signatures with cosine similarity > 0.8. Interestingly, these archetypes tend to group similar original signatures sharing either the same aetiology or similar biological processes. We believe that these findings will be useful to encourage the development of new de novo extraction methods avoiding the redundancy of information among the signatures while preserving the biological interpretation.

Game of clones: Battles in the field of carcinogenesis

Recent advances in bulk sequencing approaches as well as genomic decoding at the single-cell level have revealed surprisingly high somatic mutational burdens in normal tissues, as well as increased our understanding of the landscape of "field cancerization", that is, molecular and immune alterations in mutagen-exposed normal-appearing tissues that recapitulated those present in tumors. Charting the somatic mutational landscapes in normal tissues can have strong implications on our understanding of how tumors arise from mutagenized epithelium. Making sense of those mutations to understand the progression along the pathologic continuum of normal epithelia, preneoplasias, up to malignant tissues will help pave way for identification of ideal targets that can guide new strategies for preventing or eliminating cancers at their earliest stages of development. In this review, we will provide a brief history of field cancerization and its implications on understanding early stages of cancer pathogenesis and deviation from the pathologically "normal" state. The review will provide an overview of how mutations accumulating in normal tissues can lead to a patchwork of mutated cell clones that compete while maintaining an overall state of functional homeostasis. The review also explores the role of clonal competition in directing the fate of normal tissues and summarizes multiple mechanisms elicited in this phenomenon and which have been linked to cancer development. Finally, we highlight the importance of understanding mutations in normal tissues, as well as clonal competition dynamics (in both the epithelium and the microenvironment) and their significance in exploring new approaches to combatting cancer.

Incident Cancer Risk and Signatures Among Older MUTYH Carriers: Analysis of Population-Based and Genomic Cohorts

MUTYH carriers have an increased colorectal cancer risk in case-control studies, with loss of heterozygosity (LOH) as the presumed mechanism. We evaluated cancer risk among carriers in a prospective, population-based cohort of older adults. In addition, we assessed if cancers from carriers demonstrated mutational signatures (G:C>T:A transversions) associated with early LOH. We calculated incident risk of cancer and colorectal cancer among 13,131 sequenced study participants of the ASPirin in Reducing Events in the Elderly cohort, stratified by sex and adjusting for age, smoking, alcohol use, BMI, polyp history, history of cancer, and aspirin use. MUTYH carriers were identified among 13,033 participants in The Cancer Genome Atlas and International Cancer Genome Consortium, and somatic signatures of cancers were analyzed. Male MUTYH carriers demonstrated an increased risk for overall cancer incidence [multivariable HR, 1.66; 95% confidence interval (CI), 1.03-2.68; P = 0.038] driven by increased colorectal cancer incidence (multivariable HR, 3.55; 95% CI, 1.42-8.78; P = 0.007), as opposed to extracolonic cancer incidence (multivariable HR, 1.40; 95% CI, 0.81-2.44; P = 0.229). Female carriers did not demonstrate increased risk of cancer, colorectal cancer, or extracolonic cancers. Analysis of mutation signatures from cancers of MUTYH carriers revealed no significant contribution toward early mutagenesis from widespread G:C>T:A transversions among gastrointestinal epithelial cancers. Among cancers from carriers, somatic transversions associated with base-excision repair deficiency are uncommon, suggestive of diverse mechanisms of carcinogenesis in carriers compared with those who inherit biallelic MUTYH mutations.

PREVENTION RELEVANCE

Despite absence of loss of heterozygosity in colorectal cancers, elderly male MUTYH carriers appeared to be at increased of colorectal cancer.

The Dynamics of Somatic Mutagenesis During Life in Humans

From conception to death, human cells accumulate somatic mutations in their genomes. These mutations can contribute to the development of cancer and non-malignant diseases and have also been associated with aging. Rapid technological developments in sequencing approaches in the last few years and their application to normal tissues have greatly advanced our knowledge about the accumulation of these mutations during healthy aging. Whole genome sequencing studies have revealed that there are significant differences in mutation burden and patterns across tissues, but also that the mutation rates within tissues are surprisingly constant during adult life. In contrast, recent lineage-tracing studies based on whole-genome sequencing have shown that the rate of mutation accumulation is strongly increased early in life before birth. These early mutations, which can be shared by many cells in the body, may have a large impact on development and the origin of somatic diseases. For example, cancer driver mutations can arise early in life, decades before the detection of the malignancy. Here, we review the recent insights in mutation accumulation and mutagenic processes in normal tissues. We compare mutagenesis early and later in life and discuss how mutation rates and patterns evolve during aging. Additionally, we outline the potential impact of these mutations on development, aging and disease.

Inherited MUTYH mutations cause elevated somatic mutation rates and distinctive mutational signatures in normal human cells

Cellular DNA damage caused by reactive oxygen species is repaired by the base excision repair (BER) pathway which includes the DNA glycosylase MUTYH. Inherited biallelic MUTYH mutations cause predisposition to colorectal adenomas and carcinoma. However, the mechanistic progression from germline MUTYH mutations to MUTYH-Associated Polyposis (MAP) is incompletely understood. Here, we sequence normal tissue DNAs from 10 individuals with MAP. Somatic base substitution mutation rates in intestinal epithelial cells were elevated 2 to 4-fold in all individuals, except for one showing a 31-fold increase, and were also increased in other tissues. The increased mutation burdens were of multiple mutational signatures characterised by C > A changes. Different mutation rates and signatures between individuals are likely due to different MUTYH mutations or additional inherited mutations in other BER pathway genes. The elevated base substitution rate in normal cells likely accounts for the predisposition to neoplasia in MAP. Despite ubiquitously elevated mutation rates, individuals with MAP do not display overt evidence of premature ageing. Thus, accumulation of somatic mutations may not be sufficient to cause the global organismal functional decline of ageing.

The impact of rare germline variants on human somatic mutation processes

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Identifying colorectal cancer caused by biallelic MUTYH pathogenic variants using tumor mutational signatures

Carriers of germline biallelic pathogenic variants in the MUTYH gene have a high risk of colorectal cancer. We test 5649 colorectal cancers to evaluate the discriminatory potential of a tumor mutational signature specific to MUTYH for identifying biallelic carriers and classifying variants of uncertain clinical significance (VUS). Using a tumor and matched germline targeted multi-gene panel approach, our classifier identifies all biallelic MUTYH carriers and all known non-carriers in an independent test set of 3019 colorectal cancers (accuracy = 100% (95% confidence interval 99.87-100%)). All monoallelic MUTYH carriers are classified with the non-MUTYH carriers. The classifier provides evidence for a pathogenic classification for two VUS and a benign classification for five VUS. Somatic hotspot mutations KRAS p.G12C and PIK3CA p.Q546K are associated with colorectal cancers from biallelic MUTYH carriers compared with non-carriers (p = 2 × 10-23 and p = 6 × 10-11, respectively). Here, we demonstrate the potential application of mutational signatures to tumor sequencing workflows to improve the identification of biallelic MUTYH carriers.

Concordance of hydrogen peroxide-induced 8-oxo-guanine patterns with two cancer mutation signatures of upper GI tract tumors

Oxidative DNA damage has been linked to inflammation, cancer, and aging. Here, we have mapped two types of oxidative DNA damage, oxidized guanines produced by hydrogen peroxide and oxidized thymines created by potassium permanganate, at a single-base resolution. 8-Oxo-guanine occurs strictly dependent on the G/C sequence context and shows a pronounced peak at transcription start sites (TSSs). We determined the trinucleotide sequence pattern of guanine oxidation. This pattern shows high similarity to the cancer-associated single-base substitution signatures SBS18 and SBS36. SBS36 is found in colorectal cancers that carry mutations in MUTYH, encoding a repair enzyme that operates on 8-oxo-guanine mispairs. SBS18 is common in inflammation-associated upper gastrointestinal tract tumors including esophageal and gastric adenocarcinomas. Oxidized thymines induced by permanganate occur with a distinct dinucleotide specificity, 5'T-A/C, and are depleted at the TSS. Our data suggest that two cancer mutational signatures, SBS18 and SBS36, are caused by reactive oxygen species.

Mutational signatures are markers of drug sensitivity of cancer cells

Genomic analyses have revealed mutational footprints associated with DNA maintenance gone awry, or with mutagen exposures. Because cancer therapeutics often target DNA synthesis or repair, we asked if mutational signatures make useful markers of drug sensitivity. We detect mutational signatures in cancer cell line exomes (where matched healthy tissues are not available) by adjusting for the confounding germline mutation spectra across ancestries. We identify robust associations between various mutational signatures and drug activity across cancer cell lines; these are as numerous as associations with established genetic markers such as driver gene alterations. Signatures of prior exposures to DNA damaging agents - including chemotherapy - tend to associate with drug resistance, while signatures of deficiencies in DNA repair tend to predict sensitivity towards particular therapeutics. Replication analyses across independent drug and CRISPR genetic screening data sets reveal hundreds of robust associations, which are provided as a resource for drug repurposing guided by mutational signature markers.

A natural mutator allele shapes mutation spectrum variation in mice

Although germline mutation rates and spectra can vary within and between species, common genetic modifiers of the mutation rate have not been identified in nature¹. Here we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory strains C57BL/6J (B haplotype) and DBA/2J (D haplotype). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes². We show that mice inheriting D haplotypes at a quantitative trait locus on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily owing to the activity of a C>A-dominated mutational signature known as SBS18. The B and D quantitative trait locus haplotypes encode different alleles of Mutyh, a DNA repair gene that underlies the heritable cancer predisposition syndrome that causes colorectal tumors with a high SBS18 mutation load^3,4. Both B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.

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.

Detection and Localization of Solid Tumors Utilizing the Cancer-Type-Specific Mutational Signatures

Accurate detection and location of tumor lesions are essential for improving the diagnosis and personalized cancer therapy. However, the diagnosis of lesions with fuzzy histology is mainly dependent on experiences and with low accuracy and efficiency. Here, we developed a logistic regression model based on mutational signatures (MS) for each cancer type to trace the tumor origin. We observed MS could distinguish cancer from inflammation and healthy individuals. By collecting extensive datasets of samples from ten tumor types in the training cohort (5,001 samples) and independent testing cohort (2,580 samples), cancer-type-specific MS patterns (CTS-MS) were identified and had a robust performance in distinguishing different types of primary and metastatic solid tumors (AUC:0.76 ∼ 0.93). Moreover, we validated our model in an Asian population and found that the AUC of our model in predicting the tumor origin of the Asian population was higher than 0.7. The metastatic tumor lesions inherited the MS pattern of the primary tumor, suggesting the capability of MS in identifying the tissue-of-origin for metastatic cancers. Furthermore, we distinguished breast cancer and prostate cancer with 90% accuracy by combining somatic mutations and CTS-MS from cfDNA, indicating that the CTS-MS could improve the accuracy of cancer-type prediction by cfDNA. In summary, our study demonstrated that MS was a novel reliable biomarker for diagnosing solid tumors and provided new insights into predicting tissue-of-origin.

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

Whole-genome sequencing (WGS) permits comprehensive cancer genome analyses, revealing mutational signatures, imprints of DNA damage and repair processes that have arisen in each patient's cancer. We performed mutational signature analyses on 12,222 WGS tumor-normal matched pairs, from patients recruited via the UK National Health Service. We contrasted our results to two independent cancer WGS datasets, the International Cancer Genome Consortium (ICGC) and Hartwig Foundation, involving 18,640 WGS cancers in total. Our analyses add 40 single and 18 double substitution signatures to the current mutational signature tally. Critically, we show for each organ, that cancers have a limited number of 'common' signatures and a long tail of 'rare' signatures. We provide a practical solution for utilizing this concept of common versus rare signatures in future analyses.

A pilot study of genomic-guided induction therapy followed by immunotherapy with difluoromethylornithine maintenance for high-risk neuroblastoma

BACKGROUND

Survival for patients with high-risk neuroblastoma (HRNB) remains poor despite aggressive multimodal therapies.

AIMS

To study the feasibility and safety of incorporating a genomic-based targeted agent to induction therapy for HRNB as well as the feasibility and safety of adding difluoromethylornithine (DFMO) to anti-GD2 immunotherapy.

METHODS

Twenty newly diagnosed HRNB patients were treated on this multicenter pilot trial. Molecular tumor boards selected one of six targeted agents based on tumor-normal whole exome sequencing and tumor RNA-sequencing results. Treatment followed standard upfront HRNB chemotherapy with the addition of the selected targeted agent to cycles 3-6 of induction. Following consolidation, DFMO (750 mg/m² twice daily) was added to maintenance with dinutuximab and isotretinoin, followed by continuation of DFMO alone for 2 years. DNA methylation analysis was performed retrospectively and compared to RNA expression.

RESULTS

Of the 20 subjects enrolled, 19 started targeted therapy during cycle 3 and 1 started during cycle 5. Eighty-five percent of subjects met feasibility criteria (receiving 75% of targeted agent doses). Addition of targeted agents did not result in toxicities requiring dose reduction of chemotherapy or permanent discontinuation of targeted agent. Following standard consolidation, 15 subjects continued onto immunotherapy with DFMO. This combination was well-tolerated and resulted in no unexpected adverse events related to DFMO.

CONCLUSION

This study demonstrates the safety and feasibility of adding targeted agents to standard induction therapy and adding DFMO to immunotherapy for HRNB. This treatment regimen has been expanded to a Phase II trial to evaluate efficacy.

Monoallelic deleterious MUTYH germline variants as a driver for tumorigenesis

MUTYH encodes a glycosylase involved in the base excision repair of DNA. Biallelic pathogenic germline variants in MUTYH cause an autosomal recessive condition known as MUTYH-associated adenomatous polyposis and consequently increase the risk of colorectal cancer. However, reports of increased cancer risk in individuals carrying only one defective MUTYH allele are controversial and based on studies involving few individuals. Here, we describe a comprehensive investigation of monoallelic pathogenic MUTYH germline variants in 10,389 cancer patients across 33 different tumour types and 117,000 healthy individuals. Our results indicate that monoallelic pathogenic MUTYH germline variants can lead to tumorigenesis through a mechanism of somatic loss of heterozygosity of the functional MUTYH allele in the tumour. We confirmed that the frequency of monoallelic pathogenic MUTYH germline variants is higher in individuals with cancer than in the general population, although this frequency is not homogeneous among tumour types. We also demonstrated that the MUTYH mutational signature is present only in tumours with loss of the functional allele and found that the characteristic MUTYH base substitution (C>A) increases stop-codon generation. We identified key genes that are affected during tumorigenesis. In conclusion, we propose that carriers of the monoallelic pathogenic MUTYH germline variant are at a higher risk of developing tumours, especially those with frequent loss of heterozygosity events, such as adrenal adenocarcinoma, although the overall risk is still low. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The context-specific role of germline pathogenicity in tumorigenesis

Human cancers arise from environmental, heritable and somatic factors, but how these mechanisms interact in tumorigenesis is poorly understood. Studying 17,152 prospectively sequenced patients with cancer, we identified pathogenic germline variants in cancer predisposition genes, and assessed their zygosity and co-occurring somatic alterations in the concomitant tumors. Two major routes to tumorigenesis were apparent. In carriers of pathogenic germline variants in high-penetrance genes (5.1% overall), lineage-dependent patterns of biallelic inactivation led to tumors exhibiting mechanism-specific somatic phenotypes and fewer additional somatic oncogenic drivers. Nevertheless, 27% of cancers in these patients, and most tumors in patients with pathogenic germline variants in lower-penetrance genes, lacked particular hallmarks of tumorigenesis associated with the germline allele. The dependence of tumors on pathogenic germline variants is variable and often dictated by both penetrance and lineage, a finding with implications for clinical management.

Evaluating the utility of tumour mutational signatures for identifying hereditary colorectal cancer and polyposis syndrome carriers

OBJECTIVE

Germline pathogenic variants (PVs) in the DNA mismatch repair (MMR) genes and in the base excision repair gene MUTYH underlie hereditary colorectal cancer (CRC) and polyposis syndromes. We evaluated the robustness and discriminatory potential of tumour mutational signatures in CRCs for identifying germline PV carriers.

DESIGN

Whole-exome sequencing of formalin-fixed paraffin-embedded (FFPE) CRC tissue was performed on 33 MMR germline PV carriers, 12 biallelic MUTYH germline PV carriers, 25 sporadic MLH1 methylated MMR-deficient CRCs (MMRd controls) and 160 sporadic MMR-proficient CRCs (MMRp controls) and included 498 TCGA CRC tumours. COSMIC V3 single base substitution (SBS) and indel (ID) mutational signatures were assessed for their ability to differentiate CRCs that developed in carriers from non-carriers.

RESULTS

The combination of mutational signatures SBS18 and SBS36 contributing >30% of a CRC's signature profile was able to discriminate biallelic MUTYH carriers from all other non-carrier control CRCs with 100% accuracy (area under the curve (AUC) 1.0). SBS18 and SBS36 were associated with specific MUTYH variants p.Gly396Asp (p=0.025) and p.Tyr179Cys (p=5×10^-5), respectively. The combination of ID2 and ID7 could discriminate the 33 MMR PV carrier CRCs from the MMRp control CRCs (AUC 0.99); however, SBS and ID signatures, alone or in combination, could not provide complete discrimination (AUC 0.79) between CRCs from MMR PV carriers and sporadic MMRd controls.

CONCLUSION

Assessment of SBS and ID signatures can discriminate CRCs from biallelic MUTYH carriers and MMR PV carriers from non-carriers with high accuracy, demonstrating utility as a potential diagnostic and variant classification tool.

Novel temporal and spatial patterns of metastatic colonization from breast cancer rapid-autopsy tumor biopsies

BACKGROUND

Metastatic breast cancer is a deadly disease with a low 5-year survival rate. Tracking metastatic spread in living patients is difficult and thus poorly understood.

METHODS

Via rapid autopsy, we have collected 30 tumor samples over 3 timepoints and across 8 organs from a triple-negative metastatic breast cancer patient. The large number of sites sampled, together with deep whole-genome sequencing and advanced computational analysis, allowed us to comprehensively reconstruct the tumor's evolution at subclonal resolution.

RESULTS

The most unique, previously unreported aspect of the tumor's evolution that we observed in this patient was the presence of "subclone incubators," defined as metastatic sites where substantial tumor evolution occurs before colonization of additional sites and organs by subclones that initially evolved at the incubator site. Overall, we identified four discrete waves of metastatic expansions, each of which resulted in a number of new, genetically similar metastasis sites that also enriched for particular organs (e.g., abdominal vs bone and brain). The lung played a critical role in facilitating metastatic spread in this patient: the lung was the first site of metastatic escape from the primary breast lesion, subclones at this site were likely the source of all four subsequent metastatic waves, and multiple sites in the lung acted as subclone incubators. Finally, functional annotation revealed that many known drivers or metastasis-promoting tumor mutations in this patient were shared by some, but not all metastatic sites, highlighting the need for more comprehensive surveys of a patient's metastases for effective clinical intervention.

CONCLUSIONS

Our analysis revealed the presence of substantial tumor evolution at metastatic incubator sites in a patient, with potentially important clinical implications. Our study demonstrated that sampling of a large number of metastatic sites affords unprecedented detail for studying metastatic evolution.

Mutation signatures of pediatric acute myeloid leukemia and normal blood progenitors associated with differential patient outcomes

A subset of pediatric AML cases harbors more somatic mutations in their genomes compared to normal blood progenitors. This subset displays expression profiles that resemble more committed progenitors and associates with better patient survival.

Acquisition of oncogenic mutations with age is believed to be rate limiting for carcinogenesis. However, the incidence of leukemia in children is higher than in young adults. Here we compare somatic mutations across pediatric acute myeloid leukemia (pAML) patient-matched leukemic blasts and hematopoietic stem and progenitor cells (HSPC), as well as HSPCs from age-matched healthy donors. HSPCs in the leukemic bone marrow have limited genetic relatedness and share few somatic mutations with the cell of origin of the malignant blasts, suggesting polyclonal hematopoiesis in patients with pAML. Compared with normal HSPCs, a subset of pAML cases harbored more somatic mutations and a distinct composition of mutational process signatures. We hypothesize that these cases might have arisen from a more committed progenitor. This subset had better outcomes than pAML cases with mutation burden comparable with age-matched healthy HSPCs. Our study provides insights into the etiology and patient stratification of pAML.

POLE, POLD1, and NTHL1: the last but not the least hereditary cancer-predisposing genes

POLE, POLD1, and NTHL1 are involved in DNA replication and have recently been recognized as hereditary cancer-predisposing genes, because their alterations are associated with colorectal cancer and other tumors. POLE/POLD1-associated syndrome shows an autosomal dominant inheritance, whereas NTHL1-associated syndrome follows an autosomal recessive pattern. Although the prevalence of germline monoallelic POLE/POLD1 and biallelic NTHL1 pathogenic variants is low, they determine different phenotypes with a broad tumor spectrum overlapping that of other hereditary conditions like Lynch Syndrome or Familial Adenomatous Polyposis. Endometrial and breast cancers, and probably ovarian and brain tumors are also associated with POLE/POLD1 alterations, while breast cancer and other unusual tumors are correlated with NTHL1 pathogenic variants. POLE-mutated colorectal and endometrial cancers are associated with better prognosis and may show favorable responses to immunotherapy. Since POLE/POLD1-mutated tumors show a high tumor mutational burden producing an increase in neoantigens, the identification of POLE/POLD1 alterations could help select patients suitable for immunotherapy treatment. In this review, we will investigate the role of POLE, POLD1, and NTHL1 genetic variants in cancer predisposition, discussing the potential future therapeutic applications and assessing the utility of performing a routine genetic testing for these genes, in order to implement prevention and surveillance strategies in mutation carriers.

Mutational signatures: emerging concepts, caveats and clinical applications

Whole-genome sequencing has brought the cancer genomics community into new territory. Thanks to the sheer power provided by the thousands of mutations present in each patient's cancer, we have been able to discern generic patterns of mutations, termed 'mutational signatures', that arise during tumorigenesis. These mutational signatures provide new insights into the causes of individual cancers, revealing both endogenous and exogenous factors that have influenced cancer development. This Review brings readers up to date in a field that is expanding in computational, experimental and clinical directions. We focus on recent conceptual advances, underscoring some of the caveats associated with using the mutational signature frameworks and highlighting the latest experimental insights. We conclude by bringing attention to areas that are likely to see advancements in clinical applications.

From APC to the genetics of hereditary and familial colon cancer syndromes

Hereditary colorectal cancer (CRC) syndromes attributable to high penetrance mutations represent 9-26% of young-onset CRC cases. The clinical significance of many of these mutations is understood well enough to be used in diagnostics and as an aid in patient care. However, despite the advances made in the field, a significant proportion of familial and early-onset cases remains molecularly uncharacterized and extensive work is still needed to fully understand the genetic nature of CRC susceptibility. With the emergence of next-generation sequencing and associated methods, several predisposition loci have been unraveled, but validation is incomplete. Individuals with cancer-predisposing mutations are currently enrolled in life-long surveillance, but with the development of new treatments, such as cancer vaccinations, this might change in the not so distant future for at least some individuals. For individuals without a known cause for their disease susceptibility, prevention and therapy options are less precise. Herein, we review the progress achieved in the last three decades with a focus on how CRC predisposition genes were discovered. Furthermore, we discuss the clinical implications of these discoveries and anticipate what to expect in the next decade.

Genetic and epigenetic basis of hepatoblastoma diversity

Hepatoblastoma (HB) is the most common pediatric liver malignancy; however, hereditary predisposition and acquired molecular aberrations related to HB clinicopathological diversity are not well understood. Here, we perform an integrative genomic profiling of 163 pediatric liver tumors (154 HBs and nine hepatocellular carcinomas) based on the data acquired from a cohort study (JPLT-2). The total number of somatic mutations is precious low (0.52/Mb on exonic regions) but correlated with age at diagnosis. Telomerase reverse transcriptase (TERT) promoter mutations are prevalent in the tween HBs, selective in the transitional liver cell tumor (TLCT, > 8 years old). DNA methylation profiling reveals that classical HBs are characterized by the specific hypomethylated enhancers, which are enriched with binding sites for ASCL2, a regulatory transcription factor for definitive endoderm in Wnt-pathway. Prolonged upregulation of ASCL2, as well as fetal-liver-like methylation patterns of IGF2 promoters, suggests their "cell of origin" derived from the premature hepatoblast, similar to intestinal epithelial cells, which are highly proliferative. Systematic molecular profiling of HB is a promising approach for understanding the epigenetic drivers of hepatoblast carcinogenesis and deriving clues for risk stratification.

Defects in 8-oxo-guanine repair pathway cause high frequency of C > A substitutions in neuroblastoma

The collection of large amounts of whole-genome sequencing data allowed for identification of mutational signatures, which are characteristic combinations of substitutions in the context of neighboring bases. The clinical significance of these mutational signatures is still largely unknown. In neuroblastoma, we showed that high levels of cytosine > adenine (C > A) substitutions are associated with poor survival. We identified that these high levels of C > A substitutions result from defects in 8-oxo-guanine repair, specifically from copy number loss of the DNA glycosylases MUTYH and OGG1. The high frequency of C > A substitutions in neuroblastoma contributes to the increased adaptive capacity of these tumors. Thereby, we link basic molecular genetic mutation patterns to clinically significant tumor evolution processes.

Neuroblastomas are childhood tumors with frequent fatal relapses after induction treatment, which is related to tumor evolution with additional genomic events. Our whole-genome sequencing data analysis revealed a high frequency of somatic cytosine > adenine (C > A) substitutions in primary neuroblastoma tumors, which was associated with poor survival. We showed that increased levels of C > A substitutions correlate with copy number loss (CNL) of OGG1 or MUTYH. Both genes encode DNA glycosylases that recognize 8-oxo-guanine (8-oxoG) lesions as a first step of 8-oxoG repair. Tumor organoid models with CNL of OGG1 or MUTYH show increased 8-oxoG levels compared to wild-type cells. We used CRISPR-Cas9 genome editing to create knockout clones of MUTYH and OGG1 in neuroblastoma cells. Whole-genome sequencing of single-cell OGG1 and MUTYH knockout clones identified an increased accumulation of C > A substitutions. Mutational signature analysis of these OGG1 and MUTYH knockout clones revealed enrichment for C > A signatures 18 and 36, respectively. Clustering analysis showed that the knockout clones group together with tumors containing OGG1 or MUTYH CNL. In conclusion, we demonstrate that defects in 8-oxoG repair cause accumulation of C > A substitutions in neuroblastoma, which contributes to mutagenesis and tumor evolution.

Analyzing the Opportunities to Target DNA Double-Strand Breaks Repair and Replicative Stress Responses to Improve Therapeutic Index of Colorectal Cancer

Despite the ample improvements of CRC molecular landscape, the therapeutic options still rely on conventional chemotherapy-based regimens for early disease, and few targeted agents are recommended for clinical use in the metastatic setting. Moreover, the impact of cytotoxic, targeted agents, and immunotherapy combinations in the metastatic scenario is not fully satisfactory, especially the outcomes for patients who develop resistance to these treatments need to be improved. Here, we examine the opportunity to consider therapeutic agents targeting DNA repair and DNA replication stress response as strategies to exploit genetic or functional defects in the DNA damage response (DDR) pathways through synthetic lethal mechanisms, still not explored in CRC. These include the multiple actors involved in the repair of DNA double-strand breaks (DSBs) through homologous recombination (HR), classical non-homologous end joining (NHEJ), and microhomology-mediated end-joining (MMEJ), inhibitors of the base excision repair (BER) protein poly (ADP-ribose) polymerase (PARP), as well as inhibitors of the DNA damage kinases ataxia-telangiectasia and Rad3 related (ATR), CHK1, WEE1, and ataxia-telangiectasia mutated (ATM). We also review the biomarkers that guide the use of these agents, and current clinical trials with targeted DDR therapies.

A systematic CRISPR screen defines mutational mechanisms underpinning signatures caused by replication errors and endogenous DNA damage

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR-Cas9 knockouts (Δ) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage, and performed whole-genome sequencing of 173 subclones. ΔOGG1, ΔUNG, ΔEXO1, ΔRNF168, ΔMLH1, ΔMSH2, ΔMSH6, ΔPMS1, and ΔPMS2 produced marked mutational signatures indicative of being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-dG elimination is sequence-context-specific while uracil clearance is sequence-context-independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C>A transversions), differential misincorporation by replicative polymerases (T>C and C>T transitions), and we propose a 'reverse template slippage' model for T>A transversions. ΔMLH1, ΔMSH6, and ΔMSH2 signatures were similar to each other but distinct from ΔPMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 WGS cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies.

The diverse molecular profiles of lynch syndrome-associated colorectal cancers are (highly) dependent on underlying germline mismatch repair mutations

Lynch syndrome (LS) is a hereditary cancer syndrome that accounts for 3% of all new colorectal cancer (CRC) cases. Patients carry a germline pathogenic variant in one of the mismatch repair (MMR) genes (MLH1, MSH2, MSH6 or PMS2), which encode proteins involved in a post-replicative proofreading and editing mechanism. The clinical presentation of LS is highly heterogeneous, showing high variability in age at onset and penetrance of cancer, which may be partly attributable to the molecular profiles of carcinomas. This review discusses the frequency of alterations in the WNT/B-CATENIN, RAF/MEK/ERK and PI3K/PTEN/AKT pathways identified in all four LS subgroups and how these changes may relate to the 'three pathway model' of carcinogenesis, in which LS CRCs develop from MMR-proficient adenomas, MMR-deficient adenomas or directly from MMR-deficient crypts. Understanding the specific differences in carcinogenesis for each LS subgroup will aid in the further optimization of guidelines for diagnosis, surveillance and treatment.

Targeting DNA Damage Repair for Immune Checkpoint Inhibition: Mechanisms and Potential Clinical Applications

DNA damage repair (DDR) pathways play an essential role in maintaining genomic integrity. DDR dysfunction leads to accumulated DNA damage, predisposition to cancer, and high sensitivity to chemotherapy and radiotherapy. Recent studies have demonstrated that DDR status is associated with response to immune checkpoint inhibitors (ICIs). Among the DDR pathways, mismatch repair is one of the most recognized predictive biomarkers for ICIs. Furthermore, preclinical and early clinical studies suggest the rationale of combining agents targeting the DDR pathways, such as poly (ADP-ribose) polymerase (PARP) inhibitors, cyclin-dependent kinase 4/6 (CDK4/6) inhibitors, and ataxia telangiectasia and rad3-related (ATR) kinase inhibitors, with ICIs. In the present review, we describe the predictive role of DDR pathways in ICIs and summarize the advances in potential combination strategies of novel agents targeting DDR with ICIs for cancer treatment.

Integrated genomic analysis identifies driver genes and cisplatin-resistant progenitor phenotype in pediatric liver cancer

Pediatric liver cancers (PLCs) comprise diverse diseases affecting infants, children and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of hepatoblastoma patients, all before 3 years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between 'Hepatocytic', 'Liver Progenitor' and 'Mesenchymal' molecular subgroups of hepatoblastoma. We showed that during chemotherapy, 'Liver Progenitor' cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies.

Tumor Signature Analysis Implicates Hereditary Cancer Genes in Endometrial Cancer Development

Risk of endometrial cancer (EC) is increased ~2-fold for women with a family history of cancer, partly due to inherited pathogenic variants in mismatch repair (MMR) genes. We explored the role of additional genes as explanation for familial EC presentation by investigating germline and EC tumor sequence data from The Cancer Genome Atlas (n = 539; 308 European ancestry), and germline data from 33 suspected familial European ancestry EC patients demonstrating immunohistochemistry-detected tumor MMR proficiency. Germline variants in MMR and 26 other known/candidate EC risk genes were annotated for pathogenicity in the two EC datasets, and also for European ancestry individuals from gnomAD as a population reference set (n = 59,095). Ancestry-matched case-control comparisons of germline variant frequency and/or sequence data from suspected familial EC cases highlighted ATM, PALB2, RAD51C, MUTYH and NBN as candidates for large-scale risk association studies. Tumor mutational signature analysis identified a microsatellite-high signature for all cases with a germline pathogenic MMR gene variant. Signature analysis also indicated that germline loss-of-function variants in homologous recombination (BRCA1, PALB2, RAD51C) or base excision (NTHL1, MUTYH) repair genes can contribute to EC development in some individuals with germline variants in these genes. These findings have implications for expanded therapeutic options for EC cases.

Unraveling the genomic landscape of colorectal cancer through mutational signatures

Colorectal cancer, along with most other cancer types, is driven by somatic mutations. Characteristic patterns of somatic mutations, known as mutational signatures, arise as a result of the activities of different mutational processes. Mutational signatures have diverse origins, including exogenous and endogenous sources. In the case of colorectal cancer, the analysis of mutational signatures has elucidated specific signatures for classically associated DNA repair deficiencies, namely mismatch repair (leading to microsatellite instability), base excision repair (due to MUTYH or NTHL1 mutations), and polymerase proofreading (due to POLE and POLD1 exonuclease domain mutations). Additional signatures also play a role in colorectal cancer, including those related to normal aging and those associated with gut microbiota, as well as a number of signatures with unknown etiologies. This chapter provides an overview of the current knowledge of mutational signatures, with a focus on colorectal cancer and on the recently reported signatures in physiologically normal and inflammatory bowel disease-affected somatic colon tissues.