Measuring selection coefficients below 10(-3): method, questions, and prospects

Linking selection to demography in experimental evolution of active death in a unicellular organism

Deciphering how natural selection emerges from demographic differences among genotypes, and reciprocally how evolution affects population dynamics, is key to understanding population responses to environmental stress. This is especially true in non-trivial ecological scenarios, such as programmed cell death (PCD) in unicellular organisms, which can lead to massive population decline in response to stress. To understand how selection may operate on this trait, we exposed monocultures and mixtures of two closely related strains of the microalga Dunalliela salina, one of which induces PCD, to multiple cycles of hyper-osmotic shocks, and tracked demography and selection throughout. Population dynamics were consistent between mixtures and monocultures, suggesting that selection on PCD does not involve strong ecological interactions. The PCD-inducing strain was maintained throughout the experiment despite an initial decline, by a combination of fast population rebound following each decline, and density-dependent competition dynamics near the stationary phase that were independent of these initial population fluctuations. As a result of PCD maintenance, population decline in response to environmental stress was not counter-selected in our experiment, but persisted over 13 cycles of salinity. Our results highlight how analysing the demographic underpinnings of fitness and competition can shed light on the mechanisms underlying selection and eco-evolutionary dynamics.

Reversions mask the contribution of adaptive evolution in microbiomes

When examining bacterial genomes for evidence of past selection, the results depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by dN/dS, the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with opposite implications for dynamical intuition and applications of dN/dS. Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the dN/dS decay given only dozens of locally fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of dN/dS obtained from long timescales with caution as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short timescales.

Reversions mask the contribution of adaptive evolution in microbiomes

When examining bacterial genomes for evidence of past selection, the results obtained depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed byd N / d S , the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale-dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with exactly opposite implications for dynamical intuition and applications ofd N / d S . Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain thed N / d S decay given only dozens of locally-fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values ofd N / d S obtained from long-time scales with caution, as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short time scales.

Temperature dependence of the mutation rate towards antibiotic resistance

Objectives

Environmental conditions can influence mutation rates in bacteria. Fever is a common response to infection that alters the growth conditions of infecting bacteria. Here we examine how a temperature change, such as is associated with fever, affects the mutation rate towards antibiotic resistance.

Methods

We used a fluctuation test to assess the mutation rate towards antibiotic resistance in Escherichia coli at two different temperatures: 37°C (normal temperature) and 40°C (fever temperature). We performed this measurement for three different antibiotics with different modes of action: ciprofloxacin, rifampicin and ampicillin.

Results

In all cases, the mutation rate towards antibiotic resistance turned out to be temperature dependent, but in different ways. Fever temperatures led to a reduced mutation rate towards ampicillin resistance and an elevated mutation rate towards ciprofloxacin and rifampicin resistance.

Conclusions

This study shows that the mutation rate towards antibiotic resistance is impacted by a small change in temperature, such as associated with fever. This opens a new avenue to mitigate the emergence of antibiotic resistance by coordinating the choice of an antibiotic with the decision of whether or not to suppress fever when treating a patient. Hence, optimized combinations of antibiotics and fever suppression strategies may be a new weapon in the battle against antibiotic resistance.

Approximate Bayesian computational methods to estimate the strength of divergent selection in population genomics models

Statistical estimation of parameters in large models of evolutionary processes is often too computationally inefficient to pursue using exact model likelihoods, even with single-nucleotide polymorphism (SNP) data, which offers a way to reduce the size of genetic data while retaining relevant information. Approximate Bayesian Computation (ABC) to perform statistical inference about parameters of large models takes the advantage of simulations to bypass direct evaluation of model likelihoods. We develop a mechanistic model to simulate forward-in-time divergent selection with variable migration rates, modes of reproduction (sexual, asexual), length and number of migration-selection cycles. We investigate the computational feasibility of ABC to perform statistical inference and study the quality of estimates on the position of loci under selection and the strength of selection. To expand the parameter space of positions under selection, we enhance the model by implementing an outlier scan on summarized observed data. We evaluate the usefulness of summary statistics well-known to capture the strength of selection, and assess their informativeness under divergent selection. We also evaluate the effect of genetic drift with respect to an idealized deterministic model with single-locus selection. We discuss the role of the recombination rate as a confounding factor in estimating the strength of divergent selection, and emphasize its importance in break down of linkage disequilibrium (LD). We answer the question for which part of the parameter space of the model we recover strong signal for estimating the selection, and determine whether population differentiation-based summary statistics or LD-based summary statistics perform well in estimating selection.

Measuring differential fitness costs and interactions between genetic cassettes using fluorescent spectrophotometry

In this article, we present a method for designing, executing, and analyzing data from a microbial competition experiment. We use fluorescent reporters to label different competing strains and resolve individual growth curves using a fluorescent spectrophotometer. Our comprehensive data analysis pipeline integrates multiple experiments to simultaneously infer sources of variation, extract selection coefficients, and estimate the genetic contributions to fitness for various synthetic genetic cassettes (SGCs). To demonstrate the method, we employ a synthetic biological system based on Escherichia coli. Strains carry 1 of 10 different plasmids and one of three genomically integrated fluorescent markers. All strains are co-cultured to obtain real-time measurements of optical density (total population density) and fluorescence (sub-population densities). We identify challenges in calibrating between fluorescence and density and of fluorescent proteins maturing at different rates. To resolve these issues, we compare two methods of fluorescence calibration and correct for maturation by measuring in vivo maturation times. We provide evidence of genetic interactions occurring between our SGCs and further show how to use our statistical model to test some hypotheses about microbial growth and the costs of protein expression.IMPORTANCEFluorescently labeled co-cultures are becoming increasingly popular. The approach proposed here offers a high standard for experimental design and data analysis to measure selection coefficients and growth rates in competition. Measuring competitive differences is useful in many laboratory studies, allowing for fitness cost-correction of growth rates and ecological interactions and testing hypotheses in synthetic biology. Using time-resolved growth curves, rather than endpoint measurements, for competition assays allows us to construct a detailed scientific model that can be used to ask questions about fine-grained phenomena, such as bacterial growth dynamics, as well as higher-level phenomena, such as the interactions between synthetic cassette expression.

