DNA polymerase zeta generates clustered mutations during bypass of endogenous DNA lesions in Saccharomyces cerevisiae

Interplay of replication timing, DNA repair, and translesion synthesis in UV mutagenesis in yeast

Replication timing during S-phase impacts mutation rates in yeast and human cancers; however, the exact mechanism involved remains unclear. Here, we analyze the impact of replication timing on UV mutagenesis in Saccharomyces cerevisiae. Our analysis indicates that UV mutations are enriched in early-replicating regions of the genome in wild-type cells, but in cells deficient in global genomic-nucleotide excision repair (GG-NER), mutations are enriched in late-replicating regions. Analysis of UV damage maps revealed that cyclobutane pyrimidine dimers are enriched in late-replicating regions, but this enrichment is almost entirely due to repetitive ribosomal DNA. Complex mutations typically associated with TLS activity are also elevated in late-replicating regions in GG-NER deficient cells. We propose that UV mutagenesis is higher in early-replicating regions in repair-competent cells because there is less time to repair the lesion prior to undergoing replication. However, in the absence of GG-NER, increased TLS activity promotes UV mutagenesis in late-replicating regions.

Perspectives and opportunities in forensic human, animal, and plant integrative genomics in the Pangenome era

The Human Pangenome Reference Consortium, the Chinese Pangenome Consortium, and other plant and animal pangenome projects have announced the completion of pilot work aimed at constructing high-quality, haplotype-resolved reference graph genomes representative of global ethno-linguistically different populations or different plant and animal species. These graph-based, gapless pangenome references, which are enriched in terms of genomic diversity, completeness, and contiguity, have the potential for enhancing long-read sequencing (LRS)-based genomic research, as well as improving mappability and variant genotyping on traditional short-read sequencing platforms. We comprehensively discuss the advancements in pangenome-based genomic integrative genomic discoveries across forensic-related species (humans, animals, and plants) and summarize their applications in variant identification and forensic genomics, epigenetics, transcriptomics, and microbiome research. Recent developments in multiplexed array sequencing have introduced a highly efficient and programmable technique to overcome the limitations of short forensic marker lengths in LRS platforms. This technique enables the concatenation of short RNA transcripts and DNA fragments into LRS-optimal molecules for sequencing, assembly, and genotyping. The integration of new pangenome reference coordinates and corresponding computational algorithms will benefit forensic integrative genomics by facilitating new marker identification, accurate genotyping, high-resolution panel development, and the updating of statistical algorithms. This review highlights the necessity of integrating LRS-based platforms, pangenome-based study designs, and graph-based pangenome references in short-read mapping and LRS-based innovations to achieve precision forensic science.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

Revisiting mutagenesis at non-B DNA motifs in the human genome

Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Characterisation of the spectrum and genetic dependence of collateral mutations induced by translesion DNA synthesis

Translesion DNA synthesis (TLS) is a fundamental damage bypass pathway that utilises specialised polymerases with relaxed template specificity to achieve replication through damaged DNA. Misinsertions by low fidelity TLS polymerases may introduce additional mutations on undamaged DNA near the original lesion site, which we termed collateral mutations. In this study, we used whole genome sequencing datasets of chicken DT40 and several human cell lines to obtain evidence for collateral mutagenesis in higher eukaryotes. We found that cisplatin and UVC radiation frequently induce close mutation pairs within 25 base pairs that consist of an adduct-associated primary and a downstream collateral mutation, and genetically linked their formation to TLS activity involving PCNA ubiquitylation and polymerase κ. PCNA ubiquitylation was also indispensable for close mutation pairs observed amongst spontaneously arising base substitutions in cell lines with disrupted homologous recombination. Collateral mutation pairs were also found in melanoma genomes with evidence of UV exposure. We showed that collateral mutations frequently copy the upstream base, and extracted a base substitution signature that describes collateral mutagenesis in the presented dataset regardless of the primary mutagenic process. Using this mutation signature, we showed that collateral mutagenesis creates approximately 10–20% of non-paired substitutions as well, underscoring the importance of the process.

DNA base substitutions are the most common form of genomic mutations, formed both spontaneously and in response to environmental mutagens. One of the main mechanisms of base substitution mutagenesis is translesion synthesis, a process that relies on specialised DNA polymerases to replicate damaged DNA templates. In addition to incorrect base insertions at the site of lesions in the template, translesion polymerases may also generate ‘collateral’ mutations away from the lesion due to their lower accuracy in selecting the correct incoming nucleotide. In this study, we surveyed the whole genome sequence of experimental cell clones to examine the extent and genetic dependence of collateral mutagenesis in higher eukaryotes. Looking for close mutation pairs, we found that collateral mutations frequently occur near primary lesions generated by cisplatin or ultraviolet radiation in chicken and human cells, but are restricted to a short distance of approximately 25 base pairs. By analysing their sequence context, we showed that collateral mutations can also occur near correctly bypassed primary lesions and may be responsible for a considerable proportion of all base substitution mutations.

Tandem Substitutions in Somatic Hypermutation

