Limits to adaptation in asexual populations

In asexual populations, the rate of adaptation is basically limited by the frequency and properties of spontaneous beneficial mutations. Hence, knowledge of these mutational properties and how they are affected by particular evolutionary conditions is a precondition for understanding the process of adaptation. Here, we address how the rate of adaptation of asexual populations is limited by its population size and mutation rate, as well as by two factors affecting the fraction of mutations that confer a benefit, i.e. the initial adaptedness of the population and the variability of the environment. These factors both influence which mutations are likely to occur, as well as the probability that they will ultimately contribute to adaptation. We attempt to separate the consequences of these basic population features in terms of their effect on the rate of adaptation by using results from evolution experiments with microorganisms.     

Synchronous Waves of Failed Soft Sweeps in the Laboratory: Remarkably Rampant Clonal Interference of Alleles at a Single Locus

It has increasingly been recognized that adapting populations of microbes contain not one, but many lineages continually arising and competing at once. This process, termed “clonal interference,” alters the rate and dynamics of adaptation and biases winning mutations toward those with the largest selective effect. Here we uncovered a dramatic example of clonal interference between multiple similar mutations occurring at the same locus within replicate populations of Methylobacterium extorquens AM1. Because these mutational events involved the transposition of an insertion sequence into a narrow window of a single gene, they were both readily detectable at low frequencies and could be distinguished due to differences in insertion sites. This allowed us to detect up to 17 beneficial alleles of this type coexisting in a single population. Despite conferring a large selective benefit, the majority of these alleles rose and then fell in frequency due to other lineages emerging that were more fit. By comparing allele-frequency dynamics to the trajectories of fitness gains by these populations, we estimated the fitness values of the genotypes that contained these mutations. Collectively across all populations, these alleles arose upon backgrounds with a wide range of fitness values. Within any single population, however, multiple alleles tended to rise and fall synchronously during a single wave of multiple genotypes with nearly identical fitness values. These results suggest that alleles of large benefit arose repeatedly in failed “soft sweeps” during narrow windows of adaptation due to the combined effects of epistasis and clonal interference.     

Fitness-dependent mutation rates in finite populations

Mutation rate may be condition dependent, whereby individuals in poor condition, perhaps from high mutation load, have higher mutation rates than individuals in good condition. Agrawal (J. Evol. Biol.15, 2002, 1004) explored the basic properties of fitness-dependent mutation rate (FDMR) in infinite populations and reported some heuristic results for finite populations. The key parameter governing how infinite populations evolve under FDMR is the curvature (k) of the relationship between fitness and mutation rate. We extend Agrawal's analysis to finite populations and consider dominance and epistasis. In finite populations, the probability of long-term existence depends on k. In sexual populations, positive curvature leads to low equilibrium mutation rate, whereas negative curvature results in high mutation rate. In asexual populations, negative curvature results in rapid extinction via 'mutational meltdown', whereas positive curvature sometimes allows persistence. We speculate that fitness-dependent mutation rate may provide the conditions for genetic architecture to diverge between sexual and asexual taxa.     

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.     

The effect of spatial structure on adaptation in Escherichia coli

Populations of organisms are generally organized in a given spatial structure. However, the vast majority of population genetic studies are based on populations in which every individual competes globally. Here we use experimental evolution in Escherichia coli to directly test a recently made prediction that spatial structure slows down adaptation and that this cost increases with the mutation rate. This was studied by comparing populations of different mutation rates adapting to a liquid (unstructured) medium with populations that evolved in a Petri dish on solid (structured) medium. We find that mutators adapt faster to both environments and that adaptation is slower if there is spatial structure. We observed no significant difference in the cost of structure between mutator and wild-type populations, which suggests that clonal interference is intense in both genetic backgrounds.     

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.     

Epistasis and frequency dependence influence the fitness of an adaptive mutation in a diversifying lineage

