Adaptation, Clonal Interference, and Frequency-Dependent Interactions in a Long-Term Evolution Experiment with Escherichia coli

The balance between intrinsic and ecological fitness defines new regimes in eco-evolutionary population dynamics

Selection upon intrinsic fitness differences is one of the most basic mechanisms of evolution, fundamental to all biology. Equally, within macroscopic populations and microscopic environments, ecological interactions influence evolution. Direct experimental evidence of ecological selection between microscopic agents continues to grow. Whilst eco-evolutionary dynamics describes how interactions influence population fitness and composition, we build a model that allows ecological aspects of these interactions to fall on a spectrum independent of the intrinsic fitness of the population. With our mathematical framework, we show how ecological interactions between mutating populations modify the estimated evolutionary trajectories expected from monoculture fitnesses alone. We derive and validate analytical stationary solutions to our partial differential equations that depend on intrinsic and ecological terms, and mutation rates. We determine cases in which these interactions modify evolution in such ways as to, for example, maintain or invert existing monoculture fitness differences. This work discusses the importance of understanding ecological and intrinsic selection effects to avoid misleading conclusions from experiments and defines new ways to assess this balance from experimental results. Using published experimental data, we also show evidence that real microbiological systems can span intrinsic fitness dominant and ecological-effect dominant regimes and that ecological contributions can change with an environment to exaggerate or counteract the composite populations' intrinsic fitness differences.

Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data

In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.

Species and condition shape the mutational spectrum in experimentally evolved biofilms

IMPORTANCE

Biofilm formation is a vital factor for the survival and adaptation of bacteria in diverse environmental niches. Experimental evolution combined with the advancement of whole-population genome sequencing provides us a powerful tool to understand the genomic dynamic of evolutionary adaptation to different environments, such as during biofilm development. Previous studies described the genetic and phenotypic changes of selected clones from experimentally evolved Bacillus thuringiensis and Bacillus subtilis that were adapted under abiotic and biotic biofilm conditions. However, the full understanding of the dynamic evolutionary landscapes was lacking. Furthermore, the differences and similarities of adaptive mechanisms in B. thuringiensis and B. subtilis were not identified. To overcome these limitations, we performed longitudinal whole-population genome sequencing to study the underlying genetic dynamics at high resolution. Our study provides the first comprehensive mutational landscape of two bacterial species' biofilms that is adapted to an abiotic and biotic surface.

Metacommunity structure preserves genome diversity in the presence of gene-specific selective sweeps under moderate rates of horizontal gene transfer

The horizontal transfer of genes is fundamental for the eco-evolutionary dynamics of microbial communities, such as oceanic plankton, soil, and the human microbiome. In the case of an acquired beneficial gene, classic population genetics would predict a genome-wide selective sweep, whereby the genome spreads clonally within the community and together with the beneficial gene, removing genome diversity. Instead, several sources of metagenomic data show the existence of "gene-specific sweeps", whereby a beneficial gene spreads across a bacterial community, maintaining genome diversity. Several hypotheses have been proposed to explain this process, including the decreasing gene flow between ecologically distant populations, frequency-dependent selection from linked deleterious allelles, and very high rates of horizontal gene transfer. Here, we propose an additional possible scenario grounded in eco-evolutionary principles. Specifically, we show by a mathematical model and simulations that a metacommunity where species can occupy multiple patches, acting together with a realistic (moderate) HGT rate, helps maintain genome diversity. Assuming a scenario of patches dominated by single species, our model predicts that diversity only decreases moderately upon the arrival of a new beneficial gene, and that losses in diversity can be quickly restored. We explore the generic behaviour of diversity as a function of three key parameters, frequency of insertion of new beneficial genes, migration rates and horizontal transfer rates.Our results provides a testable explanation for how diversity can be maintained by gene-specific sweeps even in the absence of high horizontal gene transfer rates.

Evolution of a cross-feeding interaction following a key innovation in a long-term evolution experiment with Escherichia coli

The evolution of a novel trait can profoundly change an organism's effects on its environment, which can in turn affect the further evolution of that organism and any coexisting organisms. We examine these effects and feedbacks following the evolution of a novel function in the Long-Term Evolution Experiment (LTEE) with Escherichia coli. A characteristic feature of E. coli is its inability to grow aerobically on citrate (Cit-). Nonetheless, a Cit+ variant with this capacity evolved in one LTEE population after 31 000 generations. The Cit+ clade then coexisted stably with another clade that retained the ancestral Cit- phenotype. This coexistence was shaped by the evolution of a cross-feeding relationship based on C4-dicarboxylic acids, particularly succinate, fumarate, and malate, that the Cit+ variants release into the medium. Both the Cit- and Cit+ cells evolved to grow on these excreted resources. The evolution of aerobic growth on citrate thus led to a transition from an ecosystem based on a single limiting resource, glucose, to one with at least five resources that were either shared or partitioned between the two coexisting clades. Our findings show that evolutionary novelties can change environmental conditions in ways that facilitate diversity by altering ecosystem structure and the evolutionary trajectories of coexisting lineages.

Long-term evolution of proliferating cells using the eVOLVER platform

Experimental evolution using fast-growing unicellular organisms is a unique strategy for deciphering the principles and mechanisms underlying evolutionary processes as well as the architecture and wiring of basic biological functions. Over the past decade, this approach has benefited from the development of powerful systems for the continuous control of the growth of independently evolving cultures. While the first devices compatible with multiplexed experimental evolution remained challenging to implement and required constant user intervention, the recently developed eVOLVER framework represents a fully automated closed-loop system for laboratory evolution assays. However, it remained difficult to maintain and compare parallel evolving cultures in tightly controlled environments over long periods of time using eVOLVER. Furthermore, a number of tools were lacking to cope with the various issues that inevitably occur when conducting such long-term assays. Here we present a significant upgrade of the eVOLVER framework, providing major modifications of the experimental methodology, hardware and software as well as a new stand-alone protocol. Altogether, these adaptations and improvements make the eVOLVER a versatile and unparalleled set-up for long-term experimental evolution.

Revisiting the Design of the Long-Term Evolution Experiment with Escherichia coli

