The genome as a record of environmental exposure

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.

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.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Malignant features of minipig melanomas prior to spontaneous regression

In MeLiM minipigs, melanomas develop around birth, can metastasize, and have histopathologic characteristics similar to humans. Interestingly, MeLiM melanomas eventually regress. This favorable outcome raises the question of their malignancy, which we investigated. We clinically followed tens of tumors from onset to first signs of regression. Transcriptome analysis revealed an enrichment of all cancer hallmarks in melanomas, although no activating or suppressing somatic mutation were found in common driver genes. Analysis of tumor cell genomes revealed high mutation rates without UV signature. Canonical proliferative, survival and angiogenic pathways were detected in MeLiM tumor cells all along progression stages. Functionally, we show that MeLiM melanoma cells are capable to grow in immunocompromised mice, with serial passages and for a longer time than in MeLiM pigs. Pigs set in place an immune response during progression with dense infiltration by myeloid cells while melanoma cells are deficient in B2M expression. To conclude, our data on MeLiM melanomas reveal several malignancy characteristics. The combination of these features with the successful spontaneous regression of these tumors make it an outstanding model to study an efficient anti-tumor immune response.

Lung cancer in patients who have never smoked - an emerging disease

Lung cancer is the most common cause of cancer-related deaths globally. Although smoking-related lung cancers continue to account for the majority of diagnoses, smoking rates have been decreasing for several decades. Lung cancer in individuals who have never smoked (LCINS) is estimated to be the fifth most common cause of cancer-related deaths worldwide in 2023, preferentially occurring in women and Asian populations. As smoking rates continue to decline, understanding the aetiology and features of this disease, which necessitate unique diagnostic and treatment paradigms, will be imperative. New data have provided important insights into the molecular and genomic characteristics of LCINS, which are distinct from those of smoking-associated lung cancers and directly affect treatment decisions and outcomes. Herein, we review the emerging data regarding the aetiology and features of LCINS, particularly the genetic and environmental underpinnings of this disease as well as their implications for treatment. In addition, we outline the unique diagnostic and therapeutic paradigms of LCINS and discuss future directions in identifying individuals at high risk of this disease for potential screening efforts.

The chemotherapeutic drug CX-5461 is a potent mutagen in cultured human cells

The chemotherapeutic agent CX-5461, or pidnarulex, has been fast-tracked by the United States Food and Drug Administration for early-stage clinical studies of BRCA1-, BRCA2- and PALB2-mutated cancers. It is under investigation in phase I and II trials. Here, we find that, although CX-5461 exhibits synthetic lethality in BRCA1-/BRCA2-deficient cells, it also causes extensive, nonselective, collateral mutagenesis in all three cell lines tested, to magnitudes that exceed known environmental carcinogens.

Molecular landscapes of oral cancers of unknown etiology

The incidence of the mobile tongue cancer in young patients has been rising. This oral cancer (OC) type has no identified risk factors (NIRF), no established molecular markers and is not yet recognized as a distinct clinical entity. To understand this emerging malignancy, we innovatively analyzed the public head and neck cancer multi-omics data. We identified mutational signatures that successfully stratified 307 OC and 109 laryngeal cancer cases according to their clinico-pathological characteristics. The NIRF OCs exhibited significantly increased activities of endogenous clock-like and APOBEC-associated mutagenesis, alongside specific cancer driver gene mutations, distinct methylome patterns and prominent antimicrobial transcriptomic responses. Furthermore, we show that mutational signature SBS16 in OCs reflects the combined effects of alcohol drinking and tobacco smoking. Our study characterizes the unique disease histories and molecular programs of the NIRF OCs revealing that this emerging cancer subtype is likely driven by increased endogenous mutagenesis correlated with responses to microbial insults.

Arsenic is a potent co-mutagen of ultraviolet light

Arsenic enhances the carcinogenicity of ultraviolet radiation (UVR). However, the mechanisms of arsenic-driven oncogenesis are not well understood. Here, we utilize experimental systems to investigate the carcinogenic and mutagenic properties of co-exposure to arsenic and UVR. In vitro and in vivo exposures indicate that, by itself, arsenic is not mutagenic. However, in combination with UVR, arsenic exposure has a synergistic effect leading to an accelerated mouse skin carcinogenesis and to more than 2-fold enrichment of UVR mutational burden. Notably, mutational signature ID13, previously found only in UVR-associated human skin cancers, is observed exclusively in mouse skin tumors and cell lines jointly exposed to arsenic and UVR. This signature was not observed in any model system exposed purely to arsenic or purely to UVR, making ID13, to the best of our knowledge, the first co-exposure signature to be reported using controlled experimental conditions. Analysis of existing skin cancer genomics data reveals that only a subset of cancers harbor ID13 and these exhibit an elevated UVR mutagenesis. Our results report a unique mutational signature caused by a co-exposure to two environmental carcinogens and provide comprehensive evidence that arsenic is a potent co-mutagen and co-carcinogen of UVR.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Characterization of highly active mutational signatures in tumors from a large Chinese population

The majority of mutational signatures have been characterized in tumors from Western countries and the degree to which mutational signatures are similar or different in Eastern populations has not been fully explored. We leveraged a large-scale clinical sequencing cohort of tumors from a Chinese population containing 25 tumor types and found that the highly active mutational signatures were similar to those previously characterized1,2. The aristolochic acid signature SBS22 was observed in four soft tissue sarcomas and the POLE-associated signature SBS10 was observed in a gallbladder carcinoma. In lung adenocarcinoma, the polycyclic aromatic hydrocarbon (PAH) signature SBS4 was significantly higher in males compared to females but not associated with smoking status. The UV-associated signature SBS7 was significantly lower in cutaneous melanomas from the Chinese population compared to a similar American cohort. Overall, these results add to our understanding of the mutational processes that contribute to tumors from the Chinese population.

Profiling endogenous, environmental, and infectious disease mutational signatures in blastic plasmacytoid dendritic cell neoplasms

BACKGROUND

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare hematopoietic disease derived from plasmacytoid dendritic lineage cells. The disease typically shows skin as well as frequent bone marrow and peripheral blood involvement. However, the pathogenesis of this disease is still not well understood. While somatic point mutations and genetic rearrangements have been described in BPDCN, the types and origins of these mutations and relationships to other cancer types is not well understood.

MATERIALS AND METHODS

