Spontaneous mutations in diploid Saccharomyces cerevisiae: more beneficial than expected

We performed a 1012-generation mutation-accumulation (MA) experiment in the yeast, Saccharomyces cerevisiae. The MA lines exhibited a significant reduction in mean fitness and a significant increase in variance in fitness. We found that 5.75% of the fitness-altering mutations accumulated were beneficial. This finding contradicts the widely held belief that nearly all fitness-altering mutations are deleterious. The mutation rate was estimated as 6.3 x 10(-5) mutations per haploid genome per generation and the average heterozygous fitness effect of a mutation as 0.061. These estimates are compatible with previous estimates in yeast.     

Mutator phenotypes conferred by MLH1 overexpression and by heterozygosity for mlh1 mutations

Loss of DNA mismatch repair due to mutation or diminished expression of the MLH1 gene is associated with genome instability and cancer. In this study, we used a yeast model system to examine three circumstances relevant to modulation of MLH1 function. First, overexpression of wild-type MLH1 was found to cause a strong elevation of mutation rates at three different loci, similar to the mutator effect of MLH1 gene inactivation. Second, haploid yeast strains with any of six mlh1 missense mutations that mimic germ line mutations found in human cancer patients displayed a strong mutator phenotype consistent with loss of mismatch repair function. Five of these mutations affect amino acids that are homologous to residues suggested by recent crystal structure and biochemical analysis of Escherichia coli MutL to participate in ATP binding and hydrolysis. Finally, using a highly sensitive reporter gene, we detected a mutator phenotype of diploid yeast strains that are heterozygous for mlh1 mutations. Evidence suggesting that this mutator effect results not from reduced mismatch repair in the MLH1/mlh1 cells but rather from loss of the wild-type MLH1 allele in a fraction of cells is presented. Exposure to bleomycin or to UV irradiation strongly enhanced mutagenesis in the heterozygous strain but had little effect on the mutation rate in the wild-type strain. This damage-induced hypermutability may be relevant to cancer in humans with germ line mutations in only one MLH1 allele.     

Direct estimate of the mutation rate and the distribution of fitness effects in the yeast Saccharomyces cerevisiae

Estimates of the rate and frequency distribution of deleterious effects were obtained for the first time by direct scoring and characterization of individual mutations. This was achieved by applying tetrad analysis to a large number of yeast clones. The genomic rate of spontaneous mutation deleterious to a basic fitness-related trait, that of growth rate, was U = 1.1 x 10(-3) per diploid cell division. Extrapolated to the fruit fly and humans, the per generation rate would be 0.074 and 0.92, respectively. This is likely to be an underestimate because single mutations with selection coefficients s < 0.01 could not be detected. The distribution of s > or = 0.01 was studied both for spontaneous and induced mutations. The latter were induced by ethyl methanesulfonate (EMS) or resulted from defective mismatch repair. Lethal changes accounted for approximately 30-40% of the scored mutations. The mean s of nonlethal mutations was fairly high, but most frequently its value was between 0.01 and 0.05. Although the rate and distribution of very small effects could not be determined, the joint share of such mutations in decreasing average fitness was probably no larger than approximately 1%.     

Fitness evolution and the rise of mutator alleles in experimental Escherichia coli populations

We studied the evolution of high mutation rates and the evolution of fitness in three experimental populations of Escherichia coli adapting to a glucose-limited environment. We identified the mutations responsible for the high mutation rates and show that their rate of substitution in all three populations was too rapid to be accounted for simply by genetic drift. In two of the populations, large gains in fitness relative to the ancestor occurred as the mutator alleles rose to fixation, strongly supporting the conclusion that mutator alleles fixed by hitchhiking with beneficial mutations at other loci. In one population, no significant gain in fitness relative to the ancestor occurred in the population as a whole while the mutator allele rose to fixation, but a substantial and significant gain in fitness occurred in the mutator subpopulation as the mutator neared fixation. The spread of the mutator allele from rarity to fixation took >1000 generations in each population. We show that simultaneous adaptive gains in both the mutator and wild-type subpopulations (clonal interference) retarded the mutator fixation in at least one of the populations. We found little evidence that the evolution of high mutation rates accelerated adaptation in these populations.     

Contrasting dynamics of a mutator allele in asexual populations of differing size