The long-term evolution experiment (LTEE) with Escherichia coli began in 1988 and it continues to this day, with its 12 populations having recently reached 75,000 generations of evolution in a simple, well-controlled environment. The LTEE was designed to explore open-ended questions about the dynamics and repeatability of phenotypic and genetic evolution. Here I discuss various aspects of the LTEE's experimental design that have enabled its stability and success, including the choices of the culture regime, growth medium, ancestral strain, and statistical replication. I also discuss some of the challenges associated with a long-running project, such as handling procedural errors (e.g., cross-contamination) and managing the expanding collection of frozen samples. The simplicity of the experimental design and procedures have supported the long-term stability of the LTEE. That stability-along with the inherent creativity of the evolutionary process and the emergence of new genomic technologies-provides a platform that has allowed talented students and collaborators to pose questions, collect data, and make discoveries that go far beyond anything I could have imagined at the start of the LTEE.

Long-term evolution of proliferating cells using the eVOLVER platform

Experimental evolution using fast-growing unicellular organisms is a unique strategy for deciphering the principles and mechanisms underlying evolutionary processes as well as the architecture and wiring of basic biological functions. Over the past decade, this approach has benefited from the development of powerful systems for the continuous control of the growth of independently evolving cultures. While the first devices compatible with multiplexed experimental evolution remained challenging to implement and required constant user intervention, the recently-developed eVOLVER framework represents a fully automated closed-loop system for laboratory evolution assays. However, it remained difficult to maintain and compare parallel evolving cultures in tightly controlled environments over long periods of time using eVOLVER. Furthermore, a number of tools were lacking to cope with the various issues that inevitably occur when conducting such long-term assays. Here we present a significant upgrade of the eVOLVER framework, providing major modifications of the experimental methodology, hardware and software as well as a new standalone protocol. Altogether, these adaptations and improvements make the eVOLVER a versatile and unparalleled setup for long-term experimental evolution.

Eco-evolutionary modelling of microbial syntrophy indicates the robustness of cross-feeding over cross-facilitation

Syntrophic cooperation among prokaryotes is ubiquitous and diverse. It relies on unilateral or mutual aid that may be both catalytic and metabolic in nature. Hypotheses of eukaryotic origins claim that mitochondrial endosymbiosis emerged from mutually beneficial syntrophy of archaeal and bacterial partners. However, there are no other examples of prokaryotic syntrophy leading to endosymbiosis. One potential reason is that when externalized products become public goods, they incite social conflict due to selfish mutants that may undermine any mutualistic interactions. To rigorously evaluate these arguments, here we construct a general mathematical framework of the ecology and evolution of different types of syntrophic partnerships. We do so both in a general microbial and in a eukaryogenetic context. Studying the case where partners cross-feed on each other's self-inhibiting waste, we show that cooperative partnerships will eventually dominate over selfish mutants. By contrast, systems where producers actively secrete enzymes that cross-facilitate their partners' resource consumption are not robust against cheaters over evolutionary time. We conclude that cross-facilitation is unlikely to provide an adequate syntrophic origin for endosymbiosis, but that cross-feeding mutualisms may indeed have played that role.

Experimental test of the contributions of initial variation and new mutations to adaptive evolution in a novel environment

Experimental evolution is an approach that allows researchers to study organisms as they evolve in controlled environments. Despite the growing popularity of this approach, there are conceptual gaps among projects that use different experimental designs. One such gap concerns the contributions to adaptation of genetic variation present at the start of an experiment and that of new mutations that arise during an experiment. The primary source of genetic variation has historically depended largely on the study organisms. In the long-term evolution experiment (LTEE) using Escherichia coli, for example, each population started from a single haploid cell, and therefore, adaptation depended entirely on new mutations. Most other microbial evolution experiments have followed the same strategy. By contrast, evolution experiments using multicellular, sexually reproducing organisms typically start with preexisting variation that fuels the response to selection. New mutations may also come into play in later generations of these experiments, but it is generally difficult to quantify their contribution in these studies. Here, we performed an experiment using E. coli to compare the contributions of initial genetic variation and new mutations to adaptation in a new environment. Our experiment had four treatments that varied in their starting diversity, with 18 populations in each treatment. One treatment depended entirely on new mutations, while the other three began with mixtures of clones, whole-population samples, or mixtures of whole-population samples from the LTEE. We tracked a genetic marker associated with different founders in two treatments. These data revealed significant variation in fitness among the founders, and that variation impacted evolution in the early generations of our experiment. However, there were no differences in fitness among the treatments after 500 or 2,000 generations in the new environment, despite the variation in fitness among the founders. These results indicate that new mutations quickly dominated, and eventually they contributed more to adaptation than did the initial variation. Our study thus shows that preexisting genetic variation can have a strong impact on early evolution in a new environment, but new beneficial mutations may contribute more to later evolution and can even drive some initially beneficial variants to extinction.

Discovery of positive and purifying selection in metagenomic time series of hypermutator microbial populations

A general method to infer both positive and purifying selection during the real-time evolution of hypermutator pathogens would be broadly useful. To this end, we introduce a Simple Test to Infer Mode of Selection (STIMS) from metagenomic time series of evolving microbial populations. We test STIMS on metagenomic data generated by simulations of bacterial evolution, and on metagenomic data spanning 62,750 generations of Lenski’s long-term evolution experiment with Escherichia coli (LTEE). This benchmarking shows that STIMS detects positive selection in both nonmutator and hypermutator populations, and purifying selection in hypermutator populations. Using STIMS, we find strong evidence of ongoing positive selection on key regulators of the E. coli gene regulatory network, even in some hypermutator populations. STIMS also detects positive selection on regulatory genes in hypermutator populations of Pseudomonas aeruginosa that adapted to subinhibitory concentrations of colistin–an antibiotic of last resort–for just twenty-six days of laboratory evolution. Our results show that the fine-tuning of gene regulatory networks is a general mechanism for rapid and ongoing adaptation. The simplicity of STIMS, together with its intuitive visual interpretation, make it a useful test for positive and purifying selection in metagenomic data sets that track microbial evolution in real-time.

