High-order interactions distort the functional landscape of microbial consortia

Custom scoring based on ecological topology of gut microbiota associated with cancer immunotherapy outcome

The gut microbiota influences the clinical responses of cancer patients to immunecheckpoint inhibitors (ICIs). However, there is no consensus definition of detrimental dysbiosis. Based on metagenomics (MG) sequencing of 245 non-small cell lung cancer (NSCLC) patient feces, we constructed species-level co-abundance networks that were clustered into species-interacting groups (SIGs) correlating with overall survival. Thirty-seven and forty-five MG species (MGSs) were associated with resistance (SIG1) and response (SIG2) to ICIs, respectively. When combined with the quantification of Akkermansia species, this procedure allowed a person-based calculation of a topological score (TOPOSCORE) that was validated in an additional 254 NSCLC patients and in 216 genitourinary cancer patients. Finally, this TOPOSCORE was translated into a 21-bacterial probe set-based qPCR scoring that was validated in a prospective cohort of NSCLC patients as well as in colorectal and melanoma patients. This approach could represent a dynamic diagnosis tool for intestinal dysbiosis to guide personalized microbiota-centered interventions.

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.

Bioactive exometabolites drive maintenance competition in simple bacterial communities

During prolonged resource limitation, bacterial cells can persist in metabolically active states of non-growth. These maintenance periods, such as those experienced in stationary phase, can include upregulation of secondary metabolism and release of exometabolites into the local environment. As resource limitation is common in many environmental microbial habitats, we hypothesized that neighboring bacterial populations employ exometabolites to compete or cooperate during maintenance and that these exometabolite-facilitated interactions can drive community outcomes. Here, we evaluated the consequences of exometabolite interactions over the stationary phase among three environmental strains: Burkholderia thailandensis E264, Chromobacterium subtsugae ATCC 31532, and Pseudomonas syringae pv. tomato DC3000. We assembled them into synthetic communities that only permitted chemical interactions. We compared the responses (transcripts) and outputs (exometabolites) of each member with and without neighbors. We found that transcriptional dynamics were changed with different neighbors and that some of these changes were coordinated between members. The dominant competitor B. thailandensis consistently upregulated biosynthetic gene clusters to produce bioactive exometabolites for both exploitative and interference competition. These results demonstrate that competition strategies during maintenance can contribute to community-level outcomes. It also suggests that the traditional concept of defining competitiveness by growth outcomes may be narrow and that maintenance competition could be an additional or alternative measure.

IMPORTANCE

Free-living microbial populations often persist and engage in environments that offer few or inconsistently available resources. Thus, it is important to investigate microbial interactions in this common and ecologically relevant condition of non-growth. This work investigates the consequences of resource limitation for community metabolic output and for population interactions in simple synthetic bacterial communities. Despite non-growth, we observed active, exometabolite-mediated competition among the bacterial populations. Many of these interactions and produced exometabolites were dependent on the community composition but we also observed that one dominant competitor consistently produced interfering exometabolites regardless. These results are important for predicting and understanding microbial interactions in resource-limited environments.

Modeling Microbial Community Networks: Methods and Tools for Studying Microbial Interactions

Microbial interactions function as a fundamental unit in complex ecosystems. By characterizing the type of interaction (positive, negative, neutral) occurring in these dynamic systems, one can begin to unravel the role played by the microbial species. Towards this, various methods have been developed to decipher the function of the microbial communities. The current review focuses on the various qualitative and quantitative methods that currently exist to study microbial interactions. Qualitative methods such as co-culturing experiments are visualized using microscopy-based techniques and are combined with data obtained from multi-omics technologies (metagenomics, metabolomics, metatranscriptomics). Quantitative methods include the construction of networks and network inference, computational models, and development of synthetic microbial consortia. These methods provide a valuable clue on various roles played by interacting partners, as well as possible solutions to overcome pathogenic microbes that can cause life-threatening infections in susceptible hosts. Studying the microbial interactions will further our understanding of complex less-studied ecosystems and enable design of effective frameworks for treatment of infectious diseases.

A generalized simplicial model and its application

Higher-order structures, consisting of more than two individuals, provide a new perspective to reveal the missed non-trivial characteristics under pairwise networks. Prior works have researched various higher-order networks, but research for evaluating the effects of higher-order structures on network functions is still scarce. In this paper, we propose a framework to quantify the effects of higher-order structures (e.g., 2-simplex) and vital functions of complex networks by comparing the original network with its simplicial model. We provide a simplicial model that can regulate the quantity of 2-simplices and simultaneously fix the degree sequence. Although the algorithm is proposed to control the quantity of 2-simplices, results indicate it can also indirectly control simplexes more than 2-order. Experiments on spreading dynamics, pinning control, network robustness, and community detection have shown that regulating the quantity of 2-simplices changes network performance significantly. In conclusion, the proposed framework is a general and effective tool for linking higher-order structures with network functions. It can be regarded as a reference object in other applications and can deepen our understanding of the correlation between micro-level network structures and global network functions.

High immigration rates critical for establishing emigration-driven diversity in microbial communities

Unraveling the mechanisms governing the diversity of ecological communities is a central goal in ecology. Although microbial dispersal constitutes an important ecological process, the effect of dispersal on microbial diversity is poorly understood. Here, we sought to fill this gap by combining a generalized Lotka-Volterra model with experimental investigations. Our model showed that emigration increases the diversity of the community when the immigration rate crosses a defined threshold, which we identified as Ineutral. We also found that at high immigration rates, emigration weakens the relative abundance of fast-growing species and thus enhances the mass effect and increases the diversity. We experimentally confirmed this finding using co-cultures of 20 bacterial strains isolated from the soil. Our model further showed that Ineutral decreases with the increase of species pool size, growth rate, and interspecies interaction. Our work deepens the understanding of the effects of dispersal on the diversity of natural communities.

