Mutators in space: the dynamics of high-mutability clones in a two-patch model

Clones of bacteria possessing high-mutability rates (or mutators) are being observed in an increasing number of species. In a constant environment most mutations are deleterious, and hence the spontaneous mutation rate is generally low. However, mutators may play an important role in the adaptation of organisms to changing environments. To date, theoretical work has focused on temporal variability in the environment, implicitly assuming that environmental conditions are constant through space. Here, we develop a two-patch model to investigate how spatiotemporal environmental variability and dispersal might influence mutator dynamics. Environmental conditions in each patch fluctuate between two states; the rate of fluctuation varies in each patch at differing phase angles. We find that at low and intermediate rates of fluctuation, an increase in dispersal results in a decrease in the density of mutators. However, at high rates of environmental change, dispersal causes an increase in mutator density. For all frequencies of environmental fluctuation these trends are enhanced as the phase angle approaches 180 degrees. We argue that future work, both empirical and theoretical, is needed to improve our understanding of how spatiotemporal variability impacts on mutator densities and dynamics.     

The evolution of mutator genes in bacterial populations: the roles of environmental change and timing

Recent studies have found high frequencies of bacteria with increased genomic rates of mutation in both clinical and laboratory populations. These observations may seem surprising in light of earlier experimental and theoretical studies. Mutator genes (genes that elevate the genomic mutation rate) are likely to induce deleterious mutations and thus suffer an indirect selective disadvantage; at the same time, bacteria carrying them can increase in frequency only by generating beneficial mutations at other loci. When clones carrying mutator genes are rare, however, these beneficial mutations are far more likely to arise in members of the much larger nonmutator population. How then can mutators become prevalent? To address this question, we develop a model of the population dynamics of bacteria confronted with ever-changing environments. Using analytical and simulation procedures, we explore the process by which initially rare mutator alleles can rise in frequency. We demonstrate that subsequent to a shift in environmental conditions, there will be relatively long periods of time during which the mutator subpopulation can produce a beneficial mutation before the ancestral subpopulations are eliminated. If the beneficial mutation arises early enough, the overall frequency of mutators will climb to a point higher than when the process began. The probability of producing a subsequent beneficial mutation will then also increase. In this manner, mutators can increase in frequency over successive selective sweeps. We discuss the implications and predictions of these theoretical results in relation to antibiotic resistance and the evolution of mutation rates.     

Dynamics of Adaptation During Three Years of Evolution Under Long-Term Stationary Phase

Many bacterial species that cannot sporulate, such as the model bacterium Escherichia coli, can nevertheless survive for years, following exhaustion of external resources, in a state termed long-term stationary phase (LTSP). Here we describe the dynamics of E. coli adaptation during the first three years spent under LTSP. We show that during this time, E. coli continuously adapts genetically through the accumulation of mutations. For nonmutator clones, the majority of mutations accumulated appear to be adaptive under LTSP, reflected in an extremely convergent pattern of mutation accumulation. Despite the rapid and convergent manner in which populations adapt under LTSP, they continue to harbor extensive genetic variation. The dynamics of evolution of mutation rates under LTSP are particularly interesting. The emergence of mutators affects overall mutation accumulation rates as well as the mutational spectra and the ultimate spectrum of adaptive alleles acquired under LTSP. With time, mutators can evolve even higher mutation rates through the acquisition of additional mutation rate–enhancing mutations. Different mutator and nonmutator clones within a single population and time point can display extreme variation in their mutation rates, resulting in differences in both the dynamics of adaptation and their associated deleterious burdens. Despite these differences, clones that vary greatly in their mutation rates tend to coexist within their populations for many years, under LTSP.     

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.     

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.     

Ploidy controls the success of mutators and nature of mutations during budding yeast evolution

