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

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

Microsatellite evidence for high clonality and limited genetic diversity in Ziziphus celata (Rhamnaceae), an endangered, self-incompatible shrub endemic to the Lake Wales Ridge, Florida, USA

Genetic data are often crucial for designing management strategies for rare and endangered species. Ziziphus celata is an endangered sandhill shrub endemic to the Lake Wales Ridge of central Florida. This self-incompatible clonal species is known from only 14 wild populations, most of which are small (under 100 plants). Focusing on the five populations discovered in 2007, we evaluate the level of genetic diversity and identify clonal lineages within the wild populations of the species with a set of microsatellite loci. To account for somatic mutations and genotyping errors, we identified clonal lineages using a threshold cutoff for pair-wise genetic distances among samples. The microsatellites had up to 18 alleles/locus, and, consistent with outcrossing, samples were highly heterozygous (average population level H _o  = 0.69). Most populations of Z. celata consist of a single clone, and the most diverse population has only 10 clones. Overall Z. celata comprises 41 multi-locus genotypes, and 30 clonal lineages. With nearly 1,000 recorded plants (595 genotyped) and only 30 clonal lineages, Ziziphus celata is highly clonal: clonal richness, R = 0.049. The pair-wise distance method facilitates identification of clonal lineages, avoiding overestimation of clonal diversity. In most cases, the samples that grouped into a lineage were one to four plants differing from a surrounding genotype by a single microsatellite repeat insertion/deletion mutation, consistent with these having arisen via somatic mutations. Our data will enable managers to incorporate extant diversity from wild populations into ex situ collections. Additionally, our research demonstrates the utility of microsatellites for conservation of imperiled species, identifying genotypes of high priority for preservation.     

The Dynamics of Genetic Draft in Rapidly Adapting Populations

The accumulation of beneficial mutations on competing genetic backgrounds in rapidly adapting populations has a striking impact on evolutionary dynamics. This effect, known as clonal interference, causes erratic fluctuations in the frequencies of observed mutations, randomizes the fixation times of successful mutations, and leaves distinct signatures on patterns of genetic variation. Here, we show how this form of “genetic draft” affects the forward-time dynamics of site frequencies in rapidly adapting asexual populations. We calculate the probability that mutations at individual sites shift in frequency over a characteristic timescale, extending Gillespie’s original model of draft to the case where many strongly selected beneficial mutations segregate simultaneously. We then derive the sojourn time of mutant alleles, the expected fixation time of successful mutants, and the site frequency spectrum of beneficial and neutral mutations. Finally, we show how this form of draft affects inferences in the McDonald–Kreitman test and how it relates to recent observations that some aspects of genetic diversity are described by the Bolthausen–Sznitman coalescent in the limit of very rapid adaptation.     

EVOLUTIONARY REVERSALS OF ANTIBIOTIC RESISTANCE IN EXPERIMENTAL POPULATIONS OF PSEUDOMONAS AERUGINOSA

Antibiotic resistance mutations are accompanied by a fitness cost, and two mechanisms allow bacteria to adapt to this cost once antibiotic use is halted. First, it is possible for resistance to revert; second, it is possible for bacteria to adapt to the cost of resistance by compensatory mutations. Unfortunately, reversion to antibiotic sensitivity is rare, but the underlying factors that prevent reversion remain obscure. Here, we directly study the evolutionary dynamics of reversion by experimentally mimicking reversion mutations—sensitives—in populations of rifampicin‐resistant Pseudomonas aeruginosa. We show that, in our populations, most sensitives are lost due to genetic drift when they are rare. However, clonal interference from lineages carrying compensatory mutations causes a dramatic increase in the time to fixation of sensitives that escape genetic drift, and mutations surpassing the sensitives’ fitness are capable of driving transiently common sensitive lineages to extinction. Crucially, we show that the constraints on reversion arising from clonal interference are determined by the potential for compensatory adaptation of the resistant population. Although the cost of resistance provides the incentive for reversion, our study demonstrates that both the cost of resistance and the intrinsic evolvability of resistant populations interact to determine the rate and likelihood of reversion.     

Rate of Adaptation in Sexuals and Asexuals: A Solvable Model of the Fisher–Muller Effect

