# American Institute of Mathematical Sciences

## Moran process and Wright-Fisher process favor low variability

 Department of Mathematics and Applied Mathematics, Virginia Commonwealth University, 1015 Floyd Avenue, Richmond, VA 23284-2014, USA

* Corresponding author: Jan Rychtář

Received  October 2019 Revised  May 2020 Published  August 2020

We study evolutionary dynamics in finite populations. We assume the individuals are one of two competing genotypes, $A$ or $B$. The genotypes have the same average fitness but different variances and/or third central moments. We focus on two frequency-independent stochastic processes: (1) Wright-Fisher process and (2) Moran process. Both processes have two absorbing states corresponding to homogeneous populations of all $A$ or all $B$. Despite the fact that types $A$ and $B$ have the same average fitness, both stochastic dynamics differ from a random drift. In both processes, the selection favors $A$ replacing $B$ and opposes $B$ replacing $A$ if the fitness variance for $A$ is smaller than the fitness variance for $B$. In the case the variances are equal, the selection favors $A$ replacing $B$ and opposes $B$ replacing $A$ if the third central moment of $A$ is larger than the third central moment of $B$. We show that these results extend to structured populations and other dynamics where the selection acts at birth. We also demonstrate that the selection favors a larger variance in fitness if the selection acts at death.

Citation: Jan Rychtář, Dewey T. Taylor. Moran process and Wright-Fisher process favor low variability. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2020242
##### References:

show all references

