Yeast origins establish a strand bias for replicational mutagenesis

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. The biased use of synonymous codons has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. The Saccharomycotina, the fungal subphylum containing the yeasts Saccharomyces cerevisiae and Candida albicans , has been a model system for studying codon usage. We characterized codon usage across 1,154 strains from 1,051 species to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns across the subphylum. We found evidence of a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is also distinct between the 12 orders within the subphylum to such a degree that yeasts can be classified into orders with an accuracy greater than 90% using a machine learning algorithm trained on codon usage. We also characterized the degree to which codon usage bias is impacted by translational selection. Interestingly, the degree of translational selection was influenced by a combination of genome features and assembly metrics that included the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs. The order contains 24 species, and 23 are computationally predicted to lack tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that extreme avoidance of the CGN codons is associated with a decline in arginine tRNA function. Codon usage bias within the Saccharomycotina is generally consistent with previous investigations in fungi, which show a role for both genomic features and GC bias in shaping codon usage. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

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.

Spontaneous deamination of cytosine to uracil is biased to the non-transcribed DNA strand in yeast

Transcription in Saccharomyces cerevisiae is associated with elevated mutation and this partially reflects enhanced damage of the corresponding DNA. Spontaneous deamination of cytosine to uracil leads to CG>TA mutations that provide a strand-specific read-out of damage in strains that lack the ability to remove uracil from DNA. Using the CAN1 forward mutation reporter, we found that C>T and G>A mutations, which reflect deamination of the non-transcribed and transcribed DNA strands, respectively, occurred at similar rates under low-transcription conditions. By contrast, the rate of C>T mutations was 3-fold higher than G>A mutations under high-transcription conditions, demonstrating biased deamination of the non-transcribed strand (NTS). The NTS is transiently single-stranded within the ∼15 bp transcription bubble, or a more extensive region of the NTS can be exposed as part of an R-loop that can form behind RNA polymerase. Neither the deletion of genes whose products restrain R-loop formation nor the over-expression of RNase H1, which degrades R-loops, reduced the biased deamination of the NTS, and no transcription-associated R-loop formation at CAN1 was detected. These results suggest that the NTS within the transcription bubble is a target for spontaneous deamination and likely other types of DNA damage.

Mismatch repair operates at the replication fork in direct competition with mismatch extension by DNA polymerase δ

DNA mismatch repair (MMR) in eukaryotes is believed to occur post-replicatively, wherein nicks or gaps in the nascent DNA strand are suggested to serve as strand discrimination signals. However, how such signals are generated in the nascent leading strand has remained unclear. Here we examine the alternative possibility that MMR occurs in conjunction with the replication fork. To this end, we utilize mutations in the PCNA interacting peptide (PIP) domain of the Pol3 or Pol32 subunit of DNA polymerase δ (Polδ) and show that these pip mutations suppress the greatly elevated mutagenesis in yeast strains harboring the pol3-01 mutation defective in Polδ proofreading activity. And strikingly, they suppress the synthetic lethality of pol3-01 pol2-4 double mutant strains, which arises from the vastly enhanced mutability due to defects in the proofreading functions of both Polδ and Polε. Our finding that suppression of elevated mutagenesis in pol3-01 by the Polδ pip mutations requires intact MMR supports the conclusion that MMR operates at the replication fork in direct competition with other mismatch removal processes and with extension of synthesis from the mispair by Polδ. Furthermore, the evidence that Polδ pip mutations eliminate almost all the mutability of pol2-4 msh2Δ or pol3-01 pol2-4 adds strong support for a major role of Polδ in replication of both the leading and lagging DNA strands.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

Viral genome sequence datasets display pervasive evidence of strand-specific substitution biases that are best described using non-reversible nucleotide substitution models

Background

The vast majority of phylogenetic trees are inferred from molecular sequence data (nucleotides or amino acids) using time-reversible evolutionary models which assume that, for any pair of nucleotide or amino acid characters, the relative rate of X to Y substitution is the same as the relative rate of Y to X substitution. However, this reversibility assumption is unlikely to accurately reflect the actual underlying biochemical and/or evolutionary processes that lead to the fixation of substitutions. Here, we use empirical viral genome sequence data to reveal that evolutionary non-reversibility is pervasive among most groups of viruses. Specifically, we consider two non-reversible nucleotide substitution models: (1) a 6-rate non-reversible model (NREV6) in which Watson-Crick complementary substitutions occur at identical relative rates and which might therefor be most applicable to analyzing the evolution of genomes where both complementary strands are subject to the same mutational processes (such as might be expected for double-stranded (ds) RNA or dsDNA genomes); and (2) a 12-rate non-reversible model (NREV12) in which all relative substitution types are free to occur at different rates and which might therefore be applicable to analyzing the evolution of genomes where the complementary genome strands are subject to different mutational processes (such as might be expected for viruses with single-stranded (ss) RNA or ssDNA genomes).

Results

Using likelihood ratio and Akaike Information Criterion-based model tests, we show that, surprisingly, NREV12 provided a significantly better fit to 21/31 dsRNA and 20/30 dsDNA datasets than did the general time reversible (GTR) and NREV6 models with NREV6 providing a better fit than NREV12 and GTR in only 5/30 dsDNA and 2/31 dsRNA datasets. As expected, NREV12 provided a significantly better fit to 24/33 ssDNA and 40/47 ssRNA datasets. Next, we used simulations to show that increasing degrees of strand-specific substitution bias decrease the accuracy of phylogenetic inference irrespective of whether GTR or NREV12 is used to describe mutational processes. However, in cases where strand-specific substitution biases are extreme (such as in SARS-CoV-2 and Torque teno sus virus datasets) NREV12 tends to yield more accurate phylogenetic trees than those obtained using GTR.

Conclusion

We show that NREV12 should, be seriously considered during the model selection phase of phylogenetic analyses involving viral genomic sequences.

Partners in crime: Tbf1 and Vid22 promote expansions of long human telomeric repeats at an interstitial chromosome position in yeast

In humans, telomeric repeats (TTAGGG)_n are known to be present at internal chromosomal sites. These interstitial telomeric sequences (ITSs) are an important source of genomic instability, including repeat length polymorphism, but the molecular mechanisms responsible for this instability remain to be understood. Here, we studied the mechanisms responsible for expansions of human telomeric (Htel) repeats that were artificially inserted inside a yeast chromosome. We found that Htel repeats in an interstitial chromosome position are prone to expansions. The propensity of Htel repeats to expand depends on the presence of a complex of two yeast proteins: Tbf1 and Vid22. These two proteins are physically bound to an interstitial Htel repeat, and together they slow replication fork progression through it. We propose that slow progression of the replication fork through the protein complex formed by the Tbf1 and Vid22 partners at the Htel repeat cause DNA strand slippage, ultimately resulting in repeat expansions.

2-Dimensional Genetic Algorithm Exhibited an Essentiality of Gene Interaction for Evolution

