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.     

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.     

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.     

ASYMMETRIC PATCH SIZE DISTRIBUTION LEADS TO DISRUPTIVE SELECTION ON DISPERSAL

Numerous models have been designed to understand how dispersal ability evolves when organisms live in a fragmented landscape. Most of them predict a single dispersal rate at evolutionary equilibrium, and when diversification of dispersal rates has been predicted, it occurs as a response to perturbation or environmental fluctuation regimes. Yet abundant variation in dispersal ability is observed in natural populations and communities, even in relatively stable environments. We show that this diversification can operate in a simple island model without temporal variability: disruptive selection on dispersal occurs when the environment consists of many small and few large patches, a common feature in natural spatial systems. This heterogeneity in patch size results in a high variability in the number of related patch mates by individual, which, in turn, triggers disruptive selection through a high per capita variance of inclusive fitness. Our study provides a likely, parsimonious and testable explanation for the diversity of dispersal rates encountered in nature. It also suggests that biological conservation policies aiming at preserving ecological communities should strive to keep the distribution of patch size sufficiently asymmetric and variable.     

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.     

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.     

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.     

Scaling population responses to spatial environmental variability in advection-dominated systems

We model the spatial dynamics of an open population of organisms that disperse solely through advection in order to understand responses to multiscale environmental variability. We show that the distance over which a population responds to a localized perturbation, called the response length, can be characterized as an organisms average lifetime dispersal distance, unless there is strong density‐dependence in demographic or dispersal rates. Continuous spatial fluctuations in demographic rates at scales smaller than the response length will be largely averaged in the population distribution, whereas those in per capita emigration rates will be strongly tracked. We illustrate these results using a parameterized example to show how responses to environmental variability may differ in streams with different average current velocities. Our model suggests an approach to linking local dynamics dominated by dispersal processes to larger‐scale dynamics dominated by births and deaths.     

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.     

A toad more traveled: the heterogeneous invasion dynamics of cane toads in Australia

To predict the spread of invasive species, we need to understand the mechanisms that underlie their range expansion. Assuming random diffusion through homogeneous environments, invasions are expected to progress at a constant rate. However, environmental heterogeneity is expected to alter diffusion rates, especially by slowing invasions as populations encounter suboptimal environmental conditions. Here, we examine how environmental and landscape factors affect the local invasion speeds of cane toads (Chaunus [Bufo] marinus) in Australia. Using high-resolution cane toad data, we demonstrate heterogeneous regional invasion dynamics that include both decelerating and accelerating range expansions. Toad invasion speed increased in regions characterized by high temperatures, heterogeneous topography, low elevations, dense road networks, and high patch connectivity. Regional increases in the toad invasion rate might be caused by environmental conditions that facilitate toad reproduction and movement, by the evolution of long-distance dispersal ability, or by some combination of these factors. In any case, theoretical predictions that neglect environmental influences on dispersal at multiple spatial scales may prove to be inaccurate. Early predictions of cane toad range expansion rates that assumed constant diffusion across homogeneous landscapes already have been proved wrong. Future attempts to predict range dynamics for invasive species should consider heterogeneity in (1) the environmental factors that determine dispersal rates and (2) the mobility of invasive populations because dispersal-relevant traits can evolve in exotic habitats. As an invasive species spreads, it is likely to encounter conditions that influence dispersal rates via one or both of these mechanisms.     

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.     

Dispersal promotes compensatory dynamics and stability in forced metacommunities