The adaptation of large asexual populations is hampered by the competition between independently arising beneficial mutations in different individuals, which is known as clonal interference. In classic work, Fisher and Muller proposed that recombination provides an evolutionary advantage in large populations by alleviating this competition. Based on recent progress in quantifying the speed of adaptation in asexual populations undergoing clonal interference, we present a detailed analysis of the Fisher–Muller mechanism for a model genome consisting of two loci with an infinite number of beneficial alleles each and multiplicative (nonepistatic) fitness effects. We solve the deterministic, infinite population dynamics exactly and show that, for a particular, natural mutation scheme, the speed of adaptation in sexuals is twice as large as in asexuals. This result is argued to hold for any nonzero value of the rate of recombination. Guided by the infinite population result and by previous work on asexual adaptation, we postulate an expression for the speed of adaptation in finite sexual populations that agrees with numerical simulations over a wide range of population sizes and recombination rates. The ratio of the sexual to asexual adaptation speed is a function of population size that increases in the clonal interference regime and approaches 2 for extremely large populations. The simulations also show that the imbalance between the numbers of accumulated mutations at the two loci is strongly suppressed even by a small amount of recombination. The generalization of the model to an arbitrary number L of loci is briefly discussed. If each offspring samples the alleles at each locus from the gene pool of the whole population rather than from two parents, the ratio of the sexual to asexual adaptation speed is approximately equal to L in large populations. A possible realization of this scenario is the reassortment of genetic material in RNA viruses with L genomic segments.     

Fitness evolution and the rise of mutator alleles in experimental Escherichia coli populations

We studied the evolution of high mutation rates and the evolution of fitness in three experimental populations of Escherichia coli adapting to a glucose-limited environment. We identified the mutations responsible for the high mutation rates and show that their rate of substitution in all three populations was too rapid to be accounted for simply by genetic drift. In two of the populations, large gains in fitness relative to the ancestor occurred as the mutator alleles rose to fixation, strongly supporting the conclusion that mutator alleles fixed by hitchhiking with beneficial mutations at other loci. In one population, no significant gain in fitness relative to the ancestor occurred in the population as a whole while the mutator allele rose to fixation, but a substantial and significant gain in fitness occurred in the mutator subpopulation as the mutator neared fixation. The spread of the mutator allele from rarity to fixation took >1000 generations in each population. We show that simultaneous adaptive gains in both the mutator and wild-type subpopulations (clonal interference) retarded the mutator fixation in at least one of the populations. We found little evidence that the evolution of high mutation rates accelerated adaptation in these populations.     

The First Steps of Adaptation of Escherichia coli to the Gut Are Dominated by Soft Sweeps

The accumulation of adaptive mutations is essential for survival in novel environments. However, in clonal populations with a high mutational supply, the power of natural selection is expected to be limited. This is due to clonal interference - the competition of clones carrying different beneficial mutations - which leads to the loss of many small effect mutations and fixation of large effect ones. If interference is abundant, then mechanisms for horizontal transfer of genes, which allow the immediate combination of beneficial alleles in a single background, are expected to evolve. However, the relevance of interference in natural complex environments, such as the gut, is poorly known. To address this issue, we have developed an experimental system which allows to uncover the nature of the adaptive process as Escherichia coli adapts to the mouse gut. This system shows the invasion of beneficial mutations in the bacterial populations and demonstrates the pervasiveness of clonal interference. The observed dynamics of change in frequency of beneficial mutations are consistent with soft sweeps, where different adaptive mutations with similar phenotypes, arise repeatedly on different haplotypes without reaching fixation. Despite the complexity of this ecosystem, the genetic basis of the adaptive mutations revealed a striking parallelism in independently evolving populations. This was mainly characterized by the insertion of transposable elements in both coding and regulatory regions of a few genes. Interestingly, in most populations we observed a complete phenotypic sweep without loss of genetic variation. The intense clonal interference during adaptation to the gut environment, here demonstrated, may be important for our understanding of the levels of strain diversity of E. coli inhabiting the human gut microbiota and of its recombination rate.     

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

Twelve replicate populations of Escherichia coli have been evolving in the laboratory for >25 years and 60,000 generations. We analyzed bacteria from whole-population samples frozen every 500 generations through 20,000 generations for one well-studied population, called Ara−1. By tracking 42 known mutations in these samples, we reconstructed the history of this population’s genotypic evolution over this period. The evolutionary dynamics of Ara−1 show strong evidence of selective sweeps as well as clonal interference between competing lineages bearing different beneficial mutations. In some cases, sets of several mutations approached fixation simultaneously, often conveying no information about their order of origination; we present several possible explanations for the existence of these mutational cohorts. Against a backdrop of rapid selective sweeps both earlier and later, two genetically diverged clades coexisted for >6000 generations before one went extinct. In that time, many additional mutations arose in the clade that eventually prevailed. We show that the clades evolved a frequency-dependent interaction, which prevented the immediate competitive exclusion of either clade, but which collapsed as beneficial mutations accumulated in the clade that prevailed. Clonal interference and frequency dependence can occur even in the simplest microbial populations. Furthermore, frequency dependence may generate dynamics that extend the period of coexistence that would otherwise be sustained by clonal interference alone.     