Organisms often evolve elevated mutation rates as they adapt to new environments. That said, the genomic basis of adaptation in asexual hypermutator populations has been difficult to discern, because the few mutations driving adaptation are often associated with large numbers of non-adaptive “hitchhiker” mutations elsewhere in the genome. To address this research gap, we present a Simple Test to Infer Mode of Selection (STIMS) which aggregates the overall signature of selection seen over hundreds of genes. Using STIMS, we find strong positive selection on key regulators of the E. coli gene regulatory network across nonmutator and hypermutator populations of Lenski’s long-term evolution experiment with Escherichia coli. We also find positive selection on regulatory genes in populations of Pseudomonas aeruginosa that evolved hypermutability as they evolved to resist colistin, an antibiotic of last resort. Our work shows that positive and purifying selection can be resolved in hypermutator populations, given sufficient genomic data, and show that regulatory mutations are a major driver of adaptation in experimental populations of microbes.

Idiosyncratic Purifying Selection on Metabolic Enzymes in the Long-Term Evolution Experiment with Escherichia coli

Bacteria, Archaea, and Eukarya all share a common set of metabolic reactions. This implies that the function and topology of central metabolism has been evolving under purifying selection over deep time. Central metabolism may similarly evolve under purifying selection during long-term evolution experiments, although it is unclear how long such experiments would have to run (decades, centuries, millennia) before signs of purifying selection on metabolism appear. I hypothesized that central and superessential metabolic enzymes would show evidence of purifying selection in the long-term evolution experiment with Escherichia coli (LTEE). I also hypothesized that enzymes that specialize on single substrates would show stronger evidence of purifying selection in the LTEE than generalist enzymes that catalyze multiple reactions. I tested these hypotheses by analyzing metagenomic time series covering 62,750 generations of the LTEE. I find mixed support for these hypotheses, because the observed patterns of purifying selection are idiosyncratic and population-specific. To explain this finding, I propose the Jenga hypothesis, named after a children's game in which blocks are removed from a tower until it falls. The Jenga hypothesis postulates that loss-of-function mutations degrade costly, redundant, and non-essential metabolic functions. Replicate populations can therefore follow idiosyncratic trajectories of lost redundancies, despite purifying selection on overall function. I tested the Jenga hypothesis by simulating the evolution of 1,000 minimal genomes under strong purifying selection. As predicted, the minimal genomes converge to different metabolic networks. Strikingly, the core genes common to all 1,000 minimal genomes show consistent signatures of purifying selection in the LTEE.

Evolution of Habitat-Dependent Antibiotic Resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that usually causes difficult-to-treat infections due to its low intrinsic antibiotic susceptibility and outstanding capacity for becoming resistant to antibiotics. In addition, it has a remarkable metabolic versatility, being able to grow in different habitats, from natural niches to different and changing inpatient environments. Study of the environmental conditions that shape genetic and phenotypic changes of P. aeruginosa toward antibiotic resistance supposes a novelty, since experimental evolution assays are usually performed with well-defined antibiotics in regular laboratory growth media. Therefore, in this work we address the extent to which the nutrients’ availability may constrain the evolution of antibiotic resistance. We determined that P. aeruginosa genetic trajectories toward resistance to tobramycin, ceftazidime, and ceftazidime-avibactam are different when evolving in laboratory rich medium, urine, or synthetic sputum. Furthermore, our study, linking genotype with phenotype, showed a clear impact of each analyzed environment on both the fitness and resistance level associated with particular resistance mutations. This indicates that the phenotype associated with specific resistance mutations is variable and dependent on the bacterial metabolic state in each particular habitat. Our results support that the design of evolution-based strategies to tackle P. aeruginosa infections should be based on robust patterns of evolution identified within each particular infection and body location.

IMPORTANCE Predicting evolution toward antibiotic resistance (AR) and its associated trade-offs, such as collateral sensitivity, is important to design evolution-based strategies to tackle AR. However, the effect of nutrients' availability on such evolution, particularly those that can be found under in vivo infection conditions, has been barely addressed. We analyzed the evolutionary patterns of P. aeruginosa in the presence of antibiotics in different media, including urine and synthetic sputum, whose compositions are similar to the ones in infections, finding that AR evolution differs, depending on growth conditions. Furthermore, the representative mutants isolated under each condition tested render different AR levels and fitness costs, depending on nutrients’ availability, supporting the idea that environmental constraints shape the phenotypes associated with specific AR mutations. Consequently, the selection of AR mutations that render similar phenotypes is environment dependent. The analysis of evolution patterns toward AR requires studying growth conditions mimicking those that bacteria face during in vivo evolution.

Overdominant and partially dominant mutations drive clonal adaptation in diploid Saccharomyces cerevisiae

Identification of adaptive targets in experimental evolution typically relies on extensive replication and genetic reconstruction. An alternative approach is to directly assay all mutations in an evolved clone by generating pools of segregants that contain random combinations of evolved mutations. Here, we apply this method to six Saccharomyces cerevisiae clones isolated from four diploid populations that were clonally evolved for 2,000 generations in rich glucose medium. Each clone contains 17-26 mutations relative to the ancestor. We derived intermediate genotypes between the founder and the evolved clones by bulk mating sporulated cultures of the evolved clones to a barcoded haploid version of the ancestor. We competed the resulting barcoded diploids en masse and quantified fitness in the experimental and alternative environments by barcode sequencing. We estimated average fitness effects of evolved mutations using barcode-based fitness assays and whole genome sequencing for a subset of segregants. In contrast to our previous work with haploid evolved clones, we find that diploids carry fewer beneficial mutations, with modest fitness effects (up to 5.4%) in the environment in which they arose. In agreement with theoretical expectations, reconstruction experiments show that all mutations with a detectable fitness effect manifest some degree of dominance over the ancestral allele, and most are overdominant. Genotypes with lower fitness effects in alternative environments allowed us to identify conditions that drive adaptation in our system.

Rapid decline of ceftazidime resistance in antibiotic-free and sub-lethal environments is contingent on genetic background

Trade-offs of antibiotic resistance evolution, such as fitness cost and collateral sensitivity (CS), could be exploited to drive evolution toward antibiotic susceptibility. Decline of resistance may occur when resistance to other drug leads to CS to the first one and when compensatory mutations, or genetic reversion of the original ones, reduce fitness cost. Here we describe the impact of antibiotic-free and sublethal environments on declining ceftazidime resistance in different Pseudomonas aeruginosa resistant mutants. We determined that decline of ceftazidime resistance occurs within 450 generations, which is caused by newly acquired mutations and not by reversion of the original ones, and that the original CS of these mutants is preserved. In addition, we observed that the frequency and degree of this decline is contingent on genetic background. Our results are relevant to implement evolution-based therapeutic approaches, as well as to redefine global policies of antibiotic use, such as drug cycling.

