Mutation rate dynamics in a bacterial population reflect tension between adaptation and genetic load

Culture Volume Influences the Dynamics of Adaptation under Long-Term Stationary Phase

Escherichia coli and many other bacterial species, which are incapable of sporulation, can nevertheless survive within resource exhausted media by entering a state termed long-term stationary phase (LTSP). We have previously shown that E. coli populations adapt genetically under LTSP in an extremely convergent manner. Here, we examine how the dynamics of LTSP genetic adaptation are influenced by varying a single parameter of the experiment-culture volume. We find that culture volume affects survival under LTSP, with viable counts decreasing as volumes increase. Across all volumes, mutations accumulate with time, and the majority of mutations accumulated demonstrate signals of being adaptive. However, positive selection appears to affect mutation accumulation more strongly at higher, compared with lower volumes. Finally, we find that several similar genes are likely involved in adaptation across volumes. However, the specific mutations within these genes that contribute to adaptation can vary in a consistent manner. Combined, our results demonstrate how varying a single parameter of an evolutionary experiment can substantially influence the dynamics of observed adaptation.

Evolutionary epidemiology of Streptococcus iniae: Linking mutation rate dynamics with adaptation to novel immunological landscapes

Pathogens continuously adapt to changing host environments where variation in their virulence and antigenicity is critical to their long-term evolutionary success. The emergence of novel variants is accelerated in microbial mutator strains (mutators) deficient in DNA repair genes, most often from mismatch repair and oxidized-guanine repair systems (MMR and OG respectively). Bacterial MMR/OG mutants are abundant in clinical samples and show increased adaptive potential in experimental infection models, yet the role of mutators in the epidemiology and evolution of infectious disease is not well understood. Here we investigated the role of mutation rate dynamics in the evolution of a broad host range pathogen, Streptococcus iniae, using a set of 80 strains isolated globally over 40 years. We have resolved phylogenetic relationships using non-recombinant core genome variants, measured in vivo mutation rates by fluctuation analysis, identified variation in major MMR/OG genes and their regulatory regions, and phenotyped the major traits determining virulence in streptococci. We found that both mutation rate and MMR/OG genotype are remarkably conserved within phylogenetic clades but significantly differ between major phylogenetic lineages. Further, variation in MMR/OG loci correlates with occurrence of atypical virulence-associated phenotypes, infection in atypical hosts (mammals), and atypical (osseous) tissue of a vaccinated primary host. These findings suggest that mutators are likely to facilitate adaptations preceding major diversification events and may promote emergence of variation permitting colonization of a novel host tissue, novel host taxa (host jumps), and immune-escape in the vaccinated host.

Effect of arsenite and growth in biofilm conditions on the evolution of Thiomonas sp. CB2

Thiomonas bacteria are ubiquitous at acid mine drainage sites and play key roles in the remediation of water at these locations by oxidizing arsenite to arsenate, favouring the sorption of arsenic by iron oxides and their coprecipitation. Understanding the adaptive capacities of these bacteria is crucial to revealing how they persist and remain active in such extreme conditions. Interestingly, it was previously observed that after exposure to arsenite, when grown in a biofilm, some strains of Thiomonas bacteria develop variants that are more resistant to arsenic. Here, we identified the mechanisms involved in the emergence of such variants in biofilms. We found that the percentage of variants generated increased in the presence of high concentrations of arsenite (5.33 mM), especially in the detached cells after growth under biofilm-forming conditions. Analysis of gene expression in the parent strain CB2 revealed that genes involved in DNA repair were upregulated in the conditions where variants were observed. Finally, we assessed the phenotypes and genomes of the subsequent variants generated to evaluate the number of mutations compared to the parent strain. We determined that multiple point mutations accumulated after exposure to arsenite when cells were grown under biofilm conditions. Some of these mutations were found in what is referred to as ICE19, a genomic island (GI) carrying arsenic-resistance genes, also harbouring characteristics of an integrative and conjugative element (ICE). The mutations likely favoured the excision and duplication of this GI. This research aids in understanding how Thiomonas bacteria adapt to highly toxic environments, and, more generally, provides a window to bacterial genome evolution in extreme environments.

The target of selection matters: An established resistance-development-time negative genetic trade-off is not found when selecting on development time

Trade-offs are fundamental to evolutionary outcomes and play a central role in eco-evolutionary theory. They are often examined by experimentally selecting on one life-history trait and looking for negative correlations in other traits. For example, populations of the moth Plodia interpunctella selected to resist viral infection show a life-history cost with longer development times. However, we rarely examine whether the detection of such negative genetic correlations depends on the trait on which we select. Here, we examine a well-characterized negative genotypic trade-off between development time and resistance to viral infection in the moth Plodia interpunctella and test whether selection on a phenotype known to be a cost of resistance (longer development time) leads to the predicted correlated increase in resistance. If there is tight pleiotropic relationship between genes that determine development time and resistance underpinning this trade-off, we might expect increased resistance when we select on longer development time. However, we show that selecting for longer development time in this system selects for reduced resistance when compared to selection for shorter development time. This shows how phenotypes typically characterized by a trade-off can deviate from that trade-off relationship, and suggests little genetic linkage between the genes governing viral resistance and those that determine response to selection on the key life-history trait. Our results are important for both selection strategies in applied biological systems and for evolutionary modelling of host-parasite interactions.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

The adaptive potential of circular DNA accumulation in ageing cells

Carefully maintained and precisely inherited chromosomal DNA provides long-term genetic stability, but eukaryotic cells facing environmental challenges can benefit from the accumulation of less stable DNA species. Circular DNA molecules lacking centromeres segregate randomly or asymmetrically during cell division, following non-Mendelian inheritance patterns that result in high copy number instability and massive heterogeneity across populations. Such circular DNA species, variously known as extrachromosomal circular DNA (eccDNA), microDNA, double minutes or extrachromosomal DNA (ecDNA), are becoming recognised as a major source of the genetic variation exploited by cancer cells and pathogenic eukaryotes to acquire drug resistance. In budding yeast, circular DNA molecules derived from the ribosomal DNA (ERCs) have been long known to accumulate with age, but it is now clear that aged yeast also accumulate other high-copy protein-coding circular DNAs acquired through both random and environmentally-stimulated recombination processes. Here, we argue that accumulation of circular DNA provides a reservoir of heterogeneous genetic material that can allow rapid adaptation of aged cells to environmental insults, but avoids the negative fitness impacts on normal growth of unsolicited gene amplification in the young population.

Resist or perish: Fate of a microbial population subjected to a periodic presence of antimicrobial

