Structured community transitions explain the switching capacity of microbial systems.

Microbial systems appear to exhibit a relatively high switching capacity of moving back and forth among few dominant communities (taxon memberships). While this switching behavior has been mainly attributed to random environmental factors, it remains unclear the extent to which internal community dynamics affect the switching capacity of microbial systems. Here, we integrate ecological theory and empirical data to demonstrate that structured community transitions increase the dependency of future communities on the current taxon membership, enhancing the switching capacity of microbial systems. Following a structuralist approach, we propose that each community is feasible within a unique domain in environmental parameter space. Then, structured transitions between any two communities can happen with probability proportional to the size of their feasibility domains and inversely proportional to their distance in environmental parameter space-which can be treated as a special case of the gravity model. We detect two broad classes of systems with structured transitions: one class where switching capacity is high across a wide range of community sizes and another class where switching capacity is high only inside a narrow size range. We corroborate our theory using temporal data of gut and oral microbiota (belonging to class 1) as well as vaginal and ocean microbiota (belonging to class 2). These results reveal that the topology of feasibility domains in environmental parameter space is a relevant property to understand the changing behavior of microbial systems. This knowledge can be potentially used to understand the relevant community size at which internal dynamics can be operating in microbial systems.

Rapid monitoring of ecological persistence.

Effective conservation of ecological communities requires accurate and up-to-date information about whether species are persisting or declining to extinction. The persistence of an ecological community is supported by its underlying network of species interactions. While the persistence of the network supporting the whole community is the most relevant scale for conservation, in practice, only small subsets of these networks can be monitored. There is therefore an urgent need to establish links between the small snapshots of data conservationists can collect, and the "big picture" conclusions about ecosystem health demanded by policymakers, scientists, and societies. Here, we show that the persistence of small subnetworks (motifs) in isolation-that is, their persistence when considered separately from the larger network of which they are a part-is a reliable probabilistic indicator of the persistence of the network as a whole. Our methods show that it is easier to detect if an ecological community is not persistent than if it is persistent, allowing for rapid detection of extinction risk in endangered systems. Our results also justify the common practice of predicting ecological persistence from incomplete surveys by simulating the population dynamics of sampled subnetworks. Empirically, we show that our theoretical predictions are supported by data on invaded networks in restored and unrestored areas, even in the presence of environmental variability. Our work suggests that coordinated action to aggregate information from incomplete sampling can provide a means to rapidly assess the persistence of entire ecological networks and the expected success of restoration strategies.

Multispecies coexistence in fragmented landscapes.

Spatial dynamics have long been recognized as an important driver of biodiversity. However, our understanding of species' coexistence under realistic landscape configurations has been limited by lack of adequate analytical tools. To fill this gap, we develop a spatially explicit metacommunity model of multiple competing species and derive analytical criteria for their coexistence in fragmented heterogeneous landscapes. Specifically, we propose measures of niche and fitness differences for metacommunities, which clarify how spatial dynamics and habitat configuration interact with local competition to determine coexistence of species. We parameterize our model with a Bayesian approach using a 36-y time-series dataset of three Daphnia species in a rockpool metacommunity covering >500 patches. Our results illustrate the emergence of interspecific variation in extinction and recolonization processes, including their dependencies on habitat size and environmental temperature. We find that such interspecific variation contributes to the coexistence of Daphnia species by reducing fitness differences and increasing niche differences. Additionally, our parameterized model allows separating the effects of habitat destruction and temperature change on species extinction. By integrating coexistence theory and metacommunity theory, our study provides platforms to increase our understanding of species' coexistence in fragmented heterogeneous landscapes and the response of biodiversity to environmental changes.

On the limits to invasion prediction using coexistence outcomes.

The dynamics of ecological communities in nature are typically characterized by probabilistic processes involving invasion dynamics. Because of technical challenges, however, the majority of theoretical and experimental studies have focused on coexistence dynamics. Therefore, it has become central to understand the extent to which coexistence outcomes can be used to predict analogous invasion outcomes relevant to systems in nature. Here, we study the limits to this predictability under a geometric and probabilistic Lotka-Volterra framework. We show that while individual survival probability in coexistence dynamics can be fairly closely translated into invader colonization probability in invasion dynamics, the translation is less precise between community persistence and community augmentation, and worse between exclusion probability and replacement probability. These results provide a guiding and testable theoretical framework regarding the translatability of outcomes between coexistence and invasion outcomes when communities are represented by Lotka-Volterra dynamics under environmental uncertainty.