BACKGROUND: We used the budding yeast Saccharomyces cerevisiae to ask how elevated mutation rates affect the evolution of asexual eukaryotic populations. Mismatch repair defective and nonmutator strains were competed during adaptation to four laboratory environments (rich medium, low glucose, high salt, and a nonfermentable carbon source). RESULTS: In diploids, mutators have an advantage over nonmutators in all conditions, and mutators that win competitions are on average fitter than nonmutator winners. In contrast, haploid mutators have no advantage when competed against haploid nonmutators, and haploid mutator winners are less fit than nonmutator winners. The diploid mutator winners were all superior to their ancestors both in the condition they had adapted to, and in two of the other conditions. This phenotype was due to a mutation or class of mutations that confers a large growth advantage during the respiratory phase of yeast cultures that precedes stationary phase. This generalist mutation(s) was not selected in diploid nonmutator strains or in haploid strains, which adapt primarily by fixing specialist (condition-specific) mutations. In diploid mutators, such mutations also occur, and the majority accumulates after the fixation of the generalist mutation. CONCLUSIONS: We conclude that the advantage of mutators depends on ploidy and that diploid mutators can give rise to beneficial mutations that are inaccessible to nonmutators and haploid mutators.     

FIXATION OF MUTATORS IN ASEXUAL POPULATIONS: THE ROLE OF GENETIC DRIFT AND EPISTASIS

We study the evolutionary dynamics of an asexual population of nonmutators and mutators on a class of epistatic fitness landscapes. We consider the situation in which all mutations are deleterious and mutators are produced from nonmutators continually at a constant rate. We find that in an infinitely large population, a minimum nonmutator‐to‐mutator conversion rate is required to fix the mutators but an arbitrarily small conversion rate results in the fixation of mutators in a finite population. We calculate analytical expressions for the mutator fraction at mutation‐selection balance and fixation time for mutators in a finite population when the difference between the mutation rate for mutator and nonmutator is smaller (regime I) and larger (regime II) than the selection coefficient. Our main result is that in regime I, the mutator fraction and the fixation time are independent of epistasis but in regime II, mutators are rarer and take longer to fix when the decrease in fitness with the number of deleterious mutations occurs at an accelerating rate (synergistic epistasis) than at a diminishing rate (antagonistic epistasis). Our analytical results are compared with numerics and their implications are discussed.     

Interplay between Pleiotropy and Secondary Selection Determines Rise and Fall of Mutators in Stress Response

Mutators are clones whose mutation rate is about two to three orders of magnitude higher than the rate of wild-type clones and their roles in adaptive evolution of asexual populations have been controversial. Here we address this problem by using an ab initio microscopic model of living cells, which combines population genetics with a physically realistic presentation of protein stability and protein-protein interactions. The genome of model organisms encodes replication controlling genes (RCGs) and genes modeling the mismatch repair (MMR) complexes. The genotype-phenotype relationship posits that the replication rate of an organism is proportional to protein copy numbers of RCGs in their functional form and there is a production cost penalty for protein overexpression. The mutation rate depends linearly on the concentration of homodimers of MMR proteins. By simulating multiple runs of evolution of populations under various environmental stresses—stationary phase, starvation or temperature-jump—we find that adaptation most often occurs through transient fixation of a mutator phenotype, regardless of the nature of stress. By contrast, the fixation mechanism does depend on the nature of stress. In temperature jump stress, mutators take over the population due to loss of stability of MMR complexes. In contrast, in starvation and stationary phase stresses, a small number of mutators are supplied to the population via epigenetic stochastic noise in production of MMR proteins (a pleiotropic effect), and their net supply is higher due to reduced genetic drift in slowly growing populations under stressful environments. Subsequently, mutators in stationary phase or starvation hitchhike to fixation with a beneficial mutation in the RCGs, (second order selection) and finally a mutation stabilizing the MMR complex arrives, returning the population to a non-mutator phenotype. Our results provide microscopic insights into the rise and fall of mutators in adapting finite asexual populations.     

Enrichment and elimination of mutY mutators in Escherichia coli populations