Organisms consist of several genetic factors differing between species. However, the evolutionary effects of gene interactions on the evolutionary rate, adaptation, and divergence of organisms remain unknown. In a previous study, the 2-dimensional genetic algorithm (2DGA) program, including a gene interaction parameter, could simulate punctuated equilibrium under the disparity mode. Following this, we verified the effect of the number of gene interactions (gene cluster size) on evolution speed, adaptation, and divergence using the advanced 2DGA program. In this program, the population was replicated, mutated, and selected for 200,000 generations, and the fitness score, divergence, number of population, and genotype were output and plotted. The genotype data were used for evaluating the phylogenetic relations among the population. The gene cluster size 1) affected the disparity and parity mutagenesis modes differently, 2) determined the growth/exclusion rate and error threshold, and 3) accelerated or decelerated the population's speed of evolutionary advancement. In particular, when the gene cluster size expanded, the rate of increase in fitness scores decreased independently of the mutation rate and mode of mutation (disparity mode/parity mode). The mutation rate at the error threshold was also decreased by expanding the gene cluster size. Dendrograms traced the genotypes of the simulated population, indicating that the disparity mode caused the evolutionary process to enter 1) a stun mode, 2) an evolution mode, or 3) a divergence mode based on the mutation rate and gene cluster size, while the parity mode did not cause the population to enter a stun mode. Based on the above findings, we compared the predictions of the present study with evolution observed in the laboratory or the natural world and the processes of ongoing virus evolution, suggesting that our findings possibly explained the real evolution.

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.

Recombination and Pol ζ Rescue Defective DNA Replication upon Impaired CMG Helicase-Pol ε Interaction

The CMG complex (Cdc45, Mcm2–7, GINS (Psf1, 2, 3, and Sld5)) is crucial for both DNA replication initiation and fork progression. The CMG helicase interaction with the leading strand DNA polymerase epsilon (Pol ε) is essential for the preferential loading of Pol ε onto the leading strand, the stimulation of the polymerase, and the modulation of helicase activity. Here, we analyze the consequences of impaired interaction between Pol ε and GINS in Saccharomyces cerevisiae cells with the psf1-100 mutation. This significantly affects DNA replication activity measured in vitro, while in vivo, the psf1-100 mutation reduces replication fidelity by increasing slippage of Pol ε, which manifests as an elevated number of frameshifts. It also increases the occurrence of single-stranded DNA (ssDNA) gaps and the demand for homologous recombination. The psf1-100 mutant shows elevated recombination rates and synthetic lethality with rad52Δ. Additionally, we observe increased participation of DNA polymerase zeta (Pol ζ) in DNA synthesis. We conclude that the impaired interaction between GINS and Pol ε requires enhanced involvement of error-prone Pol ζ, and increased participation of recombination as a rescue mechanism for recovery of impaired replication forks.

Mismatch repair and DNA polymerase δ proofreading prevent catastrophic accumulation of leading strand errors in cells expressing a cancer-associated DNA polymerase ϵ variant

Substitutions in the exonuclease domain of DNA polymerase ϵ cause ultramutated human tumors. Yeast and mouse mimics of the most common variant, P286R, produce mutator effects far exceeding the effect of Polϵ exonuclease deficiency. Yeast Polϵ-P301R has increased DNA polymerase activity, which could underlie its high mutagenicity. We aimed to understand the impact of this increased activity on the strand-specific role of Polϵ in DNA replication and the action of extrinsic correction systems that remove Polϵ errors. Using mutagenesis reporters spanning a well-defined replicon, we show that both exonuclease-deficient Polϵ (Polϵ-exo⁻) and Polϵ-P301R generate mutations in a strictly strand-specific manner, yet Polϵ-P301R is at least ten times more mutagenic than Polϵ-exo⁻ at each location analyzed. Thus, the cancer variant remains a dedicated leading-strand polymerase with markedly low accuracy. We further show that P301R substitution is lethal in strains lacking Polδ proofreading or mismatch repair (MMR). Heterozygosity for pol2-P301R is compatible with either defect but causes strong synergistic increases in the mutation rate, indicating that Polϵ-P301R errors are corrected by Polδ proofreading and MMR. These data reveal the unexpected ease with which polymerase exchange occurs in vivo, allowing Polδ exonuclease to prevent catastrophic accumulation of Polϵ-P301R-generated errors on the leading strand.

Defects in the GINS complex increase the instability of repetitive sequences via a recombination-dependent mechanism

Faithful replication and repair of DNA lesions ensure genome maintenance. During replication in eukaryotic cells, DNA is unwound by the CMG helicase complex, which is composed of three major components: the Cdc45 protein, Mcm2-7, and the GINS complex. The CMG in complex with DNA polymerase epsilon (CMG-E) participates in the establishment and progression of the replisome. Impaired functioning of the CMG-E was shown to induce genomic instability and promote the development of various diseases. Therefore, CMG-E components play important roles as caretakers of the genome. In Saccharomyces cerevisiae, the GINS complex is composed of the Psf1, Psf2, Psf3, and Sld5 essential subunits. The Psf1-1 mutant form fails to interact with Psf3, resulting in impaired replisome assembly and chromosome replication. Here, we show increased instability of repeat tracts (mononucleotide, dinucleotide, trinucleotide and longer) in yeast psf1-1 mutants. To identify the mechanisms underlying this effect, we analyzed repeated sequence instability using derivatives of psf1-1 strains lacking genes involved in translesion synthesis, recombination, or mismatch repair. Among these derivatives, deletion of RAD52, RAD51, MMS2, POL32, or PIF1 significantly decreased DNA repeat instability. These results, together with the observed increased amounts of single-stranded DNA regions and Rfa1 foci suggest that recombinational mechanisms make important contributions to repeat tract instability in psf1-1 cells. We propose that defective functioning of the CMG-E complex in psf1-1 cells impairs the progression of DNA replication what increases the contribution of repair mechanisms such as template switch and break-induced replication. These processes require sequence homology search which in case of a repeated DNA tract may result in misalignment leading to its expansion or contraction.

Processes that ensure genome stability are crucial for all organisms to avoid mutations and decrease the risk of diseases. The coordinated activity of mechanisms underlying the maintenance of high-fidelity DNA duplication and repair is critical to deal with the malfunction of replication forks or DNA damage. Repeated sequences in DNA are particularly prone to instability; these sequences undergo expansions or contractions, leading in humans to various neurological, neurodegenerative, and neuromuscular disorders. A mutant form of one of the noncatalytic subunits of active DNA helicase complex impairs DNA replication. Here, we show that this form also significantly increases the instability of mononucleotide, dinucleotide, trinucleotide and longer repeat tracts. Our results suggest that in cells that harbor a mutated variant of the helicase complex, continuation of DNA replication is facilitated by recombination processes, and this mechanism can be highly mutagenic during repair synthesis through repetitive regions, especially regions that form secondary structures. Our results indicate that proper functioning of the DNA helicase complex is crucial for maintenance of the stability of repeated DNA sequences, especially in the context of recently described disorders in which mutations or deregulation of the human homologs of genes encoding DNA helicase subunits were observed.

Comment on "A commensal strain of Staphylococcus epidermidis protects against skin neoplasia" by Nakatsuji et al

A recent article in Science Advances described the striking discovery that the commensal Staphylococcus epidermidis strain MO34 displays antimicrobial and antitumor activities by producing a small molecule, identified as the nucleobase analog 6-N-hydroxylaminopurine (6-HAP). However, in contradiction to the literature, the authors claimed that 6-HAP is nonmutagenic and proposed that the toxic effect of 6-HAP results from its ability to inhibit, in its base form, DNA synthesis. To resolve the discrepancy, we proved by genetic experiments with bacteria and yeast that extracts of MO34 do contain a mutagenic compound whose effects are identical to chemically synthesized 6-HAP. The MO34 extract induced the same mutation spectrum as authentic 6-HAP. Notably, the toxic and mutagenic effects of both synthetic and MO34-derived 6-HAP depended on conversion to the corresponding nucleotide. The nucleobase 6-HAP does not inhibit DNA synthesis in vitro, and we conclude that 6-HAP exerts its biological activity when incorporated into DNA.