Data-driven prediction of colonization outcomes for complex microbial communities

Microbial interactions can lead to different colonization outcomes of exogenous species, be they pathogenic or beneficial in nature. Predicting the colonization of exogenous species in complex communities remains a fundamental challenge in microbial ecology, mainly due to our limited knowledge of the diverse mechanisms governing microbial dynamics. Here, we propose a data-driven approach independent of any dynamics model to predict colonization outcomes of exogenous species from the baseline compositions of microbial communities. We systematically validate this approach using synthetic data, finding that machine learning models can predict not only the binary colonization outcome but also the post-invasion steady-state abundance of the invading species. Then we conduct colonization experiments for commensal gut bacteria species Enterococcus faecium and Akkermansia muciniphila in hundreds of human stool-derived in vitro microbial communities, confirming that the data-driven approaches can predict the colonization outcomes in experiments. Furthermore, we find that while most resident species are predicted to have a weak negative impact on the colonization of exogenous species, strongly interacting species could significantly alter the colonization outcomes, e.g., Enterococcus faecalis inhibits the invasion of E. faecium invasion. The presented results suggest that the data-driven approaches are powerful tools to inform the ecology and management of microbial communities.

Structured community transitions explain the switching capacity of microbial systems

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

Richness and density jointly determine context dependence in bacterial interactions

Pairwise interactions are often used to predict features of complex microbial communities due to the challenge of measuring multi-species interactions in high dimensional contexts. This assumes that interactions are unaffected by community context. Here, we used synthetic bacterial communities to investigate that assumption by observing how interactions varied across contexts. Interactions were most often weakly negative and showed a phylogenetic signal among genera. Community richness and total density emerged as strong predictors of interaction strength and contributed to an attenuation of interactions as richness increased. Population level and per-capita measures of interactions both displayed such attenuation, suggesting factors beyond systematic changes in population size were involved; namely, changes to the interactions themselves. Nevertheless, pairwise interactions retained some explanatory power across contexts, provided those contexts were not substantially divergent in richness. These results suggest that understanding the emergent properties of microbial interactions can improve our ability to predict the features of microbial communities.

Stable States of a Microbial Community Are Formed by Dynamic Metabolic Networks with Members Functioning to Achieve Both Robustness and Plasticity

A more detailed understanding of the mechanisms underlying the formation of microbial communities is essential for the efficient management of microbial ecosystems. The stable states of microbial communities are commonly perceived as static and, thus, have not been extensively examined. The present study investigated stabilizing mechanisms, minority functions, and the reliability of quantitative ana-lyses, emphasizing a metabolic network perspective. A bacterial community, formed by batch transferred cultures supplied with phenol as the sole carbon and energy source and paddy soil as the inoculum, was analyzed using a principal coordinate ana-lysis (PCoA), mathematical models, and quantitative parameters defined as growth activity, community-changing activity, community-forming activity, vulnerable force, and resilience force depending on changes in the abundance of operational taxonomic units (OTUs) using 16S rRNA gene amplicon sequences. PCoA showed succession states until the 3rd transferred cultures and stable states from the 5th to 10th transferred cultures. Quantitative parameters indicated that the bacterial community was dynamic irrespective of the succession and stable states. Three activities fluctuated under stable states. Vulnerable and resilience forces were detected under the succession and stable states, respectively. Mathematical models indicated the construction of metabolic networks, suggesting the stabilizing mechanism of the community structure. Thirteen OTUs coexisted during stable states, and were recognized as core OTUs consisting of majorities, middle-class, and minorities. The abundance of the middle-class changed, whereas that of the others did not, which indicated that core OTUs maintained metabolic networks. Some extremely low abundance OTUs were consistently exchanged, suggesting a role for scavengers. These results indicate that stable states were formed by dynamic metabolic networks with members functioning to achieve robustness and plasticity.

Searching for Principles of Microbial Ecology Across Levels of Biological Organization

Microbial communities play pivotal roles in ecosystems across different scales, from global elemental cycles to household food fermentations. These complex assemblies comprise hundreds or thousands of microbial species whose abundances vary over time and space. Unraveling the principles that guide their dynamics at different levels of biological organization, from individual species, their interactions, to complex microbial communities, is a major challenge. To what extent are these different levels of organization governed by separate principles, and how can we connect these levels to develop predictive models for the dynamics and function of microbial communities? Here, we will discuss recent advances that point towards principles of microbial communities, rooted in various disciplines from physics, biochemistry, and dynamical systems. By considering the marine carbon cycle as a concrete example, we demonstrate how the integration of levels of biological organization can offer deeper insights into the impact of increasing temperatures, such as those associated with climate change, on ecosystem-scale processes. We argue that by focusing on principles that transcend specific microbiomes, we can pave the way for a comprehensive understanding of microbial community dynamics and the development of predictive models for diverse ecosystems.

Master regulators of biological systems in higher dimensions

A longstanding goal of biology is to identify the key genes and species that critically impact evolution, ecology, and health. Network analysis has revealed keystone species that regulate ecosystems and master regulators that regulate cellular genetic networks. Yet these studies have focused on pairwise biological interactions, which can be affected by the context of genetic background and other species present, generating higher-order interactions. The important regulators of higher-order interactions are unstudied. To address this, we applied a high-dimensional geometry approach that quantifies epistasis in a fitness landscape to ask how individual genes and species influence the interactions in the rest of the biological network. We then generated and also reanalyzed 5-dimensional datasets (two genetic, two microbiome). We identified key genes (e.g., the rbs locus and pykF) and species (e.g., Lactobacilli) that control the interactions of many other genes and species. These higher-order master regulators can induce or suppress evolutionary and ecological diversification by controlling the topography of the fitness landscape. Thus, we provide a method and mathematical justification for exploration of biological networks in higher dimensions.

