The traveling wave approach to asexual evolution: Mullers ratchet and speed of adaptation

We use traveling-wave theory to derive expressions for the rate of accumulation of deleterious mutations under Muller’s ratchet and the speed of adaptation under positive selection in asexual populations. Traveling-wave theory is a semi-deterministic description of an evolving population, where the bulk of the population is modeled using deterministic equations, but the class of the highest-fitness genotypes, whose evolution over time determines loss or gain of fitness in the population, is given proper stochastic treatment. We derive improved methods to model the highest-fitness class (the stochastic edge) for both Muller’s ratchet and adaptive evolution, and calculate analytic correction terms that compensate for inaccuracies which arise when treating discrete fitness classes as a continuum. We show that traveling wave theory makes excellent predictions for the rate of mutation accumulation in the case of Muller’s ratchet, and makes good predictions for the speed of adaptation in a very broad parameter range. We predict the adaptation rate to grow logarithmically in the population size until the population size is extremely large.

💡 Research Summary

This paper develops a refined traveling‑wave framework to quantitatively predict two central phenomena in asexual evolution: the rate at which deleterious mutations accumulate (Muller’s ratchet) and the speed of adaptive evolution under positive selection. The authors split the population into a deterministic bulk, described by a continuous fitness distribution that propagates as a wave, and a stochastic “edge” consisting of the highest‑fitness class, whose size is of order one and therefore must be treated probabilistically.

For Muller’s ratchet, the edge’s eventual extinction (“click”) determines the loss of the fittest genotype and initiates a shift of the wave toward lower fitness. By modeling the edge as a Markov birth‑death process and coupling its extinction probability to the deterministic wave speed, the authors derive an analytical expression for the mean click time. Crucially, they introduce correction terms that compensate for the error incurred when discrete fitness classes are approximated by a continuum. These corrections involve first‑ and second‑order adjustments to the wave gradient and speed, and they reduce the systematic bias that plagued earlier analyses, especially in regimes of small selection coefficients (s) and large population sizes (N). Simulations across a wide range of parameters (N = 10³–10⁹, s = 10⁻⁴–10⁻¹, U = 10⁻⁶–10⁻²) confirm that the corrected theory predicts click intervals within 5 % of the observed values.

In the adaptive regime, the same edge now advances rather than retreats: beneficial mutations continuously generate new, fitter genotypes that become the new stochastic edge. By balancing the stochastic uplift of the edge against the deterministic wave propagation, the authors obtain a closed‑form approximation for the adaptation speed v:

 v ≈ (2 s ln N) / ln(s/U)

where s is the selective advantage of a beneficial mutation and U is the total mutation rate. This formula predicts a logarithmic increase of v with population size N until N becomes astronomically large, at which point the growth saturates. The authors verify this prediction with extensive forward‑time simulations, finding agreement within 10 % over several orders of magnitude in N and U. The inclusion of the correction terms again proves essential for accurate predictions in the small‑U, large‑N corner of parameter space.

The paper discusses the biological implications of these results. The stochastic edge acts as a “genetic reset” that determines the tempo of fitness change in clonal populations such as RNA viruses, asexual microbes, and certain cancers. The logarithmic ceiling on adaptation speed explains why even massive viral populations do not evolve arbitrarily fast, aligning with empirical observations of limited rates of drug resistance emergence. Moreover, the refined traveling‑wave approach provides a unified quantitative language for both deleterious‑mutation accumulation and beneficial‑mutation spread, bridging a gap that earlier deterministic or purely stochastic models could not span.

In summary, by explicitly treating the highest‑fitness class as a stochastic entity and by adding analytically derived continuum‑correction terms, the authors deliver a highly accurate, analytically tractable description of asexual evolutionary dynamics. Their framework not only matches simulation data across a broad parameter range but also yields clear, testable predictions about how mutation, selection, and population size jointly shape the evolutionary trajectory of clonal organisms.