Mutators have been shown to hitchhike in asexual populations when the anticipated beneficial mutation supply rate of the mutator subpopulation, NU(b) (for subpopulation of size N and beneficial mutation rate U(b)) exceeds that of the wild-type subpopulation. Here, we examine the effect of total population size on mutator dynamics in asexual experimental populations of Saccharomyces cerevisiae. Although mutators quickly hitchhike to fixation in smaller populations, mutator fixation requires more and more time as population size increases; this observed delay in mutator hitchhiking is consistent with the expected effect of clonal interference. Interestingly, despite their higher beneficial mutation supply rate, mutators are supplanted by the wild type in very large populations. We postulate that this striking reversal in mutator dynamics is caused by an interaction between clonal interference, the fitness cost of the mutator allele, and infrequent large-effect beneficial mutations in our experimental populations. Our work thus identifies a potential set of circumstances under which mutator hitchhiking can be inhibited in natural asexual populations, despite recent theoretical predictions that such populations should have a net tendency to evolve ever-higher genomic mutation rates.     

King''s Formula for the Mutation Load with Epistasis

A formula by J. L. King gives the equilibrium mutation load as L = 2 sigma ui(1 - qi)/z - x) in which ui is the mutation rate to deleterious alleles at the ith locus, qi is the frequency of mutant alleles at this locus, x is the mean number of such mutant genes per individual before selection, z is the mean number in individuals eliminated by selection, and the summation is over all relevant loci. We show that this rule is inaccurate for intense selection and that a correct formula is L = 2 sigma ui(1 - qi) w/(z - x) = 2U w/(z - x) = 2U/(z - x + 2U) in which U is the mean number of new mutations per haploid genome in the population and w is the mean relative fitness before selection. If w/(z - x) less than 1/2, the mutation load is less than the Haldane value (U less than or equal to L less than or equal to 2U) and can be considerably less. In a diploid asexual population, however, with independent occurrence of mutations, L = 1 - e-2U regardless of the mode of selection.     

Mutational specificity of the dnaE173 mutator associated with a defect in the catalytic subunit of DNA polymerase III of Escherichia coli.

We developed a system to examine forward mutations that occurred in the rpsL gene of Escherichia coli placed on a multicopy plasmid. Using this system we determined the mutational specificity for a dnaE173 mutator strain in which the editing function of DNA polymerase III is impeded. The frequency of rpsL- mutations increased 32,000-fold, due to the dnaE173 mutator, and 87 independent rpsL- mutations in the mutator strain were analyzed by DNA sequencing, together with 100 mutants recovered from dnaE+ strain, as the control. While half the number of mutations that occurred in the wild-type strain were caused by insertion elements, no such mutations were recovered from the mutator strain. A novel class of mutation, named "sequence substitution" was present in mutants raised in the dnaE173 strain; seven sequence substitutions induced in the mutator strain occurred at six sites, and all were located in quasipalindromic sequences, carrying the GTG or CAC sequence at one or both endpoints. While other types of mutation were found in both strains, single-base frameshifts were the most frequent events in the mutator strain. Thus, the mutator effect on this class of mutation was 175,000-fold. A total of 95% of the single-base frameshifts in the mutator strain were additions, most of which occurred at runs of A or C bases so as to increase the number of identical residues. Base substitutions, the frequency of which was enhanced 25,000-fold by the mutator effect, occurred primarily at several hotspots in the mutator strain, whereas those induced in the wild-type strain were more randomly distributed throughout the rpsL sequence. The dnaE173 mutator also increased the frequency of duplications 28,000-fold. Of the three duplications recovered from the mutator strain, one was a simple duplication, the region of which was flanked by direct repeats. The other duplications were complex, one half part of which was in the inverted orientation of a region containing two sets of inverted repeats. The same duplications were also recovered from the wild-type strain. The present data suggest that dnaE173 is a novel class of mutator that sharply induces sequence-directed mutagenesis, yielding high frequencies of single base frameshifts, duplications with inversions, sequence substitutions and base substitutions at hotspots.     

Spontaneously arising mutL mutators in evolving Escherichia coli populations are the result of changes in repeat length