Statistical prediction of microbial metabolic traits from genomes

The metabolic activity of microbial communities is central to their role in biogeochemical cycles, human health, and biotechnology. Despite the abundance of sequencing data characterizing these consortia, it remains a serious challenge to predict microbial metabolic traits from sequencing data alone. Here we culture 96 bacterial isolates individually and assay their ability to grow on 10 distinct compounds as a sole carbon source. Using these data as well as two existing datasets, we show that statistical approaches can accurately predict bacterial carbon utilization traits from genomes. First, we show that classifiers trained on gene content can accurately predict bacterial carbon utilization phenotypes by encoding phylogenetic information. These models substantially outperform predictions made by constraint-based metabolic models automatically constructed from genomes. This result solidifies our current knowledge about the strong connection between phylogeny and metabolic traits. However, phylogeny-based predictions fail to predict traits for taxa that are phylogenetically distant from any strains in the training set. To overcome this we train improved models on gene presence/absence to predict carbon utilization traits from gene content. We show that models that predict carbon utilization traits from gene presence/absence can generalize to taxa that are phylogenetically distant from the training set either by exploiting biochemical information for feature selection or by having sufficiently large datasets. In the latter case, we provide evidence that a statistical approach can identify putatively mechanistic genes involved in metabolic traits. Our study demonstrates the potential power for predicting microbial phenotypes from genotypes using statistical approaches.

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.

A rugged yet easily navigable fitness landscape

Fitness landscape theory predicts that rugged landscapes with multiple peaks impair Darwinian evolution, but experimental evidence is limited. In this study, we used genome editing to map the fitness of >260,000 genotypes of the key metabolic enzyme dihydrofolate reductase in the presence of the antibiotic trimethoprim, which targets this enzyme. The resulting landscape is highly rugged and harbors 514 fitness peaks. However, its highest peaks are accessible to evolving populations via abundant fitness-increasing paths. Different peaks share large basins of attraction that render the outcome of adaptive evolution highly contingent on chance events. Our work shows that ruggedness need not be an obstacle to Darwinian evolution but can reduce its predictability. If true in general, the complexity of optimization problems on realistic landscapes may require reappraisal.

Statistically learning the functional landscape of microbial communities

Microbial consortia exhibit complex functional properties in contexts ranging from soils to bioreactors to human hosts. Understanding how community composition determines function is a major goal of microbial ecology. Here we address this challenge using the concept of community-function landscapes-analogues to fitness landscapes-that capture how changes in community composition alter collective function. Using datasets that represent a broad set of community functions, from production/degradation of specific compounds to biomass generation, we show that statistically inferred landscapes quantitatively predict community functions from knowledge of species presence or absence. Crucially, community-function landscapes allow prediction without explicit knowledge of abundance dynamics or interactions between species and can be accurately trained using measurements from a small subset of all possible community compositions. The success of our approach arises from the fact that empirical community-function landscapes appear to be not rugged, meaning that they largely lack high-order epistatic contributions that would be difficult to fit with limited data. Finally, we show that this observation holds across a wide class of ecological models, suggesting community-function landscapes can be efficiently inferred across a broad range of ecological regimes. Our results open the door to the rational design of consortia without detailed knowledge of abundance dynamics or interactions.

Microbial inoculants with higher capacity to colonize soils improved wheat drought tolerance

Microbial inoculants have gained increasing attention worldwide as an eco-friendly solution for improving agriculture productivity. Several studies have demonstrated their potential benefits, such as enhanced resistance to drought, salinity, and pathogens. However, the beneficial impacts of inoculants remain inconsistent. This variability is attributed to limited knowledge of the mechanisms by which microbial inoculants affect crop growth and a lack of ecological characteristics of these inoculants that limit our ability to predict their beneficial effects. The first important step is believed to be the evaluation of the inoculant's ability to colonize new habitats (soils and plant roots), which could provide crops with beneficial functions and improve the consistency and efficiency of the inoculants. In this study, we aimed to investigate the impact of three microbial inoculants (two bacterial: P1 and P2, and one fungal: P3) on the growth and stress responses of three wheat varieties in two different soil types under drought conditions. Furthermore, we investigated the impact of microbial inoculants on soil microbial communities. Plant biomass and traits were measured, and high-throughput sequencing was used to characterize bulk and rhizosphere soil microbiomes after exposure to drought stress. Under drought conditions, plant shoot weight significantly increased (11.37%) under P1 treatments compared to uninoculated controls. In addition, total nitrogen enzyme activity increased significantly under P1 in sandy soil but not in clay soil. Importantly, network analyses revealed that P1, consisting of Bacillus paralicheniformis and Bacillus subtilis, emerged as the keystone taxa in sandy soil. Conversely, P2 and P3 failed to establish as keystone taxa, which may explain their insignificant impact on wheat performance under drought conditions. In conclusion, our study emphasizes the importance of effective colonization by microbial inoculants in promoting crop growth under drought conditions. Our findings support the development of microbial inoculants that robustly colonize plant roots for improved agricultural productivity.

Deconstructing taxa x taxa xenvironment interactions in the microbiota: A theoretical examination

A major objective of microbial ecology is to identify how the composition of microbial taxa shapes host phenotypes. However, most studies focus on pairwise interactions and ignore the potentially significant effects of higher-order microbial interactions.Here, we quantify the effects of higher-order interactions among taxa on host infection risk. We apply our approach to an in silico dataset that is built to resemble a population of insect hosts with gut-associated microbial communities at risk of infection from an intestinal parasite across a breadth of nutrient environmental contexts.We find that the effect of higher-order interactions is considerable and can change appreciably across environmental contexts. Furthermore, we show that higher-order interactions can stabilize community structure thereby reducing host susceptibility to parasite invasion.Our approach illustrates how incorporating the effects of higher-order interactions among gut microbiota across environments can be essential for understanding their effects on host phenotypes.

