Higher-Order Interaction between Species Inhibits Bacterial Invasion of a Phototroph-Predator Microbial Community

Developing the 'omic toolkit of comparative physiologists

Typical 'omic analyses reduce complex biological systems to simple lists of supposedly independent variables, failing to account for changes in the wider transcriptional landscape. In this commentary, we discuss the utility of network approaches for incorporating this wider context into the study of physiological phenomena. We highlight opportunities to build on traditional network tools by utilising cutting-edge techniques to account for higher order interactions (i.e. beyond pairwise associations) within datasets, allowing for more accurate models of complex 'omic systems. Finally, we show examples of previous works utilising network approaches to gain additional insight into their organisms of interest. As 'omics grow in both their popularity and breadth of application, so does the requirement for flexible analytical tools capable of interpreting and synthesising complex datasets.

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.

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.

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.

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.

Manipulating the physical distance between cells during soil colonization reveals the importance of biotic interactions in microbial community assembly

BACKGROUND

Microbial communities are of tremendous importance for ecosystem functioning and yet we know little about the ecological processes driving the assembly of these communities in the environment. Here, we used an unprecedented experimental approach based on the manipulation of physical distance between neighboring cells during soil colonization to determine the role of bacterial interactions in soil community assembly. We hypothesized that experimentally manipulating the physical distance between bacterial cells will modify the interaction strengths leading to differences in microbial community composition, with increasing distance between neighbors favoring poor competitors.

RESULTS

We found significant differences in both bacterial community diversity, composition and co-occurrence networks after soil colonization that were related to physical distancing. We show that reducing distances between cells resulted in a loss of bacterial diversity, with at least 41% of the dominant OTUs being significantly affected by physical distancing. Our results suggest that physical distancing may differentially modulate competitiveness between neighboring species depending on the taxa present in the community. The mixing of communities that assembled at high and low cell densities did not reveal any "home field advantage" during coalescence. This confirms that the observed differences in competitiveness were due to biotic rather than abiotic filtering.

CONCLUSIONS

Our study demonstrates that the competitiveness of bacteria strongly depends on cell density and community membership, therefore highlighting the fundamental role of microbial interactions in the assembly of soil communities.

Learning beyond-pairwise interactions enables the bottom-up prediction of microbial community structure

Understanding the assembly of multispecies microbial communities represents a significant challenge in ecology and has wide applications in agriculture, wastewater treatment, and human healthcare domains. Traditionally, studies on the microbial community assembly focused on analyzing pairwise relationships among species; however, neglecting higher-order interactions, i.e., the change of pairwise relationships in the community context, may lead to substantial deviation from reality. Herein, we have proposed a simple framework that incorporates higher-order interactions into a bottom-up prediction of the microbial community assembly and examined its accuracy using a seven-member synthetic bacterial community on a host plant, duckweed. Although the synthetic community exhibited emergent properties that cannot be predicted from pairwise coculturing results, our results demonstrated that incorporating information from three-member combinations allows the acceptable prediction of the community structure and actual interaction forces within it. This reflects that the occurrence of higher-order effects follows consistent patterns, which can be predicted even from trio combinations, the smallest unit of higher-order interactions. These results highlight the possibility of predicting, explaining, and understanding the microbial community structure from the bottom-up by learning interspecies interactions from simple beyond-pairwise combinations.

Programmed microalgae-gel promotes chronic wound healing in diabetes

Chronic diabetic wounds are at lifelong risk of developing diabetic foot ulcers owing to severe hypoxia, excessive reactive oxygen species (ROS), a complex inflammatory microenvironment, and the potential for bacterial infection. Here we develop a programmed treatment strategy employing live Haematococcus (HEA). By modulating light intensity, HEA can be programmed to perform a variety of functions, such as antibacterial activity, oxygen supply, ROS scavenging, and immune regulation, suggesting its potential for use in programmed therapy. Under high light intensity (658 nm, 0.5 W/cm2), green HEA (GHEA) with efficient photothermal conversion mediate wound surface disinfection. By decreasing the light intensity (658 nm, 0.1 W/cm2), the photosynthetic system of GHEA can continuously produce oxygen, effectively resolving the problems of hypoxia and promoting vascular regeneration. Continuous light irradiation induces astaxanthin (AST) accumulation in HEA cells, resulting in a gradual transformation from a green to red hue (RHEA). RHEA effectively scavenges excess ROS, enhances the expression of intracellular antioxidant enzymes, and directs polarization to M2 macrophages by secreting AST vesicles via exosomes. The living HEA hydrogel can sterilize and enhance cell proliferation and migration and promote neoangiogenesis, which could improve infected diabetic wound healing in female mice.

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.