The opportunity for a mutation to invade a population can dramatically vary depending on the context in which this mutation occurs. Such context dependence is difficult to document as it requires the ability to measure how a mutation affects phenotypes and fitness and to manipulate the context in which the mutation occurs. We identified a mutation in a gene encoding a global regulator in one of two ecotypes that diverged from a common ancestor during 1200 generations of experimental evolution. We replaced the ancestral allele by the mutant allele, and vice versa, in several clones isolated during the time course of the evolution experiment, and compared the phenotype and fitness of clones isogenic except for the focal mutation. We show that the fitness and phenotype of the mutation are strongly affected by epistatic interactions between genes in the same genome, as well as by frequency dependent selection resulting from biotic interactions between individuals in the same population. We conclude that amongst the replicate population in which it spread, the mutation we identified is only adaptive when occurring in specific genomes and competing with specific individuals. This study thus demonstrates that the opportunity for an adaptive mutation to spread in an evolutionary lineage can only be understood in the light of its genomic and competitive environments.     

Upper-limit mutation rate estimation for a plant RNA virus

It is generally accepted that mutation rates of RNA viruses are inherently high due to the lack of proofreading mechanisms. However, direct estimates of mutation rate are surprisingly scarce, in particular for plant viruses. Here, based on the analysis of in vivo mutation frequencies in tobacco etch virus, we calculate an upper-bound mutation rate estimation of 3×10⁻⁵ per site and per round of replication; a value which turns out to be undistinguishable from the methodological error. Nonetheless, the value is barely on the lower side of the range accepted for RNA viruses, although in good agreement with the only direct estimate obtained for other plant viruses. These observations suggest that, perhaps, differences in the selective pressures operating during plant virus evolution may have driven their mutation rates towards values lower than those characteristic of other RNA viruses infecting bacteria or animals.     

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.     

Decay of unused characters by selection and drift

The reduction and loss of redundant phenotypic characters is a common feature of evolution. However, the mechanisms that drive deterioration of unused characters remain unclear. Here, we outline a simple framework where the relative importance of selective and neutral processes varies with environmental factors, because of variation in the fitness costs associated with unused traits. We tested our hypotheses using experimental evolution of the bacterium Pseudomonas fluorescens in spatially uniform environments. Results show that an unused character, swimming motility, decayed over evolutionary time and the rate of this decay varied among selection environments with different levels of resource availability. This is explained in the context of an environment-specific genetic correlation between motility and fitness, which is negative when resources are limited but neutral at higher resource levels. Thus, selection against an unused character was most effective in environments where the fitness cost was the greatest. This suggests that the same character can decay by different mechanisms depending upon environmental factors and supports previous evidence to show that resource availability can critically affect the outcomes of evolution.     

Statistical inference of sequence-dependent mutation rates

Several lines of research are now converging towards an integrated understanding of mutational mechanisms and their evolutionary implications. Experimentally, crystal structures reveal the effect of sequence context on polymerase fidelity; large-scale sequencing projects generate vast amounts of sequence polymorphism data; and locus-specific databases are being constructed. Computationally, software and analytical tools have been developed to analyze mutational data, to identify mutational hot spots, and to compare the signatures of mutagenic agents.     

Testicular melanization has evolved in birds with high mtDNA mutation rates

Melanin is mainly found in the integument of animals, but it also appears in several extracutaneous tissues. The presence of melanin in testes has been anecdotally reported in all vertebrate groups, but the causes and functions of this melanin remain unknown. Similar to other extracutaneous melanins, testicular melanin may protect male germ cells from oxidative stress. Given the high respiratory activity of spermatozoa, oxidative stress generated by mitochondrial dysfunction as a consequence of mtDNA mutations directly affects sperm viability. Thus, natural selection may favour testicular melanization in males of species with high historical mutation rates in the mitochondrial genome. Here, we tested this hypothesis using information on occurrence of testicular melanization and mutation accumulation as reflected by cytochrome b mtDNA base pair substitution rates in a large set of 134 species of birds, controlling for the confounding effects of body mass, reproductive activity and phylogeny. We found that testicular melanization has evolved in species with high rates of accumulated mitochondrial mutations and propose that this is an adaptive response related to the protective capacity of melanin against oxidative stress. In support of this hypothesis, testicular melanization was more frequently observed during the breeding season of birds (i.e. when spermatogenesis is likely to occur) than during reproductive inactivity. In contrast to other extracutaneous melanins whose abundance seems to reflect skin and coat colour, we did not find a correlation between the proportion of plumage coloured by melanins and occurrence of testicular melanization. Whereas future experimental studies should test these hypotheses, our study highlights for the first time that melanization patterns in animals may evolve as a response to historical mutation rates.     