Protection promotes energetically efficient structures in marine communities.

The sustainability of marine communities is critical for supporting many biophysical processes that provide ecosystem services that promote human well-being. It is expected that anthropogenic disturbances such as climate change and human activities will tend to create less energetically-efficient ecosystems that support less biomass per unit energy flow. It is debated, however, whether this expected development should translate into bottom-heavy (with small basal species being the most abundant) or top-heavy communities (where more biomass is supported at higher trophic levels with species having larger body sizes). Here, we combine ecological theory and empirical data to demonstrate that full marine protection promotes shifts towards top-heavy energetically-efficient structures in marine communities. First, we use metabolic scaling theory to show that protected communities are expected to display stronger top-heavy structures than disturbed communities. Similarly, we show theoretically that communities with high energy transfer efficiency display stronger top-heavy structures than communities with low transfer efficiency. Next, we use empirical structures observed within fully protected marine areas compared to disturbed areas that vary in stress from thermal events and adjacent human activity. Using a nonparametric causal-inference analysis, we find a strong, positive, causal effect between full marine protection and stronger top-heavy structures. Our work corroborates ecological theory on community development and provides a quantitative framework to study the potential restorative effects of different candidate strategies on protected areas.

Experimental evidence of the importance of multitrophic structure for species persistence.

Ecological theory predicts that species interactions embedded in multitrophic networks shape the opportunities for species to persist. However, the lack of experimental support of this prediction has limited our understanding of how species interactions occurring within and across trophic levels simultaneously regulate the maintenance of biodiversity. Here, we integrate a mathematical approach and detailed experiments in plant-pollinator communities to demonstrate the need to jointly account for species interactions within and across trophic levels when estimating the ability of species to persist. Within the plant trophic level, we show that the persistence probability of plant species increases when introducing the effects of plant-pollinator interactions. Across trophic levels, we show that the persistence probabilities of both plants and pollinators exhibit idiosyncratic changes when experimentally manipulating the multitrophic structure. Importantly, these idiosyncratic effects are not recovered by traditional simulations. Our work provides tractable experimental and theoretical platforms upon which it is possible to investigate the multitrophic factors affecting species persistence in ecological communities.

Ranking species based on sensitivity to perturbations under non-equilibrium community dynamics.

Managing ecological communities requires fast detection of species that are sensitive to perturbations. Yet, the focus on recovery to equilibrium has prevented us from assessing species responses to perturbations when abundances fluctuate over time. Here, we introduce two data-driven approaches (expected sensitivity and eigenvector rankings) based on the time-varying Jacobian matrix to rank species over time according to their sensitivity to perturbations on abundances. Using several population dynamics models, we demonstrate that we can infer these rankings from time-series data to predict the order of species sensitivities. We find that the most sensitive species are not always the ones with the most rapidly changing or lowest abundance, which are typical criteria used to monitor populations. Finally, using two empirical time series, we show that sensitive species tend to be harder to forecast. Our results suggest that incorporating information on species interactions can improve how we manage communities out of equilibrium.

Understanding the impact of third-party species on pairwise coexistence.

The persistence of virtually every single species depends on both the presence of other species and the specific environmental conditions in a given location. Because in natural settings many of these conditions are unknown, research has been centered on finding the fraction of possible conditions (probability) leading to species coexistence. The focus has been on the persistence probability of an entire multispecies community (formed of either two or more species). However, the methodological and philosophical question has always been whether we can observe the entire community and, if not, what the conditions are under which an observed subset of the community can persist as part of a larger multispecies system. Here, we derive long-term (using analytical calculations) and short-term (using simulations and experimental data) system-level indicators of the effect of third-party species on the coexistence probability of a pair (or subset) of species under unknown environmental conditions. We demonstrate that the fraction of conditions incompatible with the possible coexistence of a pair of species tends to become vanishingly small within systems of increasing numbers of species. Yet, the probability of pairwise coexistence in isolation remains approximately the expected probability of pairwise coexistence in more diverse assemblages. In addition, we found that when third-party species tend to reduce (resp. increase) the coexistence probability of a pair, they tend to exhibit slower (resp. faster) rates of competitive exclusion. Long-term and short-term effects of the remaining third-party species on all possible specific pairs in a system are not equally distributed, but these differences can be mapped and anticipated under environmental uncertainty.