The evolution of antimicrobial resistance can be strongly affected by variations of antimicrobial concentration. Here, we study the impact of periodic alternations of absence and presence of antimicrobial on resistance evolution in a microbial population, using a stochastic model that includes variations of both population composition and size, and fully incorporates stochastic population extinctions. We show that fast alternations of presence and absence of antimicrobial are inefficient to eradicate the microbial population and strongly favor the establishment of resistance, unless the antimicrobial increases enough the death rate. We further demonstrate that if the period of alternations is longer than a threshold value, the microbial population goes extinct upon the first addition of antimicrobial, if it is not rescued by resistance. We express the probability that the population is eradicated upon the first addition of antimicrobial, assuming rare mutations. Rescue by resistance can happen either if resistant mutants preexist, or if they appear after antimicrobial is added to the environment. Importantly, the latter case is fully prevented by perfect biostatic antimicrobials that completely stop division of sensitive microorganisms. By contrast, we show that the parameter regime where treatment is efficient is larger for biocidal drugs than for biostatic drugs. This sheds light on the respective merits of different antimicrobial modes of action.

Mutations Beget More Mutations-Rapid Evolution of Mutation Rate in Response to the Risk of Runaway Accumulation

The rapidity with which the mutation rate evolves could greatly impact evolutionary patterns. Nevertheless, most studies simply assume a constant rate in the time scale of interest (Kimura 1983; Drake 1991; Kumar 2005; Li 2007; Lynch 2010). In contrast, recent studies of somatic mutations suggest that the mutation rate may vary by several orders of magnitude within a lifetime (Kandoth et al. 2013; Lawrence et al. 2013). To resolve the discrepancy, we now propose a runaway model, applicable to both the germline and soma, whereby mutator mutations form a positive-feedback loop. In this loop, any mutator mutation would increase the rate of acquiring the next mutator, thus triggering a runaway escalation in mutation rate. The process can be initiated more readily if there are many weak mutators than a few strong ones. Interestingly, even a small increase in the mutation rate at birth could trigger the runaway process, resulting in unfit progeny. In slowly reproducing species, the need to minimize the risk of this uncontrolled accumulation would thus favor setting the mutation rate low. In comparison, species that starts and ends reproduction sooner do not face the risk and may set the baseline mutation rate higher. The mutation rate would evolve in response to the risk of runaway mutation, in particular, when the generation time changes. A rapidly evolving mutation rate may shed new lights on many evolutionary phenomena (Elango et al. 2006; Thomas et al. 2010, 2018; Langergraber et al. 2012; Besenbacher et al. 2019).

Experimental Evolution of Legume Symbionts: What Have We Learnt?

Rhizobia, the nitrogen-fixing symbionts of legumes, are polyphyletic bacteria distributed in many alpha- and beta-proteobacterial genera. They likely emerged and diversified through independent horizontal transfers of key symbiotic genes. To replay the evolution of a new rhizobium genus under laboratory conditions, the symbiotic plasmid of Cupriavidus taiwanensis was introduced in the plant pathogen Ralstonia solanacearum, and the generated proto-rhizobium was submitted to repeated inoculations to the C. taiwanensis host, Mimosa pudica L.. This experiment validated a two-step evolutionary scenario of key symbiotic gene acquisition followed by genome remodeling under plant selection. Nodulation and nodule cell infection were obtained and optimized mainly via the rewiring of regulatory circuits of the recipient bacterium. Symbiotic adaptation was shown to be accelerated by the activity of a mutagenesis cassette conserved in most rhizobia. Investigating mutated genes led us to identify new components of R. solanacearum virulence and C. taiwanensis symbiosis. Nitrogen fixation was not acquired in our short experiment. However, we showed that post-infection sanctions allowed the increase in frequency of nitrogen-fixing variants among a non-fixing population in the M. pudica-C. taiwanensis system and likely allowed the spread of this trait in natura. Experimental evolution thus provided new insights into rhizobium biology and evolution.

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.

Low mutational load and high mutation rate variation in gut commensal bacteria

Bacteria generally live in species-rich communities, such as the gut microbiota. Yet little is known about bacterial evolution in natural ecosystems. Here, we followed the long-term evolution of commensal Escherichia coli in the mouse gut. We observe the emergence of mutation rate polymorphism, ranging from wild-type levels to 1,000-fold higher. By combining experiments, whole-genome sequencing, and in silico simulations, we identify the molecular causes and explore the evolutionary conditions allowing these hypermutators to emerge and coexist within the microbiota. The hypermutator phenotype is caused by mutations in DNA polymerase III proofreading and catalytic subunits, which increase mutation rate by approximately 1,000-fold and stabilise hypermutator fitness, respectively. Strong mutation rate variation persists for >1,000 generations, with coexistence between lineages carrying 4 to >600 mutations. The in vivo molecular evolution pattern is consistent with fitness effects of deleterious mutations s_d ≤ 10⁻⁴/generation, assuming a constant effect or exponentially distributed effects with a constant mean. Such effects are lower than typical in vitro estimates, leading to a low mutational load, an inference that is observed in in vivo and in vitro competitions. Despite large numbers of deleterious mutations, we identify multiple beneficial mutations that do not reach fixation over long periods of time. This indicates that the dynamics of beneficial mutations are not shaped by constant positive Darwinian selection but could be explained by other evolutionary mechanisms that maintain genetic diversity. Thus, microbial evolution in the gut is likely characterised by partial sweeps of beneficial mutations combined with hitchhiking of slightly deleterious mutations, which take a long time to be purged because they impose a low mutational load. The combination of these two processes could allow for the long-term maintenance of intraspecies genetic diversity, including mutation rate polymorphism. These results are consistent with the pattern of genetic polymorphism that is emerging from metagenomics studies of the human gut microbiota, suggesting that we have identified key evolutionary processes shaping the genetic composition of this community.

Weak-effect deleterious mutations and negative frequency–dependent selection, acting on beneficial mutations, shape the dynamics of molecular evolution within the mouse gut microbiota.

Adapting the engine to the fuel: mutator populations can reduce the mutational load by reorganizing their genome structure

Background

Mutators are common in bacterial populations, both in natural isolates and in the lab. The fate of these lineages, which mutation rate is increased up to 100 ×, has long been studied using population genetics models, showing that they can spread in a population following an environmental change. However in stable conditions, they suffer from the increased mutational load, hence being overcome by non-mutators. However, these results don’t take into account the fact that an elevated mutation rate can impact the genetic structure, hence changing the sensitivity of the population to mutations. Here we used Aevol, an in silico experimental evolution platform in which genomic structures are free to evolve, in order to study the fate of mutator populations evolving for a long time in constant conditions.

Results

