A structural approach for understanding multispecies coexistence

The development of ecological systems along paths of least resistance

A long-standing question in biology is whether there are common principles that characterize the development of ecological systems (the appearance of a group of taxa), regardless of organismal diversity and environmental context.1,2,3,4,5,6,7,8,9,10,11 Classic ecological theory holds that these systems develop following a sequenced, orderly process that generally proceeds from fast-growing to slow-growing taxa and depends on life-history trade-offs.2,12,13 However, it is also possible that this developmental order is simply the path with the least environmental resistance for survival of the component species and hence favored by probability alone. Here, we use theory and data to show that the order from fast- to slow-growing taxa is the most likely developmental path for diverse systems when local taxon interactions self-organize in light of environmental resistance. First, we demonstrate theoretically that a sequenced development is more likely than a simultaneous one, at least until the number of iterations becomes so large as to be ecologically implausible. We then show that greater diversity of taxa and life histories improves the likelihood of a sequenced order from fast- to slow-growing taxa. Using data from bacterial and metazoan systems,14,15,16,17,18,19 we present empirical evidence that the developmental order of ecological systems moves along the paths of least environmental resistance. The capacity of simple principles to explain the trend in the developmental order of diverse ecological systems paves the way to an enhanced understanding of collective features of life.

Soil microbial influences over coexistence potential in multispecies plant communities in a subtropical forest

Soil microbes have long been recognized to substantially affect the coexistence of pairwise plant species across terrestrial ecosystems. However, projecting their impacts on the coexistence of multispecies plant systems remains a pressing challenge. To address this challenge, we conducted a greenhouse experiment with 540 seedlings of five tree species in a subtropical forest in China and evaluated microbial effects on multispecies coexistence using the structural method, which quantifies how the structure of species interactions influences the likelihood for multiple species to persist. Specifically, we grew seedlings alone or with competitors in different microbial contexts and fitted individual biomass to a population dynamic model to calculate intra- and interspecific interaction strength with and without soil microbes. We then used these interaction structures to calculate two metrics of multispecies coexistence, structural niche differences (which promote coexistence) and structural fitness differences (which drive exclusion), for all possible communities comprising two to five plant species. We found that soil microbes generally increased both the structural niche and fitness differences across all communities, with a much stronger effect on structural fitness differences. A further examination of functional traits between plant species pairs found that trait differences are stronger predictors of structural niche differences than of structural fitness differences, and that soil microbes have the potential to change trait-mediated plant interactions. Our findings underscore that soil microbes strongly influence the coexistence of multispecies plant systems, and also add to the experimental evidence that the influence is more on fitness differences rather than on niche differences.

Consequences of Local Conspecific Density Effects for Plant Diversity and Community Dynamics

Conspecific density dependence (CDD) in plant populations is widespread, most likely caused by local-scale biotic interactions, and has potentially important implications for biodiversity, community composition, and ecosystem processes. However, progress in this important area of ecology has been hindered by differing viewpoints on CDD across subfields in ecology, lack of synthesis across CDD-related frameworks, and misunderstandings about how empirical measurements of local CDD fit within the context of broader ecological theories on community assembly and diversity maintenance. Here, we propose a conceptual synthesis of local-scale CDD and its causes, including species-specific antagonistic and mutualistic interactions. First, we compare and clarify different uses of CDD and related concepts across subfields within ecology. We suggest the use of local stabilizing/destabilizing CDD to refer to the scenario where local conspecific density effects are more negative/positive than heterospecific effects. Second, we discuss different mechanisms for local stabilizing and destabilizing CDD, how those mechanisms are interrelated, and how they cut across several fields of study within ecology. Third, we place local stabilizing/destabilizing CDD within the context of broader ecological theories and discuss implications and challenges related to scaling up the effects of local CDD on populations, communities, and metacommunities. The ultimate goal of this synthesis is to provide a conceptual roadmap for researchers studying local CDD and its implications for population and community dynamics.

Biodiversity in changing environments: An external-driver internal-topology framework to guide intervention

Accompanying the climate crisis is the more enigmatic biodiversity crisis. Rapid reorganization of biodiversity due to global environmental change has defied prediction and tested the basic tenets of conservation and restoration. Conceptual and practical innovation is needed to support decision making in the face of these unprecedented shifts. Critical questions include: How can we generalize biodiversity change at the community level? When are systems able to reorganize and maintain integrity, and when does abiotic change result in collapse or restructuring? How does this understanding provide a template to guide when and how to intervene in conservation and restoration? To this end, we frame changes in community organization as the modulation of external abiotic drivers on the internal topology of species interactions, using plant-plant interactions in terrestrial communities as a starting point. We then explore how this framing can help translate available data on species abundance and trait distributions to corresponding decisions in management. Given the expectation that community response and reorganization are highly complex, the external-driver internal-topology (EDIT) framework offers a way to capture general patterns of biodiversity that can help guide resilience and adaptation in changing environments.

Towards mechanistic integration of the causes and consequences of biodiversity

The global biodiversity crisis has stimulated decades of research on three themes: species coexistence, biodiversity-ecosystem functioning relationships (BEF), and biodiversity-ecosystem functional stability relationships (BEFS). However, studies on these themes are largely independent, creating barriers to an integrative understanding of the causes and consequences of biodiversity. Here we review recent progress towards mechanistic integration of coexistence, BEF, and BEFS. Mechanisms underlying the three themes can be linked in various ways, potentially creating either positive or negative relationships between them. That said, we generally expect positive associations between coexistence and BEF, and between BEF and BEFS. Our synthesis represents an initial step towards integrating causes and consequences of biodiversity; future developments should include more mechanistic approaches and broader ecological contexts.