The importance of an interaction network for proper DNA polymerase ζ heterotetramer activity

Precisely controlled mechanisms have been evolved to rescue impeded DNA replication resulting from encountered obstacles and involve a set of low-fidelity translesion synthesis (TLS) DNA polymerases. Studies in recent years have brought new insights into those TLS polymerases, especially concerning the structure and subunit composition of DNA polymerase zeta (Pol ζ). Pol ζ is predominantly involved in induced mutagenesis as well as the bypass of noncanonical DNA structures, and it is proficient in extending from terminal mismatched nucleotides incorporated by major replicative DNA polymerases. Two active forms of Pol ζ, heterodimeric (Pol ζ2) and heterotetrameric (Pol ζ4) ones, have been identified and studied. Here, in the light of recent publications regarding induced and spontaneous mutagenesis and diverse interactions within Pol ζ holoenzyme, combined with Pol ζ binding to the TLS polymerase Rev1p, we discuss the subunit composition of Pol ζ in various cellular physiological conditions. Available data show that it is the heterotetrameric form of Pol ζ that is involved both during spontaneous and induced mutagenesis, and underline the importance of interactions within Pol ζ when an increased Pol ζ recruitment occurs. Understanding Pol ζ function in the bypass of DNA obstacles would give a significant insight into cellular tolerance of DNA damage, genetic instability and the onset of cancer progression.

Measuring deaminated nucleotide surveillance enzyme ITPA activity with an ATP-releasing nucleotide chimera

Nucleotide quality surveillance enzymes play important roles in human health, by detecting damaged molecules in the nucleotide pool and deactivating them before they are incorporated into chromosomal DNA or adversely affect metabolism. In particular, deamination of adenine moiety in (deoxy)nucleoside triphosphates, resulting in formation of (d)ITP, can be deleterious, leading to DNA damage, mutagenesis and other harmful cellular effects. The 21.5 kDa human enzyme that mitigates this damage by conversion of (d)ITP to monophosphate, ITPA, has been proposed as a possible therapeutic and diagnostic target for multiple diseases. Measuring the activity of this enzyme is useful both in basic research and in clinical applications involving this pathway, but current methods are nonselective and are not applicable to measurement of the enzyme from cells or tissues. Here, we describe the design and synthesis of an ITPA-specific chimeric dinucleotide (DIAL) that replaces the pyrophosphate leaving group of the native substrate with adenosine triphosphate, enabling sensitive detection via luciferase luminescence signaling. The probe is shown to function sensitively and selectively to quantify enzyme activity in vitro, and can be used to measure the activity of ITPA in bacterial, yeast and human cell lysates.

Human DNA polymerase delta double-mutant D316A;E318A interferes with DNA mismatch repair in vitro

DNA mismatch repair (MMR) is a highly-conserved DNA repair mechanism, whose primary role is to remove DNA replication errors preventing them from manifesting as mutations, thereby increasing the overall genome stability. Defects in MMR are associated with increased cancer risk in humans and other organisms. Here, we characterize the interaction between MMR and a proofreading-deficient allele of the human replicative DNA polymerase delta, PolδD316A;E318A, which has a higher capacity for strand displacement DNA synthesis than wild type Polδ. Human cell lines overexpressing PolδD316A;E318A display a mild mutator phenotype, while nuclear extracts of these cells exhibit reduced MMR activity in vitro, and these defects are complemented by overexpression or addition of exogenous human Exonuclease 1 (EXO1). By contrast, another proofreading-deficient mutant, PolδD515V, which has a weaker strand displacement activity, does not decrease the MMR activity as significantly as PolδD316A;E318A. In addition, PolδD515V does not increase the mutation frequency in MMR-proficient cells. Based on our findings, we propose that the proofreading activity restricts the strand displacement activity of Polδ in MMR. This contributes to maintain the nicks required for EXO1 entry, and in this manner ensures the dominance of the EXO1-dependent MMR pathway.

The CysB motif of Rev3p involved in the formation of the four-subunit DNA polymerase ζ is required for defective-replisome-induced mutagenesis

Eukaryotic DNA replication is performed by high-fidelity multi-subunit replicative B-family DNA polymerases (Pols) α, δ and ɛ. Those complexes are composed of catalytic and accessory subunits and organized in multicomplex machinery: the replisome. The fourth B-family member, DNA polymerase zeta (Pol ζ), is responsible for a large portion of mutagenesis in eukaryotic cells. Two forms of Pol ζ have been identified, a hetero-dimeric (Pol ζ₂ ) and a hetero-tetrameric (Pol ζ₄ ) ones and recent data have demonstrated that Pol ζ₄ is responsible for damage-induced mutagenesis. Here, using yeast Pol ζ mutant defective in the assembly of the Pol ζ four-subunit form, we show in vivo that [4Fe-4S] cluster in Pol ζ catalytic subunit (Rev3p) is also required for spontaneous (wild-type cells) and defective-replisome-induced mutagenesis - DRIM (pol3-Y708A, pol2-1 or psf1-100 cells), when cells are not treated with any external damaging agents.

DNA mismatch repair and its many roles in eukaryotic cells

DNA mismatch repair (MMR) is an important DNA repair pathway that plays critical roles in DNA replication fidelity, mutation avoidance and genome stability, all of which contribute significantly to the viability of cells and organisms. MMR is widely-used as a diagnostic biomarker for human cancers in the clinic, and as a biomarker of cancer susceptibility in animal model systems. Prokaryotic MMR is well-characterized at the molecular and mechanistic level; however, MMR is considerably more complex in eukaryotic cells than in prokaryotic cells, and in recent years, it has become evident that MMR plays novel roles in eukaryotic cells, several of which are not yet well-defined or understood. Many MMR-deficient human cancer cells lack mutations in known human MMR genes, which strongly suggests that essential eukaryotic MMR components/cofactors remain unidentified and uncharacterized. Furthermore, the mechanism by which the eukaryotic MMR machinery discriminates between the parental (template) and the daughter (nascent) DNA strand is incompletely understood and how cells choose between the EXO1-dependent and the EXO1-independent subpathways of MMR is not known. This review summarizes recent literature on eukaryotic MMR, with emphasis on the diverse cellular roles of eukaryotic MMR proteins, the mechanism of strand discrimination and cross-talk/interactions between and co-regulation of MMR and other DNA repair pathways in eukaryotic cells. The main conclusion of the review is that MMR proteins contribute to genome stability through their ability to recognize and promote an appropriate cellular response to aberrant DNA structures, especially when they arise during DNA replication. Although the molecular mechanism of MMR in the eukaryotic cell is still not completely understood, increased used of single-molecule analyses in the future may yield new insight into these unsolved questions.

Context based computational analysis and characterization of ARS consensus sequences (ACS) of Saccharomyces cerevisiae genome