Understanding the factors that govern the stability of populations and communities has gained increasing importance as habitat fragmentation and environmental perturbations continue to escalate due to human activities. Dispersal is commonly viewed as essential to the maintenance of diversity in spatially subdivided communities, but few experiments have explored how dispersal interacts with the spatiotemporal components of environmental perturbations to determine community-level stability. We examined these processes using an experimental planktonic system composed of three competing species of zooplankton. We subjected zooplankton metacommunities to varying levels of dispersal and pH perturbations that varied in their degree of spatial synchrony. We show that dispersal can reverse the destabilizing effects of environmental forcing when perturbations are spatially asynchronous. Asynchrony in pH perturbations generated spatially and temporally varying species refugia that promoted source-sink dynamics and allowed prolonged persistence of zooplankton species that were otherwise extirpated in synchronously varying metacommunities. This, in turn, increased local species diversity, promoted compensatory population dynamics, and enhanced local community-level stability. Our results indicate that patterns of spatial covariation in environmental variability are critical to predicting the effects of dispersal on the dynamics and persistence of communities.     

Spatial ecology of large herbivore populations

Models of the dynamics of large herbivore populations represent density feedbacks on the population growth rate either directly or indirectly through interactions with vegetation resources. Neither approach incorporates the spatial heterogeneity that is an essential feature of most natural environments, and modifies the population dynamics generated. This is especially true for large herbivores exploiting food resources that are rooted in space but temporally variable in quantity and quality both seasonally and annually. In this review I explore how environmental variation at different spatiotemporal scales influences the abundance of herbivore populations controlled via resources, predators or social mechanisms. Changes in abundance can be spatially disparate and dependent on different resource components at different stages of the seasonal cycle, including buffer resources restricting population crashes in extremely adverse years. GPS telemetry enables movement responses generating spatial patterns to be documented in fine spatiotemporal detail, including migration and dispersal. Models incorporating spatial heterogeneity either implicitly or explicitly are outlined, exemplifying how herbivores cope with temporal variability by exploiting spatial variability in resources and conditions. Global human dominance is generating widened climatic variation while opportunities for herbivore movements are becoming constricted. Theoretical population ecologists need to shift their focus from the workings of demographic structure towards effects of changing environmental contexts, in order to project the likely trajectories of large herbivore populations through the Anthropocene.     

Effect of predator density dependent dispersal of prey on stability of a predator-prey system

This work presents a predator-prey Lotka-Volterra model in a two patch environment. The model is a set of four ordinary differential equations that govern the prey and predator population densities on each patch. Predators disperse with constant migration rates, while prey dispersal is predator density-dependent. When the predator density is large, the dispersal of prey is more likely to occur. We assume that prey and predator dispersal is faster than the local predator-prey interaction on each patch. Thus, we take advantage of two time scales in order to reduce the complete model to a system of two equations governing the total prey and predator densities. The stability analysis of the aggregated model shows that a unique strictly positive equilibrium exists. This equilibrium may be stable or unstable. A Hopf bifurcation may occur, leading the equilibrium to be a centre. If the two patches are similar, the predator density dependent dispersal of prey has a stabilizing effect on the predator-prey system.     

Dispersal evolution in temporally variable environments: implications for plant range dynamics

Dispersal—the movement of an individual from the site of birth to a different site for reproduction—is an ecological and evolutionary driver of species ranges that shapes patterns of colonization, connectivity, gene flow, and adaptation. In plants, the traits that influence dispersal often vary within and among species, are heritable, and evolve in response to the fitness consequences of moving through heterogeneous landscapes. Spatial and temporal variation in the quality and quantity of habitat are important sources of selection on dispersal strategies across species ranges. While recent reviews have evaluated the interactions between spatial variation in habitat and dispersal dynamics, the extent to which geographic variation in temporal variability can also shape range-wide patterns in dispersal traits has not been synthesized. In this paper, we summarize key predictions from metapopulation models that evaluate how dispersal evolves in response to spatial and temporal habitat variability. Next, we compile empirical data that quantify temporal variability in plant demography and patterns of dispersal trait variation across species ranges to evaluate the hypothesis that higher temporal variability favors increased dispersal at plant range limits. We found some suggestive evidence supporting this hypothesis while more generally identifying a major gap in empirical work evaluating plant metapopulation dynamics across species ranges and geographic variation in dispersal traits. To address this gap, we propose several future research directions that would advance our understanding of the interplay between spatiotemporal variability and dispersal trait variation in shaping the dynamics of current and future species ranges.     

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.     