When can higher-order interactions produce stable coexistence?

Most ecological models are based on the assumption that species interact in pairs. Diverse communities, however, can have higher-order interactions, in which two or more species jointly impact the growth of a third species. A pitfall of the common pairwise approach is that it misses the higher-order interactions potentially responsible for maintaining natural diversity. Here, we explore the stability properties of systems where higher-order interactions guarantee that a specified set of abundances is a feasible equilibrium of the dynamics. Even these higher-order interactions which lead to equilibria do not necessarily produce stable coexistence. Instead, these systems are more likely to be stable when the pairwise interactions are weak or facilitative. Correlations between the pairwise and higher-order interactions, however, do permit robust coexistence even in diverse systems. Our work not only reveals the challenges in generating stable coexistence through higher-order interactions but also uncovers interaction patterns that can enable diversity.

Parallel evolution of alternate morphotypes of Chryseobacterium gleum during experimental evolution with Caenorhabditis elegans

Microbial evolution within polymicrobial communities is a complex process. Here, we report within-species diversification within multispecies microbial communities during experimental evolution with the nematode Caenorhabditis elegans. We describe morphological diversity in the target species Chryseobacterium gleum, which developed a novel colony morphotype in a small number of replicate communities. Alternate morphotypes coexisted with original morphotypes in communities, as well as in single-species experiments using evolved isolates. We found that the original and alternate morphotypes differed in motility and in spatial expansion in the presence of C. elegans. This study provides insight into the emergence and maintenance of intraspecies diversity in the context of microbial communities.

Parameter Sensitivity of Transient Community Dynamics

AbstractTransient dynamics have always intrigued ecologists, but current rapid environmental change (inducing transients even in previously undisturbed systems) has highlighted their importance more than ever. Here, I introduce a method for analyzing the sensitivity of transient ecological dynamics to parameter perturbations. The question the method answers is: how would the community dynamics have unfolded for some time horizon had the parameters been slightly different? I apply the method to three empirically parameterized models: competition between native forbs and exotic grasses in California, a host-parasitoid system, and an experimental chemostat predator-prey model. These applications showcase the ecological insights one can gain from models using transient sensitivity analysis. First, one can find parameters and their combinations whose perturbations disproportionately affect a system. Second, one can identify particular windows of time during which the predicted deviation from the unperturbed trajectories is especially large and utilize this information for management purposes. Third, there is an inverse relationship between transient and long-term sensitivities whenever the interacting populations are ecologically similar; paradoxically, the smaller the immediate response of the system, the more extreme its long-term response will be.

Multitrophic Higher-Order Interactions Modulate Species Persistence

AbstractEcologists increasingly recognize that interactions between two species can be affected by the density of a third species. How these higher-order interactions (HOIs) affect species persistence remains poorly understood. To explore the effect of HOIs stemming from multiple trophic layers on a plant community composition, we experimentally built a mesocosm with three plants and three pollinator species arranged in a fully nested and modified network structure. We estimated pairwise interactions among plants and between plants and pollinators, as well as HOIs initiated by a plant or a pollinator affecting plant species pairs. Using a structuralist approach, we evaluated the consequences of the statistically supported HOIs on the persistence probability of each of the three competing plant species and their combinations. HOIs substantially redistribute the strength and sign of pairwise interactions between plant species, promoting the opportunities for multispecies communities to persist compared with a non-HOI scenario. However, the physical elimination of a plant-pollinator link in the modified network structure promotes changes in per capita pairwise interactions and HOIs, resulting in a single-species community. Our study provides empirical evidence of the joint importance of HOIs and network structure in determining species persistence within diverse communities.

The three-species problem: Incorporating competitive asymmetry and intransitivity in modern coexistence theory

While natural communities can contain hundreds of species, modern coexistence theory focuses primarily on species pairs. Alternatively, the structural stability approach considers the feasibility of equilibria, gaining scalability to larger communities but sacrificing information about dynamic stability. Three-species competitive communities are a bridge to more-diverse communities. They display novel phenomena while remaining amenable to mathematical analysis, but remain incompletely understood. Here, we combine these approaches to identify the key quantities that determine three-species competition outcomes. We show that pairwise niche overlap and fitness differences are insufficient to completely characterize competitive outcomes, which requires a strictly triplet-wise quantity: cyclic asymmetry, which underlies intransitivity. Low pairwise niche overlap stabilizes the triplet, while high fitness differences promote competitive exclusion. The effect of cyclic asymmetry on stability is complex and depends on pairwise niche overlap. In summary, we elucidate how pairwise niche overlap, fitness differences and cyclic asymmetry determine three-species competition outcomes.

Trophic tug-of-war: Coexistence mechanisms within and across trophic levels

Ecological communities encompass rich diversity across multiple trophic levels. While modern coexistence theory has been widely applied to understand community assembly, its traditional formalism only allows assembly within a single trophic level. Here, using an expanded definition of niche and fitness differences applicable to multitrophic communities, we study how diversity within and across trophic levels affects species coexistence. If each trophic level is analysed separately, both lower- and higher trophic levels are governed by the same coexistence mechanisms. In contrast, if the multitrophic community is analysed as a whole, different trophic levels are governed by different coexistence mechanisms: coexistence at lower trophic levels is predominantly limited by fitness differences, whereas coexistence at higher trophic levels is predominantly limited by niche differences. This dichotomy in coexistence mechanisms is supported by theoretical derivations, simulations of phenomenological and trait-based models, and a case study of a primeval forest ecosystem. Our work provides a general and testable prediction of coexistence mechanism operating in multitrophic communities.

