MutLα suppresses error-prone DNA mismatch repair and preferentially protects noncoding DNA from mutations.

The DNA mismatch repair (MMR) system promotes genome stability and protects humans from certain types of cancer. Its primary function is the correction of DNA polymerase errors. MutLα is an important eukaryotic MMR factor. We have examined the contributions of MutLα to maintaining genome stability. We show here that loss of MutLα in yeast increases the genome-wide mutation rate by ∼130-fold and generates a genome-wide mutation spectrum that consists of small indels and base substitutions. We also show that loss of yeast MutLα leads to error-prone MMR that produces T>C base substitutions in 5'-ATA-3' sequences. In agreement with this finding, our examination of human whole genome DNA sequencing data has revealed that loss of MutLα in induced pluripotent stem cells triggers error-prone MMR that leads to the formation of T>C mutations in 5'-NTN-3' sequences. Our further analysis has shown that MutLα-independent MMR plays a role in suppressing base substitutions in N3 homopolymeric runs. In addition, we describe that MutLα preferentially protects noncoding DNA from mutations. Our study defines the contributions of MutLα-dependent and independent mechanisms to genome-wide MMR.

MutLα suppresses error-prone DNA mismatch repair and preferentially protects noncoding DNA from mutations.

The DNA mismatch repair (MMR) system promotes genome stability and protects humans from certain types of cancer. Its primary function is the correction of DNA polymerase errors. MutLα is an important eukaryotic MMR factor. We have examined the contributions of MutLα to maintaining genome stability. We show here that loss of MutLα in yeast increases the genome-wide mutation rate by ~130-fold and generates a genome-wide mutation spectrum that consists of small indels and base substitutions. We also show that loss of yeast MutLα leads to error-prone MMR that produces T>C base substitutions in 5'-ATA-3' sequences. In agreement with this finding, our examination of human whole genome DNA sequencing data has revealed that loss of MutLα in induced pluripotent stem cells triggers error-prone MMR that leads to the formation of T>C mutations in 5'-NTN-3' sequences. Our further analysis has shown that MutLα-independent MMR plays a role in suppressing base substitutions in N3 homopolymeric runs. In addition, we describe that MutLα preferentially defends noncoding DNA from mutations. Our study defines the contributions of MutLα-dependent and independent mechanisms to genome-wide MMR.

Multiple DNA repair pathways prevent acetaldehyde-induced mutagenesis in yeast.

Acetaldehyde is the primary metabolite of alcohol and is present in many environmental sources including tobacco smoke. Acetaldehyde is genotoxic, whereby it can form DNA adducts and lead to mutagenesis. Individuals with defects in acetaldehyde clearance pathways have increased susceptibility to alcohol-associated cancers. Moreover, a mutation signature specific to acetaldehyde exposure is widespread in alcohol and smoking-associated cancers. However, the pathways that repair acetaldehyde-induced DNA damage and thus prevent mutagenesis are vaguely understood. Here, we used Saccharomyces cerevisiae to systematically delete genes in each of the major DNA repair pathways to identify those that alter acetaldehyde-induced mutagenesis. We found that deletion of the nucleotide excision repair (NER) genes, RAD1 or RAD14, led to an increase in mutagenesis upon acetaldehyde exposure. Acetaldehyde-induced mutations were dependent on translesion synthesis as well as DNA inter-strand crosslink (ICL) repair in Δrad1 strains. Moreover, whole genome sequencing of the mutated isolates demonstrated an increase in C→A changes coupled with an enrichment of gCn→A changes in the acetaldehyde-treated Δrad1 isolates. The gCn→A mutation signature has been shown to be diagnostic of acetaldehyde exposure in yeast and in human cancers. We also demonstrated that the deletion of the two DNA-protein crosslink (DPC) repair proteases, WSS1 and DDI1, also led to increased acetaldehyde-induced mutagenesis. Defects in base excision repair (BER) led to a mild increase in mutagenesis, while defects in mismatch repair (MMR), homologous recombination repair (HR) and post replicative repair pathways did not impact mutagenesis upon acetaldehyde exposure. Our results in yeast were further corroborated upon analysis of whole exome sequenced liver cancers, wherein, tumors with defects in ERCC1 and ERCC4 (NER), FANCD2 (ICL repair) or SPRTN (DPC repair) carried a higher gCn→A mutation load than tumors with no deleterious mutations in these genes. Our findings demonstrate that multiple DNA repair pathways protect against acetaldehyde-induced mutagenesis.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans.

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

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.

Genome-wide contributions of the MutSα- and MutSß-dependent DNA mismatch repair pathways to the maintenance of genetic stability in Saccharomyces cerevisiae.