First-order transition to oscillation death in coupled oscillators with higher-order interactions

We investigate the dynamical evolution of Stuart-Landau oscillators globally coupled through conjugate or dissimilar variables on simplicial complexes. We report a first-order explosive phase transition from an oscillatory state to oscillation death, with higher-order (2-simplex triadic) interactions, as opposed to the second-order transition with only pairwise (1-simplex) interactions. Moreover, the system displays four distinct homogeneous steady states in the presence of triadic interactions, in contrast to the two homogeneous steady states observed with dyadic interactions. We calculate the backward transition point analytically, confirming the numerical results and providing the origin of the dynamical states in the transition region. The results are robust against the application of noise. The study will be useful in understanding complex systems, such as ecological and epidemiological, having higher-order interactions and coupling through conjugate variables.

Predictability of the community-function landscape in wine yeast ecosystems

Predictively linking taxonomic composition and quantitative ecosystem functions is a major aspiration in microbial ecology, which must be resolved if we wish to engineer microbial consortia. Here, we have addressed this open question for an ecological function of major biotechnological relevance: alcoholic fermentation in wine yeast communities. By exhaustively phenotyping an extensive collection of naturally occurring wine yeast strains, we find that most ecologically and industrially relevant traits exhibit phylogenetic signal, allowing functional traits in wine yeast communities to be predicted from taxonomy. Furthermore, we demonstrate that the quantitative contributions of individual wine yeast strains to the function of complex communities followed simple quantitative rules. These regularities can be integrated to quantitatively predict the function of newly assembled consortia. Besides addressing theoretical questions in functional ecology, our results and methodologies can provide a blueprint for rationally managing microbial processes of biotechnological relevance.

More Than the Sum of Its Parts: Unlocking the Power of Network Structure for Understanding Organization and Function in Microbiomes

Plant and soil microbiomes are integral to the health and productivity of plants and ecosystems, yet researchers struggle to identify microbiome characteristics important for providing beneficial outcomes. Network analysis offers a shift in analytical framework beyond "who is present" to the organization or patterns of coexistence between microbes within the microbiome. Because microbial phenotypes are often significantly impacted by coexisting populations, patterns of coexistence within microbiomes are likely to be especially important in predicting functional outcomes. Here, we provide an overview of the how and why of network analysis in microbiome research, highlighting the ways in which network analyses have provided novel insights into microbiome organization and functional capacities, the diverse network roles of different microbial populations, and the eco-evolutionary dynamics of plant and soil microbiomes.

Integrating a tailored recurrent neural network with Bayesian experimental design to optimize microbial community functions

Microbiomes interact dynamically with their environment to perform exploitable functions such as production of valuable metabolites and degradation of toxic metabolites for a wide range of applications in human health, agriculture, and environmental cleanup. Developing computational models to predict the key bacterial species and environmental factors to build and optimize such functions are crucial to accelerate microbial community engineering. However, there is an unknown web of interactions that determine the highly complex and dynamic behavior of these systems, which precludes the development of models based on known mechanisms. By contrast, entirely data-driven machine learning models can produce physically unrealistic predictions and often require significant amounts of experimental data to learn system behavior. We develop a physically-constrained recurrent neural network that preserves model flexibility but is constrained to produce physically consistent predictions and show that it can outperform existing machine learning methods in the prediction of certain experimentally measured species abundance and metabolite concentrations. Further, we present a closed-loop, Bayesian experimental design algorithm to guide data collection by selecting experimental conditions that simultaneously maximize information gain and target microbial community functions. Using a bioreactor case study, we demonstrate how the proposed framework can be used to efficiently navigate a large design space to identify optimal operating conditions. The proposed methodology offers a flexible machine learning approach specifically tailored to optimize microbiome target functions through the sequential design of informative experiments that seek to explore and exploit community functions.

Top-down identification of keystone taxa in the microbiome

Keystone taxa in ecological communities are native taxa that play an especially important role in the stability of their ecosystem. However, we still lack an effective framework for identifying these taxa from the available high-throughput sequencing without the notoriously difficult step of reconstructing the detailed network of inter-specific interactions. In addition, while most microbial interaction models assume pair-wise relationships, it is yet unclear whether pair-wise interactions dominate the system, or whether higher-order interactions are relevant. Here we propose a top-down identification framework, which detects keystones by their total influence on the rest of the taxa. Our method does not assume a priori knowledge of pairwise interactions or any specific underlying dynamics and is appropriate to both perturbation experiments and metagenomic cross-sectional surveys. When applied to real high-throughput sequencing of the human gastrointestinal microbiome, we detect a set of candidate keystones and find that they are often part of a keystone module - multiple candidate keystone species with correlated occurrence. The keystone analysis of single-time-point cross-sectional data is also later verified by the evaluation of two-time-points longitudinal sampling. Our framework represents a necessary advancement towards the reliable identification of these key players of complex, real-world microbial communities.

Strain dropouts reveal interactions that govern the metabolic output of the gut microbiome

The gut microbiome is complex, raising questions about the role of individual strains in the community. Here, we address this question by constructing variants of a complex defined community in which we eliminate strains that occupy the bile acid 7α-dehydroxylation niche. Omitting Clostridium scindens (Cs) and Clostridium hylemonae (Ch) eliminates secondary bile acid production and reshapes the community in a highly specific manner: eight strains change in relative abundance by >100-fold. In single-strain dropout communities, Cs and Ch reach the same relative abundance and dehydroxylate bile acids to a similar extent. However, Clostridium sporogenes increases >1,000-fold in the ΔCs but not ΔCh dropout, reshaping the pool of microbiome-derived phenylalanine metabolites. Thus, strains that are functionally redundant within a niche can have widely varying impacts outside the niche, and a strain swap can ripple through the community in an unpredictable manner, resulting in a large impact on an unrelated community-level phenotype.

