---
title: Moran Process in Evolutionary Dynamics
url: https://www.emergentmind.com/topics/moran-process
type: topic
---

# Moran Process in Evolutionary Dynamics

The Moran process is a cornerstone stochastic model in evolutionary dynamics, formalizing the spread of genetic mutations within finite, often structured, populations. Its canonical form and numerous extensions provide a quantitative foundation for studying fixation probabilities, evolutionary amplification and suppression, absorption times, the effects of spatial structure and heterogeneity, and key limits of stochastic population genetics, ecology, and evolutionary game theory.

## 1. Foundations and Mathematical Formulation

The classical Moran process models a population of constant size $N$ in which each individual is associated with a position (vertex) in a graph $G$; each vertex is occupied by either a mutant or a non-mutant (resident) individual. The process is typically discrete-time and defined by two principal update steps per iteration:
1. **Birth selection:** An individual is chosen with probability proportional to its fitness ($r$ for mutants, typically $1$ for residents).
2. **Replacement:** The offspring is placed on a uniformly random neighbor (or more generally, according to the transition structure of $G$), replacing the individual previously present.

Once all individuals are mutants (*fixation*) or all are residents (*extinction*), the system is absorbed. The **fixation probability** $\rho_G(r)$ is the probability that a single mutant, introduced uniformly at random, eventually takes over the entire system.

On general graphs, the Moran process is a Markov chain on $2^n$ states (for $n$ vertices), with transitions determined by both fitness parameters and the network structure. Extensions allow for different update rules (birth-death vs. death-birth), node-specific fitnesses, mutation and selection, and non-neutral, frequency-dependent, or type-specific interactions [2410.09476][1709.03031][1012.3913][1611.01585][1512.05632][1111.3321][1311.7631].


## 2. Fixation Probability, Amplification, and Suppression

The fixation probability in the classic, well-mixed (complete graph) case is 
\[
\rho_{K_n}(r) = \frac{1 - r^{-1}}{1 - r^{-n}}
\]
which reduces to $1/n$ for the neutral case ($r=1$). Graph structure dramatically alters this baseline.

**Amplifier of selection:** A graph family is an amplifier if $\rho_{G_n}(r) > \rho_{K_n}(r)$ for $r>1$. *Strong amplifiers* satisfy $\lim_{n \to \infty} \rho_{G_n}(r) = 1$ for all $r > 1$; *strong suppressors* yield $\lim_{n \to \infty} \rho_{G_n}(r) = 0$ [1611.01585][1512.05632][1611.04209]. The existence of such undirected graphs was disproved until recently. Construction of families such as *megastars* (directed), *dense incubators* and *strong amplifier/suppressor* families (undirected) demonstrate that strong amplification and suppression can arise in undirected populations under suitable topologies:
- **Strong amplifier:** $\rho_{A_n,\epsilon}(r) = 1 - \tilde{O}(n^{-1/3})$.
- **Strong suppressor:** $\rho_{B_n}(r) = \tilde{O}(n^{-1/4})$ [1611.01585].
No undirected graph exceeds the upper bound $\rho_G(r) \le 1 - \tilde\Omega(n^{-1/3})$.

For digraphs, extinction probabilities as low as $O(n^{-1/2})$ are attainable via *megastar* constructions, which are shown to be asymptotically optimal among all strongly-connected digraphs [1512.05632][1611.04209]. For undirected graphs, the *dense incubators* achieve $O(n^{-1/3})$ decay.

The detailed structure—bottlenecks, high-degree hubs, expansion properties—determines whether a network amplifies or suppresses selection [1611.01585][1512.05632][2410.09476].


## 3. Absorption and Colonization Times: Bounds and Extremal Structures

The absorption time is the expected number of steps for the system to reach either fixation or extinction. For general undirected graphs, the absorption time is polynomial:
- $O(n^4)$ for undirected graphs with $r>1$ [1311.7631][1111.3321], with improved $O(n^3 \exp(O((\log\log n)^3)))$ bound and explicit extremal examples exhibiting $\Theta(n^3)$ scaling [1804.02293]. On regular undirected graphs, $O(n^2)$ is tight.
- For regular digraphs, $O(n^2)$ upper and $\Omega(n \log n)$ lower bounds hold.

In directed graphs, absorption time can be exponential in $n$; explicit constructions realizing such behavior exist.

A crucial variant is the *colonization regime*, where only mutants reproduce (mutant fitness $1$, resident fitness $0$, i.e., $r \to \infty$). Here, **colonization time** is always polynomial in $n$ for any underlying graph:
- **General bound:** $T(G_n) \le \frac12 n^3 - \frac12 n^2$ [2410.09476].
- **Tightness:** This is attained only by the *backward graph* $B_n$. For complete graphs $K_n$, colonization is $O(n\log n)$, the minimum possible order.
- For undirected/two-way graphs, $O(n^{2.5})$ applies; for regular/lattice-like graphs, $O(n^2)$ is tight.

These results contrast sharply with the classic case ($r > 1$, residents reproduce), where some graphs yield exponential fixation or absorption times [1311.7631][2410.09476][1804.02293].


## 4. Impact of Spatial Structure and Heterogeneity

