1
|
Solé R, Sardanyés J, Elena SF. Phase transitions in virology. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2021; 84:115901. [PMID: 34584031 DOI: 10.1088/1361-6633/ac2ab0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2021] [Accepted: 09/28/2021] [Indexed: 06/13/2023]
Abstract
Viruses have established relationships with almost every other living organism on Earth and at all levels of biological organization: from other viruses up to entire ecosystems. In most cases, they peacefully coexist with their hosts, but in most relevant cases, they parasitize them and induce diseases and pandemics, such as the AIDS and the most recent avian influenza and COVID-19 pandemic events, causing a huge impact on health, society, and economy. Viruses play an essential role in shaping the eco-evolutionary dynamics of their hosts, and have been also involved in some of the major evolutionary innovations either by working as vectors of genetic information or by being themselves coopted by the host into their genomes. Viruses can be studied at different levels of biological organization, from the molecular mechanisms of genome replication, gene expression and encapsidation, to global pandemics. All these levels are different and yet connected through the presence of threshold conditions allowing for the formation of a capsid, the loss of genetic information or epidemic spreading. These thresholds, as occurs with temperature separating phases in a liquid, define sharp qualitative types of behaviour. Thesephase transitionsare very well known in physics. They have been studied by means of simple, but powerful models able to capture their essential properties, allowing us to better understand them. Can the physics of phase transitions be an inspiration for our understanding of viral dynamics at different scales? Here we review well-known mathematical models of transition phenomena in virology. We suggest that the advantages of abstract, simplified pictures used in physics are also the key to properly understanding the origins and evolution of complexity in viruses. By means of several examples, we explore this multilevel landscape and how minimal models provide deep insights into a diverse array of problems. The relevance of these transitions in connecting dynamical patterns across scales and their evolutionary and clinical implications are outlined.
Collapse
Affiliation(s)
- Ricard Solé
- ICREA-Complex Systems Lab, Universitat Pompeu Fabra-PRBB, Dr Aiguader 80, 08003 Barcelona, Spain
- Institut de Biologia Evolutiva, CSIC-Universitat Pompeu Fabra, Passeig Maritim de la Barceloneta 37, 08003 Barcelona, Spain
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe NM 87501, United States of America
| | - Josep Sardanyés
- Centre de Recerca Matemàtica (CRM), Edifici C, Campus de Bellaterra, Cerdanyola del Vallès, 08193 Barcelona, Spain
- Dynamical Systems and Computational Virology, CSIC Associated Unit, Institute for Integrative Systems Biology (I2SysBio)-CRM, Spain
| | - Santiago F Elena
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe NM 87501, United States of America
- Evolutionary Systems Virology Lab (I2SysBio), CSIC-Universitat de València, Catedrático Agustín Escardino 9, Paterna, 46980 València, Spain
| |
Collapse
|
2
|
Soares ADA, Wardil L, Klaczko LB, Dickman R. Hidden role of mutations in the evolutionary process. Phys Rev E 2021; 104:044413. [PMID: 34781575 DOI: 10.1103/physreve.104.044413] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Accepted: 10/05/2021] [Indexed: 11/07/2022]
Abstract
Mutations not only alter allele frequencies in a genetic pool but may also determine the fate of an evolutionary process. Here we study which allele fixes in a one-step, one-way model including the wild type and two adaptive mutations. We study the effect of the four basic evolutionary mechanisms-genetic drift, natural selection, mutation, and gene flow-on mutant fixation and its kinetics. Determining which allele is more likely to fix is not simply a question of comparing fitnesses and mutation rates. For instance, if the allele of interest is less fit than the other, then not only must it have a greater mutation rate, but also its mutation rate must exceed a specific threshold for it to prevail. We find exact expressions for such conditions. Our conclusions are based on the mathematical description of two extreme but important regimes, as well as on simulations.
Collapse
Affiliation(s)
- Alexandre de Aquino Soares
- Departamento de Física, Instituto de Ciências Exatas (ICEx), Universidade Federal de Minas Gerais (UFMG), 31270-901 Belo Horizonte, Minas Gerais, Brazil
| | - Lucas Wardil
- Departamento de Física, Instituto de Ciências Exatas (ICEx), Universidade Federal de Minas Gerais (UFMG), 31270-901 Belo Horizonte, Minas Gerais, Brazil
| | - Louis Bernard Klaczko
- Departmento de Genética, Evolução, Microbiologia e Imunologia, Instituto de Biologia, Universidade Estadual de Campinas (Unicamp), C. P. 6109, 13083-970 Campinas, São Paulo, Brazil
| | - Ronald Dickman
- Departamento de Física and National Institute of Science and Technology for Complex Systems, Instituto de Ciências Exatas (ICEx), Universidade Federal de Minas Gerais (UFMG), C. P. 702, 30123-970 Belo Horizonte, Minas Gerais, Brazil
| |
Collapse
|
3
|
Abd Raman HS, Tan S, August JT, Khan AM. Dynamics of Influenza A (H5N1) virus protein sequence diversity. PeerJ 2020; 7:e7954. [PMID: 32518710 PMCID: PMC7261124 DOI: 10.7717/peerj.7954] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Accepted: 09/26/2019] [Indexed: 11/20/2022] Open
Abstract
Background Influenza A (H5N1) virus is a global concern with potential as a pandemic threat. High sequence variability of influenza A viruses is a major challenge for effective vaccine design. A continuing goal towards this is a greater understanding of influenza A (H5N1) proteome sequence diversity in the context of the immune system (antigenic diversity), the dynamics of mutation, and effective strategies to overcome the diversity for vaccine design. Methods Herein, we report a comprehensive study of the dynamics of H5N1 mutations by analysis of the aligned overlapping nonamer positions (1–9, 2–10, etc.) of more than 13,000 protein sequences of avian and human influenza A (H5N1) viruses, reported over at least 50 years. Entropy calculations were performed on 9,408 overlapping nonamer position of the proteome to study the diversity in the context of immune system. The nonamers represent the predominant length of the binding cores for peptides recognized by the cellular immune system. To further dissect the sequence diversity, each overlapping nonamer position was quantitatively analyzed for four patterns of sequence diversity motifs: index, major, minor and unique. Results Almost all of the aligned overlapping nonamer positions of each viral proteome exhibited variants (major, minor, and unique) to the predominant index sequence. Each variant motif displayed a characteristic pattern of incidence change in relation to increased total variants. The major variant exhibited a restrictive pyramidal incidence pattern, with peak incidence at 50% total variants. Post this peak incidence, the minor variants became the predominant motif for majority of the positions. Unique variants, each sequence observed only once, were present at nearly all of the nonamer positions. The diversity motifs (index and variants) demonstrated complex inter-relationships, with motif switching being a common phenomenon. Additionally, 25 highly conserved sequences were identified to be shared across viruses of both hosts, with half conserved to several other influenza A subtypes. Discussion The presence of distinct sequences (nonatypes) at nearly all nonamer positions represents a large repertoire of reported viral variants in the proteome, which influence the variability dynamics of the viral population. This work elucidated and provided important insights on the components that make up the viral diversity, delineating inherent patterns in the organization of sequence changes that function in the viral fitness-selection. Additionally, it provides a catalogue of all the mutational changes involved in the dynamics of H5N1 viral diversity for both avian and human host populations. This work provides data relevant for the design of prophylactics and therapeutics that overcome the diversity of the virus, and can aid in the surveillance of existing and future strains of influenza viruses.
Collapse
Affiliation(s)
| | - Swan Tan
- School of Data Sciences, Perdana University, Serdang, Selangor, Malaysia.,Institute for Immunology and Informatics, University of Rhode Island, Providence, RI, United States of America
| | - Joseph Thomas August
- School of Medicine, Johns Hopkins University, Baltimore, MD, United States of America
| | - Asif M Khan
- School of Data Sciences, Perdana University, Serdang, Selangor, Malaysia.,School of Medicine, Johns Hopkins University, Baltimore, MD, United States of America.,Beykoz Institute of Life Sciences and Biotechnology, Bezmialem Vakif University, Beykoz, Istanbul, Turkey
| |
Collapse
|
4
|
Saakian DB, Cheong KH. Solution of the Crow-Kimura model with a periodically changing (two-season) fitness function. Phys Rev E 2020; 100:022403. [PMID: 31574715 DOI: 10.1103/physreve.100.022403] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2019] [Indexed: 11/07/2022]
Abstract
Since the origin of life, both evolutionary dynamics and rhythms have played a key role in the functioning of living systems. The Crow-Kimura model of periodically changing fitness function has been solved exactly, using integral equation with time-ordered exponent. We also found a simple approximate solution for the two-season case. The evolutionary dynamics accompanied by the rhythms provide important insights into the properties of certain biological systems and processes.
Collapse
Affiliation(s)
- David B Saakian
- Laboratory of Applied Physics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Kang Hao Cheong
- Science and Math Cluster, Singapore University of Technology and Design (SUTD), S487372 Singapore.,SUTD-MIT International Design Centre, S487372, Singapore
| |
Collapse
|
5
|
Abstract
We consider the model of asexual evolution with migration, which was proposed by Waclaw et al. [Phys. Rev. Lett. 105, 268101 (2010)PRLTAO0031-900710.1103/PhysRevLett.105.268101]. This model setting is based on the standard mutation scheme from the quasispecies theory but with replicators moving from one habitat to another. The primary goal is to solve exactly the infinite population-genome length version of the model for the independent random distribution of fitnesses considered in the original paper. Moreover, we propose the analytical solution for the single peak and the symmetric fitness landscape. Our analytical solution is exact at the limit of large N. We found two phases-the correlated phase with the identical distributions of mutations in both habitats and the uncorrelated phase where the second habitat is choosing another peak of distribution in sequence space compared to the peak in the first habitat.
Collapse
Affiliation(s)
- David B Saakian
- Theoretical Physics Research Group, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam; Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam; and A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, 2 Alikhanian Brothers Street, Yerevan 375036, Armenia
| |
Collapse
|
6
|
Saakian DB, Hu CK. Mathematical Models of Quasi-Species Theory and Exact Results for the Dynamics. Curr Top Microbiol Immunol 2016; 392:121-39. [PMID: 26342705 DOI: 10.1007/82_2015_471] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
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.
Collapse
Affiliation(s)
- David B Saakian
- A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, 2 Alikhanian Brothers St., Yerevan, 375036, Armenia.
| | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei, 11529, Taiwan
| |
Collapse
|
7
|
Park JM, Niestemski LR, Deem MW. Quasispecies theory for evolution of modularity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:012714. [PMID: 25679649 PMCID: PMC4477872 DOI: 10.1103/physreve.91.012714] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2014] [Indexed: 06/04/2023]
Abstract
Biological systems are modular, and this modularity evolves over time and in different environments. A number of observations have been made of increased modularity in biological systems under increased environmental pressure. We here develop a quasispecies theory for the dynamics of modularity in populations of these systems. We show how the steady-state fitness in a randomly changing environment can be computed. We derive a fluctuation dissipation relation for the rate of change of modularity and use it to derive a relationship between rate of environmental changes and rate of growth of modularity. We also find a principle of least action for the evolved modularity at steady state. Finally, we compare our predictions to simulations of protein evolution and find them to be consistent.
Collapse
Affiliation(s)
- Jeong-Man Park
- Departments of Physics & Astronomy and Bioengineering, Rice University, Houston, Texas 77005-1892, USA; Department of Physical and Biological Science, Western New England University, Springfield, Massachusetts 01119, USA; and Department of Physics, The Catholic University of Korea, Bucheon 420-743, Korea
| | - Liang Ren Niestemski
- Departments of Physics & Astronomy and Bioengineering, Rice University, Houston, Texas 77005-1892, USA; Department of Physical and Biological Science, Western New England University, Springfield, Massachusetts 01119, USA; and Department of Physics, The Catholic University of Korea, Bucheon 420-743, Korea
| | - Michael W Deem
- Departments of Physics & Astronomy and Bioengineering, Rice University, Houston, Texas 77005-1892, USA; Department of Physical and Biological Science, Western New England University, Springfield, Massachusetts 01119, USA; and Department of Physics, The Catholic University of Korea, Bucheon 420-743, Korea
| |
Collapse
|
8
|
Abstract
The quasispecies concept is introduced by means of a simple theoretical model that uses as little chemical kinetics and mathematics as possible but fully in the spirit of Albert Einstein who said: "Things should be made as simple as possible but not simpler." More elaborate treatments follow in the forthcoming chapters. It is shown that the most important results of the theory, in particular the existence of error thresholds, are not dependent on simplifying assumptions concerning the distribution of fitness values. Error thresholds are regularly found on landscapes with large and irregular scatter of fitness. After the introduction to theory, it will be shown how experimental data on the evolution of molecules or viruses may be fit to the theoretical model.
Collapse
|
9
|
Getting to Know Viral Evolutionary Strategies: Towards the Next Generation of Quasispecies Models. Curr Top Microbiol Immunol 2015; 392:201-17. [PMID: 26271604 DOI: 10.1007/82_2015_457] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Viral populations are formed by complex ensembles of genomes with broad phenotypic diversity. The adaptive strategies deployed by these ensembles are multiple and often cannot be predicted a priori. Our understanding of viral dynamics is mostly based on two kinds of empirical approaches: one directed towards characterizing molecular changes underlying fitness changes and another focused on population-level responses. Simultaneously, theoretical efforts are directed towards developing a formal picture of viral evolution by means of more realistic fitness landscapes and reliable population dynamics models. New technologies, chiefly the use of next-generation sequencing and related tools, are opening avenues connecting the molecular and the population levels. In the near future, we hope to be witnesses of an integration of these still decoupled approaches, leading into more accurate and realistic quasispecies models able to capture robust generalities and endowed with a satisfactory predictive power.
Collapse
|
10
|
Saakian DB, Ghazaryan MH, Hu CK. Punctuated equilibrium and shock waves in molecular models of biological evolution. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:022712. [PMID: 25215763 DOI: 10.1103/physreve.90.022712] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2013] [Indexed: 06/03/2023]
Abstract
We consider the dynamics in infinite population evolution models with a general symmetric fitness landscape. We find shock waves, i.e., discontinuous transitions in the mean fitness, in evolution dynamics even with smooth fitness landscapes, which means that the search for the optimal evolution trajectory is more complicated. These shock waves appear in the case of positive epistasis and can be used to represent punctuated equilibria in biological evolution during long geological time scales. We find exact analytical solutions for discontinuous dynamics at the large-genome-length limit and derive optimal mutation rates for a fixed fitness landscape to send the population from the initial configuration to some final configuration in the fastest way.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan; Yerevan Physics Institute, Alikhanian Brothers Str. 2, Yerevan 375036, Armenia; and National Center for Theoretical Sciences: Physics Division, National Taiwan University, Taipei 10617, Taiwan
| | | | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| |
Collapse
|
11
|
Amor DR, Solé RV. Catastrophic shifts and lethal thresholds in a propagating front model of unstable tumor progression. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:022710. [PMID: 25215761 DOI: 10.1103/physreve.90.022710] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2014] [Indexed: 06/03/2023]
Abstract
Unstable dynamics characterizes the evolution of most solid tumors. Because of an increased failure of maintaining genome integrity, a cumulative increase in the levels of gene mutation and loss is observed. Previous work suggests that instability thresholds to cancer progression exist, defining phase transition phenomena separating tumor-winning scenarios from tumor extinction or coexistence phases. Here we present an integral equation approach to the quasispecies dynamics of unstable cancer. The model exhibits two main phases, characterized by either the success or failure of cancer tissue. Moreover, the model predicts that tumor failure can be due to either a reduced selective advantage over healthy cells or excessive instability. We also derive an approximate, analytical solution that predicts the front speed of aggressive tumor populations on the instability space.
Collapse
Affiliation(s)
- Daniel R Amor
- ICREA-Complex Systems Lab, Universitat Pompeu Fabra, Dr. Aiguader 80, 08003 Barcelona, Spain and Institut de Biologia Evolutiva, CSIC-UPF, Psg Barceloneta, Barcelona, Spain
| | - Ricard V Solé
- ICREA-Complex Systems Lab, Universitat Pompeu Fabra, Dr. Aiguader 80, 08003 Barcelona, Spain and Institut de Biologia Evolutiva, CSIC-UPF, Psg Barceloneta, Barcelona, Spain and Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
12
|
Ancliff M, Park JM. Evolution dynamics of a model for gene duplication under adaptive conflict. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:062702. [PMID: 25019809 DOI: 10.1103/physreve.89.062702] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2014] [Indexed: 06/03/2023]
Abstract
We present and solve the dynamics of a model for gene duplication showing escape from adaptive conflict. We use a Crow-Kimura quasispecies model of evolution where the fitness landscape is a function of Hamming distances from two reference sequences, which are assumed to optimize two different gene functions, to describe the dynamics of a mixed population of individuals with single and double copies of a pleiotropic gene. The evolution equations are solved through a spin coherent state path integral, and we find two phases: one is an escape from an adaptive conflict phase, where each copy of a duplicated gene evolves toward subfunctionalization, and the other is a duplication loss of function phase, where one copy maintains its pleiotropic form and the other copy undergoes neutral mutation. The phase is determined by a competition between the fitness benefits of subfunctionalization and the greater mutational load associated with maintaining two gene copies. In the escape phase, we find a dynamics of an initial population of single gene sequences only which escape adaptive conflict through gene duplication and find that there are two time regimes: until a time t single gene sequences dominate, and after t double gene sequences outgrow single gene sequences. The time t is identified as the time necessary for subfunctionalization to evolve and spread throughout the double gene sequences, and we show that there is an optimum mutation rate which minimizes this time scale.
Collapse
Affiliation(s)
- Mark Ancliff
- Department of Physics, The Catholic University of Korea, Bucheon 420-743, Korea
| | - Jeong-Man Park
- Department of Physics, The Catholic University of Korea, Bucheon 420-743, Korea
| |
Collapse
|
13
|
Variability in mutational fitness effects prevents full lethal transitions in large quasispecies populations. Sci Rep 2014; 4:4625. [PMID: 24713667 PMCID: PMC3980229 DOI: 10.1038/srep04625] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2014] [Accepted: 03/21/2014] [Indexed: 11/08/2022] Open
Abstract
The distribution of mutational fitness effects (DMFE) is crucial to the evolutionary fate of quasispecies. In this article we analyze the effect of the DMFE on the dynamics of a large quasispecies by means of a phenotypic version of the classic Eigen's model that incorporates beneficial, neutral, deleterious, and lethal mutations. By parameterizing the model with available experimental data on the DMFE of Vesicular stomatitis virus (VSV) and Tobacco etch virus (TEV), we found that increasing mutation does not totally push the entire viral quasispecies towards deleterious or lethal regions of the phenotypic sequence space. The probability of finding regions in the parameter space of the general model that results in a quasispecies only composed by lethal phenotypes is extremely small at equilibrium and in transient times. The implications of our findings can be extended to other scenarios, such as lethal mutagenesis or genomically unstable cancer, where increased mutagenesis has been suggested as a potential therapy.
Collapse
|
14
|
Stationary distribution of self-organized states and biological information generation. Sci Rep 2013; 3:3329. [PMID: 24281357 PMCID: PMC3839033 DOI: 10.1038/srep03329] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2013] [Accepted: 11/08/2013] [Indexed: 11/08/2022] Open
Abstract
Self-organization, where spontaneous orderings occur under driven conditions, is one of the hallmarks of biological systems. We consider a statistical mechanical treatment of the biased distribution of such organized states, which become favored as a result of their catalytic activity under chemical driving forces. A generalization of the equilibrium canonical distribution describes the stationary state, which can be used to model shifts in conformational ensembles sampled by an enzyme in working conditions. The basic idea is applied to the process of biological information generation from random sequences of heteropolymers, where unfavorable Shannon entropy is overcome by the catalytic activities of selected genes. The ordering process is demonstrated with the genetic distance to a genotype with high catalytic activity as an order parameter. The resulting free energy can have multiple minima, corresponding to disordered and organized phases with first-order transitions between them.
Collapse
|
15
|
Saakian DB, Hu CK. Evolutionary advantage via common action of recombination and neutrality. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:052717. [PMID: 24329306 DOI: 10.1103/physreve.88.052717] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2012] [Revised: 07/09/2013] [Indexed: 06/03/2023]
Abstract
We investigate evolution models with recombination and neutrality. We consider the Crow-Kimura (parallel) mutation-selection model with the neutral fitness landscape, in which there is a central peak with high fitness A, and some of 1-point mutants have the same high fitness A, while the fitness of other sequences is 0. We find that the effect of recombination and neutrality depends on the concrete version of both neutrality and recombination. We consider three versions of neutrality: (a) all the nearest neighbor sequences of the peak sequence have the same high fitness A; (b) all the l-point mutations in a piece of genome of length l≥1 are neutral; (c) the neutral sequences are randomly distributed among the nearest neighbors of the peak sequences. We also consider three versions of recombination: (I) the simple horizontal gene transfer (HGT) of one nucleotide; (II) the exchange of a piece of genome of length l, HGT-l; (III) two-point crossover recombination (2CR). For the case of (a), the 2CR gives a rather strong contribution to the mean fitness, much stronger than that of HGT for a large genome length L. For the random distribution of neutral sequences there is a critical degree of neutrality ν(c), and for μ<μ(c) and (μ(c)-μ) is not large, the 2CR suppresses the mean fitness while HGT increases it; for ν much larger than ν(c), the 2CR and HGT-l increase the mean fitness larger than that of the HGT. We also consider the recombination in the case of smooth fitness landscapes. The recombination gives some advantage in the evolutionary dynamics, where recombination distinguishes clearly the mean-field-like evolutionary factors from the fluctuation-like ones. By contrast, mutations affect the mean-field-like and fluctuation-like factors similarly. Consequently, recombination can accelerate the non-mean-field (fluctuation) type dynamics without considerably affecting the mean-field-like factors.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan and A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, 2 Alikhanian Brothers Street, Yerevan 375036, Armenia and Physics Division of National Center for Theoretical Sciences Taipei Branch, National Taiwan University, Taipei, Taiwan
| | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| |
Collapse
|
16
|
Baronchelli A, Chater N, Christiansen MH, Pastor-Satorras R. Evolution in a changing environment. PLoS One 2013; 8:e52742. [PMID: 23326355 PMCID: PMC3542356 DOI: 10.1371/journal.pone.0052742] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2012] [Accepted: 11/16/2012] [Indexed: 11/18/2022] Open
Abstract
We propose a simple model for genetic adaptation to a changing environment, describing a fitness landscape characterized by two maxima. One is associated with "specialist" individuals that are adapted to the environment; this maximum moves over time as the environment changes. The other maximum is static, and represents "generalist" individuals not affected by environmental changes. The rest of the landscape is occupied by "maladapted" individuals. Our analysis considers the evolution of these three subpopulations. Our main result is that, in presence of a sufficiently stable environmental feature, as in the case of an unchanging aspect of a physical habitat, specialists can dominate the population. By contrast, rapidly changing environmental features, such as language or cultural habits, are a moving target for the genes; here, generalists dominate, because the best evolutionary strategy is to adopt neutral alleles not specialized for any specific environment. The model we propose is based on simple assumptions about evolutionary dynamics and describes all possible scenarios in a non-trivial phase diagram. The approach provides a general framework to address such fundamental issues as the Baldwin effect, the biological basis for language, or the ecological consequences of a rapid climate change.
Collapse
Affiliation(s)
- Andrea Baronchelli
- Laboratory for the Modeling of Biological and Socio-technical Systems, Northeastern University, Boston, Massachusetts, United States of America
| | - Nick Chater
- Behavioural Science Group, Warwick Business School, University of Warwick, Coventry, United Kingdom
| | - Morten H. Christiansen
- Department of Psychology, Cornell University, Ithaca, New York, United States of America
- Santa Fe Institute, Santa Fe, New Mexico, United States of America
| | - Romualdo Pastor-Satorras
- Departament de Fisica i Enginyeria Nuclear, Universitat Politecnica de Catalunya, Barcelona, Spain
- * E-mail:
| |
Collapse
|
17
|
Saakian DB, Kirakosyan Z, Hu CK. Biological evolution in a multidimensional fitness landscape. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:031920. [PMID: 23030957 DOI: 10.1103/physreve.86.031920] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2012] [Revised: 05/07/2012] [Indexed: 06/01/2023]
Abstract
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.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan.
| | | | | |
Collapse
|
18
|
Abstract
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.
Collapse
Affiliation(s)
- Esteban Domingo
- Centro de Biología Molecular Severo Ochoa (CSIC-UAM), C/ Nicolás Cabrera, Universidad Autónoma de Madrid, Cantoblanco, Madrid, Spain.
| | | | | |
Collapse
|
19
|
Quasispecies as a matter of fact: viruses and beyond. Virus Res 2011; 162:203-15. [PMID: 21945638 PMCID: PMC7172439 DOI: 10.1016/j.virusres.2011.09.018] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2011] [Revised: 09/12/2011] [Accepted: 09/12/2011] [Indexed: 12/13/2022]
Abstract
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.
Collapse
|
20
|
Guttenberg N, Tabei SMA, Dinner AR. Short-time evolution in the adaptive immune system. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:031932. [PMID: 22060428 DOI: 10.1103/physreve.84.031932] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2011] [Revised: 05/25/2011] [Indexed: 05/31/2023]
Abstract
We exploit a simple model to numerically and analytically investigate the effect of enforcing a time constraint for achieving a system-wide goal during an evolutionary dynamics. This situation is relevant to finding antibody specificities in the adaptive immune response as well as to artificial situations in which an evolutionary dynamics is used to generate a desired capability in a limited number of generations. When the likelihood of finding the target phenotype is low, we find that the optimal mutation rate can exceed the error threshold, in contrast to conventional evolutionary dynamics. We also show how a logarithmic correction to the usual inverse scaling of population size with mutation rate arises. Implications for natural and artificial evolutionary situations are discussed.
Collapse
Affiliation(s)
- Nicholas Guttenberg
- James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | | | | |
Collapse
|
21
|
Saakian DB, Biebricher CK, Hu CK. Lethal mutants and truncated selection together solve a paradox of the origin of life. PLoS One 2011; 6:e21904. [PMID: 21814563 PMCID: PMC3144202 DOI: 10.1371/journal.pone.0021904] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2010] [Accepted: 06/14/2011] [Indexed: 11/25/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- David B. Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei, Taiwan
- Yerevan Physics Institute, Yerevan, Armenia
- National Center for Theoretical Sciences: Physics Division, National Taiwan University, Taipei, Taiwan
| | | | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei, Taiwan
- Department of Physics, Beijing Normal University, Beijing, China
- * E-mail:
| |
Collapse
|
22
|
Ancliff M, Park JM. Optimal mutation rates in dynamic environments: The Eigen model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:021904. [PMID: 20866834 DOI: 10.1103/physreve.82.021904] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2010] [Indexed: 05/29/2023]
Abstract
We consider the Eigen quasispecies model with a dynamic environment. For an environment with sharp-peak fitness in which the most-fit sequence moves by k spin-flips each period T we find an asymptotic stationary state in which the quasispecies population changes regularly according to the regular environmental change. From this stationary state we estimate the maximum and the minimum mutation rates for a quasispecies to survive under the changing environment and calculate the optimum mutation rate that maximizes the population growth. Interestingly we find that the optimum mutation rate in the Eigen model is lower than that in the Crow-Kimura model, and at their optimum mutation rates the corresponding mean fitness in the eigenmodel is lower than that in the Crow-Kimura model, suggesting that the mutation process which occurs in parallel to the replication process as in the Crow-Kimura model gives an adaptive advantage under changing environment.
Collapse
Affiliation(s)
- Mark Ancliff
- Department of Physics, The Catholic University of Korea, Bucheon, Gyeonggi-do, Korea
| | | |
Collapse
|
23
|
Saakian DB, Hu CK. Selection via flatness as a dynamical effect in evolution models with finite population. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:011902. [PMID: 20866643 DOI: 10.1103/physreve.82.011902] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2010] [Revised: 06/08/2010] [Indexed: 05/29/2023]
Abstract
We investigate the phenomenon of selection via flatness. In the static case, the finiteness of the population does not seriously influence the increase of mean fitness of population due to flatness around a peak. The effect is proportional to 1/square root(L), where L is the genome length. We investigated the two peak model (high peak and a flat peak). We find that the selection of flatness for long genome lengths occurs as a dynamic phenomenon in the case of evolution with small populations. We found that two factors are crucial to define the role of flatness: special initial distribution (the population is located at centers of peaks) allows flat peak to attract more population, and the large value of mutations per population per virus life cycle sometimes also increases the role of flatness. We suggested simple criteria to identify the phenomenon of dynamical arresting of population around flat peak by experiment. We infer that selection via robustness is possible in evolution as a nonequilibrium phenomenon.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan.
| | | |
Collapse
|
24
|
Kirakosyan Z, Saakian DB, Hu CK. Evolution models with lethal mutations on symmetric or random fitness landscapes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:011904. [PMID: 20866645 DOI: 10.1103/physreve.82.011904] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/01/2010] [Indexed: 05/29/2023]
Abstract
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.
Collapse
Affiliation(s)
- Zara Kirakosyan
- Yerevan Physics Institute, Alikhanian Brothers Street 2, Yerevan 375036, Armenia
| | | | | |
Collapse
|
25
|
Saakian DB, Martirosyan AS, Hu CK. Different fitnesses for in vivo and in vitro evolutions due to the finite generation-time effect. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:061913. [PMID: 20866446 DOI: 10.1103/physreve.81.061913] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2009] [Revised: 04/07/2010] [Indexed: 05/29/2023]
Abstract
We consider the finite generation-time effect in virus evolution models, introducing differential equations with delay. The suggested approach more adequately describes the evolution in case of growing populations than the popular models of population genetics, especially for the viruses with large number of offspring during one life cycle. Now the mean fitness, as a coefficient for exponential population growth, could not be defined via instant characteristics of the model. For the constant population size the finite generation-time does not affect mean fitness in the steady state. The growing virus population is characterized by different fitness than the population with a constant size.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| | | | | |
Collapse
|
26
|
Avetisyan Z, Saakian DB. Recombination in one- and two-dimensional fitness landscapes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:051916. [PMID: 20866270 DOI: 10.1103/physreve.81.051916] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2009] [Indexed: 05/29/2023]
Abstract
We consider many-site mutation-recombination models of molecular evolution, where fitness is a function of a Hamming distance from one (one-dimensional case) or two (two-dimensional case) sequences. For the one-dimensional case, we calculate the population distribution dynamics for a model with zero fitness and an arbitrary symmetric initial distribution and find an error threshold transition point in the single-peak fitness model for a given initial symmetric distribution. We calculate the recombination period in the case of a single-peak fitness function, when the original population is located at one sequence, at some Hamming distance from the peak configuration. Steady-state fitness is calculated with finite genome length corrections. We derive analytical equations for the two-dimensional mutation-recombination model.
Collapse
Affiliation(s)
- Zh Avetisyan
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany
| | | |
Collapse
|
27
|
Park JM, Muñoz E, Deem MW. Quasispecies theory for finite populations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:011902. [PMID: 20365394 PMCID: PMC4479305 DOI: 10.1103/physreve.81.011902] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2008] [Revised: 10/26/2009] [Indexed: 05/25/2023]
Abstract
We present stochastic, finite-population formulations of the Crow-Kimura and Eigen models of quasispecies theory, for fitness functions that depend in an arbitrary way on the number of mutations from the wild type. We include back mutations in our description. We show that the fluctuation of the population numbers about the average values is exceedingly large in these physical models of evolution. We further show that horizontal gene transfer reduces by orders of magnitude the fluctuations in the population numbers and reduces the accumulation of deleterious mutations in the finite population due to Muller's ratchet. Indeed, the population sizes needed to converge to the infinite population limit are often larger than those found in nature for smooth fitness functions in the absence of horizontal gene transfer. These analytical results are derived for the steady state by means of a field-theoretic representation. Numerical results are presented that indicate horizontal gene transfer speeds up the dynamics of evolution as well.
Collapse
Affiliation(s)
- Jeong-Man Park
- Department of Physics & Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | | | | |
Collapse
|
28
|
Abstract
The immune system recognizes a myriad of invading pathogens and their toxic products. It does so with a finite repertoire of antibodies and T cell receptors. We here describe theories that quantify the dynamics of the immune system. We describe how the immune system recognizes antigens by searching the large space of receptor molecules. We consider in some detail the theories that quantify the immune response to influenza and dengue fever. We review theoretical descriptions of the complementary evolution of pathogens that occurs in response to immune system pressure. Methods including bioinformatics, molecular simulation, random energy models, and quantum field theory contribute to a theoretical understanding of aspects of immunity.
Collapse
Affiliation(s)
- Michael W Deem
- Department of Bioengineering and Physics, Rice University, Houston, TX 77005, USA.
| | | |
Collapse
|
29
|
Ancliff M, Park JM. Maximum, minimum, and optimal mutation rates in dynamic environments. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:061910. [PMID: 20365193 DOI: 10.1103/physreve.80.061910] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2009] [Indexed: 05/29/2023]
Abstract
We analyze the dynamics of the parallel mutation-selection quasispecies model with a changing environment. For an environment with the sharp-peak fitness function in which the most fit sequence changes by k spin flips every period T , we find analytical expressions for the minimum and maximum mutation rates for which a quasispecies can survive, valid in the limit of large sequence size. We find an asymptotic solution in which the quasispecies population changes periodically according to the periodic environmental change. In this state we compute the mutation rate that gives the optimal mean fitness over a period. We find that the optimal mutation rate per genome, k/T , is independent of genome size, a relationship which is observed across broad groups of real organisms.
Collapse
Affiliation(s)
- Mark Ancliff
- Department of Physics, The Catholic University of Korea, Bucheon, Gyeonggi-do 420-743, Korea
| | | |
Collapse
|
30
|
Lethal mutagenesis in a structured environment. J Theor Biol 2009; 261:67-73. [PMID: 19627995 DOI: 10.1016/j.jtbi.2009.07.014] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2009] [Revised: 07/02/2009] [Accepted: 07/08/2009] [Indexed: 02/05/2023]
Abstract
We analyze how lethal mutagenesis operates in a compartmentalized host. We assume that different compartments receive different amounts of mutagen and that virions can migrate among compartments. We address two main questions: (1) To what extent can refugia, i.e., compartments that receive little mutagen, prevent extinction? (2) Does migration among compartments limit the effectiveness of refugia? We find that if there is little migration, extinction has to be achieved separately in all compartments. In this case, the total dose of mutagen administered to the host needs to be so high that the mutagen is effective even in the refugia. By contrast, if migration is extensive, then lethal mutagenesis is effective as long as the average growth in all compartments is reduced to below replacement levels. The effectiveness of migration is governed by the ratio of virion replication and death rates, R(0). The smaller R(0), the less migration is necessary to neutralize refugia and the less mutagen is necessary to achieve extinction at high migration rates.
Collapse
|
31
|
Barton N, Coe J. On the application of statistical physics to evolutionary biology. J Theor Biol 2009; 259:317-24. [DOI: 10.1016/j.jtbi.2009.03.019] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2008] [Revised: 03/07/2009] [Accepted: 03/10/2009] [Indexed: 11/16/2022]
|
32
|
Saakian DB, Biebricher CK, Hu CK. Phase diagram for the Eigen quasispecies theory with a truncated fitness landscape. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:041905. [PMID: 19518254 DOI: 10.1103/physreve.79.041905] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2008] [Revised: 11/24/2008] [Indexed: 05/27/2023]
Abstract
Using methods of statistical physics, we present rigorous theoretical calculations of Eigen's quasispecies theory with the truncated fitness landscape which dramatically limits the available sequence space of information carriers. As the mutation rate is increased from small values to large values, one can observe three phases: the first (I) selective (also known as ferromagnetic) phase, the second (II) intermediate phase with some residual order, and the third (III) completely randomized (also known as paramagnetic) phase. We calculate the phase diagram for these phases and the concentration of information carriers in the master sequence (also known as peak configuration) x0 and other classes of information carriers. As the phase point moves across the boundary between phase I and phase II, x0 changes continuously; as the phase point moves across the boundary between phase II and phase III, x0 has a large change. Our results are applicable for the general case of a fitness landscape.
Collapse
Affiliation(s)
- David B Saakian
- Yerevan Physics Institute, Alikhanian Brothers St. 2, Yerevan 375036, Armenia
| | | | | |
Collapse
|
33
|
Muñoz E, Park JM, Deem MW. Quasispecies theory for horizontal gene transfer and recombination. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:061921. [PMID: 19256882 PMCID: PMC4478466 DOI: 10.1103/physreve.78.061921] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/03/2008] [Revised: 10/03/2008] [Indexed: 05/27/2023]
Abstract
We introduce a generalization of the parallel, or Crow-Kimura, and Eigen models of molecular evolution to represent the exchange of genetic information between individuals in a population. We study the effect of different schemes of genetic recombination on the steady-state mean fitness and distribution of individuals in the population, through an analytic field theoretic mapping. We investigate both horizontal gene transfer from a population and recombination between pairs of individuals. Somewhat surprisingly, these nonlinear generalizations of quasispecies theory to modern biology are analytically solvable. For two-parent recombination, we find two selected phases, one of which is spectrally rigid. We present exact analytical formulas for the equilibrium mean fitness of the population, in terms of a maximum principle, which are generally applicable to any permutation invariant replication rate function. For smooth fitness landscapes, we show that when positive epistatic interactions are present, recombination or horizontal gene transfer introduces a mild load against selection. Conversely, if the fitness landscape exhibits negative epistasis, horizontal gene transfer or recombination introduces an advantage by enhancing selection towards the fittest genotypes. These results prove that the mutational deterministic hypothesis holds for quasispecies models. For the discontinuous single sharp peak fitness landscape, we show that horizontal gene transfer has no effect on the fitness, while recombination decreases the fitness, for both the parallel and the Eigen models. We present numerical and analytical results as well as phase diagrams for the different cases.
Collapse
Affiliation(s)
- Enrique Muñoz
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | | | | |
Collapse
|
34
|
Saakian DB, Kirakosyan Z, Hu CK. Diploid biological evolution models with general smooth fitness landscapes and recombination. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:061907. [PMID: 18643300 DOI: 10.1103/physreve.77.061907] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2007] [Revised: 04/02/2008] [Indexed: 05/26/2023]
Abstract
Using a Hamilton-Jacobi equation approach, we obtain analytic equations for steady-state population distributions and mean fitness functions for Crow-Kimura and Eigen-type diploid biological evolution models with general smooth hypergeometric fitness landscapes. Our numerical solutions of diploid biological evolution models confirm the analytic equations obtained. We also study the parallel diploid model for the simple case of recombination and calculate the variance of distribution, which is consistent with numerical results.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| | | | | |
Collapse
|
35
|
Narnhofer H, Posch HA, Thirring W. Emergence of order in quantum extensions of the classical quasispecies evolution. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:041133. [PMID: 17994962 DOI: 10.1103/physreve.76.041133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2007] [Indexed: 05/25/2023]
Abstract
We study evolution equations which model selection and mutation within the framework of quantum mechanics. The main question is to what extent order is achieved for an ensemble of typical systems. As an indicator for mixing or purification, a quadratic entropy is used which assumes values between zero for pure states and (d-1)/d for fully mixed states. Here, d is the dimension. Whereas the classical counterpart, the quasispecies dynamics, has previously been found to be predominantly mixing, the quantum quasispecies (QS) evolution surprisingly is found to be strictly purifying for all dimensions. This is also typically true for an alternative formulation (AQS) of this quantum mechanical flow. We compare this also to analogous results for the Lindblad evolution. Although the latter may be viewed as a simple linear superposition of the purifying QS and AQS evolutions, it is found to be predominantly mixing. The reason for this behavior may be explained by the fact that the two subprocesses by themselves converge to different pure states, such that the combined process is mixing. These results also apply to high-dimensional systems.
Collapse
Affiliation(s)
- Heide Narnhofer
- Faculty of Physics, Universität Wien, Boltzmanngasse 5, A-1090 Wien, Austria.
| | | | | |
Collapse
|
36
|
Jain K. Evolutionary dynamics of the most populated genotype on rugged fitness landscapes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:031922. [PMID: 17930286 DOI: 10.1103/physreve.76.031922] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2007] [Indexed: 05/25/2023]
Abstract
We consider an asexual population evolving on rugged fitness landscapes which are defined on the multidimensional genotypic space and have many local optima. We track the most populated genotype as it changes when the population jumps from a fitness peak to a better one during the process of adaptation. This is done using the dynamics of the shell model which is a simplified version of the quasispecies model for infinite populations and standard Wright-Fisher dynamics for large finite populations. We show that the population fraction of a genotype obtained within the quasispecies model and the shell model match for fit genotypes and at short times, but the dynamics of the two models are identical for questions related to the most populated genotype. We calculate exactly several properties of the jumps in infinite populations, some of which were obtained numerically in previous works. We also present our preliminary simulation results for finite populations. In particular, we measure the jump distribution in time and find that it decays as t(-2) as in the quasispecies problem.
Collapse
Affiliation(s)
- Kavita Jain
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
37
|
Izmailian NS, Papoyan VV, Priezzhev VB, Hu CK. Self-organizing behavior in a lattice model for co-evolution of virus and immune systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:041104. [PMID: 17500862 DOI: 10.1103/physreve.75.041104] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2006] [Revised: 01/12/2007] [Indexed: 05/15/2023]
Abstract
We propose a lattice model for the co-evolution of a virus population and an adaptive immune system. We show that, under some natural assumptions, both probability distribution of the virus population and the distribution of activity of the immune system tend during the evolution to a self-organized critical state.
Collapse
Affiliation(s)
- N Sh Izmailian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| | | | | | | |
Collapse
|
38
|
Park JM, Deem MW. Phase diagrams of quasispecies theory with recombination and horizontal gene transfer. PHYSICAL REVIEW LETTERS 2007; 98:058101. [PMID: 17358903 DOI: 10.1103/physrevlett.98.058101] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2006] [Indexed: 05/14/2023]
Abstract
We consider how transfer of genetic information between individuals influences the phase diagram and mean fitness of both the Eigen and the parallel, or Crow-Kimura, models of evolution. In the absence of genetic transfer, these physical models of evolution consider the replication and point mutation of the genomes of independent individuals in a large population. A phase transition occurs, such that below a critical mutation rate an identifiable quasispecies forms. We show how transfer of genetic information changes the phase diagram and mean fitness and introduces metastability in quasispecies theory, via an analytic field theoretic mapping.
Collapse
Affiliation(s)
- J-M Park
- Department of Physics & Astronomy and Department of Bioengineering, Rice University, Houston, Texas 77005-1892, USA
| | | |
Collapse
|