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.     

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.     

Metabolic Erosion Primarily Through Mutation Accumulation,and Not Tradeoffs,Drives Limited Evolution of Substrate Specificity in Escherichia coli

Evolutionary adaptation to a constant environment is often accompanied by specialization and a reduction of fitness in other environments. We assayed the ability of the Lenski Escherichia coli populations to grow on a range of carbon sources after 50,000 generations of adaptation on glucose. Using direct measurements of growth rates, we demonstrated that declines in performance were much less widespread than suggested by previous results from Biolog assays of cellular respiration. Surprisingly, there were many performance increases on a variety of substrates. In addition to the now famous example of citrate, we observed several other novel gains of function for organic acids that the ancestral strain only marginally utilized. Quantitative growth data also showed that strains with a higher mutation rate exhibited significantly more declines, suggesting that most metabolic erosion was driven by mutation accumulation and not by physiological tradeoffs. These reductions in growth by mutator strains were ameliorated by growth at lower temperature, consistent with the hypothesis that this metabolic erosion is largely caused by destabilizing mutations to the associated enzymes. We further hypothesized that reductions in growth rate would be greatest for substrates used most differently from glucose, and we used flux balance analysis to formulate this question quantitatively. To our surprise, we found no significant relationship between decreases in growth and dissimilarity to glucose metabolism. Taken as a whole, these data suggest that in a single resource environment, specialization does not mainly result as an inevitable consequence of adaptive tradeoffs, but rather due to the gradual accumulation of disabling mutations in unused portions of the genome.     

Escape from growth restriction in small colony variants of Salmonella typhimurium by gene amplification and mutation

Antibiotic resistance in bacteria is generally associated with fitness costs that often can be reduced by second-site compensatory mutations. Here, we examined how a protamine-resistant small colony variant of Salmonella typhimurium adapts to the growth reduction conferred by a resistance mutation in hemC (encoding a haem-biosynthesis enzyme). We show that adaptation occurs in a multi-step process where fitness is successively increased. Thus, the initial adaptive response was selection for an unstable gene amplification of the mutant hemC gene that provided a small fitness increase. Fitness was increased further by a mutation that restored HemC function in one gene copy, relaxing selection for the amplification. Subsequently, the amplification segregated back to the haploid state and even higher fitness. The end result was in most cases mutant strains with a hemC sequence different from that of the wild-type strain. These findings suggest that gene amplification facilitates adaptive evolution. A higher gene dosage increases the target size for compensatory mutations and improves fitness of the cell, thereby allowing an increase in the population size, further increasing the probability of a subsequent stable mutation. Our results provide a novel genetic basis for growth compensation in small colony variants.     

Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria.

Selection of mutator alleles, increasing the mutation rate up to 10, 000-fold, has been observed during in vitro experimental evolution. This spread is ascribed to the hitchhiking of mutator alleles with favorable mutations, as demonstrated by a theoretical model using selective parameters corresponding to such experiments. Observations of unexpectedly high frequencies of mutators in natural isolates suggest that the same phenomenon could occur in the wild. But it remains questionable whether realistic in natura parameter values could also result in selection of mutators. In particular, the main parameters of adaptation, the size of the adapting population and the height and steepness of the adaptive peak characterizing adaptation, are very variable in nature. By simulation approach, we studied the effect of these parameters on the selection of mutators in asexual populations, assuming additive fitness. We show that the larger the population size, the more likely the fixation of mutator alleles. At a large population size, at least four adaptive mutations are needed for mutator fixation; moreover, under stronger selection stronger mutators are selected. We propose a model based on multiple mutations to illustrate how second-order selection can optimize population fitness when few favorable mutations are required for adaptation.     

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.     

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.     

Long-term effects of inducible mutagenic DNA repair on relative fitness and phenotypic diversification in Pseudomonas cichorii 302959

Mutagenic DNA repair (MDR) employs low-fidelity DNA polymerases capable of replicating past DNA lesions resulting from exposure to high-energy ultraviolet radiation (UVR). MDR confers UVR tolerance and activation initiates a transient mutator phenotype that may provide opportunities for adaptation. To investigate the potential role of MDR in adaptation, we have propagated parallel lineages of the highly mutable epiphytic plant pathogen Pseudomonas cichorii 302959 with daily UVR activation (UVR lineages) for ~500 generations. Here we examine those lineages through the measurement of relative fitness and observation of distinct colony morphotypes that emerged. Isolates and population samples from UVR lineages displayed gains in fitness relative to the ancestor despite increased rates of inducible mutation to rifampicin resistance. Regular activation of MDR resulted in the maintenance of genetic diversity within UVR lineages, including the reproducible diversification and coexistence of “round” and “fuzzy” colony morphotypes. These results suggest that inducible mutability may present a reasonable strategy for adaptive evolution in stressful environments by contributing to gains in relative fitness and diversification.     

Estimates of the rate and distribution of fitness effects of spontaneous mutation in Saccharomyces cerevisiae

