Coexistence of many species under a random competition-colonization trade-off.

The competition-colonization (CC) trade-off is a well-studied coexistence mechanism for metacommunities. In this setting, it is believed that the coexistence of all species requires their traits to satisfy restrictive conditions limiting their similarity. To investigate whether diverse metacommunities can assemble in a CC trade-off model, we study their assembly from a probabilistic perspective. From a pool of species with parameters (corresponding to traits) sampled at random, we compute the probability that any number of species coexist and characterize the set of species that emerges through assembly. Remarkably, almost exactly half of the species in a large pool typically coexist, with no saturation as the size of the pool grows, and with little dependence on the underlying distribution of traits. Through a mix of analytical results and simulations, we show that this unlimited niche packing emerges as assembly actively moves communities toward overdispersed configurations in niche space. Our findings also apply to a realistic assembly scenario where species invade one at a time from a fixed regional pool. When diversity arises de novo in the metacommunity, richness still grows without bound, but more slowly. Together, our results suggest that the CC trade-off can support the robust emergence of diverse communities, even when coexistence of the full species pool is exceedingly unlikely.

Metapopulations with habitat modification.

Across the tree of life, organisms modify their local environment, rendering it more or less hospitable for other species. Despite the ubiquity of these processes, simple models that can be used to develop intuitions about the consequences of widespread habitat modification are lacking. Here, we extend the classic Levins metapopulation model to a setting where each of n species can colonize patches connected by dispersal, and when patches are vacated via local extinction, they retain a "memory" of the previous occupant-modeling habitat modification. While this model can exhibit a wide range of dynamics, we draw several overarching conclusions about the effects of modification and memory. In particular, we find that any number of species may potentially coexist, provided that each is at a disadvantage when colonizing patches vacated by a conspecific. This notion is made precise through a quantitative stability condition, which provides a way to unify and formalize existing conceptual models. We also show that when patch memory facilitates coexistence, it generically induces a positive relationship between diversity and robustness (tolerance of disturbance). Our simple model provides a portable, tractable framework for studying systems where species modify and react to a shared landscape.

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.

Habitat Heterogeneity, Environmental Feedbacks, and Species Coexistence across Timescales.

AbstractClassic ecological theory explains species coexistence in variable environments. While spatial variation is often treated as an intrinsic feature of a landscape, it may be shaped and even generated by the resident community. All species modify their local environment to some extent, driving changes that can feed back to affect the composition and coexistence of the community, potentially over timescales very different from population dynamics. We introduce a simple nested modeling framework for community dynamics in heterogeneous environments, including the possible evolution of heterogeneity over time due to community-environment feedbacks. We use this model to derive analytical conditions for species coexistence in environments where heterogeneity is either fixed or shaped by feedbacks. Among other results, our approach reveals how dispersal and environmental specialization interact to shape realized patterns of habitat association and demonstrates that environmental feedbacks can tune landscape conditions to allow the stable coexistence of any number of species. Our flexible modeling framework helps explain feedback dynamics that arise in a wide range of ecosystems and offers a generic platform for exploring the interplay between species and landscape diversity.

Higher-order interactions stabilize dynamics in competitive network models.

Ecologists have long sought a way to explain how the remarkable biodiversity observed in nature is maintained. On the one hand, simple models of interacting competitors cannot produce the stable persistence of very large ecological communities. On the other hand, neutral models, in which species do not interact and diversity is maintained by immigration and speciation, yield unrealistically small fluctuations in population abundance, and a strong positive correlation between a species' abundance and its age, contrary to empirical evidence. Models allowing for the robust persistence of large communities of interacting competitors are lacking. Here we show that very diverse communities could persist thanks to the stabilizing role of higher-order interactions, in which the presence of a species influences the interaction between other species. Although higher-order interactions have been studied for decades, their role in shaping ecological communities is still unclear. The inclusion of higher-order interactions in competitive network models stabilizes dynamics, making species coexistence robust to the perturbation of both population abundance and parameter values. We show that higher-order interactions have strong effects in models of closed ecological communities, as well as of open communities in which new species are constantly introduced. In our framework, higher-order interactions are completely defined by pairwise interactions, facilitating empirical parameterization and validation of our models.

Beyond pairwise mechanisms of species coexistence in complex communities.