Designing host-associated microbiomes using the consumer/resource model

A key step towards rational microbiome engineering is the in silico sampling of realistic microbial communities that correspond to desired host phenotypes, and vice versa. This remains challenging due to a lack of generative models that simultaneously model compositions of host-associated microbiomes and host phenotypes. To that end, we present a machine learning model based on the consumer/resource (C/R) framework. In the model, variation in microbial ecosystem composition arises due to differences in the availability of effective resources (latent variables) while species' resource preferences remain conserved. Variation in the same latent variables is used to model phenotypic variation across hosts. In silico microbiomes generated by our model accurately reproduce universal and dataset-specific statistics of bacterial communities. The model allows us to address two salient questions in microbiome design: (1) which host phenotypes maximally constrain the composition of the host-associated microbiome? and (2) what are plausible microbiome compositions corresponding to user-specified host phenotypes? Thus, our model aids the design and analysis of microbial communities associated with host phenotypes of interest.

Synthetic model ecosystem of 12 cryopreservable microbial species allowing for a noninvasive approach

Simultaneous understanding of both population and ecosystem dynamics is crucial in an era marked by the degradation of ecosystem services. Experimental ecosystems are a powerful tool for understanding these dynamics; however, they often face technical challenges, typically falling into two categories: "complex but with limited replicability microcosms" and "highly replicable but overly simplistic microcosms." Herein, we present a high-throughput synthetic microcosm system comprising 12 functionally and phylogenetically diverse microbial species. These species are axenically culturable, cryopreservable, and can be measured noninvasively via microscopy, aided by machine learning. This system includes prokaryotic and eukaryotic producers and decomposers, and eukaryotic consumers to ensure functional redundancy. Our model system exhibited key features of a complex ecosystem: (i) various positive and negative interspecific interactions, (ii) higher-order interactions beyond two-species dynamics, (iii) probabilistic dynamics leading to divergent outcomes, and (iv) stable nonlinear transitions. We identified several conditions under which at least one species from each of the three functional groups-producers, consumers, and decomposers-and one functionally redundant species, persisted for over six months. These conditions set the stage for detailed investigations in the future. Given its designability and experimental replicability, our model ecosystem offers a promising platform for deeper insights integrating both population and ecosystem dynamics.

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.

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.

Warming during embryogenesis induces a lasting transcriptomic signature in fishes

Exposure to elevated temperatures during embryogenesis can influence the plasticity of tissues in later life. Despite these long-term changes in plasticity, few differentially expressed genes are ever identified, suggesting that the developmental programming of later life plasticity may occur through the modulation of other aspects of transcriptomic architecture, such as gene network organisation. Here, we use network modelling approaches to demonstrate that warm temperatures during embryonic development (developmental warming) have consistent effects in later life on the organisation of transcriptomic networks across four diverse species of fishes: Scyliorhinus canicula, Danio rerio, Dicentrarchus labrax, and Gasterosteus aculeatus. The transcriptomes of developmentally warmed fishes are characterised by an increased entropy of their pairwise gene interaction networks, implying a less structured, more 'random' set of gene interactions. We also show that, in zebrafish subject to developmental warming, the entropy of an individual gene within a network is associated with that gene's probability of expression change during temperature acclimation in later life. However, this association is absent in animals reared under 'control' conditions. Thus, the thermal environment experienced during embryogenesis can alter transcriptomic organisation in later life, and these changes may influence an individual's responsiveness to future temperature challenges.

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.