Predicting bacterial interaction outcomes from monoculture growth and supernatant assays

How to derive principles of community dynamics and stability is a central question in microbial ecology. Bottom-up experiments, in which a small number of bacterial species are mixed, have become popular to address it. However, experimental setups are typically limited because co-culture experiments are labor-intensive and species are difficult to distinguish. Here, we use a four-species bacterial community to show that information from monoculture growth and inhibitory effects induced by secreted compounds can be combined to predict the competitive rank order in the community. Specifically, integrative monoculture growth parameters allow building a preliminary competitive rank order, which is then adjusted using inhibitory effects from supernatant assays. While our procedure worked for two different media, we observed differences in species rank orders between media. We then parameterized computer simulations with our empirical data to show that higher order species interactions largely follow the dynamics predicted from pairwise interactions with one important exception. The impact of inhibitory compounds was reduced in higher order communities because their negative effects were spread across multiple target species. Altogether, we formulated three simple rules of how monoculture growth and supernatant assay data can be combined to establish a competitive species rank order in an experimental four-species community.

Schrödinger's yeast: the challenge of using transformation to compare fitness among Saccharomyces cerevisiae that differ in ploidy or zygosity

How the number of genome copies modifies the effect of random mutations remains poorly known. In yeast, researchers have investigated these effects for knock-out or other large-effect mutations, but have not accounted for differences at the mating-type locus. We set out to compare fitness differences among strains that differ in ploidy and/or zygosity using a panel of spontaneously arising mutations acquired in haploid yeast from a previous study. To ensure no genetic differences, even at the mating-type locus, we embarked on a series of transformations, which first sterilized and then temporarily introduced plasmid-borne mating types. Despite these attempts to equalize the haplotypes, fitness variation introduced during transformation swamped the differences among the original mutation-accumulation lines. While colony size looked normal, we observed a bi-modality in the maximum growth rate of our transformed yeast and determined that many of the slow growing lines were respiratory deficient ("petite"). Not previously reported, we found that yeast that were TID1/RDH54 knockouts were less likely to become petite. Even for lines with the same petite status, however, we found no correlation in fitness between the two replicate transformations performed. These results pose a challenge for any study using transformation to measure the fitness effect of genetic differences among strains. By attempting to hold haplotypes constant, we introduced more mutations that overwhelmed our ability to measure fitness differences between the genetic states. In this study, we transformed over one hundred different lines of yeast, using two independent transformations, and found that this common laboratory procedure can cause large changes to the microbe studied. Our study provides a cautionary tale of the need to use multiple transformants in fitness assays.

Phenotype Design Space Provides a Mechanistic Framework Relating Molecular Parameters to Phenotype Diversity Available for Selection

Two long-standing challenges in theoretical population genetics and evolution are predicting the distribution of phenotype diversity generated by mutation and available for selection, and determining the interaction of mutation, selection and drift to characterize evolutionary equilibria and dynamics. More fundamental for enabling such predictions is the current inability to causally link genotype to phenotype. There are three major mechanistic mappings required for such a linking - genetic sequence to kinetic parameters of the molecular processes, kinetic parameters to biochemical system phenotypes, and biochemical phenotypes to organismal phenotypes. This article introduces a theoretical framework, the Phenotype Design Space (PDS) framework, for addressing these challenges by focusing on the mapping of kinetic parameters to biochemical system phenotypes. It provides a quantitative theory whose key features include (1) a mathematically rigorous definition of phenotype based on biochemical kinetics, (2) enumeration of the full phenotypic repertoire, and (3) functional characterization of each phenotype independent of its context-dependent selection or fitness contributions. This framework is built on Design Space methods that relate system phenotypes to genetically determined parameters and environmentally determined variables. It also has the potential to automate prediction of phenotype-specific mutation rate constants and equilibrium distributions of phenotype diversity in microbial populations undergoing steady-state exponential growth, which provides an ideal reference to which more realistic cases can be compared. Although the framework is quite general and flexible, the details will undoubtedly differ for different functions, organisms and contexts. Here a hypothetical case study involving a small molecular system, a primordial circadian clock, is used to introduce this framework and to illustrate its use in a particular case. The framework is built on fundamental biochemical kinetics. Thus, the foundation is based on linear algebra and reasonable physical assumptions, which provide numerous opportunities for experimental testing and further elaboration to deal with complex multicellular organisms that are currently beyond its scope. The discussion provides a comparison of results from the PDS framework with those from other approaches in theoretical population genetics.

Extreme Sensitivity of Fitness to Environmental Conditions: Lessons from #1BigBatch

The phrase "survival of the fittest" has become an iconic descriptor of how natural selection works. And yet, precisely measuring fitness, even for single-celled microbial populations growing in controlled laboratory conditions, remains a challenge. While numerous methods exist to perform these measurements, including recently developed methods utilizing DNA barcodes, all methods are limited in their precision to differentiate strains with small fitness differences. In this study, we rule out some major sources of imprecision, but still find that fitness measurements vary substantially from replicate to replicate. Our data suggest that very subtle and difficult to avoid environmental differences between replicates create systematic variation across fitness measurements. We conclude by discussing how fitness measurements should be interpreted given their extreme environment dependence. This work was inspired by the scientific community who followed us and gave us tips as we live tweeted a high-replicate fitness measurement experiment at #1BigBatch.

Distribution of fitness effects of cross-species transformation reveals potential for fast adaptive evolution