Genome-wide experimental studies in Saccharomyces cerevisiae reveal that autonomous replicating sequence (ARS) requires an essential consensus sequence (ACS) for replication activity. Computational studies identified thousands of ACS like patterns in the genome. However, only a few hundreds of these sites act as replicating sites and the rest are considered as dormant or evolving sites. In a bid to understand the sequence makeup of replication sites, a content and context-based analysis was performed on a set of replicating ACS sequences that binds to origin-recognition complex (ORC) denoted as ORC-ACS and non-replicating ACS sequences (nrACS), that are not bound by ORC. In this study, DNA properties such as base composition, correlation, sequence dependent thermodynamic and DNA structural profiles, and their positions have been considered for characterizing ORC-ACS and nrACS. Analysis reveals that ORC-ACS depict marked differences in nucleotide composition and context features in its vicinity compared to nrACS. Interestingly, an A-rich motif was also discovered in ORC-ACS sequences within its nucleosome-free region. Profound changes in the conformational features, such as DNA helical twist, inclination angle and stacking energy between ORC-ACS and nrACS were observed. Distribution of ACS motifs in the non-coding segments points to the locations of ORC-ACS which are found far away from the adjacent gene start position compared to nrACS thereby enabling an accessible environment for ORC-proteins. Our attempt is novel in considering the contextual view of ACS and its flanking region along with nucleosome positioning in the S. cerevisiae genome and may be useful for any computational prediction scheme.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

Quantifying the contributions of base selectivity, proofreading and mismatch repair to nuclear DNA replication in Saccharomyces cerevisiae

Mismatches generated during eukaryotic nuclear DNA replication are removed by two evolutionarily conserved error correction mechanisms acting in series, proofreading and mismatch repair (MMR). Defects in both processes are associated with increased susceptibility to cancer. To better understand these processes, we have quantified base selectivity, proofreading and MMR during nuclear DNA replication in Saccharomyces cerevisiae. In the absence of proofreading and MMR, the primary leading and lagging strand replicases, polymerase ɛ and polymerase δ respectively, synthesize DNA in vivo with somewhat different error rates and specificity, and with apparent base selectivity that is more than 100 times higher than measured in vitro. Moreover, leading and lagging strand replication fidelity rely on a different balance between proofreading and MMR. On average, proofreading contributes more to replication fidelity than does MMR, but their relative contributions vary from nearly all proofreading of some mismatches to mostly MMR of other mismatches. Thus accurate replication of the two DNA strands results from a non-uniform and variable balance between error prevention, proofreading and MMR.

Evidence that processing of ribonucleotides in DNA by topoisomerase 1 is leading-strand specific

Ribonucleotides incorporated during DNA replication are removed by RNase H2-dependent ribonucleotide excision repair (RER). In RER-defective yeast, topoisomerase 1 (Top1) incises DNA at unrepaired ribonucleotides, initiating their removal, but this is accompanied by RNA-DNA-damage phenotypes. Here we show that these phenotypes are incurred by a high level of ribonucleotides incorporated by a leading strand-replicase variant, DNA polymerase (Pol) ɛ, but not by orthologous variants of the lagging-strand replicases, Pols α or δ. Moreover, loss of both RNases H1 and H2 is lethal in combination with increased ribonucleotide incorporation by Pol ɛ but not by Pols α or δ. Several explanations for this asymmetry are considered, including the idea that Top1 incision at ribonucleotides relieves torsional stress in the nascent leading strand but not in the nascent lagging strand, in which preexisting nicks prevent the accumulation of superhelical tension.

Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer

Recent sequencing studies have extensively explored the somatic alterations present in the nuclear genomes of cancers. Although mitochondria control energy metabolism and apoptosis, the origins and impact of cancer-associated mutations in mtDNA are unclear. In this study, we analyzed somatic alterations in mtDNA from 1675 tumors. We identified 1907 somatic substitutions, which exhibited dramatic replicative strand bias, predominantly C > T and A > G on the mitochondrial heavy strand. This strand-asymmetric signature differs from those found in nuclear cancer genomes but matches the inferred germline process shaping primate mtDNA sequence content. A number of mtDNA mutations showed considerable heterogeneity across tumor types. Missense mutations were selectively neutral and often gradually drifted towards homoplasmy over time. In contrast, mutations resulting in protein truncation undergo negative selection and were almost exclusively heteroplasmic. Our findings indicate that the endogenous mutational mechanism has far greater impact than any other external mutagens in mitochondria and is fundamentally linked to mtDNA replication.

Proper functioning of the GINS complex is important for the fidelity of DNA replication in yeast

The role of replicative DNA polymerases in ensuring genome stability is intensively studied, but the role of other components of the replisome is still not fully understood. One of such component is the GINS complex (comprising the Psf1, Psf2, Psf3 and Sld5 subunits), which participates in both initiation and elongation of DNA replication. Until now, the understanding of the physiological role of GINS mostly originated from biochemical studies. In this article, we present genetic evidence for an essential role of GINS in the maintenance of replication fidelity in Saccharomyces cerevisiae. In our studies we employed the psf1-1 allele (Takayama et al., 2003) and a novel psf1-100 allele isolated in our laboratory. Analysis of the levels and specificity of mutations in the psf1 strains indicates that the destabilization of the GINS complex or its impaired interaction with DNA polymerase epsilon increases the level of spontaneous mutagenesis and the participation of the error-prone DNA polymerase zeta. Additionally, a synergistic mutator effect was found for the defects in Psf1p and in the proofreading activity of Pol epsilon, suggesting that proper functioning of GINS is crucial for facilitating error-free processing of terminal mismatches created by Pol epsilon.

Replicative DNA polymerase mutations in cancer

Three DNA polymerases - Pol α, Pol δ and Pol ɛ - are essential for DNA replication. After initiation of DNA synthesis by Pol α, Pol δ or Pol ɛ take over on the lagging and leading strand respectively. Pol δ and Pol ɛ perform the bulk of replication with very high fidelity, which is ensured by Watson-Crick base pairing and 3'exonuclease (proofreading) activity. Yeast models have shown that mutations in the exonuclease domain of Pol δ and Pol ɛ homologues can cause a mutator phenotype. Recently, we identified germline exonuclease domain mutations (EDMs) in human POLD1 and POLE that predispose to 'polymerase proofreading associated polyposis' (PPAP), a disease characterised by multiple colorectal adenomas and carcinoma, with high penetrance and dominant inheritance. Moreover, somatic EDMs in POLE have also been found in sporadic colorectal and endometrial cancers. Tumors with EDMs are microsatellite stable and show an 'ultramutator' phenotype, with a dramatic increase in base substitutions.

Genome-wide mutation avalanches induced in diploid yeast cells by a base analog or an APOBEC deaminase

Genetic information should be accurately transmitted from cell to cell; conversely, the adaptation in evolution and disease is fueled by mutations. In the case of cancer development, multiple genetic changes happen in somatic diploid cells. Most classic studies of the molecular mechanisms of mutagenesis have been performed in haploids. We demonstrate that the parameters of the mutation process are different in diploid cell populations. The genomes of drug-resistant mutants induced in yeast diploids by base analog 6-hydroxylaminopurine (HAP) or AID/APOBEC cytosine deaminase PmCDA1 from lamprey carried a stunning load of thousands of unselected mutations. Haploid mutants contained almost an order of magnitude fewer mutations. To explain this, we propose that the distribution of induced mutation rates in the cell population is uneven. The mutants in diploids with coincidental mutations in the two copies of the reporter gene arise from a fraction of cells that are transiently hypersensitive to the mutagenic action of a given mutagen. The progeny of such cells were never recovered in haploids due to the lethality caused by the inactivation of single-copy essential genes in cells with too many induced mutations. In diploid cells, the progeny of hypersensitive cells survived, but their genomes were saturated by heterozygous mutations. The reason for the hypermutability of cells could be transient faults of the mutation prevention pathways, like sanitization of nucleotide pools for HAP or an elevated expression of the PmCDA1 gene or the temporary inability of the destruction of the deaminase. The hypothesis on spikes of mutability may explain the sudden acquisition of multiple mutational changes during evolution and carcinogenesis.