Clonal interference, multiple mutations and adaptation in large asexual populations

Two important problems affect the ability of asexual populations to accumulate beneficial mutations and hence to adapt. First, clonal interference causes some beneficial mutations to be outcompeted by more-fit mutations that occur in the same genetic background. Second, multiple mutations occur in some individuals, so even mutations of large effect can be outcompeted unless they occur in a good genetic background that contains other beneficial mutations. In this article, we use a Monte Carlo simulation to study how these two factors influence the adaptation of asexual populations. We find that the results depend qualitatively on the shape of the distribution of the fitness effects of possible beneficial mutations. When this distribution falls off slower than exponentially, clonal interference alone reasonably describes which mutations dominate the adaptation, although it gives a misleading picture of the evolutionary dynamics. When the distribution falls off faster than exponentially, an analysis based on multiple mutations is more appropriate. Using our simulations, we are able to explore the limits of validity of both of these approaches, and we explore the complex dynamics in the regimes where neither one is fully applicable.     

Genetic Diversity in the Interference Selection Limit

Pervasive natural selection can strongly influence observed patterns of genetic variation, but these effects remain poorly understood when multiple selected variants segregate in nearby regions of the genome. Classical population genetics fails to account for interference between linked mutations, which grows increasingly severe as the density of selected polymorphisms increases. Here, we describe a simple limit that emerges when interference is common, in which the fitness effects of individual mutations play a relatively minor role. Instead, similar to models of quantitative genetics, molecular evolution is determined by the variance in fitness within the population, defined over an effectively asexual segment of the genome (a “linkage block”). We exploit this insensitivity in a new “coarse-grained” coalescent framework, which approximates the effects of many weakly selected mutations with a smaller number of strongly selected mutations that create the same variance in fitness. This approximation generates accurate and efficient predictions for silent site variability when interference is common. However, these results suggest that there is reduced power to resolve individual selection pressures when interference is sufficiently widespread, since a broad range of parameters possess nearly identical patterns of silent site variability.     

Fitness effects of fixed beneficial mutations in microbial populations

Beneficial mutations are intuitively relevant to understanding adaptation, yet not all beneficial mutations are of consequence to the long-term evolutionary outcome of adaptation. Many beneficial mutations-mostly those of small effect-are lost due either to (1) genetic drift or to (2) competition among clones carrying different beneficial mutations, a phenomenon called the "Hill-Robertson effect" for sexual populations and "clonal interference" for asexual populations. Competition among clones becomes more prevalent with increasing genetic linkage and increasing population size, and it is thus generally characteristic of microbial populations. Together, these two phenomena suggest that only those beneficial mutations of large fitness effect should achieve fixation, despite the fact that most beneficial mutations produced are predicted to have very small fitness effects. Here, we confirm this prediction-both empirically and theoretically-by showing that fitness effects of fixed beneficial mutations follow a distribution whose mode is positive.     

The influence of deleterious mutations on adaptation in asexual populations

We study the dynamics of adaptation in asexual populations that undergo both beneficial and deleterious mutations. In particular, how the deleterious mutations affect the fixation of beneficial mutations was investigated. Using extensive Monte Carlo simulations, we find that in the "strong-selection weak mutation (SSWM)" regime or in the "clonal interference (CI)" regime, deleterious mutations rarely influence the distribution of "selection coefficients of the fixed mutations (SCFM)"; while in the "multiple mutations" regime, the accumulation of deleterious mutations would lead to a decrease in fitness significantly. We conclude that the effects of deleterious mutations on adaptation depend largely on the supply of beneficial mutations. And interestingly, the lowest adaptation rate occurs for a moderate value of selection coefficient of deleterious mutations.     

Adaptation in sexuals vs. asexuals: clonal interference and the Fisher-Muller model