The DNA mismatch repair (MMR) system is a major DNA repair system that suppresses both inherited and sporadic cancers in humans. In eukaryotes, the MutSα-dependent and MutSß-dependent MMR pathways correct DNA polymerase errors. Here, we investigated these two pathways on a whole genome level in Saccharomyces cerevisiae. We found that inactivation of MutSα-dependent MMR increases the genome-wide mutation rate by ∼17-fold and loss of MutSß-dependent MMR elevates the genome-wide mutation rate by ∼4-fold. We also found that MutSα-dependent MMR does not show a preference for protecting coding or noncoding DNA from mutations, whereas MutSß-dependent MMR preferentially protects noncoding DNA from mutations. The most frequent mutations in the msh6Δ strain are C>T transitions, whereas 1- to 6-bp deletions are the most common genetic alterations in the msh3Δ strain. Strikingly, MutSα-dependent MMR is more important than MutSß-dependent MMR for protection from 1-bp insertions, while MutSß-dependent MMR has a more critical role in the defense against 1-bp deletions and 2- to 6-bp indels. We also determined that a mutational signature of yeast MSH6 loss is similar to mutational signatures of human MMR deficiency. Furthermore, our analysis showed that compared to other 5'-NCN-3' trinucleotides, 5'-GCA-3' trinucleotides are at the highest risk of accumulating C>T transitions at the central position in the msh6Δ cells and that the presence of a G/A base at the -1 position is important for the efficient MutSα-dependent suppression of C>T transitions. Our results highlight key differences between the roles of the MutSα-dependent and MutSß-dependent MMR pathways.

Acetaldehyde makes a distinct mutation signature in single-stranded DNA.

Acetaldehyde (AA), a by-product of ethanol metabolism, is acutely toxic due to its ability to react with various biological molecules including DNA and proteins, which can greatly impede key processes such as replication and transcription and lead to DNA damage. As such AA is classified as a group 1 carcinogen by the International Agency for Research on Cancer (IARC). Previous in vitro studies have shown that AA generates bulky adducts on DNA, with signature guanine-centered (GG→TT) mutations. However, due to its weak mutagenicity, short chemical half-life, and the absence of powerful genetic assays, there is considerable variability in reporting the mutagenic effects of AA in vivo. Here, we used an established yeast genetic reporter system and demonstrate that AA treatment is highly mutagenic to cells and leads to strand-biased mutations on guanines (G→T) at a high frequency on single stranded DNA (ssDNA). We further demonstrate that AA-derived mutations occur through lesion bypass on ssDNA by the translesion polymerase Polζ. Finally, we describe a unique mutation signature for AA, which we then identify in several whole-genome and -exome sequenced cancers, particularly those associated with alcohol consumption. Our study proposes a key mechanism underlying carcinogenesis by acetaldehyde-mutagenesis of single-stranded DNA.

Internode elongation in energy cane shows remarkable clues on lignocellulosic biomass biosynthesis in Saccharum hybrids.

Energy cane is a dedicated crop to high biomass production and selected during Saccharum breeding programs to fit specific industrial needs for 2G bioethanol production. Internode elongation is one of the most important characteristics in Saccharum hybrids due to its relationship with crop yield. In this study, we selected the third internode elongation of the energy cane. To characterize this process, we divided the internode into five sections and performed a detailed transcriptome analysis (RNA-Seq) and cell wall characterization. The histological analyses revealed a remarkable gradient that spans from cell division and protoxylem lignification to the internode maturation and complete vascular bundle lignification. RNA-Seq analysis revealed more than 11,000 differentially expressed genes between the sections internal. Gene ontology analyzes showed enriched categories in each section, as well as the most expressed genes in each section, presented different biological processes. We found that the internode elongation and division zones have a large number of unique genes. Evaluated the specific profile of genes related to primary and secondary cell wall formation, cellulose synthesis, hemicellulose, lignin, and growth-related genes. For each section these genes presented different profiles along the internode in elongation in energy cane. The results of this study provide an overview of the regulation of gene expression of an internode elongation in energy cane. Gene expression analysis revealed promising candidates for transcriptional regulation of energy cane lignification and evidence key genes for the regulation of internode development, which can serve as a basis for understanding the molecular regulatory mechanisms that support the growth and development of plants in the Saccahrum complex.

Noninvasive prenatal exome sequencing diagnostic utility limited by sequencing depth and fetal fraction.