Bacterial transformation, a common mechanism of horizontal gene transfer, can speed up adaptive evolution. How its costs and benefits depend on the growth environment is poorly understood. Here, we characterize the distributions of fitness effects (DFE) of transformation in different conditions and test whether they predict in which condition transformation is beneficial. To determine the DFEs, we generate hybrid libraries between the recipient Bacillus subtilis and different donor species and measure the selection coefficient of each hybrid strain. In complex medium, the donor Bacillus vallismortis confers larger fitness effects than the more closely related donor Bacillus spizizenii. For both donors, the DFEs show strong effect beneficial transfers, indicating potential for fast adaptive evolution. While some transfers of B. vallismortis DNA show pleiotropic effects, various transfers are beneficial only under a single growth condition, indicating that the recipient can benefit from a variety of donor genes to adapt to varying growth conditions. We scrutinize the predictive value of the DFEs by laboratory evolution under different growth conditions and show that the DFEs correctly predict the condition at which transformation confers a benefit. We conclude that transformation has a strong potential for speeding up adaptation to varying environments by profiting from a gene pool shared between closely related species.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

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

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

Precise measurement of the fitness effects of spontaneous mutations by droplet digital PCR in Burkholderia cenocepacia

Understanding how mutations affect survivability is a key component to knowing how organisms and complex traits evolve. However, most mutations have a minor effect on fitness and these effects are difficult to resolve using traditional molecular techniques. Therefore, there is a dire need for more accurate and precise fitness measurements methods. Here, we measured the fitness effects in Burkholderia cenocepacia HI2424 mutation accumulation (MA) lines using droplet-digital polymerase chain reaction (ddPCR). Overall, the fitness measurements from ddPCR-MA are correlated positively with fitness measurements derived from traditional phenotypic marker assays (r = 0.297, P = 0.05), but showed some differences. First, ddPCR had significantly lower measurement variance in fitness (F = 3.78, P < 2.6 × 10-13) in control experiments. Second, the mean fitness from ddPCR-MA measurements were significantly lower than phenotypic marker assays (-0.0041 vs -0.0071, P = 0.006). Consistent with phenotypic marker assays, ddPCR-MA measurements observed multiple (27/43) lineages that significantly deviated from mean fitness, suggesting that a majority of the mutations are neutral or slightly deleterious and intermixed with a few mutations that have extremely large effects. Of these mutations, we found a significant excess of mutations within DNA excinuclease and Lys R transcriptional regulators that have extreme deleterious and beneficial effects, indicating that modifications to transcription and replication may have a strong effect on organismal fitness. This study demonstrates the power of ddPCR as a ubiquitous method for high-throughput fitness measurements in both DNA- and RNA-based organisms regardless of cell type or physiology.

Plasmid fitness costs are caused by specific genetic conflicts enabling resolution by compensatory mutation

Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid carriage are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or gene expression level. By combining the results of experimental evolution with genetics and transcriptomics, we show here that fitness costs of 2 divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid carriage. We identify one of these compensatory loci, the chromosomal gene PFLU4242, as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by up-regulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained up-regulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer (HGT) due to their propensity for amelioration by single compensatory mutations, helping to explain why plasmids are so common in bacterial genomes.

Predicting population genetic change in an autocorrelated random environment: Insights from a large automated experiment

Most natural environments exhibit a substantial component of random variation, with a degree of temporal autocorrelation that defines the color of environmental noise. Such environmental fluctuations cause random fluctuations in natural selection, affecting the predictability of evolution. But despite long-standing theoretical interest in population genetics in stochastic environments, there is a dearth of empirical estimation of underlying parameters of this theory. More importantly, it is still an open question whether evolution in fluctuating environments can be predicted indirectly using simpler measures, which combine environmental time series with population estimates in constant environments. Here we address these questions by using an automated experimental evolution approach. We used a liquid-handling robot to expose over a hundred lines of the micro-alga Dunaliella salina to randomly fluctuating salinity over a continuous range, with controlled mean, variance, and autocorrelation. We then tracked the frequencies of two competing strains through amplicon sequencing of nuclear and choloroplastic barcode sequences. We show that the magnitude of environmental fluctuations (determined by their variance), but also their predictability (determined by their autocorrelation), had large impacts on the average selection coefficient. The variance in frequency change, which quantifies randomness in population genetics, was substantially higher in a fluctuating environment. The reaction norm of selection coefficients against constant salinity yielded accurate predictions for the mean selection coefficient in a fluctuating environment. This selection reaction norm was in turn well predicted by environmental tolerance curves, with population growth rate against salinity. However, both the selection reaction norm and tolerance curves underestimated the variance in selection caused by random environmental fluctuations. Overall, our results provide exceptional insights into the prospects for understanding and predicting genetic evolution in randomly fluctuating environments.

Being able to predict evolution under natural selection is important for many applied fields of biology, ranging from agriculture to medicine or conservation. However, this endeavor is complicated by factors that inherently limit our ability to predict the future, such as random fluctuations in the environment. Population genetic theory indicates that probabilistic predictions can still be made in this context, but the extent to which this holds empirically, and whether these predictions can be based on simple measurements, are still open questions. Making progress on answering these questions can be achieved by capitalizing on experiments where the environment is precisely controlled over many generations. Here, we used a pipetting robot to generate random time series of salinities with controlled patterns of fluctuations, which we imposed on a microalga, Dunaliella salina. Tracking the frequencies of two genotypes in a mixture by sequencing two short barcode sequences, we were able to show how patterns of fluctuating selection relate to the fluctuating environment. Interestingly, parts of these responses, but not all, could be predicted by simpler measurements in constant environments, allowing precise characterization the limits and prospects for predicting evolution in fluctuating environments.

Are Nonsynonymous Transversions Generally More Deleterious than Nonsynonymous Transitions?