The per-genome, per-generation rate of spontaneous mutation affecting fitness (U) and the mean fitness cost per mutation (s) are important parameters in evolutionary genetics, but have been estimated for few species. We estimated U and sh (the heterozygous effect of mutations) for two diploid yeast strains differing only in the DNA mismatch-repair deficiency used to elevate the mutation rate in one (mutator) strain. Mutations were allowed to accumulate in 50 replicate lines of each strain, during 36 transfers of randomly chosen single colonies (approximately 600 generations). Among wild-type lines, fitnesses were bimodal, with one mode showing no change in mean fitness. The other mode showed a mean 29.6% fitness decline and the petite phenotype, usually caused by partial deletion of the mitochondrial genome. Excluding petites, maximum-likelihood estimates adjusted for the effect of selection were U = 9.5 x 10(-5) and sh = 0.217 for the wild type. Among the mutator lines, the best fit was obtained with 0.005 < or = U < or = 0.94 and 0.049 > or = sh > or = 0.0003. Like other recently tested model organisms, wild-type yeast have low mutation rates, with high mean fitness costs per mutation. Inactivation of mismatch repair increases the frequency of slightly deleterious mutations by approximately two orders of magnitude.     

Genetic adaptation of Pseudomonas aeruginosa to the airways of cystic fibrosis patients is catalyzed by hypermutation