Over the course of thousands of generations of growth in a glucose-limited environment, 3 of 12 experimental populations of Escherichia coli spontaneously and independently evolved greatly increased mutation rates. In two of the populations, the mutations responsible for this increased mutation rate lie in the same region of the mismatch repair gene mutL. In this region, a 6-bp repeat is present in three copies in the gene of the wild-type ancestor of the experimental populations but is present in four copies in one of the experimental populations and two copies in the other. These in-frame mutations either add or delete the amino acid sequence LA in the MutL protein. We determined that the replacement of the wild-type sequence with either of these mutations was sufficient to increase the mutation rate of the wild-type strain to a level comparable to that of the mutator strains. Complementation of strains bearing the mutator mutations with wild-type copies of either mutL or the mismatch repair gene uvrD rescued the wild-type mutation rate. The position of the mutator mutations-in the region of MutL known as the ATP lid-suggests a possible deficiency in MutL's ATPase activity as the cause of the mutator phenotype. The similarity of the two mutator mutations (despite the independent evolutionary histories of the populations that gave rise to them) leads to a discussion of the potential adaptive role of DNA repeats.     

Mutational Spectrum Drives the Rise of Mutator Bacteria

Understanding how mutator strains emerge in bacterial populations is relevant both to evolutionary theory and to reduce the threat they pose in clinical settings. The rise of mutator alleles is understood as a result of their hitchhiking with linked beneficial mutations, although the factors that govern this process remain unclear. A prominent but underappreciated fact is that each mutator allele increases only a specific spectrum of mutational changes. This spectrum has been speculated to alter the distribution of fitness effects of beneficial mutations, potentially affecting hitchhiking. To study this possibility, we analyzed the fitness distribution of beneficial mutations generated from different mutator and wild-type Escherichia coli strains. Using antibiotic resistance as a model system, we show that mutational spectra can alter these distributions substantially, ultimately determining the competitive ability of each strain across environments. Computer simulation showed that the effect of mutational spectrum on hitchhiking dynamics follows a non-linear function, implying that even slight spectrum-dependent fitness differences are sufficient to alter mutator success frequency by several orders of magnitude. These results indicate an unanticipated central role for the mutational spectrum in the evolution of bacterial mutation rates. At a practical level, this study indicates that knowledge of the molecular details of resistance determinants is crucial for minimizing mutator evolution during antibiotic therapy.     

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

A Mutator Affecting the Region of the Iso-1-Cytochrome c Gene in Yeast

The mutator gene DEL1 in the yeast Saccharomyces cerevisiae causes a high rate of formation of multisite mutations that encompass the following three adjacent genes: CYC1, which determines the structure of iso-1-cytochrome c; RAD7, which controls UV sensitivity; and OSM1, which controls osomotic sensitivity. The simplest hypothesis is that these multisite mutations are deletions, although it has not been excluded that they may involve other types of gross chromosomal aberrations. In contrast, normal strains do not produce such multisite mutations even after mutagenic treatments. The multisite mutations arise at a rate of approximately 10(-5) to 10(-6) per cell per division in DEL1 strains, which is much higher than rates observed for mutation of genes in normal strains. For example, normal strains produce all types of cyc1 mutants at a low rate of approximately 10(-8) to 10(-9). No evidence for multisite mutations was obtained upon analysis of numerous spontaneous ade1, ade2, met2 and met15 mutants isolated in a DEL1 strain. DEL1 appears to be both cis- and trans-dominant. The location of the DEL1 gene and the lack of effect on other genes suggest that the mutator acts only on a region adjacent to itself.     

Elimination of the Yeast Rad6 Ubiquitin Conjugase Enhances Base-Pair Transitions and G.c -> T.a Transversions as Well as Transposition of the Ty Element: Implications for the Control of Spontaneous Mutation

The RAD6 gene of the yeast Saccharomyces cerevisiae encodes an enzyme that conjugates ubiquitin to other proteins. Defects in RAD6 confer a mutator phenotype due, in part, to an increased rate of transposition of the yeast Ty element. To further delineate the role of protein ubiquitination in the control of spontaneous mutagenesis in yeast, we have characterized 202 mutations that arose spontaneously in the SUP4-o gene carried on a centromere vector in a RAD6 deletion strain. The resulting mutational spectrum was compared to that for 354 spontaneous SUP4-o mutations isolated in the isogenic wild-type parent. This comparison revealed that the rad6 mutator enhanced the rate of single base-pair substitution, as well as Ty insertion, but did not affect the rates of the other mutational classes detected. Relative to the wild-type parent, Ty inserted at considerably more SUP4-o positions in the rad6 strain with a significantly smaller fraction detected at a transposition hotspot. These findings suggest that, in addition to the rate of transposition, protein ubiquitination might influence the target site specificity of Ty insertion. The increase in the substitution rate accounted for approximately 90% of the rad6 mutator effect but only the two transitions and the G. C----T.A transversion were enhanced. Analysis of the distribution of these events within SUP4-o suggested that the site specificity of the substitutions was influenced by DNA sequence context. Transformation of heteroduplex plasmid DNAs into the two strains demonstrated that the rad6 mutator did not reduce the efficiency of correcting mismatches that could give rise to the transitions or transversion nor did it bias restoration of the mismatches to the incorrect base-pairs. These results are discussed in relation to possible mechanisms that might link ubiquitination of proteins to spontaneous mutation rates.     

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.     

