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.     

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.     

The speed of evolution and maintenance of variation in asexual populations

BACKGROUND: The rate at which beneficial mutations accumulate determines how fast asexual populations evolve, but this is only partially understood. Some recent clonal-interference models suggest that evolution in large asexual populations is limited because smaller beneficial mutations are outcompeted by larger beneficial mutations that occur in different lineages within the same population. This analysis assumes that the important mutations fix one at a time; it ignores multiple beneficial mutations that occur in the lineage of an earlier beneficial mutation, before the first mutation in the series can fix. We focus on the effects of such multiple mutations. RESULTS: Our analysis predicts that the variation in fitness maintained by a continuously evolving population increases as the logarithm of the population size and logarithm of the mutation rate and thus yields a similar logarithmic increase in the speed of evolution. To test these predictions, we evolved asexual budding yeast in glucose-limited media at a range of population sizes and mutation rates. CONCLUSIONS: We find that their evolution is dominated by the accumulation of multiple mutations of moderate effect. Our results agree with our theoretical predictions and are inconsistent with the one-by-one fixation of mutants assumed by recent clonal-interference analysis.     

The degeneration of asexual haploid populations and the speed of Muller's ratchet

The accumulation of deleterious mutations due to the process known as Muller's ratchet can lead to the degeneration of nonrecombining populations. We present an analytical approximation for the rate at which this process is expected to occur in a haploid population. The approximation is based on a diffusion equation and is valid when N exp(-u/s) > 1, where N is the population size, u is the rate at which deleterious mutations occur, and s is the effect of each mutation on fitness. Simulation results are presented to show that the approximation estimates the rate of the process better than previous approximations for values of mutation rates and selection coefficients that are compatible with the biological data. Under certain conditions, the ratchet can turn at a biologically significant rate when the deterministic equilibrium number of individuals free of mutations is substantially >100. The relevance of this process for the degeneration of Y or neo-Y chromosomes is discussed.     

The fate of transposable elements in asexual populations

Sexual reproduction and recombination are important for maintaining a stable copy number of transposable elements (TEs). In sexual populations, elements can be contained by purifying selection against host carriers with higher element copy numbers; however, in the absence of sex and recombination, asexual populations could be driven to extinction by an unchecked proliferation of TEs. Here we provide a theoretical framework for analyzing TE dynamics under asexual reproduction. Analytic results show that, in an infinite asexual population, an equilibrium in copy number is achieved if no element excision is possible, but that all TEs are eliminated if there is some excision. In a finite population, computer simulations demonstrate that small populations are driven to extinction by a Muller's ratchet-like process of element accumulation, but that large populations can be cured of vertically transmitted TEs, even with excision rates well below transposition rates. These results may have important consequences for newly arisen asexual lineages and may account for the lack of deleterious retrotransposons in the putatively ancient asexual bdelloid rotifers.     

The balance between mutators and nonmutators in asexual populations

Mutator alleles, which elevate an individual's mutation rate from 10 to 10,000-fold, have been found at high frequencies in many natural and experimental populations. Mutators are continually produced from nonmutators, often due to mutations in mismatch-repair genes. These mutators gradually accumulate deleterious mutations, limiting their spread. However, they can occasionally hitchhike to high frequencies with beneficial mutations. We study the interplay between these effects. We first analyze the dynamics of the balance between the production of mutator alleles and their elimination due to deleterious mutations. We find that when deleterious mutation rates are high in mutators, there will often be many "young," recently produced mutators in the population, and the fact that deleterious mutations only gradually eliminate individuals from a population is important. We then consider how this mutator-nonmutator balance can be disrupted by beneficial mutations and analyze the circumstances in which fixation of mutator alleles is likely. We find that dynamics is crucial: even in situations where selection on average acts against mutators, so they cannot stably invade, the mutators can still occasionally generate beneficial mutations and hence be important to the evolution of the population.     

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.     

Mutability as an altruistic trait in finite asexual populations

Mutation rate (MR) is a crucial determinant of the evolutionary process. Optimal MR may enable efficient evolutionary searching and therefore increase the fitness of the population over time. Nevertheless, individuals may favor MRs that are far from being optimal for the whole population. Instead, each individual may tend to mutate at rates that selfishly increase its own relative fitness. We show that in some cases, undergoing a mutation is altruistic, i.e., it increases the expected fitness of the population, but decreases the expected fitness of the mutated individual itself. In this case, if the population is uniform (completely mixed, undivided), immutability is evolutionary stable and is probably selected for. However, our examination of a segregated population, which is divided into several groups (or patches), shows that the optimal, altruistic MR may out-compete the selfish MR if the coupling between the groups is neither too strong nor too weak. This demonstrates that the population structure is crucial for the succession of the evolutionary process itself. For example, in a uniform population, the evolutionary process may be stopped before the highest fitness is reached, as demonstrated in a one-pick fitness landscape. In addition, we show that the dichotomy between evolutionary stable and optimal MRs can be seen as a special case of a more general phenomenon in which optimal behaviors may be destabilized in finite populations, since optimal sub-populations may become extinct before the benefit of their behavior is expressed.     