Modulating Biological Rhythms: A Noncomputational Strategy Harnessing Nonlinearity and Decoupling Frequency and Amplitude

Understanding and achieving concurrent modulation of amplitude and frequency, particularly adjusting one quantity and simultaneously sustaining the other at an invariant level, are of paramount importance for complex biophysical systems, including the signal pathway where different frequency indicates different upstream signal yielding a certain downstream physiological function while different amplitude further determines different efficacy of a downstream output. However, such modulators with clearly described and universally valid mechanisms are still lacking. Here, we rigorously propose an easy-to-use control strategy containing only one or two steps, leveraging the nonlinearity in the modulated systems to decouple frequency and amplitude in a noncomputational manner. The strategy's efficacy is demonstrated using representative biochemical systems and, thus, it could be potentially applicable to modulating rhythms in experiments of biochemistry and synthetic biology.

Connecting higher-order interactions with ecological stability in experimental aquatic food webs

Community ecology is built on theories that represent the strength of interactions between species as pairwise links. Higher-order interactions (HOIs) occur when a species changes the pairwise interaction between a focal pair. Recent theoretical work has highlighted the stabilizing role of HOIs for large, simulated communities, yet it remains unclear how important higher-order effects are in real communities. Here, we used experimental communities of aquatic protists to examine the relationship between HOIs and stability (as measured by the persistence of a species in a community). We cultured a focal pair of consumers in the presence of additional competitors and a predator and collected time series data of their abundances. We then fitted competition models with and without HOIs to measure interaction strength between the focal pair across different community compositions. We used survival analysis to measure the persistence of individual species. We found evidence that additional species positively affected persistence of the focal species and that HOIs were present in most of our communities. However, persistence was only linked to HOIs for one of the focal species. Our results vindicate community ecology theory positing that species interactions may deviate from assumptions of pairwise interactions, opening avenues to consider possible consequences for coexistence and stability.

Algae drive convergent bacterial community assembly at low dilution frequency

Microbial community assembly is a complex dynamical process that determines community structure and function. The interdependence of inter-species interactions and nutrient availability presents a challenge for understanding community assembly. We sought to understand how external nutrient supply rate modulated interactions to affect the assembly process. A statistical decomposition of taxonomic structures of bacterial communities assembled with and without algae and at varying dilution frequencies allowed the separation of the effects of biotic (presence of algae) and abiotic (dilution frequency) factors on community assembly. For infrequent dilutions, the algae strongly impact community assembly, driving initially diverse bacterial consortia to converge to a common structure. Analyzing sequencing data revealed that this convergence is largely mediated by a decline in the relative abundance of specific taxa in the presence of algae. This study shows that complex phototroph-heterotroph communities can be powerful model systems for understanding assembly processes relevant to the global ecosystem functioning.

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.

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.

Experimental community coalescence sheds light on microbial interactions in soil and restores impaired functions

BACKGROUND

Microbes typically live in communities where individuals can interact with each other in numerous ways. However, knowledge on the importance of these interactions is limited and derives mainly from studies using a limited number of species grown in coculture. Here, we manipulated soil microbial communities to assess the contribution of interactions between microorganisms for assembly of the soil microbiome.

RESULTS

By combining experimental removal (taxa depletion in the community) and coalescence (mixing of manipulated and control communities) approaches, we demonstrated that interactions between microorganisms can play a key role in determining their fitness during soil recolonization. The coalescence approach not only revealed the importance of density-dependent interactions in microbial community assembly but also allowed to restore partly or fully community diversity and soil functions. Microbial community manipulation resulted in shifts in both inorganic nitrogen pools and soil pH, which were related to the proportion of ammonia-oxidizing bacteria.

CONCLUSIONS

Our work provides new insights into the understanding of the importance of microbial interactions in soil. Our top-down approach combining removal and coalescence manipulation also allowed linking community structure and ecosystem functions. Furthermore, these results highlight the potential of manipulating microbial communities for the restoration of soil ecosystems. Video Abstract.

Microbial dynamics and bioreactor performance are interlinked with organic matter removal from wastewater treatment plant effluent