It has been suggested that, due to the structure of the genetic code, nonsynonymous transitions are less likely than transversions to cause radical changes in amino acid physicochemical properties so are on average less deleterious. This view was supported by some but not all mutagenesis experiments. Because laboratory measures of fitness effects have limited sensitivities and relative frequencies of different mutations in mutagenesis studies may not match those in nature, we here revisit this issue using comparative genomics. We extend the standard codon model of sequence evolution by adding the parameter η that quantifies the ratio of the fixation probability of transitional nonsynonymous mutations to that of transversional nonsynonymous mutations. We then estimate η from the concatenated alignment of all protein-coding DNA sequences of two closely related genomes. Surprisingly, η ranges from 0.13 to 2.0 across 90 species pairs sampled from the tree of life, with 51 incidences of η < 1 and 30 incidences of η >1 that are statistically significant. Hence, whether nonsynonymous transversions are overall more deleterious than nonsynonymous transitions is species-dependent. Because the corresponding groups of amino acid replacements differ between nonsynonymous transitions and transversions, η is influenced by the relative exchangeabilities of amino acid pairs. Indeed, an extensive search reveals that the large variation in η is primarily explainable by the recently reported among-species disparity in amino acid exchangeabilities. These findings demonstrate that genome-wide nucleotide substitution patterns in coding sequences have species-specific features and are more variable among evolutionary lineages than are currently thought.

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

High-throughput analysis of adaptation using barcoded strains of Saccharomyces cerevisiae

BACKGROUND

Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge.

RESULTS

We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics.

CONCLUSIONS

Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Selective Sweep at a QTL in a Randomly Fluctuating Environment

Adaptation is mediated by phenotypic traits that are often near continuous, and undergo selective pressures that may change with the environment. The dynamics of allelic frequencies at underlying quantitative trait loci (QTL) depend on their own phenotypic effects, but also possibly on other polymorphic loci affecting the same trait, and on environmental change driving phenotypic selection. Most environments include a substantial component of random noise, characterized both by its magnitude and its temporal autocorrelation, which sets the timescale of environmental predictability. I investigate the dynamics of a mutation affecting a quantitative trait in an autocorrelated stochastic environment that causes random fluctuations of an optimum phenotype. The trait under selection may also exhibit background polygenic variance caused by many polymorphic loci of small effects elsewhere in the genome. In addition, the mutation at the QTL may affect phenotypic plasticity, the phenotypic response of given genotype to its environment of development or expression. Stochastic environmental fluctuations increase the variance of the evolutionary process, with consequences for the probability of a complete sweep at the QTL. Background polygenic variation critically alters this process, by setting an upper limit to stochastic variance of population genetics at the QTL. For a plasticity QTL, stochastic fluctuations also influences the expected selection coefficient, and alleles with the same expected trajectory can have very different stochastic variances. Finally, a mutation may be favored through its effect on plasticity despite causing a systematic mismatch with optimum, which is compensated by evolution of the mean background phenotype.

The fitness cost of mismatch repair mutators in Saccharomyces cerevisiae: partitioning the mutational load

Mutational load is the depression in a population's mean fitness that results from the continual influx of deleterious mutations. Here, we directly estimate the mutational load in a population of haploid Saccharomyces cerevisiae that are deficient for mismatch repair. We partition the load in haploids into two components. To estimate the load due to nonlethal mutations, we measure the competitive fitness of hundreds of randomly selected clones from both mismatch-repair-deficient and -proficient populations. Computation of the mean clone fitness for the mismatch-repair-deficient strain permits an estimation of the nonlethal load, and the histogram of fitness provides an interesting visualization of a loaded population. In a separate experiment, in order to estimate the load due to lethal mutations (i.e. the lethal mutation rate), we manipulate thousands of individual pairs of mother and daughter cells and track their fates. These two approaches yield point estimates for the two contributors to load, and the addition of these estimates is nearly equal to the separately measured short-term competitive fitness deficit for the mismatch-repair-deficient strain. This correspondence suggests that there is no need to invoke direct fitness effects to explain the fitness difference between mismatch-repair-deficient and -proficient strains. Assays in diploids are consistent with deleterious mutations in diploids tending towards recessivity. These results enhance our understanding of mutational load, a central population genetics concept, and we discuss their implications for the evolution of mutation rates.

Predicting microbial growth in a mixed culture from growth curve data

Determining the fitness of specific microbial genotypes has extensive application in microbial genetics, evolution, and biotechnology. While estimates from growth curves are simple and allow high throughput, they are inaccurate and do not account for interactions between costs and benefits accruing over different parts of a growth cycle. For this reason, pairwise competition experiments are the current "gold standard" for accurate estimation of fitness. However, competition experiments require distinct markers, making them difficult to perform between isolates derived from a common ancestor or between isolates of nonmodel organisms. In addition, competition experiments require that competing strains be grown in the same environment, so they cannot be used to infer the fitness consequence of different environmental perturbations on the same genotype. Finally, competition experiments typically consider only the end-points of a period of competition so that they do not readily provide information on the growth differences that underlie competitive ability. Here, we describe a computational approach for predicting density-dependent microbial growth in a mixed culture utilizing data from monoculture and mixed-culture growth curves. We validate this approach using 2 different experiments with Escherichia coli and demonstrate its application for estimating relative fitness. Our approach provides an effective way to predict growth and infer relative fitness in mixed cultures.

Inferring the distribution of fitness effects of spontaneous mutations in Chlamydomonas reinhardtii

Spontaneous mutations are the source of new genetic variation and are thus central to the evolutionary process. In molecular evolution and quantitative genetics, the nature of genetic variation depends critically on the distribution of effects of mutations on fitness and other quantitative traits. Spontaneous mutation accumulation (MA) experiments have been the principal approach for investigating the overall rate of occurrence and cumulative effect of mutations but have not allowed the phenotypic effects of individual mutations to be studied directly. Here, we crossed MA lines of the green alga Chlamydomonas reinhardtii with its unmutated ancestral strain to create haploid recombinant lines, each carrying an average of 50% of the accumulated mutations in a large number of combinations. With the aid of the genome sequences of the MA lines, we inferred the genotypes of the mutations, assayed their growth rate as a measure of fitness, and inferred the distribution of fitness effects (DFE) using a Bayesian mixture model. We infer that the DFE is highly leptokurtic (L-shaped). Of mutations with absolute fitness effects exceeding 1%, about one-sixth increase fitness in the laboratory environment. The inferred distribution of effects for deleterious mutations is consistent with a strong role for nearly neutral evolution. Specifically, such a distribution predicts that nucleotide variation and genetic variation for quantitative traits will be insensitive to change in the effective population size.

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.