Reciprocal facilitation between ants and small mammals in tidal marshes

The role of facilitation in shaping natural communities has primarily been studied in the context of plant assemblages, while its relevance for mobile animals remains less understood. Our study investigates whether reciprocal interspecific facilitation may exist between fire ants (Solenopsis richteri) and cavies (Cavia aperea), two mobile animals, in the SW Atlantic coast brackish marshes. Field samples showed a spatial association between ant mounds and cavies, and that ants prefer to use cavy runways for movement within the marsh. Through experiments involving transplanting the dominant plant, cordgrass (Spartina densiflora), and manipulating cavy presence in areas with and without ant mounds, we observed that cavies forage extensively (and defecate more) near ant mounds. The ants actively remove cavy droppings in their mound vicinity. These ant activities and interactions with cavy droppings led to reduced moisture and organic content while increasing nitrate and phosphate levels in marsh sediment. Consequently, this enhanced plant growth, indirectly facilitating the cavies, which preferred consuming vegetation near ant mounds. These cascading indirect effects persisted over time; even four months after cavies left the marshes, transplanted plants near ant mounds remained larger and exhibited more leaf senescence when exposed to cavy herbivory. Therefore, the networks of positive interactions appear to generate simultaneous selection among species (populations), promoting coexistence within the community. Although complex, these reciprocal facilitative effects among mobile animals may be more common than currently believed and should be further studied to gain a better understanding of the underlying mechanisms driving species coexistence in natural communities.

Processes governing species richness in communities exposed to temporal environmental stochasticity: A review and synthesis of modelling approaches

Research into the processes governing species richness has often assumed that the environment is fixed, whereas realistic environments are often characterised by random fluctuations over time. This temporal environmental stochasticity (TES) changes the demographic rates of species populations, with cascading effects on community dynamics and species richness. Theoretical and applied studies have used process-based mathematical models to determine how TES affects species richness, but under a variety of frameworks. Here, we critically review such studies to synthesise their findings and draw general conclusions. We first provide a broad mathematical framework encompassing the different ways in which TES has been modelled. We then review studies that have analysed models with TES under the assumption of negligible interspecific interactions, such that a community is conceptualised as the sum of independent species populations. These analyses have highlighted how TES can reduce species richness by increasing the frequency at which a species becomes rare and therefore prone to extinction. Next, we review studies that have relaxed the assumption of negligible interspecific interactions. To simplify the corresponding models and make them analytically tractable, such studies have used mean-field theory to derive fixed parameters representing the typical strength of interspecific interactions under TES. The resulting analyses have highlighted community-level effects that determine how TES affects species richness, for species that compete for a common limiting resource. With short temporal correlations of environmental conditions, a non-linear averaging effect of interspecific competition strength over time gives an increase in species richness. In contrast, with long temporal correlations of environmental conditions, strong selection favouring the fittest species between changes in environmental conditions results in a decrease in species richness. We compare such results with those from invasion analysis, which examines invasion growth rates (IGRs) instead of species richness directly. Qualitative differences sometimes arise because the IGR is the expected growth rate of a species when it is rare, which does not capture the variation around this mean or the probability of the species becoming rare. Our review elucidates key processes that have been found to mediate the negative and positive effects of TES on species richness, and by doing so highlights key areas for future research.

Plant interaction networks reveal the limits of our understanding of diversity maintenance

Species interactions are key drivers of biodiversity and ecosystem stability. Current theoretical frameworks for understanding the role of interactions make many assumptions which unfortunately, do not always hold in natural, diverse communities. This mismatch extends to annual plants, a common model system for studying coexistence, where interactions are typically averaged across environmental conditions and transitive competitive hierarchies are assumed to dominate. We quantify interaction networks for a community of annual wildflowers in Western Australia across a natural shade gradient at local scales. Whilst competition dominated, intraspecific and interspecific facilitation were widespread in all shade categories. Interaction strengths and directions varied substantially despite close spatial proximity and similar levels of local species richness, with most species interacting in different ways under different environmental conditions. Contrary to expectations, all networks were predominantly intransitive. These findings encourage us to rethink how we conceive of and categorize the mechanisms driving biodiversity in plant systems.

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.

Robustness versus productivity during evolutionary community assembly: short-term synergies and long-term trade-offs

The realization that evolutionary feedbacks need to be considered to fully grasp ecological dynamics has sparked interest in the effect of evolution on community properties like coexistence and productivity. However, little is known about the evolution of community robustness and productivity along diversification processes in species-rich systems. We leverage the recent structural approach to coexistence together with adaptive dynamics to study such properties and their relationships in a general trait-based model of competition on a niche axis. We show that the effects of coevolution on coexistence are two-fold and contrasting depending on the time scale considered. In the short term, evolution of niche differentiation strengthens coexistence, while long-term diversification leads to niche packing and decreased robustness. Moreover, we find that coevolved communities tend to be on average more robust and more productive than non-evolutionary assemblages. We illustrate how our theoretical predictions echo in observed empirical patterns and the implications of our results for empiricists and applied ecologists. We suggest that some of our results such as the improved robustness of Evolutionarily Stable Communities could be tested experimentally in suitable model systems.