To probe the origins of BPDCN, we analyzed the exome sequence data of 9 tumor-normal pair cases of BPDCN. We utilized SignatureAnalyzer, SigProfiler and a custom microbial analysis pipeline to understand the relevance of endogenous and environmental mutagenic processes.

RESULTS

Our results identified a significant tobacco exposure and aging genetic signature as well as signatures related to nucleotide excision repair deficiency, ultra violet (UV) exposure, and endogenous deamination in BPDCN. We also assessed the samples for microbial infectious disease organisms but did not find a link to a microbial etiology.

CONCLUSION

The identification of a tobacco exposure and aging genetic signature in patients with BPDCN suggests that environmental and endogenous genetic changes may be central to the oncogenesis of BPDCN.

Genome-wide maps of UVA and UVB mutagenesis in yeast reveal distinct causative lesions and mutational strand asymmetries

Ultraviolet (UV) light primarily causes C > T substitutions in lesion-forming dipyrimidine sequences. However, many of the key driver mutations in melanoma do not fit this canonical UV signature, but are instead caused by T > A, T > C, or C > A substitutions. To what extent exposure to the UVB or UVA spectrum of sunlight can induce these noncanonical mutation classes, and the molecular mechanism involved is unclear. Here, we repeatedly exposed wild-type or repair-deficient yeast (Saccharomyces cerevisiae) to UVB or UVA light and characterized the resulting mutations by whole genome sequencing. Our data indicate that UVB induces C > T and T > C substitutions in dipyrimidines, and T > A substitutions that are often associated with thymine-adenine (TA) sequences. All of these mutation classes are induced in nucleotide excision repair-deficient cells and show transcriptional strand asymmetry, suggesting they are caused by helix-distorting UV photoproducts. In contrast, UVA exposure induces orders of magnitude fewer mutations with a distinct mutation spectrum. UVA-induced mutations are elevated in Ogg1-deficient cells, and the resulting spectrum consists almost entirely of C > A/G > T mutations, indicating they are likely derived from oxidative guanine lesions. These mutations show replication asymmetry, with elevated G > T mutations on the leading strand, suggesting there is a strand bias in the removal or bypass of guanine lesions during replication. Finally, we develop a mutation reporter to show that UVA induces a G > T reversion mutation in yeast that mimics the oncogenic NRAS Q61K mutation in melanoma. Taken together, these findings indicate that UVA and UVB exposure can induce many of the noncanonical mutation classes that cause driver mutations in melanoma.

Contributions of replicative and translesion DNA polymerases to mutagenic bypass of canonical and atypical UV photoproducts

UV exposure induces a mutation signature of C > T substitutions at dipyrimidines in skin cancers. We recently identified additional UV-induced AC > TT and A > T substitutions that could respectively cause BRAF V600K and V600E oncogenic mutations. The mutagenic bypass mechanism past these atypical lesions, however, is unknown. Here, we whole genome sequenced UV-irradiated yeast and used reversion reporters to delineate the roles of replicative and translesion DNA polymerases in mutagenic bypass of UV-lesions. Our data indicates that yeast DNA polymerase eta (pol η) has varied impact on UV-induced mutations: protecting against C > T substitutions, promoting T > C and AC > TT substitutions, and not impacting A > T substitutions. Surprisingly, deletion rad30Δ increased novel UV-induced C > A substitutions at CA dinucleotides. In contrast, DNA polymerases zeta (pol ζ) and epsilon (pol ε) participated in AC > TT and A > T mutations. These results uncover lesion-specific accurate and mutagenic bypass of UV lesions, which likely contribute to key driver mutations in melanoma.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Using genomic scars to select immunotherapy beneficiaries in advanced non-small cell lung cancer

In advanced non-small cell lung cancer (NSCLC), response to immunotherapy is difficult to predict from pre-treatment information. Given the toxicity of immunotherapy and its financial burden on the healthcare system, we set out to identify patients for whom treatment is effective. To this end, we used mutational signatures from DNA mutations in pre-treatment tissue. Single base substitutions, doublet base substitutions, indels, and copy number alteration signatures were analysed in [Formula: see text] patients (the discovery set). We found that tobacco smoking signature (SBS4) and thiopurine chemotherapy exposure-associated signature (SBS87) were linked to durable benefit. Combining both signatures in a machine learning model separated patients with a progression-free survival hazard ratio of 0.40[Formula: see text] on the cross-validated discovery set and 0.24[Formula: see text] on an independent external validation set ([Formula: see text]). This paper demonstrates that the fingerprints of mutagenesis, codified through mutational signatures, select advanced NSCLC patients who may benefit from immunotherapy, thus potentially reducing unnecessary patient burden.

Excerno: Using Mutational Signatures in Sequencing Data to Filter False Variants Caused by Clinical Archival

The accurate detection of point mutations from pathology slides using sequencing data is of great importance in cancer genomics and precision oncology. Formalin-fixation paraffin-embedding (FFPE) is a widely used technique to preserve pathology tissues. The FFPE process introduces artificial C > T mutations in next-generation sequencing, so we set out to develop excerno, a method to score and filter such spurious variants. FFPE mutational artifacts follow a mutational signature. By using the FFPE signature and Bayes' formula, we can calculate the probability of a mutation resulting from the FFPE process and use this probability to filter FFPE variants. We implement this method as the excerno R package. We tested excerno by simulating mutations across all 60-baseline mutational signatures from the Catalog of Somatic Mutations in Cancer (COSMIC) and combining them with mutations following the FFPE mutational signature. The sensitivity and specificity of excerno are adversely affected by the cosine similarity between the baseline and FFPE signatures (cosFFPE). Higher percentages of FFPE mutations (pctFFPE) result in increased sensitivity and reduced specificity. The specificity and sensitivity of excerno can be predicted as linear model with an interaction term using cosFFPE and pctFFPE, with an R2of0.84 and 0.79, respectively. Finally, we tested excerno using six RNA sequencing cancer samples and observed concordant trends of specificity and sensitivity with respect to our simulated data. The excerno R package can be used to annotate and filter FFPE-induced mutations in cancer genomics. Our method is adversely affected by cosFFPE and pctFFPE.

Quantification and Mapping of Alkylation in the Human Genome Reveal Single Nucleotide Resolution Precursors of Mutational Signatures