Spatial heterogeneity of mortality and temporal fluctuation in fertility promote coexistence but not vice versa: a random-community approach

Both spatial heterogeneity and temporal fluctuation of the environment are important mechanisms promoting species coexistence, but they work in different manners. We consider many pairs of species with randomly generated survivorship and fertility in the lottery model, and examine how the variability in demographic processes affects the outcome of competition. The results are: [1] Coexistence is easier if habitat difference in mortality is greater, or if year-to-year variation in reproductive rate is larger. But neither habitat-difference in fertility nor temporal variation in mortality promotes coexistence. [2] Mean fertility does not affect the outcome if CV remains constant. In contrast, enhanced mean mortality decreases the fraction of coexisting pairs if the environment fluctuates temporally. [3] We also investigate the effect of limited dispersal of propagules between habitats. Compared with the complete mixing case, the fraction of coexisting pairs is clearly enhanced if the spatial heterogeneity is the major source of environmental variation, but shows slight increase if the temporal fluctuation is dominant. We conclude that spatial heterogeneity is likely to work more effectively in promoting species coexistence than temporal fluctuation, especially when the species suffer relatively high mortality, and disperse their propagules in a limited spatial scale.     

Spatial scaling of consumer-resource interactions in advection-dominated systems

Ecologists studying consumer-resource interactions in advection-dominated systems such as streams and rivers frequently seek to link the results of small-scale experiments with larger-scale patterns of distribution and abundance. Accomplishing this goal requires determining the characteristic scale, termed the response length, at which there is a shift from local dynamics dominated by advective dispersal to larger-scale dynamics dominated by births and deaths. Here, we model the dynamics of consumer-resource systems in a spatially variable, advective environment and show how consumer-resource interactions alter the response length relative to its single-species value. For one case involving a grazer that emigrates in response to high predator density, we quantify the changes using published data from small-scale experiments on aquatic invertebrates. Using Fourier analysis, we describe the responses of advection-dominated consumer-resource systems to spatially extended environmental variability in a way that involves explicit consideration of the response length. The patterns we derive for different consumer-resource systems exhibit important similarities in how component populations respond to spatial environmental variability affecting dispersal as opposed to demographic parameters.     

Evolution and adaptation in Pseudomonas aeruginosa biofilms driven by mismatch repair system-deficient mutators

Pseudomonas aeruginosa is an important opportunistic pathogen causing chronic airway infections, especially in cystic fibrosis (CF) patients. The majority of the CF patients acquire P. aeruginosa during early childhood, and most of them develop chronic infections resulting in severe lung disease, which are rarely eradicated despite intensive antibiotic therapy. Current knowledge indicates that three major adaptive strategies, biofilm development, phenotypic diversification, and mutator phenotypes [driven by a defective mismatch repair system (MRS)], play important roles in P. aeruginosa chronic infections, but the relationship between these strategies is still poorly understood. We have used the flow-cell biofilm model system to investigate the impact of the mutS associated mutator phenotype on development, dynamics, diversification and adaptation of P. aeruginosa biofilms. Through competition experiments we demonstrate for the first time that P. aeruginosa MRS-deficient mutators had enhanced adaptability over wild-type strains when grown in structured biofilms but not as planktonic cells. This advantage was associated with enhanced micro-colony development and increased rates of phenotypic diversification, evidenced by biofilm architecture features and by a wider range and proportion of morphotypic colony variants, respectively. Additionally, morphotypic variants generated in mutator biofilms showed increased competitiveness, providing further evidence for mutator-driven adaptive evolution in the biofilm mode of growth. This work helps to understand the basis for the specific high proportion and role of mutators in chronic infections, where P. aeruginosa develops in biofilm communities.