Cell size-dependent species sensitivity to nanoparticles underlies changes in phytoplankton diversity and productivity

Nanoparticle pollution has been shown to affect various organisms. However, the effects of nanoparticles on species interactions, and the role of species traits, such as body size, in modulating these effects, are not well-understood. We addressed this issue using competing freshwater phytoplankton species exposed to copper oxide nanoparticles. Increasing nanoparticle concentration resulted in decreased phytoplankton species growth rates and community productivity (both abundance and biomass). Importantly, we consistently found that nanoparticles had greater negative effects on species with smaller cell sizes, such that nanoparticle pollution weakened the competitive dominance of smaller species and promoted species diversity. Moreover, nanoparticles reduced the growth rate differences and competitive ability differences of competing species, while having little effect on species niche differences. Consequently, nanoparticle pollution reduced the selection effect on phytoplankton community abundance, but increased the selection effect on community biomass. Our results suggest cell size as a key functional trait to consider when predicting phytoplankton community structure and ecosystem functioning in the face of increasing nanopollution.

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.

Assembly archetypes in ecological communities

An instrumental discovery in comparative and developmental biology is the existence of assembly archetypes that synthesize the vast diversity of organisms' body plans-from legs and wings to human arms-into simple, interpretable and general design principles. Here, we combine a novel mathematical formalism based on category theory with experimental data to show that similar 'assembly archetypes' exist at the larger organization scale of ecological communities when assembling a species pool across diverse environmental contexts, particularly when species interactions are highly structured. We applied our formalism to clinical data discovering two assembly archetypes that differentiate between healthy and unhealthy human gut microbiota. The concept of assembly archetypes and the methods to synthesize them can pave the way to discovering the general assembly principles of the ecological communities we observe in nature.

Structural asymmetry in biotic interactions as a tool to understand and predict ecological persistence

A universal feature of ecological systems is that species do not interact with others with the same sign and strength. Yet, the consequences of this asymmetry in biotic interactions for the short- and long-term persistence of individual species and entire communities remains unclear. Here, we develop a set of metrics to evaluate how asymmetric interactions among species translate to asymmetries in their individual vulnerability to extinction under changing environmental conditions. These metrics, which solve previous limitations of how to independently quantify the size from the shape of the so-called feasibility domain, provide rigorous advances to understand simultaneously why some species and communities present more opportunities to persist than others. We further demonstrate that our shape-related metrics are useful to predict short-term changes in species' relative abundances during 7 years in a Mediterranean grassland. Our approach is designed to be applied to any ecological system regardless of the number of species and type of interactions. With it, we show that is possible to obtain both mechanistic and predictive information on ecological persistence for individual species and entire communities, paving the way for a stronger integration of theoretical and empirical research.

A spatially structured mathematical model of the gut microbiome reveals factors that increase community stability

Given the importance of gut microbial communities for human health, we may want to ensure their stability in terms of species composition and function. Here, we built a mathematical model of a simplified gut composed of two connected patches where species and metabolites can flow from an upstream patch, allowing upstream species to affect downstream species' growth. First, we found that communities in our model are more stable if they assemble through species invasion over time compared to combining a set of species from the start. Second, downstream communities are more stable when species invade the downstream patch less frequently than the upstream patch. Finally, upstream species that have positive effects on downstream species can further increase downstream community stability. Despite it being quite abstract, our model may inform future research on designing more stable microbial communities or increasing the stability of existing ones.

Mean species responses predict effects of environmental change on coexistence

Environmental change research is plagued by the curse of dimensionality: the number of communities at risk and the number of environmental drivers are both large. This raises the pressing question if a general understanding of ecological effects is achievable. Here, we show evidence that this is indeed possible. Using theoretical and simulation-based evidence for bi- and tritrophic communities, we show that environmental change effects on coexistence are proportional to mean species responses and depend on how trophic levels on average interact prior to environmental change. We then benchmark our findings using relevant cases of environmental change, showing that means of temperature optima and of species sensitivities to pollution predict concomitant effects on coexistence. Finally, we demonstrate how to apply our theory to the analysis of field data, finding support for effects of land use change on coexistence in natural invertebrate communities.

Periodic temporal environmental variations induce coexistence in resource competition models

Natural ecosystems, in particular on the microbial scale, are inhabited by a large number of species. The population size of each species is affected by interactions of individuals with each other and by spatial and temporal changes in environmental conditions, such as resource abundance. Here, we use a generic population dynamics model to study how, and under what conditions, a periodic temporal environmental variation can alter an ecosystem's composition and biodiversity. We demonstrate that using timescale separation allows one to qualitatively predict the long-term population dynamics of interacting species in varying environments. We show that the notion of Tilman's R* rule, a well-known principle that applies for constant environments, can be extended to periodically varying environments if the timescale of environmental changes (e.g., seasonal variations) is much faster than the timescale of population growth (doubling time in bacteria). When these timescales are similar, our analysis shows that a varying environment deters the system from reaching a steady state, and stable coexistence between multiple species becomes possible. Our results posit that biodiversity can in part be attributed to natural environmental variations.

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.

Mechanistic Models of Trophic Interactions: Opportunities for Species Richness and Challenges for Modern Coexistence Theory

