Survival of the scarcer

Does deterministic coexistence theory matter in a finite world?

Contemporary studies of species coexistence are underpinned by deterministic models that assume that competing species have continuous (i.e., noninteger) densities, live in infinitely large landscapes, and coexist over infinite time horizons. By contrast, in nature, species are composed of discrete individuals subject to demographic stochasticity and occur in habitats of finite size where extinctions occur in finite time. One consequence of these discrepancies is that metrics of species' coexistence derived from deterministic theory may be unreliable predictors of the duration of species coexistence in nature. These coexistence metrics include invasion growth rates and niche and fitness differences, which are now commonly applied in theoretical and empirical studies of species coexistence. In this study, we tested the efficacy of deterministic coexistence metrics on the duration of species coexistence in a finite world. We introduce new theoretical and computational methods to estimate coexistence times in stochastic counterparts of classic deterministic models of competition. Importantly, we parameterized this model using experimental field data for 90 pairwise combinations of 18 species of annual plants, allowing us to derive biologically informed estimates of coexistence times for a natural system. Strikingly, we found that for species expected to deterministically coexist, community sizes containing only 10 individuals had predicted coexistence times of more than 1000 years. We also found that invasion growth rates explained 60% of the variation in intrinsic coexistence times, reinforcing their general usefulness in studies of coexistence. However, only by integrating information on both invasion growth rates and species' equilibrium population sizes could most (>99%) of the variation in species coexistence times be explained. This integration was achieved with demographically uncoupled single-species models solely determined by the invasion growth rates and equilibrium population sizes. Moreover, because of a complex relationship between niche overlap/fitness differences and equilibrium population sizes, increasing niche overlap and increasing fitness differences did not always result in decreasing coexistence times, as deterministic theory would predict. Nevertheless, our results tend to support the informed use of deterministic theory for understanding the duration of species' coexistence while highlighting the need to incorporate information on species' equilibrium population sizes in addition to invasion growth rates.

Stochastic resonance of abundance fluctuations and mean time to extinction in an ecological community

Periodic environmental changes are commonly observed in nature from the amount of daylight to seasonal temperature. These changes usually affect individuals' death or birth rates, dragging the system from its previous stable states. When the fluctuation of abundance is amplified due to such changes, extinction of species may be accelerated. To see this effect, we examine how the abundance and the mean time to extinction respond to the periodic environmental changes. We consider a population wherein two species coexist together implemented by three rules-birth, spontaneous death, and death from competitions. As the interspecific interaction strength is varied, we observe the resonance behavior in both fluctuations of abundances and the mean time to extinction. Our result suggests that neither too high nor too low competition rates make the system more susceptible to environmental changes.

Effects of niche overlap on coexistence, fixation and invasion in a population of two interacting species

Synergistic and antagonistic interactions in multi-species populations-such as resource sharing and competition-result in remarkably diverse behaviours in populations of interacting cells, such as in soil or human microbiomes, or clonal competition in cancer. The degree of inter- and intra-specific interaction can often be quantified through the notion of an ecological 'niche'. Typically, weakly interacting species that occupy largely distinct niches result in stable mixed populations, while strong interactions and competition for the same niche result in rapid extinctions of some species and fixations of others. We investigate the transition of a deterministically stable mixed population to a stochasticity-induced fixation as a function of the niche overlap between the two species. We also investigate the effect of the niche overlap on the population stability with respect to external invasions. Our results have important implications for a number of experimental systems.

CHAOTIC COEXISTENCE IN A RESOURCE–CONSUMER MODEL

We study the interspecies competition in a simple resource–consumer model which includes the resource supply as a dynamic variable. In the model, an organism of each species needs to consume a certain minimum rate of resource (food) to survive; a higher rate of consumption is required for reproduction. We analyze the orbit diagrams and Lyapunov exponents in order to characterize the system dynamical behavior. We observe that under a periodic food supply, the system can exhibit coexistence in the form of limit cycle oscillations. For a certain parameter range, we observe chaotic behavior emerging from successive period doublings and quasi-periodicity.

Population size changes and extinction risk of populations driven by mutant interactors

Spontaneous random mutations are an important source of variation in populations. Many evolutionary models consider mutants with a fixed fitness, chosen from a fitness distribution without considering microscopic interactions among the residents and mutants. Here, we go beyond this and consider "mutant interactors," which lead to new interactions between the residents and invading mutants that can affect the population size and the extinction risk of populations. We model microscopic interactions between individuals by using a dynamic interaction matrix, the dimension of which increases with the emergence of a new mutant and decreases with extinction. The new interaction parameters of the mutant follow a probability distribution around the payoff entries of its ancestor. These new interactions can drive the population away from the previous equilibrium and lead to changes in the population size. Thus, the population size is an evolving property rather than an externally controlled variable. We calculate the average population size of our stochastic system over time and quantify the extinction risk of the population by the mean time to extinction.