The tremendous diversity of species in ecological communities has motivated a century of research into the mechanisms that maintain biodiversity. However, much of this work examines the coexistence of just pairs of competitors. This approach ignores those mechanisms of coexistence that emerge only in diverse competitive networks. Despite the potential for these mechanisms to create conditions under which the loss of one competitor triggers the loss of others, we lack the knowledge needed to judge their importance for coexistence in nature. Progress requires borrowing insight from the study of multitrophic interaction networks, and coupling empirical data to models of competition.

No robust multispecies coexistence in a canonical model of plant-soil feedbacks.

Plant-soil feedbacks (PSFs) are considered a key mechanism generating frequency-dependent dynamics in plant communities. Negative feedbacks, in particular, are often invoked to explain coexistence and the maintenance of diversity in species-rich communities. However, the primary modelling framework used to study PSFs considers only two plant species, and we lack clear theoretical expectations for how these complex interactions play out in communities with natural levels of diversity. Here, we extend this canonical model of PSFs to include an arbitrary number of plant species and analyse the dynamics. Surprisingly, we find that coexistence of more than two species is virtually impossible, suggesting that alternative theoretical frameworks are needed to describe feedbacks observed in diverse natural communities. Drawing on our analysis, we discuss future directions for PSF models and implications for experimental study of PSF-mediated coexistence in the field.

Tractable models of ecological assembly.

Ecological assembly is a fundamental and yet poorly understood process. Three main obstacles hinder the development of a theory of assembly, and when these issues are sidestepped by making strong assumptions, one can build an assembly graph in which nodes are ecological communities and edges are invasions shifting their composition. The graph can then be analysed directly, without the need to consider dynamics. To showcase this framework, we build and analyse assembly graphs for the competitive Lotka-Volterra model, showing that in these cases sequential assembly (in which species invade a community one at a time) can reach the same configurations found when starting the system with all species present at different initial conditions. We discuss how our results can advance our understanding of assembly both from an empirical and a theoretical point of view, informing the study of ecological restoration and the design of ecological communities.

Telling ecological networks apart by their structure: A computational challenge.

Ecologists have been compiling ecological networks for over a century, detailing the interactions between species in a variety of ecosystems. To this end, they have built networks for mutualistic (e.g., pollination, seed dispersal) as well as antagonistic (e.g., herbivory, parasitism) interactions. The type of interaction being represented is believed to be reflected in the structure of the network, which would differ substantially between mutualistic and antagonistic networks. Here, we put this notion to the test by attempting to determine the type of interaction represented in a network based solely on its structure. We find that, although it is easy to separate different kinds of nonecological networks, ecological networks display much structural variation, making it difficult to distinguish between mutualistic and antagonistic interactions. We therefore frame the problem as a challenge for the community of scientists interested in computational biology and machine learning. We discuss the features a good solution to this problem should possess and the obstacles that need to be overcome to achieve this goal.

Last name analysis of mobility, gender imbalance, and nepotism across academic systems.

In biology, last names have been used as proxy for genetic relatedness in pioneering studies of neutral theory and human migrations. More recently, analyzing the last name distribution of Italian academics has raised the suspicion of nepotism, with faculty hiring their relatives for academic posts. Here, we analyze three large datasets containing the last names of all academics in Italy, researchers from France, and those working at top public institutions in the United States. Through simple randomizations, we show that the US academic system is geographically well-mixed, whereas Italian academics tend to work in their native region. By contrasting maiden and married names, we can detect academic couples in France. Finally, we detect the signature of nepotism in the Italian system, with a declining trend. The claim that our tests detect nepotism as opposed to other effects is supported by the fact that we obtain different results for the researchers hired after 2010, when an antinepotism law was in effect.

Transcriptome resilience predicts thermotolerance in Caenorhabditis elegans.

BACKGROUND: The detrimental effects of a short bout of stress can persist and potentially turn lethal, long after the return to normal conditions. Thermotolerance, which is the capacity of an organism to withstand relatively extreme temperatures, is influenced by the response during stress exposure, as well as the recovery process afterwards. While heat-shock response mechanisms have been studied intensively, predicting thermal tolerance remains a challenge. RESULTS: Here, we use the nematode Caenorhabditis elegans to measure transcriptional resilience to heat stress and predict thermotolerance. Using principal component analysis in combination with genome-wide gene expression profiles collected in three high-resolution time series during control, heat stress, and recovery conditions, we infer a quantitative scale capturing the extent of stress-induced transcriptome dynamics in a single value. This scale provides a basis for evaluating transcriptome resilience, defined here as the ability to depart from stress-expression dynamics during recovery. Independent replication across multiple highly divergent genotypes reveals that the transcriptional resilience parameter measured after a spike in temperature is quantitatively linked to long-term survival after heat stress. CONCLUSION: Our findings imply that thermotolerance is an intrinsic property that pre-determines long-term outcome of stress and can be predicted by the transcriptional resilience parameter. Inferring the transcriptional resilience parameters of higher organisms could aid in evaluating rehabilitation strategies after stresses such as disease and trauma.