Optimizing bioreactor performance for organic matter removal can achieve sustainable and energy-efficient micropollutant removal in subsequent tertiary treatment. Bioreactor performance heavily depends on its resident microbial community; hence, a deeper understanding of community dynamics is essential. The microbial communities of three different bioreactors (biological activated carbon, moving bed biofilm reactor, sand filter), used for organic matter removal from wastewater treatment effluent, were characterized by 16S rRNA gene amplicon sequence analysis. An interdependency between bioreactor performance and microbial community profile was observed. Overall, Proteobacteria was the most predominant phylum, and Comamonadaceae was the most predominant family in all bioreactors. The relative abundance of the genus Roseococcus was positively correlated with organic matter removal. A generalized Lotka-Volterra (gLV) model was established to understand the interactions in the microbial community. By identifying microbial dynamics and their role in bioreactors, a strategy can be developed to improve bioreactor performance.

Controlling the human microbiome

We coexist with a vast number of microbes that live in and on our bodies. Those microbes and their genes are collectively known as the human microbiome, which plays important roles in human physiology and diseases. We have acquired extensive knowledge of the organismal compositions and metabolic functions of the human microbiome. However, the ultimate proof of our understanding of the human microbiome is reflected in our ability to manipulate it for health benefits. To facilitate the rational design of microbiome-based therapies, there are many fundamental questions to be addressed at the systems level. Indeed, we need a deep understanding of the ecological dynamics associated with such a complex ecosystem before we rationally design control strategies. In light of this, this review discusses progress from various fields, e.g., community ecology, network science, and control theory, that are helping us make progress toward the ultimate goal of controlling the human microbiome.

Chlamydomonas reinhardtii, a Reference Organism to Study Algal-Microbial Interactions: Why Can't They Be Friends?

The stability and harmony of ecological niches rely on intricate interactions between their members. During evolution, organisms have developed the ability to thrive in different environments, taking advantage of each other. Among these organisms, microalgae are a highly diverse and widely distributed group of major primary producers whose interactions with other organisms play essential roles in their habitats. Understanding the basis of these interactions is crucial to control and exploit these communities for ecological and biotechnological applications. The green microalga Chlamydomonas reinhardtii, a well-established model, is emerging as a model organism for studying a wide variety of microbial interactions with ecological and economic significance. In this review, we unite and discuss current knowledge that points to C. reinhardtii as a model organism for studying microbial interactions.

Effect of phenol formaldehyde-associated microplastics on soil microbial community, assembly, and functioning

Increasing investigations explore the effects of plastic pollutants on bacterial communities, diversity, and functioning in various ecosystems. However, the impact of microplastics (MPs) on the eukaryotic community, microbial assemblages, and interactions is still limited. Here, we investigated bacterial and micro-eukaryotic communities and functioning in soils with different concentrations of phenol formaldehyde-associated MPs (PF-MPs), and revealed the factors, such as soil properties, microbial community assembly, and interactions between microbes, influencing them. Our results showed that a high concentration (1%) of PF-MPs decreased the microbial interactions and the contribution of deterministic processes to the community assembly of microbes, and consequently changed the communities of bacteria, but not eukaryotes. A significant and negative relationship was determined between N₂O emission rate and functional genes related to nitrification, indicating that the competitive interactions between functional microbes would affect the nitrogen cycling of soil ecosystem. We further found that vegetable biomass weakly decreased in treatments with a higher concentration of PF-MPs and positively related to the diversity of micro-eukaryotic communities and functional diversity of bacterial communities. These results suggest that a high concentration of the PF-MPs would influence crop growth by changing microbial communities, interactions, and eukaryotic and functional diversity. Our findings provide important evidence for agriculture management of phenol formaldehyde and suggest that we must consider their threats to microbial community compositions, diversity, and assemblage in soils due to the accumulation of PF-MPs widely used in the field.

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.

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.

Resource availability drives bacteria community resistance to pathogen invasion via altering bacterial pairwise interactions