Leapfrog dynamics in phage‐bacteria coevolution revealed by joint analysis of cross‐infection phenotypes and whole genome sequencing

Viruses and their hosts can undergo coevolutionary arms races where hosts evolve increased resistance and viruses evolve counter-resistance. Given these arms race dynamics (ARD), both players are predicted to evolve along a single trajectory as more recently evolved genotypes replace their predecessors. By coupling phenotypic and genomic analyses of coevolving populations of bacteriophage λ and Escherichia coli, we find conflicting evidence for ARD. Virus-host infection phenotypes fit the ARD model, yet genomic analyses revealed fluctuating selection dynamics. Rather than coevolution unfolding along a single trajectory, cryptic genetic variation emerges and is maintained at low frequency for generations until it eventually supplants dominant lineages. These observations suggest a hybrid 'leapfrog' dynamic, revealing weaknesses in the predictive power of standard coevolutionary models. The findings shed light on the mechanisms that structure coevolving ecological networks and reveal the limits of using phenotypic or genomic data alone to differentiate coevolutionary dynamics.

Intra- and inter-host evolution of H9N2 influenza A virus in Japanese quail

Influenza A viruses (IAVs) are constantly evolving. Crucial steps in the infection cycle, such as sialic acid (SA) receptor binding on the host cell surface, can either promote or hamper the emergence of new variants. We previously assessed the relative fitness in Japanese quail of H9N2 variant viruses differing at a single amino acid position, residue 216 in the hemagglutinin (HA) viral surface protein. This site is known to modulate SA recognition. Our prior study generated a valuable set of longitudinal samples from quail transmission groups where the inoculum comprised different mixed populations of HA 216 variant viruses. Here, we leveraged these samples to examine the evolutionary dynamics of viral populations within and between inoculated and naïve contact quails. We found that positive selection dominated HA gene evolution, but fixation of the fittest variant depended on the competition mixture. Analysis of the whole genome revealed further evidence of positive selection acting both within and between hosts. Positive selection drove fixation of variants in non-HA segments within inoculated and contact quails. Importantly, transmission bottlenecks were modulated by the molecular signature at HA 216, revealing viral receptor usage as a determinant of transmitted diversity. Overall, we show that selection strongly shaped the evolutionary dynamics within and between quails. These findings support the notion that selective processes act effectively on IAV populations in poultry hosts, facilitating rapid viral evolution in this ecological niche.

Predatory Bacteria Select for Sustained Prey Diversity

Predator impacts on prey diversity are often studied among higher organisms over short periods, but microbial predator-prey systems allow examination of prey-diversity dynamics over evolutionary timescales. We previously showed that Escherichia coli commonly evolved minority mucoid phenotypes in response to predation by the bacterial predator Myxococcus xanthus by one time point of a coevolution experiment now named MyxoEE-6. Here we examine mucoid frequencies across several MyxoEE-6 timepoints to discriminate between the hypotheses that mucoids were increasing to fixation, stabilizing around equilibrium frequencies, or heading to loss toward the end of MyxoEE-6. In four focal coevolved prey populations, mucoids rose rapidly early in the experiment and then fluctuated within detectable minority frequency ranges through the end of MyxoEE-6, generating frequency dynamics suggestive of negative frequency-dependent selection. However, a competition experiment between mucoid and non-mucoid clones found a predation-specific advantage of the mucoid clone that was insensitive to frequency over the examined range, leaving the mechanism that maintains minority mucoidy unresolved. The advantage of mucoidy under predation was found to be associated with reduced population size after growth (productivity) in the absence of predators, suggesting a tradeoff between productivity and resistance to predation that we hypothesize may reverse mucoid vs non-mucoid fitness ranks within each MyxoEE-6 cycle. We also found that mucoidy was associated with diverse colony phenotypes and diverse candidate mutations primarily localized in the exopolysaccharide operon yjbEFGH. Collectively, our results show that selection from predatory bacteria can generate apparently stable sympatric phenotypic polymorphisms within coevolving prey populations and also allopatric diversity across populations by selecting for diverse mutations and colony phenotypes associated with mucoidy. More broadly, our results suggest that myxobacterial predation increases long-term diversity within natural microbial communities.

High-throughput characterization of mutations in genes that drive clonal evolution using multiplex adaptome capture sequencing

Understanding how cells are likely to evolve can guide medical interventions and bioengineering efforts that must contend with unwanted mutations. The adaptome of a cell-the neighborhood of genetic changes that are most likely to drive adaptation in a given environment-can be mapped by tracking rare beneficial variants during the early stages of clonal evolution. We used multiplex adaptome capture sequencing (mAdCap-seq), a procedure that combines unique molecular identifiers and hybridization-based enrichment, to characterize mutations in eight Escherichia coli genes known to be under selection in a laboratory environment. We tracked 301 mutations at frequencies as low as 0.01% and inferred the fitness effects of 240 of these mutations. There were distinct molecular signatures of selection on protein structure and function for the three genes with the most beneficial mutations. Our results demonstrate how mAdCap-seq can be used to deeply profile a targeted portion of a cell's adaptome.

Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and high-frequency mutation N501Y in different levels of control measures

Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus' genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1,058,020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3% of global COVID-19 cases and determined three types of mutation dynamics: high-frequency, medium-frequency, and low-frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that medium-frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of high-frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with high-frequency or not-high-frequency and both follows similar trends in different levels of stringency.

Multi-omic Characterization of Intraspecies Variation in Laboratory and Natural Environments

Investigation of microbial communities has led to many advances in our understanding of ecosystem function, whether that ecosystem is a subglacial lake or the human gut. Within these communities, much emphasis has been placed on interspecific variation and between-species relationships. However, with current advances in sequencing technology resulting in both the reduction in sequencing costs and the rise of shotgun metagenomic sequencing, the importance of intraspecific variation and within-species relationships is becoming realized. Our group conducts multi-omic analyses to understand how spatial structure and resource availability influence diversification within a species and the potential for long-term coexistence of multiple ecotypes within a microbial community. Here, we present examples of ecotypic variation observed in the lab and in the wild, current challenges faced when investigating intraspecies diversity, and future developments that we expect to define the field over the next 5 years.

Adaptive laboratory evolution of microbial co-cultures for improved metabolite secretion