In vivo bypass of 8-oxodG

8-oxoG is one of the most common and mutagenic DNA base lesions caused by oxidative damage. However, it has not been possible to study the replication of a known 8-oxoG base in vivo in order to determine the accuracy of its replication, the influence of various components on that accuracy, and the extent to which an 8-oxoG might present a barrier to replication. We have been able to place a single 8-oxoG into the Saccharomyces cerevisiae chromosome in a defined location using single-strand oligonucleotide transformation and to study its replication in a fully normal chromosome context. During replication, 8-oxoG is recognized as a lesion and triggers a switch to translesion synthesis by Pol η, which replicates 8-oxoG with an accuracy (insertion of a C opposite the 8-oxoG) of approximately 94%. In the absence of Pol η, template switching to the newly synthesized sister chromatid is observed at least one third of the time; replication of the 8-oxoG in the absence of Pol η is less than 40% accurate. The mismatch repair (MMR) system plays an important role in 8-oxoG replication. Template switching is blocked by MMR and replication accuracy even in the absence of Pol η is approximately 95% when MMR is active. These findings indicate that in light of the overlapping mechanisms by which errors in 8-oxoG replication can be avoided in the cell, the mutagenic threat of 8-oxoG is due more to its abundance than the effect of a single lesion. In addition, the methods used here should be applicable to the study of any lesion that can be stably incorporated into synthetic oligonucleotides.

Ribonucleotides are signals for mismatch repair of leading-strand replication errors

To maintain genome stability, mismatch repair of nuclear DNA replication errors must be directed to the nascent strand, likely by DNA ends and PCNA. Here we show that the efficiency of mismatch repair in Saccharomyces cerevisiae is reduced by inactivating RNase H2, which nicks DNA containing ribonucleotides incorporated during replication. In strains encoding mutator polymerases, this reduction is preferential for repair of mismatches made by leading-strand DNA polymerase ε as compared to lagging-strand DNA polymerase δ. The results suggest that RNase-H2-dependent processing of ribonucleotides transiently present in DNA after replication may direct mismatch repair to the continuously replicated nascent leading strand.

Emergence of DNA polymerase ε antimutators that escape error-induced extinction in yeast

DNA polymerases (Pols) ε and δ perform the bulk of yeast leading- and lagging-strand DNA synthesis. Both Pols possess intrinsic proofreading exonucleases that edit errors during polymerization. Rare errors that elude proofreading are extended into duplex DNA and excised by the mismatch repair (MMR) system. Strains that lack Pol proofreading or MMR exhibit a 10- to 100-fold increase in spontaneous mutation rate (mutator phenotype), and inactivation of both Pol δ proofreading (pol3-01) and MMR is lethal due to replication error-induced extinction (EEX). It is unclear whether a similar synthetic lethal relationship exists between defects in Pol ε proofreading (pol2-4) and MMR. Using a plasmid-shuffling strategy in haploid Saccharomyces cerevisiae, we observed synthetic lethality of pol2-4 with alleles that completely abrogate MMR (msh2Δ, mlh1Δ, msh3Δ msh6Δ, or pms1Δ mlh3Δ) but not with partial MMR loss (msh3Δ, msh6Δ, pms1Δ, or mlh3Δ), indicating that high levels of unrepaired Pol ε errors drive extinction. However, variants that escape this error-induced extinction (eex mutants) frequently emerged. Five percent of pol2-4 msh2Δ eex mutants encoded second-site changes in Pol ε that reduced the pol2-4 mutator phenotype between 3- and 23-fold. The remaining eex alleles were extragenic to pol2-4. The locations of antimutator amino-acid changes in Pol ε and their effects on mutation spectra suggest multiple mechanisms of mutator suppression. Our data indicate that unrepaired leading- and lagging-strand polymerase errors drive extinction within a few cell divisions and suggest that there are polymerase-specific pathways of mutator suppression. The prevalence of suppressors extragenic to the Pol ε gene suggests that factors in addition to proofreading and MMR influence leading-strand DNA replication fidelity.

Exonuclease 1 preferentially repairs mismatches generated by DNA polymerase α

The Saccharomyces cerevisiae EXO1 gene encodes a 5' exonuclease that participates in mismatch repair (MMR) of DNA replication errors. Deleting EXO1 was previously shown to increase mutation rates to a greater extent when combined with a mutator variant (pol3-L612M) of the lagging strand replicase, DNA polymerase δ (Pol δ), than when combined with a mutator variant (pol2-M644G) of the leading strand replicase, DNA polymerase ɛ (Pol ɛ). Here we confirm that result, and extend the approach to examine the effect of deleting EXO1 in a mutator variant (pol1-L868M) of Pol α, the proofreading-deficient and least accurate of the three nuclear replicases that is responsible for initiating Okazaki fragment synthesis. We find that deleting EXO1 increases the mutation rate in the Pol α mutator strain to a significantly greater extent than in the Pol δ or Pol ɛ mutator strains, thereby preferentially reducing the efficiency of MMR of replication errors generated by Pol α. Because these mismatches are closer to the 5' ends of Okazaki fragments than are mismatches made by Pol δ or Pol ɛ, the results not only support the previous suggestion that Exo1 preferentially excises lagging strand replication errors during mismatch repair, they further imply that the 5' ends serve as entry points for 5' excision of replication errors made by Pol α, and possibly as strand discrimination signals for MMR. Nonetheless, mutation rates in the Pol α mutator strain are 5- to 25-fold lower in an exo1Δ strain as compared to an msh2Δ strain completely lacking MMR, indicating that in the absence of Exo1, most replication errors made by Pol α can still be removed in an Msh2-dependent manner by other nucleases and/or by strand displacement.

DNA polymerase delta in DNA replication and genome maintenance

The eukaryotic genome is in a constant state of modification and repair. Faithful transmission of the genomic information from parent to daughter cells depends upon an extensive system of surveillance, signaling, and DNA repair, as well as accurate synthesis of DNA during replication. Often, replicative synthesis occurs over regions of DNA that have not yet been repaired, presenting further challenges to genomic stability. DNA polymerase δ (pol δ) occupies a central role in all of these processes: catalyzing the accurate replication of a majority of the genome, participating in several DNA repair synthetic pathways, and contributing structurally to the accurate bypass of problematic lesions during translesion synthesis. The concerted actions of pol δ on the lagging strand, pol ϵ on the leading strand, associated replicative factors, and the mismatch repair (MMR) proteins results in a mutation rate of less than one misincorporation per genome per replication cycle. This low mutation rate provides a high level of protection against genetic defects during development and may prevent the initiation of malignancies in somatic cells. This review explores the role of pol δ in replication fidelity and genome maintenance.

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