A three-species synthetic community model whose rapid response to antagonism allows the study of higher-order dynamics and emergent properties in minutes

Microbial communities can be considered complex adaptive systems. Understanding how these systems arise from different components and how the dynamics of microbial interactions allow for species coexistence are fundamental questions in ecology. To address these questions, we built a three-species synthetic community, called BARS (Bacillota A + S + R). Each species in this community exhibits one of three ecological roles: Antagonistic, Sensitive, or Resistant, assigned in the context of a sediment community. We show that the BARS community reproduces features of complex communities and exhibits higher-order interaction (HOI) dynamics. In paired interactions, the majority of the S species (Sutcliffiella horikoshii 20a) population dies within 5 min when paired with the A species (Bacillus pumilus 145). However, an emergent property appears upon adding the third interactor, as antagonism of species A over S is not observed in the presence of the R species (Bacillus cereus 111). For the paired interaction, within the first 5 min, the surviving population of the S species acquires tolerance to species A, and species A ceases antagonism. This qualitative change reflects endogenous dynamics leading to the expression for tolerance to an antagonistic substance. The stability reached in the triple interaction exhibits a nonlinear response, highly sensitive to the density of the R species. In summary, our HOI model allows the study of the assembly dynamics of a three-species community and evaluating the immediate outcome within a 30 min frame. The BARS has features of a complex system where the paired interactions do not predict the community dynamics. The model is amenable to mechanistic dissection and to modeling how the parts integrate to achieve collective properties.

Global epistasis on fitness landscapes

Epistatic interactions between mutations add substantial complexity to adaptive landscapes and are often thought of as detrimental to our ability to predict evolution. Yet, patterns of global epistasis, in which the fitness effect of a mutation is well-predicted by the fitness of its genetic background, may actually be of help in our efforts to reconstruct fitness landscapes and infer adaptive trajectories. Microscopic interactions between mutations, or inherent nonlinearities in the fitness landscape, may cause global epistasis patterns to emerge. In this brief review, we provide a succinct overview of recent work about global epistasis, with an emphasis on building intuition about why it is often observed. To this end, we reconcile simple geometric reasoning with recent mathematical analyses, using these to explain why different mutations in an empirical landscape may exhibit different global epistasis patterns—ranging from diminishing to increasing returns. Finally, we highlight open questions and research directions.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Construction and Modeling of a Coculture Microplate for Real-Time Measurement of Microbial Interactions

The dynamic structures of microbial communities emerge from the complex network of interactions between their constituent microorganisms. Quantitative measurements of these interactions are important for understanding and engineering ecosystem structure. Here, we present the development and application of the BioMe plate, a redesigned microplate device in which pairs of wells are separated by porous membranes. BioMe facilitates the measurement of dynamic microbial interactions and integrates easily with standard laboratory equipment. We first applied BioMe to recapitulate recently characterized, natural symbiotic interactions between bacteria isolated from the Drosophila melanogaster gut microbiome. Specifically, the BioMe plate allowed us to observe the benefit provided by two Lactobacillus strains to an Acetobacter strain. We next explored the use of BioMe to gain quantitative insight into the engineered obligate syntrophic interaction between a pair of Escherichia coli amino acid auxotrophs. We integrated experimental observations with a mechanistic computational model to quantify key parameters associated with this syntrophic interaction, including metabolite secretion and diffusion rates. This model also allowed us to explain the slow growth observed for auxotrophs growing in adjacent wells by demonstrating that, under the relevant range of parameters, local exchange between auxotrophs is essential for efficient growth. The BioMe plate provides a scalable and flexible approach for the study of dynamic microbial interactions. IMPORTANCE Microbial communities participate in many essential processes from biogeochemical cycles to the maintenance of human health. The structure and functions of these communities are dynamic properties that depend on poorly understood interactions among different species. Unraveling these interactions is therefore a crucial step toward understanding natural microbiota and engineering artificial ones. Microbial interactions have been difficult to measure directly, largely due to limitations of existing methods to disentangle the contribution of different organisms in mixed cocultures. To overcome these limitations, we developed the BioMe plate, a custom microplate-based device that enables direct measurement of microbial interactions, by detecting the abundance of segregated populations of microbes that can exchange small molecules through a membrane. We demonstrated the possible application of the BioMe plate for studying both natural and artificial consortia. BioMe is a scalable and accessible platform that can be used to broadly characterize microbial interactions mediated by diffusible molecules.

Data-driven prediction of colonization outcomes for complex microbial communities

Complex microbial interactions can lead to different colonization outcomes of exogenous species, be they pathogenic or beneficial in nature. Predicting the colonization of exogenous species in complex communities remains a fundamental challenge in microbial ecology, mainly due to our limited knowledge of the diverse physical, biochemical, and ecological processes governing microbial dynamics. Here, we proposed a data-driven approach independent of any dynamics model to predict colonization outcomes of exogenous species from the baseline compositions of microbial communities. We systematically validated this approach using synthetic data, finding that machine learning models (including Random Forest and neural ODE) can predict not only the binary colonization outcome but also the post-invasion steady-state abundance of the invading species. Then we conducted colonization experiments for two commensal gut bacteria species Enterococcus faecium and Akkermansia muciniphila in hundreds of human stool-derived in vitro microbial communities, confirming that the data-driven approach can successfully predict the colonization outcomes. Furthermore, we found that while most resident species were predicted to have a weak negative impact on the colonization of exogenous species, strongly interacting species could significantly alter the colonization outcomes, e.g., the presence of Enterococcus faecalis inhibits the invasion of E. faecium . The presented results suggest that the data-driven approach is a powerful tool to inform the ecology and management of complex microbial communities.

