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.     

Impact of mutation rate on the adaptation of gut bacteria

To study the role of mutator bacteria in the evolution of bacterial populations, we followed the impact of the mutation rate of Escherichia coli strains in the colonisation of the gut of axenic mice and the evolution of the mutation rate of bacterial populations living in the gut. We show that mutator bacteria have an advantage during the colonization. This adaptive advantage comes from their ability to generate adaptive mutations faster than wild type strains, mutations that allow their maintenance in the ecosystem. However, while mutator bacteria are becoming specialised to the environment they are living in, they accumulate mutations that may be deleterious or lethal in secondary environments. By following the evolution of the mutation rate of bacterial populations living in the gut of mice receiving antibiotics, we show that this therapy selects not only for antibiotic resistant mutants but also for mutator alleles that enhance mutation rates and are responsible for the appearance of the resistance. The costs of a high mutation rate, due to the accumulation of mutations, is seen in environments where changes are recurrent. In an ever-changing situation where every change is new, mutator bacteria might help the evolution of bacterial populations.     

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 trade-off between oxidative stress resistance and DNA repair plays a role in the evolution of elevated mutation rates in bacteria

The dominant paradigm for the evolution of mutator alleles in bacterial populations is that they spread by indirect selection for linked beneficial mutations when bacteria are poorly adapted. In this paper, we challenge the ubiquity of this paradigm by demonstrating that a clinically important stressor, hydrogen peroxide, generates direct selection for an elevated mutation rate in the pathogenic bacterium Pseudomonas aeruginosa as a consequence of a trade-off between the fidelity of DNA repair and hydrogen peroxide resistance. We demonstrate that the biochemical mechanism underlying this trade-off in the case of mutS is the elevated secretion of catalase by the mutator strain. Our results provide, to our knowledge, the first experimental evidence that direct selection can favour mutator alleles in bacterial populations, and pave the way for future studies to understand how mutation and DNA repair are linked to stress responses and how this affects the evolution of bacterial mutation rates.     

Evolution of mutation rates in bacteria

Evolutionary success of bacteria relies on the constant fine-tuning of their mutation rates, which optimizes their adaptability to constantly changing environmental conditions. When adaptation is limited by the mutation supply rate, under some conditions, natural selection favours increased mutation rates by acting on allelic variation of the genetic systems that control fidelity of DNA replication and repair. Mutator alleles are carried to high frequency through hitchhiking with the adaptive mutations they generate. However, when fitness gain no longer counterbalances the fitness loss due to continuous generation of deleterious mutations, natural selection favours reduction of mutation rates. Selection and counter-selection of high mutation rates depends on many factors: the number of mutations required for adaptation, the strength of mutator alleles, bacterial population size, competition with other strains, migration, and spatial and temporal environmental heterogeneity. Such modulations of mutation rates may also play a role in the evolution of antibiotic resistance.     

An overview of the mechanisms of mutagenesis and carcinogenesis

Cancer is a genetic disease due to the accumulation of numerous mutations rendering the tumour cell insensitive to control by the local cellular environment and by the whole organism. Analysis of the frequency of appearance of human cancer as a function of age shows that between four and seven mutations in key genes are usually necessary to produce most human cancers. Interesting debates in the literature are concerned with the idea that normal mutation rates followed by selective advantage of mutated clones are enough to produce the numerous mutations found in human cancers. Alternatively, the mutator phenotype hypothesis is based on the idea that the normal mutation rates are insufficient to account for the multiple mutations found in tumours. It is, however, difficult not only to know this exact mutation frequency in cells but also to know the total number of cell divisions giving rise to a cancer. Therefore, during at least one step in the carcinogenic process, a mutator phenotype in target cells may occur due to mutations controlling the fidelity of DNA replication or DNA repair, the apoptosis pathways or the cell cycle checkpoint regulations. Among the multiple mutations found in human cancers such as gene amplification, chromosome alterations and translocations, point mutations are very important and the molecular mechanisms of their production are well documented. I will describe in detail the various mechanisms that a cell can use to produce point mutations due to lower fidelity in the DNA polymerisation step or to inefficient repair pathways. The presence of multiple mutations in human cancer is interesting not only in terms of understanding the carcinogenesis process in humans but also in eventually promoting strategies to decrease the efficiency of this process and to increase cancer therapy regimen.     