Chemical modifications to DNA bases, including DNA adducts arising from reactions with electrophilic chemicals, are well-known to impact cell growth, miscode during replication, and influence disease etiology. However, knowledge of how genomic sequences and structures influence the accumulation of alkylated DNA bases is not broadly characterized with high resolution, nor have these patterns been linked with overall quantities of modified bases in the genome. For benzo(a) pyrene (BaP), a ubiquitous environmental carcinogen, we developed a single-nucleotide resolution damage sequencing method to map in a human lung cell line the main mutagenic adduct arising from BaP. Furthermore, we combined this analysis with quantitative mass spectrometry to evaluate the dose-response profile of adduct formation. By comparing damage abundance with DNase hypersensitive sites, transcription levels, and other genome annotation data, we found that although overall adduct levels rose with increasing chemical exposure concentration, genomic distribution patterns consistently correlated with chromatin state and transcriptional status. Moreover, due to the single nucleotide resolution characteristics of this DNA damage map, we could determine preferred DNA triad sequence contexts for alkylation accumulation, revealing a characteristic DNA damage signature. This new BaP damage signature had a profile highly similar to mutational signatures identified previously in lung cancer genomes from smokers. Thus, these data provide insight on how genomic features shape the accumulation of alkylation products in the genome and predictive strategies for linking single-nucleotide resolution in vitro damage maps with human cancer mutations.

Arsenic is a potent co-mutagen of ultraviolet light

Environmental co-exposures are widespread and are major contributors to carcinogenic mechanisms. Two well-established environmental agents causing skin cancer are ultraviolet radiation (UVR) and arsenic. Arsenic is a known co-carcinogen that enhances UVR's carcinogenicity. However, the mechanisms of arsenic co-carcinogenesis are not well understood. In this study, we utilized primary human keratinocytes and a hairless mouse model to investigate the carcinogenic and mutagenic properties of co-exposure to arsenic and UVR. In vitro and in vivo exposures revealed that, on its own, arsenic is neither mutagenic nor carcinogenic. However, in combination with UVR, arsenic exposure has a synergistic effect leading to an accelerated mouse skin carcinogenesis as well as to more than 2-fold enrichment of UVR mutational burden. Notably, mutational signature ID13, previously found only in UVR-associated human skin cancers, was observed exclusively in mouse skin tumors and cell lines jointly exposed to arsenic and UVR. This signature was not observed in any model system exposed purely to arsenic or purely to UVR, making ID13 the first co-exposure signature to be reported using controlled experimental conditions. Analysis of existing genomics data from basal cell carcinomas and melanomas revealed that only a subset of human skin cancers harbor ID13 and, consistent with our experimental observations, these cancers exhibited an elevated UVR mutagenesis. Our results provide the first report of a unique mutational signature caused by a co-exposure to two environmental carcinogens and the first comprehensive evidence that arsenic is a potent co-mutagen and co-carcinogen of UVR. Importantly, our findings suggest that a large proportion of human skin cancers are not formed purely due to UVR exposure but rather due to a co-exposure of UVR and other co-mutagens such as arsenic.

The mutational impact of Illudin S on human cells

Illudin S (ILS) is a fungal sesquiterpene secondary metabolite with potent genotoxic and cytotoxic properties. Early genetic studies and more recent genome-wide CRISPR screens showed that Illudin-induced lesions are preferentially repaired by transcription-coupled nucleotide excision repair (TC-NER) with some contribution from post-replication repair pathways. In line with these results, Irofulven, a semi-synthetic ILS analog was recently shown to be particularly effective on cell lines and patient-derived xenografts with impaired NER (e.g. ERCC2/3 mutations), raising hope that ILS-derived molecules may soon enter the clinic. Despite the therapeutic potential of ILS and its analogs, we still lack a global understanding of their mutagenic potential. Here, we characterize the mutational signatures associated with chronic exposure to ILS in human cells. ILS treatment rapidly stalls DNA replication and transcription, leading to the activation of the replication stress response and the accumulation of DNA damage. Novel single and double base substitution signatures as well as a characteristic indel signature indicate that ILS treatment preferentially alkylates purine residues and induces oxidative stress, confirming prior in vitro data. Many mutation contexts exhibit a strong transcriptional strand bias, highlighting the contribution of TC-NER to the repair of ILS lesions. Finally, collateral mutations are also observed in response to ILS, suggesting a contribution of translesion synthesis pathways to ILS tolerance. Accordingly, ILS treatment led to the rapid recruitment of the Y-family DNA polymerase kappa onto chromatin, supporting its preferential use for ILS lesion bypass. Altogether, our work provides the first global assessment of the genomic impact of ILS, demonstrating the contribution of multiple DNA repair pathways to ILS resistance and mutagenicity.

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.

Replicative Instability Drives Cancer Progression

In the past decade, defective DNA repair has been increasingly linked with cancer progression. Human tumors with markers of defective DNA repair and increased replication stress exhibit genomic instability and poor survival rates across tumor types. Seminal studies have demonstrated that genomic instability develops following inactivation of BRCA1, BRCA2, or BRCA-related genes. However, it is recognized that many tumors exhibit genomic instability but lack BRCA inactivation. We sought to identify a pan-cancer mechanism that underpins genomic instability and cancer progression in BRCA-wildtype tumors. Methods: Using multi-omics data from two independent consortia, we analyzed data from dozens of tumor types to identify patient cohorts characterized by poor outcomes, genomic instability, and wildtype BRCA genes. We developed several novel metrics to identify the genetic underpinnings of genomic instability in tumors with wildtype BRCA. Associated clinical data was mined to analyze patient responses to standard of care therapies and potential differences in metastatic dissemination. Results: Systematic analysis of the DNA repair landscape revealed that defective single-strand break repair, translesion synthesis, and non-homologous end-joining effectors drive genomic instability in tumors with wildtype BRCA and BRCA-related genes. Importantly, we find that loss of these effectors promotes replication stress, therapy resistance, and increased primary carcinoma to brain metastasis. Conclusions: Our results have defined a new pan-cancer class of tumors characterized by replicative instability (RIN). RIN is defined by the accumulation of intra-chromosomal, gene-level gain and loss events at replication stress sensitive (RSS) genome sites. We find that RIN accelerates cancer progression by driving copy number alterations and transcriptional program rewiring that promote tumor evolution. Clinically, we find that RIN drives therapy resistance and distant metastases across multiple tumor types.

Establishing Linkages Among DNA Damage, Mutagenesis, and Genetic Diseases

