Estimation of the rate and effect of new beneficial mutations in asexual populations

The rate and effect of available beneficial mutations are key parameters in determining how a population adapts to a new environment. However, these parameters are poorly known, in large part because of the difficulty of designing and interpreting experiments to examine the rare and intrinsically stochastic process of mutation occurrence. We present a new approach to estimate the rate and selective advantage of beneficial mutations that underlie the adaptation of asexual populations. We base our approach on the analysis of experiments that track the effect of newly arising beneficial mutations on the dynamics of a neutral marker in evolving bacterial populations and develop efficient estimators of mutation rate and selective advantage. Using extensive simulations, we evaluate the accuracy of our estimators and conclude that they are quite robust to the use of relatively low experimental replication. To validate the predictions of our model, we compare theoretical and experimentally determined estimates of the selective advantage of the first beneficial mutation to fix in a series of ten replicate populations. We find that our theoretical predictions are not significantly different from experimentally determined selection coefficients. Application of our method to suitably designed experiments will allow estimation of how population evolvability depends on demographic and initial fitness parameters.     

Modelling the fracture-healing process as a moving-interface problem using an interface-capturing approach

We present a novel numerical model of the fracture-healing process using interface-capturing techniques, a well-known approach from fields like fluid dynamics, to describe tissue growth. One advantage of this method is its direct connection to experimentally observable parameters, including tissue-growth velocities. In our model, osteogenesis, chondrogenesis and revascularisation are triggered by mechanical stimuli via mechano-transduction based on previously established hypothesis of Claes and Heigele. After experimentally verifying the convergence of the numerical method, we compare the predictions of our model with those of the already established Ulm bone-healing model, which serves as a benchmark, and corroborate our results with existing animal experiments. We demonstrate that the new model can predict the history of the interfragmentary movement and forecast a tissue evolution that appears similar to the experimental results. Furthermore, we compare the relative tissue concentration in the healing domain with outcomes of animal experiments. Finally, we discuss the possible application of the model to new fields, where numerical simulations could also prove beneficial.     

Eco‐evolutionary dynamics in a coevolving host–virus system

Eco‐evolutionary dynamics have been shown to be important for understanding population and community stability and their adaptive potential. However, coevolution in the framework of eco‐evolutionary theory has not been addressed directly. Combining experiments with an algal host and its viral parasite, and mathematical model analyses we show eco‐evolutionary dynamics in antagonistic coevolving populations. The interaction between antagonists initially resulted in arms race dynamics (ARD) with selective sweeps, causing oscillating host–virus population dynamics. However, ARD ended and populations stabilised after the evolution of a general resistant host, whereas a trade‐off between host resistance and growth then maintained host diversity over time (trade‐off driven dynamics). Most importantly, our study shows that the interaction between ecology and evolution had important consequences for the predictability of the mode and tempo of adaptive change and for the stability and adaptive potential of populations.     

Quasispecies Theory and the Behavior of RNA Viruses

A large number of medically important viruses, including HIV, hepatitis C virus, and influenza, have RNA genomes. These viruses replicate with extremely high mutation rates and exhibit significant genetic diversity. This diversity allows a viral population to rapidly adapt to dynamic environments and evolve resistance to vaccines and antiviral drugs. For the last 30 years, quasispecies theory has provided a population-based framework for understanding RNA viral evolution. A quasispecies is a cloud of diverse variants that are genetically linked through mutation, interact cooperatively on a functional level, and collectively contribute to the characteristics of the population. Many predictions of quasispecies theory run counter to traditional views of microbial behavior and evolution and have profound implications for our understanding of viral disease. Here, we discuss basic principles of quasispecies theory and describe its relevance for our understanding of viral fitness, virulence, and antiviral therapeutic strategy.     

Fitness recovery and compensatory evolution in natural mutant lines of C. elegans

Deleterious mutation accumulation plays a central role in evolutionary genetics, conservation biology, human health, and evolutionary medicine (e.g., methods of viral attenuation for live vaccines). It is therefore important to understand whether and how quickly populations with accumulated deleterious mutational loads can recover fitness through adaptive evolution. We used laboratory experimental evolution with four long-term mutation-accumulation (MA) lines of Caenorhabditis elegans nematodes to study the dynamics of such fitness evolution. We previously showed that when homozygous mutant populations are evolved in large population sizes, they can rapidly achieve wild-type fitness through the accumulation of new beneficial or compensatory epistatic mutations. Here, we expand this approach to demonstrate that when replicate lineages are initiated from the same mutant genotype, phenotypic evolution is only sometimes repeatable. MA genotypes that recovered ancestral fitness in the previous experiment did not always do so here. Further, the pattern of adaptive evolution in independently evolved replicates was contingent upon the MA genotype and varied among fitness-related traits. Our findings suggest that new beneficial mutations can drive rapid fitness evolution, but that the adaptive process is rendered somewhat unpredictable by its susceptibility to chance events and sensitivity to the evolutionary history of the starting population.     