Upon antigen recognition, activation-induced cytosine deaminase initiates affinity maturation of the B-cell receptor by somatic hypermutation (SHM) through error-prone DNA repair pathways. SHM typically creates single nucleotide substitutions, but tandem substitutions may also occur. We investigated incidence and sequence context of tandem substitutions by massive parallel sequencing of V(D)J repertoires in healthy human donors. Mutation patterns were congruent with SHM-derived single nucleotide mutations, delineating initiation of the tandem substitution by AID. Tandem substitutions comprised 5,7% of AID-induced mutations. The majority of tandem substitutions represents single nucleotide juxtalocations of directly adjacent sequences. These observations were confirmed in an independent cohort of healthy donors. We propose a model where tandem substitutions are predominantly generated by translesion synthesis across an apyramidinic site that is typically created by UNG. During replication, apyrimidinic sites transiently adapt an extruded configuration, causing skipping of the extruded base. Consequent strand decontraction leads to the juxtalocation, after which exonucleases repair the apyramidinic site and any directly adjacent mismatched base pairs. The mismatch repair pathway appears to account for the remainder of tandem substitutions. Tandem substitutions may enhance affinity maturation and expedite the adaptive immune response by overcoming amino acid codon degeneracies or mutating two adjacent amino acid residues simultaneously.

A modified fluctuation assay reveals a natural mutator phenotype that drives mutation spectrum variation within Saccharomyces cerevisiae

Although studies of Saccharomyces cerevisiae have provided many insights into mutagenesis and DNA repair, most of this work has focused on a few laboratory strains. Much less is known about the phenotypic effects of natural variation within S. cerevisiae’s DNA repair pathways. Here, we use natural polymorphisms to detect historical mutation spectrum differences among several wild and domesticated S. cerevisiae strains. To determine whether these differences are likely caused by genetic mutation rate modifiers, we use a modified fluctuation assay with a CAN1 reporter to measure de novo mutation rates and spectra in 16 of the analyzed strains. We measure a 10-fold range of mutation rates and identify two strains with distinctive mutation spectra. These strains, known as AEQ and AAR, come from the panel’s ‘Mosaic beer’ clade and share an enrichment for C > A mutations that is also observed in rare variation segregating throughout the genomes of several Mosaic beer and Mixed origin strains. Both AEQ and AAR are haploid derivatives of the diploid natural isolate CBS 1782, whose rare polymorphisms are enriched for C > A as well, suggesting that the underlying mutator allele is likely active in nature. We use a plasmid complementation test to show that AAR and AEQ share a mutator allele in the DNA repair gene OGG1, which excises 8-oxoguanine lesions that can cause C > A mutations if left unrepaired.

Protection of the C. elegans germ cell genome depends on diverse DNA repair pathways during normal proliferation

Maintaining genome integrity is particularly important in germ cells to ensure faithful transmission of genetic information across generations. Here we systematically describe germ cell mutagenesis in wild-type and 61 DNA repair mutants cultivated over multiple generations. ~44% of the DNA repair mutants analysed showed a >2-fold increased mutagenesis with a broad spectrum of mutational outcomes. Nucleotide excision repair deficiency led to higher base substitution rates, whereas polh-1(Polη) and rev-3(Polζ) translesion synthesis polymerase mutants resulted in 50-400 bp deletions. Signatures associated with defective homologous recombination fall into two classes: 1) brc-1/BRCA1 and rad-51/RAD51 paralog mutants showed increased mutations across all mutation classes, 2) mus-81/MUS81 and slx-1/SLX1 nuclease, and him-6/BLM, helq-1/HELQ or rtel-1/RTEL1 helicase mutants primarily accumulated structural variants. Repetitive and G-quadruplex sequence-containing loci were more frequently mutated in specific DNA repair backgrounds. Tandem duplications embedded in inverted repeats were observed in helq-1 helicase mutants, and a unique pattern of 'translocations' involving homeologous sequences occurred in rip-1 recombination mutants. atm-1/ATM checkpoint mutants harboured structural variants specifically enriched in subtelomeric regions. Interestingly, locally clustered mutagenesis was only observed for combined brc-1 and cep-1/p53 deficiency. Our study provides a global view of how different DNA repair pathways contribute to prevent germ cell mutagenesis.

XPC deficiency increases risk of hematologic malignancies through mutator phenotype and characteristic mutational signature

Recent studies demonstrated a dramatically increased risk of leukemia in patients with a rare genetic disorder, Xeroderma Pigmentosum group C (XP-C), characterized by constitutive deficiency of global genome nucleotide excision repair (GG-NER). The genetic mechanisms of non-skin cancers in XP-C patients remain unexplored. In this study, we analyze a unique collection of internal XP-C tumor genomes including 6 leukemias and 2 sarcomas. We observe a specific mutational pattern and an average of 25-fold increase of mutation rates in XP-C versus sporadic leukemia which we presume leads to its elevated incidence and early appearance. We describe a strong mutational asymmetry with respect to transcription and the direction of replication in XP-C tumors suggesting association of mutagenesis with bulky purine DNA lesions of probably endogenous origin. These findings suggest existence of a balance between formation and repair of bulky DNA lesions by GG-NER in human body cells which is disrupted in XP-C patients.

Xeroderma Pigmentosum group C (XP-C) is a rare genetic disorder characterised by deficient DNA repair leading to skin and internal cancer, but the latter is not well understood molecularly. Here the authors sequence genomes of non-skin cancers from XP-C patients to unravel its mutational patterns.

Sexually dimorphic DNA damage responses and mutation avoidance in the mouse germline

In this study, Bloom and Schimenti examine the response of primordial germ cells to DNA damage. Using both environmental and genetic stresses, the authors reveal the importance of the G1 checkpoint in preventing accumulation of complex mutations in the germline, and the differentiation of the DNA damage response during germ cell development.

