The Red Queen hypothesis, named after Lewis Carroll's character who must keep running just to stay in place, proposes that hosts and parasites are locked in a perpetual evolutionary arms race. As hosts evolve new resistance alleles, parasites evolve to overcome them, causing allele frequencies in both populations to oscillate indefinitely rather than reach a stable equilibrium. This matching-alleles model is implemented using coupled ordinary differential equations: dp/dt = p(1−p)(w_H − w̄_H) and dq/dt = q(1−q)(w_P − w̄_P), where host fitness decreases when the parasite allele frequency matches the host allele, and vice versa.
The simulation lets you adjust selection strength (the cost of being matched), mutation rate, and starting allele frequencies for both host and parasite. The time-series panel shows the cycling allele frequencies, while the phase-plane panel reveals the trajectory as a closed orbit — a hallmark of neutral oscillations driven by frequency-dependent selection.
What is the Red Queen hypothesis?
The Red Queen hypothesis (Van Valen, 1973) states that species must constantly evolve simply to maintain their current fitness relative to co-evolving antagonists such as parasites. In host–parasite systems this manifests as perpetual cycling of resistance and virulence alleles, with neither side ever gaining a permanent advantage.
Why do allele frequencies oscillate rather than reach equilibrium?
The dynamics are driven by frequency-dependent selection: a rare host genotype has high fitness because parasites are not adapted to it, but as the host genotype becomes common, parasites evolve to target it, reducing its fitness again. This negative frequency dependence creates a restoring force that keeps frequencies cycling rather than fixing at any value.
What does the matching-alleles model assume?
In the matching-alleles model, a parasite is infectious only if its allele matches the host allele exactly — like a lock and key. Host fitness is therefore 1 − s·q (where q is the parasite allele frequency) and parasite fitness is 1 + s·p. This is contrasted with the gene-for-gene model used for plant–pathogen coevolution, where each host resistance gene corresponds to a specific avirulence gene.
The selection coefficient s sets how strongly being matched by a parasite reduces host fitness (and boosts parasite fitness). Higher s produces faster and wider oscillations, while very low s produces slow, shallow cycles. When s approaches 1, allele frequencies swing dramatically between near-fixation at 0 and 1.
Mutation continuously injects both alleles at rate μ, preventing either allele from reaching fixation at 0 or 1. Without mutation the orbits can spiral inward or outward depending on the model. With mutation the system tends toward stable cycles around the interior equilibrium point at p = q = 0.5.
The phase plane plots parasite frequency q against host frequency p. In the matching-alleles model the trajectory forms a closed elliptical orbit around the equilibrium (0.5, 0.5). The colour gradient from green to red shows how the trajectory evolved through time, and the yellow dot marks the current state. Circular orbits indicate neutral oscillations; spiralling inward or outward indicates damping or amplification.
Strong evidence comes from the New Zealand snail Potamopyrgus antipodarum and its trematode parasites: sexual snail populations show much higher allele diversity than asexual clones, consistent with Red Queen coevolution. Similar cycling has been measured in Daphnia water fleas and their parasitic bacteria Pasteuria ramosa over decades of lake sediment records.
The simulation counts how many times the host frequency p crosses the midpoint (0.5) from below and above and uses these zero-crossings to estimate the period. Because the model is deterministic, the cycle length is approximately 2π / √(s²) generations — longer cycles arise from weaker selection or higher mutation that damps the oscillations.
This is one of the most debated questions in evolutionary biology. The Red Queen hypothesis proposes that sex is maintained because it shuffles allele combinations each generation, making offspring unpredictable to parasites adapted to the parent's genotype. Evidence from Potamopyrgus snails supports this: populations under intense parasite pressure maintain higher rates of sexual reproduction than parasite-free populations.
Starting at the interior equilibrium (p₀ = q₀ = 0.5) places the system exactly at the fixed point, so in the absence of noise or mutation it stays there. Any small perturbation — or the mutation term in the model — will push it off this unstable equilibrium and into oscillation. This demonstrates that the equilibrium is neutrally stable in the pure matching-alleles model.