DNA damage by chemicals, radiation, or oxidative stress leads to a mutational spectrum, which is complex because it is determined in part by lesion structure, the DNA sequence context of the lesion, lesion repair kinetics, and the type of cells in which the lesion is replicated. Accumulation of mutations may give rise to genetic diseases such as cancer and therefore understanding the process underlying mutagenesis is of immense importance to preserve human health. Chemical or physical agents that cause cancer often leave their mutational fingerprints, which can be used to back-calculate the molecular events that led to disease. To make a clear link between DNA lesion structure and the mutations a given lesion induces, the field of single-lesion mutagenesis was developed. In the last three decades this area of research has seen much growth in several directions, which we attempt to describe in this Perspective.

Aristolochic acid-associated cancers: a public health risk in need of global action

Aristolochic acids (AAs) are a group of naturally occurring compounds present in many plant species of the Aristolochiaceae family. Exposure to AA is a significant risk factor for severe nephropathy, and urological and hepatobiliary cancers (among others) that are often recurrent and characterized by the prominent mutational fingerprint of AA. However, herbal medicinal products that contain AA continue to be manufactured and marketed worldwide with inadequate regulation, and possible environmental exposure routes receive little attention. As the trade of food and dietary supplements becomes increasingly globalized, we propose that further inaction on curtailing AA exposure will have far-reaching negative effects on the disease trends of AA-associated cancers. Our Review aims to systematically present the historical and current evidence for the mutagenicity and carcinogenicity of AA, and the effect of removing sources of AA exposure on cancer incidence trends. We discuss the persisting challenges of assessing the scale of AA-related carcinogenicity, and the obstacles that must be overcome in curbing AA exposure and preventing associated cancers. Overall, this Review aims to strengthen the case for the implementation of prevention measures against AA's multifaceted, detrimental and potentially fully preventable effects on human cancer development.

Quantification of APOBEC3 Mutation Rates Affecting the VP1 Gene of BK Polyomavirus In Vivo

Mutations in the BK polyomavirus (BKPyV) capsid accumulate in kidney transplant (KTx) recipients with persistent virus replication. They are associated with neutralization escape and appear to arise as a result of cytosine deamination by host cell APOBEC3A/B enzymes. To study the mutagenic processes occurring in patients, we amplified the typing region of the VP1 gene, sequenced the amplicons to a depth of 5000–10,000×, and identified rare mutations, which were fitted to COSMIC mutational signatures. Background mutations were identified in amplicons from plasmids carrying the BKPyV genome and compared to mutations observed in 148 samples from 23 KTx recipients in France and in Vietnam. Three mutational signatures were consistently observed in urine, serum, and kidney biopsy samples, two of which, SBS2 and SBS13, corresponded to APOBEC3A/B activity. In addition, a third signature with no known etiology, SBS89, was detected both in patient samples, and in cells infected in vitro with BKPyV. Quantitatively, APOBEC3A/B mutation rates in urine samples were strongly correlated with urine viral load, and also appeared to vary between individuals. These results confirm that APOBEC3A/B is a major, but not the only, source of BKPyV genome mutations in patients.

Tobacco Use and Response to Immune Checkpoint Inhibitor Therapy in Non-Small Cell Lung Cancer

Tobacco is a known risk factor for lung cancer, and continued tobacco use is associated with poorer outcomes across multiple lung cancer treatment modalities including surgery, chemotherapy and radiation therapy. Less is known about the association of tobacco use and outcomes with immune checkpoint inhibitors (ICIs), which are becoming an important part of the treatment landscape in lung cancer, both in metastatic and curative settings. We reviewed the literature on the association of tobacco and tumor biology as it relates to immunotherapy. We also reviewed the association of tobacco use on outcomes among phase III randomized clinical trials involving ICIs in non-small cell lung cancer (NSCLC). We identified that patients with a smoking history may have a greater benefit with ICI treatment compared to never smokers in both treatment-naïve (HR 0.82, 95% CI 0.69–0.97, vs. HR 1.06, 95% CI 0.81–1.38) and pre-treated (HR 0.79, 95% CI 0.70–0.90 vs. 1.03, 95% CI 0.74–1.43) settings. In trials where smoking status was further defined, ex-smokers appear to demonstrate greater benefit with ICI therapy compared to current smokers (HR 0.78, 95% CI 0.59–1.01 vs. 0.91, 95% CI 0.72–1.14). We conclude by offering our perspective on future directions in this area of research, including implementation of standardized collection and analysis of tobacco use in clinical trials involving ICI therapy in lung cancer and other disease sites, and also evaluating how tobacco may affect toxicities related to ICI therapy. Based on our review, we believe that a patient’s history of tobacco smoking does have a role to play in guiding treatment decision making in patients with lung cancer.

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 and processes in hepatobiliary cancers

The evolutionary history of hepatobiliary cancers is embedded in their genomes. By analysing their catalogue of somatic mutations and the DNA sequence context in which they occur, it is possible to infer the mechanisms underpinning tumorigenesis. These mutational signatures reflect the exogenous and endogenous origins of genetic damage as well as the capacity of hepatobiliary cells to repair and replicate DNA. Genomic analysis of thousands of patients with hepatobiliary cancers has highlighted the diversity of mutagenic processes active in these malignancies, highlighting a prominent source of the inter-cancer-type, inter-patient, intertumour and intratumoural heterogeneity that is observed clinically. However, a substantial proportion of mutational signatures detected in hepatocellular carcinoma and biliary tract cancer remain of unknown cause, emphasizing the important contribution of processes yet to be identified. Exploiting mutational signatures to retrospectively understand hepatobiliary carcinogenesis could advance preventative management of these aggressive tumours as well as potentially predict treatment response and guide the development of therapies targeting tumour evolution.

The Cooked Meat Carcinogen 2-Amino-1-methyl-6-phenylimidazo[4,5-b]pyridine Hair Dosimeter, DNA Adductomics Discovery, and Associations with Prostate Cancer Pathology Biomarkers

