Robust permanence for ecological equations with internal and external feedbacks

Building modern coexistence theory from the ground up: The role of community assembly

Modern coexistence theory (MCT) is one of the leading methods to understand species coexistence. It uses invasion growth rates-the average, per-capita growth rate of a rare species-to identify when and why species coexist. Despite significant advances in dissecting coexistence mechanisms when coexistence occurs, MCT relies on a 'mutual invasibility' condition designed for two-species communities but poorly defined for species-rich communities. Here, we review well-known issues with this component of MCT and propose a solution based on recent mathematical advances. We propose a clear framework for expanding MCT to species-rich communities and for understanding invasion resistance as well as coexistence, especially for communities that could not be analysed with MCT so far. Using two data-driven community models from the literature, we illustrate the utility of our framework and highlight the opportunities for bridging the fields of community assembly and species coexistence.

Permanence via invasion graphs: incorporating community assembly into modern coexistence theory

To understand the mechanisms underlying species coexistence, ecologists often study invasion growth rates of theoretical and data-driven models. These growth rates correspond to average per-capita growth rates of one species with respect to an ergodic measure supporting other species. In the ecological literature, coexistence often is equated with the invasion growth rates being positive. Intuitively, positive invasion growth rates ensure that species recover from being rare. To provide a mathematically rigorous framework for this approach, we prove theorems that answer two questions: (i) When do the signs of the invasion growth rates determine coexistence? (ii) When signs are sufficient, which invasion growth rates need to be positive? We focus on deterministic models and equate coexistence with permanence, i.e., a global attractor bounded away from extinction. For models satisfying certain technical assumptions, we introduce invasion graphs where vertices correspond to proper subsets of species (communities) supporting an ergodic measure and directed edges correspond to potential transitions between communities due to invasions by missing species. These directed edges are determined by the signs of invasion growth rates. When the invasion graph is acyclic (i.e. there is no sequence of invasions starting and ending at the same community), we show that permanence is determined by the signs of the invasion growth rates. In this case, permanence is characterized by the invasibility of all [Formula: see text] communities, i.e., communities without species i where all other missing species have negative invasion growth rates. To illustrate the applicability of the results, we show that dissipative Lotka-Volterra models generically satisfy our technical assumptions and computing their invasion graphs reduces to solving systems of linear equations. We also apply our results to models of competing species with pulsed resources or sharing a predator that exhibits switching behavior. Open problems for both deterministic and stochastic models are discussed. Our results highlight the importance of using concepts about community assembly to study coexistence.

Complex community-wide consequences of consumer sexual dimorphism

Sexual dimorphism is a ubiquitous source of within‐species variation, yet the community‐level consequences of sex differences remain poorly understood.

Here, we analyse a bitrophic model of two competing resource species and a sexually reproducing consumer species.

We show that consumer sex differences in resource acquisition can have striking consequences for consumer‐resource coexistence, abundance and dynamics.

Under both direct interspecific competition and apparent competition between two resource species, sexual dimorphism in consumers' attack rates can mediate coexistence of the resource species, while in other cases can lead to exclusion when stable coexistence is typically expected. Slight sex differences in total resource acquisition also can reverse competitive outcomes and lead to density cycles. These effects are expected whenever both consumer sexes require different amounts or types of resources to reproduce.

Our results suggest that consumer sexual dimorphism, which is common, has wide‐reaching implications for the assembly and dynamics of natural communities.

Sick of eating: Eco-evo-immuno dynamics of predators and their trophically acquired parasites

When predators consume prey, they risk becoming infected with their prey's parasites, which can then establish the predator as a secondary host. A predator population's diet therefore influences what parasites it is exposed to, as has been repeatedly shown in many species such as threespine stickleback (Gasterosteus aculeatus) (more benthic‐feeding individuals obtain nematodes from oligocheate prey, whereas limnetic‐feeding individuals catch cestodes from copepod prey). These differing parasite encounters, in turn, determine how natural selection acts on the predator's immune system. We might therefore expect that ecoevolutionary dynamics of a predator's diet (as determined by its ecomorphology) should drive correlated evolution of its immune traits. Conversely, the predator's immunity to certain parasites might alter the relative costs and benefits of different prey, driving evolution of its ecomorphology. To evaluate the potential for ecological morphology to drive evolution of immunity, and vice versa, we use a quantitative genetics framework coupled with an ecological model of a predator and two prey species (the diet options). Our analysis reveals fundamental asymmetries in the evolution of ecomorphology and immunity. When ecomorphology rapidly evolves, it determines how immunity evolves, but not vice versa. Weak trade‐offs in ecological morphology select for diet generalists despite strong immunological trade‐offs, but not vice versa. Only weak immunological trade‐offs can explain negative diet‐infection correlations across populations. The analysis also reveals that eco‐evo‐immuno feedbacks destabilize population dynamics when trade‐offs are sufficiently weak and heritability is sufficiently high. Collectively, these results highlight the delicate interplay between multivariate trait evolution and the dynamics of ecological communities.

