Speciation due to hybrid necrosis in plant-pathogen models

We develop a model for speciation due to postzygotic incompatibility generated by autoimmune reactions. The model is based on predator-prey interactions between a host plants and their pathogens. Such interactions are often frequency-dependent, so that pathogen attack is focused on the most abundant plant phenotype, while rare plant types may escape pathogen attack. Thus, frequency dependence can generate disruptive selection, which can give rise to speciation if distant phenotypes become reproductively isolated. Based on recent experimental evidence from {\it Arabidopsis}, we assume that at the molecular level, incompatibility between strains is caused by epistatic interactions between two proteins in the plant immune system, the guard and the guardee. Within each plant strain, immune reactions occur when the guardee protein is modified by a pathogen effector, and the guard subsequently binds to the guardee, thus precipitating an immune response. However, when guard and guardee proteins come from phenotypically distant parents, a hybrid’s immune system can be triggered by erroneous interactions between these proteins even in the absence of pathogen attack, leading to severe autoimmune reactions in hybrids. Our model shows how phenotypic variation generated by frequency-dependent host-pathogen interactions can lead to postzygotic incompatibility between extremal types, and hence to speciation.

💡 Research Summary

The paper presents a theoretical framework that links frequency‑dependent host‑pathogen interactions with post‑zygotic reproductive isolation driven by autoimmune incompatibility, ultimately leading to speciation. The authors begin by noting that hybrid necrosis—severe auto‑immune reactions observed in inter‑population hybrids of Arabidopsis—represents a concrete example of post‑zygotic incompatibility. Recent molecular work has identified the “guard–guardee” system as the mechanistic basis: a guard protein monitors a guardee that is modified by pathogen effectors; when modification occurs, the guard binds the guardee and triggers immunity. However, when guard and guardee alleles derived from phenotypically distant parents are combined in a hybrid, erroneous binding can occur even in the absence of a pathogen, causing lethal auto‑immunity.

To explore how such molecular incompatibility can arise from ecological dynamics, the authors construct a predator‑prey (host‑pathogen) model in which pathogen attack is frequency‑dependent: the pathogen preferentially targets the most common plant phenotype, while rare phenotypes experience reduced pressure. Mathematically, plant density N(x,t) is a continuous function of a phenotypic trait x, and pathogen attack intensity a(x) is modeled as a saturating function of N(x). Plant growth follows logistic dynamics with a mortality term proportional to a(x)·N(x). Pathogen dynamics depend on the average host phenotype, creating a feedback loop that drives the plant trait distribution.

Simulation of the coupled equations reveals a two‑stage process. First, negative frequency‑dependence generates disruptive selection, splitting the plant population into two distinct phenotypic peaks. Each peak evolves a distinct set of guard and guardee alleles that are co‑adapted to the local pathogen pressure. Second, when individuals from the two peaks interbreed, hybrids inherit mismatched guard–guardee combinations. Because the binding affinity between guard and guardee is highly sensitive to phenotypic distance, the hybrid’s immune system can be erroneously activated without any pathogen, producing hybrid necrosis. This creates a strong post‑zygotic barrier that dramatically reduces hybrid fitness.

Parameter sweeps show that speciation is most likely when (i) pathogen attack is broad enough to maintain strong frequency‑dependent pressure, and (ii) the guard–guardee interaction exhibits steep dependence on phenotypic divergence. If mutation rates are too low, phenotypic divergence never reaches the threshold for incompatibility; if they are too high, the system becomes globally unstable, with many genotypes suffering auto‑immunity. The model therefore predicts a “Goldilocks” zone of genetic variation that maximizes the likelihood of speciation via hybrid necrosis.

The authors validate the model qualitatively against Arabidopsis data, where specific NLR (nucleotide‑binding leucine‑rich repeat) and RIN4 variants from different ecotypes cause severe auto‑immune phenotypes when combined. This empirical observation aligns with the model’s prediction that divergent guard–guardee alleles generate hybrid incompatibility.

In conclusion, the study integrates ecological frequency‑dependent selection with molecular immunogenetics to propose a plausible pathway from host‑pathogen coevolution to reproductive isolation and speciation. It extends classic ecological speciation theory by incorporating a concrete genetic mechanism, offering testable predictions for future work in plant systems and potentially in other host‑parasite interactions.