Phenotypic variability promotes diversity and stability in competitive communities.

Intraspecific variation is at the core of evolutionary theory, and yet, from an ecological perspective, we have few robust expectations for how this variation should affect the dynamics of large communities. Here, by adapting an approach from evolutionary game theory, we show that the incorporation of phenotypic variability into competitive networks dramatically alters the dynamics across ecological timescales, stabilising the systems and buffering the communities against demographic perturbations. The beneficial effects of phenotypic variability are strongest when there are substantial differences among phenotypes and when phenotypes are inherited with moderately high fidelity; yet even low levels of variation lead to significant increases in diversity, stability, and robustness. By identifying a simple and ubiquitous stabilising force in competitive communities, this work contributes to our core understanding of how biological diversity is maintained in natural systems.

Reconciling empirical interactions and species coexistence.

Coexistence in ecological communities is governed largely by the nature and intensity of species interactions. Countless studies have proposed methods to infer these interactions from empirical data, yet models parameterised using such data often fail to recover observed coexistence patterns. Here, we propose a method to reconcile empirical parameterisations of community dynamics with species-abundance data, ensuring that the predicted equilibrium is consistent with the observed abundance distribution. To illustrate the approach, we explore two case studies: an experimental freshwater algal community and a long-term time series of displacement in an intertidal community. We demonstrate how our method helps recover observed coexistence patterns, capture the core dynamics of the system, and, in the latter case, predict the impacts of experimental extinctions. Collectively, these results demonstrate an intuitive approach for reconciling observed and empirical data, improving our ability to explore the links between species interactions and coexistence in natural systems.

Network spandrels reflect ecological assembly.

Ecological networks that exhibit stable dynamics should theoretically persist longer than those that fluctuate wildly. Thus, network structures which are over-represented in natural systems are often hypothesised to be either a cause or consequence of ecological stability. Rarely considered, however, is that these network structures can also be by-products of the processes that determine how new species attempt to join the community. Using a simulation approach in tandem with key results from random matrix theory, we illustrate how historical assembly mechanisms alter the structure of ecological networks. We demonstrate that different community assembly scenarios can lead to the emergence of structures that are often interpreted as evidence of 'selection for stability'. However, by controlling for the underlying selection pressures, we show that these assembly artefacts-or spandrels-are completely unrelated to stability or selection, and are instead by-products of how new species are introduced into the system. We propose that these network-assembly spandrels are critically overlooked aspects of network theory and stability analysis, and we illustrate how a failure to adequately account for historical assembly can lead to incorrect inference about the causes and consequences of ecological stability.

Understanding the role of parasites in food webs using the group model.

Parasites are ubiquitous and have been shown to influence macroscopic measures of ecological network structure, such as connectance and robustness, as well as local structure, such as subgraph frequencies. Nevertheless, they are often under-represented in ecological studies due to their small size and often complex life cycles. We consider whether or not parasites play structurally unique roles in ecological networks; that is, can we distinguish parasites from other species using network structure alone? We partition the species in a community statistically using the group model, and we test whether or not parasites tend to cluster in their own groups, using a measure of "imbalance." We find that parasites form highly imbalanced groups, and that concomitant predation, in which a predator consumes a prey and its parasites, but not the number of interactions, improves the group model's ability to distinguish parasites from non-parasites. This work demonstrates that parasites and non-parasites interact in networks in statistically distinct ways, and that these differences are partly, but not entirely, due to the existence of concomitant predation.

Stability criteria for complex ecosystems.

Forty years ago, May proved that sufficiently large or complex ecological networks have a probability of persisting that is close to zero, contrary to previous expectations. May analysed large networks in which species interact at random. However, in natural systems pairs of species have well-defined interactions (for example predator-prey, mutualistic or competitive). Here we extend May's results to these relationships and find remarkable differences between predator-prey interactions, which are stabilizing, and mutualistic and competitive interactions, which are destabilizing. We provide analytic stability criteria for all cases. We use the criteria to prove that, counterintuitively, the probability of stability for predator-prey networks decreases when a realistic food web structure is imposed or if there is a large preponderance of weak interactions. Similarly, stability is negatively affected by nestedness in bipartite mutualistic networks. These results are found by separating the contribution of network structure and interaction strengths to stability. Stable predator-prey networks can be arbitrarily large and complex, provided that predator-prey pairs are tightly coupled. The stability criteria are widely applicable, because they hold for any system of differential equations.