In previous work (E. E. Smith, D. G. Buckley, Z. Wu, C. Saenphimmachack, L. R. Hoffman, D. A. D'Argenio, S. I. Miller, B. W. Ramsey, D. P. Speert, S. M. Moskowitz, J. L. Burns, R. Kaul, and M. V. Olson, Proc. Natl. Acad. Sci. USA 103:8487-8492, 2006) it was shown that Pseudomonas aeruginosa undergoes intense genetic adaptation during chronic respiratory infection (CRI) in cystic fibrosis (CF) patients. We used the same collection of isolates to explore the role of hypermutation in this process, since one of the hallmarks of CRI is the high prevalence of DNA mismatch repair (MMR) system-deficient mutator strains. The presence of mutations in 34 genes (many of them positively linked to adaptation in CF patients) in the study collection of 90 P. aeruginosa isolates obtained longitudinally from 29 CF patients was not homogeneous; on the contrary, mutations were significantly concentrated in the mutator lineages, which represented 17% of the isolates (87% MMR deficient). While sequential nonmutator lineages acquired a median of only 0.25 mutation per year of infection, mutator lineages accumulated more than 3 mutations per year. On the whole-genome scale, data for the first fully sequenced late CF isolate, which was also shown to be an MMR-deficient mutator, also support these findings. Moreover, for the first time the predicted amplification of mutator populations due to hitchhiking with adaptive mutations in the course of natural human infections is clearly documented. Interestingly, increased accumulation of mutations in mutator lineages was not a consequence of overrepresentation of mutations in genes involved in antimicrobial resistance, the only adaptive trait linked so far to hypermutation in CF patients, demonstrating that hypermutation also plays a major role in P. aeruginosa genome evolution and adaptation during CRI.     

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.     

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.     

Beneficial mutation selection balance and the effect of linkage on positive selection

When beneficial mutations are rare, they accumulate by a series of selective sweeps. But when they are common, many beneficial mutations will occur before any can fix, so there will be many different mutant lineages in the population concurrently. In an asexual population, these different mutant lineages interfere and not all can fix simultaneously. In addition, further beneficial mutations can accumulate in mutant lineages while these are still a minority of the population. In this article, we analyze the dynamics of such multiple mutations and the interplay between multiple mutations and interference between clones. These result in substantial variation in fitness accumulating within a single asexual population. The amount of variation is determined by a balance between selection, which destroys variation, and beneficial mutations, which create more. The behavior depends in a subtle way on the population parameters: the population size, the beneficial mutation rate, and the distribution of the fitness increments of the potential beneficial mutations. The mutation-selection balance leads to a continually evolving population with a steady-state fitness variation. This variation increases logarithmically with both population size and mutation rate and sets the rate at which the population accumulates beneficial mutations, which thus also grows only logarithmically with population size and mutation rate. These results imply that mutator phenotypes are less effective in larger asexual populations. They also have consequences for the advantages (or disadvantages) of sex via the Fisher-Muller effect; these are discussed briefly.     

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.     

Mutator Mutations Enhance Tumorigenic Efficiency across Fitness Landscapes

Background

Tumorigenesis requires multiple genetic changes. Mutator mutations are mutations that increase genomic instability, and according to the mutator hypothesis, accelerate tumorigenesis by facilitating oncogenic mutations. Alternatively, repeated lineage selection and expansion without increased mutation frequency may explain observed cancer incidence. Mutator lineages also risk increased deleterious mutations, leading to extinction, thus providing another counterargument to the mutator hypothesis. Both selection and extinction involve changes in lineage fitness, which may be represented as “trajectories” through a “fitness landscape” defined by genetics and environment.

Methodology/Principal Findings

Here I systematically analyze the relative efficiency of tumorigenesis with and without mutator mutations by evaluating archetypal fitness trajectories using deterministic and stochastic mathematical models. I hypothesize that tumorigenic mechanisms occur clinically in proportion to their relative efficiency. This work quantifies the relative importance of mutator pathways as a function of experimentally measurable parameters, demonstrating that mutator pathways generally enhance efficiency of tumorigenesis. An optimal mutation rate for tumor evolution is derived, and shown to differ from that for species evolution.

Conclusions/Significance

The models address the major counterarguments to the mutator hypothesis, confirming that mutator mechanisms are generally more efficient routes to tumorigenesis than non-mutator mechanisms. Mutator mutations are more likely to occur early, and to occur when more oncogenic mutations are required to create a tumor. Mutator mutations likely occur in a minority of premalignant lesions, but these mutator premalignant lesions are disproportionately likely to develop into malignant tumors. Tumor heterogeneity due to mutator mutations may contribute to therapeutic resistance, and the degree of heterogeneity of tumors may need to be considered when therapeutic strategies are devised. The model explains and predicts important biological observations in bacterial and mouse systems, as well as clinical observations.     

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.     

Compensatory adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium

Most chromosomal mutations that cause antibiotic resistance impose fitness costs on the bacteria. This biological cost can often be reduced by compensatory mutations. In Salmonella typhimurium, the nucleotide substitution AAA42 --> AAC in the rpsL gene confers resistance to streptomycin. The resulting amino acid substitution (K42N) in ribosomal protein S12 causes an increased rate of ribosomal proofreading and, as a result, the rate of protein synthesis, bacterial growth and virulence are decreased. Eighty-one independent lineages of the low-fitness, K42N mutant were evolved in the absence of antibiotic to ameliorate the costs. From the rate of fixation of compensated mutants and their fitness, the rate of compensatory mutations was estimated to be > or = 10-7 per cell per generation. The size of the population bottleneck during evolution affected fitness of the adapted mutants: a larger bottleneck resulted in higher average fitness. Only four of the evolved lineages contained streptomycin-sensitive revertants. The remaining 77 lineages contained mutants that were still fully streptomycin resistant, had retained the original resistance mutation and also acquired compensatory mutations. Most of the compensatory mutations, resulting in at least 35 different amino acid substitutions, were novel single-nucleotide substitutions in the rpsD, rpsE, rpsL or rplS genes encoding the ribosomal proteins S4, S5, S12 and L19 respectively. Our results show that the deleterious effects of a resistance mutation can be compensated by an unexpected variety of mutations.     

Mutator dynamics in fluctuating environments

Populations with high mutation rates (mutator clones) are being found in increasing numbers of species, and a clear link is being established between the presence of mutator clones and drug resistance. Mutator clones exist despite the fact that in a constant environment most mutations are deleterious, with the spontaneous mutation rate generally held at a low value. This implies that mutator clones have an important role in the adaptation of organisms to changing environments. Our study examines how mutator dynamics vary according to the frequency of environmental fluctuations. Although recent studies have considered a single environmental switch, here we investigate mutator dynamics in a regularly varying environment, seeking to mimic conditions present, for example, under certain drug or pesticide regimes. Our model provides four significant new insights. First, the results demonstrate that mutators are most prevalent under intermediate rates of environmental change. When the environment oscillates more rapidly, mutators are unable to provide sufficient adaptability to keep pace with the frequent changes in selection pressure and, instead, a population of 'environmental generalists' dominates. Second, our findings reveal that mutator dynamics may be complex, exhibiting limit cycles and chaos. Third, we demonstrate that when each beneficial mutation provides a greater gain in fitness, mutators achieve higher densities in more rapidly fluctuating environments. Fourth, we find that mutators of intermediate strength reach higher densities than very weak or strong mutators.     

Hypermutation and stress adaptation in bacteria

Hypermutability is a phenotype characterized by a moderate to high elevation of spontaneous mutation rates and could result from DNA replication errors, defects in error correction mechanisms and many other causes. The elevated mutation rates are helpful to organisms to adapt to sudden and unforeseen threats to survival. At the same time hypermutability also leads to the generation of many deleterious mutations which offset its adaptive value and therefore disadvantageous. Nevertheless, it is very common in nature, especially among clinical isolates of pathogens. Hypermutability is inherited by indirect (second order) selection along with the beneficial mutations generated. At large population sizes and high mutation rates many cells in the population could concurrently acquire beneficial mutations of varying adaptive (fitness) values. These lineages compete with the ancestral cells and also among themselves for fixation. The one with the ‘fittest’ mutation gets fixed ultimately while the others are lost. This has been called ‘clonal interference’ which puts a speed limit on adaptation. The original clonal interference hypothesis has been modified recently. Nonheritable (transient) hypermtability conferring significant adaptive benefits also occur during stress response although its molecular basis remains controversial. The adaptive benefits of heritable hypermutability are discussed with emphasis on host–pathogen interactions.