A general theory for the evolutionary dynamics of virulence

Most theory on the evolution of virulence is based on a game-theoretic approach. One potential shortcoming of this approach is that it does not allow the prediction of the evolutionary dynamics of virulence. Such dynamics are of interest for several reasons: for experimental tests of theory, for the development of useful virulence management protocols, and for understanding virulence evolution in situations where the epidemiological dynamics never reach equilibrium and/or when evolutionary change occurs on a timescale comparable to that of the epidemiological dynamics. Here we present a general theory similar to that of quantitative genetics in evolutionary biology that allows for the easy construction of models that include both within-host mutation as well as superinfection and that is capable of predicting both the short- and long-term evolution of virulence. We illustrate the generality and intuitive appeal of the theory through a series of examples showing how it can lead to transparent interpretations of the selective forces governing virulence evolution. It also leads to novel predictions that are not possible using the game-theoretic approach. The general theory can be used to model the evolution of other pathogen traits as well.     

Evolution and population dynamics in stochastic environments

Inter-generational temporal variability of the environment is important in the evolution and adaptation of phenotypic traits. We discuss a population-dynamic approach which plays a central role in the analysis of evolutionary processes. The basic principle is that the phenotypes with the greatest long-term average growth rate will dominate the entire population. The calculation of longterm average growth rates for populations under temporal stochasticity can be highly cumbersome. However, for a discrete non-overlapping population, it is identical to the geometric mean of the growth rates (geometric mean fitness), which is usually different from the simple arithmetic mean of growth rates. Evolutionary outcomes based on geometric mean fitness are often very different from the predictions based on the usual arithmetic mean fitness. In this paper we illustrate the concept of geometric mean fitness in a few simple models. We discuss its implications for the adaptive evolution of phenotypes, e.g. foraging under predation risks and clutch size. Next, we present an application: the risk-spreading egg-laying behaviour of the cabbage white butterfly, and develop a two-patch population dynamic model to show how the optimal solution diverges from the ssual arithmetic mean approach. The dynamics of these stochastic models cannot be predicted from the dynamics of simple deterministic models. Thus the inclusion of stochastic factors in the analyses of populations is essential to the understanding of not only population dynamics, but also their evolutionary dynamics.     

Non-Selective Evolution of Growing Populations

Non-selective effects, like genetic drift, are an important factor in modern conceptions of evolution, and have been extensively studied for constant population sizes (Kimura, 1955; Otto and Whitlock, 1997). Here, we consider non-selective evolution in the case of growing populations that are of small size and have varying trait compositions (e.g. after a population bottleneck). We find that, in these conditions, populations never fixate to a trait, but tend to a random limit composition, and that the distribution of compositions “freezes” to a steady state. This final state is crucially influenced by the initial conditions. We obtain these findings from a combined theoretical and experimental approach, using multiple mixed subpopulations of two Pseudomonas putida strains in non-selective growth conditions (Matthijs et al, 2009) as model system. The experimental results for the population dynamics match the theoretical predictions based on the Pólya urn model (Eggenberger and Pólya, 1923) for all analyzed parameter regimes. In summary, we show that exponential growth stops genetic drift. This result contrasts with previous theoretical analyses of non-selective evolution (e.g. genetic drift), which investigated how traits spread and eventually take over populations (fixate) (Kimura, 1955; Otto and Whitlock, 1997). Moreover, our work highlights how deeply growth influences non-selective evolution, and how it plays a key role in maintaining genetic variability. Consequently, it is of particular importance in life-cycles models (Melbinger et al, 2010; Cremer et al, 2011; Cremer et al, 2012) of periodically shrinking and expanding populations.     

Evolution of RNA virus in spatially structured heterogeneous environments

A hallmark of the infectious cycle for many RNA viruses parasitizing multicellular hosts is the need to invade and successfully replicate in tissues that comprise a variety of cell types. Thus, multicellular hosts represent a heterogeneous environment to evolving viral populations. To understand viral adaptation to multicellular hosts, we took a double approach. First, we developed a mathematical model that served to make predictions concerning the dynamics of viral populations evolving in heterogeneous environments. Second, the predictions were tested by evolving vesicular stomatitis virus in vitro on a spatially structured environment formed by three different cell types. In the absence of gene flow, adaptation was tissue-specific, but fitness in all tissues decreased with migration rate. The performance in a given tissue was negatively correlated with its distance to the tissue hosting the population. This correlation decreased with migration rate.