Starting from wild-types that were pre-evolved for 300,000 generations, we let 100 mutator populations (point mutation rate ×100) evolve for 100,000 further generations in constant conditions. As expected all populations initially undergo a fitness loss. However, after that the mutator populations started to recover. Most populations ultimately recovered their ancestors fitness, and a significant fraction became even fitter than the non-mutator control clones that evolved in parallel. By analyzing the genomes of the mutators, we show that the fitness recovery is due to two mechanisms: i. an increase in robustness through compaction of the coding part of the mutator genomes, ii. an increase of the selection coefficient that decreases the mean-fitness of the population. Strikingly the latter is due to the accumulation of non-coding sequences in the mutators genomes.

Conclusion

Our results show that the mutational burden that is classically thought to be associated with mutator phenotype is escapable. On the long run mutators adapted their genomes and reshaped the distribution of mutation effects. Therewith the lineage is able to recover fitness even though the population still suffers the elevated mutation rate. Overall these results change our view of mutator dynamics: by being able to reduce the deleterious effect of the elevated mutation rate, mutator populations may be able to last for a very long time; A situation commonly observed in nature.

On the feasibility of saltational evolution

In evolutionary biology, it is generally assumed that evolution occurs in the weak mutation limit, that is, the frequency of multiple mutations simultaneously occurring in the same genome and the same generation is negligible. We employ mathematical modeling to show that, although under the typical parameter values of the evolutionary process the probability of multimutational leaps is indeed low, they might become substantially more likely under stress, when the mutation rate is dramatically elevated. We hypothesize that stress-induced mutagenesis in microbes is an evolvable adaptive strategy. Multimutational leaps might matter also in other cases of substantially increased mutation rate, such as growing tumors or evolution of primordial replicators.

Is evolution always gradual or can it make leaps? We examine a mathematical model of an evolutionary process on a fitness landscape and obtain analytic solutions for the probability of multimutation leaps, that is, several mutations occurring simultaneously, within a single generation in 1 genome, and being fixed all together in the evolving population. The results indicate that, for typical, empirically observed combinations of the parameters of the evolutionary process, namely, effective population size, mutation rate, and distribution of selection coefficients of mutations, the probability of a multimutation leap is low, and accordingly the contribution of such leaps is minor at best. However, we show that, taking sign epistasis into account, leaps could become an important factor of evolution in cases of substantially elevated mutation rates, such as stress-induced mutagenesis in microbes. We hypothesize that stress-induced mutagenesis is an evolvable adaptive strategy.

The prokaryotic Argonaute proteins enhance homology sequence-directed recombination in bacteria

Argonaute proteins are present and conserved in all domains of life. Recently characterized prokaryotic Argonaute proteins (pAgos) participates in host defense by DNA interference. Here, we report that the Natronobacterium gregoryi Argonaute (NgAgo) enhances gene insertions or deletions in Pasteurella multocida and Escherichia coli at efficiencies of 80–100%. Additionally, the effects are in a homologous arms-dependent but guide DNA- and potential enzyme activity-independent manner. Interestingly, such effects were also observed in other pAgos fragments including Thermus thermophilus Argonaute (TtAgo), Aquifex aeolicus Argonaute (AaAgo) and Pyrococcus furiosus Argonaute (PfAgo). The underlying mechanism of the NgAgo system is a positive selection process mainly through its PIWI-like domain interacting with recombinase A (recA) to enhance recA-mediated DNA strand exchange. Our study reveals a novel system for enhancing homologous sequence-guided gene editing in bacteria.

SpCLUST: Towards a fast and reliable clustering for potentially divergent biological sequences

This paper presents SpCLUST, a new C++ package that takes a list of sequences as input, aligns them with MUSCLE, computes their similarity matrix in parallel and then performs the clustering. SpCLUST extends a previously released software by integrating additional scoring matrices which enables it to cover the clustering of amino-acid sequences. The similarity matrix is now computed in parallel according to the master/slave distributed architecture, using MPI. Performance analysis, realized on two real datasets of 100 nucleotide sequences and 1049 amino-acids ones, show that the resulting library substantially outperforms the original Python package. The proposed package was also intensively evaluated on simulated and real genomic and protein data sets. The clustering results were compared to the most known traditional tools, such as UCLUST, CD-HIT and DNACLUST. The comparison showed that SpCLUST outperforms the other tools when clustering divergent sequences, and contrary to the others, it does not require any user intervention or prior knowledge about the input sequences.

Long-term evolution of the natural isolate of Escherichia coli 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon

In vitro experimental evolution has taught us many lessons on the molecular bases of adaptation. To move towards more natural settings, evolution in the mice gut has been successfully performed. Yet, these experiments suffered from the use of laboratory strains as well as the use of axenic or streptomycin-treated mice to maintain the inoculated strains. To circumvent these limitations, we conducted a one-year experimental evolution in vivo using a natural isolate of E. coli, strain 536, in conditions mimicking as much as possible natural environment with mother-to-offspring microbiota transmission. Mice were then distributed in 24 independent cages and separated into two different diets: a regular one (chow diet, CD) and high-fat and high-sugar one (Western Diet, WD). Genome sequences revealed an early and rapid selection during the breastfeeding period that selected the constitutive expression of the well-characterized lactose operon. E. coli was lost significantly more in CD than WD; however, we could not detect any genomic signature of selection, nor any diet specificities during the later part of the experiments. The apparently neutral evolution presumably due to low population size maintained nevertheless at high frequency the early selected mutations affecting lactose regulation. The rapid loss of lactose operon regulation challenges the idea that plastic gene expression is both optimal and stable in the wild.

Inactivation of a Mismatch-Repair System Diversifies Genotypic Landscape of Escherichia coli During Adaptive Laboratory Evolution

Adaptive laboratory evolution (ALE) is used to find causal mutations that underlie improved strain performance under the applied selection pressure. ALE studies have revealed that mutator populations tend to outcompete their non-mutator counterparts following the evolutionary trajectory. Among them, mutS-inactivated mutator cells, characterize d by a dysfunctional methyl-mismatch repair system, are frequently found in ALE experiments. Here, we examined mutS inactivation as an approach to facilitate ALE of Escherichia coli. The wild-type E. coli MG1655 and mutS knock-out derivative (ΔmutS) were evolved in parallel for 800 generations on lactate or glycerol minimal media in a serial-transfer experiment. Whole-genome re-sequencing of each lineage at 100-generation intervals revealed that (1) mutations emerge rapidly in the ΔmutS compared to in the wild-type strain; (2) mutations were more than fourfold higher in the ΔmutS strain at the end-point populations compared to the wild-type strain; and (3) a significant number of random mutations accumulated in the ΔmutS strains. We then measured the fitness of the end-point populations on an array of non-adaptive carbon sources. Interestingly, collateral fitness increases on non-adaptive carbon sources were more pronounced in the ΔmutS strains than the parental strain. Fitness measurement of single mutants revealed that the collateral fitness increase seen in the mutator lineages can be attributed to a pool of random mutations. Together, this study demonstrates that short-term mutator ALE extensively expands possible genotype space, resulting in versatile bacteria with elevated fitness levels across various carbon sources.

