Hypermutability impedes cooperation in pathogenic bacteria

When the supply of beneficial mutations limits adaptation, bacterial mutator alleles can reach high frequencies by hitchhiking with advantageous mutations. However, when populations are well adapted to their environments, the increased rate of deleterious mutations makes hypermutability selectively disadvantageous. Here, we consider a further cost of hypermutability: its potential to break down cooperation (group-beneficial behavior that is costly to the individual). This probably occurs for three reasons. First, an increased rate at which 'cheating' genotypes are generated; second, an increased probability of producing efficient cheats; and third, a decrease in relatedness (not addressed in the present study). We used Pseudomonas aeruginosa's production of extracellular iron-scavenging molecules, siderophores, to determine if cheating evolved more readily in mutator populations. Siderophore production is costly to individual bacteria but benefits all nearby cells. Siderophore-deficient cheats therefore have a selective advantage within populations. We observed the de novo evolution and subsequent increase in frequency of siderophore cheats within both wild-type and mutator populations for 200 generations. Cheats appeared and increased in frequency more rapidly in mutator populations. The presence of cheats was costly to the group, as shown by a negative correlation between cheat frequency and population density.     

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.     

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.     

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.     

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.     

Detection of mutator subpopulations in Salmonella typhimurium LT2 by reversion of his alleles

Defects in the methyl-directed mismatch repair lead to both the hypermutability phenotype and removal of a barrier to genetic exchange between species. Mutator bacteria carrying such defects occur frequently among bacterial pathogens, suggesting that subpopulations of mutators are contained within pathogen clones and give rise to the genetic variants that are acted upon by selective forces to allow survival or successful infection. We report here on the detection of the mutator subpopulation in Salmonella typhimurium and determination of its frequency in laboratory cultures. The analysis involved screening for mutators among revertants of S. typhimurium histidine auxotrophs selected for the His⁺ phenotype, since the frequency of mutators is expected to be increased in the selected mutant population they helped to spawn. The increases in spontaneous reversion of histidine mutations were first measured in isogenic strains carrying mismatch repair-defective mutH, mutL, mutS, or uvrD alleles, relative to their mismatch repair-proficient counterparts. Screening for the mutator phenotype in nearly 12,000 revertants of repair-proficient strains carrying his mutations highly stimulated for reversion in mutator backgrounds, the base substitution in hisG428 and frameshift in hisC3076, yielded five mutator strains (0.04%). the his⁺ reversion mutations contained within the newly-arisen mutator strains were characteristic of the predominant nucleotide changes expected in such mutators, as assessed by comparison with the spectra for reversion events in wild-type and mismatch correction-defective backgrounds. The results show that subpopulations of mutators, residing in normal populations at a finite frequency, can be culled from the culture by strong selection for a required phenotype. We calculate that the frequency of mutators in the unselected population of S. typhimurium is 1–4×10⁻⁶, an incidence of 10-fold lower than that expected based on studies of laboratory cultures of Escherichia coli.     

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.     

Wider access to genotypic space facilitates loss of cooperation in a bacterial mutator

Understanding the ecological, evolutionary and genetic factors that affect the expression of cooperative behaviours is a topic of wide biological significance. On a practical level, this field of research is useful because many pathogenic microbes rely on the cooperative production of public goods (such as nutrient scavenging molecules, toxins and biofilm matrix components) in order to exploit their hosts. Understanding the evolutionary dynamics of cooperation is particularly relevant when considering long-term, chronic infections where there is significant potential for intra-host evolution. The impact of responses to non-social selection pressures on social evolution is arguably an under-examined area. In this paper, we consider how the evolution of a non-social trait--hypermutability--affects the cooperative production of iron-scavenging siderophores by the opportunistic human pathogen Pseudomonas aeruginosa. We confirm an earlier prediction that hypermutability accelerates the breakdown of cooperation due to increased sampling of genotypic space, allowing mutator lineages to generate non-cooperative genotypes with the ability to persist at high frequency and dominate populations. This may represent a novel cost of hypermutability.     

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.     

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 bacterial mutation rates under simultaneous selection by interspecific and social parasitism

Many bacterial populations harbour substantial numbers of hypermutable bacteria, in spite of hypermutation being associated with deleterious mutations. One reason for the persistence of hypermutators is the provision of novel mutations, enabling rapid adaptation to continually changing environments, for example coevolving virulent parasites. However, hypermutation also increases the rate at which intraspecific parasites (social cheats) are generated. Interspecific and intraspecific parasitism are therefore likely to impose conflicting selection pressure on mutation rate. Here, we combine theory and experiments to investigate how simultaneous selection from inter- and intraspecific parasitism affects the evolution of bacterial mutation rates in the plant-colonizing bacterium Pseudomonas fluorescens. Both our theoretical and experimental results suggest that phage presence increases and selection for public goods cooperation (the production of iron-scavenging siderophores) decreases selection for mutator bacteria. Moreover, phages imposed a much greater growth cost than social cheating, and when both selection pressures were imposed simultaneously, selection for cooperation did not affect mutation rate evolution. Given the ubiquity of infectious phages in the natural environment and clinical infections, our results suggest that phages are likely to be more important than social interactions in determining mutation rate evolution.     

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 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.     

Negotiation and appeasement can be more effective drivers of sociality than kin selection