Well-done cooked red meat consumption is linked to aggressive prostate cancer (PC) risk. Identifying mutation-inducing DNA adducts in the prostate genome can advance our understanding of chemicals in meat that may contribute to PC. 2-Amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP), a heterocyclic aromatic amine (HAA) formed in cooked meat, is a potential human prostate carcinogen. PhIP was measured in the hair of PC patients undergoing prostatectomy, bladder cancer patients under treatment for cystoprostatectomy, and patients treated for benign prostatic hyperplasia (BPH). PhIP hair levels were above the quantification limit in 123 of 205 subjects. When dichotomizing prostate pathology biomarkers, the geometric mean PhIP hair levels were higher in patients with intermediate and elevated-risk prostate-specific antigen values than lower-risk values <4 ng/mL (p = 0.03). PhIP hair levels were also higher in patients with intermediate and high-risk Gleason scores ≥7 compared to lower-risk Gleason score 6 and BPH patients (p = 0.02). PC patients undergoing prostatectomy had higher PhIP hair levels than cystoprostatectomy or BPH patients (p = 0.02). PhIP-DNA adducts were detected in 9.4% of the patients assayed; however, DNA adducts of other carcinogenic HAAs, and benzo[a]pyrene formed in cooked meat, were not detected. Prostate specimens were also screened for 10 oxidative stress-associated lipid peroxidation (LPO) DNA adducts. Acrolein 1,N²-propano-2'-deoxyguanosine adducts were detected in 54.5% of the patients; other LPO adducts were infrequently detected. Acrolein adducts were not associated with prostate pathology biomarkers, although DNA adductomic profiles differed between PC patients with low and high-grade Gleason scores. Many DNA adducts are of unknown origin; however, dG adducts of formaldehyde and a series of purported 4-hydroxy-2-alkenals were detected at higher abundance in a subset of patients with elevated Gleason scores. The PhIP hair biomarker and DNA adductomics data support the paradigm of well-done cooked meat and oxidative stress in aggressive PC risk.

Discriminating Spontaneous From Cigarette Smoke and THS 2.2 Aerosol Exposure-Related Proliferative Lung Lesions in A/J Mice by Using Gene Expression and Mutation Spectrum Data

Mice, especially A/J mice, have been widely employed to elucidate the underlying mechanisms of lung tumor formation and progression and to derive human-relevant modes of action. Cigarette smoke (CS) exposure induces tumors in the lungs; but, non-exposed A/J mice will also develop lung tumors spontaneously with age, which raises the question of discriminating CS-related lung tumors from spontaneous ones. However, the challenge is that spontaneous tumors are histologically indistinguishable from the tumors occurring in CS-exposed mice. We conducted an 18-month inhalation study in A/J mice to assess the impact of lifetime exposure to Tobacco Heating System (THS) 2.2 aerosol relative to exposure to 3R4F cigarette smoke (CS) on toxicity and carcinogenicity endpoints. To tackle the above challenge, a 13-gene gene signature was developed based on an independent A/J mouse CS exposure study, following by a one-class classifier development based on the current study. Identifying gene signature in one data set and building classifier in another data set addresses the feature/gene selection bias which is a well-known problem in literature. Applied to data from this study, this gene signature classifier distinguished tumors in CS-exposed animals from spontaneous tumors. Lung tumors from THS 2.2 aerosol-exposed mice were significantly different from those of CS-exposed mice but not from spontaneous tumors. The signature was also applied to human lung adenocarcinoma gene expression data (from The Cancer Genome Atlas) and discriminated cancers in never-smokers from those in ever-smokers, suggesting translatability of our signature genes from mice to humans. A possible application of this gene signature is to discriminate lung cancer patients who may benefit from specific treatments (i.e., EGFR tyrosine kinase inhibitors). Mutational spectra from a subset of samples were also utilized for tumor classification, yielding similar results. “Landscaping” the molecular features of A/J mouse lung tumors highlighted, for the first time, a number of events that are also known to play a role in human lung tumorigenesis, such as Lrp1b mutation and Ros1 overexpression. This study shows that omics and computational tools provide useful means of tumor classification where histopathological evaluation alone may be unsatisfactory to distinguish between age- and exposure-related lung tumors.

Signatures of TOP1 transcription-associated mutagenesis in cancer and germline

The mutational landscape is shaped by many processes. Genic regions are vulnerable to mutation but are preferentially protected by transcription-coupled repair1. In microorganisms, transcription has been demonstrated to be mutagenic2,3; however, the impact of transcription-associated mutagenesis remains to be established in higher eukaryotes4. Here we show that ID4-a cancer insertion-deletion (indel) mutation signature of unknown aetiology5 characterized by short (2 to 5 base pair) deletions -is due to a transcription-associated mutagenesis process. We demonstrate that defective ribonucleotide excision repair in mammals is associated with the ID4 signature, with mutations occurring at a TNT sequence motif, implicating topoisomerase 1 (TOP1) activity at sites of genome-embedded ribonucleotides as a mechanistic basis. Such TOP1-mediated deletions occur somatically in cancer, and the ID-TOP1 signature is also found in physiological settings, contributing to genic de novo indel mutations in the germline. Thus, although topoisomerases protect against genome instability by relieving topological stress6, their activity may also be an important source of mutations in the human genome.

Clinical feature-related single-base substitution sequence signatures identified with an unsupervised machine learning approach

BACKGROUND

Mutation processes leave different signatures in genes. For single-base substitutions, previous studies have suggested that mutation signatures are not only reflected in mutation bases but also in neighboring bases. However, because of the lack of a method to identify features of long sequences next to mutation bases, the understanding of how flanking sequences influence mutation signatures is limited.

METHODS

We constructed a long short-term memory-self organizing map (LSTM-SOM) unsupervised neural network. By extracting mutated sequence features via LSTM and clustering similar features with the SOM, single-base substitutions in The Cancer Genome Atlas database were clustered according to both their mutation site and flanking sequences. The relationship between mutation sequence signatures and clinical features was then analyzed. Finally, we clustered patients into different classes according to the composition of the mutation sequence signatures by the K-means method and then studied the differences in clinical features and survival between classes.

RESULTS

Ten classes of mutant sequence signatures (mutation blots, MBs) were obtained from 2,141,527 single-base substitutions via LSTM-SOM machine learning approach. Different features in mutation bases and flanking sequences were revealed among MBs. MBs reflect both the site and pathological features of cancers. MBs were related to clinical features, including age, sex, and cancer stage. The class of an MB in a given gene was associated with survival. Finally, patients were clustered into 7 classes according to the MB composition. Significant differences in survival and clinical features were observed among different patient classes.

CONCLUSIONS

We provided a method for analyzing the characteristics of mutant sequences. Result of this study showed that flanking sequences, together with mutation bases, shape the signatures of SBSs. MBs were shown related to clinical features and survival of cancer patients. Composition of MBs is a feasible predictive factor of clinical prognosis. Further study of the mechanism of MBs related to cancer characteristics is suggested.

The Mutagenic Impact of Environmental Exposures in Human Cells and Cancer: Imprints Through Time