Adaptive laboratory evolution has proven highly effective for obtaining microorganisms with enhanced capabilities. Yet, this method is inherently restricted to the traits that are positively linked to cell fitness, such as nutrient utilization. Here, we introduce coevolution of obligatory mutualistic communities for improving secretion of fitness‐costly metabolites through natural selection. In this strategy, metabolic cross‐feeding connects secretion of the target metabolite, despite its cost to the secretor, to the survival and proliferation of the entire community. We thus co‐evolved wild‐type lactic acid bacteria and engineered auxotrophic Saccharomyces cerevisiae in a synthetic growth medium leading to bacterial isolates with enhanced secretion of two B‐group vitamins, viz., riboflavin and folate. The increased production was specific to the targeted vitamin, and evident also in milk, a more complex nutrient environment that naturally contains vitamins. Genomic, proteomic and metabolomic analyses of the evolved lactic acid bacteria, in combination with flux balance analysis, showed altered metabolic regulation towards increased supply of the vitamin precursors. Together, our findings demonstrate how microbial metabolism adapts to mutualistic lifestyle through enhanced metabolite exchange.

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR- and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

Experimental-Evolution-Driven Identification of Arabidopsis Rhizosphere Competence Genes in Pseudomonas protegens

Beneficial plant root-associated microorganisms carry out a range of functions that are essential for plant performance. Establishment of a bacterium on plant roots, however, requires overcoming several challenges, including competition with neighboring microorganisms and host immunity. Forward and reverse genetics have led to the identification of mechanisms that are used by beneficial microorganisms to overcome these challenges, such as the production of iron-chelating compounds, the formation of strong biofilms, or the concealment of characteristic microbial molecular patterns that trigger the host immune system. However, how such mechanisms arose from an evolutionary perspective is much less understood. To study bacterial adaptation in the rhizosphere, we employed experimental evolution to track the physiological and genetic dynamics of root-dwelling Pseudomonas protegens in the Arabidopsis thaliana rhizosphere under axenic conditions. This simplified binary one plant/one bacterium system allows for the amplification of key adaptive mechanisms for bacterial rhizosphere colonization. We identified 35 mutations, including single-nucleotide polymorphisms, insertions, and deletions, distributed over 28 genes. We found that mutations in genes encoding global regulators and in genes for siderophore production, cell surface decoration, attachment, and motility accumulated in parallel, underlining the finding that bacterial adaptation to the rhizosphere follows multiple strategies. Notably, we observed that motility increased in parallel across multiple independent evolutionary lines. All together, these results underscore the strength of experimental evolution in identifying key genes, pathways, and processes for bacterial rhizosphere colonization and a methodology for the development of elite beneficial microorganisms with enhanced root-colonizing capacities that can support sustainable agriculture in the future. IMPORTANCE Beneficial root-associated microorganisms carry out many functions that are essential for plant performance. Establishment of a bacterium on plant roots, however, requires overcoming many challenges. Previously, diverse mechanisms that are used by beneficial microorganisms to overcome these challenges were identified. However, how such mechanisms have developed from an evolutionary perspective is much less understood. Here, we employed experimental evolution to track the evolutionary dynamics of a root-dwelling pseudomonad on the root of Arabidopsis. We found that mutations in global regulators, as well as in genes for siderophore production, cell surface decoration, attachment, and motility, accumulate in parallel, emphasizing these strategies for bacterial adaptation to the rhizosphere. We identified 35 mutations distributed over 28 genes. All together, our results demonstrate the power of experimental evolution in identifying key pathways for rhizosphere colonization and a methodology for the development of elite beneficial microorganisms that can support sustainable agriculture.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

Universal Constraints on Protein Evolution in the Long-Term Evolution Experiment with Escherichia coli

Although it is well known that abundant proteins evolve slowly across the tree of life, there is little consensus for why this is true. Here, I report that abundant proteins evolve slowly in the hypermutator populations of Lenski's long-term evolution experiment with Escherichia coli (LTEE). Specifically, the density of all observed mutations per gene, as measured in metagenomic time series covering 60,000 generations of the LTEE, significantly anticorrelates with mRNA abundance, protein abundance, and degree of protein-protein interaction. The same pattern holds for nonsynonymous mutation density. However, synonymous mutation density, measured across the LTEE hypermutator populations, positively correlates with protein abundance. These results show that universal constraints on protein evolution are visible in data spanning three decades of experimental evolution. Therefore, it should be possible to design experiments to answer why abundant proteins evolve slowly.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

Maintenance of Metabolic Plasticity despite Relaxed Selection in a Long-Term Evolution Experiment with Escherichia coli

AbstractTraits that are unused in a given environment are subject to processes that tend to erode them, leading to reduced fitness in other environments. Although this general tendency is clear, we know much less about why some traits are lost while others are retained and about the roles of mutation and selection in generating different responses. We addressed these issues by examining populations of a facultative anaerobe, Escherichia coli, that have evolved for >30 years in the presence of oxygen, with relaxed selection for anaerobic growth and the associated metabolic plasticity. We asked whether evolution led to the loss, improvement, or maintenance of anaerobic growth, and we analyzed gene expression and mutational data sets to understand the outcomes. We identified genomic signatures of both positive and purifying selection on aerobic-specific genes, while anaerobic-specific genes showed clear evidence of relaxed selection. We also found parallel evolution at two interacting loci that regulate anaerobic growth. We competed the ancestor and evolved clones from each population in an anoxic environment, and we found that anaerobic fitness had not decayed, despite relaxed selection. In summary, relaxed selection does not necessarily reduce an organism's fitness in other environments. Instead, the genetic architecture of the traits under relaxed selection and their correlations with traits under positive and purifying selection may sometimes determine evolutionary outcomes.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment

BACKGROUND

Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit.

RESULTS

We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences.

CONCLUSIONS

Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Whose trait is it anyways? Coevolution of joint phenotypes and genetic architecture in mutualisms

