Wright-Fisher Intra-Patch Reproduction

• Wright-Fisher Intra-Patch Reproduction is a model capturing stochastic allele dynamics in fixed-size populations through genetic drift and natural selection.
• It employs Fisher’s angular transformation to convert coordinate-dependent diffusion into tractable closed-form Green’s functions for both neutral and selective scenarios.
• The approach yields robust short-time approximations validated by simulations, providing precise estimates of fixation and extinction probabilities.

Wright-Fisher intra-patch reproduction describes the stochastic evolutionary dynamics of competing alleles or variants within a fixed-size population (“patch”), integrating the influences of genetic drift and natural selection. The formalism is anchored in the Wright-Fisher Fokker-Planck equation, enabling precise quantification of transition probabilities for allele frequencies over short timescales. Fisher’s angular transformation provides a natural way to regularize these dynamics, leading to closed-form Green’s functions that unify both neutral and selective scenarios and have direct applications in population genetics, evolutionary inference, and deep sequencing datasets (Khatri, 2015).

1. The Wright-Fisher Fokker-Planck Equation in Fixed-Size Populations

The intra-patch regime assumes a fixed population size $N$, within which two or more self-reproducing variants compete. Let $G(x,x_{0};t)$ denote the transition density of the allele frequency $x \in [0,1]$ at time $t$ given initial frequency $x_0$. In the diffusion approximation, the Wright-Fisher Fokker-Planck equation is

$$\frac{\partial}{\partial t}\,G(x,x_{0};t) = -\frac{\partial}{\partial x}\!\Bigl[s\,x(1-x)\,G\Bigr] + \frac{1}{2N}\,\frac{\partial^{2}}{\partial x^{2}}\!\Bigl[x(1-x)\,G\Bigr],$$

with $G(x,x_{0};0)=\delta(x-x_{0})$. The “drift” term ($s$) encodes directional selection, while the “diffusion” term models genetic drift with variance $x(1-x)/N$. The equation is subject to absorbing boundaries at $x=0$ (extinction) and $G(x,x_{0};t)$0 (fixation).

2. Fisher’s Angular Transformation and the Natural Length Scale

To simplify the coordinate-dependent diffusion, Fisher introduced the angular transformation: $G(x,x_{0};t)$1 In $G(x,x_{0};t)$2-space, the Jacobian $G(x,x_{0};t)$3 regularizes the diffusion, mapping the frequency-dependent $G(x,x_{0};t)$4 to a constant $G(x,x_{0};t)$5 in $G(x,x_{0};t)$6. This recasting allows direct analogy with Brownian motion, where the dynamics become coordinate-independent and more mathematically tractable.

3. Stochastic Differential Equation and the Effective Potential Landscape

In the angular coordinate, under neutrality ($G(x,x_{0};t)$7), the transformation yields the SDE: $G(x,x_{0};t)$8 This drift can be expressed via the potential $G(x,x_{0};t)$9. The potential features a maximum at $x \in [0,1]$0 with divergences at boundaries ($x \in [0,1]$1 or $x \in [0,1]$2), corresponding to extinction or fixation. The allele frequency behaves as a Brownian particle in an unstable potential, highlighting the drift-driven tendency for rare variants to be lost and common variants to fix.

4. Short-Time Asymptotic Green’s Functions: Neutral and Selective Cases

Neutral Case ($x \in [0,1]$3)

Expanding the drift around $x \in [0,1]$4 linearizes the process, leading to: $x \in [0,1]$5 The Green’s function in $x \in [0,1]$6 and then $x \in [0,1]$7 space is

$x \in [0,1]$8

Arbitrary Selection ($x \in [0,1]$9)

Here, the SDE in $t$0 space is

$t$1

Using a heuristic Gaussian approximation:

• The mean $t$2 solves the nonlinear deterministic ODE,

$t$3

with $t$4, $t$5.

• The variance is determined by linearizing the drift around $t$6,

$t$7

where $t$8.

The resulting $t$9 retains a closed Gaussian form.

5. Validity, Accuracy, and Regime Limitations

The heuristic Gaussian approximation assumes the process in $x_0$0-space remains approximately Gaussian for short times ($x_0$1 mean fixation time), with the mean given by the full nonlinear ODE and variance by the local linearization. This regime applies for $x_0$2 (neutral drift time scale) or $x_0$3 (selected), and for frequencies not near absorbing boundaries. The approximation yields excellent agreement with stochastic simulation, accurately capturing fixation/extinction tails in these windows (Khatri, 2015). For longer times or frequencies near boundaries, accuracy necessarily degrades.

6. Implications and Applications in Evolutionary Genetics

Within a patch, allele frequencies evolve as Brownian particles in an unstable potential, providing an explicit physical interpretation for the eventual fate (fixation or extinction) of variants. The closed-form, short-time Green’s functions allow explicit computation of transition densities and time-dependent probabilities of fixation/extinction, unifying neutral and selective models in a single analytic framework. These results are particularly applicable to inference of selection from allele-frequency time series in high-throughput sequencing, especially in experimental microbial or viral populations with short generation times and large $x_0$4, directly informing models of evolutionary dynamics under strong genetic drift and selection (Khatri, 2015).

7. Connections and Extensions

The Fisher angular transformation not only regularizes the mathematical treatment of the Wright-Fisher process but also clarifies connections to other potential well models in stochastic processes. The asymptotic Green’s function framework developed by Hallatschek and colleagues establishes a rigorous foundation for time-dependent population-genetic inference. A plausible implication is that similar coordinate transformations could prove fruitful in other classes of stochastic, finite-population processes, particularly where strong boundary effects and short-time dynamics are central.

For researchers analyzing deep sequencing or experimental evolution datasets, these intra-patch Wright-Fisher results provide the foundation for robust estimation of selection coefficients and evolutionary forecasting in finite populations (Khatri, 2015).