Fisher and Muller's theory that recombination speeds adaptation by eliminating competition among beneficial mutations has proved a popular explanation for the advantage of sex. Recent theoretical studies have attempted to quantify the speed of adaptation under the Fisher-Muller model, partly in an attempt to understand the role of "clonal interference" in microbial experimental evolution. We reexamine adaptation in sexuals vs. asexuals, using a model of DNA sequence evolution. In this model, a modest number of sites can mutate to beneficial alleles and the fitness effects of these mutations are unequal. We study (1) transition probabilities to different beneficial mutations; (2) waiting times to the first and the last substitutions of beneficial mutations; and (3) trajectories of mean fitness through time. We find that some of these statistics are surprisingly similar between sexuals and asexuals. These results highlight the importance of the choice of substitution model in assessing the Fisher-Muller advantage of sex.     

The Properties of Adaptive Walks in Evolving Populations of Fungus

The rarity of beneficial mutations has frustrated efforts to develop a quantitative theory of adaptation. Recent models of adaptive walks, the sequential substitution of beneficial mutations by selection, make two compelling predictions: adaptive walks should be short, and fitness increases should become exponentially smaller as successive mutations fix. We estimated the number and fitness effects of beneficial mutations in each of 118 replicate lineages of Aspergillus nidulans evolving for approximately 800 generations at two population sizes using a novel maximum likelihood framework, the results of which were confirmed experimentally using sexual crosses. We find that adaptive walks do indeed tend to be short, and fitness increases become smaller as successive mutations fix. Moreover, we show that these patterns are associated with a decreasing supply of beneficial mutations as the population adapts. We also provide empirical distributions of fitness effects among mutations fixed at each step. Our results provide a first glimpse into the properties of multiple steps in an adaptive walk in asexual populations and lend empirical support to models of adaptation involving selection towards a single optimum phenotype. In practical terms, our results suggest that the bulk of adaptation is likely to be accomplished within the first few steps.     

Centennial clonal stability of asexual Daphnia in Greenland lakes despite climate variability

Climate and environmental condition drive biodiversity at many levels of biological organization, from populations to ecosystems. Combined with paleoecological reconstructions, palaeogenetic information on resident populations provides novel insights into evolutionary trajectories and genetic diversity driven by environmental variability. While temporal observations of changing genetic structure are often made of sexual populations, little is known about how environmental change affects the long‐term fate of asexual lineages. Here, we provide information on obligately asexual, triploid Daphnia populations from three Arctic lakes in West Greenland through the past 200–300 years to test the impact of environmental change on the temporal and spatial population genetic structure. The contrasting ecological state of the lakes, specifically regarding salinity and habitat structure may explain the observed lake‐specific clonal composition over time. Palaeolimnological reconstructions show considerable regional environmental fluctuations since 1,700 (the end of the Little Ice Age), but the population genetic structure in two lakes was almost unchanged with at most two clones per time period. Their local populations were strongly dominated by a single clone that has persisted for 250–300 years. We discuss possible explanations for the apparent population genetic stability: (a) persistent clones are general‐purpose genotypes that thrive under broad environmental conditions, (b) clonal lineages evolved subtle genotypic differences unresolved by microsatellite markers, or (c) epigenetic modifications allow for clonal adaptation to changing environmental conditions. Our results motivate research into the mechanisms of adaptation in these populations, as well as their evolutionary fate in the light of accelerating climate change in the polar regions.     

Diminishing returns of population size in the rate of RNA virus adaptation

Whenever an asexual viral population evolves by adapting to new environmental conditions, beneficial mutations, the ultimate cause of adaptation, are randomly produced and then fixed in the population. The larger the population size and the higher the mutation rate, the more beneficial mutations can be produced per unit time. With the usually high mutation rate of RNA viruses and in a large enough population, several beneficial mutations could arise at the same time but in different genetic backgrounds, and if the virus is asexual, they will never be brought together through recombination. Thus, the best of these genotypes must outcompete each other on their way to fixation. This competition among beneficial mutations has the effect of slowing the overall rate of adaptation. This phenomenon is known as clonal interference. Clonal interference predicts a speed limit for adaptation as the population size increases. In the present report, by varying the size of evolving vesicular stomatitis virus populations, we found evidence clearly demonstrating this speed limit and thus indicating that clonal interference might be an important factor modulating the rate of adaptation to an in vitro cell system. Several evolutionary and epidemiological implications of the clonal interference model applied to RNA viruses are discussed.     

Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria.