Untangling the complexity of priority effects in multispecies communities.

The history of species immigration can dictate how species interact in local communities, thereby causing historical contingency in community assembly. Since immigration history is rarely known, these historical influences, or priority effects, pose a major challenge in predicting community assembly. Here, we provide a graph-based, non-parametric, theoretical framework for understanding the predictability of community assembly as affected by priority effects. To develop this framework, we first show that the diversity of possible priority effects increases super-exponentially with the number of species. We then point out that, despite this diversity, the consequences of priority effects for multispecies communities can be classified into four basic types, each of which reduces community predictability: alternative stable states, alternative transient paths, compositional cycles and the lack of escapes from compositional cycles to stable states. Using a neural network, we show that this classification of priority effects enables accurate explanation of community predictability, particularly when each species immigrates repeatedly. We also demonstrate the empirical utility of our theoretical framework by applying it to two experimentally derived assembly graphs of algal and ciliate communities. Based on these analyses, we discuss how the framework proposed here can help guide experimental investigation of the predictability of history-dependent community assembly.

Understanding the emergence of contingent and deterministic exclusion in multispecies communities.

Competitive exclusion can be classified as deterministic or as historically contingent. While competitive exclusion is common in nature, it has remained unclear when multispecies communities formed by more than two species should be dominated by deterministic or contingent exclusion. Here, we take a fully parameterised model of an empirical competitive system between invasive annual and native perennial plant species to explain both the emergence and sources of competitive exclusion in multispecies communities. Using a structural approach to understand the range of parameters promoting deterministic and contingent exclusions, we then find heuristic theoretical support for the following three general conclusions. First, we find that the life-history of perennial species increases the probability of observing contingent exclusion by increasing their effective intrinsic growth rates. Second, we find that the probability of observing contingent exclusion increases with weaker intraspecific competition, and not with the level of hierarchical competition. Third, we find a shift from contingent exclusion to deterministic exclusion with increasing numbers of competing species. Our work provides a heuristic framework to increase our understanding about the predictability of species persistence within multispecies communities.

Observed Ecological Communities Are Formed by Species Combinations That Are among the Most Likely to Persist under Changing Environments.

AbstractDespite the rich biodiversity found in nature, it is unclear to what extent some combinations of interacting species, while conceivable in a given place and time, may never be realized. Yet solving this problem is important for understanding the role of randomness and predictability in the assembly of ecological communities. Here we show that the specific combinations of interacting species that emerge from the ecological dynamics within regional species pools are not all equally likely to be seen; rather, they are among the most likely to persist under changing environments. First, we use niche-based competition matrices and Lotka-Volterra models to demonstrate that realized combinations of interacting species are more likely to persist under random parameter perturbations than the majority of potential combinations with the same number of species that could have been formed from the regional pool. We then corroborate our theoretical results using a 10-year observational study, recording 88 plant-herbivore communities across three different forest successional stages. By inferring and validating plant-mediated communities of competing herbivore species, we find that observed combinations of herbivores have an expected probability of species persistence higher than half of all potential combinations. Our findings open up the opportunity to establish a formal probabilistic and predictive understanding of the composition of ecological communities.

Telling ecological networks apart by their structure: An environment-dependent approach.

The network architecture of an ecological community describes the structure of species interactions established in a given place and time. It has been suggested that this architecture presents unique features for each type of ecological interaction: e.g., nested and modular architectures would correspond to mutualistic and antagonistic interactions, respectively. Recently, Michalska-Smith and Allesina (2019) proposed a computational challenge to test whether it is indeed possible to differentiate ecological interactions based on network architecture. Contrary to the expectation, they found that this differentiation is practically impossible, moving the question to why it is not possible to differentiate ecological interactions based on their network architecture alone. Here, we show that this differentiation becomes possible by adding the local environmental information where the networks were sampled. We show that this can be explained by the fact that environmental conditions are a confounder of ecological interactions and network architecture. That is, the lack of association between network architecture and type of ecological interactions changes by conditioning on the local environmental conditions. Additionally, we find that environmental conditions are linked to the stability of ecological networks, but the direction of this effect depends on the type of interaction network. This suggests that the association between ecological interactions and network architectures exists, but cannot be fully understood without attention to the environmental conditions acting upon them.