Two alternative frameworks explain the evolution of cooperation in the face of conflicting interests. Conflicts can be alleviated by kinship, the alignment of interests by virtue of shared genes, or by negotiation strategies, allowing mutually beneficial trading of services or commodities. Although negotiation often occurs in kin-structured populations, the interplay of kin- and negotiation-based mechanisms in the evolution of cooperation remains an unresolved issue. Inspired by the biology of a cooperatively breeding fish, we developed an individual-based simulation model to study the evolution of negotiation-based cooperation in relation to different levels of genetic relatedness. We show that the evolution of negotiation strategies leads to an equilibrium where subordinates appease dominants by conditional cooperation, resulting in high levels of help and low levels of aggression. This negotiation-based equilibrium can be reached both in the absence of relatedness and in a kin-structured population. However, when relatedness is high, evolution often ends up in an alternative equilibrium where subordinates help their kin unconditionally. The level of help at this kin-selected equilibrium is considerably lower than at the negotiation-based equilibrium, and it corresponds to a level reached when responsiveness is prevented from evolving in the simulations. A mathematical invasion analysis reveals that, quite generally, the alignment of payoffs due to the relatedness of interaction partners tends to impede selection for harsh but effective punishment of defectors. Hence kin structure will often hamper rather than facilitate the evolution of productive cooperation.     

A simple rule for the evolution of contingent cooperation in large groups

Humans cooperate in large groups of unrelated individuals, and many authors have argued that such cooperation is sustained by contingent reward and punishment. However, such sanctioning systems can also stabilize a wide range of behaviours, including mutually deleterious behaviours. Moreover, it is very likely that large-scale cooperation is derived in the human lineage. Thus, understanding the evolution of mutually beneficial cooperative behaviour requires knowledge of when strategies that support such behaviour can increase when rare. Here, we derive a simple formula that gives the relatedness necessary for contingent cooperation in n-person iterated games to increase when rare. This rule applies to a wide range of pay-off functions and assumes that the strategies supporting cooperation are based on the presence of a threshold fraction of cooperators. This rule suggests that modest levels of relatedness are sufficient for invasion by strategies that make cooperation contingent on previous cooperation by a small fraction of group members. In contrast, only high levels of relatedness allow the invasion by strategies that require near universal cooperation. In order to derive this formula, we introduce a novel methodology for studying evolution in group structured populations including local and global group-size regulation and fluctuations in group size.     

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.     

Proliferation of mutators in A cell population.

A Lac- strain of Escherichia coli that reverts by the addition of a G to a G-G-G-G-G-G sequence was used to study the proliferation of mutators in a bacterial culture. Selection for the Lac+ phenotype, which is greatly stimulated in mismatch repair-deficient strains, results in an increase in the percentage of mutators in the selected population from less than 1 per 100,000 cells to 1 per 200 cells. All the mutators detected were deficient in the mismatch repair system. Mutagenesis results in a similar increase in the percentage of mutators. Mutagenesis combined with a single selection can result in a population of more than 50% mutators when a sample of several thousand cells is grown out and selected. Mutagenesis combined with two or more successive selections can generate a population that is 100% mutator. These experiments are discussed in relation to ideas that an early step in carcinogenesis is the creation of a mutator phenotype.     

Coexistence and Within-Host Evolution of Diversified Lineages of Hypermutable Pseudomonas aeruginosa in Long-term Cystic Fibrosis Infections

The advent of high-throughput sequencing techniques has made it possible to follow the genomic evolution of pathogenic bacteria by comparing longitudinally collected bacteria sampled from human hosts. Such studies in the context of chronic airway infections by Pseudomonas aeruginosa in cystic fibrosis (CF) patients have indicated high bacterial population diversity. Such diversity may be driven by hypermutability resulting from DNA mismatch repair system (MRS) deficiency, a common trait evolved by P. aeruginosa strains in CF infections. No studies to date have utilized whole-genome sequencing to investigate within-host population diversity or long-term evolution of mutators in CF airways. We sequenced the genomes of 13 and 14 isolates of P. aeruginosa mutator populations from an Argentinian and a Danish CF patient, respectively. Our collection of isolates spanned 6 and 20 years of patient infection history, respectively. We sequenced 11 isolates from a single sample from each patient to allow in-depth analysis of population diversity. Each patient was infected by clonal populations of bacteria that were dominated by mutators. The in vivo mutation rate of the populations was ∼100 SNPs/year–∼40-fold higher than rates in normo-mutable populations. Comparison of the genomes of 11 isolates from the same sample showed extensive within-patient genomic diversification; the populations were composed of different sub-lineages that had coexisted for many years since the initial colonization of the patient. Analysis of the mutations identified genes that underwent convergent evolution across lineages and sub-lineages, suggesting that the genes were targeted by mutation to optimize pathogenic fitness. Parallel evolution was observed in reduction of overall catabolic capacity of the populations. These findings are useful for understanding the evolution of pathogen populations and identifying new targets for control of chronic infections.     

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.     

Direct selection for mutators in Escherichia coli

We have constructed strains that allow a direct selection for mutators of Escherichia coli on a single plate medium. The plate selection is based on using two different markers whose reversion is enhanced by a given mutator. Plates containing limiting amounts of each respective nutrient allow the growth of ghost colonies or microcolonies that give rise to full-size colonies only if a reversion event occurs. Because two successive mutational events are required, mutator cells are favored to generate full-size colonies. Reversion of a third marker allows direct visualization of the mutator phenotype by the large number of blue papillae in the full-size colonies. We also describe plate selections involving three successive nutrient markers followed by a fourth papillation step. Different frameshift or base substitution mutations are used to select for mismatch-repair-defective strains (mutHLS and uvrD). We can detect and monitor mutator cells arising spontaneously, at frequencies lower than 10(-5) in the population. Also, we can measure a mutator cascade, in which one type of mutator (mutT) generates a second mutator (mutHLS) that then allows stepwise frameshift mutations. We discuss the relevance of mutators arising on a single medium as a result of cells overcoming successive growth barriers to the development and progression of cancerous tumors, some of which are mutator cell lines.