Microbial interactions within resident communities are a major determinant of resistance to pathogen invasion. Yet, interactions vary with environmental conditions, raising the question of how community composition and environments interactively shape invasion resistance. Here, we use resource availability (RA) as a model parameter altering the resistance of model bacterial communities to invasion by the plant pathogenic bacterium Ralstonia solanacearum. We found that at high RA, interactions between resident bacterial species were mainly driven by the direct antagonism, in terms of the means of invader inhibition. Consequently, the competitive resident communities with a higher production of antibacterial were invaded to a lesser degree than facilitative communities. At low RA, bacteria produced little direct antagonist potential, but facilitative communities reached a relatively higher community productivity, which showed higher resistance to pathogen invasion than competitive communities with lower productivities. This framework may lay the basis to understand complex microbial interactions and biological invasion as modulated by the dynamic changes of environmental resource availability.

The microbial communities (bacteria, algae, zooplankton, and fungi) improved biofloc technology including the nitrogen-related material cycle in Litopenaeus vannamei farms

Microbes are essential in biofloc technology for controlling nitrogen levels in water. The composition and function of microorganisms with biofloc systems were reported; however, data on microorganisms other than bacteria, such as algae (which are essential in the nitrogen cycle) and zooplankton (which are bacterial and algal predators), remain limited. The microbial communities (including bacteria, algae, zooplankton, and fungi) were investigated in shrimp farms using biofloc technology. Using Illumina MiSeq sequencing, the V4 region of 18S rRNA and the V3-V4 region of 16S rRNA were utilized for the analysis of the eukaryotic and prokaryotic microbial communities. As a result, it was found that the biofloc in the shrimp farm consisted of 48.73%-73.04% eukaryotic organisms and 26.96%-51.27% prokaryotic organisms. In these shrimp farms, prokaryotic microbial communities had higher specie richness and diversity than eukaryotic microbial communities. However, the eukaryotic microbial communities were more abundant than their prokaryotic counterparts, while algae and zooplankton dominated them. It was discovered that the structures of the microbial communities in the shrimp farms seemed to depend on the effects of predation by zooplankton and other related organisms. The results provided the nitrogen cycle in biofloc systems by the algal and bacterial groups in microbial communities.

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.

Obtaining Bioproducts from the Studies of Signals and Interactions between Microalgae and Bacteria

The applications of microalgae biomass have been widely studied worldwide. The classical processes used in outdoor cultivations of microalgae, in closed or open photobioreactors, occur in the presence of bacteria. Understanding how communication between cells occurs through quorum sensing and evaluating co-cultures allows the production of microalgae and cyanobacteria to be positively impacted by bacteria, in order to guarantee safety and profitability in the production process. In addition, the definition of the effects that occur during an interaction, promotes insights to improve the production of biomolecules, and to develop innovative products. This review presents the interactions between microalgae and bacteria, including compounds exchanges and communication, and addresses the development of new pharmaceutical, cosmetic and food bioproducts from microalgae based on these evaluations, such as prebiotics, vegan skincare products, antimicrobial compounds, and culture media with animal free protein for producing vaccines and other biopharmaceutical products. The use of microalgae as raw biomass or in biotechnological platforms is in line with the fulfillment of the 2030 Agenda related to the Sustainable Development Goals (SDGs).

Disaggregation as an interaction mechanism among intestinal bacteria

The gut microbiome contains hundreds of interacting species that together influence host health and development. The mechanisms by which intestinal microbes can interact, however, remain poorly mapped and are often modeled as spatially unstructured competitions for chemical resources. Recent imaging studies examining the zebrafish gut have shown that patterns of aggregation are central to bacterial population dynamics. In this study, we focus on bacterial species of genera Aeromonas and Enterobacter. Two zebrafish gut-derived isolates, Aeromonas ZOR0001 (AE) and Enterobacter ZOR0014 (EN), when mono-associated with the host, are highly aggregated and located primarily in the intestinal midgut. An Aeromonas isolate derived from the commensal strain, Aeromonas-MB4 (AE-MB4), differs from the parental strain in that it is composed mostly of planktonic cells localized to the anterior gut. When challenged by AE-MB4, clusters of EN rapidly fragment into non-motile, slow-growing, dispersed individual cells with overall abundance two orders of magnitude lower than the mono-association value. In the presence of a certain set of additional gut bacterial species, these effects on EN are dampened. In particular, if AE-MB4 invades an already established multi-species community, EN persists in the form of large aggregates. These observations reveal an unanticipated competition mechanism based on manipulation of bacterial spatial organization, namely dissolution of aggregates, and provide evidence that multi-species communities may facilitate stable intestinal co-existence.