Germ cells specified during fetal development form the foundation of the mammalian germline. These primordial germ cells (PGCs) undergo rapid proliferation, yet the germline is highly refractory to mutation accumulation compared with somatic cells. Importantly, while the presence of endogenous or exogenous DNA damage has the potential to impact PGCs, there is little known about how these cells respond to stressors. To better understand the DNA damage response (DDR) in these cells, we exposed pregnant mice to ionizing radiation (IR) at specific gestational time points and assessed the DDR in PGCs. Our results show that PGCs prior to sex determination lack a G1 cell cycle checkpoint. Additionally, the response to IR-induced DNA damage differs between female and male PGCs post-sex determination. IR of female PGCs caused uncoupling of germ cell differentiation and meiotic initiation, while male PGCs exhibited repression of piRNA metabolism and transposon derepression. We also used whole-genome single-cell DNA sequencing to reveal that genetic rescue of DNA repair-deficient germ cells (Fancm^−/−) leads to increased mutation incidence and biases. Importantly, our work uncovers novel insights into how PGCs exposed to DNA damage can become developmentally defective, leaving only those genetically fit cells to establish the adult germline.

Landscape of multi-nucleotide variants in 125,748 human exomes and 15,708 genomes

Multi-nucleotide variants (MNVs), defined as two or more nearby variants existing on the same haplotype in an individual, are a clinically and biologically important class of genetic variation. However, existing tools typically do not accurately classify MNVs, and understanding of their mutational origins remains limited. Here, we systematically survey MNVs in 125,748 whole exomes and 15,708 whole genomes from the Genome Aggregation Database (gnomAD). We identify 1,792,248 MNVs across the genome with constituent variants falling within 2 bp distance of one another, including 18,756 variants with a novel combined effect on protein sequence. Finally, we estimate the relative impact of known mutational mechanisms - CpG deamination, replication error by polymerase zeta, and polymerase slippage at repeat junctions - on the generation of MNVs. Our results demonstrate the value of haplotype-aware variant annotation, and refine our understanding of genome-wide mutational mechanisms of MNVs.

Crossing fitness valleys via double substitutions within codons

BACKGROUND

Single nucleotide substitutions in protein-coding genes can be divided into synonymous (S), with little fitness effect, and non-synonymous (N) ones that alter amino acids and thus generally have a greater effect. Most of the N substitutions are affected by purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases potentially could alleviate the deleterious effect of single substitutions, making them subject to positive selection. To elucidate the effects of selection on double substitutions in all codons, it is critical to differentiate selection from mutational biases.

RESULTS

We addressed the evolutionary regimes of within-codon double substitutions in 37 groups of closely related prokaryotic genomes from diverse phyla by comparing the fractions of double substitutions within codons to those of the equivalent double S substitutions in adjacent codons. Under the assumption that substitutions occur one at a time, all within-codon double substitutions can be represented as "ancestral-intermediate-final" sequences (where "intermediate" refers to the first single substitution and "final" refers to the second substitution) and can be partitioned into four classes: (1) SS, S intermediate-S final; (2) SN, S intermediate-N final; (3) NS, N intermediate-S final; and (4) NN, N intermediate-N final. We found that the selective pressure on the second substitution markedly differs among these classes of double substitutions. Analogous to single S (synonymous) substitutions, SS double substitutions evolve neutrally, whereas analogous to single N (non-synonymous) substitutions, SN double substitutions are subject to purifying selection. In contrast, NS show positive selection on the second step because the original amino acid is recovered. The NN double substitutions are heterogeneous and can be subject to either purifying or positive selection, or evolve neutrally, depending on the amino acid similarity between the final or intermediate and the ancestral states.

CONCLUSIONS

The results of the present, comprehensive analysis of the evolutionary landscape of within-codon double substitutions reaffirm the largely conservative regime of protein evolution. However, the second step of a double substitution can be subject to positive selection when the first step is deleterious. Such positive selection can result in frequent crossing of valleys on the fitness landscape.

Pol V-Mediated Translesion Synthesis Elicits Localized Untargeted Mutagenesis during Post-replicative Gap Repair

In vivo, replication forks proceed beyond replication-blocking lesions by way of downstream repriming, generating daughter strand gaps that are subsequently processed by post-replicative repair pathways such as homologous recombination and translesion synthesis (TLS). The way these gaps are filled during TLS is presently unknown. The structure of gap repair synthesis was assessed by sequencing large collections of single DNA molecules that underwent specific TLS events in vivo. The higher error frequency of specialized relative to replicative polymerases allowed us to visualize gap-filling events at high resolution. Unexpectedly, the data reveal that a specialized polymerase, Pol V, synthesizes stretches of DNA both upstream and downstream of a site-specific DNA lesion. Pol V-mediated untargeted mutations are thus spread over several hundred nucleotides, strongly eliciting genetic instability on either side of a given lesion. Consequently, post-replicative gap repair may be a source of untargeted mutations critical for gene diversification in adaptation and evolution.

The absence of the catalytic domains of Saccharomyces cerevisiae DNA polymerase ϵ strongly reduces DNA replication fidelity

The four B-family DNA polymerases α, δ, ϵ and ζ cooperate to accurately replicate the eukaryotic nuclear genome. Here, we report that a Saccharomyces cerevisiae strain encoding the pol2-16 mutation that lacks Pol ϵ's polymerase and exonuclease activities has increased dNTP concentrations and an increased mutation rate at the CAN1 locus compared to wild type yeast. About half of this mutagenesis disappears upon deleting the REV3 gene encoding the catalytic subunit of Pol ζ. The remaining, still strong, mutator phenotype is synergistically elevated in an msh6Δ strain and has a mutation spectrum characteristic of mistakes made by Pol δ. The results support a model wherein slow-moving replication forks caused by the lack of Pol ϵ's catalytic domains result in greater involvement of mutagenic DNA synthesis by Pol ζ as well as diminished proofreading by Pol δ during replication.