Massively Parallel Assays and Quantitative Sequence-Function Relationships

Over the last decade, a rich variety of massively parallel assays have revolutionized our understanding of how biological sequences encode quantitative molecular phenotypes. These assays include deep mutational scanning, high-throughput SELEX, and massively parallel reporter assays. Here, we review these experimental methods and how the data they produce can be used to quantitatively model sequence-function relationships. In doing so, we touch on a diverse range of topics, including the identification of clinically relevant genomic variants, the modeling of transcription factor binding to DNA, the functional and evolutionary landscapes of proteins, and cis-regulatory mechanisms in both transcription and mRNA splicing. We further describe a unified conceptual framework and a core set of mathematical modeling strategies that studies in these diverse areas can make use of. Finally, we highlight key aspects of experimental design and mathematical modeling that are important for the results of such studies to be interpretable and reproducible.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

Slower environmental change hinders adaptation from standing genetic variation

Evolutionary responses to environmental change depend on the time available for adaptation before environmental degradation leads to extinction. Explicit tests of this relationship are limited to microbes where adaptation usually depends on the sequential fixation of de novo mutations, excluding standing variation for genotype-by-environment fitness interactions that should be key for most natural species. For natural species evolving from standing genetic variation, adaptation at slower rates of environmental change may be impeded since the best genotypes at the most extreme environments can be lost during evolution due to genetic drift or founder effects. To address this hypothesis, we perform experimental evolution with self-fertilizing populations of the nematode Caenorhabditis elegans and develop an inference model to describe natural selection on extant genotypes under environmental change. Under a sudden environmental change, we find that selection rapidly increases the frequency of genotypes with high fitness in the most extreme environment. In contrast, under a gradual environmental change selection first favors genotypes that are worse at the most extreme environment. We demonstrate with a second set of evolution experiments that, as a consequence of slower environmental change and thus longer periods to reach the most extreme environments, genetic drift and founder effects can lead to the loss of the most beneficial genotypes. We further find that maintenance of standing genetic variation can retard the fixation of the best genotypes in the most extreme environment because of interference between them. Taken together, these results show that slower environmental change can hamper adaptation from standing genetic variation and they support theoretical models indicating that standing variation for genotype-by-environment fitness interactions critically alters the pace and outcome of adaptation under environmental change.

Expansion of the metazoan virosphere: progress, pitfalls, and prospects