During life, the DNA of our cells is continuously exposed to external damaging processes. Despite the activity of various repair mechanisms, DNA damage eventually results in the accumulation of mutations in the genomes of our cells. Oncogenic mutations are at the root of carcinogenesis, and carcinogenic agents are often highly mutagenic. Over the past decade, whole genome sequencing data of healthy and tumor tissues have revealed how cells in our body gradually accumulate mutations because of exposure to various mutagenic processes. Dissection of mutation profiles based on the type and context specificities of the altered bases has revealed a variety of signatures that reflect past exposure to environmental mutagens, ranging from chemotherapeutic drugs to genotoxic gut bacteria. In this review, we discuss the latest knowledge on somatic mutation accumulation in human cells, and how environmental mutagenic factors further shape the mutation landscapes of tissues. In addition, not all carcinogenic agents induce mutations, which may point to alternative tumor-promoting mechanisms, such as altered clonal selection dynamics. In short, we provide an overview of how environmental factors induce mutations in the DNA of our healthy cells and how this contributes to carcinogenesis. A better understanding of how environmental mutagens shape the genomes of our cells can help to identify potential preventable causes of cancer.

Dissemination of Mycobacterium abscessus via global transmission networks

Mycobacterium abscessus, a multidrug-resistant nontuberculous mycobacterium, has emerged as a major pathogen affecting people with cystic fibrosis (CF). Although originally thought to be acquired independently from the environment, most individuals are infected with one of several dominant circulating clones (DCCs), indicating the presence of global transmission networks of M. abscessus. How and when these clones emerged and spread globally is unclear. Here, we use evolutionary analyses of isolates from individuals both with and without CF to reconstruct the population history, spatiotemporal spread and recent transmission networks of the DCCs. We demonstrate synchronous expansion of six unrelated DCCs in the 1960s, a period associated with major changes in CF care and survival. Each of these clones has spread globally as a result of rare intercontinental transmission events. We show that the DCCs, but not environmentally acquired isolates, exhibit a specific smoking-associated mutational signature and that current transmission networks include individuals both with and without CF. We therefore propose that the DCCs initially emerged in non-CF populations but were then amplified and spread through the CF community. While individuals with CF are probably the most permissive host, non-CF individuals continue to play a key role in transmission networks and may facilitate long-distance transmission.

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.

Co-Exposure to Aristolochic Acids I and II Increases DNA Adduct Formation Responsible for Aristolochic Acid I-Mediated Carcinogenicity in Rats

The plant extract aristolochic acid (AA), containing aristolochic acids I (AAI) and II (AAII) as major components, causes aristolochic acid nephropathy (AAN) and Balkan endemic nephropathy (BEN), unique renal diseases associated with upper urothelial cancer. Recently (Chemical Research in Toxicology 33(11), 2804–2818, 2020), we showed that the in vivo metabolism of AAI and AAII in Wistar rats is influenced by their co-exposure (i.e., AAI/AAII mixture). Using the same rat model, we investigated how exposure to the AAI/AAII mixture can influence AAI and AAII DNA adduct formation (i.e., AA-mediated genotoxicity). Using ³²P-postlabelling, we found that AA-DNA adduct formation was increased in the livers and kidneys of rats treated with AAI/AAII mixture compared to rats treated with AAI or AAII alone. Measuring the activity of enzymes involved in AA metabolism, we showed that enhanced AA-DNA adduct formation might be caused partially by both decreased AAI detoxification as a result of hepatic CYP2C11 inhibition during treatment with AAI/AAII mixture and by hepatic or renal NQO1 induction, the key enzyme predominantly activating AA to DNA adducts. Moreover, our results indicate that AAII might act as an inhibitor of AAI detoxification in vivo. Consequently, higher amounts of AAI might remain in liver and kidney tissues, which can be reductively activated, resulting in enhanced AAI DNA adduct formation. Collectively, these results indicate that AAII present in the plant extract AA enhances the genotoxic properties of AAI (i.e., AAI DNA adduct formation). As patients suffering from AAN and BEN are always exposed to the plant extract (i.e., AAI/AAII mixture), our findings are crucial to better understanding host factors critical for AAN- and BEN-associated urothelial malignancy.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

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.

A comparative analysis of the mutagenicity of platinum-containing chemotherapeutic agents reveals direct and indirect mutagenic mechanisms

Platinum-based drugs are a mainstay of cancer chemotherapy. However, their mutagenic effect can increase tumour heterogeneity, contribute to the evolution of treatment resistance and also induce secondary malignancies. We coupled whole genome sequencing with phenotypic investigations on two cell line models to compare the magnitude and examine the mechanism of mutagenicity of cisplatin, carboplatin and oxaliplatin. Cisplatin induced significantly more base substitution mutations than carboplatin or oxaliplatin when used at equitoxic concentrations on human TK6 or chicken DT40 cells, and also induced the highest number of short insertions and deletions. The analysis of base substitution spectra revealed that all three tested platinum drugs elicit both a direct mutagenic effect at purine dinucleotides, and an indirect effect of accelerating endogenous mutagenic processes, whereas the direct mutagenic effect appeared to correlate with the level of DNA damage caused as assessed through histone H2AX phosphorylation and single-cell agarose gel electrophoresis, the indirect mutagenic effects were equal. The different mutagenicity and DNA-damaging effect of equitoxic platinum drug treatments suggest that DNA damage independent mechanisms significantly contribute to their cytotoxicity. Thus, the comparatively high mutagenicity of cisplatin should be taken into account in the design of chemotherapeutic regimens.

Direct quantification of in vivo mutagenesis and carcinogenesis using duplex sequencing

Error-corrected next-generation sequencing (ecNGS) can be used to rapidly detect and quantify the in vivo mutagenic impact of environmental exposures or endogenous processes in any tissue, from any species, at any genomic location. The greater speed, higher scalability, richer data outputs, and cross-species and cross-locus applicability of ecNGS compared to existing methods make it a powerful new tool for mutational research, regulatory safety testing, and emerging clinical applications.