##### References:
Fixation probabilities for the Wright-Fisher process (left) and the Moran process (right). Darker areas correspond to higher fixation probabilities; $0.1$ corresponds to a random drift. For $N = 10$, we considered genotypes $G_i; i = 1, \ldots, 11$ that have a fitness $i+1$ with probability $1/i$ and $1$ with probability $1-1/i$. The average fitness of $G_i$ is $2$, the variance is $i-1$. For every pair of $i,j$, we run $10^5$ simulations starting with a single $G_i$ individual among $G_j$ individuals. Note that the individuals' fitness is sometimes more than double the expected value
Left: A graphical representation of the genotypes $G_p$ for $p\in\{0.1, 0.2, \ldots, 0.9\}$. A fitness of $G_p$ is $\mu-x_p$ with probability $p$ and $\mu+y_p$ with probability $1-p$ where $\mu = 4$, $x_p = (1-p)y_p/p$ and $y_p = \sqrt{((1-p)^2/p + (1-p))^{-1}}$. The average fitness of $G_i$ is $\mu$ represented by the horizontal dotted line, the variance is $1$. The genotypes are color coded by $p$. The center of a disc corresponds to the fitness, the area of a disc corresponds to the probability of attaining such a fitness. The thick black curve is the third central moment $s_p$ of the genotype $G_p$ with values of the right $y$ axis. Right: Fixation probabilities for the Moran process. The darker the color, the larger the fixation probability; $0.2$ corresponds to a random drift. For $N = 5$ and every pair of $p_i,p_j$, we run $10^7$ simulations starting with a single $G_{p_i}$ individual among $G_{p_j}$ individuals
Fixation probabilities for the Db-Moran process. At every step, an individual is selected to die with the probability inversely proportional to their fitness and is then replaced by a copy of a randomly selected remaining individuals. For $N = 10$, we considered genotypes $G_i; i = 1, \ldots, 11$ as in Figure 1; a fitness of $G_i$ is $i+1$ with probability $1/i$ and $1$ with probability $1-1/i$. The average fitness of $G_i$ is $2$, the variance is $i-1$. For every pair of $i,j$, we run $10^5$ simulations starting with a single $G_i$ individual among $G_j$ individuals
Errors in approximation of $1/y$ by Taylor polynomials $P_n(y)$ at $y_0$. The approximation by a higher degree polynomial is better only on $(0,2y_0)$ and worse on $(2y_0,\infty)$
The mean square errors (MSE) of the estimate of $\mathbb{E}[\pi_A'|\pi_A]-\pi_A$ by $E_2$ (left) and by $E_3$ (right). The darker the color, the larger the error. For $N = 10$, we considered genotypes $G_i; i = 1, \ldots, 21$; a fitness of genotype $G_i$ is $i+1$ with probability $1/i$ and $1$ with probability $1-1/i$. The average fitness is $2$, the variance is $i-1$. For every $\pi_A\in\{1/N, \ldots, (N-1)/N\}$, we run $10^6$ simulations to estimate $\mathbb{E}[\pi_A'|\pi_A]-\pi_A$ numerically by the average different $\widehat{E}$. We then calculated the MSE as $\frac{1}{N-1}\sum_{\pi_A\in\left \{\frac1N, \ldots, \frac{N-1}{N}\right\}} \left(\widehat{E}-E_2\right)^2$
 [1] Xiaohong Li, Mingxin Sun, Zhaohua Gong, Enmin Feng. Multistage optimal control for microbial fed-batch fermentation process. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021040 [2] Meiqiao Ai, Zhimin Zhang, Wenguang Yu. First passage problems of refracted jump diffusion processes and their applications in valuing equity-linked death benefits. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021039 [3] Shangzhi Li, Shangjiang Guo. Permanence and extinction of a stochastic SIS epidemic model with three independent Brownian motions. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2693-2719. doi: 10.3934/dcdsb.2020201 [4] Juan Manuel Pastor, Javier García-Algarra, Javier Galeano, José María Iriondo, José J. Ramasco. A simple and bounded model of population dynamics for mutualistic networks. Networks & Heterogeneous Media, 2015, 10 (1) : 53-70. doi: 10.3934/nhm.2015.10.53 [5] Linlin Li, Bedreddine Ainseba. Large-time behavior of matured population in an age-structured model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2561-2580. doi: 10.3934/dcdsb.2020195 [6] Jianping Gao, Shangjiang Guo, Wenxian Shen. Persistence and time periodic positive solutions of doubly nonlocal Fisher-KPP equations in time periodic and space heterogeneous media. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2645-2676. doi: 10.3934/dcdsb.2020199 [7] Brandy Rapatski, James Yorke. Modeling HIV outbreaks: The male to female prevalence ratio in the core population. Mathematical Biosciences & Engineering, 2009, 6 (1) : 135-143. doi: 10.3934/mbe.2009.6.135 [8] Arseny Egorov. Morse coding for a Fuchsian group of finite covolume. Journal of Modern Dynamics, 2009, 3 (4) : 637-646. doi: 10.3934/jmd.2009.3.637 [9] Murat Uzunca, Ayşe Sarıaydın-Filibelioǧlu. Adaptive discontinuous galerkin finite elements for advective Allen-Cahn equation. Numerical Algebra, Control & Optimization, 2021, 11 (2) : 269-281. doi: 10.3934/naco.2020025 [10] Hakan Özadam, Ferruh Özbudak. A note on negacyclic and cyclic codes of length $p^s$ over a finite field of characteristic $p$. Advances in Mathematics of Communications, 2009, 3 (3) : 265-271. doi: 10.3934/amc.2009.3.265 [11] Lunji Song, Wenya Qi, Kaifang Liu, Qingxian Gu. A new over-penalized weak galerkin finite element method. Part Ⅱ: Elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2581-2598. doi: 10.3934/dcdsb.2020196 [12] Zengyun Wang, Jinde Cao, Zuowei Cai, Lihong Huang. Finite-time stability of impulsive differential inclusion: Applications to discontinuous impulsive neural networks. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2677-2692. doi: 10.3934/dcdsb.2020200 [13] Marita Holtmannspötter, Arnd Rösch, Boris Vexler. A priori error estimates for the space-time finite element discretization of an optimal control problem governed by a coupled linear PDE-ODE system. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021014 [14] Jan Prüss, Laurent Pujo-Menjouet, G.F. Webb, Rico Zacher. Analysis of a model for the dynamics of prions. Discrete & Continuous Dynamical Systems - B, 2006, 6 (1) : 225-235. doi: 10.3934/dcdsb.2006.6.225 [15] Johannes Kellendonk, Lorenzo Sadun. Conjugacies of model sets. Discrete & Continuous Dynamical Systems - A, 2017, 37 (7) : 3805-3830. doi: 10.3934/dcds.2017161 [16] Didier Bresch, Thierry Colin, Emmanuel Grenier, Benjamin Ribba, Olivier Saut. A viscoelastic model for avascular tumor growth. Conference Publications, 2009, 2009 (Special) : 101-108. doi: 10.3934/proc.2009.2009.101 [17] Ondrej Budáč, Michael Herrmann, Barbara Niethammer, Andrej Spielmann. On a model for mass aggregation with maximal size. Kinetic & Related Models, 2011, 4 (2) : 427-439. doi: 10.3934/krm.2011.4.427 [18] Martin Bohner, Sabrina Streipert. Optimal harvesting policy for the Beverton--Holt model. Mathematical Biosciences & Engineering, 2016, 13 (4) : 673-695. doi: 10.3934/mbe.2016014 [19] Chin-Chin Wu. Existence of traveling wavefront for discrete bistable competition model. Discrete & Continuous Dynamical Systems - B, 2011, 16 (3) : 973-984. doi: 10.3934/dcdsb.2011.16.973 [20] Michael Grinfeld, Amy Novick-Cohen. Some remarks on stability for a phase field model with memory. Discrete & Continuous Dynamical Systems - A, 2006, 15 (4) : 1089-1117. doi: 10.3934/dcds.2006.15.1089

2019 Impact Factor: 1.27