Exome-wide assessment of the functional impact and pathogenicity of multinucleotide mutations

Approximately 2% of de novo single-nucleotide variants (SNVs) appear as part of clustered mutations that create multinucleotide variants (MNVs). MNVs are an important source of genomic variability as they are more likely to alter an encoded protein than a SNV, which has important implications in disease as well as evolution. Previous studies of MNVs have focused on their mutational origins and have not systematically evaluated their functional impact and contribution to disease. We identified 69,940 MNVs and 91 de novo MNVs in 6688 exome-sequenced parent–offspring trios from the Deciphering Developmental Disorders Study comprising families with severe developmental disorders. We replicated the previously described MNV mutational signatures associated with DNA polymerase zeta, an error-prone translesion polymerase, and the APOBEC family of DNA deaminases. We estimate the simultaneous MNV germline mutation rate to be 1.78 × 10⁻¹⁰ mutations per base pair per generation. We found that most MNVs within a single codon create a missense change that could not have been created by a SNV. MNV-induced missense changes were, on average, more physicochemically divergent, were more depleted in highly constrained genes (pLI ≥ 0.9), and were under stronger purifying selection compared with SNV-induced missense changes. We found that de novo MNVs were significantly enriched in genes previously associated with developmental disorders in affected children. This shows that MNVs can be more damaging than SNVs even when both induce missense changes, and are an important variant type to consider in relation to human disease.

R-loop formation by dCas9 is mutagenic in Saccharomyces cerevisiae

Cas9 binds and cleaves specific DNA sequences by inducing the formation of an R-loop between the guide RNA and its genomic target. While targeting of active Cas9 to a genomic locus is highly mutagenic because Cas9 creates DNA double strand breaks, targeting of dead Cas9 (dCas9) is presumed not to be mutagenic, as dCas9 lacks DNA endonuclease activity. Here, we show that dCas9 targeting induces mutations in yeast, particularly when targeted to the non-transcribed strand of a gene. dCas9-induced mutations cluster near the guide RNA target region and are comprised of single nucleotide substitutions, small insertions and deletions, and even complex mutations, depending upon the particular guide RNA target. We show that many of these mutations are a consequence of cytosine deamination events occurring on the non-target strand of the dCas9-induced R-loop, while others are associated with homopolymer instability or translesion DNA synthesis. Targeting of dCas9 by a mismatch-containing guide RNA also increases CAN1 mutation frequency, particularly in an ung1Δ mutant strain, suggesting that dCas9 induces mutations through similar mechanisms at off-target sites. These findings indicate that DNA binding by dCas9 is mutagenic in yeast, likely because dCas9 induces the formation of an R-loop at its target site.

Lack of G1/S control destabilizes the yeast genome via replication stress-induced DSBs and illegitimate recombination

The protein Swi6 in Saccharomyces cerevisiae is a cofactor in two complexes that regulate the transcription of the genes controlling the G1/S transition. It also ensures proper oxidative and cell wall stress responses. Previously, we found that Swi6 was crucial for the survival of genotoxic stress. Here, we show that a lack of Swi6 causes replication stress leading to double-strand break (DSB) formation, inefficient DNA repair and DNA content alterations, resulting in high cell mortality. Comparative genome hybridization experiments revealed that there was a random genome rearrangement in swi6Δ cells, whereas in diploid swi6Δ/swi6Δ cells, chromosome V is duplicated. SWI4 and PAB1, which are located on chromosome V and are known multicopy suppressors of swi6Δ phenotypes, partially reverse swi6Δ genome instability when overexpressed. Another gene on chromosome V, RAD51, also supports swi6Δ survival, but at a high cost; Rad51-dependent illegitimate recombination in swi6Δ cells appears to connect DSBs, leading to genome rearrangement and preventing cell death.This article has an associated First Person interview with the first author of the paper.

Genome Instability Is Promoted by the Chromatin-Binding Protein Spn1 in Saccharomyces cerevisiae

Cells expend a large amount of energy to maintain their DNA sequence. DNA repair pathways, cell cycle checkpoint activation, proofreading polymerases, and chromatin structure are ways in which the cell minimizes changes to the genome. During replication, the DNA-damage tolerance pathway allows the replication forks to bypass damage on the template strand. This avoids prolonged replication fork stalling, which can contribute to genome instability. The DNA-damage tolerance pathway includes two subpathways: translesion synthesis and template switch. Post-translational modification of PCNA and the histone tails, cell cycle phase, and local DNA structure have all been shown to influence subpathway choice. Chromatin architecture contributes to maintaining genome stability by providing physical protection of the DNA and by regulating DNA-processing pathways. As such, chromatin-binding factors have been implicated in maintaining genome stability. Using Saccharomyces cerevisiae, we examined the role of Spn1 (Suppresses postrecruitment gene number 1), a chromatin-binding and transcription elongation factor, in DNA-damage tolerance. Expression of a mutant allele of SPN1 results in increased resistance to the DNA-damaging agent methyl methanesulfonate, lower spontaneous and damage-induced mutation rates, along with increased chronological life span. We attribute these effects to an increased usage of the template switch branch of the DNA-damage tolerance pathway in the spn1 strain. This provides evidence for a role of wild-type Spn1 in promoting genome instability, as well as having ties to overcoming replication stress and contributing to chronological aging.

Multinucleotide mutations cause false inferences of lineage-specific positive selection