Multiple sequence changes that are simultaneously introduced in a single DNA transaction have a higher probability of altering gene function than do single base substitutions. DNA polymerase zeta (Pol ζ) has been shown to introduce such clustered mutations under specific selective and/or DNA damage-producing conditions. In this study, a forward mutation assay was used to determine the specificity of spontaneous mutations generated in Saccharomyces cerevisiae when either wild-type Pol ζ or a mutator Pol ζ variant (rev3-L979F) bypasses endogenous lesions. Mutagenesis in strains proficient for nucleotide excision repair (NER) was compared to mutagenesis in NER-deficient strains that retain unrepaired endogenous DNA lesions in the genome. Compared to NER-proficient strains, NER-deficient rad14Δ strains have elevated mutation rates that depend on Pol ζ. Rates are most strongly elevated for tandem base pair substitutions and clusters of multiple, closely spaced mutations. Both types of mutations depend on Pol ζ, but not on Pol η. Rates of each are further elevated in yeast strains bearing the rev3-979F allele. The results indicate that when Pol ζ performs mutagenic bypass of endogenous, helix-distorting lesions, it catalyzes a short track of processive, error-prone synthesis. We discuss the implications of this unique catalytic property of Pol ζ to its evolutionary conservation and possibly to multistage carcinogenesis.

Modulation of mutagenesis in eukaryotes by DNA replication fork dynamics and quality of nucleotide pools

The rate of mutations in eukaryotes depends on a plethora of factors and is not immediately derived from the fidelity of DNA polymerases (Pols). Replication of chromosomes containing the anti-parallel strands of duplex DNA occurs through the copying of leading and lagging strand templates by a trio of Pols α, δ and ϵ, with the assistance of Pol ζ and Y-family Pols at difficult DNA template structures or sites of DNA damage. The parameters of the synthesis at a given location are dictated by the quality and quantity of nucleotides in the pools, replication fork architecture, transcription status, regulation of Pol switches, and structure of chromatin. The result of these transactions is a subject of survey and editing by DNA repair.

Asymmetry indices for analysis and prediction of replication origins in eukaryotic genomes

DNA replication was recently shown to induce the formation of compositional skews in the genomes of the yeasts Saccharomyces cerevisiae and Kluyveromyces lactis. In this work, I have characterized further GC and TA skew variations in the vicinity of S. cerevisiae replication origins and termination sites, and defined asymmetry indices for origin analysis and prediction. The presence of skew jumps at some termination sites in the S. cerevisiae genome was established. The majority of S. cerevisiae replication origins are marked by an oriented consensus sequence called ACS, but no evidence could be found for asymmetric origin firing that would be linked to ACS orientation. Asymmetry indices related to GC and TA skews were defined, and a global asymmetry index I_GC,TA was described. I_GC,TA was found to strongly correlate with origin efficiency in S. cerevisiae and to allow the determination of sets of intergenes significantly enriched in origin loci. The generalized use of asymmetry indices for origin prediction in naive genomes implies the determination of the direction of the skews, i.e. the identification of which strand, leading or lagging, is enriched in G and which one is enriched in T. Recent work indicates that in Candida albicans and in several related species, centromeres contain early and efficient replication origins. It has been proposed that the skew jumps observed at these positions would reflect the activity of these origins, thus allowing to determine the direction of the skews in these genomes. However, I show here that the skew jumps at C. albicans centromeres are not related to replication and that replication-associated GC and TA skews in C. albicans have in fact the opposite directions of what was proposed.

Defect of Dpb2p, a noncatalytic subunit of DNA polymerase ɛ, promotes error prone replication of undamaged chromosomal DNA in Saccharomyces cerevisiae

The Saccharomyces cerevisiae DNA polymerase epsilon holoenzyme (Pol ɛ HE) is composed of four subunits: Pol2p, Dpb2p, Dpb3p and Dpb4p. The biological functions of Pol2p, the catalytic subunit of Pol ɛ, are subject of active investigation, while the role of the other three, noncatalytic subunits, is not well defined. We showed previously that mutations in Dpb2p, a noncatalytic but essential subunit of Pol ɛ HE, influence the fidelity of DNA replication in yeast cells. The strength of the mutator phenotype due to the different dpb2 alleles was inversely proportional to the strength of protein-protein interactions between Pol2p and the mutated forms of Dpb2p. To understand better the mechanisms of the contribution of Dpb2p to the controlling of the level of spontaneous mutagenesis we undertook here a further genetic analysis of the mutator phenotype observed in dpb2 mutants. We demonstrate that the presence of mutated forms of Dpb2p in the cell not only influences the intrinsic fidelity of Pol ɛ but also facilitates more frequent participation of error-prone DNA polymerase zeta (Pol ζ) in DNA replication. The obtained results suggest that the structural integrity of Pol ɛ HE is a crucial contributor to accurate chromosomal DNA replication and, when compromised, favors participation of error prone DNA Pol ζ in this process.

Proofreading of ribonucleotides inserted into DNA by yeast DNA polymerase ɛ

We have investigated the ability of the 3' exonuclease activity of Saccharomyces cerevisiae DNA polymerase ɛ (Pol ɛ) to proofread newly inserted ribonucleotides (rNMPs). During DNA synthesis in vitro, Pol ɛ proofreads ribonucleotides with apparent efficiencies that vary from none at some locations to more than 90% at others, with rA and rU being more efficiently proofread than rC and rG. Previous studies show that failure to repair ribonucleotides in the genome of rnh201Δ strains that lack RNase H2 activity elevates the rate of short deletions in tandem repeat sequences. Here we show that this rate is increased by 2-4-fold in pol2-4 rnh201Δ strains that are also defective in Pol ɛ proofreading. In comparison, defective proofreading in these same strains increases the rate of base substitutions by more than 100-fold. Collectively, the results indicate that although proofreading of an 'incorrect' sugar is less efficient than is proofreading of an incorrect base, Pol ɛ does proofread newly inserted rNMPs to enhance genome stability.

DNA replication fidelity in Escherichia coli: a multi-DNA polymerase affair

High accuracy (fidelity) of DNA replication is important for cells to preserve the genetic identity and to prevent the accumulation of deleterious mutations. The error rate during DNA replication is as low as 10(-9) to 10(-11) errors per base pair. How this low level is achieved is an issue of major interest. This review is concerned with the mechanisms underlying the fidelity of the chromosomal replication in the model system Escherichia coli by DNA polymerase III holoenzyme, with further emphasis on participation of the other, accessory DNA polymerases, of which E. coli contains four (Pols I, II, IV, and V). Detailed genetic analysis of mutation rates revealed that (1) Pol II has an important role as a back-up proofreader for Pol III, (2) Pols IV and V do not normally contribute significantly to replication fidelity, but can readily do so under conditions of elevated expression, (3) participation of Pols IV and V, in contrast to that of Pol II, is specific to the lagging strand, and (4) Pol I also makes a lagging-strand-specific fidelity contribution, limited, however, to the faithful filling of the Okazaki fragment gaps. The fidelity role of the Pol III τ subunit is also reviewed.

Structure and function of eukaryotic DNA polymerase δ

DNA polymerase δ (Pol δ) is a member of the B-family DNA polymerases and is one of the major replicative DNA polymerases in eukaryotes. In addition to chromosomal DNA replication it is also involved in DNA repair and recombination. Pol δ is a multi-subunit complex comprised of a catalytic subunit and accessory subunits. The latter subunits play a critical role in the regulation of Pol δ functions. Recent progress in the structural characterization of Pol δ, together with a vast number of biochemical and functional studies, provides the basis for understanding the intriguing mechanisms of its regulation during DNA replication, repair and recombination. In this chapter we review the current state of the Pol δ structure-function relationship with an emphasis on the role of its accessory subunits.

Mutator suppression and escape from replication error-induced extinction in yeast