The evolution of low mutation rates in experimental mutator populations of Saccharomyces cerevisiae

Mutation is the source of both beneficial adaptive variation and deleterious genetic load, fueling the opposing selective forces than shape mutation rate evolution. This dichotomy is well illustrated by the evolution of the mutator phenotype, a genome-wide 10- to 100-fold increase in mutation rate. This phenotype has often been observed in clonally expanding populations exposed to novel or frequently changing conditions. Although studies of both experimental and natural populations have shed light on the evolutionary forces that lead to the spread of the mutator allele through a population, significant gaps in our understanding of mutator evolution remain. Here we use an experimental evolution approach to investigate the conditions required for the evolution of a reduction in mutation rate and the mechanisms by which populations tolerate the accumulation of deleterious mutations. We find that after ～6,700 generations, four out of eight experimental mutator lines had evolved a decreased mutation rate. We provide evidence that the accumulation of deleterious mutations leads to selection for reduced mutation rate clones in populations of mutators. Finally, we test the long-term consequences of the mutator phenotype, finding that mutator lines follow different evolutionary trajectories, some of which lead to drug resistance.     

Evolutionary Consequences of Mutation and Selection within an Individual

Whether in sexual or asexual organisms, selection among cell lineages during development is an effective way of eliminating deleterious mutations. Using a mathematical analysis, we find that relatively small differences in cell replication rates during development can translate into large differences in the proportion of mutant cells within the adult, especially when development involves a large number of cell divisions. Consequently, intraorganismal selection can substantially reduce the deleterious mutation rate observed among offspring as well as the mutation load within a population, because cells rather than individuals provide the selective ``deaths' necessary to stem the tide of deleterious mutations. The reduction in mutation rate among offspring is more pronounced in organisms with plastic development than in those with structured development. It is also more pronounced in asexual organisms that produce multicellular rather than unicellular offspring. By effecting the mutation rate, intraorganismal selection may have broad evolutionary implications; as an example, we consider its influence on the evolution of ploidy levels, finding that cell-lineage selection is more effective in haploids and tends to favor their evolution.     

Fixation probability of rare nonmutator and evolution of mutation rates

Although mutations drive the evolutionary process, the rates at which the mutations occur are themselves subject to evolutionary forces. Our purpose here is to understand the role of selection and random genetic drift in the evolution of mutation rates, and we address this question in asexual populations at mutation‐selection equilibrium neglecting selective sweeps. Using a multitype branching process, we calculate the fixation probability of a rare nonmutator in a large asexual population of mutators and find that a nonmutator is more likely to fix when the deleterious mutation rate of the mutator population is high. Compensatory mutations in the mutator population are found to decrease the fixation probability of a nonmutator when the selection coefficient is large. But, surprisingly, the fixation probability changes nonmonotonically with increasing compensatory mutation rate when the selection is mild. Using these results for the fixation probability and a drift‐barrier argument, we find a novel relationship between the mutation rates and the population size. We also discuss the time to fix the nonmutator in an adapted population of asexual mutators, and compare our results with experiments.     

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.     

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.     

Conditional coalescent trees with two mutation rates and their application to genomic instability

Humans have invested several genes in DNA repair and fidelity replication. To account for the disparity between the rarity of mutations in normal cells and the large number of mutations present in cancer, an hypothesis is that cancer cells must exhibit a mutator phenotype (genomic instability) during tumor progression, with the initiation of abnormal mutation rates caused by the loss of mismatch repair. In this study we introduce a stochastic model of mutation in tumor cells with the aim of estimating the amount of genomic instability due to the alteration of DNA repair genes. Our approach took into account the difficulties generated by sampling within tumoral clones and the fact that these clones must be difficult to isolate. We provide corrections to two classical statistics to obtain unbiased estimators of the raised mutation rate, and we show that large statistical errors may be associated with such estimators. The power of these new statistics to reject genomic instability is assessed and proved to increase with the intensity of mutation rates. In addition, we show that genomic instability cannot be detected unless the raised mutation rates exceed the normal rates by a factor of at least 1000.     