Evolutionary biologists typically envision a trait's genetic basis and fitness effects occurring within a single species. However, traits can be determined by and have fitness consequences for interacting species, thus evolving in multiple genomes. This is especially likely in mutualisms, where species exchange fitness benefits and can associate over long periods of time. Partners may experience evolutionary conflict over the value of a multi-genomic trait, but such conflicts may be ameliorated by mutualism's positive fitness feedbacks. Here, we develop a simulation model of a host-microbe mutualism to explore the evolution of a multi-genomic trait. Coevolutionary outcomes depend on whether hosts and microbes have similar or different optimal trait values, strengths of selection and fitness feedbacks. We show that genome-wide association studies can map joint traits to loci in multiple genomes and describe how fitness conflict and fitness feedback generate different multi-genomic architectures with distinct signals around segregating loci. Partner fitnesses can be positively correlated even when partners are in conflict over the value of a multi-genomic trait, and conflict can generate strong mutualistic dependency. While fitness alignment facilitates rapid adaptation to a new optimum, conflict maintains genetic variation and evolvability, with implications for applied microbiome science.

Increased time sampling in an evolve-and-resequence experiment with outcrossing Saccharomyces cerevisiae reveals multiple paths of adaptive change

"Evolve and resequence" (E&R) studies combine experimental evolution and whole-genome sequencing to interrogate the genetics underlying adaptation. Due to ease of handling, E&R work with asexual organisms such as bacteria can employ optimized experimental design, with large experiments and many generations of selection. By contrast, E&R experiments with sexually reproducing organisms are more difficult to implement, and design parameters vary dramatically among studies. Thus, efforts have been made to assess how these differences, such as number of independent replicates, or size of experimental populations, impact inference. We add to this work by investigating the role of time sampling-the number of discrete time points sequence data are collected from evolving populations. Using data from an E&R experiment with outcrossing Saccharomyces cerevisiae in which populations were sequenced 17 times over ~540 generations, we address the following questions: (a) Do more time points improve the ability to identify candidate regions underlying selection? And (b) does high-resolution sampling provide unique insight into evolutionary processes driving adaptation? We find that while time sampling does not improve the ability to identify candidate regions, high-resolution sampling does provide valuable opportunities to characterize evolutionary dynamics. Increased time sampling reveals three distinct trajectories for adaptive alleles: one consistent with classic population genetic theory (i.e., models assuming constant selection coefficients), and two where trajectories suggest more context-dependent responses (i.e., models involving dynamic selection coefficients). We conclude that while time sampling has limited impact on candidate region identification, sampling eight or more time points has clear benefits for studying complex evolutionary dynamics.

Environmental and Physiological Factors Affecting High-Throughput Measurements of Bacterial Growth

Bacterial growth under nutrient-rich and starvation conditions is intrinsically tied to the environmental history and physiological state of the population. While high-throughput technologies have enabled rapid analyses of mutant libraries, technical and biological challenges complicate data collection and interpretation. Here, we present a framework for the execution and analysis of growth measurements with improved accuracy over that of standard approaches. Using this framework, we demonstrate key biological insights that emerge from consideration of culturing conditions and history. We determined that quantification of the background absorbance in each well of a multiwell plate is critical for accurate measurements of maximal growth rate. Using mathematical modeling, we demonstrated that maximal growth rate is dependent on initial cell density, which distorts comparisons across strains with variable lag properties. We established a multiple-passage protocol that alleviates the substantial effects of glycerol on growth in carbon-poor media, and we tracked growth rate-mediated fitness increases observed during a long-term evolution of Escherichia coli in low glucose concentrations. Finally, we showed that growth of Bacillus subtilis in the presence of glycerol induces a long lag in the next passage due to inhibition of a large fraction of the population. Transposon mutagenesis linked this phenotype to the incorporation of glycerol into lipoteichoic acids, revealing a new role for these envelope components in resuming growth after starvation. Together, our investigations underscore the complex physiology of bacteria during bulk passaging and the importance of robust strategies to understand and quantify growth.IMPORTANCE How starved bacteria adapt and multiply under replete nutrient conditions is intimately linked to their history of previous growth, their physiological state, and the surrounding environment. While automated equipment has enabled high-throughput growth measurements, data interpretation and knowledge gaps regarding the determinants of growth kinetics complicate comparisons between strains. Here, we present a framework for growth measurements that improves accuracy and attenuates the effects of growth history. We determined that background absorbance quantification and multiple passaging cycles allow for accurate growth rate measurements even in carbon-poor media, which we used to reveal growth-rate increases during long-term laboratory evolution of Escherichia coli Using mathematical modeling, we showed that maximum growth rate depends on initial cell density. Finally, we demonstrated that growth of Bacillus subtilis with glycerol inhibits the future growth of most of the population, due to lipoteichoic acid synthesis. These studies highlight the challenges of accurate quantification of bacterial growth behaviors.

Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment

Evolutionary innovations allow populations to colonize new ecological niches. We previously reported that aerobic growth on citrate (Cit⁺) evolved in an Escherichia coli population during adaptation to a minimal glucose medium containing citrate (DM25). Cit⁺ variants can also grow in citrate-only medium (DM0), a novel environment for E. coli. To study adaptation to this niche, we founded two sets of Cit⁺ populations and evolved them for 2500 generations in DM0 or DM25. The evolved lineages acquired numerous parallel mutations, many mediated by transposable elements. Several also evolved amplifications of regions containing the maeA gene. Unexpectedly, some evolved populations and clones show apparent declines in fitness. We also found evidence of substantial cell death in Cit⁺ clones. Our results thus demonstrate rapid trait refinement and adaptation to the new citrate niche, while also suggesting a recalcitrant mismatch between E. coli physiology and growth on citrate.

Controlling the Implementation of Transgenic Microbes: Are We Ready for What Synthetic Biology Has to Offer?

Synthetic biology has promised and delivered on an impressive array of applications based on genetically modified microorganisms. While novel biotechnology undoubtedly offers benefits, like all new technology, precautions should be considered during implementation to reduce the risk of both known and unknown adverse effects. To achieve containment of transgenic microorganisms, confidence to a near-scientific certainty that they cannot transfer their transgenic genes to other organisms, and that they cannot survive to propagate in unintended environments, is a priority. Here, we present an in-depth summary of biological containment systems for micro-organisms published to date, including the production of a genetic firewall through genome recoding and physical containment of microbes using auxotrophies, regulation of essential genes, and expression of toxic genes. The level of containment required to consider a transgenic organism suitable for deployment is discussed, as well as standards of practice for developing new containment systems.

Predicting Drug Resistance Using Deep Mutational Scanning