AbstractMany potential mechanisms promote species coexistence, but we know little about their relative importance. To compare multiple mechanisms, we modeled a two-trophic planktonic food web based on mechanistic species interactions and empirically measured species traits. We simulated thousands of possible communities under realistic and altered interaction strengths to assess the relative importance of three potential drivers of phytoplankton and zooplankton species richness: resource-mediated coexistence mechanisms, predator-prey interactions, and trait trade-offs. Next, we computed niche and fitness differences of competing zooplankton to obtain a deeper understanding of how these mechanisms determine species richness. We found that predator-prey interactions were the most important driver of phytoplankton and zooplankton species richness and that large zooplankton fitness differences were associated with low species richness, but zooplankton niche differences were not associated with species richness. However, for many communities we could not apply modern coexistence theory to compute niche and fitness differences of zooplankton because of conceptual issues with the invasion growth rates arising from trophic interactions. We therefore need to expand modern coexistence theory to fully investigate multitrophic-level communities.

Contextualizing the ecology of plant-plant interactions and constructive networks

Botanical concepts have traditionally viewed the environment as a static box containing plants. In this box, plants compete with one another and act as passive resource consumers subjected to the environment in a top-down manner. This entails that plants have only negative effects on other plants and have no influence on the environment. By contrast, there is increasing evidence that plants have positive, bottom-up engineering effects and diversity effects on other plants and on the environment. Here, to overcome the limitations of top-down environmental control, antagonistic-only and pairwise interactions, I propose the concept of constructive networks. Constructive networks unify niche construction and network theory recognizing that (i) plants have manifold ecological functions and impacts on their neighbours, and (ii) the environment shapes and is shaped by diverse organisms, primarily plants. Constructive networks integrate both plant-environment and plant-plant interactions in a relational context. They address how plants influence the environment and support or inhibit other plant species by physically, biochemically and ecologically shaping environmental conditions. Constructive networks acknowledge the fact that diverse plants change and create novel environmental conditions and co-produce, share and transform resources, thereby influencing biological communities and the environment in constructive ways. Different interaction types are considered simultaneously in constructive networks. Yet, the main limitation to understanding constructive networks is the identification of plant links. This barrier may be overcome by applying complexity theory and statistical mechanics to comparative data and experimental field botany. Considering multiple interaction types and feedback between plants and the environment may improve our understanding of mechanisms responsible for biodiversity maintenance and help us to better anticipate the response of plant systems to global change.

Relationships of temperature and biodiversity with stability of natural aquatic food webs

Temperature and biodiversity changes occur in concert, but their joint effects on ecological stability of natural food webs are unknown. Here, we assess these relationships in 19 planktonic food webs. We estimate stability as structural stability (using the volume contraction rate) and temporal stability (using the temporal variation of species abundances). Warmer temperatures were associated with lower structural and temporal stability, while biodiversity had no consistent effects on either stability property. While species richness was associated with lower structural stability and higher temporal stability, Simpson diversity was associated with higher temporal stability. The responses of structural stability were linked to disproportionate contributions from two trophic groups (predators and consumers), while the responses of temporal stability were linked both to synchrony of all species within the food web and distinctive contributions from three trophic groups (predators, consumers, and producers). Our results suggest that, in natural ecosystems, warmer temperatures can erode ecosystem stability, while biodiversity changes may not have consistent effects.

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.

Navigation between initial and desired community states using shortcuts

Ecological management problems often involve navigating from an initial to a desired community state. We ask whether navigation without brute-force additions and deletions of species is possible via: adding/deleting a small number of individuals of a species, changing the environment, and waiting. Navigation can yield direct paths (single sequence of actions) or shortcut paths (multiple sequences of actions with lower cost than a direct path). We ask (1) when is non-brute-force navigation possible?; (2) do shortcuts exist and what are their properties?; and (3) what heuristics predict shortcut existence? Using a state diagram framework applied to several empirical datasets, we show that (1) non-brute-force navigation is only possible between some state pairs, (2) shortcuts exist between many state pairs; and (3) changes in abundance and richness are the strongest predictors of shortcut existence, independent of dataset and algorithm choices. State diagrams thus unveil hidden strategies for manipulating species coexistence and efficiently navigating between states.

Non-random interactions within and across guilds shape the potential to coexist in multi-trophic ecological communities

Theory posits that the persistence of species in ecological communities is shaped by their interactions within and across trophic guilds. However, we lack empirical evaluations of how the structure, strength and sign of biotic interactions drive the potential to coexist in diverse multi-trophic communities. Here, we model community feasibility domains, a theoretically informed measure of multi-species coexistence probability, from grassland communities comprising more than 45 species on average from three trophic guilds (plants, pollinators and herbivores). Contrary to our hypothesis, increasing community complexity, measured either as the number of guilds or community richness, did not decrease community feasibility. Rather, we observed that high degrees of species self-regulation and niche partitioning allow for maintaining larger levels of community feasibility and higher species persistence in more diverse communities. Our results show that biotic interactions within and across guilds are not random in nature and both structures significantly contribute to maintaining multi-trophic diversity.

Nitrogen enrichment and foliar fungal pathogens affect the mechanisms of multispecies plant coexistence