The Mutator Phenotype: Adapting Microbial Evolution to Cancer Biology

The mutator phenotype hypothesis was postulated almost 40 years ago to reconcile the observation that while cancer cells display widespread mutational burden, acquisition of mutations in non-transformed cells is a rare event. Moreover, it also suggested that cancer evolution could be fostered by increased genome instability. Given the evolutionary conservation throughout the tree of life and the genetic tractability of model organisms, yeast and bacterial species pioneered studies to dissect the functions of genes required for genome maintenance (caretaker genes) or for cell growth control (gatekeeper genes). In this review, we first provide an overview of what we learned from model organisms about the roles of these genes and the genome instability that arises as a consequence of their dysregulation. We then discuss our current understanding of how mutator phenotypes shape the evolution of bacteria and yeast species. We end by bringing clinical evidence that lessons learned from single-cell organisms can be applied to tumor evolution.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

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.

IAMBEE: a web-service for the identification of adaptive pathways from parallel evolved clonal populations

IAMBEE is a web server designed for the Identification of Adaptive Mutations in Bacterial Evolution Experiments (IAMBEE). Input data consist of genotype information obtained from independently evolved clonal populations or strains that show the same adapted behavior (phenotype). To distinguish adaptive from passenger mutations, IAMBEE searches for neighborhoods in an organism-specific interaction network that are recurrently mutated in the adapted populations. This search for recurrently mutated network neighborhoods, as proxies for pathways is driven by additional information on the functional impact of the observed genetic changes and their dynamics during adaptive evolution. In addition, the search explicitly accounts for the differences in mutation rate between the independently evolved populations. Using this approach, IAMBEE allows exploiting parallel evolution to identify adaptive pathways. The web-server is freely available at http://bioinformatics.intec.ugent.be/iambee/ with no login requirement.

Interference Effects of Deleterious and Beneficial Mutations in Large Asexual Populations

Linked beneficial and deleterious mutations are known to decrease the fixation probability of a favorable mutation in large asexual populations. While the hindering effect of strongly deleterious mutations on adaptive evolution has been well studied, how weakly deleterious mutations, either in isolation or with superior beneficial mutations, influence the rate of adaptation has not been fully explored. When the selection against the deleterious mutations is weak, the beneficial mutant can fix in many genetic backgrounds, besides the one it arose on. Here, taking this factor into account, I obtain an accurate analytical expression for the fixation probability of a beneficial mutant in an asexual population at mutation-selection balance. I then exploit this result along with clonal interference theory to investigate the joint effect of linked beneficial and deleterious mutations on the rate of adaptation, and identify parameter regions where it is reduced due to interference by either beneficial or deleterious or both types of mutations. I also study the evolution of mutation rates in adapting asexual populations, and find that linked beneficial mutations have a stronger influence than the deleterious mutations on mutator fixation.

Changes in Intrinsic Antibiotic Susceptibility during a Long-Term Evolution Experiment with Escherichia coli

High-level resistance often evolves when populations of bacteria are exposed to antibiotics, by either mutations or horizontally acquired genes. There is also variation in the intrinsic resistance levels of different bacterial strains and species that is not associated with any known history of exposure. In many cases, evolved resistance is costly to the bacteria, such that resistant types have lower fitness than their progenitors in the absence of antibiotics. Some longer-term studies have shown that bacteria often evolve compensatory changes that overcome these tradeoffs, but even those studies have typically lasted only a few hundred generations. In this study, we examine changes in the susceptibilities of 12 populations of Escherichia coli to 15 antibiotics after 2,000 and 50,000 generations without exposure to any antibiotic. On average, the evolved bacteria were more susceptible to most antibiotics than was their ancestor. The bacteria at 50,000 generations tended to be even more susceptible than after 2,000 generations, although most of the change occurred during the first 2,000 generations. Despite the general trend toward increased susceptibility, we saw diverse outcomes with different antibiotics. For streptomycin, which was the only drug to which the ancestral strain was highly resistant, none of the evolved lines showed any increased susceptibility. The independently evolved lineages often exhibited correlated responses to the antibiotics, with correlations usually corresponding to their modes of action. On balance, our study shows that bacteria with low levels of intrinsic resistance often evolve to become even more susceptible to antibiotics in the absence of corresponding selection.IMPORTANCE Resistance to antibiotics often evolves when bacteria encounter antibiotics. However, bacterial strains and species without any known exposure to these drugs also vary in their intrinsic susceptibility. In many cases, evolved resistance has been shown to be costly to the bacteria, such that resistant types have reduced competitiveness relative to their sensitive progenitors in the absence of antibiotics. In this study, we examined changes in the susceptibilities of 12 populations of Escherichia coli to 15 antibiotics after 2,000 and 50,000 generations without exposure to any drug. The evolved bacteria tended to become more susceptible to most antibiotics, with most of the change occurring during the first 2,000 generations, when the bacteria were undergoing rapid adaptation to their experimental conditions. On balance, our findings indicate that bacteria with low levels of intrinsic resistance can, in the absence of relevant selection, become even more susceptible to antibiotics.

Directed evolution of Escherichia coli with lower-than-natural plasmid mutation rates

Unwanted evolution of designed DNA sequences limits metabolic and genome engineering efforts. Engineered functions that are burdensome to host cells and slow their replication are rapidly inactivated by mutations, and unplanned mutations with unpredictable effects often accumulate alongside designed changes in large-scale genome editing projects. We developed a directed evolution strategy, Periodic Reselection for Evolutionarily Reliable Variants (PResERV), to discover mutations that prolong the function of a burdensome DNA sequence in an engineered organism. Here, we used PResERV to isolate Escherichia coli cells that replicate ColE1-type plasmids with higher fidelity. We found mutations in DNA polymerase I and in RNase E that reduce plasmid mutation rates by 6- to 30-fold. The PResERV method implicitly selects to maintain the growth rate of host cells, and high plasmid copy numbers and gene expression levels are maintained in some of the evolved E. coli strains, indicating that it is possible to improve the genetic stability of cellular chassis without encountering trade-offs in other desirable performance characteristics. Utilizing these new antimutator E. coli and applying PResERV to other organisms in the future promises to prevent evolutionary failures and unpredictability to provide a more stable genetic foundation for synthetic biology.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

A Rapid and Precise Mutation-Activated Fluorescence Reporter for Analyzing Acute Mutagenesis Frequency