The Effect of Intra- and Interspecific Competition on Coexistence in Multispecies Communities.

For two competing species, intraspecific competition must exceed interspecific competition for coexistence. To generalize this well-known criterion to multiple competing species, one must take into account both the distribution of interaction strengths and community structure. Here we derive a multispecies generalization of the two-species rule in the context of symmetric Lotka-Volterra competition and obtain explicit stability conditions for random competitive communities. We then explore the influence of community structure on coexistence. Results show that both the most and least stabilized cases have striking global structures, with a nested pattern emerging in both cases. The distribution of intraspecific coefficients leading to the most and least stabilized communities also follows a predictable pattern that can be justified analytically. In addition, we show that the size of the parameter space allowing for feasible communities always increases with the strength of intraspecific effects in a characteristic way that is independent of the interspecific interaction structure. We conclude by discussing possible extensions of our results to nonsymmetric competition.

Ocean acidification affects competition for space: projections of community structure using cellular automata.

Historical ecological datasets from a coastal marine community of crustose coralline algae (CCA) enabled the documentation of ecological changes in this community over 30 years in the Northeast Pacific. Data on competitive interactions obtained from field surveys showed concordance between the 1980s and 2013, yet also revealed a reduction in how strongly species interact. Here, we extend these empirical findings with a cellular automaton model to forecast ecological dynamics. Our model suggests the emergence of a new dominant competitor in a global change scenario, with a reduced role of herbivory pressure, or trophic control, in regulating competition among CCA. Ocean acidification, due to its energetic demands, may now instead play this role in mediating competitive interactions and thereby promote species diversity within this guild.

Metapopulation persistence in random fragmented landscapes.

Habitat destruction and land use change are making the world in which natural populations live increasingly fragmented, often leading to local extinctions. Although local populations might undergo extinction, a metapopulation may still be viable as long as patches of suitable habitat are connected by dispersal, so that empty patches can be recolonized. Thus far, metapopulations models have either taken a mean-field approach, or have modeled empirically-based, realistic landscapes. Here we show that an intermediate level of complexity between these two extremes is to consider random landscapes, in which the patches of suitable habitat are randomly arranged in an area (or volume). Using methods borrowed from the mathematics of Random Geometric Graphs and Euclidean Random Matrices, we derive a simple, analytic criterion for the persistence of the metapopulation in random fragmented landscapes. Our results show how the density of patches, the variability in their value, the shape of the dispersal kernel, and the dimensionality of the landscape all contribute to determining the fate of the metapopulation. Using this framework, we derive sufficient conditions for the population to be spatially localized, such that spatially confined clusters of patches act as a source of dispersal for the whole landscape. Finally, we show that a regular arrangement of the patches is always detrimental for persistence, compared to the random arrangement of the patches. Given the strong parallel between metapopulation models and contact processes, our results are also applicable to models of disease spread on spatial networks.

What Can Interaction Webs Tell Us About Species Roles?

The group model is a useful tool to understand broad-scale patterns of interaction in a network, but it has previously been limited in use to food webs, which contain only predator-prey interactions. Natural populations interact with each other in a variety of ways and, although most published ecological networks only include information about a single interaction type (e.g., feeding, pollination), ecologists are beginning to consider networks which combine multiple interaction types. Here we extend the group model to signed directed networks such as ecological interaction webs. As a specific application of this method, we examine the effects of including or excluding specific interaction types on our understanding of species roles in ecological networks. We consider all three currently available interaction webs, two of which are extended plant-mutualist networks with herbivores and parasitoids added, and one of which is an extended intertidal food web with interactions of all possible sign structures (+/+, -/0, etc.). Species in the extended food web grouped similarly with all interactions, only trophic links, and only nontrophic links. However, removing mutualism or herbivory had a much larger effect in the extended plant-pollinator webs. Species removal even affected groups that were not directly connected to those that were removed, as we found by excluding a small number of parasitoids. These results suggest that including additional species in the network provides far more information than additional interactions for this aspect of network structure. Our methods provide a useful framework for simplifying networks to their essential structure, allowing us to identify generalities in network structure and better understand the roles species play in their communities.