Metagenomic sequencing has led to a recent and rapid expansion of the animal virome. It has uncovered a multitude of new virus lineages from under-sampled host groups, including many that break up long branches in the virus tree, and many that display unexpected genome sizes and structures. Although there are challenges to inferring the existence of a virus from a `virus-like sequence', in the absence of an isolate the analysis of nucleic acid (including small RNAs) and sequence data can provide considerable confidence. As a consequence, this period of molecular natural history is helping to reshape our views of deep virus evolution.

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.

Fitness effects of altering gene expression noise in Saccharomyces cerevisiae

Gene expression noise is an evolvable property of biological systems that describes differences in expression among genetically identical cells in the same environment. Prior work has shown that expression noise is heritable and can be shaped by selection, but the impact of variation in expression noise on organismal fitness has proven difficult to measure. Here, we quantify the fitness effects of altering expression noise for the TDH3 gene in Saccharomyces cerevisiae. We show that increases in expression noise can be deleterious or beneficial depending on the difference between the average expression level of a genotype and the expression level maximizing fitness. We also show that a simple model relating single-cell expression levels to population growth produces patterns consistent with our empirical data. We use this model to explore a broad range of average expression levels and expression noise, providing additional insight into the fitness effects of variation in expression noise.

Experimental demonstration of the impact of hard and soft selection regimes on polymorphism maintenance in spatially heterogeneous environments

Predicting and managing contemporary adaption requires a proper understanding of the determinants of genetic variation. Spatial heterogeneity of the environment may stably maintain polymorphism when habitat contribution to the next generation can be considered independent of the degree of adaptation of local populations within habitats (i.e., under soft selection). In contrast, when habitats contribute proportionally to the mean fitness of the populations they host (hard selection), polymorphism is not expected to be maintained by selection. Although mathematically established decades ago, this prediction had never been demonstrated experimentally. Here, we provide an experimental test in which polymorphic populations of Escherichia coli growing in heterogeneous habitats were exposed to hard and soft selection regimes. As predicted by theory, polymorphism was preserved longer under soft selection. Complementary tests established that soft selection slowed fixation processes and could even protect polymorphism in the long term by providing a systematic advantage to rare genotypes.

Evolution of bacteria specialization along an antibiotic dose gradient

Antibiotic and pesticide resistance of pathogens are major and pressing worldwide issues. Resistance evolution is often considered in simplified ecological contexts: treated versus nontreated environments. In contrast, antibiotic usually present important dose gradients: from ecosystems to hospitals to polluted soils, in treated patients across tissues. However, we do not know whether adaptation to low or high doses involves different phenotypic traits, and whether these traits trade‐off with each other. In this study, we investigated the occurrence of such fitness trade‐offs along a dose gradient by evolving experimentally resistant lines of Escherichia coli at different antibiotic concentrations for ∼400 generations. Our results reveal fast evolution toward specialization following the first mutational step toward resistance, along with pervasive trade‐offs among different evolution doses. We found clear and regular fitness patterns of specialization, which converged rapidly from different initial starting points. These findings are consistent with a simple fitness peak shift model as described by the classical evolutionary ecology theory of adaptation across environmental gradients. We also found that the fitness costs of resistance tend to be compensated through time at low doses whereas they increase through time at higher doses. This cost evolution follows a linear trend with the log‐dose of antibiotic along the gradient. These results suggest a general explanation for the variability of the fitness costs of resistance and their evolution. Overall, these findings call for more realistic models of resistance management incorporating dose‐specialization.

Genetic Adaptation to Growth Under Laboratory Conditions in Escherichia coli and Salmonella enterica

Experimental evolution under controlled laboratory conditions is becoming increasingly important to address various evolutionary questions, including, for example, the dynamics and mechanisms of genetic adaptation to different growth and stress conditions. In such experiments, mutations typically appear that increase the fitness under the conditions tested (medium adaptation), but that are not necessarily of interest for the specific research question. Here, we have identified mutations that appeared during serial passage of E. coli and S. enterica in four different and commonly used laboratory media and measured the relative competitive fitness and maximum growth rate of 111 genetically re-constituted strains, carrying different single and multiple mutations. Little overlap was found between the mutations that were selected in the two species and the different media, implying that adaptation occurs via different genetic pathways. Furthermore, we show that commonly occurring adaptive mutations can generate undesired genetic variation in a population and reduce the accuracy of competition experiments. However, by introducing media adaptation mutations with large effects into the parental strain that was used for the evolution experiment, the variation (standard deviation) was decreased 10-fold, and it was possible to measure fitness differences between two competitors as small as |s| < 0.001.

The evolution of bacterial cell size: the internal diffusion-constraint hypothesis

Size is one of the most important biological traits influencing organismal ecology and evolution. However, we know little about the drivers of body size evolution in unicellulars. A long-term evolution experiment (Lenski's LTEE) in which Escherichia coli adapts to a simple glucose medium has shown that not only the growth rate and the fitness of the bacterium increase over time but also its cell size. This increase in size contradicts prominent 'external diffusion' theory (EDC) predicting that cell size should have evolved toward smaller cells. Among several scenarios, we propose and test an alternative 'internal diffusion-constraint' (IDC) hypothesis for cell size evolution. A change in cell volume affects metabolite concentrations in the cytoplasm. The IDC states that a higher metabolism can be achieved by a reduction in the molecular traffic time inside of the cell, by increasing its volume. To test this hypothesis, we studied a population from the LTEE. We show that bigger cells with greater growth and CO₂ production rates and lower mass-to-volume ratio were selected over time in the LTEE. These results are consistent with the IDC hypothesis. This novel hypothesis offers a promising approach for understanding the evolutionary constraints on cell size.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.

Fisher's geometrical model and the mutational patterns of antibiotic resistance across dose gradients

Fisher's geometrical model (FGM) has been widely used to depict the fitness effects of mutations. It is a general model with few underlying assumptions that gives a large and comprehensive view of adaptive processes. It is thus attractive in several situations, for example adaptation to antibiotics, but comes with limitations, so that more mechanistic approaches are often preferred to interpret experimental data. It might be possible however to extend FGM assumptions to better account for mutational data. This is theoretically challenging in the context of antibiotic resistance because resistance mutations are assumed to be rare. In this article, we show with Escherichia coli how the fitness effects of resistance mutations screened at different doses of nalidixic acid vary across a dose-gradient. We found experimental patterns qualitatively consistent with the basic FGM (rate of resistance across doses, gamma distributed costs) but also unexpected patterns such as a decreasing mean cost of resistance with increasing screen dose. We show how different extensions involving mutational modules and variations in trait covariance across environments, can be discriminated based on these data. Overall, simple extensions of the FGM accounted well for complex mutational effects of resistance mutations across antibiotic doses.

Role of epistasis on the fixation probability of a non-mutator in an adapted asexual population

The mutation rate of a well adapted population is prone to reduction so as to have a lower mutational load. We aim to understand the role of epistatic interactions between the fitness affecting mutations in this process. Using a multitype branching process, the fixation probability of a single non-mutator emerging in a large asexual mutator population is analytically calculated here. The mutator population undergoes deleterious mutations at constant, but at a much higher rate than that of the non-mutator. We find that antagonistic epistasis lowers the chances of mutation rate reduction, while synergistic epistasis enhances it. Below a critical value of epistasis, the fixation probability behaves non-monotonically with variation in the mutation rate of the background population. Moreover, the variation of this critical value of the epistasis parameter with the strength of the mutator is discussed in the appendix. For synergistic epistasis, when selection is varied, the fixation probability reduces overall, with damped oscillations.

Relating fitness to long-term environmental variations in natura

Quantifying links between ecological processes and adaptation dynamics in natura remains a crucial challenge. Many studies have documented the strength, form and direction of selection, and its variations in space and time, but only a few managed to link these variations to their proximal causes. This step is, however, crucial, if we are to understand how the variation in selective pressure affects adaptive allele dynamics in natural settings. We used data from a long-term survey (about 30 years) monitoring the adaptation to insecticides of Culex pipiens mosquitoes in Montpellier area (France), focusing on three resistance alleles of the Ester locus. We used a population genetics model taking temporal and spatial variations in selective pressure into account, to assess the quantitative relationships between variations in the proximal agent of selection (amounts of insecticide sprayed) and the fitness of resistance alleles. The response to variations in selective pressure was fast, and the alleles displayed different fitness-to-environment relationships: the analyses revealed that even slight changes in insecticide doses could induce changes in the strength and direction of selection, thus changing the fitness ranking of the adaptive alleles. They also revealed that selective pressures other than the insecticides used for mosquito control affected the resistance allele dynamics. These fitness-to-environment relationships, fast responses and continuous evolution limit our ability to predict the outcome of adaptive allele dynamics in a changing environment, but they clearly contribute to the maintenance of polymorphism in natural populations. Our study also emphasizes the necessity of long-term surveys in evolutionary ecology.

The Fitness Effects of Spontaneous Mutations Nearly Unseen by Selection in a Bacterium with Multiple Chromosomes

Mutation accumulation (MA) experiments employ the strategy of minimizing the population size of evolving lineages to greatly reduce effects of selection on newly arising mutations. Thus, most mutations fix within MA lines independently of their fitness effects. This approach, more recently combined with genome sequencing, has detailed the rates, spectra, and biases of different mutational processes. However, a quantitative understanding of the fitness effects of mutations virtually unseen by selection has remained an untapped opportunity. Here, we analyzed the fitness of 43 sequenced MA lines of the multi-chromosome bacterium Burkholderia cenocepacia that had each undergone 5554 generations of MA and accumulated an average of 6.73 spontaneous mutations. Most lineages exhibited either neutral or deleterious fitness in three different environments in comparison with their common ancestor. The only mutational class that was significantly overrepresented in lineages with reduced fitness was the loss of the plasmid, though nonsense mutations, missense mutations, and coding insertion-deletions were also overrepresented in MA lineages whose fitness had significantly declined. Although the overall distribution of fitness effects was similar between the three environments, the magnitude and even the sign of the fitness of a number of lineages changed with the environment, demonstrating that the fitness of some genotypes was environmentally dependent. These results present an unprecedented picture of the fitness effects of spontaneous mutations in a bacterium with multiple chromosomes and provide greater quantitative support for the theory that the vast majority of spontaneous mutations are neutral or deleterious.

The Effects of Population Size Histories on Estimates of Selection Coefficients from Time-Series Genetic Data

Many approaches have been developed for inferring selection coefficients from time series data while accounting for genetic drift. These approaches have been motivated by the intuition that properly accounting for the population size history can significantly improve estimates of selective strengths. However, the improvement in inference accuracy that can be attained by modeling drift has not been characterized. Here, by comparing maximum likelihood estimates of selection coefficients that account for the true population size history with estimates that ignore drift by assuming allele frequencies evolve deterministically in a population of infinite size, we address the following questions: how much can modeling the population size history improve estimates of selection coefficients? How much can mis-inferred population sizes hurt inferences of selection coefficients? We conduct our analysis under the discrete Wright–Fisher model by deriving the exact probability of an allele frequency trajectory in a population of time-varying size and we replicate our results under the diffusion model. For both models, we find that ignoring drift leads to estimates of selection coefficients that are nearly as accurate as estimates that account for the true population history, even when population sizes are small and drift is high. This result is of interest because inference methods that ignore drift are widely used in evolutionary studies and can be many orders of magnitude faster than methods that account for population sizes.

Environmental variation alters the fitness effects of rifampicin resistance mutations in Pseudomonas aeruginosa

The fitness effects of antibiotic resistance mutations in antibiotic-free conditions play a key role in determining the long-term maintenance of resistance. Although resistance is usually associated with a cost, the impact of environmental variation on the cost of resistance is poorly understood. Here, we test the impact of heterogeneity in temperature and resource availability on the fitness effects of antibiotic resistance using strains of the pathogenic bacterium Pseudomonas aeruginosa carrying clinically important rifampicin resistance mutations. Although the rank order of fitness was generally maintained across environments, fitness effects relative to the wild type differed significantly. Changes in temperature had a profound impact on the fitness effects of resistance, whereas changes in carbon substrate had only a weak impact. This suggests that environmental heterogeneity may influence whether the costs of resistance are likely to be ameliorated by second-site compensatory mutations or by reversion to wild-type rpoB. Our results highlight the need to consider environmental heterogeneity and genotype-by-environment interactions for fitness in models of resistance evolution.

Predation and selection for antibiotic resistance in natural environments

Genes encoding resistance to antibiotics appear, like the antibiotics themselves, to be ancient, originating long before the rise of the era of anthropogenic antibiotics. However, detailed understanding of the specific biological advantages of antibiotic resistance in natural environments is still lacking, thus limiting our efforts to prevent environmental influx of resistance genes. Here, we propose that antibiotic-resistant cells not only evade predation from antibiotic producers but also take advantage of nutrients released from cells that are killed by the antibiotic-producing bacteria. Thus, predation is potentially an important mechanism for driving antibiotic resistance during slow or stationary phase of growth when nutrients are deprived. This adds to explain the ancient nature and widespread occurrence of antibiotic resistance in natural environments unaffected by anthropogenic antibiotics. In particular, we suggest that nutrient-poor environments including indoor environments, for example, clean rooms and intensive care units may serve as a reservoir and source for antibiotic-producing as well as antibiotic-resistant bacteria.

The Impact of Selection, Gene Conversion, and Biased Sampling on the Assessment of Microbial Demography

Recent studies have linked demographic changes and epidemiological patterns in bacterial populations using coalescent-based approaches. We identified 26 studies using skyline plots and found that 21 inferred overall population expansion. This surprising result led us to analyze the impact of natural selection, recombination (gene conversion), and sampling biases on demographic inference using skyline plots and site frequency spectra (SFS). Forward simulations based on biologically relevant parameters from Escherichia coli populations showed that theoretical arguments on the detrimental impact of recombination and especially natural selection on the reconstructed genealogies cannot be ignored in practice. In fact, both processes systematically lead to spurious interpretations of population expansion in skyline plots (and in SFS for selection). Weak purifying selection, and especially positive selection, had important effects on skyline plots, showing patterns akin to those of population expansions. State-of-the-art techniques to remove recombination further amplified these biases. We simulated three common sampling biases in microbiological research: uniform, clustered, and mixed sampling. Alone, or together with recombination and selection, they further mislead demographic inferences producing almost any possible skyline shape or SFS. Interestingly, sampling sub-populations also affected skyline plots and SFS, because the coalescent rates of populations and their sub-populations had different distributions. This study suggests that extreme caution is needed to infer demographic changes solely based on reconstructed genealogies. We suggest that the development of novel sampling strategies and the joint analyzes of diverse population genetic methods are strictly necessary to estimate demographic changes in populations where selection, recombination, and biased sampling are present.

Diverse phenotypic and genetic responses to short-term selection in evolving Escherichia coli populations

Beneficial mutations fuel adaptation by altering phenotypes that enhance the fit of organisms to their environment. However, the phenotypic effects of mutations often depend on ecological context, making the distribution of effects across multiple environments essential to understanding the true nature of beneficial mutations. Studies that address both the genetic basis and ecological consequences of adaptive mutations remain rare. Here, we characterize the direct and pleiotropic fitness effects of a collection of 21 first-step beneficial mutants derived from naïve and adapted genotypes used in a long-term experimental evolution of Escherichia coli. Whole-genome sequencing was able to identify the majority of beneficial mutations. In contrast to previous studies, we find diverse fitness effects of mutations selected in a simple environment and few cases of genetic parallelism. The pleiotropic effects of these mutations were predominantly positive but some mutants were highly antagonistic in alternative environments. Further, the fitness effects of mutations derived from the adapted genotypes were dramatically reduced in nearly all environments. These findings suggest that many beneficial variants are accessible from a single point on the fitness landscape, and the fixation of alternative beneficial mutations may have dramatic consequences for niche breadth reduction via metabolic erosion.

Constraints on adaptation of Escherichia coli to mixed-resource environments increase over time

Can a population evolved in two resources reach the same fitness in both as specialist populations evolved in each of the individual resources? This question is central to theories of ecological specialization, the maintenance of genetic variation, and sympatric speciation, yet relatively few experiments have examined costs of generalism over long-term adaptation. We tested whether selection in environments containing two resources limits a population's ability to adapt to the individual resources by comparing the fitness of replicate Escherichia coli populations evolved for 6000 generations in the presence of glucose or lactose alone (specialists), or in varying presentations of glucose and lactose together (generalists). We found that all populations had significant fitness increases in both resources, though the magnitude and rate of these increases differed. For the first 4000 generations, most generalist populations increased in fitness as quickly in the individual resources as the corresponding specialist populations. From 5000 generations, however, a widespread cost of adaptation affected all generalists, indicating a growing constraint on their abilities to adapt to two resources simultaneously. Our results indicate that costs of generalism are prevalent, but may influence evolutionary trajectories only after a period of cost-free adaptation.

Malthusian Parameters as Estimators of the Fitness of Microbes: A Cautionary Tale about the Low Side of High Throughput

The maximum exponential growth rate, the Malthusian parameter (MP), is commonly used as a measure of fitness in experimental studies of adaptive evolution and of the effects of antibiotic resistance and other genes on the fitness of planktonic microbes. Thanks to automated, multi-well optical density plate readers and computers, with little hands-on effort investigators can readily obtain hundreds of estimates of MPs in less than a day. Here we compare estimates of the relative fitness of antibiotic susceptible and resistant strains of E. coli, Pseudomonas aeruginosa and Staphylococcus aureus based on MP data obtained with automated multi-well plate readers with the results from pairwise competition experiments. This leads us to question the reliability of estimates of MP obtained with these high throughput devices and the utility of these estimates of the maximum growth rates to detect fitness differences.

The spectrum of adaptive mutations in experimental evolution

A primary goal of recent work in experimental evolution is to probe the molecular basis of adaptation. This requires an understanding of the individual mutations in evolving populations: their identity, their physiological and fitness effects, and the interactions between them. The combination of high-throughput methods for laboratory evolution and next-generation sequencing methods now makes it possible to identify and quantify mutations in hundreds of replicate populations over thousands of generations, and to directly measure fitness effects and epistatic interactions. Many laboratories are now leveraging these tools to study the molecular basis of adaptation and the reproducibility of evolutionary outcomes across a variety of model systems. Genetic analyses on evolved populations are shedding light on the statistics of epistasis between evolved mutations. Here we review the current understanding of the spectrum of mutations observed across these systems, with a focus on epistatic interactions between beneficial mutations and constraints on evolutionary outcomes. We emphasize evolution in asexual microbes, where next generation sequencing methods have been widely applied.

Neutral null models for diversity in serial transfer evolution experiments

Evolution experiments with microorganisms coupled with genome-wide sequencing now allow for the systematic study of population genetic processes under a wide range of conditions. In learning about these processes in natural, sexual populations, neutral models that describe the behavior of diversity and divergence summaries have played a pivotal role. It is therefore natural to ask whether neutral models, suitably modified, could be useful in the context of evolution experiments. Here, we introduce coalescent models for polymorphism and divergence under the most common experimental evolution assay, a serial transfer experiment. This relatively simple setting allows us to address several issues that could affect diversity patterns in evolution experiments, whether selection is operating or not: the transient behavior of neutral polymorphism in an experiment beginning from a single clone, the effects of randomness in the timing of cell division and noisiness in population size in the dilution stage. In our analyses and discussion, we emphasize the implications for experiments aimed at measuring diversity patterns and making inferences about population genetic processes based on these measurements.

The dynamics of niche evolution upon abrupt environmental change

Abrupt environmental changes are of particular interest to understand how species can quickly evolve at the boundary of their current niche. In particular the "sliding niche" model, wherein a niche shifts globally toward the new condition, has been used in understanding and modeling this process. Here, we investigate the dynamics of relative fitness change in four evolutionary replicates of Escherichia coli populations exposed to an extreme pH shift. We analyzed these changes at generations 500, 1000, and 2000 to determine whether niche global deformations fully capture the temporal dynamics of niche evolution. Strikingly, this analysis reveals that fitness variations can indeed be attributed to simple and global deformation of an underlying simple niche template. Analysis from two experimental replicates displays a transient increase in niche width, consistent with recent theory considering plasticity evolution in the context of an abrupt environmental change. We term this scenario the "sidestep niche model."