Phylogenetic tests of adaptive evolution, such as the widely used branch-site test, assume that nucleotide substitutions occur singly and independently. But recent research has shown that errors at adjacent sites often occur during DNA replication, and the resulting multinucleotide mutations (MNMs) are overwhelmingly likely to be nonsynonymous. We evaluated whether the branch-site test (BST) might misinterpret sequence patterns produced by MNMs as false support for positive selection. We analyzed two genome-scale datasets– one from mammals and one from flies – and found that codons with multiple differences account for virtually all the support for lineage-specific positive selection in the BST. Simulations under conditions derived from these alignments but without positive selection show that realistic rates of MNMs cause a strong and systematic bias towards false inferences of selection. This bias is sufficient under empirically derived conditions to produce false positive inferences as often as the branch-site test infers positive selection from the empirical data. Although some genes with BST-positive results may have evolved adaptively, the test cannot distinguish sequence patterns produced by authentic positive selection from those caused by neutral fixation of MNMs. Many published inferences of adaptive evolution using this technique may therefore be artifacts of model violation caused by unincorporated neutral mutational processes. We introduce a model that incorporates MNMs and may help to ameliorate this bias.

Reading the genome like a history book

Mathematical modeling sheds light on the evolution of human genetic variation

Hypermutation signature reveals a slippage and realignment model of translesion synthesis by Rev3 polymerase in cisplatin-treated yeast

Gene-gene or gene-drug interactions are typically quantified using fitness as a readout because the data are continuous and easily measured in high throughput. However, to what extent fitness captures the range of other phenotypes that show synergistic effects is usually unknown. Using Saccharomyces cerevisiae and focusing on a matrix of DNA repair mutants and genotoxic drugs, we quantify 76 gene-drug interactions based on both mutation rate and fitness and find that these parameters are not connected. Independent of fitness defects, we identified six cases of synthetic hypermutation, where the combined effect of the drug and mutant on mutation rate was greater than predicted. One example occurred when yeast lacking RAD1 were exposed to cisplatin, and we characterized this interaction using whole-genome sequencing. Our sequencing results indicate mutagenesis by cisplatin in rad1Δ cells appeared to depend almost entirely on interstrand cross-links at GpCpN motifs. Interestingly, our data suggest that the following base on the template strand dictates the addition of the mutated base. This result differs from cisplatin mutation signatures in XPF-deficient Caenorhabditis elegans and supports a model in which translesion synthesis polymerases perform a slippage and realignment extension across from the damaged base. Accordingly, DNA polymerase ζ activity was essential for mutagenesis in cisplatin-treated rad1Δ cells. Together these data reveal the potential to gain new mechanistic insights from nonfitness measures of gene-drug interactions and extend the use of mutation accumulation and whole-genome sequencing analysis to define DNA repair mechanisms.

Multi-nucleotide de novo Mutations in Humans

Mutation of the DNA molecule is one of the most fundamental processes in biology. In this study, we use 283 parent-offspring trios to estimate the rate of mutation for both single nucleotide variants (SNVs) and short length variants (indels) in humans and examine the mutation process. We found 17812 SNVs, corresponding to a mutation rate of 1.29 × 10-8 per position per generation (PPPG) and 1282 indels corresponding to a rate of 9.29 × 10-10 PPPG. We estimate that around 3% of human de novo SNVs are part of a multi-nucleotide mutation (MNM), with 558 (3.1%) of mutations positioned less than 20kb from another mutation in the same individual (median distance of 525bp). The rate of de novo mutations is greater in late replicating regions (p = 8.29 × 10-19) and nearer recombination events (p = 0.0038) than elsewhere in the genome.

Clusters of Multiple Mutations: Incidence and Molecular Mechanisms

It has been long understood that mutation distribution is not completely random across genomic space and in time. Indeed, recent surprising discoveries identified multiple simultaneous mutations occurring in tiny regions within chromosomes while the rest of the genome remains relatively mutation-free. Mechanistic elucidation of these phenomena, called mutation showers, mutation clusters, or kataegis, in parallel with findings of abundant clustered mutagenesis in cancer genomes, is ongoing. So far, the combination of factors most important for clustered mutagenesis is the induction of DNA lesions within unusually long and persistent single-strand DNA intermediates. In addition to being a fascinating phenomenon, clustered mutagenesis also became an indispensable tool for identifying a previously unrecognized major source of mutation in cancer, APOBEC cytidine deaminases. Future research on clustered mutagenesis may shed light onto important mechanistic details of genome maintenance, with potentially profound implications for human health.

DNA polymerase ι functions in the generation of tandem mutations during somatic hypermutation of antibody genes

DNA polymerase ι (Pol ι) is an attractive candidate for somatic hypermutation in antibody genes because of its low fidelity. To identify a role for Pol ι, we analyzed mutations in two strains of mice with deficiencies in the enzyme: 129 mice with negligible expression of truncated Pol ι, and knock-in mice that express full-length Pol ι that is catalytically inactive. Both strains had normal frequencies and spectra of mutations in the variable region, indicating that loss of Pol ι did not change overall mutagenesis. We next examined if Pol ι affected tandem mutations generated by another error-prone polymerase, Pol ζ. The frequency of contiguous mutations was analyzed using a novel computational model to determine if they occur during a single DNA transaction or during two independent events. Analyses of 2,000 mutations from both strains indicated that Pol ι-compromised mice lost the tandem signature, whereas C57BL/6 mice accumulated significant amounts of double mutations. The results support a model where Pol ι occasionally accesses the replication fork to generate a first mutation, and Pol ζ extends the mismatch with a second mutation.

Histone H2B mono-ubiquitylation maintains genomic integrity at stalled replication forks