The approach to mutation-selection balance in an infinite asexual population, and the evolution of mutation rates

A method is described for calculating the dynamics of the distribution of fitness in an infinite asexual population which is subject to unconditionally deleterious mutations with independent effects. This method is applied to the problem of calculating the frequency of a mutator subpopulation, at equilibrium between mutation and indirect selection due to association with deleterious mutations. Many mutator alleles are produced by loss-of-function mutations in polymerase or mismatch repair genes. Previous calculations have ignored the fact that this creates a flux of higher fitness individuals into the mutator subpopulation. This flux raises the mean fitness of the mutator subpopulation, and when this factor is taken into account, the frequency of the mutator may be more than an order of magnitude greater than recent theoretical work has suggested.     

Experimental adaptation of Salmonella typhimurium to mice

Experimental evolution is a powerful approach to study the dynamics and mechanisms of bacterial niche specialization. By serial passage in mice, we evolved 18 independent lineages of Salmonella typhimurium LT2 and examined the rate and extent of adaptation to a mainly reticuloendothelial host environment. Bacterial mutation rates and population sizes were varied by using wild-type and DNA repair-defective mutator (mutS) strains with normal and high mutation rates, respectively, and by varying the number of bacteria intraperitoneally injected into mice. After <200 generations of adaptation all lineages showed an increased fitness as measured by a faster growth rate in mice (selection coefficients 0.11-0.58). Using a generally applicable mathematical model we calculated the adaptive mutation rate for the wild-type bacterium to be >10(-6)/cell/generation, suggesting that the majority of adaptive mutations are not simple point mutations. For the mutator lineages, adaptation to mice was associated with a loss of fitness in secondary environments as seen by a reduced metabolic capability. During adaptation there was no indication that a high mutation rate was counterselected. These data show that S. typhimurium can rapidly and extensively increase its fitness in mice but this niche specialization is, at least in mutators, associated with a cost.     

Estimating numbers of EMS-induced mutations affecting life history traits in Caenorhabditis elegans in crosses between inbred sublines

Inbred lines of the nematode Caenorhabditis elegans containing independent EMS-induced mutations were crossed to the ancestral wild-type strain (N2). Replicated inbred sublines were generated from the F1 offspring under conditions of minimal selection and, along with the N2 and mutant progenitor lines, were assayed for several fitness correlates including relative fitness (w). A modification of the Castle-Wright estimator and a maximum-likelihood (ML) method were used to estimate the numbers and effects of detectable mutations affecting these characters. The ML method allows for variation in mutational effects by fitting either one or two classes of mutational effect, and uses a Box-Cox power transformation of residual values to account for a skewed distribution of residuals. Both the Castle-Wright and the ML analyses suggest that most of the variation among sublines was due to a few (approximately 1.5-2.5 on average) large-effect mutations. Under ML, a model with two classes of mutational effects, including a class with small effects, fitted better than a single mutation class model, although not significantly better. Nonetheless, given that we expect there to be many mutations induced per line, our results support the hypothesis that mutations vary widely in their effects.     

Low impact of germline transposition on the rate of mildly deleterious mutation in Caenorhabditis elegans

Little is known about the role of transposable element (TE) insertion in the production of mutations with mild effects on fitness, the class of mutations thought to be central to the evolution of many basic features of natural populations. We propagated mutation-accumulation (MA) lines of two RNAi-deficient strains of Caenorhabditis elegans that exhibit germline transposition. We show here that the impact of TE activity was to raise the level of mildly deleterious mutation by 2- to 8.5-fold, as estimated from fecundity, longevity, and body length measurements, compared to that observed in a parallel MA experiment with a control strain characterized by a lack of germline transposition. Despite this increase, the rate of mildly deleterious mutation was between one and two orders of magnitude lower than the rate of TE accumulation, which was approximately two new insertions per genome per generation. This study suggests that high rates of TE activity do not necessarily translate into high rates of detectable nonlethal mutation.     