Merging dynamical and structural indicators to measure resilience in multispecies systems.

Resilience is broadly understood as the ability of an ecological system to resist and recover from perturbations acting on species abundances and on the system's structure. However, one of the main problems in assessing resilience is to understand the extent to which measures of recovery and resistance provide complementary information about a system. While recovery from abundance perturbations has a strong tradition under the analysis of dynamical stability, it is unclear whether this same formalism can be used to measure resistance to structural perturbations (e.g. perturbations to model parameters). Here, we provide a framework grounded on dynamical and structural stability in Lotka-Volterra systems to link recovery from small perturbations on species abundances (i.e. dynamical indicators) with resistance to parameter perturbations of any magnitude (i.e. structural indicators). We use theoretical and experimental multispecies systems to show that the faster the recovery from abundance perturbations, the higher the resistance to parameter perturbations. We first use theoretical systems to show that the return rate along the slowest direction after a small random abundance perturbation (what we call full recovery) is negatively correlated with the largest random parameter perturbation that a system can withstand before losing any species (what we call full resistance). We also show that the return rate along the second fastest direction after a small random abundance perturbation (what we call partial recovery) is negatively correlated with the largest random parameter perturbation that a system can withstand before at most one species survives (what we call partial resistance). Then, we use a dataset of experimental microbial systems to confirm our theoretical expectations and to demonstrate that full and partial components of resilience are complementary. Our findings reveal that we can obtain the same level of information about resilience by measuring either a dynamical (i.e. recovery) or a structural (i.e. resistance) indicator. Irrespective of the chosen indicator (dynamical or structural), our results show that we can obtain additional information by separating the indicator into its full and partial components. We believe these results can motivate new theoretical approaches and empirical analyses to increase our understanding about risk in ecological systems.

Structural forecasting of species persistence under changing environments.

The persistence of a species in a given place not only depends on its intrinsic capacity to consume and transform resources into offspring, but also on how changing environmental conditions affect its growth rate. However, the complexity of factors has typically taken us to choose between understanding and predicting the persistence of species. To tackle this limitation, we propose a probabilistic approach rooted on the statistical concepts of ensemble theory applied to statistical mechanics and on the mathematical concepts of structural stability applied to population dynamics models - what we call structural forecasting. We show how this new approach allows us to estimate a probability of persistence for single species in local communities; to understand and interpret this probability conditional on the information we have concerning a system; and to provide out-of-sample predictions of species persistence as good as the best experimental approaches without the need of extensive amounts of data.

On the Consequences of the Interdependence of Stabilizing and Equalizing Mechanisms.

We present an overlooked but important property of modern coexistence theory (MCT), along with two key new results and their consequences. The overlooked property is that stabilizing mechanisms (increasing species' niche differences) and equalizing mechanisms (reducing species' fitness differences) have two distinct sets of meanings within MCT: one in a two-species context and another in a general multispecies context. We demonstrate that the two-species framework is not a special case of the multispecies one, and therefore these two parallel frameworks must be studied independently. Our first result is that, using the two-species framework and mechanistic consumer-resource models, stabilizing and equalizing mechanisms exhibit complex interdependence, such that changing one will simultaneously change the other. Furthermore, the nature and direction of this simultaneous change sensitively depend on model parameters. The second result states that while MCT is often seen as bridging niche and neutral modes of coexistence by building a niche-neutrality continuum, the interdependence between stabilizing and equalizing mechanisms acts to break this continuum under almost any biologically relevant circumstance. We conclude that the complex entanglement of stabilizing and equalizing terms makes their impact on coexistence difficult to understand, but by seeing them as aggregated effects (rather than underlying causes) of coexistence, we may increase our understanding of ecological dynamics.

Beware z-scores.

Song, Rohr, and Saavedra (2017) have proposed a methodology to compare network properties across systems with different sizes and constraints, in response to the fact that z-scores cannot be used for such purposes. Simmons, Hoeppke, and Sutherland (2019) have shown that part of the methodology can be improved. Here, we show that all previous results hold and are strengthened by the new methodology.