Histone modifications play an important role in regulating access to DNA for transcription, DNA repair and DNA replication. A central player in these events is the mono-ubiquitylation of histone H2B (H2Bub1), which has been shown to regulate nucleosome dynamics. Previously, it was shown that H2Bub1 was important for nucleosome assembly onto nascent DNA at active replication forks. In the absence of H2Bub1, incomplete chromatin structures resulted in several replication defects. Here, we report new evidence, which shows that loss of H2Bub1 contributes to genomic instability in yeast. Specifically, we demonstrate that H2Bub1-deficient yeast accumulate mutations at a high frequency under conditions of replicative stress. This phenotype is due to an aberrant DNA Damage Tolerance (DDT) response upon fork stalling. We show that H2Bub1 normally functions to promote error-free translesion synthesis (TLS) mediated by DNA polymerase eta (Polη). Without H2Bub1, DNA polymerase zeta (Polζ) is responsible for a highly mutagenic alternative mechanism. While H2Bub1 does not appear to regulate other DDT pathways, error-free DDT mechanisms are employed by H2Bub1-deficient cells as another means for survival. However, in these instances, the anti-recombinase, Srs2, is essential to prevent the accumulation of toxic HR intermediates that arise in an unconstrained chromatin environment.

Concurrent nucleotide substitution mutations in the human genome are characterized by a significantly decreased transition/transversion ratio

There is accumulating evidence that the number of multiple-nucleotide substitutions (MNS) occurring in closely spaced sites in eukaryotic genomes is significantly higher than would be predicted from the random accumulation of independently generated single-nucleotide substitutions (SNS). Although this excess can in principle be accounted for by the concept of transient hypermutability, a general mutational signature of concurrent MNS mutations has not so far been evident. Employing a dataset (N = 449) of "concurrent" double MNS mutations causing human inherited disease, we have identified just such a mutational signature: concurrently generated double MNS mutations exhibit a >twofold lower transition/transversion ratio (termed RTs/Tv ) than independently generated de novo SNS mutations (<0.80 vs. 2.10; P = 2.69 × 10(-14) ). We replicated this novel finding through a similar analysis employing two double MNS variant datasets with differing abundances of concurrent events (150,521 variants with both substitutions on the same haplotypic lineage vs. 94,875 variants whose component substitutions were on different haplotypic lineages) plus 5,430,874 SNS variants, all being derived from the whole-genome sequencing of seven Chinese individuals. Evaluation of the newly observed mutational signature in diverse contexts provides solid support for the postulated role of translesion synthesis DNA polymerases in transient hypermutability.

Polymerase ζ Activity Is Linked to Replication Timing in Humans: Evidence from Mutational Signatures

Replication timing is an important determinant of germline mutation patterns, with a higher rate of point mutations in late replicating regions. Mechanisms underlying this association remain elusive. One of the suggested explanations is the activity of error-prone DNA polymerases in late-replicating regions. Polymerase zeta (pol ζ), an essential error-prone polymerase biased toward transversions, also has a tendency to produce dinucleotide mutations (DNMs), complex mutational events that simultaneously affect two adjacent nucleotides. Experimental studies have shown that pol ζ is strongly biased toward GC→AA/TT DNMs. Using primate divergence data, we show that the GC→AA/TT pol ζ mutational signature is the most frequent among DNMs, and its rate exceeds the mean rate of other DNM types by a factor of approximately 10. Unlike the overall rate of DNMs, the pol ζ signature drastically increases with the replication time in the human genome. Finally, the pol ζ signature is enriched in transcribed regions, and there is a strong prevalence of GC→TT over GC→AA DNMs on the nontemplate strand, indicating association with transcription. A recurrently occurring GC→TT DNM in HRAS and SOD1 genes causes the Costello syndrome and amyotrophic lateral sclerosis correspondently; we observe an approximately 1 kb long mutation hotspot enriched by transversions near these DNMs in both cases, suggesting a link between these diseases and pol ζ activity. This study uncovers the genomic preferences of pol ζ, shedding light on a novel cause of mutational heterogeneity along the genome.

Extensive de novo mutation rate variation between individuals and across the genome of Chlamydomonas reinhardtii

Describing the process of spontaneous mutation is fundamental for understanding the genetic basis of disease, the threat posed by declining population size in conservation biology, and much of evolutionary biology. Directly studying spontaneous mutation has been difficult, however, because new mutations are rare. Mutation accumulation (MA) experiments overcome this by allowing mutations to build up over many generations in the near absence of natural selection. Here, we sequenced the genomes of 85 MA lines derived from six genetically diverse strains of the green alga Chlamydomonas reinhardtii. We identified 6843 new mutations, more than any other study of spontaneous mutation. We observed sevenfold variation in the mutation rate among strains and that mutator genotypes arose, increasing the mutation rate approximately eightfold in some replicates. We also found evidence for fine-scale heterogeneity in the mutation rate, with certain sequence motifs mutating at much higher rates, and clusters of multiple mutations occurring at closely linked sites. There was little evidence, however, for mutation rate heterogeneity between chromosomes or over large genomic regions of 200 kbp. We generated a predictive model of the mutability of sites based on their genomic properties, including local GC content, gene expression level, and local sequence context. Our model accurately predicted the average mutation rate and natural levels of genetic diversity of sites across the genome. Notably, trinucleotides vary 17-fold in rate between the most and least mutable sites. Our results uncover a rich heterogeneity in the process of spontaneous mutation both among individuals and across the genome.

Genetic instability in budding and fission yeast-sources and mechanisms