Directional mutation pressure,mutator mutations,and dynamics of molecular evolution

Using a general form of the directional mutation theory, this paper analyzes the effect of mutations in mutator genes on the G + C content of DNA, the frequency of substitution mutations, and evolutionary changes (cumulative mutations) under various degrees of selective constraints. Directional mutation theory predicts that when the mutational bias between A/T and G/C nucleotide pairs is equilibrated with the base composition of a neutral set of DNA nucleotides, the mutation frequency per gene will be much lower than the frequency immediately after the mutator mutation takes place. This prediction explains the wide variation of the DNA G + C content among unicellular organisms and possibly also the wide intragenomic heterogeneity of third codon positions for the genes of multicellular eukaryotes. The present analyses lead to several predictions that are not consistent with a number of the frequently held assumptions in the field of molecular evolution, including belief in a constant rate of evolution, symmetric branching of phylogenetic trees, the generality of higher mutation frequency for neutral sets of nucleotides, the notion that mutator mutations are generally deleterious because of their high mutation rates, and teleological explanations of DNA base composition. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Combining Magnetic Sorting of Mother Cells and Fluctuation Tests to Analyze Genome Instability During Mitotic Cell Aging in Saccharomyces cerevisiae

Saccharomyces cerevisiae has been an excellent model system for examining mechanisms and consequences of genome instability. Information gained from this yeast model is relevant to many organisms, including humans, since DNA repair and DNA damage response factors are well conserved across diverse species. However, S. cerevisiae has not yet been used to fully address whether the rate of accumulating mutations changes with increasing replicative (mitotic) age due to technical constraints. For instance, measurements of yeast replicative lifespan through micromanipulation involve very small populations of cells, which prohibit detection of rare mutations. Genetic methods to enrich for mother cells in populations by inducing death of daughter cells have been developed, but population sizes are still limited by the frequency with which random mutations that compromise the selection systems occur. The current protocol takes advantage of magnetic sorting of surface-labeled yeast mother cells to obtain large enough populations of aging mother cells to quantify rare mutations through phenotypic selections. Mutation rates, measured through fluctuation tests, and mutation frequencies are first established for young cells and used to predict the frequency of mutations in mother cells of various replicative ages. Mutation frequencies are then determined for sorted mother cells, and the age of the mother cells is determined using flow cytometry by staining with a fluorescent reagent that detects bud scars formed on their cell surfaces during cell division. Comparison of predicted mutation frequencies based on the number of cell divisions to the frequencies experimentally observed for mother cells of a given replicative age can then identify whether there are age-related changes in the rate of accumulating mutations. Variations of this basic protocol provide the means to investigate the influence of alterations in specific gene functions or specific environmental conditions on mutation accumulation to address mechanisms underlying genome instability during replicative aging.     

The Fixation Probability of Rare Mutators in Finite Asexual Populations

A mutator is an allele that increases the mutation rate throughout the genome by disrupting some aspect of DNA replication or repair. Mutators that increase the mutation rate by the order of 100-fold have been observed to spontaneously emerge and achieve high frequencies in natural populations and in long-term laboratory evolution experiments with Escherichia coli. In principle, the fixation of mutator alleles is limited by (i) competition with mutations in wild-type backgrounds, (ii) additional deleterious mutational load, and (iii) random genetic drift. Using a multiple-locus model and employing both simulation and analytic methods, we investigate the effects of these three factors on the fixation probability P_fix of an initially rare mutator as a function of population size N, beneficial and deleterious mutation rates, and the strength of mutations s. Our diffusion-based approximation for P_fix successfully captures effects ii and iii when selection is fast compared to mutation (). This enables us to predict the conditions under which mutators will be evolutionarily favored. Surprisingly, our simulations show that effect i is typically small for strong-effect mutators. Our results agree semiquantitatively with existing laboratory evolution experiments and suggest future experimental directions.     

Development of spontaneous mutations in somatic mamalian cells and DNA replication. II. Expression of gene mutations in cultured cells]

