Mutational spectrum drives the rise of mutator bacteria

Quantitative analysis of relationship between mutation rate and speed of adaptation under antibiotic exposure in Escherichia coli

Mutations are the ultimate source of biological evolution that creates genetic variation in populations. Mutations can create new advantageous traits but can also potentially interfere with pre-existing organismal functions. Therefore, organisms may have evolved mutation rates to appropriate levels to maintain or improve their fitness. In this study, we aimed to experimentally quantify the relationship between the mutation rate and evolution of antibiotic resistance. We conducted an evolution experiment using 12 Escherichia coli mutator strains with increased mutation rates and five antibiotics. Our results demonstrated that the rate of adaptation generally increased with higher mutation rates, except in a single mutator strain with the highest mutation rate, which exhibited a significant decline in evolutionary speed. To further elucidate these findings, we developed a simple population dynamics model that successfully recapitulated the observed dependence of adaptation speed on mutation rate. These findings provide important insights into the evolution of mutation rate accompanied by the evolution.

Experimental evolution of evolvability

Evolvability-the capacity to generate adaptive variation-is a trait that can itself evolve through natural selection. However, the idea that mutation can become biased toward adaptive outcomes remains controversial. In this work, we report the evolution of enhanced evolvability through localized hypermutation in experimental populations of bacteria. The evolved mechanism is analogous to the mutation-prone sequences of contingency loci observed in pathogenic bacteria. Central to this outcome was a lineage-level selection process, where success depended on the capacity to evolve between two phenotypic states. Subsequent evolution showed that the hypermutable locus is itself evolvable with respect to alterations in the frequency of environmental change. Lineages with localized hypermutability were more likely to acquire additional adaptive mutations, revealing an unanticipated benefit.

Evolution of evolvability in rapidly adapting populations

Mutations can alter the short-term fitness of an organism, as well as the rates and benefits of future mutations. While numerous examples of these evolvability modifiers have been observed in rapidly adapting microbial populations, existing theory struggles to predict when they will be favoured by natural selection. Here we develop a mathematical framework for predicting the fates of genetic variants that modify the rates and benefits of future mutations in linked genomic regions. We derive analytical expressions showing how the fixation probabilities of these variants depend on the size of the population and the diversity of competing mutations. We find that competition between linked mutations can dramatically enhance selection for modifiers that increase the benefits of future mutations, even when they impose a strong direct cost on fitness. However, we also find that modest direct benefits can be sufficient to drive evolutionary dead ends to fixation. Our results suggest that subtle differences in evolvability could play an important role in shaping the long-term success of genetic variants in rapidly evolving microbial populations.

Direct and indirect selection in a proofreading polymerase

The traits that affect evolvability are subject to indirect selection, as these traits affect the course of evolution over many generations rather than the direct replicative fitness of an individual. However, the evolution of evolvability-determining traits is often difficult to study because putative evolvability alleles often have confounding direct fitness effects of unknown origin and size. Here, we study theoretically and experimentally the evolution of mutation rates in proofreading polymerases with orthogonal control of direct and indirect selection. Mutagenic DNA polymerases enjoy a long-time fitness advantage by enhancing the rate of acquiring beneficial mutations. However, this is offset by a short-time fitness penalty, which we trace to a counterintuitive trade-off between mutation rates and activity in proofreading polymerases. Since these fitness effects act on different timescales, no one number characterizes the fitness of a mutator allele. We find unusual dynamic features in the resulting evolutionary dynamics, such as kinetic exclusion, selection by dynamic environments, and Rock-Paper-Scissors dynamics in the absence of ecology. Our work has implications for the evolution of mutation rates and more broadly, evolution in the context of an anti-correlation between mutation rates and short term fitness.

Shifts in Mutation Bias Promote Mutators by Altering the Distribution of Fitness Effects

AbstractRecent experimental evidence demonstrates that shifts in mutational biases-for example, increases in transversion frequency-can change the distribution of fitness effects of mutations (DFE). In particular, reducing or reversing a prevailing bias can increase the probability that a de novo mutation is beneficial. It has also been shown that mutator bacteria are more likely to emerge if the beneficial mutations they generate have a larger effect size than observed in the wild type. Here, we connect these two results, demonstrating that mutator strains that reduce or reverse a prevailing bias have a positively shifted DFE, which in turn can dramatically increase their emergence probability. Since changes in mutation rate and bias are often coupled through the gain and loss of DNA repair enzymes, our results predict that the invasion of mutator strains will be facilitated by shifts in mutation bias that offer improved access to previously undersampled beneficial mutations.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Hypermutator emergence in experimental Escherichia coli populations is stress-type dependent

