Adaptation dynamics of the quasispecies model

We study the adaptation dynamics of an initially maladapted population evolving via the elementary processes of mutation and selection. The evolution occurs on rugged fitness landscapes which are defined on the multi-dimensional genotypic space and have many local peaks separated by low fitness valleys. We mainly focus on the Eigen’s model that describes the deterministic dynamics of an infinite number of self-replicating molecules. In the stationary state, for small mutation rates such a population forms a {\it quasispecies} which consists of the fittest genotype and its closely related mutants. The quasispecies dynamics on rugged fitness landscape follow a punctuated (or step-like) pattern in which a population jumps from a low fitness peak to a higher one, stays there for a considerable time before shifting the peak again and eventually reaches the global maximum of the fitness landscape. We calculate exactly several properties of this dynamical process within a simplified version of the quasispecies model.

💡 Research Summary

The paper investigates how an initially maladapted population evolves on rugged fitness landscapes when only mutation and selection act, using the deterministic Eigen quasispecies model for an infinite number of self‑replicating molecules. The authors first formulate the classic Eigen equation, where the concentration of genotype i evolves as

 dx_i/dt = Σ_j Q_{ij} f(j) x_j – Φ(t) x_i,

with Q_{ij} the mutation matrix (dependent on Hamming distance and per‑site mutation rate μ) and Φ(t) the mean replication rate. In the stationary regime and for sufficiently small μ, the population collapses into a “quasispecies”: a cloud of sequences centered on the fittest genotype and its close mutants.

When the fitness landscape contains many local peaks separated by low‑fitness valleys, the dynamics become punctuated. Starting from a low‑fitness peak, the population remains trapped for a long “search phase” while random mutations explore the surrounding genotype space. Once a mutant reaches a higher peak, a rapid “transition phase” follows: the mutant’s higher replication rate causes it to sweep the population, producing a step‑like jump in average fitness.

To quantify this process, the authors diagonalize the mutation‑selection operator, obtaining eigenvalues λ_k and eigenvectors v_k. The dominant eigenvalue λ_1 equals the mean fitness Φ, while the gap λ_1 – λ_2 controls the speed of peak shifts. They derive closed‑form expressions for the probability of moving from peak i to peak j after time t, showing that the transition rate scales roughly as

 r_{i→j} ≈ A exp