Drug resistance is a major healthcare challenge, resulting in a continuous need to develop new inhibitors. The development of these inhibitors requires an understanding of the mechanisms of resistance for a critical mass of occurrences. Recent genome editing technologies based on high-throughput DNA synthesis and sequencing may help to predict mutations resulting in resistance by testing large mutagenesis libraries. Here we describe the rationale of this approach, with examples and relevance to drug development and resistance in malaria.

Fitness and Productivity Increase with Ecotypic Diversity among Escherichia coli Strains That Coevolved in a Simple, Constant Environment

Polymicrobial consortia occur in both environmental and clinical settings. In many cases, diversity and productivity correlate in these consortia, especially when sustained by positive, density-dependent interactions. However, the evolutionary history of such entities is typically obscure, making it difficult to establish the relative fitness of consortium partners and to use those data to illuminate the diversity-productivity relationship. Here, we dissect an Escherichia coli consortium that evolved under continuous glucose limitation in the laboratory from a single common ancestor. We show that a partnership consisting of cross-feeding ecotypes is better able to secure primary and secondary resources and to convert those resources to offspring than the ancestral clone. Such interactions may be a prelude to a special form of syntrophy and are likely determinants of microbial community structure in nature, including those having clinical significance such as chronic infections.

The productivity of a biological community often correlates with its diversity. In the microbial world this phenomenon can sometimes be explained by positive, density-dependent interactions such as cross-feeding and syntrophy. These metabolic interactions help account for the astonishing variety of microbial life and drive many of the biogeochemical cycles without which life as we know it could not exist. While it is difficult to recapitulate experimentally how these interactions evolved among multiple taxa, we can explore in the laboratory how they arise within one. These experiments provide insight into how different bacterial ecotypes evolve and from these, possibly new “species.” We have previously shown that in a simple, constant environment a single clone of Escherichia coli can give rise to a consortium of genetically and phenotypically differentiated strains, in effect, a set of ecotypes, that coexist by cross-feeding. We marked these different ecotypes and their shared ancestor by integrating fluorescent protein into their genomes and then used flow cytometry to show that each evolved strain is more fit than the shared ancestor, that pairs of evolved strains are fitter still, and that the entire consortium is the fittest of all. We further demonstrate that the rank order of fitness values agrees with estimates of yield, indicating that an experimentally evolved consortium more efficiently converts primary and secondary resources to offspring than its ancestor or any member acting in isolation.

IMPORTANCE Polymicrobial consortia occur in both environmental and clinical settings. In many cases, diversity and productivity correlate in these consortia, especially when sustained by positive, density-dependent interactions. However, the evolutionary history of such entities is typically obscure, making it difficult to establish the relative fitness of consortium partners and to use those data to illuminate the diversity-productivity relationship. Here, we dissect an Escherichia coli consortium that evolved under continuous glucose limitation in the laboratory from a single common ancestor. We show that a partnership consisting of cross-feeding ecotypes is better able to secure primary and secondary resources and to convert those resources to offspring than the ancestral clone. Such interactions may be a prelude to a special form of syntrophy and are likely determinants of microbial community structure in nature, including those having clinical significance such as chronic infections.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Interpreting the Tape of Life: Ancestry-Based Analyses Provide Insights and Intuition about Evolutionary Dynamics

Fine-scale evolutionary dynamics can be challenging to tease out when focused on the broad brush strokes of whole populations over long time spans. We propose a suite of diagnostic analysis techniques that operate on lineages and phylogenies in digital evolution experiments, with the aim of improving our capacity to quantitatively explore the nuances of evolutionary histories in digital evolution experiments. We present three types of lineage measurements: lineage length, mutation accumulation, and phenotypic volatility. Additionally, we suggest the adoption of four phylogeny measurements from biology: phylogenetic richness, phylogenetic divergence, phylogenetic regularity, and depth of the most-recent common ancestor. In addition to quantitative metrics, we also discuss several existing data visualizations that are useful for understanding lineages and phylogenies: state sequence visualizations, fitness landscape overlays, phylogenetic trees, and Muller plots. We examine the behavior of these metrics (with the aid of data visualizations) in two well-studied computational contexts: (1) a set of two-dimensional, real-valued optimization problems under a range of mutation rates and selection strengths, and (2) a set of qualitatively different environments in the Avida digital evolution platform. These results confirm our intuition about how these metrics respond to various evolutionary conditions and indicate their broad value.

Increasing growth rate slows adaptation when genotypes compete for diffusing resources

The rate at which a species responds to natural selection is a central predictor of the species’ ability to adapt to environmental change. It is well-known that spatially-structured environments slow the rate of adaptation due to increased intra-genotype competition. Here, we show that this effect magnifies over time as a species becomes better adapted and grows faster. Using a reaction-diffusion model, we demonstrate that growth rates are inextricably coupled with effective spatial scales, such that higher growth rates cause more localized competition. This has two effects: selection requires more generations for beneficial mutations to fix, and spatially-caused genetic drift increases. Together, these effects diminish the value of additional growth rate mutations in structured environments.

What determines how quickly a beneficial mutation will spread through a population? The intuitive answer is that mutations that confer faster growth rates will spread at a rate that is relative to the size of the growth-rate benefit. Indeed, this is true in a well-mixed environment where all genotypes compete globally. But most organisms don’t live in a simple well-mixed environment. Many organisms, like bacteria, live in a structured environment, such as on the surface of a solid substrate. Does life on a surface change the expectation about the spread of faster-growing mutants? We developed a mathematical model to answer this question, and found that on a surface, the actual growth rates—not just the relative growth rates—were critical to determining how fast a faster-growing mutant spread through a population. When the simulated organisms grew slowly, competition was basically global and a faster-growing mutant could pre-empt resources from far-away competitors. In contrast, when organisms grew more quickly, competition became much more localized, and the faster-growing mutant could only steal resources from neighboring competitors. This result means that there are diminishing returns to series of mutations which confer growth-rate benefits. This idea will help us predict and understand future and past evolutionary trajectories.

Antibiotic Resistance Evolution Is Contingent on the Quorum-Sensing Response in Pseudomonas aeruginosa