Genotypes exhibiting an increased mutation rate, called hypermutators, can propagate in microbial populations because they can have an advantage due to the higher supply of beneficial mutations needed for adaptation. Although this is a frequently observed phenomenon in natural and laboratory populations, little is known about the influence of parameters such as the degree of maladaptation, stress intensity, and the genetic architecture for adaptation on the emergence of hypermutators. To address this knowledge gap, we measured the emergence of hypermutators over ~1,000 generations in experimental Escherichia coli populations exposed to different levels of osmotic or antibiotic stress. Our stress types were chosen based on the assumption that the genetic architecture for adaptation differs between them. Indeed, we show that the size of the genetic basis for adaptation is larger for osmotic stress compared to antibiotic stress. During our experiment, we observed an increased emergence of hypermutators in populations exposed to osmotic stress but not in those exposed to antibiotic stress, indicating that hypermutator emergence rates are stress type dependent. These results support our hypothesis that hypermutator emergence is linked to the size of the genetic basis for adaptation. In addition, we identified other parameters that covaried with stress type (stress level and IS transposition rates) that might have contributed to an increased hypermutator provision and selection. Our results provide a first comparison of hypermutator emergence rates under varying stress conditions and point towards complex interactions of multiple stress-related factors on the evolution of mutation rates.

A Multimodal Approach towards Genomic Identification of Protein Inhibitors of Uracil-DNA Glycosylase

DNA-mimicking proteins encoded by viruses can modulate processes such as innate cellular immunity. An example is Ung-family uracil-DNA glycosylase inhibition, which prevents Ung-mediated degradation via the stoichiometric protein blockade of the Ung DNA-binding cleft. This is significant where uracil-DNA is a key determinant in the replication and distribution of virus genomes. Unrelated protein folds support a common physicochemical spatial strategy for Ung inhibition, characterised by pronounced sequence plasticity within the diverse fold families. That, and the fact that relatively few template sequences are biochemically verified to encode Ung inhibitor proteins, presents a barrier to the straightforward identification of Ung inhibitors in genomic sequences. In this study, distant homologs of known Ung inhibitors were characterised via structural biology and structure prediction methods. A recombinant cellular survival assay and in vitro biochemical assay were used to screen distant variants and mutants to further explore tolerated sequence plasticity in motifs supporting Ung inhibition. The resulting validated sequence repertoire defines an expanded set of heuristic sequence and biophysical signatures shared by known Ung inhibitor proteins. A computational search of genome database sequences and the results of recombinant tests of selected output sequences obtained are presented here.

Shifts in mutation spectra enhance access to beneficial mutations

Biased mutation spectra are pervasive, with wide variation in the magnitude of mutational biases that influence genome evolution and adaptation. How do such diverse biases evolve? Our experiments show that changing the mutation spectrum allows populations to sample previously undersampled mutational space, including beneficial mutations. The resulting shift in the distribution of fitness effects is advantageous: Beneficial mutation supply and beneficial pleiotropy both increase, while deleterious load reduces. More broadly, simulations indicate that reducing or reversing the direction of a long-term bias is always selectively favored. Such changes in mutation bias can occur easily via altered function of DNA repair genes. A phylogenetic analysis shows that these genes are repeatedly gained and lost in bacterial lineages, leading to frequent bias shifts in opposite directions. Thus, shifts in mutation spectra may evolve under selection and can directly alter the outcome of adaptive evolution by facilitating access to beneficial mutations.

Transient mutation bias increases the predictability of evolution on an empirical genotype-phenotype landscape

Predicting how a population will likely navigate a genotype-phenotype landscape requires consideration of selection in combination with mutation bias, which can skew the likelihood of following a particular trajectory. Strong and persistent directional selection can drive populations to ascend toward a peak. However, with a greater number of peaks and more routes to reach them, adaptation inevitably becomes less predictable. Transient mutation bias, which operates only on one mutational step, can influence landscape navigability by biasing the mutational trajectory early in the adaptive walk. This sets an evolving population upon a particular path, constraining the number of accessible routes and making certain peaks and routes more likely to be realized than others. In this work, we employ a model system to investigate whether such transient mutation bias can reliably and predictably place populations on a mutational trajectory to the strongest selective phenotype or usher populations to realize inferior phenotypic outcomes. For this we use motile mutants evolved from ancestrally non-motile variants of the microbe Pseudomonas fluorescens SBW25, of which one trajectory exhibits significant mutation bias. Using this system, we elucidate an empirical genotype-phenotype landscape, where the hill-climbing process represents increasing strength of the motility phenotype, to reveal that transient mutation bias can facilitate rapid and predictable ascension to the strongest observed phenotype in place of equivalent and inferior trajectories. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Evolution of the SARS-CoV-2 Mutational Spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here, we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single-nucleotide mutations at 4-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly 2-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al. (2022. Mutational spectra distinguish SARS-CoV-2 replication niches. bioRxiv, doi:10.1101/2022.09.27.509649). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors, although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution and suggests that human SARS-CoV-2 may be trending toward a mutation spectrum more similar to that of other animal sarbecoviruses.