Cells rely on a network of conserved pathways to govern DNA replication fidelity. Loss of polymerase proofreading or mismatch repair elevates spontaneous mutation and facilitates cellular adaptation. However, double mutants are inviable, suggesting that extreme mutation rates exceed an error threshold. Here we combine alleles that affect DNA polymerase δ (Pol δ) proofreading and mismatch repair to define the maximal error rate in haploid yeast and to characterize genetic suppressors of mutator phenotypes. We show that populations tolerate mutation rates 1,000-fold above wild-type levels but collapse when the rate exceeds 10⁻³ inactivating mutations per gene per cell division. Variants that escape this error-induced extinction (eex) rapidly emerge from mutator clones. One-third of the escape mutants result from second-site changes in Pol δ that suppress the proofreading-deficient phenotype, while two-thirds are extragenic. The structural locations of the Pol δ changes suggest multiple antimutator mechanisms. Our studies reveal the transient nature of eukaryotic mutators and show that mutator phenotypes are readily suppressed by genetic adaptation. This has implications for the role of mutator phenotypes in cancer.

Organisms strike a balance between genetic continuity and change. Most cells are well adapted to their niches and therefore invest heavily in mechanisms that maintain accurate DNA replication. When cell populations are confronted with changing environmental conditions, “mutator” clones with high mutation rates emerge and readily adapt to the new conditions by rapidly acquiring beneficial mutations. However, deleterious mutations also accumulate, raising the question: what level of mutational burden can cell populations sustain before collapsing? Here we experimentally determine the maximal mutation rate in haploid yeast. We observe that yeast can withstand a 1,000-fold increase in mutation rate without losing colony forming capacity. Yet no strains survive a 10,000-fold increase in mutation rate. Escape mutants with an “anti-mutator” phenotype frequently emerge from cell populations undergoing this error-induced extinction. The diversity of antimutator changes suggests that strong mutator phenotypes in nature may be inherently transient, ensuring that rapid adaptation is followed by genetic attenuation which preserves the beneficial, adaptive mutations. These observations are relevant to microbial populations during infection as well as the somatic evolution of cancer cells.

DNA replication induces compositional biases in yeast

Asymmetries intrinsic to the process of DNA replication are expected to cause differences in the substitution patterns of the leading and the lagging strands and to induce compositional biases. These biases have been detected in the majority of eubacterial genomes but rarely in eukaryotes. Only in the human genome, the activity of a minority of replication origins seems to generate compositional biases. In this work, we provide evidence for replication-associated GC and TA skews in the genomes of two yeast species, Saccharomyces cerevisiae and Kluyveromyces lactis, whereas the data for the Schizosaccharomyces pombe genome are less conclusive. In contrast with the genomes of Homo sapiens and of the majority of eubacteria, the leading strand is enriched in cytosine and adenine in both S. cerevisiae and K. lactis. We observed significant variations across the interorigin intervals of several substitution rates in the S. cerevisiae lineage since its divergence from S. paradoxus. We also found that the S. cerevisiae genome is far from compositional equilibrium and that its present compositional biases are due to substitution rates operating before its divergence from S. paradoxus. Finally, we observed that replication and transcription tend to be cooriented in the S. cerevisiae genome, especially for genes encoding subunits of protein complexes. Taken together, our results suggest that replication-related compositional biases may be a feature of many eukaryotic genomes despite the stochastic nature of the firing of replication origins in these genomes.

G-quadruplex-induced instability during leading-strand replication

G-quadruplexes are four-stranded nucleic acid structures whose biological functions remain poorly understood. In the yeast S. cerevisiae, we report that G-quadruplexes form and, if not properly processed, pose a specific challenge to replication. We show that the G-quadruplex-prone CEB1 tandem array is tolerated when inserted near ARS305 replication origin in wild-type cells but is very frequently destabilized upon treatment with the potent Phen-DC(3) G-quadruplex ligand, or in the absence of the G-quadruplex-unwinding Pif1 helicase, only when the G-rich strand is the template of leading-strand replication. The orientation-dependent instability is associated with the formation of Rad51-Rad52-dependent X-shaped intermediates during replication detected by two-dimensional (2D) gels, and relies on the presence of intact G-quadruplex motifs in CEB1 and on the activity of ARS305. The asymmetrical behaviour of G-quadruplex prone sequences during replication has implications for their evolutionary dynamics within genomes, including the maintenance of G-rich telomeres.

Inaccurate DNA synthesis in cell extracts of yeast producing active human DNA polymerase iota

Mammalian Pol ι has an unusual combination of properties: it is stimulated by Mn²⁺ ions, can bypass some DNA lesions and misincorporates “G” opposite template “T” more frequently than incorporates the correct “A.” We recently proposed a method of detection of Pol ι activity in animal cell extracts, based on primer extension opposite the template T with a high concentration of only two nucleotides, dGTP and dATP (incorporation of “G” versus “A” method of Gening, abbreviated as “misGvA”). We provide unambiguous proof of the “misGvA” approach concept and extend the applicability of the method for the studies of variants of Pol ι in the yeast model system with different cation cofactors. We produced human Pol ι in baker's yeast, which do not have a POLI ortholog. The “misGvA” activity is absent in cell extracts containing an empty vector, or producing catalytically dead Pol ι, or Pol ι lacking exon 2, but is robust in the strain producing wild-type Pol ι or its catalytic core, or protein with the active center L62I mutant. The signature pattern of primer extension products resulting from inaccurate DNA synthesis by extracts of cells producing either Pol ι or human Pol η is different. The DNA sequence of the template is critical for the detection of the infidelity of DNA synthesis attributed to DNA Pol ι. The primer/template and composition of the exogenous DNA precursor pool can be adapted to monitor replication fidelity in cell extracts expressing various error-prone Pols or mutator variants of accurate Pols. Finally, we demonstrate that the mutation rates in yeast strains producing human DNA Pols ι and η are not elevated over the control strain, despite highly inaccurate DNA synthesis by their extracts.

The role of DNA polymerase alpha in the control of mutagenesis in Saccharomyces cerevisiae cells starved for nutrients

In nature, micro organisms experience numerous environmental stresses and generally grow poorly most of the time. In the last two decades it has become evident that mutations arise not only in actively dividing cells but also in non-replicating or slowly replicating cells starved for nutrients. In yeast, precise base selection and proofreading by replicative DNA polymerases δ and ε keep starvation-associated mutagenesis (SAM) at basal levels. Less is known about the role of replicative DNA polymerase α (Pol α). Here we provide evidence that Pol α is involved in the control of SAM in yeast cells starved for adenine by participation in sporadic replication and/or DNA repair under these conditions.

Differential correction of lagging-strand replication errors made by DNA polymerases {alpha} and {delta}

Mismatch repair (MMR) of replication errors requires DNA ends that can direct repair to the newly synthesized strand containing the error. For all but those organisms that use adenine methylation to generate nicks, the source of these ends in vivo is unknown. One possibility is that MMR may have a "special relation to the replication complex" [Wagner R, Jr., Meselson M (1976) Proc Natl Acad Sci USA 73:4135-4139], perhaps one that allows 5' or 3' DNA ends associated with replication to act as strand discrimination signals. Here we examine this hypothesis, based on the logic that errors made by yeast DNA polymerase α (Pol α), which initiates Okazaki fragments during lagging-strand replication, will always be closer to a 5' end than will be more internal errors generated by DNA polymerase δ (Pol δ), which takes over for Pol α to complete lagging-strand replication. When we compared MMR efficiency for errors made by variant forms of these two polymerases, Msh2-dependent repair efficiencies for mismatches made by Pol α were consistently higher than for those same mismatches when made by Pol δ. Thus, one special relationship between MMR and replication is that MMR is more efficient for the least accurate of the major replicative polymerases, exonuclease-deficient Pol α. This observation is consistent with the close proximity and possible use of 5' ends of Okazaki fragments for strand discrimination, which could increase the probability of Msh2-dependent MMR by 5' excision, by a Msh2-dependent strand displacement mechanism, or both.