Mutagenesis reporters are critical for quantifying genome stability. However, current methods rely on cell survival/death to report mutation, which takes weeks and prevents evaluation of acute or time-dependent changes. Existing methods also have other limitations, such as cell type restrictions. Using our discovery that mCherryFP fluorescence depends on residue Trp98, we replaced this codon with a stop codon to generate a mutation biosensor (termed CherryOFF), with a green fluorescence protein (GFP) as an internal control. We found that the red fluorescence of this biosensor is activated by a specific A/T-G/C nucleotide transition. Compared with the established hypoxanthine phosphoribosyl transferase assay, our reporter has similar or better ability to detect changes of mutation frequency induced by physical/chemical mutagens or manipulation of mutation-related genes. Furthermore, CherryOFF-GFP can report mutagenesis independently of cell-death events, can be adapted to many cell types, and can generate readouts within 1 day for the measurement of acute or time-dependent events.

Spatial Vulnerabilities of the Escherichia coli Genome to Spontaneous Mutations Revealed with Improved Duplex Sequencing

Investigation of spontaneous mutations by next-generation sequencing technology has attracted extensive attention lately due to the fundamental roles of spontaneous mutations in evolution and pathological processes. However, these studies only focused on the mutations accumulated through many generations during long-term (possibly be years of) culturing, but not the freshly generated mutations that occur at very low frequencies. In this study, we established a molecularly barcoded deep sequencing strategy to detect low abundant spontaneous mutations in genomes of bacteria cell cultures. Genome-wide spontaneous mutations in 15 Escherichia coli cell culture samples were defined with a high confidence (P < 0.01). We also developed a hotspot-calling approach based on the run-length encoding algorithm to find the genomic regions that are vulnerable to the spontaneous mutations. The hotspots for the mutations appeared to be highly conserved across the bacteria samples. Further biological annotation of these regions indicated that most of the spontaneous mutations were located at the repeat domains or nonfunctional domains of the genomes, suggesting the existence of mechanisms that could somehow prevent the occurrence of mutations in crucial genic areas. This study provides a more faithful picture of mutation occurrence and spectra in a single expansion process without long-term culturing.

High mutation rates limit evolutionary adaptation in Escherichia coli

Mutation is fundamental to evolution, because it generates the genetic variation on which selection can act. In nature, genetic changes often increase the mutation rate in systems that range from viruses and bacteria to human tumors. Such an increase promotes the accumulation of frequent deleterious or neutral alleles, but it can also increase the chances that a population acquires rare beneficial alleles. Here, we study how up to 100-fold increases in Escherichia coli's genomic mutation rate affect adaptive evolution. To do so, we evolved multiple replicate populations of asexual E. coli strains engineered to have four different mutation rates for 3000 generations in the laboratory. We measured the ability of evolved populations to grow in their original environment and in more than 90 novel chemical environments. In addition, we subjected the populations to whole genome population sequencing. Although populations with higher mutation rates accumulated greater genetic diversity, this diversity conveyed benefits only for modestly increased mutation rates, where populations adapted faster and also thrived better than their ancestors in some novel environments. In contrast, some populations at the highest mutation rates showed reduced adaptation during evolution, and failed to thrive in all of the 90 alternative environments. In addition, they experienced a dramatic decrease in mutation rate. Our work demonstrates that the mutation rate changes the global balance between deleterious and beneficial mutational effects on fitness. In contrast to most theoretical models, our experiments suggest that this tipping point already occurs at the modest mutation rates that are found in the wild.

Population heterogeneity in microbial bioprocesses: origin, analysis, mechanisms, and future perspectives

Population heterogeneity is omnipresent in all bioprocesses even in homogenous environments. Its origin, however, is only so well understood that potential strategies like bet-hedging, noise in gene expression and division of labour that lead to population heterogeneity can be derived from experimental studies simulating the dynamics in industrial scale bioprocesses. This review aims at summarizing the current state of the different parts of single cell studies in bioprocesses. This includes setups to visualize different phenotypes of single cells, computational approaches connecting single cell physiology with environmental influence and special cultivation setups like scale-down reactors that have been proven to be useful to simulate large-scale conditions. A step in between investigation of populations and single cells is studying subpopulations with distinct properties that differ from the rest of the population with sub-omics methods which are also presented here. Moreover, the current knowledge about population heterogeneity in bioprocesses is summarized for relevant industrial production hosts and mixed cultures, as they provide the unique opportunity to distribute metabolic burden and optimize production processes in a way that is impossible in traditional monocultures. In the end, approaches to explain the underlying mechanism of population heterogeneity and the evidences found to support each hypothesis are presented. For instance, population heterogeneity serving as a bet-hedging strategy that is used as coordinated action against bioprocess-related stresses while at the same time spreading the risk between individual cells as it ensures the survival of least a part of the population in any environment the cells encounter.

Evolution of a Dominant Natural Isolate of Escherichia coli in the Human Gut over the Course of a Year Suggests a Neutral Evolution with Reduced Effective Population Size

In vitro and in vivo evolution experiments on Escherichia coli revealed several principles of bacterial adaptation. However, few data are available in the literature describing the behavior of E. coli in its natural environment. We attempted here to study the evolution in the human gut of a commensal dominant E. coli clone, ED1a belonging to the B2 phylogroup, through a longitudinal genomic study. We sequenced 24 isolates sampled at three different time points within a healthy individual over almost a year. We computed a mutation rate of 6.90 × 10^-7 mutations per base per year of the chromosome for E. coli ED1a in healthy human gut. We observed very limited genomic diversity and could not detect any evidence of selection, in contrast to what is observed in experimental evolution over a similar length of time. We therefore suggest that ED1a, being well adapted to the healthy human gut, evolves mostly neutrally with a low effective population size (N_e of ≈500 to 1,700).IMPORTANCE In this study, we follow the genomic fate of a dominant clone of Escherichia coli in the human gut of a healthy individual over about a year. We could compute a low annual mutation rate that supports low diversity, and we could not retrieve any clear signature of selection. These observations support a neutral evolution of E. coli in the human gut, compatible with a very limited effective population size that deviates drastically with the observations made previously in experimental evolution.

Adaptive evolution of genomically recoded Escherichia coli

