Ecological landscapes guide the assembly of optimal microbial communities

Global epistasis and the emergence of function in microbial consortia

The many functions of microbial communities emerge from a complex web of interactions between organisms and their environment. This poses a significant obstacle to engineering microbial consortia, hindering our ability to harness the potential of microorganisms for biotechnological applications. In this study, we demonstrate that the collective effect of ecological interactions between microbes in a community can be captured by simple statistical models that predict how adding a new species to a community will affect its function. These predictive models mirror the patterns of global epistasis reported in genetics, and they can be quantitatively interpreted in terms of pairwise interactions between community members. Our results illuminate an unexplored path to quantitatively predicting the function of microbial consortia from their composition, paving the way to optimizing desirable community properties and bringing the tasks of predicting biological function at the genetic, organismal, and ecological scales under the same quantitative formalism.

Sparsity of higher-order landscape interactions enables learning and prediction for microbiomes

Microbiome engineering offers the potential to leverage microbial communities to improve outcomes in human health, agriculture, and climate. To translate this potential into reality, it is crucial to reliably predict community composition and function. But a brute force approach to cataloging community function is hindered by the combinatorial explosion in the number of ways we can combine microbial species. An alternative is to parameterize microbial community outcomes using simplified, mechanistic models, and then extrapolate these models beyond where we have sampled. But these approaches remain data-hungry, as well as requiring an a priori specification of what kinds of mechanisms are included and which are omitted. Here, we resolve both issues by introducing a mechanism-agnostic approach to predicting microbial community compositions and functions using limited data. The critical step is the identification of a sparse representation of the community landscape. We then leverage this sparsity to predict community compositions and functions, drawing from techniques in compressive sensing. We validate this approach on in silico community data, generated from a theoretical model. By sampling just [Formula: see text]1% of all possible communities, we accurately predict community compositions out of sample. We then demonstrate the real-world application of our approach by applying it to four experimental datasets and showing that we can recover interpretable, accurate predictions on composition and community function from highly limited data.

Universal abundance fluctuations across microbial communities, tropical forests, and urban populations

The growth of complex populations, such as microbial communities, forests, and cities, occurs over vastly different spatial and temporal scales. Although research in different fields has developed detailed, system-specific models to understand each individual system, a unified analysis of different complex populations is lacking; such an analysis could deepen our understanding of each system and facilitate cross-pollination of tools and insights across fields. Here, we use a shared framework to analyze time-series data of the human gut microbiome, tropical forest, and urban employment. We demonstrate that a single, three-parameter model of stochastic population dynamics can reproduce the empirical distributions of population abundances and fluctuations in all three datasets. The three parameters characterizing a species measure its mean abundance, deterministic stability, and stochasticity. Our analysis reveals that, despite the vast differences in scale, all three systems occupy a similar region of parameter space when time is measured in generations. In other words, although the fluctuations observed in these systems may appear different, this difference is primarily due to the different physical timescales associated with each system. Further, we show that the distribution of temporal abundance fluctuations is described by just two parameters and derive a two-parameter functional form for abundance fluctuations to improve risk estimation and forecasting.

Functional convergence in slow-growing microbial communities arises from thermodynamic constraints

The dynamics of microbial communities is complex, determined by competition for metabolic substrates and cross-feeding of byproducts. Species in the community grow by harvesting energy from chemical reactions that transform substrates to products. In many anoxic environments, these reactions are close to thermodynamic equilibrium and growth is slow. To understand the community structure in these energy-limited environments, we developed a microbial community consumer-resource model incorporating energetic and thermodynamic constraints on an interconnected metabolic network. The central element of the model is product inhibition, meaning that microbial growth may be limited not only by depletion of metabolic substrates but also by accumulation of products. We demonstrate that these additional constraints on microbial growth cause a convergence in the structure and function of the community metabolic network-independent of species composition and biochemical details-providing a possible explanation for convergence of community function despite taxonomic variation observed in many natural and industrial environments. Furthermore, we discovered that the structure of community metabolic network is governed by the thermodynamic principle of maximum free energy dissipation. Our results predict the decrease of functional convergence in faster growing communities, which we validate by analyzing experimental data from anaerobic digesters. Overall, the work demonstrates how universal thermodynamic principles may constrain community metabolism and explain observed functional convergence in microbial communities.

Pulsed, continuous or somewhere in between? Resource dynamics matter in the optimisation of microbial communities

The optimisation of synthetic and natural microbial communities has vast potential for emerging applications in medicine, agriculture and industry. Realising this goal is contingent on a close correlation between theory, experiments, and the real world. Although the temporal pattern of resource supply can play a major role in microbial community assembly, resource dynamics are commonly treated inconsistently in theoretical and experimental research. Here we explore how the composition of communities varies under continuous resource supply, typical of theoretical approaches, versus pulsed resource supply, typical of experiments. Using simulations of classical resource competition models, we show that community composition diverges rapidly between the two regimes, with almost zero overlap in composition once the pulsing interval stretches beyond just four hours. The implication for the rapidly growing field of microbial community optimisation is that the resource supply regime must be tailored to the community being optimised. As such, we argue that resource supply dynamics should be considered both a constraint in the design of novel microbial communities and as a tuning mechanism for the optimisation of pre-existing communities like those found in the human gut.