The kinetics of mutator sweeps was followed in two independent populations of Escherichia coli grown for up to 350 generations in glucose-limited continuous culture. A rapid elevation of mutation rates was observed in both populations within 120-150 generations, as was apparent from major increases in the proportion of the populations with unselected mutations in fhuA. The increase in mutation rates was due to sweeps by mutY mutators. In both cultures, the enrichment of mutators resulted from hitchhiking with identified beneficial mutations increasing fitness under glucose limitation; mutY hitchhiked with mgl mutations in one culture and ptsG in the other. In both cases, mutators were enriched to constitute close to 100% of the population before a periodic selection event reduced the frequency of unselected mutations and mutators in the cultures. The high proportion of mutators persisted for 150 generations in one population but began to be eliminated within 50 generations in the other. The persistence of mutator, as well as experimental data showing that mutY bacteria were as fit as near-isogenic mutY(+) bacteria in competition experiments, suggest that mutator load by deleterious mutations did not explain the rapidly diminishing proportion of mutators in the populations. The nonmutators sweeping out mutators were also unlikely to have arisen by reversion or antimutator mutations; the mutY mutations were major deletions in each case and the bacteria sweeping out mutators contained intact mutY. By following mgl allele frequencies in one population, we discovered that mutators were outcompeted by bacteria that had rare mgl mutations previously as well as additional beneficial mutation(s). The pattern of appearance of mutY, but not its elimination, conforms to current models of mutator sweeps in bacterial populations. A mutator with a narrow mutational spectrum like mutY may be lost if the requirement for beneficial mutations is for changes other than GC --> TA transversions. Alternatively, epistatic interactions between mutator mutation and beneficial mutations need to be postulated to explain mutator elimination.     

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.     

The balance between mutators and nonmutators in asexual populations

Mutator alleles, which elevate an individual's mutation rate from 10 to 10,000-fold, have been found at high frequencies in many natural and experimental populations. Mutators are continually produced from nonmutators, often due to mutations in mismatch-repair genes. These mutators gradually accumulate deleterious mutations, limiting their spread. However, they can occasionally hitchhike to high frequencies with beneficial mutations. We study the interplay between these effects. We first analyze the dynamics of the balance between the production of mutator alleles and their elimination due to deleterious mutations. We find that when deleterious mutation rates are high in mutators, there will often be many "young," recently produced mutators in the population, and the fact that deleterious mutations only gradually eliminate individuals from a population is important. We then consider how this mutator-nonmutator balance can be disrupted by beneficial mutations and analyze the circumstances in which fixation of mutator alleles is likely. We find that dynamics is crucial: even in situations where selection on average acts against mutators, so they cannot stably invade, the mutators can still occasionally generate beneficial mutations and hence be important to the evolution of the population.     

IMPACT OF BACTERIAL MUTATION RATE ON COEVOLUTIONARY DYNAMICS BETWEEN BACTERIA AND PHAGES

Mutator bacteria are frequently found in natural populations of bacteria and although coevolution with parasitic viruses (phages) is thought to be one reason for their persistence, it remains unclear how the presence of mutators affects coevolutionary dynamics. We hypothesized that phages must themselves adapt more rapidly or go extinct, in the face of rapidly evolving mutator bacteria. We compared the coevolutionary dynamics of wild‐type Pseudomonas fluorescens SBW25 with a lytic phage to the dynamics of an isogenic mutator of P. fluorescens SBW25 together with the same phage. At the beginning of the experiment both wild‐type bacteria and mutator bacteria coevolved with phages. However, mutators rapidly evolved higher levels of sympatric resistance to phages. The phages were unable to “keep‐up” with the mutator bacteria, and these rates of coevolution declined to less than the rates of coevolution between the phages and wild‐type bacteria. By the end of the experiment, the sympatric resistance of the mutator bacteria was not significantly different to the sympatric resistance of the wild‐type bacteria. This suggests that the importance of mutators in the coevolutionary interactions with a particular phage population is likely to be short‐lived. More generally, the results demonstrate that coevolving enemies may escape from Red‐Queen dynamics.     

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.     

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.     