Graph structure, spatial heterogeneity in connectivity, and environmental effects are central determinants of evolutionary outcome:
- **Environmental heterogeneity**: Variation in node-specific fitnesses for mutants and/or residents has asymmetrical effects. In large populations, variance in mutant fitness suppresses fixation probability, while resident heterogeneity has no effect (but can amplify fixation in small systems). The selection condition for mutants is $\overline{a} > \overline{b}$, where $\overline{a}$ and $\overline{b}$ are arithmetic means over environments [1709.03031].
- **Resource heterogeneity**: Represented as node colorings (e.g., two-colorings for resource-rich/poor sites), determines fixation probabilities. All properly two-colored, undirected, regular graphs are evolutionarily equivalent; fixation probability depends only on the coloring parameters, not the details of graph connectivity. Dynamically fluctuating resources scramble these effects, tending toward neutrality [2009.10693].
- **Random graphs**: On $G_{n,p}$ at the connectivity threshold, the degree of the initial mutant is the dominant predictor of fixation probability. Extreme degree variability can drive fixation outcomes away from classical predictions [2409.11615].
- **Spatial and genetic constraints**: In structured populations, mating restrictions and spatial grids increase effective mutation rates. Explicit thresholds for speciation by pattern formation—as a balance of mutation, spatial, and genetic constraints—can be derived [1012.3913].


## 5. Computational and Algorithmic Aspects

Exact computation of fixation probabilities is generally infeasible for arbitrary graphs due to exponential state space. However, **fully polynomial randomized approximation schemes (FPRAS)** for both fixation ($r \ge 1$) and extinction probabilities ($r > 0$) exist for undirected graphs, enabled by polynomial absorption time bounds [1111.3321][1311.7631][1804.02293]:
- The key Monte Carlo method is to simulate the process many times, estimating probabilities by frequency of absorption outcomes.
- Early stopping rules and improved lower bounds (depending on minimum degree, drift) sharply reduce simulation requirements.
- In the generalized and heterogeneous Moran processes (e.g., node- or type-specific fitness), submodularity in mutant-biased settings enables greedy approximation algorithms with $(1-1/e)$-guarantees for seed set selection [2404.15986][2201.02248].
- For processes with positional advantage, optimization is NP-hard but tractable in special cases (strong/weak selection), with efficient schemes for undirected graphs and weak selection limit [2201.02248].
- With type-specific graphs (mutants/residents have different connectivity), FPRAS exists only in restricted scenarios (e.g., mutant clique) [1710.07365].

Analytic approaches for evolutionary game theory and complex dynamics include path-integral formulations, yielding exact transition probabilities as sums over evolutionary trajectories, thus capturing the full stochastic dynamics beyond mean-field or diffusive approximations [2209.01060].


## 6. Extensions: Mutation, Genealogy, Frequency-Dependence, and Replacer Phenotype

The basic Moran process has been extended in multiple directions:
- **Mutation:** Allowing for mutation leads to stationary distributions or recurrent Markov processes. Metastable switching times and invariant measures can be analyzed via diffusion/WKB approximations; finite populations exhibit skew and rare transitions absent in deterministic limits [1512.06886].
- **Multivariate/frequency-dependent selection:** The multitype Moran process accommodates arbitrary $S$ types with frequency-dependent selection, yielding direct correspondences (under appropriate scaling) to the deterministic Lotka-Volterra equations; stability and coexistence in the deterministic model translate to features of the Moran process's quasi-stationary distribution [1110.4444].
- **Sampled Genealogy Process:** Incorporates explicit genealogical trees under Moran demography, capturing the full distribution of possible sample genealogies, stationary under Kingman's coalescent, and extending to arbitrary finite $n$ and arbitrary sampling patterns [2002.11184].
- **Replacer mutants:** A novel "neighborhood-aware" phenotype, which targets only the opposite type for replacement, achieves substantially higher fixation probability for neutral and even deleterious mutants (scaling as $1/\sqrt{N}$ rather than $1/N$ in well-mixed populations), and dramatically enhances evolutionary stability once established [2511.04417].


## 7. Measuring and Interpreting Evolutionary Timescales

Time in the Moran process is not unique: **step counting** (number of discrete updates) can yield paradoxes in spatial populations (e.g., starting with more mutants can paradoxically increase time to fixation in certain lollipop graphs). The **real duration** should reflect the cumulative reproductive rate; in the colonization regime, each update's "weight" is $1/k$ with $k$ mutants present, aligning with a continuous-time Markov process. This measure reveals the true scaling of evolutionary dynamics across different spatial structures and avoids artifacts of classic step-counting [2410.09476].

The key quantitative relationships are:
- Colonization time (steps): $T(G_n)$
- Real evolutionary duration: $R(G_n,v) = T(G_n,v)/n$
- Bounds (general/undirected/regular): $O(n^3), O(n^{2.5}), O(n^2)$ in steps, translating to $n^2, n^{1.5}, n$ in real time units.

In summary, the Moran process and its myriad extensions reveal that the structure of populations, the heterogeneity and dynamism of environments, and the details of reproductive strategy and updating all exert quantifiable, sometimes extreme, influence on evolutionary fate. Upper and lower bounds—attained or tightly approximated by explicit extremal graphs—shape both the simulation feasibility and conceptual understanding of selection, drift, and invasion in structured populations.

Source: https://www.emergentmind.com/topics/moran-process