Can deleterious mutations surf deterministic population waves?

Technical Summary: Can Deleterious Mutations Surf Deterministic Population Waves?

Problem Statement
The paper addresses the phenomenon of "gene surfing," where rare mutations spread across large spatial regions during a range expansion. While gene surfing of neutral mutations is well-established, the behavior of deleterious mutations remains contentious. Experimental and simulation studies offer conflicting evidence: some suggest deleterious mutations fail to surf, others suggest transient surfing, and some indicate they can spread extensively, leading to "expansion load" (a progressive decline in fitness at the front). The central question is whether deleterious mutations can surf deterministic population waves, or if their presence at the expansion front is merely a result of mutation-selection equilibrium.

Methodology
The authors analyze a deterministic version of the spatial Muller's ratchet, modeling an asexual population subject to mutation, migration, and density-dependent reproduction and death. The system is described by an infinite set of coupled reaction-diffusion partial differential equations (PDEs), where $u_k(t, x)$ represents the population density of individuals carrying exactly $k$ deleterious mutations.

Key methodological components include:

1. Mathematical Framework: The study establishes the well-posedness (existence and uniqueness) of mild solutions to the infinite system of PDEs in the space of summable sequences ($\ell^1$).
2. Feynman–Kac Representation: To overcome the challenge of non-local interactions in type space (where the reproduction/death rates depend on the total population density $\|u\|_{\ell^1}$), the authors derive a Feynman–Kac representation for each coordinate $u_k$. This probabilistic representation allows for quantitative control over the evolution of the mutational profile.
3. Tracer Dynamics: To distinguish between surfing and equilibrium, the authors employ a "tracer dynamics" approach. They label a subset of the initial population carrying deleterious mutations and track the density of these labeled descendants ($u^*$) over time. If surfing occurs, the labeled mutants should propagate with the wave; if not, their density should decay relative to the wild type.
4. Regime Analysis: The analysis focuses on the "monostable" regime, specifically under Fisher–KPP conditions (where competition dominates cooperation), and derives bounds on the ratios of subpopulations with different mutation loads.

Key Contributions and Results

• Well-Posedness: The paper rigorously proves the existence and uniqueness of global non-negative mild solutions to the infinite system of PDEs (Theorem 2.1). It establishes that solutions are regular ($C^{1,2}$) and uniformly bounded.
• Quantitative Bounds on Mutational Profiles: In the monostable regime, the authors derive quantitative bounds on the ratio between the density of individuals carrying $k$ mutations ($u_k$) and the density of mutation-free individuals ($u_0$). They show that as time $T \to \infty$, this ratio is uniformly bounded between two sequences that converge to specific limits determined by the mutation rate and fitness parameters (Theorem 2.4). This confirms that the mutational profile remains well-behaved during expansion.
• Spreading Speed: Under Fisher–KPP conditions, the authors rigorously determine the spreading speed $c^*$ of the population into an empty habitat. They confirm the non-rigorous computations of Foutel-Rodier and Etheridge [15], showing that $c^* = \sqrt{2m((1-\mu)q_+(0) - q_-(0))}$ (Theorem 2.6). Furthermore, they demonstrate that all subpopulations (regardless of mutation load) spread at this same speed.
• No Surfing of Deleterious Mutations: The primary result (Theorem 2.7) utilizes tracer dynamics to show that deleterious mutations cannot surf deterministic waves. While the fraction of the population carrying deleterious mutations remains positive at the expansion front (due to the mutation-selection equilibrium established in Theorem 2.4), the labeled descendants of the initial mutant subpopulation decay exponentially. The expansion front is colonized by new mutations arising from the wild-type (mutation-free) subpopulation, not by the surfing of pre-existing deleterious lineages.

Significance and Claims
The paper claims to provide a rigorous mathematical justification for the conjecture that deleterious mutations do not surf in the deterministic limit of the spatial Muller's ratchet.

• Resolution of Contradictions: The authors suggest their results help explain contradictory findings in biological data regarding expansion load. They posit that while surfing may occur in finite, stochastic populations (due to demographic noise), the deterministic limit (representing very large populations) precludes it. Thus, any surfing effect in finite populations must be limited.
• Validation of Previous Work: The results rigorously validate the non-rigorous derivations and conjectures of Foutel-Rodier and Etheridge [15] regarding the spreading speed and the equilibrium distribution of mutations.
• Methodological Advancement: The work demonstrates the utility of Feynman–Kac representations in analyzing non-local, infinite-dimensional reaction-diffusion systems where standard comparison principles fail.

The paper concludes that in the deterministic monostable regime, the presence of deleterious mutations at the expansion front is a consequence of the continuous generation of new mutations from the wild type (mutation-selection equilibrium), rather than the successful surfing of existing deleterious lineages.