Phase diagram for the Eigen quasispecies theory with a truncated fitness landscape

Virus Evolution on Fitness Landscapes

The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.

Phase transitions in evolutionary dynamics

Sharp changes in state, such as transitions from survival to extinction, are hallmarks of evolutionary dynamics in biological systems. These transitions can be explored using the techniques of statistical physics and the physics of nonlinear and complex systems. For example, a survival-to-extinction transition can be characterized as a non-equilibrium phase transition to an absorbing state. Here, we review the literature on phase transitions in evolutionary dynamics. We discuss directed percolation transitions in cellular automata and evolutionary models, and models that diverge from the directed percolation universality class. We explore in detail an example of an absorbing phase transition in an agent-based model of evolutionary dynamics, including previously unpublished data demonstrating similarity to, but also divergence from, directed percolation, as well as evidence for phase transition behavior at multiple levels of the model system's evolutionary structure. We discuss phase transition models of the error catastrophe in RNA virus dynamics and phase transition models for transition from chemistry to biochemistry, i.e., the origin of life. We conclude with a review of phase transition dynamics in models of natural selection, discuss the possible role of phase transitions in unraveling fundamental unresolved questions regarding multilevel selection and the major evolutionary transitions, and assess the future outlook for phase transitions in the investigation of evolutionary dynamics.

Distinct critical behaviors from the same state in quantum spin and population dynamics perspectives

There is a deep connection between the ground states of transverse-field spin systems and the late-time distributions of evolving viral populations-within simple models, both are obtained from the principal eigenvector of the same matrix. However, that vector is the wave-function amplitude in the quantum spin model, whereas it is the probability itself in the population model. We show that this seemingly minor difference has significant consequences: Phase transitions that are discontinuous in the spin system become continuous when viewed through the population perspective, and transitions that are continuous become governed by new critical exponents. We introduce a more general class of models that encompasses both cases and that can be solved exactly in a mean-field limit. Numerical results are also presented for a number of one-dimensional chains with power-law interactions. We see that well-worn spin models of quantum statistical mechanics can contain unexpected new physics and insights when treated as population-dynamical models and beyond, motivating further studies.

Mathematical Models of Quasi-Species Theory and Exact Results for the Dynamics

We formulate the Crow-Kimura, discrete-time Eigen model, and continuous-time Eigen model. These models are interrelated and we established an exact mapping between them. We consider the evolutionary dynamics for the single-peak fitness and symmetric smooth fitness. We applied the quantum mechanical methods to find the exact dynamics of the evolution model with a single-peak fitness. For the smooth symmetric fitness landscape, we map exactly the evolution equations into Hamilton-Jacobi equation (HJE). We apply the method to the Crow-Kimura (parallel) and Eigen models. We get simple formulas to calculate the dynamics of the maximum of distribution and the variance. We review the existing mathematical tools of quasi-species theory.

Phylogeny and evolution of RNA structure

Darwin's conviction that all living beings on Earth are related and the graph of relatedness is tree-shaped has been essentially confirmed by phylogenetic reconstruction first from morphology and later from data obtained by molecular sequencing. Limitations of the phylogenetic tree concept were recognized as more and more sequence information became available. The other path-breaking idea of Darwin, natural selection of fitter variants in populations, is cast into simple mathematical form and extended to mutation-selection dynamics. In this form the theory is directly applicable to RNA evolution in vitro and to virus evolution. Phylogeny and population dynamics of RNA provide complementary insights into evolution and the interplay between the two concepts will be pursued throughout this chapter. The two strategies for understanding evolution are ultimately related through the central paradigm of structural biology: sequence ⇒ structure ⇒ function. We elaborate on the state of the art in modeling both phylogeny and evolution of RNA driven by reproduction and mutation. Thereby the focus will be laid on models for phylogenetic sequence evolution as well as evolution and design of RNA structures with selected examples and notes on simulation methods. In the perspectives an attempt is made to combine molecular structure, population dynamics, and phylogeny in modeling evolution.

Stochastic simulations suggest that HIV-1 survives close to its error threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, , however, is not known. Application of the quasispecies theory to determine poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and . We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated to be substitutions/site/replication, ∼2–6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Currently available antiretroviral drugs curtail HIV infection but fail to eradicate the virus. A strategy of intervention radically different from that employed by current drugs has been proposed by the molecular quasispecies theory. The theory predicts that increasing the viral mutation rate beyond a critical value, called the error threshold, would cause a severe loss of genetic information, potentially leading to viral clearance. Several chemical mutagens are now being developed that can increase the mutation rate of HIV-1. Their success depends on reliable estimates of the error threshold of HIV-1, which are currently lacking. The quasispecies theory cannot be applied directly to HIV-1: the theory considers an infinitely large population of asexually reproducing haploid individuals, whereas HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We employed detailed stochastic simulations that overcome the limitations of the quasispecies theory and accurately mimic HIV-1 evolution in vivo. With these simulations, we estimated the error threshold of HIV-1 to be ∼2–6-fold higher than its natural mutation rate, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs.