Conceptualizing microbe-plasmid communities as complex adaptive systems

Plasmids shape microbial communities' diversity, structure, and function. Nevertheless, we lack a mechanistic understanding of how community structure and dynamics emerge from local microbe-plasmid interactions and coevolution. Addressing this gap is challenging because multiple processes operate simultaneously at multiple levels of organization. For example, immunity operates between a plasmid and a cell, but incompatibility mechanisms regulate coexistence between plasmids. Conceptualizing microbe-plasmid communities as complex adaptive systems is a promising approach to overcoming these challenges. I illustrate how agent-based evolutionary modeling, extended by network analysis, can be used to quantify the relative importance of local processes governing community dynamics. These theoretical developments can advance our understanding of plasmid ecology and evolution, especially when combined with empirical data.

The community-function landscape of microbial consortia

Quantitatively linking the composition and function of microbial communities is a major aspiration of microbial ecology. Microbial community functions emerge from a complex web of molecular interactions between cells, which give rise to population-level interactions among strains and species. Incorporating this complexity into predictive models is highly challenging. Inspired by a similar problem in genetics of predicting quantitative phenotypes from genotypes, an ecological community-function (or structure-function) landscape could be defined that maps community composition and function. In this piece, we present an overview of our current understanding of these community landscapes, their uses, limitations, and open questions. We argue that exploiting the parallels between both landscapes could bring powerful predictive methodologies from evolution and genetics into ecology, providing a boost to our ability to engineer and optimize microbial consortia.

Engineering plant microbiomes by integrating eco-evolutionary principles into current strategies

Engineering plant microbiomes has the potential to improve plant health in a rapid and sustainable way. Rapidly changing climates and relatively long timelines for plant breeding make microbiome engineering an appealing approach to improving food security. However, approaches that have shown promise in the lab have not resulted in wide-scale implementation in the field. Here, we suggest the use of an integrated approach, combining mechanistic molecular and genetic knowledge, with ecological and evolutionary theory, to target knowledge gaps in plant microbiome engineering that may facilitate translatability of approaches into the field. We highlight examples where understanding microbial community ecology is essential for a holistic understanding of the efficacy and consequences of microbiome engineering. We also review examples where understanding plant-microbe evolution could facilitate the design of plants able to recruit specific microbial communities. Finally, we discuss possible trade-offs in plant-microbiome interactions that should be considered during microbiome engineering efforts so as not to introduce off-target negative effects. We include classic and emergent approaches, ranging from microbial inoculants to plant breeding to host-driven microbiome engineering, and address areas that would benefit from multidisciplinary approaches.

Ecological landscapes guide the assembly of optimal microbial communities

Assembling optimal microbial communities is key for various applications in biofuel production, agriculture, and human health. Finding the optimal community is challenging because the number of possible communities grows exponentially with the number of species, and so an exhaustive search cannot be performed even for a dozen species. A heuristic search that improves community function by adding or removing one species at a time is more practical, but it is unknown whether this strategy can discover an optimal or nearly optimal community. Using consumer-resource models with and without cross-feeding, we investigate how the efficacy of search depends on the distribution of resources, niche overlap, cross-feeding, and other aspects of community ecology. We show that search efficacy is determined by the ruggedness of the appropriately-defined ecological landscape. We identify specific ruggedness measures that are both predictive of search performance and robust to noise and low sampling density. The feasibility of our approach is demonstrated using experimental data from a soil microbial community. Overall, our results establish the conditions necessary for the success of the heuristic search and provide concrete design principles for building high-performing microbial consortia.

Higher-order interactions shape microbial interactions as microbial community complexity increases

Non-pairwise interactions, or higher-order interactions (HOIs), in microbial communities have been described as significant drivers of emergent features in microbiomes. Yet, the re-organization of microbial interactions between pairwise cultures and larger communities remains largely unexplored from a molecular perspective but is central to our understanding and further manipulation of microbial communities. Here, we used a bottom-up approach to investigate microbial interaction mechanisms from pairwise cultures up to 4-species communities from a simple microbiome (Hafnia alvei, Geotrichum candidum, Pencillium camemberti and Escherichia coli). Specifically, we characterized the interaction landscape for each species combination involving E. coli by identifying E. coli's interaction-associated mutants using an RB-TnSeq-based interaction assay. We observed a deep reorganization of the interaction-associated mutants, with very few 2-species interactions conserved all the way up to a 4-species community and the emergence of multiple HOIs. We further used a quantitative genetics strategy to decipher how 2-species interactions were quantitatively conserved in higher community compositions. Epistasis-based analysis revealed that, of the interactions that are conserved at all levels of complexity, 82% follow an additive pattern. Altogether, we demonstrate the complex architecture of microbial interactions even within a simple microbiome, and provide a mechanistic and molecular explanation of HOIs.

A Transcriptomic Analysis of Higher-Order Ecological Interactions in a Eukaryotic Model Microbial Ecosystem