Changes in resources (e.g. nitrogen) and enemies (e.g. foliar pathogens) are key drivers of plant diversity and composition. However, their effects have not been connected to the niche and fitness differences that determine multispecies coexistence. Here, we combined a structuralist theoretical approach with a detailed grassland experiment factorially applying nitrogen addition and foliar fungal pathogen suppression to evaluate the joint effect of nitrogen and pathogens on niche and fitness differences, across a gradient from two to six interacting species. Nitrogen addition and pathogen suppression modified species interaction strengths and intrinsic growth rates, leading to reduced multispecies fitness differences. However, contrary to expected, we also observed that they promote stabilising niche differences. Although these modifications did not substantially alter species richness, they predicted major changes in community composition. Indirect interactions between species explained these community changes in smaller assemblages (three and four species) but lost importance in favour of direct pairwise interactions when more species were involved (five and six). Altogether, our work shows that explicitly considering the number of interacting species is critical for better understanding the direct and indirect processes by which nitrogen enrichment and pathogen communities shape coexistence in grasslands.

Niche theory for positive plant-soil feedbacks

Interactions between plants and the soil are an important ecological process in terrestrial ecosystems as they affect plant community structure: when and where we find different plant species. Those interactions are typically thought of as one-directional: local soil conditions filter through dispersing species to produce a community of locally adapted plants. However, plants can modify local physicochemical soil conditions via their roots and associations with soil microbes. These may in turn affect the local fitness of other plants, making plant-soil interactions bidirectional. In order to understand how they differ from other ecological processes that structure plant communities, we need a theory connecting these individual-level plant-soil feedbacks to community-level patterns. Here, we build this theory with a mathematical model of plant community dynamics in which soil conditioning is explicitly modeled over time and depends on the density of the plants. We analyze this model to describe the long-term composition and spatial distribution of the plant community. Our main result is that positive plant-soil feedbacks will create clustering of species with similar soil preferences. The composition of these clusters is further influenced by niche width and conditioning strength. In contrast with competitive dynamics driven by niche overlap, only species belonging to the same cluster can maintain high relative abundance in the community. Spatial heterogeneity in the form of an environmental gradient generates patches, each representing a single cluster. However, such patchiness is disfavored when species differ in dispersal ability. We show that stronger dispersers cannot take over the habitat as long as an exogenous driver favors soil conditions that benefit the other species. If exogenous drivers supersede soil conditioning by plants, we retrieve classic habitat filtering, where species are selected based on their suitability to the local environment. Overall, we provide a novel mathematical model for positive plant-soil feedback that we use to describe the spatial patterns of plant abundance and traits related to soil preference and conditioning ability.

Heterogeneity of interaction strengths and its consequences on ecological systems

Ecosystems are formed by networks of species and their interactions. Traditional models of such interactions assume a constant interaction strength between a given pair of species. However, there is often significant trait variation among individual organisms even within the same species, causing heterogeneity in their interaction strengths with other species. The consequences of such heterogeneous interactions for the ecosystem have not been studied systematically. As a theoretical exploration, we analyze a simple ecosystem with trophic interactions between two predators and a shared prey, which would exhibit competitive exclusion in models with homogeneous interactions. We consider several scenarios where individuals of the prey species differentiate into subpopulations with different interaction strengths. We show that in all these cases, whether the heterogeneity is inherent, reversible, or adaptive, the ecosystem can stabilize at a new equilibrium where all three species coexist. Moreover, the prey population that has heterogeneous interactions with its predators reaches a higher density than it would without heterogeneity, and can even reach a higher density in the presence of two predators than with just one. Our results suggest that heterogeneity may be a naturally selected feature of ecological interactions that have important consequences for the stability and diversity of ecosystems.

Effects of multi-resource addition on grassland plant productivity and biodiversity along a resource gradient

The change of plant biodiversity caused by resource-enhancing global changes has greatly affected grassland productivity. However, it remains unclear how multi-resource enrichment induces the effects of multifaceted biodiversity on grassland productivity under different site resource constraints. We conducted a multiple resource addition (MRA) experiment of water and nutrients at three sites located along a resource gradient in northern China. This allowed us to assess the response of aboveground net primary productivity (ANPP), species (species richness and plant density), functional (functional richness and community-weighted mean of traits) and phylogenetic (phylogenetic richness) diversity to increasing number of MRA. We used structural equation model (SEM) to examine the direct and indirect effects of MRA and multifaceted biodiversity on ANPP. The combined addition of the four resources increased ANPP at all three sites. But with increasing number of MRA, biodiversity varied at the three sites. At the high resource constraint site, species richness, plant density and leaf nitrogen concentration (LNC) increased. At the medium resource constraint site, plant height and LNC increased, leaf dry matter content (LDMC) decreased. At the low resource constraint site, species, functional and phylogenetic richness decreased, and height increased. The SEM showed that MRA increased ANPP directly at all three sites, and indirectly by increasing plant density at the high constraint site and height at the medium constraint site. Independent of MRA, ANPP was affected by height at the high resource constraint site and LNC at the low resource constraint site. Our results illustrate that multi-resource addition positively affects productivity, while affects biodiversity depending on site resource constraint. The study highlights that site resource constraint conditions need to be taken into consideration to better predict grassland structure and function, particularly under the future multifaceted global change scenarios.

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.

EcoEvoApps: Interactive apps for theoretical models in ecology and evolutionary biology