Clostridium beijerinckii strain degeneration is driven by the loss of Spo0A activity

Solventogenic clostridia represent a diverse group of anaerobic, spore-forming bacteria capable of producing acetone, butanol and ethanol through their unique biphasic metabolism. An intrinsic problem with these organisms however is their tendency to degenerate when repeatedly subcultured or when grown continuously. This phenomenon sees cells lose their ability to produce solvents and spores, posing a significant problem for industrial applications. To investigate the mechanistic and evolutionary basis of degeneration we combined comparative genomics, ultra-deep sequencing, and concepts of sociomicrobiology using Clostridium beijerinckii NCIMB 8052 as our model organism. These approaches revealed spo0A, the master regulator gene involved in spore and solvent formation, to be key to the degeneration process in this strain. Comparative genomics of 71 degenerate variants revealed four distinct hotspot regions that contained considerably more mutations than the rest of the genome. These included spo0A as well as genes suspected to regulate its expression and activity. Ultra-deep sequencing of populations during the subculturing process showed transient increases in mutations we believe linked to the spo0A network, however, these were ultimately dominated by mutations in the master regulator itself. Through frequency-dependent fitness assays, we found that spo0A mutants gained a fitness advantage, relative to the wild type, presumably allowing for propagation throughout the culture. Combined, our data provides new insights into the phenomenon of clostridial strain degeneration and the C. beijerinckii NCIMB 8052 solvent and spore regulation network.

Evolution of the SARS-CoV-2 mutational spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single nucleotide mutations at four-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly two-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al (2022). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors-although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution, and suggests that human SARS-CoV-2 may be trending towards a mutation spectrum more similar to that of other animal sarbecoviruses.

Hypermutator strains of Pseudomonas aeruginosa reveal novel pathways of resistance to combinations of cephalosporin antibiotics and beta-lactamase inhibitors

Hypermutation due to DNA mismatch repair (MMR) deficiencies can accelerate the development of antibiotic resistance in Pseudomonas aeruginosa. Whether hypermutators generate resistance through predominantly similar molecular mechanisms to wild-type (WT) strains is not fully understood. Here, we show that MMR-deficient P. aeruginosa can evolve resistance to important broad-spectrum cephalosporin/beta-lactamase inhibitor combination antibiotics through novel mechanisms not commonly observed in WT lineages. Using whole-genome sequencing (WGS) and transcriptional profiling of isolates that underwent in vitro adaptation to ceftazidime/avibactam (CZA), we characterized the detailed sequence of mutational and transcriptional changes underlying the development of resistance. Surprisingly, MMR-deficient lineages rapidly developed high-level resistance (>256 μg/mL) largely without corresponding fixed mutations or transcriptional changes in well-established resistance genes. Further investigation revealed that these isolates had paradoxically generated an early inactivating mutation in the mexB gene of the MexAB-OprM efflux pump, a primary mediator of CZA resistance in P. aeruginosa, potentially driving an evolutionary search for alternative resistance mechanisms. In addition to alterations in a number of genes not known to be associated with resistance, 2 mutations were observed in the operon encoding the RND efflux pump MexVW. These mutations resulted in a 4- to 6-fold increase in resistance to ceftazidime, CZA, cefepime, and ceftolozane-tazobactam when engineered into a WT strain, demonstrating a potentially important and previously unappreciated mechanism of resistance to these antibiotics in P. aeruginosa. Our results suggest that MMR-deficient isolates may rapidly evolve novel resistance mechanisms, sometimes with complex dynamics that reflect gene inactivation that occurs with hypermutation. The apparent ease with which hypermutators may switch to alternative resistance mechanisms for which antibiotics have not been developed may carry important clinical implications.

Co-option of stress mechanisms in the origin of evolutionary novelties

