Experimental analysis of molecular events during mutational periodic selections in bacterial evolution

A fundamental feature of bacterial evolution is a succession of adaptive mutational sweeps when fitter mutants take over a population. To understand the processes involved in mutational successions, Escherichia coli continuous cultures were analyzed for changes at two loci where mutations provide strong transport advantages to fitness under steady-state glucose limitation. Three separate sweeps, observed as classic periodic selection events causing a change in the frequency of neutral mutations (in fhuA causing phage T5 resistance), were identified with changes at particular loci. Two of the sweeps were associated with a reduction in the frequency of neutral mutations and the concurrent appearance of at least 13 alleles at the mgl or mlc loci, respectively. These mgl and mlc polymorphisms were of many mutational types, so were not the result of a mutator or directed mutation event. The third sweep observed was altogether distinct and involved hitchhiking between T5 resistance and advantageous mgl mutations. Moreover, the hitchhiking event coincided with an increase in mutation rates, due to the transient appearance of a strong mutator in the population. The spectrum of mgl mutations among mutator isolates was distinct and due to mutS. The mutator-associated periodic selection also resulted in mgl and fhuA polymorphism in the sweeping population. These examples of periodic selections maintained significant genotypic diversity even in a rapidly evolving culture, with no individual "winner clone" or genotype purging the population.     

Bypass of genetic constraints during mutator evolution to antibiotic resistance

Genetic constraints can block many mutational pathways to optimal genotypes in real fitness landscapes, yet the extent to which this can limit evolution remains to be determined. Interestingly, mutator bacteria elevate only specific types of mutations, and therefore could be very sensitive to genetic constraints. Testing this possibility is not only clinically relevant, but can also inform about the general impact of genetic constraints in adaptation. Here, we evolved 576 populations of two mutator and one wild-type Escherichia coli to doubling concentrations of the antibiotic cefotaxime. All strains carried TEM-1, a β-lactamase enzyme well known by its low availability of mutational pathways. Crucially, one of the mutators does not elevate any of the relevant first-step mutations known to improve cefatoximase activity. Despite this, both mutators displayed a similar ability to evolve more than 1000-fold resistance. Initial adaptation proceeded in parallel through general multi-drug resistance mechanisms. High-level resistance, in contrast, was achieved through divergent paths; with the a priori inferior mutator exploiting alternative mutational pathways in PBP3, the target of the antibiotic. These results have implications for mutator management in clinical infections and, more generally, illustrate that limits to natural selection in real organisms are alleviated by the existence of multiple loci contributing to fitness.     

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.     

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.     

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.     

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.     

Mutational neighbourhood and mutation supply rate constrain adaptation in Pseudomonas aeruginosa

Understanding adaptation by natural selection requires understanding the genetic factors that determine which beneficial mutations are available for selection. Here, using experimental evolution of rifampicin-resistant Pseudomonas aeruginosa, we show that different genotypes vary in their capacity for adaptation to the cost of antibiotic resistance. We then use sequence data to show that the beneficial mutations associated with fitness recovery were specific to particular genetic backgrounds, suggesting that genotypes had access to different sets of beneficial mutations. When we manipulated the supply rate of beneficial mutations, by altering effective population size during evolution, we found that it constrained adaptation in some selection lines by restricting access to rare beneficial mutations, but that the effect varied among the genotypes in our experiment. These results suggest that mutational neighbourhood varies even among genotypes that differ by a single amino acid change, and this determines their capacity for adaptation as well as the influence of population biology processes that alter mutation supply rate.     

Role of hypermutability on bacterial fitness and emergence of resistance in experimental osteomyelitis due to Staphylococcus aureus