Efforts are underway to construct several recoded genomes anticipated to exhibit multivirus resistance, enhanced nonstandard amino acid (nsAA) incorporation, and capability for synthetic biocontainment. Although our laboratory pioneered the first genomically recoded organism (Escherichia coli strain C321.∆A), its fitness is far lower than that of its nonrecoded ancestor, particularly in defined media. This fitness deficit severely limits its utility for nsAA-linked applications requiring defined media, such as live cell imaging, metabolic engineering, and industrial-scale protein production. Here, we report adaptive evolution of C321.∆A for more than 1,000 generations in independent replicate populations grown in glucose minimal media. Evolved recoded populations significantly exceeded the growth rates of both the ancestral C321.∆A and nonrecoded strains. We used next-generation sequencing to identify genes mutated in multiple independent populations, and we reconstructed individual alleles in ancestral strains via multiplex automatable genome engineering (MAGE) to quantify their effects on fitness. Several selective mutations occurred only in recoded evolved populations, some of which are associated with altering the translation apparatus in response to recoding, whereas others are not apparently associated with recoding, but instead correct for off-target mutations that occurred during initial genome engineering. This report demonstrates that laboratory evolution can be applied after engineering of recoded genomes to streamline fitness recovery compared with application of additional targeted engineering strategies that may introduce further unintended mutations. In doing so, we provide the most comprehensive insight to date into the physiology of the commonly used C321.∆A strain.

The Odyssey of the Ancestral Escherich Strain through Culture Collections: an Example of Allopatric Diversification