The ability to accurately measure mutations is critical for basic research and identifying potential drug and chemical carcinogens. Current methods for in vivo quantification of mutagenesis are limited because they rely on transgenic rodent systems that are low-throughput, expensive, prolonged, and do not fully represent other species such as humans. Next-generation sequencing (NGS) is a conceptually attractive alternative for detecting mutations in the DNA of any organism; however, the limit of resolution for standard NGS is poor. Technical error rates (∼1 × 10⁻³) of NGS obscure the true abundance of somatic mutations, which can exist at per-nucleotide frequencies ≤1 × 10⁻⁷. Using duplex sequencing, an extremely accurate error-corrected NGS (ecNGS) technology, we were able to detect mutations induced by three carcinogens in five tissues of two strains of mice within 31 d following exposure. We observed a strong correlation between mutation induction measured by duplex sequencing and the gold-standard transgenic rodent mutation assay. We identified exposure-specific mutation spectra of each compound through trinucleotide patterns of base substitution. We observed variation in mutation susceptibility by genomic region, as well as by DNA strand. We also identified a primordial marker of carcinogenesis in a cancer-predisposed strain of mice, as evidenced by clonal expansions of cells carrying an activated oncogene, less than a month after carcinogen exposure. These findings demonstrate that ecNGS is a powerful method for sensitively detecting and characterizing mutagenesis and the early clonal evolutionary hallmarks of carcinogenesis. Duplex sequencing can be broadly applied to basic mutational research, regulatory safety testing, and emerging clinical applications.

Tethering-facilitated DNA 'opening' and complementary roles of β-hairpin motifs in the Rad4/XPC DNA damage sensor protein

XPC/Rad4 initiates eukaryotic nucleotide excision repair on structurally diverse helix-destabilizing/distorting DNA lesions by selectively 'opening' these sites while rapidly diffusing along undamaged DNA. Previous structural studies showed that Rad4, when tethered to DNA, could also open undamaged DNA, suggesting a 'kinetic gating' mechanism whereby lesion discrimination relied on efficient opening versus diffusion. However, solution studies in support of such a mechanism were lacking and how 'opening' is brought about remained unclear. Here, we present crystal structures and fluorescence-based conformational analyses on tethered complexes, showing that Rad4 can indeed 'open' undamaged DNA in solution and that such 'opening' can largely occur without one or the other of the β-hairpin motifs in the BHD2 or BHD3 domains. Notably, the Rad4-bound 'open' DNA adopts multiple conformations in solution notwithstanding the DNA's original structure or the β-hairpins. Molecular dynamics simulations reveal compensatory roles of the β-hairpins, which may render robustness in dealing with and opening diverse lesions. Our study showcases how fluorescence-based studies can be used to obtain information complementary to ensemble structural studies. The tethering-facilitated DNA 'opening' of undamaged sites and the dynamic nature of 'open' DNA may shed light on how the protein functions within and beyond nucleotide excision repair in cells.

Bedrock radioactivity influences the rate and spectrum of mutation

All organisms on Earth are exposed to low doses of natural radioactivity but some habitats are more radioactive than others. Yet, documenting the influence of natural radioactivity on the evolution of biodiversity is challenging. Here, we addressed whether organisms living in naturally more radioactive habitats accumulate more mutations across generations using 14 species of waterlice living in subterranean habitats with contrasted levels of radioactivity. We found that the mitochondrial and nuclear mutation rates across a waterlouse species’ genome increased on average by 60% and 30%, respectively, when radioactivity increased by a factor of three. We also found a positive correlation between the level of radioactivity and the probability of G to T (and complementary C to A) mutations, a hallmark of oxidative stress. We conclude that even low doses of natural bedrock radioactivity influence the mutation rate possibly through the accumulation of oxidative damage, in particular in the mitochondrial genome.

Atypical UV Photoproducts Induce Non-canonical Mutation Classes Associated with Driver Mutations in Melanoma

Somatic mutations in skin cancers and other ultraviolet (UV)-exposed cells are typified by C>T and CC>TT substitutions at dipyrimidine sequences; however, many oncogenic “driver” mutations in melanoma do not fit this UV signature. Here, we use genome sequencing to characterize mutations in yeast repeatedly irradiated with UV light. Analysis of ~50,000 UV-induced mutations reveals abundant non-canonical mutations, including T>C, T>A, and AC>TT substitutions. These mutations display transcriptional asymmetry that is modulated by nucleotide excision repair (NER), indicating that they are caused by UV photoproducts. Using a sequencing method called UV DNA endonuclease sequencing (UVDE-seq), we confirm the existence of an atypical thymine-adenine photoproduct likely responsible for UV-induced T>A substitutions. Similar non-canonical mutations are present in skin cancers, which also display transcriptional asymmetry and dependence on NER. These include multiple driver mutations, most prominently the recurrent BRAF V600E and V600K substitutions, suggesting that mutations arising from rare, atypical UV photoproducts may play a role in melanomagenesis.

UV mutagenesis has been well studied, but many driver mutations in melanoma do not fit the canonical UV signature. Using whole-genome sequencing, Laughery et al. show that UV induces a broader spectrum of mutations than anticipated. Non-canonical UV mutations are likely caused by atypical photoproducts, which may contribute to melanomagenesis.

The mutational signature profile of known and suspected human carcinogens in mice

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analyzing tumor genomes from mice chronically exposed to 1 of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue-specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-trichloropropane, a haloalkane and pollutant of drinking water, and find these and other signatures in human tumor genomes. We define the cross-species genomic landscape of tumors induced by an important compendium of agents with relevance to human health.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella-another positive single stranded (ss) RNA virus. Each of the rubella virus isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella virus, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2-likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Mutagenicity of acrylamide and glycidamide in human TP53 knock-in (Hupki) mouse embryo fibroblasts

Acrylamide is a suspected human carcinogen formed during high-temperature cooking of starch-rich foods. It is metabolised by cytochrome P450 2E1 to its reactive metabolite glycidamide, which forms pre-mutagenic DNA adducts. Using the human TP53 knock-in (Hupki) mouse embryo fibroblasts (HUFs) immortalisation assay (HIMA), acrylamide- and glycidamide-induced mutagenesis was studied in the tumour suppressor gene TP53. Selected immortalised HUF clones were also subjected to next-generation sequencing to determine mutations across the whole genome. The TP53-mutant frequency after glycidamide exposure (1.1 mM for 24 h, n = 198) was 9% compared with 0% in cultures treated with acrylamide [1.5 (n = 24) or 3 mM (n = 6) for 48 h] and untreated vehicle (water) controls (n = 36). Most glycidamide-induced mutations occurred at adenines with A > T/T > A and A > G/T > C mutations being the most common types. Mutations induced by glycidamide occurred at specific TP53 codons that have also been found to be mutated in human tumours (i.e., breast, ovary, colorectal, and lung) previously associated with acrylamide exposure. The spectrum of TP53 mutations was further reflected by the mutations detected by whole-genome sequencing (WGS) and a distinct WGS mutational signature was found in HUF clones treated with glycidamide that was again characterised by A > G/T > C and A > T/T > A mutations. The WGS mutational signature showed similarities with COSMIC mutational signatures SBS3 and 25 previously found in human tumours (e.g., breast and ovary), while the adenine component was similar to COSMIC SBS4 found mostly in smokers’ lung cancer. In contrast, in acrylamide-treated HUF clones, only culture-related background WGS mutational signatures were observed. In summary, the results of the present study suggest that glycidamide may be involved in the development of breast, ovarian, and lung cancer.