Different works have explored independently the evolution toward antibiotic resistance and the role of eco-adaptive mutations in the adaptation to a new habitat (as the infected host) of bacterial pathogens. However, knowledge about the connection between both processes is still limited. We address this issue by comparing the evolutionary trajectories toward antibiotic resistance of a Pseudomonas aeruginosa lasR defective mutant and its parental wild-type strain, when growing in presence of two ribosome-targeting antibiotics. Quorum-sensing lasR defective mutants are selected in P. aeruginosa populations causing chronic infections. Further, we observed they are also selected in vitro as a first adaptation for growing in culture medium. By using experimental evolution and whole-genome sequencing, we found that the evolutionary trajectories of P. aeruginosa in presence of these antibiotics are different in lasR defective and in wild-type backgrounds, both at the phenotypic and the genotypic levels. Recreation of a set of mutants in both genomic backgrounds (either wild type or lasR defective) allowed us to determine the existence of negative epistatic interactions between lasR and antibiotic resistance determinants. These epistatic interactions could lead to mutual contingency in the evolution of antibiotic resistance when P. aeruginosa colonizes a new habitat in presence of antibiotics. If lasR mutants are selected first, this would constraint antibiotic resistance evolution. Conversely, when resistance mutations (at least those studied in the present work) are selected, lasR mutants may not be selected in presence of antibiotics. These results underlie the importance of contingency and epistatic interactions in modulating antibiotic resistance evolution.

Evolutionary Rescue from a Wave of Biological Invasion

Evolution can potentially rescue populations from being driven extinct by biological invasions, but predictions for this occurrence are generally lacking. Here I derive theoretical predictions for evolutionary rescue of a resident population experiencing invasion from an introduced competitor that spreads over its introduced range as a traveling spatial wave that displaces residents. I compare the likelihood of evolutionary rescue from invasion for two modes of competition: exploitation and interference competition. I find that, all else equal, evolutionary rescue is less effective at preventing extinction caused by interference-driven invasions than by exploitation-driven invasions. Rescue from interference-driven invasions is, surprisingly, independent of invader dispersal rate or the speed of invasion and is more weakly dependent on range size than in the exploitation-driven case. In contrast, rescue from exploitation-driven invasions strongly depends on range size and is less likely during fast invasions. The results presented here have potential applications for conserving endemic species from nonnative invaders and for ensuring extinction of pests using intentionally introduced biocontrol agents.

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.

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.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

Designing Automated, High-throughput, Continuous Cell Growth Experiments Using eVOLVER

Continuous culture methods enable cells to be grown under quantitatively controlled environmental conditions, and are thus broadly useful for measuring fitness phenotypes and improving our understanding of how genotypes are shaped by selection. Extensive recent efforts to develop and apply niche continuous culture devices have revealed the benefits of conducting new forms of cell culture control. This includes defining custom selection pressures and increasing throughput for studies ranging from long-term experimental evolution to genome-wide library selections and synthetic gene circuit characterization. The eVOLVER platform was recently developed to meet this growing demand: a continuous culture platform with a high degree of scalability, flexibility, and automation. eVOLVER provides a single standardizing platform that can be (re)-configured and scaled with minimal effort to perform many different types of high-throughput or multi-dimensional growth selection experiments. Here, a protocol is presented to provide users of the eVOLVER framework a description for configuring the system to conduct a custom, large-scale continuous growth experiment. Specifically, the protocol guides users on how to program the system to multiplex two selection pressures - temperature and osmolarity - across many eVOLVER vials in order to quantify fitness landscapes of Saccharomyces cerevisiae mutants at fine resolution. We show how the device can be configured both programmatically, through its open-source web-based software, and physically, by arranging fluidic and hardware layouts. The process of physically setting up the device, programming the culture routine, monitoring and interacting with the experiment in real-time over the internet, sampling vials for subsequent offline analysis, and post experiment data analysis are detailed. This should serve as a starting point for researchers across diverse disciplines to apply eVOLVER in the design of their own complex and high-throughput cell growth experiments to study and manipulate biological systems.

Nonlinear frequency-dependent selection promotes long-term coexistence between bacteria species

Negative frequency-dependent selection (NFDS) is an important mechanism for species coexistence and for the maintenance of genetic polymorphism. Long-term coexistence nevertheless requires NFDS interactions to be resilient to further evolution of the interacting species or genotypes. For closely related genotypes, NFDS interactions have been shown to be preserved through successive rounds of evolution in coexisting lineages. On the contrary, the evolution of NFDS interactions between distantly related species has received less attention. Here, we tracked the co-evolution of Escherichia coli and Citrobacter freundii that initially differ in their ecological characteristics. We showed that these two bacterial species engaged in an NFDS interaction particularly resilient to further evolution: despite a very strong asymmetric rate of adaptation, their coexistence was maintained owing to an NFDS pattern where fitness increases steeply as the frequency decreases towards zero. Using a model, we showed how and why such NFDS pattern can emerge. These findings provide a robust explanation for the long-term maintenance of species at very low frequencies.

Fibroblasts and alectinib switch the evolutionary games played by non-small cell lung cancer

Heterogeneity in strategies for survival and proliferation among the cells which constitute a tumour is a driving force behind the evolution of resistance to cancer therapy. The rules mapping the tumour’s strategy distribution to the fitness of individual strategies can be represented as an evolutionary game. We develop a game assay to measure effective evolutionary games in co-cultures of non-small cell lung cancer cells which are sensitive and resistant to the anaplastic lymphoma kinase inhibitor Alectinib. The games are not only quantitatively different between different environments, but targeted therapy and cancer associated fibroblasts qualitatively switch the type of game being played by the in-vitro population from Leader to Deadlock. This observation provides empirical confirmation of a central theoretical postulate of evolutionary game theory in oncology: we can treat not only the player, but also the game. Although we concentrate on measuring games played by cancer cells, the measurement methodology we develop can be used to advance the study of games in other microscopic systems by providing a quantitative description of non-cell-autonomous effects.

Multiscale Evolutionary Dynamics of Host-Associated Microbiomes

The composite members of the microbiota face a range of selective pressures and must adapt to persist in the host. We highlight recent work characterizing the evolution and transfer of genetic information across nested scales of host-associated microbiota, which enable resilience to biotic and abiotic perturbations. At the strain level, we consider the preservation and diversification of adaptive information in progeny lineages. At the community level, we consider genetic exchange between distinct microbes in the ecosystem. Finally, we frame microbiomes as open systems subject to acquisition of novel information from foreign ecosystems through invasion by outsider microbes.