OBJECTIVE: Sequencing cell-free DNA now allows detection of large chromosomal abnormalities and dominant Mendelian disorders in the prenatal period. Improving upon these methods would allow newborn screening programs to begin with prenatal genetics, ultimately improving the management of rare genetic disorders. METHODS: As a pilot study, we performed exome sequencing on the cell-free DNA from three mothers with singleton pregnancies to assess the viability of broad sequencing modalities in a noninvasive prenatal setting. RESULTS: We found poor resolution of maternal and fetal genotypes due to both sampling and technical issues. CONCLUSION: We find broad sequencing modalities inefficient for noninvasive prenatal applications. Alternatively, we suggest a more targeted path forward for noninvasive prenatal genotyping.

The Shu complex prevents mutagenesis and cytotoxicity of single-strand specific alkylation lesions.

Three-methyl cytosine (3meC) are toxic DNA lesions, blocking base pairing. Bacteria and humans express members of the AlkB enzymes family, which directly remove 3meC. However, other organisms, including budding yeast, lack this class of enzymes. It remains an unanswered evolutionary question as to how yeast repairs 3meC, particularly in single-stranded DNA. The yeast Shu complex, a conserved homologous recombination factor, aids in preventing replication-associated mutagenesis from DNA base damaging agents such as methyl methanesulfonate (MMS). We found that MMS-treated Shu complex-deficient cells exhibit a genome-wide increase in A:T and G:C substitutions mutations. The G:C substitutions displayed transcriptional and replicational asymmetries consistent with mutations resulting from 3meC. Ectopic expression of a human AlkB homolog in Shu-deficient yeast rescues MMS-induced growth defects and increased mutagenesis. Thus, our work identifies a novel homologous recombination-based mechanism mediated by the Shu complex for coping with alkylation adducts.

Pre-capture multiplexing provides additional power to detect copy number variation in exome sequencing.

BACKGROUND: As exome sequencing (ES) integrates into clinical practice, we should make every effort to utilize all information generated. Copy-number variation can lead to Mendelian disorders, but small copy-number variants (CNVs) often get overlooked or obscured by under-powered data collection. Many groups have developed methodology for detecting CNVs from ES, but existing methods often perform poorly for small CNVs and rely on large numbers of samples not always available to clinical laboratories. Furthermore, methods often rely on Bayesian approaches requiring user-defined priors in the setting of insufficient prior knowledge. This report first demonstrates the benefit of multiplexed exome capture (pooling samples prior to capture), then presents a novel detection algorithm, mcCNV ("multiplexed capture CNV"), built around multiplexed capture. RESULTS: We demonstrate: (1) multiplexed capture reduces inter-sample variance; (2) our mcCNV method, a novel depth-based algorithm for detecting CNVs from multiplexed capture ES data, improves the detection of small CNVs. We contrast our novel approach, agnostic to prior information, with the the commonly-used ExomeDepth. In a simulation study mcCNV demonstrated a favorable false discovery rate (FDR). When compared to calls made from matched genome sequencing, we find the mcCNV algorithm performs comparably to ExomeDepth. CONCLUSION: Implementing multiplexed capture increases power to detect single-exon CNVs. The novel mcCNV algorithm may provide a more favorable FDR than ExomeDepth. The greatest benefits of our approach derive from (1) not requiring a database of reference samples and (2) not requiring prior information about the prevalance or size of variants.

Evaluating the analytical validity of circulating tumor DNA sequencing assays for precision oncology.

Circulating tumor DNA (ctDNA) sequencing is being rapidly adopted in precision oncology, but the accuracy, sensitivity and reproducibility of ctDNA assays is poorly understood. Here we report the findings of a multi-site, cross-platform evaluation of the analytical performance of five industry-leading ctDNA assays. We evaluated each stage of the ctDNA sequencing workflow with simulations, synthetic DNA spike-in experiments and proficiency testing on standardized, cell-line-derived reference samples. Above 0.5% variant allele frequency, ctDNA mutations were detected with high sensitivity, precision and reproducibility by all five assays, whereas, below this limit, detection became unreliable and varied widely between assays, especially when input material was limited. Missed mutations (false negatives) were more common than erroneous candidates (false positives), indicating that the reliable sampling of rare ctDNA fragments is the key challenge for ctDNA assays. This comprehensive evaluation of the analytical performance of ctDNA assays serves to inform best practice guidelines and provides a resource for precision oncology.

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.

Characterization of systemic genomic instability in budding yeast.

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Adaptation and selection shape clonal evolution of tumors during residual disease and recurrence.