Tumour predisposition and cancer syndromes as models to study gene-environment interactions

Cell division and organismal development are exquisitely orchestrated and regulated processes. The dysregulation of the molecular mechanisms underlying these processes may cause cancer, a consequence of cell-intrinsic and/or cell-extrinsic events. Cellular DNA can be damaged by spontaneous hydrolysis, reactive oxygen species, aberrant cellular metabolism or other perturbations that cause DNA damage. Moreover, several environmental factors may damage the DNA, alter cellular metabolism or affect the ability of cells to interact with their microenvironment. While some environmental factors are well established as carcinogens, there remains a large knowledge gap of others owing to the difficulty in identifying them because of the typically long interval between carcinogen exposure and cancer diagnosis. DNA damage increases in cells harbouring mutations that impair their ability to correctly repair the DNA. Tumour predisposition syndromes in which cancers arise at an accelerated rate and in different organs - the equivalent of a sensitized background - provide a unique opportunity to examine how gene-environment interactions influence cancer risk when the initiating genetic defect responsible for malignancy is known. Understanding the molecular processes that are altered by specific germline mutations, environmental exposures and related mechanisms that promote cancer will allow the design of novel and effective preventive and therapeutic strategies.

Next-Generation Sequencing Analysis Identified Genomic Alterations in Pathological Morphologies of 3 Cases of Pulmonary Carcinosarcoma

Background

Pulmonary carcinosarcomas (PCSs) are a heterogeneous group of non-small-cell lung carcinomas (NSCLCs) with aggressiveness and a poor prognosis. Although genetic mutations of some common lung cancer subtypes have been extensively studied, the molecular characteristics of PCSs and the existence of abnormal target genes are unknown.

Methods

In this study, the clinical and molecular characterization in 3 pulmonary sarcomatoid carcinomas (PSCs) were presented using microscope analysis and next-generation sequencing (NGS) analysis.

Results

The results revealed a carcinosarcomas subtype presenting squamous cell carcinoma and sarcoma components in all 3 cases. NGS analysis showed that 182, 316 and 230 shared mutations were detected between sarcoma and lung carcinoma from 3 patients. Two identical alterations in two genes (CSMD3 and RYR3) that were all shared by the two components in 3 patients. Tumor suppressor gene TP53 (5/6, 83%) showed the highest mutation frequency for driver genes here. Additionally, we focused on an LYST mutation which was mainly present in the sarcoma components. Moreover, the clonal evolution and signature analysis confirm that lung squamous cell carcinoma and sarcoma in each PCS patient may have come from a common ancestor, and mutagenesis was possibly related to indirect effects of tobacco, age or other unknown factors.

Conclusion

Our results indicate that genetic analysis and molecular targeted therapy are necessary for the identification and treatment of these rare lung tumors. CSMD3 and LYST, as common mutation genes, may be a potential therapeutic target in PCS.

Mutational Landscape of Esophageal Squamous Cell Carcinoma in an Indian Cohort

Esophageal squamous cell carcinoma (ESCC) is the most common histological subtype of esophageal cancer in India. Cigarette smoking and chewing tobacco are known risk factors associated with ESCC. However, genomic alterations associated with ESCC in India are not well-characterized. In this study, we carried out exome sequencing to characterize the mutational landscape of ESCC tumors from subjects with a varied history of tobacco usage. Whole exome sequence analysis of ESCC from an Indian cohort revealed several genes that were mutated or had copy number changes. ESCC from tobacco chewers had a higher frequency of C:G > A:T transversions and 2-fold enrichment for mutation signature 4 compared to smokers and non-users of tobacco. Genes, such as TP53, CSMD3, SYNE1, PIK3CA, and NOTCH1 were found to be frequently mutated in Indian cohort. Mutually exclusive mutation patterns were observed in PIK3CA-NOTCH1, DNAH5-ZFHX4, MUC16-FAT1, and ZFHX4-NOTCH1 gene pairs. Recurrent amplifications were observed in 3q22-3q29, 11q13.3-q13.4, 7q22.1-q31.1, and 8q24 regions. Approximately 53% of tumors had genomic alterations in PIK3CA making this pathway a promising candidate for targeted therapy. In conclusion, we observe enrichment of mutation signature 4 in ESCC tumors from patients with a history of tobacco chewing. This is likely due to direct exposure of esophagus to tobacco carcinogens when it is chewed and swallowed. Genomic alterations were frequently observed in PIK3CA-AKT pathway members independent of the history of tobacco usage. PIK3CA pathway can be potentially targeted in ESCC which currently has no effective targeted therapeutic options.

Chemically induced mutations in a MutaMouse reporter gene inform mechanisms underlying human cancer mutational signatures

Transgenic rodent (TGR) models use bacterial reporter genes to quantify in vivo mutagenesis. Pairing TGR assays with next-generation sequencing (NGS) enables comprehensive mutation pattern analysis to inform mutational mechanisms. We used this approach to identify 2751 independent lacZ mutations in the bone marrow of MutaMouse animals exposed to four chemical mutagens: benzo[a]pyrene, N-ethyl-N-nitrosourea, procarbazine, and triethylenemelamine. We also collected published data for 706 lacZ mutations from eight additional environmental mutagens. We report that lacZ gene sequencing generates chemical-specific mutation signatures observed in human cancers with established environmental causes. For example, the mutation signature of benzo[a]pyrene, a carcinogen present in tobacco smoke, matched the signature associated with tobacco-induced lung cancers. Our results suggest that the analysis of chemically induced mutations in the lacZ gene shortly after exposure provides an effective approach to characterize human-relevant mechanisms of carcinogenesis and propose novel environmental causes of mutation signatures observed in human cancers.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella - another positive single stranded (ss) RNA virus. Each of the rubella isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2 - likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.