Nonlinear ecological interactions within microbial ecosystems and their contribution to ecosystem functioning remain largely unexplored. Higher-order interactions, or interactions in systems comprised of more than two members that cannot be explained by cumulative pairwise interactions, are particularly understudied, especially in eukaryotic microorganisms. The wine fermentation ecosystem presents an ideal model to study yeast ecosystem establishment and functioning. Some pairwise ecological interactions between wine yeast species have been characterized, but very little is known about how more complex, multispecies systems function. Here, we evaluated nonlinear ecosystem properties by determining the transcriptomic response of Saccharomyces cerevisiae to pairwise versus tri-species culture. The transcriptome revealed that genes expressed during pairwise coculture were enriched in the tri-species data set but also that just under half of the data set comprised unique genes attributed to a higher-order response. Through interactive protein-association network visualizations, a holistic cell-wide view of the gene expression data was generated, which highlighted known stress response and metabolic adaptation mechanisms which were specifically activated during tri-species growth. Further, extracellular metabolite data corroborated that the observed differences were a result of a biotic stress response. This provides exciting new evidence showing the presence of higher-order interactions within a model microbial ecosystem. IMPORTANCE Higher-order interactions are one of the major blind spots in our understanding of microbial ecosystems. These systems remain largely unpredictable and are characterized by nonlinear dynamics, in particular when the system is comprised of more than two entities. By evaluating the transcriptomic response of S. cerevisiae to an increase in culture complexity from a single species to two- and three-species systems, we were able to confirm the presence of a unique response in the more complex setting that could not be explained by the responses observed at the pairwise level. This is the first data set that provides molecular targets for further analysis to explain unpredictable ecosystem dynamics in yeast.

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.

Engineering microbial consortia with rationally designed cellular interactions

Synthetic microbial consortia represent a frontier of synthetic biology that promises versatile engineering of cellular functions. They are primarily developed through the design and construction of cellular interactions that coordinate individual dynamics and generate collective behaviors. Here we review recent advances in the engineering of synthetic communities through cellular-interaction programming. We first examine fundamental building blocks for intercellular communication and unidirectional positive and negative interactions. We then recap the assembly of the building blocks for creating bidirectional interactions in two-species ecosystems, which is followed by the discussion of engineering toward complex communities with increasing species numbers, under spatial contexts, and via model-guided design. We conclude by summarizing major challenges and future opportunities of engineered microbial ecosystems.

Ecological modelling approaches for predicting emergent properties in microbial communities

Recent studies have brought forward the critical role of emergent properties in shaping microbial communities and the ecosystems of which they are a part. Emergent properties-patterns or functions that cannot be deduced linearly from the properties of the constituent parts-underlie important ecological characteristics such as resilience, niche expansion and spatial self-organization. While it is clear that emergent properties are a consequence of interactions within the community, their non-linear nature makes mathematical modelling imperative for establishing the quantitative link between community structure and function. As the need for conservation and rational modulation of microbial ecosystems is increasingly apparent, so is the consideration of the benefits and limitations of the approaches to model emergent properties. Here we review ecosystem modelling approaches from the viewpoint of emergent properties. We consider the scope, advantages and limitations of Lotka-Volterra, consumer-resource, trait-based, individual-based and genome-scale metabolic models. Future efforts in this research area would benefit from capitalizing on the complementarity between these approaches towards enabling rational modulation of complex microbial ecosystems.

Recurrent neural networks enable design of multifunctional synthetic human gut microbiome dynamics

Predicting the dynamics and functions of microbiomes constructed from the bottom-up is a key challenge in exploiting them to our benefit. Current models based on ecological theory fail to capture complex community behaviors due to higher order interactions, do not scale well with increasing complexity and in considering multiple functions. We develop and apply a long short-term memory (LSTM) framework to advance our understanding of community assembly and health-relevant metabolite production using a synthetic human gut community. A mainstay of recurrent neural networks, the LSTM learns a high dimensional data-driven non-linear dynamical system model. We show that the LSTM model can outperform the widely used generalized Lotka-Volterra model based on ecological theory. We build methods to decipher microbe-microbe and microbe-metabolite interactions from an otherwise black-box model. These methods highlight that Actinobacteria, Firmicutes and Proteobacteria are significant drivers of metabolite production whereas Bacteroides shape community dynamics. We use the LSTM model to navigate a large multidimensional functional landscape to design communities with unique health-relevant metabolite profiles and temporal behaviors. In sum, the accuracy of the LSTM model can be exploited for experimental planning and to guide the design of synthetic microbiomes with target dynamic functions.

Competition for nutrients increases invasion resistance during assembly of microbial communities

The assembly of microbial communities through sequential invasions of microbial species is challenging to study experimentally. Here, I used genome‐scale metabolic models of multiple species to model community assembly. Each such model represents all known biochemical reactions that a species uses to build biomass from nutrients in the environment. Species interactions in such models emerge from first biochemical principles, either through competition for environmental nutrients, or through cross‐feeding on metabolic by‐products excreted by resident species. I used these models to study 250 community assembly sequences. In each such sequence, a community changes through successive species invasions. During the 250 assembly sequences, communities become more species‐rich and invasion‐resistant. Resistance against both constructive and destructive invasions – those that entail species extinction – is associated with high community productivity, high biomass, and low concentrations of unused carbon. Competition for nutrients outweighs the influence of cross‐feeding on the growth rate of individual species. In a community assembly network of all communities that arise during the 250 assembly sequences, some communities occur more often than expected by chance. These include invasion resistant “attractor” communities with high biomass that arise late in community assembly and persist preferentially because of their invasion resistance. Genome‐scale metabolic models can reveal generic properties of microbial communities that are independent of the resident species and the environment.

Dynamic proteome allocation regulates the profile of interaction of auxotrophic bacterial consortia