The dependence of spontaneous rate and phenotypic expression of 6-mercaptopurine resistance mutations of DNA replication synthesis was studied in cultured Chinese hamster cells. Spontaneous mutations arising with a constant rate per cell per time unit independently on DNA replication rate were shown to be expressed only in the course of subsequent cell divisions. The frequency of N-nitrosomethylurea induced mutations in cells with reduced and normal DNA replication rate is approximately the same. However, DNA replication synthesis is necessary for the phenotypic expression of both induced and spontaneous mutations. The causes of the phenotypic lag are discussed.     

Implications of genetic heterogeneity in cancer

DNA sequencing studies have established that many cancers contain tens of thousands of clonal mutations throughout their genomes, which is difficult to reconcile with the very low rate of mutation in normal human cells. This observation provides strong evidence for the mutator phenotype hypothesis, which proposes that a genome-wide elevation in the spontaneous mutation rate is an early step in carcinogenesis. An elevated mutation rate implies that cancers undergo continuous evolution, generating multiple subpopulations of cells that differ from one another in DNA sequence. The extensive heterogeneity in DNA sequence and continual tumor evolution that would occur in the context of a mutator phenotype have important implications for cancer diagnosis and therapy.     

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.     

Strong and weak male mutation bias at different sites in the primate genomes: insights from the human-chimpanzee comparison

Male mutation bias is a higher mutation rate in males than in females thought to result from the greater number of germ line cell divisions in males. If errors in DNA replication cause most mutations, then the magnitude of male mutation bias, measured as the male-to-female mutation rate ratio (alpha), should reflect the relative excess of male versus female germ line cell divisions. Evolutionary rates averaged among all sites in a sequence and compared between mammalian sex chromosomes were shown to be indeed higher in males than in females. However, it is presently unknown whether individual classes of substitutions exhibit such bias. To address this issue, we investigated male mutation bias separately at non-CpG and CpG sites using human-chimpanzee whole-genome alignments. We observed strong male mutation bias at non-CpG sites: alpha in the X-autosome comparison was approximately 6-7, which was similar to the male-to-female ratio in the number of germ line cell divisions. In contrast, mutations at CpG sites exhibited weak male mutation bias: alpha in the X-autosome comparison was only approximately 2-3. This is consistent with the methylation-induced and replication-independent mechanism of CpG transitions, which constitute the majority of mutations at CpG sites. Interestingly, our study also indicated weak male mutation bias for transversions at CpG sites, implying a spontaneous mechanism largely not associated with replication. Male mutation bias was equally strong at CpG and non-CpG sites located within unmethylated "CpG islands," suggesting the replication-dependent origin of these mutations. Thus, we found that the strength of male mutation bias is nonuniform in the primate genomes. Importantly, we discovered that male mutation bias depends on the proportion of CpG sites in the loci compared. This might explain the differences in the magnitude of primate male mutation bias observed among studies.     

Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations

We examined rates of DNA sequence evolution in 12 populations of Escherichia coli propagated in a glucose minimal medium for 20,000 generations. Previous work saw mutations mediated by mobile elements in these populations, but the extent of other genomic changes was not investigated. Four of the populations evolved defects in DNA repair and became mutators. Some 500 bp was sequenced in each of 36 genes for 50 clones, including 2 ancestral variants, 2 clones from each population at generation 10,000, and 2 from each at generation 20,000. Ten mutations were found in total, all point mutations including mostly synonymous substitutions and nonsynonymous polymorphisms; all 10 were found in mutator populations. We compared the observed sequence evolution to predictions based on different scenarios. The number of synonymous substitutions is lower than predicted from measured mutation rates in E. coli, but the number is higher than rates based on comparing E. coli and Salmonella genomes. Extrapolating to the entire genome, these data predict about 250 synonymous substitutions on average per mutator population, but only about 3 synonymous substitutions per nonmutator population, during 20,000 generations. These data illustrate the challenge of finding sequence variation among bacterial isolates that share such a recent ancestor. However, this limited variation also provides a useful baseline for research aimed at finding the beneficial substitutions in these populations.