This study was designed to investigate the role of hypermutability of Staphylococcus aureus on bacterial fitness and antibiotic resistance in a model of chronic bone infection. An isogenic pair of strains, S. aureus RN4220 and its mutator counterpart inactivated in the mutL gene were used in a rat model of osteomyelitis of the tibia. The effect of the mutator phenotype in the emergence of antibiotic resistance was assessed in rats infected by each strain separately and treated with rifampicin for 5 days. No difference between the two strains was found in bacterial growth in vitro and in bacterial survival in the animal model, indicating no fitness defect in the mutator strain. In competition studies performed in rats coinfected with the two strains at a same ratio and sacrificed at different times from day 3 to day 42 postinoculation, the mutator strain was clearly disadvantaged because it was found in all rats and at all study times at lower counts (P<0.05 from day 14 to day 42). Two of the 16 rats infected by the mutator strain and none of the 14 rats infected by the wild-type strain had acquired rifampicin-resistant mutants (P=0.4). Data suggest that the S. aureus mutator phenotype was associated with a decreased bacterial fitness in vivo and did not confer significant advantage in the acquisition of antibiotic resistance in a model of chronic bone infection.     

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.     

Genetic Architecture of Intrinsic Antibiotic Susceptibility

Background

Antibiotic exposure rapidly selects for more resistant bacterial strains, and both a drug''s chemical structure and a bacterium''s cellular network affect the types of mutations acquired.

Methodology/Principal Findings

To better characterize the genetic determinants of antibiotic susceptibility, we exposed a transposon-mutagenized library of Escherichia coli to each of 17 antibiotics that encompass a wide range of drug classes and mechanisms of action. Propagating the library for multiple generations with drug concentrations that moderately inhibited the growth of the isogenic parental strain caused the abundance of strains with even minor fitness advantages or disadvantages to change measurably and reproducibly. Using a microarray-based genetic footprinting strategy, we then determined the quantitative contribution of each gene to E. coli''s intrinsic antibiotic susceptibility. We found both loci whose removal increased general antibiotic tolerance as well as pathways whose down-regulation increased tolerance to specific drugs and drug classes. The beneficial mutations identified span multiple pathways, and we identified pairs of mutations that individually provide only minor decreases in antibiotic susceptibility but that combine to provide higher tolerance.

Conclusions/Significance

Our results illustrate that a wide-range of mutations can modulate the activity of many cellular resistance processes and demonstrate that E. coli has a large mutational target size for increasing antibiotic tolerance. Furthermore, the work suggests that clinical levels of antibiotic resistance might develop through the sequential accumulation of chromosomal mutations of small individual effect.     

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.     

Genotype‐by‐environment interactions due to antibiotic resistance and adaptation in Escherichia coli

Mutations that are beneficial in one environment can have different fitness effects in other environments. In the context of antibiotic resistance, the resulting genotype‐by‐environment interactions potentially make selection on resistance unpredictable in heterogeneous environments. Furthermore, resistant bacteria frequently fix additional mutations during evolution in the absence of antibiotics. How do these two types of mutations interact to determine the bacterial phenotype across different environments? To address this, I used Escherichia coli as a model system, measuring the effects of nine different rifampicin resistance mutations on bacterial growth in 31 antibiotic‐free environments. I did this both before and after approximately 200 generations of experimental evolution in antibiotic‐free conditions (LB medium), and did the same for the antibiotic‐sensitive wild type after adaptation to the same environment. The following results were observed: (i) bacteria with and without costly resistance mutations adapted to experimental conditions and reached similar levels of competitive fitness; (ii) rifampicin resistance mutations and adaptation to LB both indirectly altered growth in other environments; and (iii) resistant‐evolved genotypes were more phenotypically different from the ancestor and from each other than resistant‐nonevolved and sensitive‐evolved genotypes. This suggests genotype‐by‐environment interactions generated by antibiotic resistance mutations, observed previously in short‐term experiments, are more pronounced after adaptation to other types of environmental variation, making it difficult to predict long‐term selection on resistance mutations from fitness effects in a single environment.     

A trade-off between oxidative stress resistance and DNA repair plays a role in the evolution of elevated mutation rates in bacteria

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

Effects of antibiotic resistance alleles on bacterial evolutionary responses to viral parasites