Microbial ecosystems are composed of multiple species in constant metabolic exchange. A pervasive interaction in microbial communities is metabolic cross-feeding and occurs when the metabolic burden of producing costly metabolites is distributed between community members, in some cases for the benefit of all interacting partners. In particular, amino acid auxotrophies generate obligate metabolic inter-dependencies in mixed populations and have been shown to produce a dynamic profile of interaction that depends upon nutrient availability. However, identifying the key components that determine the pair-wise interaction profile remains a challenging problem, partly because metabolic exchange has consequences on multiple levels, from allocating proteomic resources at a cellular level to modulating the structure, function and stability of microbial communities. To evaluate how ppGpp-mediated resource allocation drives the population-level profile of interaction, here we postulate a multi-scale mathematical model that incorporates dynamics of proteome partition into a population dynamics model. We compare our computational results with experimental data obtained from co-cultures of auxotrophic Escherichia coli K12 strains under a range of amino acid concentrations and population structures. We conclude by arguing that the stringent response promotes cooperation by inhibiting the growth of fast-growing strains and promoting the synthesis of metabolites essential for other community members.

Linking Reactive Oxygen Species (ROS) to Abiotic and Biotic Feedbacks in Plant Microbiomes: The Dose Makes the Poison

In nature, plants develop in complex, adaptive environments. Plants must therefore respond efficiently to environmental stressors to maintain homeostasis and enhance their fitness. Although many coordinated processes remain integral for achieving homeostasis and driving plant development, reactive oxygen species (ROS) function as critical, fast-acting orchestrators that link abiotic and biotic responses to plant homeostasis and development. In addition to the suite of enzymatic and non-enzymatic ROS processing pathways that plants possess, they also rely on their microbiota to buffer and maintain the oxidative window needed to balance anabolic and catabolic processes. Strong evidence has been communicated recently that links ROS regulation to the aggregated function(s) of commensal microbiota and plant-growth-promoting microbes. To date, many reports have put forth insightful syntheses that either detail ROS regulation across plant development (independent of plant microbiota) or examine abiotic–biotic feedbacks in plant microbiomes (independent of clear emphases on ROS regulation). Here we provide a novel synthesis that incorporates recent findings regarding ROS and plant development in the context of both microbiota regulation and plant-associated microbes. Specifically, we discuss various roles of ROS across plant development to strengthen the links between plant microbiome functioning and ROS regulation for both basic and applied research aims.

Evaluation of Existing Methods for High-Order Epistasis Detection

Finding epistatic interactions among loci when expressing a phenotype is a widely employed strategy to understand the genetic architecture of complex traits in GWAS. The abundance of methods dedicated to the same purpose, however, makes it increasingly difficult for scientists to decide which method is more suitable for their studies. This work compares the different epistasis detection methods published during the last decade in terms of runtime, detection power and type I error rate, with a special emphasis on high-order interactions. Results show that in terms of detection power, the only methods that perform well across all experiments are the exhaustive methods, although their computational cost may be prohibitive in large-scale studies. Regarding non-exhaustive methods, not one could consistently find epistasis interactions when marginal effects are absent. If marginal effects are present, there are methods that perform well for high-order interactions, such as BADTrees, FDHE-IW, SingleMI or SNPHarvester. As for false-positive control, only SNPHarvester, FDHE-IW and DCHE show good results. The study concludes that there is no single epistasis detection method to recommend in all scenarios. Authors should prioritize exhaustive methods when sufficient computational resources are available considering the data set size, and resort to non-exhaustive methods when the analysis time is prohibitive.

An ensemble approach to the structure-function problem in microbial communities

The metabolic activity of microbial communities plays a primary role in the flow of essential nutrients throughout the biosphere. Molecular genetics has revealed the metabolic pathways that model organisms utilize to generate energy and biomass, but we understand little about how the metabolism of diverse, natural communities emerges from the collective action of its constituents. We propose that quantifying and mapping metabolic fluxes to sequencing measurements of genomic, taxonomic, or transcriptional variation across an ensemble of diverse communities, either in the laboratory or in the wild, can reveal low-dimensional descriptions of community structure that can explain or predict their emergent metabolic activity. We survey the types of communities for which this approach might be best suited, review the analytical techniques available for quantifying metabolite fluxes in communities, and discuss what types of data analysis approaches might be lucrative for learning the structure-function mapping in communities from these data.

Exploring the roles of microbes in facilitating plant adaptation to climate change

Plants benefit from their close association with soil microbes which assist in their response to abiotic and biotic stressors. Yet much of what we know about plant stress responses is based on studies where the microbial partners were uncontrolled and unknown. Under climate change, the soil microbial community will also be sensitive to and respond to abiotic and biotic stressors. Thus, facilitating plant adaptation to climate change will require a systems-based approach that accounts for the multi-dimensional nature of plant-microbe-environment interactions. In this perspective, we highlight some of the key factors influencing plant-microbe interactions under stress as well as new tools to facilitate the controlled study of their molecular complexity, such as fabricated ecosystems and synthetic communities. When paired with genomic and biochemical methods, these tools provide researchers with more precision, reproducibility, and manipulability for exploring plant-microbe-environment interactions under a changing climate.

Genomic structure predicts metabolite dynamics in microbial communities

The metabolic activities of microbial communities play a defining role in the evolution and persistence of life on Earth, driving redox reactions that give rise to global biogeochemical cycles. Community metabolism emerges from a hierarchy of processes, including gene expression, ecological interactions, and environmental factors. In wild communities, gene content is correlated with environmental context, but predicting metabolite dynamics from genomes remains elusive. Here, we show, for the process of denitrification, that metabolite dynamics of a community are predictable from the genes each member of the community possesses. A simple linear regression reveals a sparse and generalizable mapping from gene content to metabolite dynamics for genomically diverse bacteria. A consumer-resource model correctly predicts community metabolite dynamics from single-strain phenotypes. Our results demonstrate that the conserved impacts of metabolic genes can predict community metabolite dynamics, enabling the prediction of metabolite dynamics from metagenomes, designing denitrifying communities, and discovering how genome evolution impacts metabolism.