Cells are constantly confronted with endogenous and exogenous factors that affect their genomes. Eons of evolution have allowed the cellular mechanisms responsible for preserving the genome to adjust for achieving contradictory objectives: to maintain the genome unchanged and to acquire mutations that allow adaptation to environmental changes. One evolutionary mechanism that has been refined for survival is genetic variation. In this review, we describe the mechanisms responsible for two biological processes: genome maintenance and mutation tolerance involved in generations of genetic variations in mitotic cells of both Saccharomyces cerevisiae and Schizosaccharomyces pombe. These processes encompass mechanisms that ensure the fidelity of replication, DNA lesion sensing and DNA damage response pathways, as well as mechanisms that ensure precision in chromosome segregation during cell division. We discuss various factors that may influence genome stability, such as cellular ploidy, the phase of the cell cycle, transcriptional activity of a particular region of DNA, the proficiency of DNA quality control systems, the metabolic stage of the cell and its respiratory potential, and finally potential exposure to endogenous or environmental stress.

The stability of budding and fission yeast genomes is influenced by two contradictory factors: (1) the need to be fully functional, which is ensured through the replication fidelity pathways of nuclear and mitochondrial genomes through sensing and repairing DNA damage, through precise chromosome segregation during cell division; and (2) the need to acquire changes for adaptation to environmental challenges.

DNA polymerase ζ-dependent lesion bypass in Saccharomyces cerevisiae is accompanied by error-prone copying of long stretches of adjacent DNA

Translesion synthesis (TLS) helps cells to accomplish chromosomal replication in the presence of unrepaired DNA lesions. In eukaryotes, the bypass of most lesions involves a nucleotide insertion opposite the lesion by either a replicative or a specialized DNA polymerase, followed by extension of the resulting distorted primer terminus by DNA polymerase ζ (Polζ). The subsequent events leading to disengagement of the error-prone Polζ from the primer terminus and its replacement with an accurate replicative DNA polymerase remain largely unknown. As a first step toward understanding these events, we aimed to determine the length of DNA stretches synthesized in an error-prone manner during the Polζ-dependent lesion bypass. We developed new in vivo assays to identify the products of mutagenic TLS through a plasmid-borne tetrahydrofuran lesion and a UV-induced chromosomal lesion. We then surveyed the region downstream of the lesion site (in respect to the direction of TLS) for the presence of mutations indicative of an error-prone polymerase activity. The bypass of both lesions was associated with an approximately 300,000-fold increase in the mutation rate in the adjacent DNA segment, in comparison to the mutation rate during normal replication. The hypermutated tract extended 200 bp from the lesion in the plasmid-based assay and as far as 1 kb from the lesion in the chromosome-based assay. The mutation rate in this region was similar to the rate of errors produced by purified Polζ during copying of undamaged DNA in vitro. Further, no mutations downstream of the lesion were observed in rare TLS products recovered from Polζ-deficient cells. This led us to conclude that error-prone Polζ synthesis continues for several hundred nucleotides after the lesion bypass is completed. These results provide insight into the late steps of TLS and show that error-prone TLS tracts span a substantially larger region than previously appreciated.

Genomic instability is associated with multiple genetic diseases. Endogenous and exogenous DNA-damaging factors constitute a major source of genomic instability. Mutations occur when DNA lesions are bypassed by specialized translesion synthesis (TLS) DNA polymerases that are less accurate than the normal replicative polymerases. The discovery of the remarkable infidelity of the TLS enzymes at the turn of the century immediately suggested that their contribution to replication must be tightly restricted to sites of DNA damage to avoid excessive mutagenesis. The actual extent of error-prone synthesis that accompanies TLS in vivo has never been estimated. We describe a novel genetic approach to measure the length of DNA synthesized by TLS polymerases upon their recruitment to sites of DNA damage. We show that stretches of error-prone synthesis associated with the bypass of a single damaged nucleotide span at least 200 and sometimes up to 1,000 nucleotide-long segments, resulting in more than a 300,000-fold increase in mutagenesis in the surrounding region. We speculate that processive synthesis of long DNA stretches by error-prone polymerases could contribute to clustered mutagenesis, a phenomenon that allows for rapid genome changes without significant loss of fitness and plays an important role in tumorigenesis, the immune response and adaptation.

Phosphorylation of PCNA by EGFR inhibits mismatch repair and promotes misincorporation during DNA synthesis

Proliferating cell nuclear antigen (PCNA) plays essential roles in eukaryotic cells during DNA replication, DNA mismatch repair (MMR), and other events at the replication fork. Earlier studies show that PCNA is regulated by posttranslational modifications, including phosphorylation of tyrosine 211 (Y211) by the epidermal growth factor receptor (EGFR). However, the functional significance of Y211-phosphorylated PCNA remains unknown. Here, we show that PCNA phosphorylation by EGFR alters its interaction with mismatch-recognition proteins MutSα and MutSβ and interferes with PCNA-dependent activation of MutLα endonuclease, thereby inhibiting MMR at the initiation step. Evidence is also provided that Y211-phosphorylated PCNA induces nucleotide misincorporation during DNA synthesis. These findings reveal a novel mechanism by which Y211-phosphorylated PCNA promotes cancer development and progression via facilitating error-prone DNA replication and suppressing the MMR function.

Eukaryotic DNA polymerase ζ