It is widely accepted that stressful conditions can facilitate evolutionary change. The mechanisms elucidated thus far accomplish this with a generic increase in heritable variation that facilitates more rapid adaptive evolution, often via plastic modifications of existing characters. Through scrutiny of different meanings of stress in biological research, and an explicit recognition that stressors must be characterized relative to their effect on capacities for maintaining functional integrity, we distinguish between: (1) previously identified stress‐responsive mechanisms that facilitate evolution by maintaining an adaptive fit with the environment, and (2) the co‐option of stress‐responsive mechanisms that are specific to stressors leading to the origin of novelties via compensation. Unlike standard accounts of gene co‐option that identify component sources of evolutionary change, our model documents the cost‐benefit trade‐offs and thereby explains how one mechanism—an immediate response to acute stress—is transformed evolutionarily into another—routine protection from recurring stressors. We illustrate our argument with examples from cell type origination as well as processes and structures at higher levels of organization. These examples suggest a general principle of evolutionary origination based on the capacity to switch between regulatory states related to reproduction and proliferation versus survival and differentiation.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Bacterial phenotypic heterogeneity in DNA repair and mutagenesis

Genetically identical cells frequently exhibit striking heterogeneity in various phenotypic traits such as their morphology, growth rate, or gene expression. Such non-genetic diversity can help clonal bacterial populations overcome transient environmental challenges without compromising genome stability, while genetic change is required for long-term heritable adaptation. At the heart of the balance between genome stability and plasticity are the DNA repair pathways that shield DNA from lesions and reverse errors arising from the imperfect DNA replication machinery. In principle, phenotypic heterogeneity in the expression and activity of DNA repair pathways can modulate mutation rates in single cells and thus be a source of heritable genetic diversity, effectively reversing the genotype-to-phenotype dogma. Long-standing evidence for mutation rate heterogeneity comes from genetics experiments on cell populations, which are now complemented by direct measurements on individual living cells. These measurements are increasingly performed using fluorescence microscopy with a temporal and spatial resolution that enables localising, tracking, and counting proteins with single-molecule sensitivity. In this review, we discuss which molecular processes lead to phenotypic heterogeneity in DNA repair and consider the potential consequences on genome stability and dynamics in bacteria. We further inspect these concepts in the context of DNA damage and mutation induced by antibiotics.

Point mutations in topoisomerase I alter the mutation spectrum in E. coli and impact the emergence of drug resistance genotypes

Identifying the molecular mechanisms that give rise to genetic variation is essential for the understanding of evolutionary processes. Previously, we have used adaptive laboratory evolution to enable biomass synthesis from CO2 in Escherichia coli. Genetic analysis of adapted clones from two independently evolving populations revealed distinct enrichment for insertion and deletion mutational events. Here, we follow these observations to show that mutations in the gene encoding for DNA topoisomerase I (topA) give rise to mutator phenotypes with characteristic mutational spectra. Using genetic assays and mutation accumulation lines, we find that point mutations in topA increase the rate of sequence deletion and duplication events. Interestingly, we observe that a single residue substitution (R168C) results in a high rate of head-to-tail (tandem) short sequence duplications, which are independent of existing sequence repeats. Finally, we show that the unique mutation spectrum of topA mutants enhances the emergence of antibiotic resistance in comparison to mismatch-repair (mutS) mutators, and leads to new resistance genotypes. Our findings highlight a potential link between the catalytic activity of topoisomerases and the fundamental question regarding the emergence of de novo tandem repeats, which are known modulators of bacterial evolution.

Prevalence of hypermutator isolates of Achromobacter spp. from cystic fibrosis patients

Bacteria colonising the lungs of cystic fibrosis (CF) patients encounter high selective pressures. Hypermutation facilitates adaptation to fluctuating environments, and hypermutator strains are frequently isolated from CF patients. We investigated the prevalence of hypermutator isolates of Achromobacter spp. among patients affiliated with the CF Centre in Aarhus, Denmark. By exposure to rifampicin, the mutation frequency was determined for 90 isolates of Achromobacter spp. cultured from 42 CF patients; 20 infections were categorised as chronic, 22 as intermittent. The genetic mechanisms of hypermutation were examined by comparing DNA repair gene sequences from hypermutator and normomutator isolates. Achromobacter spp. cultured from 11 patients were categorised as hypermutators, and this phenotype was exclusively associated with chronic infections. Isolates of the Danish epidemic strain (DES) of Achromobacter ruhlandii cultured from patients from both Danish CF centres showed elevated mutation frequencies. The hypermutator state of Achromobacter spp. was most commonly associated with nonsynonymous mutations in the DNA mismatch repair gene mutS; a single clone had developed a substitution in the S-adenosyl-L-methionine-dependent methyltransferase putatively involved in DNA repair mechanisms, but not previously linked to the hypermutator phenotype. Hypermutation is prevalent among clinical isolates of Achromobacter spp. and could be a key determinant for the extraordinary adaptation and persistence of DES.