Fixation probabilities in populations under demographic fluctuations

We study the fixation probability of a mutant type when introduced into a resident population. We implement a stochastic competitive Lotka-Volterra model with two types and intra- and interspecific competition. The model further allows for stochastically varying population sizes. The competition coefficients are interpreted in terms of inverse payoffs emerging from an evolutionary game. Since our study focuses on the impact of the competition values, we assume the same net growth rate for both types. In this general framework, we derive a formula for the fixation probability [Formula: see text] of the mutant type under weak selection. We find that the most important parameter deciding over the invasion success of the mutant is its death rate due to competition with the resident. Furthermore, we compare our approximation to results obtained by implementing population size changes deterministically in order to explore the parameter regime of validity of our method. Finally, we put our formula in the context of classical evolutionary game theory and observe similarities and differences to the results obtained in that constant population size setting.

Extinction dynamics from metastable coexistences in an evolutionary game

Deterministic evolutionary game dynamics can lead to stable coexistences of different types. Stochasticity, however, drives the loss of such coexistences. This extinction is usually accompanied by population size fluctuations. We investigate the most probable extinction trajectory under such fluctuations by mapping a stochastic evolutionary model to a problem of classical mechanics using the Wentzel-Kramers-Brillouin (WKB) approximation. Our results show that more abundant types in a coexistence may be more likely to go extinct first, in good agreement with previous results. The distance between the coexistence and extinction points is not a good predictor of extinction either. Instead, the WKB method correctly predicts the type going extinct first.

The advantage of being slow: The quasi-neutral contact process

According to the competitive exclusion principle, in a finite ecosystem, extinction occurs naturally when two or more species compete for the same resources. An important question that arises is: when coexistence is not possible, which mechanisms confer an advantage to a given species against the other(s)? In general, it is expected that the species with the higher reproductive/death ratio will win the competition, but other mechanisms, such as asymmetry in interspecific competition or unequal diffusion rates, have been found to change this scenario dramatically. In this work, we examine competitive advantage in the context of quasi-neutral population models, including stochastic models with spatial structure as well as macroscopic (mean-field) descriptions. We employ a two-species contact process in which the “biological clock” of one species is a factor of αslower than that of the other species. Our results provide new insights into how stochasticity and competition interact to determine extinction in finite spatial systems. We find that a species with a slower biological clock has an advantage if resources are limited, winning the competition against a species with a faster clock, in relatively small systems. Periodic or stochastic environmental variations also favor the slower species, even in much larger systems.

Extinction of oscillating populations

Established populations often exhibit oscillations in their sizes that, in the deterministic theory, correspond to a limit cycle in the space of population sizes. If a population is isolated, the intrinsic stochasticity of elemental processes can ultimately bring it to extinction. Here we study extinction of oscillating populations in a stochastic version of the Rosenzweig-MacArthur predator-prey model. To this end we develop a WKB (Wentzel, Kramers and Brillouin) approximation to the master equation, employing the characteristic population size as the large parameter. Similar WKB theories have been developed previously in the context of population extinction from an attracting multipopulation fixed point. We evaluate the extinction rates and find the most probable paths to extinction from the limit cycle by applying Floquet theory to the dynamics of an effective four-dimensional WKB Hamiltonian. We show that the entropic barriers to extinction change in a nonanalytic way as the system passes through the Hopf bifurcation. We also study the subleading pre-exponential factors of the WKB approximation.

Spatial complementarity and the coexistence of species

Coexistence of apparently similar species remains an enduring paradox in ecology. Spatial structure has been predicted to enable coexistence even when population-level models predict competitive exclusion if it causes each species to limit its own population more than that of its competitor. Nevertheless, existing hypotheses conflict with regard to whether clustering favours or precludes coexistence. The spatial segregation hypothesis predicts that in clustered populations the frequency of intra-specific interactions will be increased, causing each species to be self-limiting. Alternatively, individuals of the same species might compete over greater distances, known as heteromyopia, breaking down clusters and opening space for a second species to invade. In this study we create an individual-based model in homogeneous two-dimensional space for two putative sessile species differing only in their demographic rates and the range and strength of their competitive interactions. We fully characterise the parameter space within which coexistence occurs beyond population-level predictions, thereby revealing a region of coexistence generated by a previously-unrecognised process which we term the triadic mechanism. Here coexistence occurs due to the ability of a second generation of offspring of the rarer species to escape competition from their ancestors. We diagnose the conditions under which each of three spatial coexistence mechanisms operates and their characteristic spatial signatures. Deriving insights from a novel metric - ecological pressure - we demonstrate that coexistence is not solely determined by features of the numerically-dominant species. This results in a common framework for predicting, given any pair of species and knowledge of the relevant parameters, whether they will coexist, the mechanism by which they will do so, and the resultant spatial pattern of the community. Spatial coexistence arises from complementary combinations of traits in each species rather than solely through self-limitation.