Biological evolution in a multidimensional fitness landscape

We considered a multiblock molecular model of biological evolution, in which fitness is a function of the mean types of alleles located at different parts (blocks) of the genome. We formulated an infinite population model with selection and mutation, and calculated the mean fitness. For the case of recombination, we formulated a model with a multidimensional fitness landscape (the dimension of the space is equal to the number of blocks) and derived a theorem about the dynamics of initially narrow distribution. We also considered the case of lethal mutations. We also formulated the finite population version of the model in the case of lethal mutations. Our models, derived for the virus evolution, are interesting also for the statistical mechanics and the Hamilton-Jacobi equation as well.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Quasispecies as a matter of fact: viruses and beyond

We review the origins of the quasispecies concept and its relevance for RNA virus evolution, viral pathogenesis and antiviral treatment strategies. We emphasize a critical point of quasispecies that refers to genome collectivities as the unit of selection, and establish parallels between RNA viruses and some cellular systems such as bacteria and tumor cells. We refer also to tantalizing new observations that suggest quasispecies behavior in prions, perhaps as a result of the same quantum-mechanical indeterminations that underlie protein conformation and error-prone replication in genetic systems. If substantiated, these observations with prions could lead to new research on the structure-function relationship of non-nucleic acid biological molecules.

Dynamics of alternative modes of RNA replication for positive-sense RNA viruses

We propose and study nonlinear mathematical models describing the intracellular time dynamics of viral RNA accumulation for positive-sense single-stranded RNA viruses. Our models consider different replication modes ranging between two extremes represented by the geometric replication (GR) and the linear stamping machine replication (SMR). We first analyse a model that quantitatively reproduced experimental data for the accumulation dynamics of both polarities of turnip mosaic potyvirus RNAs. We identify a non-degenerate transcritical bifurcation governing the extinction of both strands depending on three key parameters: the mode of replication (α), the replication rate (r) and the degradation rate (δ) of viral strands. Our results indicate that the bifurcation associated with α generically takes place when the replication mode is closer to the SMR, thus suggesting that GR may provide viral strands with an increased robustness against degradation. This transcritical bifurcation, which is responsible for the switching from an active to an absorbing regime, suggests a smooth (i.e. second-order), absorbing-state phase transition. Finally, we also analyse a simplified model that only incorporates asymmetry in replication tied to differential replication modes.

Lethal mutants and truncated selection together solve a paradox of the origin of life

Background

Many attempts have been made to describe the origin of life, one of which is Eigen's cycle of autocatalytic reactions [Eigen M (1971) Naturwissenschaften 58, 465–523], in which primordial life molecules are replicated with limited accuracy through autocatalytic reactions. For successful evolution, the information carrier (either RNA or DNA or their precursor) must be transmitted to the next generation with a minimal number of misprints. In Eigen's theory, the maximum chain length that could be maintained is restricted to nucleotides, while for the most primitive genome the length is around . This is the famous error catastrophe paradox. How to solve this puzzle is an interesting and important problem in the theory of the origin of life.

Methodology/Principal Findings

We use methods of statistical physics to solve this paradox by carefully analyzing the implications of neutral and lethal mutants, and truncated selection (i.e., when fitness is zero after a certain Hamming distance from the master sequence) for the critical chain length. While neutral mutants play an important role in evolution, they do not provide a solution to the paradox. We have found that lethal mutants and truncated selection together can solve the error catastrophe paradox. There is a principal difference between prebiotic molecule self-replication and proto-cell self-replication stages in the origin of life.

Conclusions/Significance

We have applied methods of statistical physics to make an important breakthrough in the molecular theory of the origin of life. Our results will inspire further studies on the molecular theory of the origin of life and biological evolution.

Evolution models with lethal mutations on symmetric or random fitness landscapes

We calculate the mean fitness for evolution models, when the fitness is a function of the Hamming distance from a reference sequence, and there is a probability that this fitness is nullified (Eigen model case) or tends to the negative infinity (Crow-Kimura model case). We calculate the mean fitness of these models. The mean fitness is calculated also for the random fitnesses with logarithmic-normal distribution, reasonably describing sometimes the situation with RNA viruses.