Effective coastal Escherichia coli monitoring by unmanned aerial vehicles (UAV) thermal infrared images

Coastal Escherichia coli (E. coli) significantly influence ocean safety and public health, thus requiring an effective E. coli pollution monitoring. However conventional detection relying on manual field sampling is time-consuming. Here, this study established an E. coli estimation model based on thermal remote sensing of unmanned aerial vehicles (UAV). This model was developed against one-year comprehensive field work in a representative sandy beach and further validated against 50 beaches in Hong Kong to evaluate its applicability. The estimated E. coli concentrations were in a reliable agreement with direct measurements. For this model, this study deployed the radon-222 (²²²Rn) as a bridging tracer to couple UAV thermal images and coastal E. coli concentrations. Coastal ²²²Rn can be reflected on the UAV thermal images, and there was a good positive correlation between the ²²²Rn activity and coastal E. coli concentration via one-year field data. Hence, coupling the ²²²Rn activity estimated from UAV thermal images and the relationship between ²²²Rn and E. coli, this study can readily monitor coastal E. coli by UAV. These findings highlighted that UAV technology is an effective approach to measure the E. coli concentrations and can further pave the way for an efficient coastal E. coli monitoring and public health risk warning.

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.

Detecting and interpreting higher-order interactions in ecological communities

When species simultaneously compete with two or more species of competitor, higher-order interactions (HOIs) can lead to emergent properties not present when species interact in isolated pairs. To extend ecological theory to multi-competitor communities, ecologists must confront the challenges of measuring and interpreting HOIs in models of competition fit to data from nature. Such efforts are hindered by the fact that different studies use different definitions, and these definitions have unclear relationships to one another. Here, we propose a distinction between 'soft' HOIs, which identify possible interaction modification by competitors, and 'hard' HOIs, which identify interactions uniquely emerging in systems with three or more competitors. We show how these two classes of HOI differ in their motivation and interpretation, as well as the tests one uses to identify them in models fit to data. We then show how to operationalise this structure of definitions by analysing the results of a simulated competition experiment underlain by a consumer resource model. In the course of doing so, we clarify the challenges of interpreting HOIs in nature, and suggest a more precise framing of this research endeavour to catalyse further investigations.

Synthetic nonlinear computation for genetic circuit design

Computation frameworks have been studied in synthetic biology to achieve biosignals integration and processing, for biosensing and therapeutics applications. Biological systems exhibit nonlinearity across scales from the molecular level, to biochemical network and intercellular systems. At the molecular level, cooperative bindings contribute to nonlinear molecular signal processing in a way similar to weight variables. At the intracellular network level, feedback and feedforward regulations result in cell behaviors such as multistability and adaptation. When biochemical networks are distributed in different cell groups, intercelluar networks can generate population dynamics. Here, we review works that highlight nonlinear computations in synthetic biology. We group the works according to the scale of implementations, from the cis-transcription level, to biochemical circuit level and cellular networks.

Protist Predation Influences the Temperature Response of Bacterial Communities

Temperature strongly influences microbial community structure and function, in turn contributing to global carbon cycling that can fuel further warming. Recent studies suggest that biotic interactions among microbes may play an important role in determining the temperature responses of these communities. However, how predation regulates these microbiomes under future climates is still poorly understood. Here, we assess whether predation by a key global bacterial consumer—protists—influences the temperature response of the community structure and function of a freshwater microbiome. To do so, we exposed microbial communities to two cosmopolitan protist species—Tetrahymena thermophila and Colpidium sp.—at two different temperatures, in a month-long microcosm experiment. While microbial biomass and respiration increased with temperature due to community shifts, these responses changed over time and in the presence of protists. Protists influenced microbial biomass and respiration rate through direct and indirect effects on bacterial community structure, and predator presence actually reduced microbial respiration at elevated temperature. Indicator species analyses showed that these predator effects were mostly determined by phylum-specific bacterial responses to protist density and cell size. Our study supports previous findings that temperature is an important driver of microbial communities but also demonstrates that the presence of a large predator can mediate these responses to warming.

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.

