Towards a mathematical understanding of invasion resistance in multispecies communities

Multispecies community composition and dynamics are key to health and disease across biological systems, a prominent example being microbial ecosystems. Explaining the forces that govern diversity and resilience in the microbial consortia making up our body’s defences remains a challenge. In this, theoretical models are crucial, to bridge the gap between species dynamics and underlying mechanisms and to develop analytic insight. Here we propose a replicator equation framework to model multispecies dynamics where an explicit notion of invasion resistance of a system emerges and can be studied explicitly. For illustration, we derive the conceptual link between such replicator equation and N microbial species’ growth and interaction traits, stemming from micro-scale environmental modification. Within this replicator framework, mean invasion fitness arises, evolves dynamically, and may undergo critical predictable shifts with global environmental changes. This mathematical approach clarifies the key role of this resident system trait for invader success, and highlights interaction principles among N species that optimize their collective resistance to invasion. We propose this model based on the replicator equation as a powerful new avenue to study, test and validate mechanisms of invasion resistance and colonization in multispecies microbial ecosystems and beyond.