The survival and recurrence of residual tumor cells following therapy constitutes one of the biggest obstacles to obtaining cures in breast cancer, but it remains unclear how the clonal composition of tumors changes during relapse. We use cellular barcoding to monitor clonal dynamics during tumor recurrence in vivo. We find that clonal diversity decreases during tumor regression, residual disease, and recurrence. The recurrence of dormant residual cells follows several distinct routes. Approximately half of the recurrent tumors exhibit clonal dominance with a small number of subclones comprising the vast majority of the tumor; these clonal recurrences are frequently dependent upon Met gene amplification. A second group of recurrent tumors comprises thousands of subclones, has a clonal architecture similar to primary tumors, and is dependent upon the Jak/Stat pathway. Thus the regrowth of dormant tumors proceeds via multiple routes, producing recurrent tumors with distinct clonal composition, genetic alterations, and drug sensitivities.

Tn-Seq Analysis Identifies Genes Important for Yersinia pestis Adherence during Primary Pneumonic Plague.

Following inhalation, Yersinia pestis rapidly colonizes the lung to establish infection during primary pneumonic plague. Although several adhesins have been identified in Yersinia spp., the factors mediating early Y. pestis adherence in the lung remain unknown. To identify genes important for Y. pestis adherence during primary pneumonic plague, we used transposon insertion sequencing (Tn-seq). Wild-type and capsule mutant (Δcaf1) Y. pestis transposon mutant libraries were serially passaged in vivo to enrich for nonadherent mutants in the lung using a mouse model of primary pneumonic plague. Sequencing of the passaged libraries revealed six mutants that were significantly enriched in both the wild-type and Δcaf1Y. pestis backgrounds. The enriched mutants had insertions in genes that encode transcriptional regulators, chaperones, an endoribonuclease, and YPO3903, a hypothetical protein. Using single-strain infections and a transcriptional analysis, we identified a significant role for YPO3903 in Y. pestis adherence in the lung and showed that YPO3903 regulated transcript levels of psaA, which encodes a fimbria previously implicated in Y. pestis adherence in vitro Deletion of psaA had a minor effect on Y. pestis adherence in the lung, suggesting that YPO3903 regulates other adhesins in addition to psaA By enriching for mutations in genes that regulate the expression or assembly of multiple genes or proteins, we obtained screen results indicating that there may be not just one dominant adhesin but rather several factors that contribute to early Y. pestis adherence during primary pneumonic plague.IMPORTANCE Colonization of the lung by Yersinia pestis is a critical first step in establishing infection during primary pneumonic plague, a disease characterized by high lethality. However, the mechanisms by which Y. pestis adheres in the lung after inhalation remain elusive. Here, we used Tn-seq to identify Y. pestis genes important for adherence early during primary pneumonic plague. Our mutant enrichment strategy resulted in the identification of genes important for regulation and assembly of genes and proteins rather than adhesin genes themselves. These results reveal that there may be multiple Y. pestis adhesins or redundancy among adhesins. Identifying the adhesins regulated by the genes identified in our enrichment screen may reveal novel therapeutic targets for preventing Y. pestis adherence and the subsequent development of pneumonic plague.

Long transposon-rich centromeres in an oomycete reveal divergence of centromere features in Stramenopila-Alveolata-Rhizaria lineages.

Centromeres are chromosomal regions that serve as platforms for kinetochore assembly and spindle attachments, ensuring accurate chromosome segregation during cell division. Despite functional conservation, centromere DNA sequences are diverse and often repetitive, making them challenging to assemble and identify. Here, we describe centromeres in an oomycete Phytophthora sojae by combining long-read sequencing-based genome assembly and chromatin immunoprecipitation for the centromeric histone CENP-A followed by high-throughput sequencing (ChIP-seq). P. sojae centromeres cluster at a single focus at different life stages and during nuclear division. We report an improved genome assembly of the P. sojae reference strain, which enabled identification of 15 enriched CENP-A binding regions as putative centromeres. By focusing on a subset of these regions, we demonstrate that centromeres in P. sojae are regional, spanning 211 to 356 kb. Most of these regions are transposon-rich, poorly transcribed, and lack the histone modification H3K4me2 but are embedded within regions with the heterochromatin marks H3K9me3 and H3K27me3. Strikingly, we discovered a Copia-like transposon (CoLT) that is highly enriched in the CENP-A chromatin. Similar clustered elements are also found in oomycete relatives of P. sojae, and may be applied as a criterion for prediction of oomycete centromeres. This work reveals a divergence of centromere features in oomycetes as compared to other organisms in the Stramenopila-Alveolata-Rhizaria (SAR) supergroup including diatoms and Plasmodium falciparum that have relatively short and simple regional centromeres. Identification of P. sojae centromeres in turn also advances the genome assembly.

Correction: Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

[This corrects the article DOI: 10.1371/journal.pbio.3000464.].

Mutation signatures specific to DNA alkylating agents in yeast and cancers.

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.