More than a century ago, Theodor Escherich isolated the bacterium that was to become Escherichia coli, one of the most studied organisms. Not long after, the strain began an odyssey and landed in many laboratories across the world. As laboratory culture conditions could be responsible for major changes in bacterial strains, we conducted a genome analysis of isolates of this emblematic strain from different culture collections (England, France, the United States, Germany). Strikingly, many discrepancies between the isolates were observed, as revealed by multilocus sequence typing (MLST), the presence of virulence-associated genes, core genome MLST, and single nucleotide polymorphism/indel analyses. These differences are correlated with the phylogeographic history of the strain and were due to an unprecedented number of mutations in coding DNA repair functions such as mismatch repair (MutL) and oxidized guanine nucleotide pool cleaning (MutT), conferring a specific mutational spectrum and leading to a mutator phenotype. The mutator phenotype was probably acquired during subculturing and corresponded to second-order selection. Furthermore, all of the isolates exhibited hypersusceptibility to antibiotics due to mutations in efflux pump- and porin-encoding genes, as well as a specific mutation in the sigma factor-encoding gene rpoS. These defects reflect a self-preservation and nutritional competence tradeoff allowing survival under the starvation conditions imposed by storage. From a clinical point of view, dealing with such mutator strains can lead microbiologists to draw false conclusions about isolate relatedness and may impact therapeutic effectiveness. IMPORTANCE Mutator phenotypes have been described in laboratory-evolved bacteria, as well as in natural isolates. Several genes can be impacted, each of them being associated with a typical mutational spectrum. By studying one of the oldest strains available, the ancestral Escherich strain, we were able to identify its mutator status leading to tremendous genetic diversity among the isolates from various collections and allowing us to reconstruct the phylogeographic history of the strain. This mutator phenotype was probably acquired during the storage of the strain, promoting adaptation to a specific environment. Other mutations in rpoS and efflux pump- and porin-encoding genes highlight the acclimatization of the strain through self-preservation and nutritional competence regulation. This strain history can be viewed as unintentional experimental evolution in culture collections all over the word since 1885, mimicking the long-term experimental evolution of E. coli of Lenski et al. (O. Tenaillon, J. E. Barrick, N. Ribeck, D. E. Deatherage, J. L. Blanchard, A. Dasgupta, G. C. Wu, S. Wielgoss, S. Cruveiller, C. Médigue, D. Schneider, and R. E. Lenski, Nature 536:165-170, 2016, https://doi.org/10.1038/nature18959) that shares numerous molecular features.

Analysis of bacterial genomes from an evolution experiment with horizontal gene transfer shows that recombination can sometimes overwhelm selection

Few experimental studies have examined the role that sexual recombination plays in bacterial evolution, including the effects of horizontal gene transfer on genome structure. To address this limitation, we analyzed genomes from an experiment in which Escherichia coli K-12 Hfr (high frequency recombination) donors were periodically introduced into 12 evolving populations of E. coli B and allowed to conjugate repeatedly over the course of 1000 generations. Previous analyses of the evolved strains from this experiment showed that recombination did not accelerate adaptation, despite increasing genetic variation relative to asexual controls. However, the resolution in that previous work was limited to only a few genetic markers. We sought to clarify and understand these puzzling results by sequencing complete genomes from each population. The effects of recombination were highly variable: one lineage was mostly derived from the donors, while another acquired almost no donor DNA. In most lineages, some regions showed repeated introgression and others almost none. Regions with high introgression tended to be near the donors' origin of transfer sites. To determine whether introgressed alleles imposed a genetic load, we extended the experiment for 200 generations without recombination and sequenced whole-population samples. Beneficial alleles in the recipient populations were occasionally driven extinct by maladaptive donor-derived alleles. On balance, our analyses indicate that the plasmid-mediated recombination was sufficiently frequent to drive donor alleles to fixation without providing much, if any, selective advantage.

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.

One Year Genome Evolution of Lausannevirus in Allopatric versus Sympatric Conditions

Amoeba-resisting microorganisms raised a great interest during the last decade. Among them, some large DNA viruses present huge genomes up to 2.5 Mb long, exceeding the size of small bacterial genomes. The rate of genome evolution in terms of mutation, deletion, and gene acquisition in these genomes is yet unknown. Given the suspected high plasticity of viral genomes, the microevolution of the 346 kb genome of Lausannevirus, a member of Megavirales, was studied. Hence, Lausannevirus was co-cultured within the amoeba Acanthamoeba castellanii over one year. Despite a low number of mutations, the virus showed a genome reduction of 3.7% after 12 months. Lausannevirus genome evolution in sympatric conditions was investigated by its co-culture with Estrella lausannensis, an obligate intracellular bacterium, in the amoeba A. castellanii during one year. Cultures were split every 3 months. Genome sequencing revealed that in these conditions both, Lausannevirus and E. lausannensis, show stable genome, presenting no major rearrangement. In fact, after one year they acquired from 2 to 7 and from 4 to 10 mutations per culture for Lausannevirus and E. lausannensis, respectively. Interestingly, different mutations in the endonuclease encoding genes of Lausannevirus were observed in different subcultures, highlighting the importance of this gene product in the replication of Lausannevirus. Conversely, mutations in E. lausannensis were mainly located in a gene encoding for a phosphoenolpyruvate–protein phosphotransferase (PtsI), implicated in sugar metabolism. Moreover, in our conditions and with our analyses we detected no horizontal gene transfer during one year of co-culture.

Rapid Genetic Adaptation during the First Four Months of Survival under Resource Exhaustion

Many bacteria, including the model bacterium Escherichia coli can survive for years within spent media, following resource exhaustion. We carried out evolutionary experiments, followed by whole genome sequencing of hundreds of evolved clones to study the dynamics by which E. coli adapts during the first 4 months of survival under resource exhaustion. Our results reveal that bacteria evolving under resource exhaustion are subject to intense selection, manifesting in rapid mutation accumulation, enrichment in functional mutation categories and extremely convergent adaptation. In the most striking example of convergent adaptation, we found that across five independent populations adaptation to conditions of resource exhaustion occurs through mutations to the three same specific positions of the RNA polymerase core enzyme. Mutations to these three sites are strongly antagonistically pleiotropic, in that they sharply reduce exponential growth rates in fresh media. Such antagonistically pleiotropic mutations, combined with the accumulation of additional mutations, severely reduce the ability of bacteria surviving under resource exhaustion to grow exponentially in fresh media. We further demonstrate that the three positions at which these resource exhaustion mutations occur are conserved for the ancestral E. coli allele, across bacterial phyla, with the exception of nonculturable bacteria that carry the resource exhaustion allele at one of these positions, at very high frequencies. Finally, our results demonstrate that adaptation to resource exhaustion is not limited by mutational input and that bacteria are able to rapidly adapt under resource exhaustion in a temporally precise manner through allele frequency fluctuations.

The impact of bottlenecks on microbial survival, adaptation, and phenotypic switching in host-pathogen interactions

Microbial pathogens and viruses can often maintain sufficient population diversity to evade a wide range of host immune responses. However, when populations experience bottlenecks, as occurs frequently during initiation of new infections, pathogens require specialized mechanisms to regenerate diversity. We address the evolution of such mechanisms, known as stochastic phenotype switches, which are prevalent in pathogenic bacteria. We analyze a model of pathogen diversification in a changing host environment that accounts for selective bottlenecks, wherein different phenotypes have distinct transmission probabilities between hosts. We show that under stringent bottlenecks, such that only one phenotype can initiate new infections, there exists a threshold stochastic switching rate below which all pathogen lineages go extinct, and above which survival is a near certainty. We determine how quickly stochastic switching rates can evolve by computing a fitness landscape for the evolutionary dynamics of switching rates, and analyzing its dependence on both the stringency of bottlenecks and the duration of within-host growth periods. We show that increasing the stringency of bottlenecks or decreasing the period of growth results in faster adaptation of switching rates. Our model provides strong theoretical evidence that bottlenecks play a critical role in accelerating the evolutionary dynamics of pathogens.

The dynamics of molecular evolution over 60,000 generations

The outcomes of evolution are determined by a stochastic dynamical process that governs how mutations arise and spread through a population. However, it is difficult to observe these dynamics directly over long periods and across entire genomes. Here we analyse the dynamics of molecular evolution in twelve experimental populations of Escherichia coli, using whole-genome metagenomic sequencing at five hundred-generation intervals through sixty thousand generations. Although the rate of fitness gain declines over time, molecular evolution is characterized by signatures of rapid adaptation throughout the duration of the experiment, with multiple beneficial variants simultaneously competing for dominance in each population. Interactions between ecological and evolutionary processes play an important role, as long-term quasi-stable coexistence arises spontaneously in most populations, and evolution continues within each clade. We also present evidence that the targets of natural selection change over time, as epistasis and historical contingency alter the strength of selection on different genes. Together, these results show that long-term adaptation to a constant environment can be a more complex and dynamic process than is often assumed.

Mutator genomes decay, despite sustained fitness gains, in a long-term experiment with bacteria

Understanding the extreme variation among bacterial genomes remains an unsolved challenge in evolutionary biology, despite long-standing debate about the relative importance of natural selection, mutation, and random drift. A potentially important confounding factor is the variation in mutation rates between lineages and over evolutionary history, which has been documented in several species. Mutation accumulation experiments have shown that hypermutability can erode genomes over short timescales. These results, however, were obtained under conditions of extremely weak selection, casting doubt on their general relevance. Here, we circumvent this limitation by analyzing genomes from mutator populations that arose during a long-term experiment with Escherichia coli, in which populations have been adaptively evolving for >50,000 generations. We develop an analytical framework to quantify the relative contributions of mutation and selection in shaping genomic characteristics, and we validate it using genomes evolved under regimes of high mutation rates with weak selection (mutation accumulation experiments) and low mutation rates with strong selection (natural isolates). Our results show that, despite sustained adaptive evolution in the long-term experiment, the signature of selection is much weaker than that of mutational biases in mutator genomes. This finding suggests that relatively brief periods of hypermutability can play an outsized role in shaping extant bacterial genomes. Overall, these results highlight the importance of genomic draft, in which strong linkage limits the ability of selection to purge deleterious mutations. These insights are also relevant to other biological systems evolving under strong linkage and high mutation rates, including viruses and cancer cells.

Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii

Pathogenic microbes confront an evolutionary conflict between the pressure to maintain genome stability and the need to adapt to mounting external stresses. Bacteria often respond with elevated mutation rates, but little evidence exists of stable eukaryotic hypermutators in nature. Whole genome resequencing of the human fungal pathogen Cryptococcus deuterogattii identified an outbreak lineage characterized by a nonsense mutation in the mismatch repair component MSH2. This defect results in a moderate mutation rate increase in typical genes, and a larger increase in genes containing homopolymer runs. This allows facile inactivation of genes with coding homopolymer runs including FRR1, which encodes the target of the immunosuppresive antifungal drugs FK506 and rapamycin. Our study identifies a eukaryotic hypermutator lineage spread over two continents and suggests that pathogenic eukaryotic microbes may experience similar selection pressures on mutation rate as bacterial pathogens, particularly during long periods of clonal growth or while expanding into new environments.

As humans, we often think of genetic mutations as being bad. Over the past several decades we have seen health warnings issued on a variety of environmental exposures, from cigarettes to tanning beds, and with good reason because they cause mutations. For multicellular organisms like humans, these mutations are strongly associated with cancer. But in bacteria, this is not true. In fact, the rate at which mutations occur sometimes increases to help bacteria cope with stressful environments.

Unlike bacteria, humans are eukaryotes – the name given to organisms whose cells contain different compartments separated by membranes, such as the nucleus of the cell. For years, we have assumed that eukaryotic microbes, like fungi and parasites, act more like humans than like bacteria because work in budding yeast (another eukaryote) has suggested this to be the case. However, recent work in disease-causing fungi has shown that, much like bacteria, elevated mutation rates may help them to respond to stress. This could also enable fungi to become resistant to drugs used to treat fungal infections.

Cryptococcus deuterogattii is a fungus that causes human diseases including meningoencephalitis and a lung infection called pulmonary cryptococcosis. An ongoing outbreak of the fungus began in the Pacific Northwest of Canada in the late 1990s and emerged in the United States in 2006/2007. Among isolates closely related to those fungi causing the outbreak, three were found that appear to have a specific mutation in their DNA mismatch repair pathway, meaning that they may also experience a higher mutation rate. These strains are also less able to cause disease than others.

Billmyre et al. now demonstrate experimentally that all three isolates have a specific DNA mismatch repair defect, and show that these fungi experience elevated mutation rates, resulting in what is known as a hypermutator state. Furthermore, whole genome sequencing and phylogenetic analysis showed that these hypermutator strains are derived from the outbreak-causing fungi, and that their reduced ability to cause disease is likely a result of accumulating mutations and the loss of the ability to grow at the higher temperatures found in the human body.

Fungal infections are difficult to treat, in part because there are a limited number of available drugs. Elevated mutation rates will likely increase how often and how rapidly fungi develop resistance to these drugs. Understanding how commonly fungi exhibit a hypermutator state that could impact the development of drug resistance will therefore be important for treating patients with fungal infections, which account for millions of infections and hundreds of thousands of deaths annually worldwide.

Fixation probability of a nonmutator in a large population of asexual mutators

In an adapted population of mutators in which most mutations are deleterious, a nonmutator that lowers the mutation rate is under indirect selection and can sweep to fixation. Using a multitype branching process, we calculate the fixation probability of a rare nonmutator in a large population of asexual mutators. We show that when beneficial mutations are absent, the fixation probability is a nonmonotonic function of the mutation rate of the mutator: it first increases sublinearly and then decreases exponentially. We also find that beneficial mutations can enhance the fixation probability of a nonmutator. Our analysis is relevant to an understanding of recent experiments in which a reduction in the mutation rates has been observed.

Genetic variation in adaptability and pleiotropy in budding yeast

Evolution can favor organisms that are more adaptable, provided that genetic variation in adaptability exists. Here, we quantify this variation among 230 offspring of a cross between diverged yeast strains. We measure the adaptability of each offspring genotype, defined as its average rate of adaptation in a specific environmental condition, and analyze the heritability, predictability, and genetic basis of this trait. We find that initial genotype strongly affects adaptability and can alter the genetic basis of future evolution. Initial genotype also affects the pleiotropic consequences of adaptation for fitness in a different environment. This genetic variation in adaptability and pleiotropy is largely determined by initial fitness, according to a rule of declining adaptability with increasing initial fitness, but several individual QTLs also have a significant idiosyncratic role. Our results demonstrate that both adaptability and pleiotropy are complex traits, with extensive heritable differences arising from naturally occurring variation.

Discovering Complete Quasispecies in Bacterial Genomes

Mobile genetic elements can be found in almost all genomes. Possibly the most common nonautonomous mobile genetic elements in bacteria are repetitive extragenic palindromic doublets forming hairpins (REPINs) that can occur hundreds of times within a genome. The sum of all REPINs in a genome can be viewed as an evolving population because REPINs replicate and mutate. In contrast to most other biological populations, we know the exact composition of the REPIN population and the sequence of each member of the population. Here, we model the evolution of REPINs as quasispecies. We fit our quasispecies model to 10 different REPIN populations from 10 different bacterial strains and estimate effective duplication rates. Our estimated duplication rates range from ∼5 × 10-9 to 15 × 10-9 duplications per bacterial generation per REPIN. The small range and the low level of the REPIN duplication rates suggest a universal trade-off between the survival of the REPIN population and the reduction of the mutational load for the host genome. The REPIN populations we investigated also possess features typical of other natural populations. One population shows hallmarks of a population that is going extinct, another population seems to be growing in size, and we also see an example of competition between two REPIN populations.

Coordinated Changes in Mutation and Growth Rates Induced by Genome Reduction

Genome size is determined during evolution, but it can also be altered by genetic engineering in laboratories. The systematic characterization of reduced genomes provides valuable insights into the cellular properties that are quantitatively described by the global parameters related to the dynamics of growth and mutation. In the present study, we analyzed a small collection of W3110 Escherichia coli derivatives containing either the wild-type genome or reduced genomes of various lengths to examine whether the mutation rate, a global parameter representing genomic plasticity, was affected by genome reduction. We found that the mutation rates of these cells increased with genome reduction. The correlation between genome length and mutation rate, which has been reported for the evolution of bacteria, was also identified, intriguingly, for genome reduction. Gene function enrichment analysis indicated that the deletion of many of the genes encoding membrane and transport proteins play a role in the mutation rate changes mediated by genome reduction. Furthermore, the increase in the mutation rate with genome reduction was highly associated with a decrease in the growth rate in a nutrition-dependent manner; thus, poorer media showed a larger change that was of higher significance. This negative correlation was strongly supported by experimental evidence that the serial transfer of the reduced genome improved the growth rate and reduced the mutation rate to a large extent. Taken together, the global parameters corresponding to the genome, growth, and mutation showed a coordinated relationship, which might be an essential working principle for balancing the cellular dynamics appropriate to the environment.

Genome reduction is a powerful approach for investigating the fundamental rules for living systems. Whether genetically disturbed genomes have any specific properties that are different from or similar to those of natively evolved genomes has been under investigation. In the present study, we found that Escherichia coli cells with reduced genomes showed accelerated nucleotide substitution errors (mutation rates), although these cells retained the normal DNA mismatch repair systems. Intriguingly, this finding of correlation between reduced genome size and a higher mutation rate was consistent with the reported evolution of mutation rates. Furthermore, the increased mutation rate was quantitatively associated with a decreased growth rate, indicating that the global parameters related to the genome, growth, and mutation, which represent the amount of genetic information, the efficiency of propagation, and the fidelity of replication, respectively, are dynamically coordinated.

Evolution of organismal stoichiometry in a long-term experiment with Escherichia coli

Organismal stoichiometry refers to the relative proportion of chemical elements in the biomass of organisms, and it can have important effects on ecological interactions from population to ecosystem scales. Although stoichiometry has been studied extensively from an ecological perspective, much less is known about the rates and directions of evolutionary changes in elemental composition. We measured carbon, nitrogen and phosphorus content of 12 Escherichia coli populations that evolved under controlled carbon-limited, serial-transfer conditions for 50 000 generations. The bacteria evolved higher relative nitrogen and phosphorus content, consistent with selection for increased use of the more abundant elements. Total carbon assimilated also increased, indicating more efficient use of the limiting element. We also measured stoichiometry in one population repeatedly through time. Stoichiometry changed more rapidly in early generations than later on, similar to the trajectory seen for competitive fitness. Altogether, our study shows that stoichiometry evolved over long time periods, and that it did so in a predictable direction, given the carbon-limited environment.