Limited Pairwise Synergistic and Antagonistic Interactions Impart Stability to Microbial Communities

One of the central goals of ecology is to explain and predict coexistence of species. In this context, microbial communities provide a model system where community structure can be studied in environmental niches and in laboratory conditions. A community of microbial population is stabilized by interactions between participating species. However, the nature of these stabilizing interactions has remained largely unknown. Theory and experiments have suggested that communities are stabilized by antagonistic interactions between member species, and destabilized by synergistic interactions. However, experiments have also revealed that a large fraction of all the interactions between species in a community are synergistic in nature. To understand the relative significance of the two types of interactions (synergistic vs. antagonistic) between species, we perform simulations of microbial communities with a small number of participating species using two frameworks—a replicator equation and a Lotka-Volterra framework. Our results demonstrate that synergistic interactions between species play a critical role in maintaining diversity in cultures. These interactions are critical for the ability of the communities to survive perturbations and maintain diversity. We follow up the simulations with quantification of the extent to which synergistic and antagonistic interactions are present in a bacterial community present in a soil sample. Overall, our results show that community stability is largely achieved with the help of synergistic interactions between participating species. However, we perform experiments to demonstrate that antagonistic interactions, in specific circumstances, can also contribute toward community stability.

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.

Interactions of E. coli with algae and aquatic vegetation in natural waters

Both algae and bacteria are essential inhabitants of surface waters. Their presence is of ecological significance and sometimes of public health concern triggering various control actions. Interactions of microalgae, macroalgae, submerged aquatic vegetation, and bacteria appear to be important phenomena necessitating a deeper understanding by those involved in research and management of microbial water quality. Given the long-standing reliance on Escherichia coli as an indicator of the potential presence of pathogens in natural waters, understanding its biology in aquatic systems is necessary. The major effects of algae and aquatic vegetation on E. coli growth and survival, including changes in the nutrient supply, modification of water properties and constituents, impact on sunlight radiation penetration, survival as related to substrate attachment, algal mediation of secondary habitats, and survival inhibition due to the release of toxic substances and antibiotics, are discussed in this review. An examination of horizontal gene transfer and antibiotic resistance potential, strain-specific interactions, effects on the microbial, microalgae, and grazer community structure, and hydrodynamic controls is given. Outlooks due to existing and expected consequences of climate change and advances in observation technologies via high-resolution satellite imaging, unmanned aerial vehicles (drones), and mathematical modeling are additionally covered. The multiplicity of interactions among bacteria, algae, and aquatic vegetation as well as multifaceted impacts of these interactions, create a wide spectrum of research opportunities and technology developments.

High-order correlations in species interactions lead to complex diversity-stability relationships for ecosystems

How ecosystems maintain stability is an active area of research. Inspired by applications of random matrix theory in nuclear physics, May showed decades ago that in an ecosystem model with many randomly interacting species, increasing species diversity decreases the stability of the ecosystem. There have since been many additions to May's efforts, one being an improved understanding the effect of mutualistic, competitive, or predator-prey-like correlations between pairs of species. Here we extend a random matrix technique developed in the context of spin-glass theory to study the effect of high-order correlations among species interactions. The resulting analytically solvable models include next-to-nearest-neighbor correlations in the species interaction network, such as the enemy of my enemy is my friend, as well as higher-order correlations. We find qualitative differences from May and others' models, including nonmonotonic diversity-stability relationships. Furthermore, inclusion of particular next-to-nearest-neighbor correlations in predator-prey as opposed to mutualist-competitive networks causes the former to transition to being more stable at higher species diversity. We discuss potential applicability of our results to microbiota engineering and to the ecology of interpredator interactions, such as cub predation between lions and hyenas as well as companionship between humans and dogs.