Irregularities in genetic variation and mutation rates with environmental stresses

The appearance of new mutations is determined by the equilibrium between DNA error formation and repair. In bacteria like Escherichia coli, stresses are thought shift this balance towards increased mutagenesis. Recent findings, however, suggest a very uneven relationship between stress and mutations. Only a subset of stressful environments increase the net rate of mutation and different forms of nutritional stress (such as oxygen, carbon or phosphorus limitations) result in markedly different mutation rates after similar reductions in growth rate. Moreover, different stresses result in altered mutational spectra, with some increasing transposition and others increasing indel formation. Single-base substitution rates are lower with some stresses than in unstressed bacteria. Indeed, changes to the mix of mutations with stress are more widespread than a marked increase in net mutation rate. Much remains to be learned on how environments have unique mutational signatures and why some stresses are more mutagenic than others. Even beyond stress-induced genetic variation, the fundamental unresolved question in the stress-mutation relationship is the adaptive value of different types of mutations and mutation rates; is transposition, for example, more advantageous under anaerobic conditions? It remains to be investigated whether stress-specific genetic variation impacts on evolvability differentially in distinct environments.

Mutation bias and GC content shape antimutator invasions

Mutators represent a successful strategy in rapidly adapting asexual populations, but theory predicts their eventual extinction due to their unsustainably large deleterious load. While antimutator invasions have been documented experimentally, important discrepancies among studies remain currently unexplained. Here we show that a largely neglected factor, the mutational idiosyncrasy displayed by different mutators, can play a major role in this process. Analysing phylogenetically diverse bacteria, we find marked and systematic differences in the protein-disruptive effects of mutations caused by different mutators in species with different GC compositions. Computer simulations show that these differences can account for order-of-magnitude changes in antimutator fitness for a realistic range of parameters. Overall, our results suggest that antimutator dynamics may be highly dependent on the specific genetic, ecological and evolutionary history of a given population. This context-dependency further complicates our understanding of mutators in clinical settings, as well as their role in shaping bacterial genome size and composition.

Opposing effects of final population density and stress on Escherichia coli mutation rate

Evolution depends on mutations. For an individual genotype, the rate at which mutations arise is known to increase with various stressors (stress-induced mutagenesis-SIM) and decrease at high final population density (density-associated mutation-rate plasticity-DAMP). We hypothesised that these two forms of mutation-rate plasticity would have opposing effects across a nutrient gradient. Here we test this hypothesis, culturing Escherichia coli in increasingly rich media. We distinguish an increase in mutation rate with added nutrients through SIM (dependent on error-prone polymerases Pol IV and Pol V) and an opposing effect of DAMP (dependent on MutT, which removes oxidised G nucleotides). The combination of DAMP and SIM results in a mutation rate minimum at intermediate nutrient levels (which can support 7 × 108 cells ml-1). These findings demonstrate a strikingly close and nuanced relationship of ecological factors-stress and population density-with mutation, the fuel of all evolution.

Parallel Evolution of High-Level Aminoglycoside Resistance in Escherichia coli Under Low and High Mutation Supply Rates

Antibiotic resistance is a major concern in public health worldwide, thus there is much interest in characterizing the mutational pathways through which susceptible bacteria evolve resistance. Here we use experimental evolution to explore the mutational pathways toward aminoglycoside resistance, using gentamicin as a model, under low and high mutation supply rates. Our results show that both normo and hypermutable strains of Escherichia coli are able to develop resistance to drug dosages > 1,000-fold higher than the minimal inhibitory concentration for their ancestors. Interestingly, such level of resistance was often associated with changes in susceptibility to other antibiotics, most prominently with increased resistance to fosfomycin. Whole-genome sequencing revealed that all resistant derivatives presented diverse mutations in five common genetic elements: fhuA, fusA and the atpIBEFHAGDC, cyoABCDE, and potABCD operons. Despite the large number of mutations acquired, hypermutable strains did not pay, apparently, fitness cost. In contrast to recent studies, we found that the mutation supply rate mainly affected the speed (tempo) but not the pattern (mode) of evolution: both backgrounds acquired the mutations in the same order, although the hypermutator strain did it faster. This observation is compatible with the adaptive landscape for high-level gentamicin resistance being relatively smooth, with few local maxima; which might be a common feature among antibiotics for which resistance involves multiple loci.

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.