This review focuses on eukaryotic DNA polymerase ζ (Pol ζ), the enzyme responsible for the bulk of mutagenesis in eukaryotic cells in response to DNA damage. Pol ζ is also responsible for a large portion of mutagenesis during normal cell growth, in response to spontaneous damage or to certain DNA structures and other blocks that stall DNA replication forks. Novel insights in mutagenesis have been derived from recent advances in the elucidation of the subunit structure of Pol ζ. The lagging strand DNA polymerase δ shares the small Pol31 and Pol32 subunits with the Rev3-Rev7 core assembly giving a four subunit Pol ζ complex that is the active form in mutagenesis. Furthermore, Pol ζ forms essential interactions with the mutasome assembly factor Rev1 and with proliferating cell nuclear antigen (PCNA). These interactions are modulated by posttranslational modifications such as ubiquitination and phosphorylation that enhance translesion synthesis (TLS) and mutagenesis.

Do mutator mutations fuel tumorigenesis?

The mutator phenotype hypothesis proposes that the mutation rate of normal cells is insufficient to account for the large number of mutations found in human cancers. Consequently, human tumors exhibit an elevated mutation rate that increases the likelihood of a tumor acquiring advantageous mutations. The hypothesis predicts that tumors are composed of cells harboring hundreds of thousands of mutations, as opposed to a small number of specific driver mutations, and that malignant cells within a tumor therefore constitute a highly heterogeneous population. As a result, drugs targeting specific mutated driver genes or even pathways of mutated driver genes will have only limited anticancer potential. In addition, because the tumor is composed of such a diverse cell population, tumor cells harboring drug-resistant mutations will exist prior to the administration of any chemotherapeutic agent. We present recent evidence in support of the mutator phenotype hypothesis, major arguments against this concept, and discuss the clinical consequences of tumor evolution fueled by an elevated mutation rate. We also consider the therapeutic possibility of altering the rate of mutation accumulation. Most significantly, we contend that there is a need to fundamentally reconsider current approaches to personalized cancer therapy. We propose that targeting cellular pathways that alter the rate of mutation accumulation in tumors will ultimately prove more effective than attempting to identify and target mutant driver genes or driver pathways.

Error-prone polymerase activity causes multinucleotide mutations in humans

About 2% of human genetic polymorphisms have been hypothesized to arise via multinucleotide mutations (MNMs), complex events that generate SNPs at multiple sites in a single generation. MNMs have the potential to accelerate the pace at which single genes evolve and to confound studies of demography and selection that assume all SNPs arise independently. In this paper, we examine clustered mutations that are segregating in a set of 1092 human genomes, demonstrating that the signature of MNM becomes enriched as large numbers of individuals are sampled. We estimate the percentage of linked SNP pairs that were generated by simultaneous mutation as a function of the distance between affected sites and show that MNMs exhibit a high percentage of transversions relative to transitions, findings that are reproducible in data from multiple sequencing platforms and cannot be attributed to sequencing error. Among tandem mutations that occur simultaneously at adjacent sites, we find an especially skewed distribution of ancestral and derived alleles, with GC → AA, GA → TT, and their reverse complements making up 27% of the total. These mutations have been previously shown to dominate the spectrum of the error-prone polymerase Pol ζ, suggesting that low-fidelity DNA replication by Pol ζ is at least partly responsible for the MNMs that are segregating in the human population. We develop statistical estimates of MNM prevalence that can be used to correct phylogenetic and population genetic inferences for the presence of complex mutations.

Human Pol ζ purified with accessory subunits is active in translesion DNA synthesis and complements Pol η in cisplatin bypass

DNA polymerase ζ (Pol ζ) is a eukaryotic B-family DNA polymerase that specializes in translesion synthesis and is essential for normal embryogenesis. At a minimum, Pol ζ consists of a catalytic subunit Rev3 and an accessory subunit Rev7. Mammalian Rev3 contains >3,000 residues and is twice as large as the yeast homolog. To date, no vertebrate Pol ζ has been purified for biochemical characterization. Here we report purification of a series of human Rev3 deletion constructs expressed in HEK293 cells and identification of a minimally catalytically active human Pol ζ variant. With a tagged form of an active Pol ζ variant, we isolated two additional accessory subunits of human Pol ζ, PolD2 and PolD3. The purified four-subunit Pol ζ4 (Rev3-Rev7-PolD2-PolD3) is much more efficient and more processive at bypassing a 1,2-intrastrand d(GpG)-cisplatin cross-link than the two-subunit Pol ζ2 (Rev3-Rev7). We show that complete bypass of cisplatin lesions requires Pol η to insert dCTP opposite the 3' guanine and Pol ζ4 to extend the primers.

Patterns and mutational signatures of tandem base substitutions causing human inherited disease

Tandem base substitutions (TBSs) are multiple mutations that comprise two or more contiguous nucleotide substitutions without any net gain or loss of bases. They have recently become recognized as a distinct category of human genomic variant. However, their role in causing human inherited disease so far has not been studied methodically. Here, using data from the Human Gene Mutation Database (http://www.hgmd.org), we identified 477 events to be TBSs (doublets, 448; triplets, 16; and quadruplets to octuplets, 13). A comprehensive sequence pattern and context analysis implied the likely fundamental importance of translesion synthesis (TLS) DNA polymerases in generating these diverse TBSs but revealed that TLS polymerases may operate differently in generating TBSs of ≤ 3 bases (bypass of endogenous DNA lesions) than those of ≥ 4 bases (serial replication slippage). Moreover, GC was found to be the most frequently affected dinucleotide with GC/GC>AA/TT being the most frequent double TBS. Comparison with cancer genome mutational spectra allowed us to conclude that human germline TBSs arise predominantly through the action of endogenous mechanisms of mutagenesis rather than through exposure to exogenous mutagens. Finally, the rates of double and triple TBSs were estimated to be 0.2-1.2 × 10(-10) and 0.8-4.8 × 10(-12) per base per generation, respectively.