Scaling down the microbial loop: data-driven modelling of growth interactions in a diatom-bacterium co-culture

An intricate set of interactions characterizes marine ecosystems. One of the most important is represented by the microbial loop, which includes the exchange of dissolved organic matter (DOM) from phototrophic organisms to heterotrophic bacteria. Here, it can be used as the major carbon and energy source. This interaction is one of the foundations of the entire ocean food-web. The carbon fixed by phytoplankton can be redirected to bacteria in two main ways; either (i) bacteria feed on dead phytoplankton cells or (ii) DOM is actively released by phytoplankton (a process resulting in up to 50% of the fixed carbon leaving the cell). Here, we have set up a co-culture of the diatom Phaeodactylum tricornutum and the chemoheterotrophic bacterium Pseudoalteromonas haloplanktis TAC125 and used this system to study the interactions between these two representatives of the microbial loop. We show that the bacterium can thrive on diatom-derived carbon and that this growth can be sustained by both diatom dead cells and diatom-released compounds. These observations were formalized in a network of putative interactions between P. tricornutum and P. haloplanktis and implemented in a model that reproduces the observed co-culture dynamics, revealing an overall accuracy of our hypotheses in explaining the experimental data.

Closed microbial communities self-organize to persistently cycle carbon

Life on Earth depends on ecologically driven nutrient cycles to regenerate resources. Understanding how nutrient cycles emerge from a complex web of ecological processes is a central challenge in ecology. However, we lack model ecosystems that can be replicated, manipulated, and quantified in the laboratory, making it challenging to determine how changes in composition and the environment impact cycling. Enabled by a new high-precision method to quantify carbon cycling, we show that materially closed microbial ecosystems (CES) provided with only light self-organize to robustly cycle carbon. Studying replicate CES that support a carbon cycle reveals variable community composition but a conserved set of metabolic capabilities. Our study helps establish CES as model biospheres for studying how ecosystems persistently cycle nutrients.

Cycles of nutrients (N, P, etc.) and resources (C) are a defining emergent feature of ecosystems. Cycling plays a critical role in determining ecosystem structure at all scales, from microbial communities to the entire biosphere. Stable cycles are essential for ecosystem persistence because they allow resources and nutrients to be regenerated. Therefore, a central problem in ecology is understanding how ecosystems are organized to sustain robust cycles. Addressing this problem quantitatively has proved challenging because of the difficulties associated with manipulating ecosystem structure while measuring cycling. We address this problem using closed microbial ecosystems (CES), hermetically sealed microbial consortia provided with only light. We develop a technique for quantifying carbon cycling in hermetically sealed microbial communities and show that CES composed of an alga and diverse bacterial consortia self-organize to robustly cycle carbon for months. Comparing replicates of diverse CES, we find that carbon cycling does not depend strongly on the taxonomy of the bacteria present. Moreover, despite strong taxonomic differences, self-organized CES exhibit a conserved set of metabolic capabilities. Therefore, an emergent carbon cycle enforces metabolic but not taxonomic constraints on ecosystem organization. Our study helps establish closed microbial communities as model ecosystems to study emergent function and persistence in replicate systems while controlling community composition and the environment.

Open challenges for microbial network construction and analysis

Microbial network construction is a popular explorative data analysis technique in microbiome research. Although a large number of microbial network construction tools has been developed to date, there are several issues concerning the construction and interpretation of microbial networks that have received less attention. The purpose of this perspective is to draw attention to these underexplored challenges of microbial network construction and analysis.

Addressing context dependence in ecology

Context dependence is widely invoked to explain disparate results in ecology. It arises when the magnitude or sign of a relationship varies due to the conditions under which it is observed. Such variation, especially when unexplained, can lead to spurious or seemingly contradictory conclusions, which can limit understanding and our ability to transfer findings across studies, space, and time. Using examples from biological invasions, we identify two types of context dependence resulting from four sources: mechanistic context dependence arises from interaction effects; and apparent context dependence can arise from the presence of confounding factors, problems of statistical inference, and methodological differences among studies. Addressing context dependence is a critical challenge in ecology, essential for increased understanding and prediction.