Mutation accumulation under UV radiation in Escherichia coli

Mutations are induced by not only intrinsic factors such as inherent molecular errors but also by extrinsic mutagenic factors such as UV radiation. Therefore, identifying the mutational properties for both factors is necessary to achieve a comprehensive understanding of evolutionary processes both in nature and in artificial situations. Although there have been extensive studies on intrinsic factors, the mutational profiles of extrinsic factors are poorly understood on a genomic scale. Here, we explored the mutation profiles of UV radiation, a ubiquitous mutagen, in Escherichia coli on the genomic scale. We performed an evolution experiment under periodic UV radiation for 28 days. The accumulation speed of the mutations was found to increase so that it exceeded that of a typical mutator strain with deficient mismatch repair processes. The huge contribution of the extrinsic factors to all mutations consequently increased the risk of the destruction of inherent error correction systems. The spectrum of the UV-induced mutations was broader than that of the spontaneous mutations in the mutator. The broad spectrum and high upper limit of the frequency of occurrence suggested ubiquitous roles for UV radiation in accelerating the evolutionary process.

The road not taken: Could stress-specific mutations lead to different evolutionary paths?

Organisms often encounter stressful conditions, some of which damage their DNA. In response, some organisms show a high expression of error-prone DNA repair machinery, causing a temporary increase in the genome-wide mutation rate. Although we now have a detailed map of the molecular mechanisms underlying such stress-induced mutagenesis (SIM), it has been hotly debated whether SIM alters evolutionary dynamics. Key to this controversy is our poor understanding about which stresses increase mutagenesis and their long-term consequences for adaptation. In a new study with Escherichia coli, Maharjan and Ferenci show that while only some nutritional stresses (phosphorous and carbon limitation) increase total mutation rates, each stress generates a unique spectrum of mutations. Their results suggest the potential for specific stresses to shape evolutionary dynamics and highlight the necessity for explicit tests of the long-term evolutionary impacts of SIM.

Adaptive tuning of mutation rates allows fast response to lethal stress in Escherichia coli

While specific mutations allow organisms to adapt to stressful environments, most changes in an organism's DNA negatively impact fitness. The mutation rate is therefore strictly regulated and often considered a slowly-evolving parameter. In contrast, we demonstrate an unexpected flexibility in cellular mutation rates as a response to changes in selective pressure. We show that hypermutation independently evolves when different Escherichia coli cultures adapt to high ethanol stress. Furthermore, hypermutator states are transitory and repeatedly alternate with decreases in mutation rate. Specifically, population mutation rates rise when cells experience higher stress and decline again once cells are adapted. Interestingly, we identified cellular mortality as the major force driving the quick evolution of mutation rates. Together, these findings show how organisms balance robustness and evolvability and help explain the prevalence of hypermutation in various settings, ranging from emergence of antibiotic resistance in microbes to cancer relapses upon chemotherapy.

Mutation tendency of mutator Plasmodium berghei with proofreading-deficient DNA polymerase δ

In this study, we investigated the mutation tendency of a mutator rodent malaria parasite, Plasmodium berghei, with proofreading-deficient DNA polymerase δ. Wild-type and mutator parasites were maintained in mice for over 24 weeks, and the genome-wide accumulated mutations were determined by high-throughput sequencing. The mutator P. berghei had a significant preference for C/G to A/T substitutions; thus, its genome had a trend towards a higher AT content. The mutation rate was influenced by the sequence context, and mutations were markedly elevated at TCT. Some genes mutated repeatedly in replicate passage lines. In particular, knockout mutations of the AP2-G gene were frequent, which conferred strong growth advantages on parasites during the blood stage but at the cost of losing the ability to form gametocytes. This is the first report to demonstrate a biased mutation tendency in malaria parasites, and its results help to promote our basic understanding of Plasmodium genetics.

Genome-Wide Biases in the Rate and Molecular Spectrum of Spontaneous Mutations in Vibrio cholerae and Vibrio fischeri

