Stabilizing Interactions in Assemblages with Weak Shared Evolutionary History

Do multispecies communities need deep coevolutionary history to generate stabilizing interactions that maintain species coexistence? We investigate through computational analysis of Lotka-Volterra communities across a coevolutionary gradient.

Category: q-bio.PE — Populations and Evolution | Reference: Reynebeau et al. (2026) arXiv:2601.14213

NFD at θ=0

+0.0728
Mean invasion growth rate in novel assemblages with no shared evolutionary history

Species Invading from Rarity

65.7%
Fraction of species with positive invasion growth rate at θ=0

Stabilizing Niche Diff.

0.525
Mean (1-ρ) at θ=0 -- substantial without coevolution

Eco-Evo α Reduction

-8.0%
Mean interspecific competition declines from 0.618 to 0.569 over 500 generations

Methods

Three complementary computational analyses addressing the open problem

Experiment 1 Coevolutionary Gradient Sweep

Sweep θ from 0 (random assembly) to 1 (coevolved community) across Lotka-Volterra communities with S=10 species, 30 replicates per θ value, measuring NFD via invasion-from-rarity analysis.

Experiment 2 Eco-Evolutionary Rescue

Starting from θ=0, allow interaction coefficients to evolve via mutation and selection over 500 generations, testing whether NFD can emerge de novo through niche differentiation.

Experiment 3 MCT Pairwise Decomposition

Decompose coexistence into stabilizing niche differences (1-ρ) and fitness ratios under Modern Coexistence Theory across the coevolutionary gradient (S=8, 50 replicates).

Coevolutionary Gradient Sweep

How does shared evolutionary history (θ) affect NFD strength and species persistence?

NFD Strength vs. Coevolutionary History

Species Persistence vs. Coevolutionary History

Fraction of Species with Positive Invasion Growth Rate

NFD Strength: θ=0 vs θ=1 (Summary)

Eco-Evolutionary Rescue Dynamics

Can novel assemblages develop NFD through rapid evolution?

NFD Emergence over Eco-Evolutionary Generations

Mean Interspecific Competition (αij) over Time

Modern Coexistence Theory Decomposition

Stabilizing niche differences and pairwise coexistence across θ

Stabilizing Niche Difference (1-ρ) and Niche Overlap (ρ)

Fraction of Coexisting Species Pairs

Summary Results

Key metrics at selected θ values from the gradient sweep (S=10, 30 replicates)

θ Mean NFD SD NFD Surviving (of 10) % Positive Invasion 1-ρ (MCT) Coexisting Pairs
0.00+0.07280.03186.6765.7%0.52597.0%
0.10+0.08960.01797.5074.3%0.51897.2%
0.25+0.07420.02797.3773.7%0.52598.3%
0.50+0.06280.01787.2771.3%0.52898.2%
0.75+0.04260.01116.6765.0%0.53593.9%
1.00+0.00270.00234.0040.0%0.54585.6%

Eco-Evolutionary Rescue Summary

Metrics at start and end of 500-generation simulation (S=10, θinitial=0)

Metric Generation 0 Generation 499 Change
Mean NFD+0.0859+0.0866+0.8%
Mean αij0.6180.569-8.0%
Surviving species910+1
Fraction positive invasion0.900.80-0.10

Key Conclusions

Finding 1 NFD Does Not Require Coevolution

Novel assemblages (θ=0) exhibit positive NFD with mean invasion growth rate +0.0728 and 65.7% of species achieving positive invasion rates. Stabilizing niche differences (1-ρ=0.525) are nearly as large as in coevolved communities.

Finding 2 Coevolution Reshapes but Does Not Uniformly Strengthen NFD

Increasing θ reduces mean NFD and fraction of coexisting pairs because structured niche partitioning introduces fitness asymmetries between neighboring species on the niche axis.

Finding 3 Eco-Evolutionary Dynamics Maintain NFD

Novel assemblages maintain positive NFD over 500 generations while mean interspecific competition decreases by 8.0% through evolved niche differentiation, suggesting rapid coevolution can stabilize novel communities.