Local adaptation in rotifer populations

The adaptation of organisms to their environment has been a subject of study for a long time. One method to study adaptations in populations involves comparing contemporary populations of the same species under different selective regimes, in what is known as a ??local adaptation?? study. A previous study of the cyclically parthenogenetic rotifer Brachionus plicatilis found high heritabilities for some life-history traits. Some of these life-history traits significantly differed among six populations from Eastern Spain and data suggested some traits to have higher evolutionary rates than neutral genetic markers. Here, by studying the same B. plicatilis populations, we examine the variation and possible local adaptation of their main life-history traits, closely related to fitness, in relation to habitat salinity and temperature. These environmental factors have been shown to play a key role in the ecological differentiation among co-generic species of B. plicatilis. The results obtained in this study show that: (1) the seasonality of rotifer populations from Eastern Spain has profoundly influenced sexual reproduction strategies; (2) salinity is probably a key factor in the ecological specialization of some populations; and (3) rotifer populations harbour high variability in their fitness components.     

Genetic variation and the fate of beneficial mutations in asexual populations

The fate of a newly arising beneficial mutation depends on many factors, such as the population size and the availability and fitness effects of other mutations that accumulate in the population. It has proved difficult to understand how these factors influence the trajectories of particular mutations, since experiments have primarily focused on characterizing successful clones emerging from a small number of evolving populations. Here, we present the results of a massively parallel experiment designed to measure the full spectrum of possible fates of new beneficial mutations in hundreds of experimental yeast populations, whether these mutations are ultimately successful or not. Using strains in which a particular class of beneficial mutation is detectable by fluorescence, we followed the trajectories of these beneficial mutations across 592 independent populations for 1000 generations. We find that the fitness advantage provided by individual mutations plays a surprisingly small role. Rather, underlying "background" genetic variation is quickly generated in our initially clonal populations and plays a crucial role in determining the fate of each individual beneficial mutation in the evolving population.     

Parameters adaptation in the populations models

Ecology-evolutionary models of low dimensions were developed on the basis of competitive selection criteria. Dynamics of variables (number of individuals) and the search of evolutionary-stable values of parameters (biological characterictics of populations) were monitored in the suggested models. If the environmental temperature is changing periodically, the average (a) and width (d) of temperature tolerance range appears to be the important parameters. By model experiments it was established that stable values of temperature (a), favorable for development of highly specialized algae (d is low) were close to minimum and maximum of temperature curve. And for the low specialized algae (d is high) this values were close to the average temperature of environment. In a similar manner, a set of evolutionally stable parameters (a, d) was established for either of the two interacted populations (competitors and "predator-prey"). The hypotheses concerning it's geometric structure and the process of coevolution is formulated.     

Intriguing asexual life in marginal populations of the brown seaweed Fucus vesiculosus

Reproduction of attached large brown algae is known to occur only by sexual zygotes. Using microsatellites we show evolution of asexual reproduction in the bladder wrack promoting population persistence in the brackish water Baltic Sea (< 6 psu). Here a dwarf morph of Fucus vesiculosus is dominated by a single clone but clonal reproduction is also present in the common form of the species. We describe a possible mechanism for vegetative reproduction of attached algae, and conclude that clonality plays an important role in persistence and dispersal of these marginal populations, in which sexual reproduction is impaired by low salinity.     

Loss of least-loaded class in asexual populations due to drift and epistasis

We consider the dynamics of a nonrecombining haploid population of finite size that accumulates deleterious mutations irreversibly. This ratchet-like process occurs at a finite speed in the absence of epistasis, but it has been suggested that synergistic epistasis can halt the ratchet. Using a diffusion theory, we find explicit analytical expressions for the typical time between successive clicks of the ratchet for both nonepistatic and epistatic fitness functions. Our calculations show that the interclick time is of a scaling form that in the absence of epistasis gives a speed that is determined by size of the least-loaded class and the selection coefficient. With synergistic interactions, the ratchet speed is found to approach zero rapidly for arbitrary epistasis. Our analytical results are in good agreement with the numerical simulations.