DNA replication fidelity and cancer

Cancer is fueled by mutations and driven by adaptive selection. Normal cells avoid deleterious mutations by replicating their genomes with extraordinary accuracy. Here we review the pathways governing DNA replication fidelity and discuss evidence implicating replication errors (point mutation instability or PIN) in carcinogenesis.

Epigenetically-inherited centromere and neocentromere DNA replicates earliest in S-phase

Eukaryotic centromeres are maintained at specific chromosomal sites over many generations. In the budding yeast Saccharomyces cerevisiae, centromeres are genetic elements defined by a DNA sequence that is both necessary and sufficient for function; whereas, in most other eukaryotes, centromeres are maintained by poorly characterized epigenetic mechanisms in which DNA has a less definitive role. Here we use the pathogenic yeast Candida albicans as a model organism to study the DNA replication properties of centromeric DNA. By determining the genome-wide replication timing program of the C. albicans genome, we discovered that each centromere is associated with a replication origin that is the first to fire on its respective chromosome. Importantly, epigenetic formation of new ectopic centromeres (neocentromeres) was accompanied by shifts in replication timing, such that a neocentromere became the first to replicate and became associated with origin recognition complex (ORC) components. Furthermore, changing the level of the centromere-specific histone H3 isoform led to a concomitant change in levels of ORC association with centromere regions, further supporting the idea that centromere proteins determine origin activity. Finally, analysis of centromere-associated DNA revealed a replication-dependent sequence pattern characteristic of constitutively active replication origins. This strand-biased pattern is conserved, together with centromere position, among related strains and species, in a manner independent of primary DNA sequence. Thus, inheritance of centromere position is correlated with a constitutively active origin of replication that fires at a distinct early time. We suggest a model in which the distinct timing of DNA replication serves as an epigenetic mechanism for the inheritance of centromere position.

Centromeres form at the same chromosomal position from generation to generation, yet in most species this inheritance occurs in a DNA sequence–independent manner that is not well understood. Here, we determine the timing of DNA replication across the genome of the human fungal pathogen Candida albicans and find that centromeric DNA is the first locus to replicate on each chromosome. Furthemore, this unique replication timing may be important for centromere inheritance, based on several observations. First, DNA sequence patterns at centromeres indicate that, despite high levels of primary sequence divergence, the region has served as a replication origin for millions of years; second, formation of a neocentromere (a new centromere formed at an ectopic locus following deletion of the native centromere DNA) results in the establishment of a new, early-firing origin of replication; and third, a centromere-specific protein, Cse4p, recruits origin replication complex proteins in a concentration-dependent manner. Thus, centromere position is inherited by an epigenetic mechanism that appears to be defined by a distinctively early firing DNA replication origin.

In vivo mutation analysis using the ΦX174 transgenic mouse and comparisons with other transgenes and endogenous genes

The ΦX174 transgenic mouse was first developed as an in vivo Ames test, detecting base pair substitution (bps) at a single bp in a reversion assay. A forward mutational assay was also developed, which is a gain of function assay that also detects bps exclusively. Later work with both assays focused on establishing that a mutation was fixed in vivo using single-burst analysis: determining the number of mutant progeny virus from an electroporated cell by dividing the culture into aliquots before scoring mutants. We review results obtained from single-burst analysis, including testing the hypothesis that high mutant frequencies (MFs) of G:C to A:T mutation recovered by transgenic targets include significant numbers of unrepaired G:T mismatches. Comparison between the ΦX174 and lacI transgenes in mouse spleen indicates that the spontaneous bps mutation frequency per nucleotide (mf(n)) is not significantly lower for ΦX174 than for lacI; the response to ENU is also comparable. For the lacI transgene, the spontaneous bps mf(n) is highly age-dependent up to 12 weeks of age and the linear trend extrapolates at conception to a frequency close to the human bps mf(n) per generation of 1.7 × 10(-8). Unexpectedly, we found that the lacI somatic (spleen) bps mf(n) per cell division at early ages was estimated to be the same as for the human germ-line. The bps mf(n) in bone marrow for the gpt transgene is comparable to spleen for the lacI and ΦX174 transgenes. We conclude that the G:C to A:T transition is characteristic of spontaneous in vivo mutation and that the MFs measured in these transgenes at early ages reflect the expected accumulation of in vivo mutation typical of endogenous mammalian mutation rates. However, spontaneous and induced mf(n)s per nucleotide for the cII gene in spleen are 5-10 times higher than for these other transgenes.

Participation of DNA polymerase zeta in replication of undamaged DNA in Saccharomyces cerevisiae

Translesion synthesis DNA polymerases contribute to DNA damage tolerance by mediating replication of damaged templates. Due to the low fidelity of these enzymes, lesion bypass is often mutagenic. We have previously shown that, in Saccharomyces cerevisiae, the contribution of the error-prone DNA polymerase zeta (Polzeta) to replication and mutagenesis is greatly enhanced if the normal replisome is defective due to mutations in replication genes. Here we present evidence that this defective-replisome-induced mutagenesis (DRIM) results from the participation of Polzeta in the copying of undamaged DNA rather than from mutagenic lesion bypass. First, DRIM is not elevated in strains that have a high level of endogenous DNA lesions due to defects in nucleotide excision repair or base excision repair pathways. Second, DRIM remains unchanged when the level of endogenous oxidative DNA damage is decreased by using anaerobic growth conditions. Third, analysis of the spectrum of mutations occurring during DRIM reveals the characteristic error signature seen during replication of undamaged DNA by Polzeta in vitro. These results extend earlier findings in Escherichia coli indicating that Y-family DNA polymerases can contribute to the copying of undamaged DNA. We also show that exposure of wild-type yeast cells to the replication inhibitor hydroxyurea causes a Polzeta-dependent increase in mutagenesis. This suggests that DRIM represents a response to replication impediment per se rather than to specific defects in the replisome components.

DNA polymerases at the eukaryotic fork-20 years later

Function of the eukaryotic genome depends on efficient and accurate replication of anti-parallel DNA strands. Eukaryotic DNA polymerases have different properties adapted to perform a wide spectrum of DNA transactions. Here we focus on major players in the bulk replication, DNA polymerases of the B-family. We review the organization of the replication fork in eukaryotes in a historical perspective, analyze contemporary models and propose a new integrative model of the fork.

Single-molecule studies of DNA replisome function

Fast and accurate replication of DNA is accomplished by the interactions of multiple proteins in the dynamic DNA replisome. The DNA replisome effectively coordinates the leading and lagging strand synthesis of DNA. These complex, yet elegantly organized, molecular machines have been studied extensively by kinetic and structural methods to provide an in-depth understanding of the mechanism of DNA replication. Owing to averaging of observables, unique dynamic information of the biochemical pathways and reactions is concealed in conventional ensemble methods. However, recent advances in the rapidly expanding field of single-molecule analyses to study single biomolecules offer opportunities to probe and understand the dynamic processes involved in large biomolecular complexes such as replisomes. This review will focus on the recent developments in the biochemistry and biophysics of DNA replication employing single-molecule techniques and the insights provided by these methods towards a better understanding of the intricate mechanisms of DNA replication.