ABSTRACT
Structure, composition, and stability of ecological populations are shaped by the inter- and intraspecies interactions within their communities. It remains to be fully understood how the interplay of these interactions with other factors, such as immigration, controls the structure, the diversity, and the long-term stability of ecological systems in the presence of noise and fluctuations. We address this problem using a minimal model of interacting multispecies ecological communities that incorporates competition, immigration, and demographic noise. We find that a complete phase diagram exhibits rich behavior with multiple regimes that go beyond the classical "niche" and "neutral" regimes, extending and modifying the "rare biosphere" or "niche-like" dichotomy. In particular, we observe regimes that cannot be characterized as either niche or neutral where a multimodal species abundance distribution is observed. We characterize the transitions between the different regimes and show how these arise from the underlying kinetics of the species turnover, extinction, and invasion. Our model serves as a minimal null model of noisy competitive ecological systems, against which more complex models that include factors such as mutations and environmental noise can be compared.
Subject(s)
Ecosystem , Models, Biological , Biodiversity , Biota , Kinetics , Population DynamicsABSTRACT
Quantifying biochemical reaction rates within complex cellular processes remains a key challenge of systems biology even as high-throughput single-cell data have become available to characterize snapshots of population variability. That is because complex systems with stochastic and non-linear interactions are difficult to analyze when not all components can be observed simultaneously and systems cannot be followed over time. Instead of using descriptive statistical models, we show that incompletely specified mechanistic models can be used to translate qualitative knowledge of interactions into reaction rate functions from covariability data between pairs of components. This promises to turn a globally intractable problem into a sequence of solvable inference problems to quantify complex interaction networks from incomplete snapshots of their stochastic fluctuations.