1
|
Gilpin W, Feldman MW. Cryptic selection forces and dynamic heritability in generalized phenotypic evolution. Theor Popul Biol 2018; 125:20-29. [PMID: 30528351 DOI: 10.1016/j.tpb.2018.11.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Revised: 11/10/2018] [Accepted: 11/14/2018] [Indexed: 11/26/2022]
Abstract
Individuals with different phenotypes can have widely-varying responses to natural selection, yet many classical approaches to evolutionary dynamics emphasize only how a population's average phenotype increases in fitness over time. However, recent experimental results have produced examples of populations that have multiple fitness peaks, or that experience frequency-dependence that affects the direction and strength of selection on certain individuals. Here, we extend classical fitness gradient formulations of natural selection in order to describe the dynamics of a phenotype distribution in terms of its moments-such as the mean, variance, and skewness. The number of governing equations in our model can be adjusted in order to capture different degrees of detail about the population. We compare our simplified model to direct Wright-Fisher simulations of evolution in several canonical fitness landscapes, and we find that our model provides a low-dimensional description of complex dynamics not typically explained by classical theory, such as cryptic selection forces due to selection on trait ranges, time-variation of the heritability, and nonlinear responses to stabilizing or disruptive selection due to asymmetric trait distributions. In addition to providing a framework for extending general understanding of common qualitative concepts in phenotypic evolution - such as fitness gradients, selection pressures, and heritability - our approach has practical importance for studying evolution in contexts in which genetic analysis is infeasible.
Collapse
Affiliation(s)
- William Gilpin
- Department of Applied Physics, Stanford University, Stanford, CA, United States.
| | - Marcus W Feldman
- Department of Biology, Stanford University, Stanford, CA, United States
| |
Collapse
|
2
|
Huang GR, Saakian DB, Hu CK. Accurate analytic solution of chemical master equations for gene regulation networks in a single cell. Phys Rev E 2018; 97:012412. [PMID: 29448337 DOI: 10.1103/physreve.97.012412] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2017] [Indexed: 12/21/2022]
Abstract
Studying gene regulation networks in a single cell is an important, interesting, and hot research topic of molecular biology. Such process can be described by chemical master equations (CMEs). We propose a Hamilton-Jacobi equation method with finite-size corrections to solve such CMEs accurately at the intermediate region of switching, where switching rate is comparable to fast protein production rate. We applied this approach to a model of self-regulating proteins [H. Ge et al., Phys. Rev. Lett. 114, 078101 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.078101] and found that as a parameter related to inducer concentration increases the probability of protein production changes from unimodal to bimodal, then to unimodal, consistent with phenotype switching observed in a single cell.
Collapse
Affiliation(s)
- Guan-Rong Huang
- Physics Division, National Center for Theoretical Sciences, Hsinchu 30013, Taiwan
| | - 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.,Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| | - Chin-Kun Hu
- Physics Division, National Center for Theoretical Sciences, Hsinchu 30013, Taiwan.,Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan.,Department of Systems Science, University of Shanghai for Science and Technology, Shanghai 200093, China.,Department of Physics, National Dong Hwa University, Hualien 97401, Taiwan
| |
Collapse
|
3
|
Bravetti A, Padilla P. An optimal strategy to solve the Prisoner's Dilemma. Sci Rep 2018; 8:1948. [PMID: 29386635 PMCID: PMC5792647 DOI: 10.1038/s41598-018-20426-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Accepted: 11/06/2017] [Indexed: 11/25/2022] Open
Abstract
Cooperation is a central mechanism for evolution. It consists of an individual paying a cost in order to benefit another individual. However, natural selection describes individuals as being selfish and in competition among themselves. Therefore explaining the origin of cooperation within the context of natural selection is a problem that has been puzzling researchers for a long time. In the paradigmatic case of the Prisoner's Dilemma (PD), several schemes for the evolution of cooperation have been proposed. Here we introduce an extension of the Replicator Equation (RE), called the Optimal Replicator Equation (ORE), motivated by the fact that evolution acts not only at the level of individuals of a population, but also among competing populations, and we show that this new model for natural selection directly leads to a simple and natural rule for the emergence of cooperation in the most basic version of the PD. Contrary to common belief, our results reveal that cooperation can emerge among selfish individuals because of selfishness itself: if the final reward for being part of a society is sufficiently appealing, players spontaneously decide to cooperate.
Collapse
Affiliation(s)
- Alessandro Bravetti
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, México City, 04510, Mexico.
| | - Pablo Padilla
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, México City, 04510, Mexico
- Fitzwilliam College, University of Cambridge, Storey's Way, CB3 ODG, UK
| |
Collapse
|
4
|
Saakian DB, Bratus AS, Hu CK. Crossing fitness canyons by a finite population. Phys Rev E 2017; 95:062405. [PMID: 28709354 DOI: 10.1103/physreve.95.062405] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2016] [Indexed: 11/07/2022]
Abstract
We consider the Wright-Fisher model of the finite population evolution on a fitness landscape defined in the sequence space by a path of nearly neutral mutations. We study a specific structure of the fitness landscape: One of the intermediate mutations on the mutation path results in either a large fitness value (climbing up a fitness hill) or a low fitness value (crossing a fitness canyon), the rest of the mutations besides the last one are neutral, and the last sequence has much higher fitness than any intermediate sequence. We derive analytical formulas for the first arrival time of the mutant with two point mutations. For the first arrival problem for the further mutants in the case of canyon crossing, we analytically deduce how the mean first arrival time scales with the population size and fitness difference. The location of the canyon on the path of sequences has a crucial role. If the canyon is at the beginning of the path, then it significantly prolongs the first arrival time; otherwise it just slightly changes it. Furthermore, the fitness hill at the beginning of the path strongly prolongs the arrival time period; however, the hill located near the end of the path shortens it. We optimize the first arrival time by applying a nonzero selection to the intermediate sequences. We extend our results and provide a scaling for the valley crossing time via the depth of the canyon and population size in the case of a fitness canyon at the first position. Our approach is useful for understanding some complex evolution systems, e.g., the evolution of cancer.
Collapse
Affiliation(s)
- David B Saakian
- Theoretical Physics Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
| | - Alexander S Bratus
- Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow 119992, Russia
| | - Chin-Kun Hu
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan.,National Center for Theoretical Sciences, National Tsing Hua University, Hsinchu 30013, Taiwan.,Business School, University of Shanghai for Science and Technology, Shanghai 200093, China
| |
Collapse
|
5
|
Saakian DB, Hu CK. Solution of classical evolutionary models in the limit when the diffusion approximation breaks down. Phys Rev E 2016; 94:042422. [PMID: 27841654 DOI: 10.1103/physreve.94.042422] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2016] [Indexed: 06/06/2023]
Abstract
The discrete time mathematical models of evolution (the discrete time Eigen model, the Moran model, and the Wright-Fisher model) have many applications in complex biological systems. The discrete time Eigen model rather realistically describes the serial passage experiments in biology. Nevertheless, the dynamics of the discrete time Eigen model is solved in this paper. The 90% of results in population genetics are connected with the diffusion approximation of the Wright-Fisher and Moran models. We considered the discrete time Eigen model of asexual virus evolution and the Wright-Fisher model from population genetics. We look at the logarithm of probabilities and apply the Hamilton-Jacobi equation for the models. We define exact dynamics for the population distribution for the discrete time Eigen model. For the Wright-Fisher model, we express the exact steady state solution and fixation probability via the solution of some nonlocal equation then give the series expansion of the solution via degrees of selection and mutation rates. The diffusion theories result in the zeroth order approximation in our approach. The numeric confirms that our method works in the case of strong selection, whereas the diffusion method fails there. Although the diffusion method is exact for the mean first arrival time, it provides incorrect approximation for the dynamics of the tail of distribution.
Collapse
Affiliation(s)
- David B Saakian
- Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan
- 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
- National Center for Theoretical Sciences, National Tsing Hua University, Hsinchu 30013, Taiwan
- Business School, University of Shanghai for Science and Technology, Shanghai 200093, China
| |
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
|