The integration of theory and data drives progress in science, but a persistent barrier to such integration in ecology and evolutionary biology is that theory is often developed and expressed in the form of mathematical models that can feel daunting and inaccessible for students and empiricists with variable quantitative training and attitudes towards math. A promising way to make mathematical models more approachable is to embed them into interactive tools with which one can visually evaluate model structures and directly explore model outcomes through simulation. To promote such interactive learning of quantitative models, we developed EcoEvoApps, a collection of free, open-source, and multilingual R/Shiny apps that include model overviews, interactive model simulations, and code to implement these models directly in R. The package currently focuses on canonical models of population dynamics, species interactions, and landscape ecology. These apps help illustrate fundamental results from theoretical ecology and can serve as valuable teaching tools in classroom settings. We present data from student surveys which show that students rate these apps as useful learning tools, and that using interactive apps leads to substantial gains in students' interest and confidence in working with mathematical models. This points to the potential for interactive activities to make theoretical models more accessible to a wider audience, and thus facilitate the feedback between theory and data across ecology and evolutionary biology.

Small rainfall changes drive substantial changes in plant coexistence

Although precipitation patterns have long been known to shape plant distributions1, the effect of changing climate on the interactions of species and therefore community composition is far less understood2,3. Here, we explored how changes in precipitation alter competitive dynamics via direct effects on individual species, as well as by the changing strength of competitive interactions between species, using an annual grassland community in California. We grew plants under ambient and reduced precipitation in the field to parameterize a competition model4 with which we quantified the stabilizing niche and fitness differences that determine species coexistence in each rainfall regime. We show that reduced precipitation had little direct effect on species grown alone, but it qualitatively shifted predicted competitive outcomes for 10 of 15 species pairs. In addition, species pairs that were functionally more similar were less likely to experience altered outcomes, indicating that functionally diverse communities may be most threatened by changing interactions. Our results highlight how important it is to account for changes to species interactions when predicting species and community response to global change.

Coexistence in diverse communities with higher-order interactions

A central assumption in most ecological models is that the interactions in a community operate only between pairs of species. However, two species may interactively affect the growth of a focal species. Although interactions among three or more species, called higher-order interactions, have the potential to modify our theoretical understanding of coexistence, ecologists lack clear expectations for how these interactions shape community structure. Here we analytically predict and numerically confirm how the variability and strength of higher-order interactions affect species coexistence. We found that as higher-order interaction strengths became more variable across species, fewer species could coexist, echoing the behavior of pairwise models. If interspecific higher-order interactions became too harmful relative to self-regulation, coexistence in diverse communities was destabilized, but coexistence was also lost when these interactions were too weak and mutualistic higher-order effects became prevalent. This behavior depended on the functional form of the interactions as the destabilizing effects of the mutualistic higher-order interactions were ameliorated when their strength saturated with species' densities. Last, we showed that more species-rich communities structured by higher-order interactions lose species more readily than their species-poor counterparts, generalizing classic results for community stability. Our work provides needed theoretical expectations for how higher-order interactions impact species coexistence in diverse communities.

Application of modern coexistence theory to rare plant restoration provides early indication of restoration trajectories

Restoration ecology commonly seeks to re-establish species of interest in degraded habitats. Despite a rich understanding of how succession influences re-establishment, there are several outstanding questions that remain unaddressed: are short-term abundances sufficient to determine long-term re-establishment success, and what factors contribute to unpredictable restorations outcomes? In other words, when restoration fails, is it because the restored habitat is substandard, because of strong competition with invasive species, or alternatively due to changing environmental conditions that would equally impact established populations? Here, we re-purpose tools developed from modern coexistence theory to address these questions, and apply them to an effort to restore the endangered Contra Costa goldfields (Lasthenia conjugens) in constructed ("restored") California vernal pools. Using 16 years of data, we construct a population model of L. conjugens, a species of conservation concern due primarily to habitat loss and invasion of exotic grasses. We show that initial, short-term appearances of restoration success from population abundances is misleading, as year-to-year fluctuations cause long-term population growth rates to fall below zero. The failure of constructed pools is driven by lower maximum growth rates compared with reference ("natural") pools, coupled with a stronger negative sensitivity to annual fluctuations in abiotic conditions that yield decreased maximum growth rates. Nonetheless, our modeling shows that fluctuations in competition (mainly with exotic grasses) benefit L. conjugens through periods of competitive release, especially in constructed pools of intermediate pool depth. We therefore show how reductions in invasives and seed addition in pools of particular depths could change the outcome of restoration for L. conjugens. By applying a largely theoretical framework to the urgent goal of ecological restoration, our study provides a blueprint for predicting restoration success, and identifies future actions to reverse species loss.

Integrating eco-evolutionary dynamics and modern coexistence theory

Community ecology typically assumes that competitive exclusion and species coexistence are unaffected by evolution on the time scale of ecological dynamics. However, recent studies suggest that rapid evolution operating concurrently with competition may enable species coexistence. Such findings necessitate general theory that incorporates the coexistence contributions of eco-evolutionary processes in parallel with purely ecological mechanisms and provides metrics for quantifying the role of evolution in shaping competitive outcomes in both modelling and empirical contexts. To foster the development of such theory, here we extend the interpretation of the two principal metrics of modern coexistence theory-niche and competitive ability differences-to systems where competitors evolve. We define eco-evolutionary versions of these metrics by considering how invading and resident species adapt to conspecific and heterospecific competitors. We show that the eco-evolutionary niche and competitive ability differences are sums of ecological and evolutionary processes, and that they accurately predict the potential for stable coexistence in previous theoretical studies of eco-evolutionary dynamics. Finally, we show how this theory frames recent empirical assessments of rapid evolution effects on species coexistence, and how empirical work and theory on species coexistence and eco-evolutionary dynamics can be further integrated.

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.

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.

Stability criteria for the consumption and exchange of essential resources