Disruption of the Rad52 Gene Alters the Spectrum of Spontaneous Sup4-O Mutations in Saccharomyces Cerevisiae

Defects in the RAD52 gene of the yeast Saccharomyces cerevisiae confer a mutator phenotype. To characterize this effect in detail, a collection of 238 spontaneous SUP4-o mutations arising in a strain having a disrupted RAD52 gene was analyzed by DNA sequencing. The resulting mutational spectrum was compared to that derived from an examination of 222 spontaneous mutations selected in a nearisogenic wild-type (RAD52) strain. This comparison revealed that the mutator phenotype was associated with an increase in the frequency of base-pair substitutions. All possible types of substitution were detected but there was a reduction in the relative fraction of A.T----G.C transitions and an increase in the proportion of G.C----C.G transversions. These changes were sufficient to cause a twofold greater preference for substitutions at G.C sites in the rad52 strain despite a decrease in the fraction of G.C----T.A transversions. There were also considerable differences between the distributions of substitutions within the SUP4-o gene. Base-pair changes occurred at fewer sites in the rad52 strain but the mutated sites included several that were not detected in the RAD52 background. Only two of the four sites that were mutated most frequently in the rad52 strain were also prominent in the wild-type strain and mutation frequencies at almost all sites common to both strains were greater for the rad52 derivative. Although single base-pair deletions occurred in the two strains with similar frequencies, several classes of mutation that were recovered in the wild-type background including multiple base-pair deletions, insertions of the yeast transposable element Ty, and more complex changes, were not detected in the rad52 strain.(ABSTRACT TRUNCATED AT 250 WORDS)     

Environment dependence of mutational parameters for viability in Drosophila melanogaster

The genomic rate of mildly deleterious mutations (U) figures prominently in much evolutionary and ecological theory. In Drosophila melanogaster, estimates of U have varied widely, from <0.1 to nearly 1 per zygote. The source of this variation is unknown, but could include differences in the conditions used for assaying fitness traits. We examined how assay conditions affect estimates of the rates and effects of viability-depressing mutations in two sets of lines with accumulated spontaneous mutations on the second chromosome. In each set, the among-line variance in egg-to-adult viability was significantly greater when viability was assayed using a high parental density than when it was assayed using a low density. In contrast, the proportional decline in viability due to new mutations did not differ between densities. Two other manipulations, lowering the temperature and adding ethanol to the medium, had no significant effects on either the mean decline or among-line variance. Cross-environment genetic correlations in viability were generally close to one, implying that most mutations reduced viability in all environments. Using data from the low-density, lower-bound estimates of U approached the classic, high values of Mukai and Ohnishi; at the high density, U estimates were similar to recently reported low values. The difference in estimated mutation rates, taken at face value, would imply that many mutations affected fitness at low density but not at high density, but this is shown to be incompatible with the observed high cross-environment correlations. Possible reasons for this discrepancy are discussed. Regardless of the interpretation, the results show that assay conditions can have a large effect on estimates of mutational parameters for fitness traits.     

Rates and fitness consequences of new mutations in humans

The human mutation rate per nucleotide site per generation (μ) can be estimated from data on mutation rates at loci causing Mendelian genetic disease, by comparing putatively neutrally evolving nucleotide sequences between humans and chimpanzees and by comparing the genome sequences of relatives. Direct estimates from genome sequencing of relatives suggest that μ is about 1.1 × 10(-8), which is about twofold lower than estimates based on the human-chimp divergence. This implies that an average of ~70 new mutations arise in the human diploid genome per generation. Most of these mutations are paternal in origin, but the male:female mutation rate ratio is currently uncertain and might vary even among individuals within a population. On the basis of a method proposed by Kondrashov and Crow, the genome-wide deleterious mutation rate (U) can be estimated from the product of the number of nucleotide sites in the genome, μ, and the mean selective constraint per site. Although the presence of many weakly selected mutations in human noncoding DNA makes this approach somewhat problematic, estimates are U ≈ 2.2 for the whole diploid genome per generation and 0.35 for mutations that change an amino acid of a protein-coding gene. A genome-wide deleterious mutation rate of 2.2 seems higher than humans could tolerate if natural selection is "hard," but could be tolerated if selection acts on relative fitness differences between individuals or if there is synergistic epistasis. I argue that in the foreseeable future, an accumulation of new deleterious mutations is unlikely to lead to a detectable decline in fitness of human populations.