Structural changes within trophic levels are constrained by within-family assembly rules at lower trophic levels.

Historical contingency broadly refers to the proposition that even random historical events can constrain the ecological and evolutionary pathways of organisms and that of entire communities. Focusing on communities, these pathways can be reflected into specific structural changes within and across trophic levels - how species interact with and affect each other - which has important consequences for species coexistence. Using the registry of the last 2000 years of plant introductions and their novel herbivores encountered in Central Europe, we find that the order of arrival of closely related (but not of distantly related) plant species constrained the structural changes within the trophic level formed by herbivore species across the observation period. Because it is difficult for field and lab experiments to be conducted over hundreds of years to record and replay the assembly history of a community, our study provides an alternative to understand how structural changes have occurred across extensive periods of time.

Structural stability as a consistent predictor of phenological events.

The timing of the first and last seasonal appearance of a species in a community typically follows a pattern that is governed by temporal factors. While it has been shown that changes in the environment are linked to phenological changes, the direction of this link appears elusive and context-dependent. Thus, finding consistent predictors of phenological events is of central importance for a better assessment of expected changes in the temporal dynamics of ecological communities. Here we introduce a measure of structural stability derived from species interaction networks as an estimator of the expected range of environmental conditions compatible with the existence of a community. We test this measure as a predictor of changes in species richness recorded on a daily basis in a high-arctic plant-pollinator community during two spring seasons. We find that our measure of structural stability is the only consistent predictor of changes in species richness among different ecological and environmental variables. Our findings suggest that measures based on the notion of structural stability can synthesize the expected variation of environmental conditions tolerated by a community, and explain more consistently the phenological changes observed in ecological communities.

Will a small randomly assembled community be feasible and stable?

How likely is it that few species can randomly assemble into a feasible and stable community? Some studies have answered that as long as the community is feasible, it will nearly always be stable. In contrast, other studies have answered that the likelihood is almost null. Here, we show that the origin of this debate has been the underestimation of the association of the parameter space of intrinsic growth rates with the feasibility and stability properties of small randomly-assembled communities. In particular, we demonstrate that not all parameterizations and sampling distributions of intrinsic growth rates lead to the same probabilities of stability and feasibility, which could mistakenly lead to under- or overestimate the stability properties of feasible communities. Additionally, we find that stability imposes a filtering of species abundances "towards" more even distributions in small feasible randomly-assembled communities. This indicates that the stability of feasible communities is inherently linked to the starting distribution of species abundances, a characteristic that many times has been ignored, but should be incorporated in manageable lab and field experiments. Overall, the return to this debate is a central reminder that a more systematic exploration of the feasible parameter space is necessary to derive general conclusions about the stability properties of ecological communities.

A guideline to study the feasibility domain of multi-trophic and changing ecological communities.

The feasibility domain of an ecological community can be described by the set of environmental abiotic and biotic conditions under which all co-occurring and interacting species in a given site and time can have positive abundances. Mathematically, the feasibility domain corresponds to the parameter space compatible with positive (feasible) solutions at equilibrium for all the state variables in a system under a given model of population dynamics. Under specific dynamics, the existence of a feasible equilibrium is a necessary condition for species persistence regardless of whether the feasible equilibrium is dynamically stable or not. Thus, the size of the feasibility domain can also be used as an indicator of the tolerance of a community to random environmental variations. This has motivated a rich research agenda to estimate the feasibility domain of ecological communities. However, these methodologies typically assume that species interactions are static, or that input and output energy flows on each trophic level are unconstrained. Yet, this is different to how communities behave in nature. Here, we present a step-by-step quantitative guideline providing illustrative examples, computational code, and mathematical proofs to study systematically the feasibility domain of ecological communities under changes of interspecific interactions and subject to different constraints on the trophic energy flows. This guideline covers multi-trophic communities that can be formed by any type of interspecific interactions. Importantly, we show that the relative size of the feasibility domain can significantly change as a function of the biological information taken into consideration. We believe that the availability of these methods can allow us to increase our understanding about the limits at which ecological communities may no longer tolerate further environmental perturbations, and can facilitate a stronger integration of theoretical and empirical research.