The vast diversity in nucleotide composition and architecture among bacterial genomes may be partly explained by inherent biases in the rates and spectra of spontaneous mutations. Bacterial genomes with multiple chromosomes are relatively unusual but some are relevant to human health, none more so than the causative agent of cholera, Vibrio cholerae Here, we present the genome-wide mutation spectra in wild-type and mismatch repair (MMR) defective backgrounds of two Vibrio species, the low-%GC squid symbiont V. fischeri and the pathogen V. cholerae, collected under conditions that greatly minimize the efficiency of natural selection. In apparent contrast to their high diversity in nature, both wild-type V. fischeri and V. cholerae have among the lowest rates for base-substitution mutations (bpsms) and insertion-deletion mutations (indels) that have been measured, below 10^-3/genome/generation. Vibrio fischeri and V. cholerae have distinct mutation spectra, but both are AT-biased and produce a surprising number of multi-nucleotide indels. Furthermore, the loss of a functional MMR system caused the mutation spectra of these species to converge, implying that the MMR system itself contributes to species-specific mutation patterns. Bpsm and indel rates varied among genome regions, but do not explain the more rapid evolutionary rates of genes on chromosome 2, which likely result from weaker purifying selection. More generally, the very low mutation rates of Vibrio species correlate inversely with their immense population sizes and suggest that selection may not only have maximized replication fidelity but also optimized other polygenic traits relative to the constraints of genetic drift.

Optimization of Polymyxin B in Combination with Doripenem To Combat Mutator Pseudomonas aeruginosa

Development of spontaneous mutations in Pseudomonas aeruginosa has been associated with antibiotic failure, leading to high rates of morbidity and mortality. Our objective was to evaluate the pharmacodynamics of polymyxin B combinations against rapidly evolving P. aeruginosa mutator strains and to characterize the time course of bacterial killing and resistance via mechanism-based mathematical models. Polymyxin B or doripenem alone and in combination were evaluated against six P. aeruginosa strains: wild-type PAO1, mismatch repair (MMR)-deficient (mutS and mutL) strains, and 7,8-dihydro-8-oxo-deoxyguanosine system (GO) base excision repair (BER)-deficient (mutM, mutT, and mutY) strains over 48 h. Pharmacodynamic modeling was performed using S-ADAPT and facilitated by SADAPT-TRAN. Mutator strains displayed higher mutation frequencies than the wild type (>600-fold). Exposure to monotherapy was followed by regrowth, even at high polymyxin B concentrations of up to 16 mg/liter. Polymyxin B and doripenem combinations displayed enhanced killing activity against all strains where complete eradication was achieved for polymyxin B concentrations of >4 mg/liter and doripenem concentrations of 8 mg/liter. Modeling suggested that the proportion of preexisting polymyxin B-resistant subpopulations influenced the pharmacodynamic profiles for each strain uniquely (fraction of resistance values are -8.81 log10 for the wild type, -4.71 for the mutS mutant, and -7.40 log10 for the mutM mutant). Our findings provide insight into the optimization of polymyxin B and doripenem combinations against P. aeruginosa mutator strains.

The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa

Background

Natural genetic variation ultimately arises from the process of mutation. Knowledge of how the raw material for evolution is produced is necessary for a full understanding of several fundamental evolutionary concepts. We performed a mutation accumulation experiment with wild-type and mismatch-repair deficient, mutator lines of the pathogenic bacterium Pseudomonas aeruginosa, and used whole-genome sequencing to reveal the genome-wide rate, spectrum, distribution, leading/lagging bias, and context-dependency of spontaneous mutations.

Results

Wild-type base-pair mutation and indel rates were ~10⁻¹⁰ and ~10⁻¹¹ per nucleotide per generation, respectively, and deficiencies in the mismatch-repair system caused rates to increase by over two orders of magnitude. A universal bias towards AT was observed in wild-type lines, but was reversed in mutator lines to a bias towards GC. Biases for which types of mutations occurred during replication of the leading versus lagging strand were detected reciprocally in both replichores. The distribution of mutations along the chromosome was non-random, with peaks near the terminus of replication and at positions intermediate to the replication origin and terminus. A similar distribution bias was observed along the chromosome in natural populations of P. aeruginosa. Site-specific mutation rates were higher when the focal nucleotide was immediately flanked by C:G pairings.

Conclusions

Whole-genome sequencing of mutation accumulation lines allowed the comprehensive identification of mutations and revealed what factors of molecular and genomic architecture affect the mutational process. Our study provides a more complete view of how several mechanisms of mutation, mutation repair, and bias act simultaneously to produce the raw material for evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2244-3) contains supplementary material, which is available to authorized users.

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.

Fitness is strongly influenced by rare mutations of large effect in a microbial mutation accumulation experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.

The advantage of arriving first: characteristic times in finite size populations of error-prone replicators