Selection of mutator alleles, increasing the mutation rate up to 10, 000-fold, has been observed during in vitro experimental evolution. This spread is ascribed to the hitchhiking of mutator alleles with favorable mutations, as demonstrated by a theoretical model using selective parameters corresponding to such experiments. Observations of unexpectedly high frequencies of mutators in natural isolates suggest that the same phenomenon could occur in the wild. But it remains questionable whether realistic in natura parameter values could also result in selection of mutators. In particular, the main parameters of adaptation, the size of the adapting population and the height and steepness of the adaptive peak characterizing adaptation, are very variable in nature. By simulation approach, we studied the effect of these parameters on the selection of mutators in asexual populations, assuming additive fitness. We show that the larger the population size, the more likely the fixation of mutator alleles. At a large population size, at least four adaptive mutations are needed for mutator fixation; moreover, under stronger selection stronger mutators are selected. We propose a model based on multiple mutations to illustrate how second-order selection can optimize population fitness when few favorable mutations are required for adaptation.     

Long-term effects of inducible mutagenic DNA repair on relative fitness and phenotypic diversification in Pseudomonas cichorii 302959

Mutagenic DNA repair (MDR) employs low-fidelity DNA polymerases capable of replicating past DNA lesions resulting from exposure to high-energy ultraviolet radiation (UVR). MDR confers UVR tolerance and activation initiates a transient mutator phenotype that may provide opportunities for adaptation. To investigate the potential role of MDR in adaptation, we have propagated parallel lineages of the highly mutable epiphytic plant pathogen Pseudomonas cichorii 302959 with daily UVR activation (UVR lineages) for ~500 generations. Here we examine those lineages through the measurement of relative fitness and observation of distinct colony morphotypes that emerged. Isolates and population samples from UVR lineages displayed gains in fitness relative to the ancestor despite increased rates of inducible mutation to rifampicin resistance. Regular activation of MDR resulted in the maintenance of genetic diversity within UVR lineages, including the reproducible diversification and coexistence of “round” and “fuzzy” colony morphotypes. These results suggest that inducible mutability may present a reasonable strategy for adaptive evolution in stressful environments by contributing to gains in relative fitness and diversification.     

Adaptation to High Ethanol Reveals Complex Evolutionary Pathways

Tolerance to high levels of ethanol is an ecologically and industrially relevant phenotype of microbes, but the molecular mechanisms underlying this complex trait remain largely unknown. Here, we use long-term experimental evolution of isogenic yeast populations of different initial ploidy to study adaptation to increasing levels of ethanol. Whole-genome sequencing of more than 30 evolved populations and over 100 adapted clones isolated throughout this two-year evolution experiment revealed how a complex interplay of de novo single nucleotide mutations, copy number variation, ploidy changes, mutator phenotypes, and clonal interference led to a significant increase in ethanol tolerance. Although the specific mutations differ between different evolved lineages, application of a novel computational pipeline, PheNetic, revealed that many mutations target functional modules involved in stress response, cell cycle regulation, DNA repair and respiration. Measuring the fitness effects of selected mutations introduced in non-evolved ethanol-sensitive cells revealed several adaptive mutations that had previously not been implicated in ethanol tolerance, including mutations in PRT1, VPS70 and MEX67. Interestingly, variation in VPS70 was recently identified as a QTL for ethanol tolerance in an industrial bio-ethanol strain. Taken together, our results show how, in contrast to adaptation to some other stresses, adaptation to a continuous complex and severe stress involves interplay of different evolutionary mechanisms. In addition, our study reveals functional modules involved in ethanol resistance and identifies several mutations that could help to improve the ethanol tolerance of industrial yeasts.     

Bacterial recombination promotes the evolution of multi-drug-resistance in functionally diverse populations

Bacterial recombination is believed to be a major factor explaining the prevalence of multi-drug-resistance (MDR) among pathogenic bacteria. Despite extensive evidence for exchange of resistance genes from retrospective sequence analyses, experimental evidence for the evolutionary benefits of bacterial recombination is scarce. We compared the evolution of MDR between populations of Acinetobacter baylyi in which we manipulated both the recombination rate and the initial diversity of strains with resistance to single drugs. In populations lacking recombination, the initial presence of multiple strains resistant to different antibiotics inhibits the evolution of MDR. However, in populations with recombination, the inhibitory effect of standing diversity is alleviated and MDR evolves rapidly. Moreover, only the presence of DNA harbouring resistance genes promotes the evolution of resistance, ruling out other proposed benefits for recombination. Together, these results provide direct evidence for the fitness benefits of bacterial recombination and show that this occurs by mitigation of functional interference between genotypes resistant to single antibiotics. Although analogous to previously described mechanisms of clonal interference among alternative beneficial mutations, our results actually highlight a different mechanism by which interactions among co-occurring strains determine the benefits of recombination for bacterial evolution.