Competition between the dam mutator and the isogenic wild-type of Escherichia coli

Previous studies have shown that the mutT, mutH, mutL and mutS mutators of Escherichia coli confer a marked selective advantage on their respective hosts in competition with otherwise isogenic wild-type strains. We have conducted competition experiments between dam- and dam+ strains of Escherichia coli and have found that dam mutator strains are negatively selected. Although dam- is the first mutator to have a lower fitness than wild-type under chemostat conditions our result does not contradict the hypothesis that increased mutation rates are of evolutionary advantage under environmental stress conditions. Only in the special case of dam- does the advantage of higher mutation rates not outweigh the disadvantage due to the dam- -caused heavy pleiotropic effects.     

Mutation Rate Evolution in Replicator Dynamics

The mutation rate of an organism is itself evolvable. In stable environments, if faithful replication is costless, theory predicts that mutation rates will evolve to zero. However, positive mutation rates can evolve in novel or fluctuating environments, as analytical and empirical studies have shown. Previous work on this question has focused on environments that fluctuate independently of the evolving population. Here we consider fluctuations that arise from frequency-dependent selection in the evolving population itself. We investigate how the dynamics of competing traits can induce selective pressure on the rates of mutation between these traits. To address this question, we introduce a theoretical framework combining replicator dynamics and adaptive dynamics. We suppose that changes in mutation rates are rare, compared to changes in the traits under direct selection, so that the expected evolutionary trajectories of mutation rates can be obtained from analysis of pairwise competition between strains of different rates. Depending on the nature of frequency-dependent trait dynamics, we demonstrate three possible outcomes of this competition. First, if trait frequencies are at a mutation–selection equilibrium, lower mutation rates can displace higher ones. Second, if trait dynamics converge to a heteroclinic cycle—arising, for example, from “rock-paper-scissors” interactions—mutator strains succeed against non-mutators. Third, in cases where selection alone maintains all traits at positive frequencies, zero and nonzero mutation rates can coexist indefinitely. Our second result suggests that relatively high mutation rates may be observed for traits subject to cyclical frequency-dependent dynamics.     

Migration promotes mutator alleles in subdivided populations

Mutator alleles that elevate the genomic mutation rate may invade nonrecombining populations by hitchhiking with beneficial mutations. Mutators have been repeatedly observed to take over adapting laboratory populations and have been found at high frequencies in both microbial pathogen and cancer populations in nature. Recently, we have shown that mutators are only favored by selection in sufficiently large populations and transition to being disfavored as population size decreases. This population size‐dependent sign inversion in selective effect suggests that population structure may also be an important determinant of mutation rate evolution. Although large populations may favor mutators, subdividing such populations into sufficiently small subpopulations (demes) might effectively inhibit them. On the other hand, migration between small demes that otherwise inhibit hitchhiking may promote mutator fixation in the whole metapopulation. Here, we use stochastic, agent‐based simulations and evolution experiments with the yeast Saccharomyces cerevisiae to show that mutators can, indeed, be favored by selection in subdivided metapopulations composed of small demes connected by sufficient migration. In fact, we show that population structure plays a previously unsuspected role in promoting mutator success in subdivided metapopulations when migration is rare.     

The effect of spatial structure on adaptation in Escherichia coli

Populations of organisms are generally organized in a given spatial structure. However, the vast majority of population genetic studies are based on populations in which every individual competes globally. Here we use experimental evolution in Escherichia coli to directly test a recently made prediction that spatial structure slows down adaptation and that this cost increases with the mutation rate. This was studied by comparing populations of different mutation rates adapting to a liquid (unstructured) medium with populations that evolved in a Petri dish on solid (structured) medium. We find that mutators adapt faster to both environments and that adaptation is slower if there is spatial structure. We observed no significant difference in the cost of structure between mutator and wild-type populations, which suggests that clonal interference is intense in both genetic backgrounds.     

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.