We study the evolution of a finite size population formed by mutationally isolated lineages of error-prone replicators in a two-peak fitness landscape. Computer simulations are performed to gain a stochastic description of the system dynamics. More specifically, for different population sizes, we compute the probability of each lineage being selected in terms of their mutation rates and the amplification factors of the fittest phenotypes. We interpret the results as the compromise between the characteristic time a lineage takes to reach its fittest phenotype by crossing the neutral valley and the selective value of the sequences that form the lineages. A main conclusion is drawn: for finite population sizes, the survival probability of the lineage that arrives first to the fittest phenotype rises significantly.

Genome destabilizing mutator alleles drive specific mutational trajectories in Saccharomyces cerevisiae

In addition to environmental factors and intrinsic variations in base substitution rates, specific genome-destabilizing mutations can shape the mutational trajectory of genomes. How specific alleles influence the nature and position of accumulated mutations in a genomic context is largely unknown. Understanding the impact of genome-destabilizing alleles is particularly relevant to cancer genomes where biased mutational signatures are identifiable. We first created a more complete picture of cellular pathways that impact mutation rate using a primary screen to identify essential Saccharomyces cerevisiae gene mutations that cause mutator phenotypes. Drawing primarily on new alleles identified in this resource, we measure the impact of diverse mutator alleles on mutation patterns directly by whole-genome sequencing of 68 mutation-accumulation strains derived from wild-type and 11 parental mutator genotypes. The accumulated mutations differ across mutator strains, displaying base-substitution biases, allele-specific mutation hotspots, and break-associated mutation clustering. For example, in mutants of POLα and the Cdc13–Stn1–Ten1 complex, we find a distinct subtelomeric bias for mutations that we show is independent of the target sequence. Together our data suggest that specific genome-instability mutations are sufficient to drive discrete mutational signatures, some of which share properties with mutation patterns seen in tumors. Thus, in a population of cells, genome-instability mutations could influence clonal evolution by establishing discrete mutational trajectories for genomes.

Endogenous stress caused by faulty oxidation reactions fosters evolution of 2,4-dinitrotoluene-degrading bacteria

Environmental strain Burkholderia sp. DNT mineralizes the xenobiotic compound 2,4-dinitrotoluene (DNT) owing to the catabolic dnt genes borne by plasmid DNT, but the process fails to promote significant growth. To investigate this lack of physiological return of such an otherwise complete metabolic route, cells were exposed to DNT under various growth conditions and the endogenous formation of reactive oxygen species (ROS) monitored in single bacteria. These tests revealed the buildup of a strong oxidative stress in the population exposed to DNT. By either curing the DNT plasmid or by overproducing the second activity of the biodegradation route (DntB) we could trace a large share of ROS production to the first reaction of the route, which is executed by the multicomponent dioxygenase encoded by the dntA gene cluster. Naphthalene, the ancestral substrate of the dioxygenase from which DntA has evolved, also caused significant ROS formation. That both the old and the new substrate brought about a considerable cellular stress was indicative of a still-evolving DntA enzyme which is neither optimal any longer for naphthalene nor entirely advantageous yet for growth of the host strain on DNT. We could associate endogenous production of ROS with likely error-prone repair mechanisms of DNA damage, and the ensuing stress-induced mutagenesis in cells exposed to DNT. It is thus plausible that the evolutionary roadmap for biodegradation of xenobiotic compounds like DNT was largely elicited by mutagenic oxidative stress caused by faulty reactions of precursor enzymes with novel but structurally related substrates-to-be.

Many bacteria have acquired the capacity of metabolizing chemical compounds that have never been in the Biosphere before the onset of contemporary synthetic chemistry. However, the factors that shape the new metabolic properties of such microorganisms remain obscure. We examined the performance of a still-evolving metabolic pathway for biodegradation of 2,4-dinitrotoluene (DNT, an archetypal xenobiotic compound) borne by a Burkholderia strain isolated from soil in an ammunition plant. The biodegradation pathway likely arose from a precursor set of genes for catabolism of naphthalene (although Burkholderia does not degrade this compound any longer), and is now advancing towards the new substrate, DNT. We found that the action of the first enzyme of the biodegradation pathway, a Rieske-type dioxygenase, on the still-suboptimal substrate (DNT) generates a high level of endogenous reactive oxygen species. This, in turn, damages DNA and increases mutagenesis, ultimately resulting in the creation of novelty that may foster evolution of xenobiotic-degrading variants of the strain hosting the biodegradation pathway. The very metabolic problem thus somehow seems to stimulate the exploration of the solution space. Our data is fully consistent with the notion that stress caused by faulty dioxygenation of DNT accelerates the rate of bacterial evolution.