Models of consumer effects on a shared resource environment have helped clarify how the interplay of consumer traits and resource supply impact stable coexistence. Recent models generalize this picture to include the exchange of resources alongside resource competition. These models exemplify the fact that although consumers shape the resource environment, the outcome of consumer interactions is context-dependent: such models can have either stable or unstable equilibria, depending on the resource supply. However, these recent models focus on a simplified version of microbial metabolism where the depletion of resources always leads to consumer growth. Here, we model an arbitrarily large system of consumers governed by Liebig’s law, where species require and deplete multiple resources, but each consumer’s growth rate is only limited by a single one of these resources. Resources that are taken up but not incorporated into new biomass are leaked back into the environment, possibly transformed by intracellular reactions, thereby tying the mismatch between depletion and growth to cross-feeding. For this set of dynamics, we show that feasible equilibria can be either stable or unstable, again depending on the resource environment. We identify special consumption and production networks which protect the community from instability when resources are scarce. Using simulations, we demonstrate that the qualitative stability patterns derived analytically apply to a broader class of network structures and resource inflow profiles, including cases where multiple species coexist on only one externally supplied resource. Our stability criteria bear some resemblance to classic stability results for pairwise interactions, but also demonstrate how environmental context can shape coexistence patterns when resource limitation and exchange are modeled directly.

Gut microbiomes of cyprinid fish exhibit host-species symbiosis along gut trait and diet

Teleost omnivorous fish that coexist partially sharing resources are likely to modify their gut traits and microbiome as a feedback mechanism between ecological processes and evolution. However, we do not understand how the core gut microbiome supports the metabolic capacity of the host and regulates digestive functions in specialized omnivorous fish gut traits. Therefore, we evaluated the gut microbiome of eight omnivorous fish from a single family (i.e., Cyprinidae) in the current study. We examined the correlation between host phylogeny, diet composition, and intestinal morphological traits related to the intestinal microbiome. The results indicated that cyprinid fish with similar relative gut lengths had considerable gut microbiome similarity. Notably, the SL (short relative gut length) group, as zoobenthos and zooplankton specialists, was abundant in Proteobacteria and was less abundant in Firmicutes than in the ML (medium relative gut length) and LL (long relative gut length) groups. These fish could extract nutrients from aquatic plants and algae. Additionally, we found the relative abundance of Clostridium and Romboutsia to be positively correlated with host relative gut length but negatively correlated with the relative abundance of Cetobacterium, Plesiomonas, Bacteroides, and Lactobacillus, and host-relative gut length. We also show a positive linear relationship between host gut microbiome carbohydrate metabolism and relative gut length, while the amino acid and lipid metabolism of the gut microbiome was negatively correlated with host-relative gut length. In addition, omnivorous species competing for resources improve their ecological adaptability through the specialization of gut length, which is closely related to variation in the synergy of the gut microbiome. Above all, specialized gut microbiota and associated gut morphologies enable fish to variably tolerate resource fluctuation and improve the utilization efficiency of nutrient extraction from challenging food resources.

Quantifying mechanisms of coexistence in disease ecology

Pathogen coexistence depends on ecological processes operating at both within and between-host scales, making it difficult to quantify which processes may promote or prevent coexistence. Here, we propose that adapting modern coexistence theory-traditionally applied in plant communities-to pathogen systems provides an exciting approach for examining mechanisms of coexistence operating across different spatial scales. We first overview modern coexistence theory and its mechanistic decomposition; we subsequently adapt the framework to quantify how spatial variation in pathogen density, host resources and immunity, and their interaction may promote pathogen coexistence. We apply this derivation to an example two pathogen, multi-scale model comparing two scenarios with generalist and strain-specific immunity: one with demographic equivalency among pathogens and one with demographic trade-offs among pathogens. We then show how host-pathogen feedbacks generate spatial heterogeneity that promote pathogen coexistence and decompose those mechanisms to quantify how each spatial heterogeneity contributes to that coexistence. Specifically, coexistence of demographically equivalent pathogens occurs due to spatial variation in host resources, immune responses, and pathogen aggregation. With a competition-colonization trade-off, the superior colonizer requires spatial heterogeneity to coexist, whereas the superior competitor does not. Finally, we suggest ways forward for linking theory and empirical tests of coexistence in disease systems.

Competition for pollinators destabilizes plant coexistence

Mounting concern over the global decline of pollinators has fuelled calls for investigating their role in maintaining plant diversity^1,2. Theory predicts that competition for pollinators can stabilize interactions between plant species by providing opportunities for niche differentiation³, while at the same time can drive competitive imbalances that favour exclusion⁴. Here we empirically tested these contrasting effects by manipulating competition for pollinators in a way that predicts its long-term implications for plant coexistence. We subjected annual plant individuals situated across experimentally imposed gradients in neighbour density to either ambient insect pollination or a pollen supplementation treatment alleviating competition for pollinators. The vital rates of these individuals informed plant population dynamic models predicting the key theoretical metrics of species coexistence. Competition for pollinators generally destabilized the interactions between plant species, reducing the proportion of pairs expected to coexist. Interactions with pollinators also influenced the competitive imbalances between plant species, effects that are expected to strengthen with pollinator decline, potentially disrupting plant coexistence. Indeed, results from an experiment simulating pollinator decline showed that plant species experiencing greater reductions in floral visitation also suffered greater declines in population growth rate. Our results reveal that competition for pollinators may weaken plant coexistence by destabilizing interactions and contributing to competitive imbalances, information critical for interpreting the impacts of pollinator decline.