Low abundance of Escherichia coli microsatellites is associated with an extremely low mutation rate

It is widely assumed that microsatellites are generated by replication slippage, a mutation process specific to repetitive DNA. Consistent with their high mutation rate, microsatellites are highly abundant in most eukaryotic genomes. In Escherichia coli, however, microsatellites are rare and short despite the fact that a high microsatellite mutation rate was described. We show that this high microsatellite instability depends on the presence of the F-plasmid. E. coli cells lacking the F-plasmid have extremely low microsatellite mutation rates. This result provides a possible explanation for the genome-wide low density of microsatellites in E. coli. Furthermore, we show that the F-plasmid induced microsatellite instability is independent of the mismatch repair pathway.     

The distribution of fitness effects of new beneficial mutations in Pseudomonas fluorescens

Theoretical studies of adaptation emphasize the importance of understanding the distribution of fitness effects (DFE) of new mutations. We report the isolation of 100 adaptive mutants—without the biasing influence of natural selection—from an ancestral genotype whose fitness in the niche occupied by the derived type is extremely low. The fitness of each derived genotype was determined relative to a single reference type and the fitness effects found to conform to a normal distribution. When fitness was measured in a different environment, the rank order changed, but not the shape of the distribution. We argue that, even with detailed knowledge of the genetic architecture underpinning the adaptive types (as is the case here), the DFEs remain unpredictable, and we discuss the possibility that general explanations for the shape of the DFE might not be possible in the absence of organism-specific biological details.     

Genic mutation rates in mammals: local similarity,chromosomal heterogeneity,and X-versus-autosome disparity

The reduction of mutation rates on the mammalian X chromosome relative to autosomes is most often explained in the literature as evidence of male-driven evolution. This hypothesis attributes lowered mutation rates on the X chromosome to the fact that this chromosome spends less time in the germline of males than in the germline of females. In contrast to this majority view, two articles argued that the patterns of mutation rates across chromosomes are inconsistent with male-driven evolution. One article reported a 40% reduction in synonymous substitution rates (Ks) for X-linked genes relative to autosomes in the mouse-rat lineage. The authors argued that this reduction is too dramatic to be explained by male-driven evolution and concluded that selection has systematically reduced mutation rate on the X chromosome to a level optimal for this male-hemizygous chromosome. More recently, a second article found that chromosomal mutation rates in both the human-mouse and mouse-rat lineages were so heterogeneous that the X chromosome was not an outlier. Here again, the authors argued that this is at odds with male-driven evolution and suggested that selection has modulated chromosomal mutation rates to locally optimal levels, thus extending the argument of the first mentioned article to include autosomes. Here, we reexamine these conclusions using mouse-rat and human-mouse coding-region data. We find a more modest reduction of Ks on the X chromosome, but our results contradict the finding that the X chromosome is not distinct from autosomes. Multiple statistical tests show that Ks rates on the X chromosome differ systematically from the autosomes in both lineages. We conclude that the moderate reduction of mutation rate on the X chromosome of both lineages is consistent with male-driven evolution; however, the large variance in mutation rates across chromosomes suggests that mutation rates are affected by additional factors besides male-driven evolution. Investigation of mutation rates by synteny reveals that synteny blocks, rather than entire chromosomes, might represent the unit of mutation rate variation.