Antibiotic resistance has wide-ranging effects on bacterial phenotypes and evolution. However, the influence of antibiotic resistance on bacterial responses to parasitic viruses remains unclear, despite the ubiquity of such viruses in nature and current interest in therapeutic applications. We experimentally investigated this by exposing various Escherichia coli genotypes, including eight antibiotic-resistant genotypes and a mutator, to different viruses (lytic bacteriophages). Across 960 populations, we measured changes in population density and sensitivity to viruses, and tested whether variation among bacterial genotypes was explained by their relative growth in the absence of parasites, or mutation rate towards phage resistance measured by fluctuation tests for each phage. We found that antibiotic resistance had relatively weak effects on adaptation to phages, although some antibiotic-resistance alleles impeded the evolution of resistance to phages via growth costs. By contrast, a mutator allele, often found in antibiotic-resistant lineages in pathogenic populations, had a relatively large positive effect on phage-resistance evolution and population density under parasitism. This suggests costs of antibiotic resistance may modify the outcome of phage therapy against pathogenic populations previously exposed to antibiotics, but the effects of any co-occurring mutator alleles are likely to be stronger.     

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.     

The cost of multiple drug resistance in Pseudomonas aeruginosa

The spread of bacterial antibiotic resistance mutations is thought to be constrained by their pleiotropic fitness costs. Here we investigate the fitness costs of resistance in the context of the evolution of multiple drug resistance (MDR), by measuring the cost of acquiring streptomycin resistance mutations (StrepR) in independent strains of the bacterium Pseudomonas aeruginosa carrying different rifampicin resistance (RifR) mutations. In the absence of antibiotics, StrepR mutations are associated with similar fitness costs in different RifR genetic backgrounds. The cost of StrepR mutations is greater in a rifampicin‐sensitive (RifS) background, directly demonstrating antagonistic epistasis between resistance mutations. In the presence of rifampicin, StrepR mutations have contrasting effects in different RifR backgrounds: StrepR mutations have no detectable costs in some RifR backgrounds and massive fitness costs in others. Our results clearly demonstrate the importance of epistasis and genotype‐by‐environment interactions for the evolution of MDR.     

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.     

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.     

Mutations in dnaA and a cryptic interaction site increase drug resistance in Mycobacterium tuberculosis

Genomic dissection of antibiotic resistance in bacterial pathogens has largely focused on genetic changes conferring growth above a single critical concentration of drug. However, reduced susceptibility to antibiotics—even below this breakpoint—is associated with poor treatment outcomes in the clinic, including in tuberculosis. Clinical strains of Mycobacterium tuberculosis exhibit extensive quantitative variation in antibiotic susceptibility but the genetic basis behind this spectrum of drug susceptibility remains ill-defined. Through a genome wide association study, we show that non-synonymous mutations in dnaA, which encodes an essential and highly conserved regulator of DNA replication, are associated with drug resistance in clinical M. tuberculosis strains. We demonstrate that these dnaA mutations specifically enhance M. tuberculosis survival during isoniazid treatment via reduced expression of katG, the activator of isoniazid. To identify DnaA interactors relevant to this phenotype, we perform the first genome-wide biochemical mapping of DnaA binding sites in mycobacteria which reveals a DnaA interaction site that is the target of recurrent mutation in clinical strains. Reconstructing clinically prevalent mutations in this DnaA interaction site reproduces the phenotypes of dnaA mutants, suggesting that clinical strains of M. tuberculosis have evolved mutations in a previously uncharacterized DnaA pathway that quantitatively increases resistance to the key first-line antibiotic isoniazid. Discovering genetic mechanisms that reduce drug susceptibility and support the evolution of high-level drug resistance will guide development of biomarkers capable of prospectively identifying patients at risk of treatment failure in the clinic.     

Competition between high- and higher-mutating strains of Escherichia coli

Experimental studies have shown that a mutator allele can readily hitchhike to fixation with beneficial mutations in an asexual population having a low, wild-type mutation rate. Here, we show that a genotype bearing two mutator alleles can supplant a population already fixed for one mutator allele. Our results provide experimental support for recent theory predicting that mutator alleles will tend to accumulate in asexual populations by hitchhiking with beneficial mutations, causing an ever-higher genomic mutation rate.