Destabilizing evolutionary and eco-evolutionary feedbacks drive empirical eco-evolutionary cycles

We develop a method to identify how ecological, evolutionary, and eco-evolutionary feedbacks influence system stability. We apply our method to nine empirically parametrized eco-evolutionary models of exploiter-victim systems from the literature and identify which particular feedbacks cause some systems to converge to a steady state or to exhibit sustained oscillations. We find that ecological feedbacks involving the interactions between all species and evolutionary and eco-evolutionary feedbacks involving only the interactions between exploiter species (predators or pathogens) are typically stabilizing. In contrast, evolutionary and eco-evolutionary feedbacks involving the interactions between victim species (prey or hosts) are destabilizing more often than not. We also find that while eco-evolutionary feedbacks rarely altered system stability from what would be predicted from just ecological and evolutionary feedbacks, eco-evolutionary feedbacks have the potential to alter system stability at faster or slower speeds of evolution. As the number of empirical studies demonstrating eco-evolutionary feedbacks increases, we can continue to apply these methods to determine whether the patterns we observe are common in other empirical communities.

Persistence and extinction for stochastic ecological models with internal and external variables

The dynamics of species' densities depend both on internal and external variables. Internal variables include frequencies of individuals exhibiting different phenotypes or living in different spatial locations. External variables include abiotic factors or non-focal species. These internal or external variables may fluctuate due to stochastic fluctuations in environmental conditions. The interplay between these variables and species densities can determine whether a particular population persists or goes extinct. To understand this interplay, we prove theorems for stochastic persistence and exclusion for stochastic ecological difference equations accounting for internal and external variables. Specifically, we use a stochastic analog of average Lyapunov functions to develop sufficient and necessary conditions for (i) all population densities spending little time at low densities i.e. stochastic persistence, and (ii) population trajectories asymptotically approaching the extinction set with positive probability. For (i) and (ii), respectively, we provide quantitative estimates on the fraction of time that the system is near the extinction set, and the probability of asymptotic extinction as a function of the initial state of the system. Furthermore, in the case of persistence, we provide lower bounds for the expected time to escape neighborhoods of the extinction set. To illustrate the applicability of our results, we analyze stochastic models of evolutionary games, Lotka-Volterra dynamics, trait evolution, and spatially structured disease dynamics. Our analysis of these models demonstrates environmental stochasticity facilitates coexistence of strategies in the hawk-dove game, but inhibits coexistence in the rock-paper-scissors game and a Lotka-Volterra predator-prey model. Furthermore, environmental fluctuations with positive auto-correlations can promote persistence of evolving populations and persistence of diseases in patchy landscapes. While our results help close the gap between the persistence theories for deterministic and stochastic systems, we highlight several challenges for future research.

Stochastic population growth in spatially heterogeneous environments: the density-dependent case

This work is devoted to studying the dynamics of a structured population that is subject to the combined effects of environmental stochasticity, competition for resources, spatio-temporal heterogeneity and dispersal. The population is spread throughout n patches whose population abundances are modeled as the solutions of a system of nonlinear stochastic differential equations living on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[0,\infty )^n$$\end{document}[0,∞)n. We prove that r, the stochastic growth rate of the total population in the absence of competition, determines the long-term behaviour of the population. The parameter r can be expressed as the Lyapunov exponent of an associated linearized system of stochastic differential equations. Detailed analysis shows that if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ r>0$$\end{document}r>0, the population abundances converge polynomially fast to a unique invariant probability measure on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(0,\infty )^n$$\end{document}(0,∞)n, while when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ r<0$$\end{document}r<0, the population abundances of the patches converge almost surely to 0 exponentially fast. This generalizes and extends the results of Evans et al. (J Math Biol 66(3):423–476, 2013) and proves one of their conjectures. Compared to recent developments, our model incorporates very general density-dependent growth rates and competition terms. Furthermore, we prove that persistence is robust to small, possibly density dependent, perturbations of the growth rates, dispersal matrix and covariance matrix of the environmental noise. We also show that the stochastic growth rate depends continuously on the coefficients. Our work allows the environmental noise driving our system to be degenerate. This is relevant from a biological point of view since, for example, the environments of the different patches can be perfectly correlated. We show how one can adapt the nondegenerate results to the degenerate setting. As an example we fully analyze the two-patch case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2$$\end{document}n=2, and show that the stochastic growth rate is a decreasing function of the dispersion rate. In particular, coupling two sink patches can never yield persistence, in contrast to the results from the non-degenerate setting treated by Evans et al. which show that sometimes coupling by dispersal can make the system persistent.