Metabolic resource allocation in individual microbes determines ecosystem interactions and spatial dynamics

Moving from genome-scale to community-scale metabolic models for the human gut microbiome

Metabolic models of individual microorganisms or small microbial consortia have become standard research tools in the bioengineering and systems biology fields. However, extending metabolic modelling to diverse microbial communities, such as those in the human gut, remains a practical challenge from both modelling and experimental validation perspectives. In complex communities, metabolic models accounting for community dynamics, or those that consider multiple objectives, may provide optimal predictions over simpler steady-state models, but require a much higher computational cost. Here we describe some of the strengths and limitations of microbial community-scale metabolic models and argue for a robust validation framework for developing personalized, mechanistic and accurate predictions of microbial community metabolic behaviours across environmental contexts. Ultimately, quantitatively accurate microbial community-scale metabolic models could aid in the design and testing of personalized prebiotic, probiotic and dietary interventions that optimize for translationally relevant outcomes.

Replicating community dynamics reveals how initial composition shapes the functional outcomes of bacterial communities

Bacterial communities play key roles in global biogeochemical cycles, industry, agriculture, human health, and animal husbandry. There is therefore great interest in understanding bacterial community dynamics so that they can be controlled and engineered to optimise ecosystem services. We assess the reproducibility and predictability of bacterial community dynamics by creating a frozen archive of hundreds of naturally-occurring bacterial communities that we repeatedly revive and track in a standardised, complex resource environment. Replicate communities follow reproducible trajectories and the community dynamics closely map to ecosystem functioning. However, even under standardised conditions, the communities exhibit tipping-points, where small differences in initial community composition create divergent compositional and functional outcomes. The predictability of community trajectories therefore requires detailed knowledge of rugged compositional landscapes where ecosystem properties are not the inevitable result of prevailing environmental conditions but can be tilted toward different outcomes depending on the initial community composition. Our results shed light on the relationship between composition and function, opening new avenues to understand the feasibility and limitations of function prediction in complex microbial communities.

Cross-feeding of amino acid pathway intermediates is common in co-cultures of auxotrophic Escherichia coli

Amino acid auxotrophy refers to an organism's inability to synthesize one or more amino acids that are required for cell growth. In microbiome research, co-cultures of amino acid auxotrophs are often used to investigate metabolite cross-feeding interactions and model community dynamics. Thus far, it has been implicitly assumed that amino acids are mainly cross-fed between these auxotrophs. However, this assumption has not been fully verified. For example, it could be that intermediates of amino acid biosynthesis pathways are exchanged instead, or in addition to amino acids. If true, this would significantly increase the complexity of metabolic interactions that needs to be considered. Here, we show that metabolic pathway intermediates are indeed exchanged in many co-cultures of amino acid auxotrophs. To demonstrate this, we selected 25 E. coli single gene knockouts that are auxotrophic for five different amino acids: arginine, histidine, isoleucine, proline, and tryptophan. In co-culture experiments, we paired strains that shared the same amino acid auxotrophy and monitored cell growth. We observed growth in 23 out of 55 strain pairings, indicating that pathway intermediates were exchanged between the strains. To provide further support for cross-feeding of pathway intermediates, auxotrophic E. coli strains were cultured in media supplemented with commercially available metabolic pathway intermediates at different concentrations. Supplementing media with these metabolites recovered cell growth as was predicted from the co-culture experiments. Most of these metabolites supported high growth rates, even when present at low concentrations (10 μM), suggesting the presence of high affinity transporters for these metabolites. In total, we identified eight metabolic pathway intermediates that were likely exchanged between the auxotrophic E. coli strains and verified six of these, including histidinol, N-acetyl-L-ornithine, L-ornithine, L-citrulline, keto-isoleucine and anthranilate. Taken together, this work demonstrates that exchange of metabolic pathway intermediates is more common than has been assumed so far. In future, these exchanges must be explicitly considered when constructing models of metabolite cross-feeding interactions in microbial communities and when interpreting results from microbiome studies involving auxotrophic strains.

Integrated use of electrochemical anaerobic reactors and genomic based modeling for characterizing methanogenic activity in microbial communities exposed to BTEX contamination

In soil polluted with benzene, toluene, ethylbenzene, and xylenes (BTEX), oxygen is rapidly depleted by aerobic respiration, creating a redox gradient across the plume. Under anaerobic conditions, BTEX biodegradation is then coupled with fermentation and methanogenesis. This study aimed to characterize this multi-step process, focusing on the interactions and functional roles of key microbial groups involved. A reactor system, comprising an Anaerobic Bioreactor (AB) and two Microbial Electrolysis Cell (MEC) chambers, designed to represent different spatial zones along the redox gradient, operated for 160 days with intermittent exposure to BTEX. The functional differentiation of each chamber was reflected by the gas emission profiles: 50%, 12% and 84% methane in the AB, anode and cathode chambers, respectively. The taxonomic profiling, assessed using 16S amplicon sequencing, led to the identification chamber-characteristic taxonomic groups. To translate the taxonomic shift into a functional shift, community dynamics was transformed into a simulative platform based on genome scale metabolic models constructed for 21 species that capture both key functionalities and taxonomies. Representatives include BTEX degraders, fermenters, iron reducers acetoclastic and hydrogenotrophic methanogens. Functionality was inferred according to the identification of the functional gene bamA as a biomarker for anaerobic BTEX degradation, taxonomy and literature support. Comparison of the predicted performances of the reactor-specific communities confirmed that the simulation successfully captured the experimentally recorded functional variation. Variations in the predicted exchange profiles between chambers capture reported and novel competitive and cooperative interactions between methanogens and non-methanogens. Examples include the exchange profiles of hypoxanthine (HYXN) and acetate between fermenters and methanogens, suggesting mechanisms underlying the supportive/repressive effect of taxonomic divergence on methanogenesis. Hence, the platform represents a pioneering attempt to capture the full spectrum of community activity in methanogenic hydrocarbon biodegradation while supporting the future design of optimization strategies.

The role of bacterial metabolism in human gut colonization

Can we anticipate the emergence of the next pandemic antibiotic-resistant bacterial clone? Addressing such an ambitious question relies on our ability to comprehensively understand the ecological and epidemiological factors fostering the evolution of high-risk clones. Among these factors, the ability to persistently colonize and thrive in the human gut is crucial for most high-risk clones. Nonetheless, the causes and mechanisms facilitating successful gut colonization remain obscure. Here, we review recent evidence that suggests that bacterial metabolism plays a pivotal role in determining the ability of high-risk clones to colonize the human gut. Subsequently, we outline novel approaches that enable the exploration of microbial metabolism at an unprecedented scale and level of detail. A thorough understanding of the constraints and opportunities of bacterial metabolism in gut colonization will foster our ability to predict the emergence of high-risk clones and take appropriate containment strategies.

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

Harnessing microbial division of labor for biomanufacturing: a review of laboratory and formal modeling approaches

Bioprocess industries aim to meet the increasing demand for product complexity by designing enhanced cellular and metabolic capabilities for the host. Monocultures, standard biomanufacturing workhorses, are often restricted in their capability to meet these demands, and the solution often involves the genetic modification of the host. Synthetic microbial communities are a promising alternative to monocultures because they exhibit division of labor, enabling efficient resource utilization and pathway modularity. This specialization minimizes metabolic burden and enhances robustness to perturbations, providing a competitive advantage. Despite this potential, their utilization in biotechnological or bioprocessing applications remains limited. The recent emergence of new and innovative community design tools and strategies, particularly those harnessing the division of labor, holds promise to change this outlook. Understanding the microbial interactions governing natural microbial communities can be used to identify complementary partners, informing synthetic community design. Therefore, we particularly consider engineering division of labor in synthetic microbial communities as a viable solution to accelerate progress in the field. This review presents the current understanding of how microbial interactions enable division of labor and how this information can be used to design synthetic microbial communities to perform tasks otherwise unfeasible to individual organisms. We then evaluate laboratory and formal modeling approaches specifically developed to: elucidate microbial community physiology, guide experimental design, and improve our understanding of complex community interactions assisting synthetic community design. By synthesizing these insights, we aim to present a comprehensive framework that advances the use of microbial communities in biomanufacturing applications.

Genome-scale models in human metabologenomics

Metabologenomics integrates metabolomics with other omics data types to comprehensively study the genetic and environmental factors that influence metabolism. These multi-omics data can be incorporated into genome-scale metabolic models (GEMs), which are highly curated knowledge bases that explicitly account for genes, transcripts, proteins and metabolites. By including all known biochemical reactions catalysed by enzymes and transporters encoded in the human genome, GEMs analyse and predict the behaviour of complex metabolic networks. Continued advancements to the scale and scope of GEMs - from cells and tissues to microbiomes and the whole body - have helped to design effective treatments and develop better diagnostic tools for metabolic diseases. Furthermore, increasing amounts of multi-omics data are incorporated into GEMs to better identify the underlying mechanisms, biomarkers and potential drug targets of metabolic diseases.

Lysis of Escherichia coli by colicin Ib contributes to bacterial cross-feeding by releasing active β-galactosidase

The diffusible toxin ColIb produced by Salmonella enterica serovar Typhimurium SL1344 is a potent inhibitor of Escherichia coli growth. To identify and parameterize metabolic cross-feeding in states of competition, we established defined communities in which E. coli was the only species able to access a sole carbon source, lactose. Although ColIb was predicted to undermine cross-feeding by killing the lactose-converting E. coli, S. enterica populations thrived in co-culture. We discovered that ColIb caused the release of active β-galactosidase from E. coli cells, which induced galactose uptake by S. enterica. Although iron limitation stimulates ColIb production and makes E. coli more sensitive to the toxin, ColIb killing in iron-limited conditions did not enhance iron acquisition or siderophore scavenging by S. enterica. Also unexpected was the rapid rate at which resistance to ColIb evolved in E. coli through spontaneous mutation of the ColIb receptor gene cirA or horizontal acquisition of the S. enterica colicin immunity gene imm. Mathematical modelling effectively predicted the growth kinetics of E. coli and S. enterica populations, revealing a tractable system in which ColIb can shrink a competitor population while simultaneously amplifying the metabolic contributions of the suppressed population.

Metabolic Flux Modeling in Marine Ecosystems

Ocean metabolism constitutes a complex, multiscale ensemble of biochemical reaction networks harbored within and between the boundaries of a myriad of organisms. Gaining a quantitative understanding of how these networks operate requires mathematical tools capable of solving in silico the resource allocation problem each cell faces in real life. Toward this goal, stoichiometric modeling of metabolism, such as flux balance analysis, has emerged as a powerful computational tool for unraveling the intricacies of metabolic processes in microbes, microbial communities, and multicellular organisms. Here, we provide an overview of this approach and its applications, future prospects, and practical considerations in the context of marine sciences. We explore how flux balance analysis has been employed to study marine organisms, help elucidate nutrient cycling, and predict metabolic capabilities within diverse marine environments, and highlight future prospects for this field in advancing our knowledge of marine ecosystems and their sustainability.

Azotobacter vinelandii N2 fixation increases in co-culture with the PGPR Bacillus subtilis in a nitrogen concentration-dependent manner

Biological nitrogen fixation (BNF) is an essential source of new nitrogen (N) for terrestrial ecosystems. The abiotic factors regulating BNF have been extensively studied in various ecosystems and laboratory settings. Despite this, our understanding of the impact of neighboring bacteria on N2 fixer activity remains limited. Here, we explored this question using a co-culture of the two model species: the free-living diazotroph Azotobacter vinelandii and the non-fixing plant growth-promoting rhizobacteria Bacillus subtilis. We observed that the interaction between the two bacteria was modulated by N availability. Under N-replete conditions, B. subtilis outcompeted A. vinelandii in the co-culture. Under N-limiting conditions, BNF activity by A. vinelandii was enhanced in the presence of B. subtilis. Reciprocally, the presence of A. vinelandii repressed sporulation by B. subtilis and supported its growth likely through N transfer. N inputs by A. vinelandii were doubled in the presence of B. subtilis compared to the monoculture, primarily due to the retention of a robust N2 fixation activity in the stationary phase. A proteomic analysis revealed that A. vinelandii N metabolism, particularly the molybdenum nitrogenase isoform protein levels (NifK and NifD), was upregulated during the stationary growth phase in the presence of B. subtilis. This study revealed that N stress drives bacterial interactions and activity in a two-species community, especially in the stationary phase.

IMPORTANCE

Reducing inputs of chemical N fertilizers is essential to develop a more sustainable agriculture. The stimulation of biological nitrogen fixation by N2 fixers in multispecies cultures, here the plant growth-promoting rhizobacteria Azotobacter vinelandii and Bacillus subtilis, opens opportunities for the formulation of biofertilizers consortia. While most research on N2 fixation historically focussed on the exponential growth phase of microorganisms, we observed that Bacillus subtilis stimulated Azotobacter vinelandii N2 fixation mostly during the stationary phase. This result highlights that more research on the factors controlling N2 fixation repression during the stationary growth phase, especially bacteria-bacteria interactions, is eagerly needed.

A systematic discussion and comparison of the construction methods of synthetic microbial community

Synthetic microbial community has widely concerned in the fields of agriculture, food and environment over the past few years. However, there is little consensus on the method to synthetic microbial community from construction to functional verification. Here, we review the concept, characteristics, history and applications of synthetic microbial community, summarizing several methods for synthetic microbial community construction, such as isolation culture, core microbiome mining, automated design, and gene editing. In addition, we also systematically summarized the design concepts, technological thresholds, and applicable scenarios of various construction methods, and highlighted their advantages and limitations. Ultimately, this review provides four efficient, detailed, easy-to-understand and -follow steps for synthetic microbial community construction, with major implications for agricultural practices, food production, and environmental governance.

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Exploring interspecific interaction variability in microbiota: A review

Interspecific interactions are an important component and a strong selective force in microbial communities. Over the past few decades, there has been a growing awareness of the variability in microbial interactions, and various studies are already unraveling the inner working dynamics in microbial communities. This has prompted scientists to develop novel techniques for characterizing the varying interspecific interactions among microbes. Here, we review the precise definitions of pairwise and high-order interactions, summarize the key concepts related to interaction variability, and discuss the strengths and weaknesses of emerging characterization techniques. Specifically, we found that most methods can accurately predict or provide direct information about microbial pairwise interactions. However, some of these methods inevitably mask the underlying high-order interactions in the microbial community. Making reasonable assumptions and choosing a characterization method to explore varying microbial interactions should allow us to better understand and engineer dynamic microbial systems.

Algorithm of spatial-temporal simulation for environment-strain interactions in strain-strain consortia based on resource competition mechanism

Interaction simulation for co-culture systems is important for optimizing culture conditions and improving yields. For industrial production, the environment significantly affects the spatial-temporal microbial interactions. However, the current research on polymicrobial interactions mainly focuses on interaction patterns among strains, and neglects the environment influence. Based on the resource competition relationship between two strains, this research set up the modules of cellular physicochemical properties, nutrient uptake and metabolite release, cellular survival, cell swimming and substrate diffusion, and investigated the spatial-temporal strain-environment interactions through module coupling and data mining. Furthermore, in an Escherichia coli-Saccharomyces cerevisiae consortium, the total net reproduction rate decreased as glucose was consumed. E. coli gradually dominated favorable positions due to its higher glucose utilization capacity, reaching 100 % abundance with a competitive strength of 0.86 for glucose. Conversely, S. cerevisiae decreased to 0 % abundance with a competitive strength of 0.14. The simulation results of environment influence on strain competitiveness showed that inoculation ratio and dissolved oxygen strongly influenced strain competitiveness. Specifically, strain competitiveness increased with higher inoculation ratio, whereas E. coli competitiveness increased as dissolved oxygen increased, in contrast to S. cerevisiae. On the other hand, substrate diffusion condition, micronutrients and toxins had minimal influence on strain competitiveness. This method offers a straightforward procedure without featured downscaling and provides novel insights into polymicrobial interaction simulation.

Assembly of functional microbial ecosystems: from molecular circuits to communities

Microbes compete and cooperate with each other via a variety of chemicals and circuits. Recently, to decipher, simulate, or reconstruct microbial communities, many researches have been engaged in engineering microbiomes with bottom-up synthetic biology approaches for diverse applications. However, they have been separately focused on individual perspectives including genetic circuits, communications tools, microbiome engineering, or promising applications. The strategies for coordinating microbial ecosystems based on different regulation circuits have not been systematically summarized, which calls for a more comprehensive framework for the assembly of microbial communities. In this review, we summarize diverse cross-talk and orthogonal regulation modules for de novo bottom-up assembling functional microbial ecosystems, thus promoting further consortia-based applications. First, we review the cross-talk communication-based regulations among various microbial communities from intra-species and inter-species aspects. Then, orthogonal regulations are summarized at metabolites, transcription, translation, and post-translation levels, respectively. Furthermore, to give more details for better design and optimize various microbial ecosystems, we propose a more comprehensive design-build-test-learn procedure including function specification, chassis selection, interaction design, system build, performance test, modeling analysis, and global optimization. Finally, current challenges and opportunities are discussed for the further development and application of microbial ecosystems.

Dynamic metabolic modeling of ATP allocation during viral infection

Viral pathogens, like SARS-CoV-2, hijack the host's macromolecular production machinery, imposing an energetic burden that is distributed across cellular metabolism. To explore the dynamic metabolic tension between the host's survival and viral replication, we developed a computational framework that uses genome-scale models to perform dynamic Flux Balance Analysis of human cell metabolism during virus infections. Relative to previous models, our framework addresses the physiology of viral infections of non-proliferating host cells through two new features. First, by incorporating the lipid content of SARS-CoV-2 biomass, we discovered activation of previously overlooked pathways giving rise to new predictions of possible drug targets. Furthermore, we introduce a dynamic model that simulates the partitioning of resources between the virus and the host cell, capturing the extent to which the competition depletes the human cells from essential ATP. By incorporating viral dynamics into our COMETS framework for spatio-temporal modeling of metabolism, we provide a mechanistic, dynamic and generalizable starting point for bridging systems biology modeling with viral pathogenesis. This framework could be extended to broadly incorporate phage dynamics in microbial systems and ecosystems.

Metabolic Trade-Offs Can Reverse the Resource-Diversity Relationship

AbstractFor species that partition resources, the classic expectation is that increasing resource diversity allows for increased species diversity. On the other hand, for neutral species, such as those competing equally for a single resource, diversity reflects a balance between the rate of introduction of novelty (e.g., by immigration or speciation) and the rate of extinction. Recent models of microbial metabolism have identified scenarios where metabolic trade-offs among species partitioning multiple resources can produce emergent neutral-like dynamics. In this hybrid scenario, one might expect that both resource diversity and immigration will act to boost species diversity. We show, however, that the reverse may be true: when metabolic trade-offs hold and population sizes are sufficiently large, increasing resource diversity can act to reduce species diversity, sometimes drastically. This reversal is explained by a generic transition between neutral- and niche-like dynamics, driven by the diversity of resources. The inverted resource-diversity relationship that results may be a signature of consumer-resource systems with strong metabolic trade-offs.

Effect of microbial interactions on performance of community metabolic modeling algorithms: flux balance analysis (FBA), community FBA (cFBA) and SteadyCom

To explore the impact of microbial interactions on outcomes from three prevalent algorithms (Flux Balance Analysis (FBA), community FBA (cFBA), and SteadyCom) analyzing microbial community metabolic networks, five toy community models representing common microbial interactions were designed. These include commensalism, mutualism, competition, mutualism-competition, and commensalism-competition. Various scenarios, considering different biomass yields and substrate constraints, were examined for each type. In commensal communities, all algorithms consistently produced similar results. However, changes in biomass yields and substrate constraints led to variable abundances (0.33-0.8) and community growth rates (2-5 1/h) within a broad range. For competitive communities, all algorithms predicted growth of fastest-growing member. To comply with the natural coexistence of members, suboptimal solutions over optimal point are recommended. FBA faced challenges in modeling mutualism, consistently predicting growth of only one member. Although cFBA and SteadyCom resulted in a lower community growth rate, coexistence of both members were satisfied. In toy models with dual interactions, more realistic outcomes were achieved contrary to purely competitive model as the dependency fosters the coexistence which was missing in the competitive only scenarios. These findings emphasize the importance of algorithm choice based on specific microbial interaction types for reliable community behavior predictions.​.

Game-Theoretic Flux Balance Analysis Model for Predicting Stable Community Composition

Models for microbial interactions attempt to understand and predict the steady state network of inter-species relationships in a community, e.g. competition for shared metabolites, and cooperation through cross-feeding. Flux balance analysis (FBA) is an approach that was introduced to model the interaction of a particular microbial species with its environment. This approach has been extended to analyzing interactions in a community of microbes; however, these approaches have two important drawbacks: first, one has to numerically solve a differential equation to identify the steady state, and second, there are no methods available to analyze the stability of the steady state. We propose a game theory based community FBA model wherein species compete to maximize their individual growth rate, and the state of the community is given by the resulting Nash equilibrium. We develop a computationally efficient method for directly computing the steady state biomasses and fluxes without solving a differential equation. We also develop a method to determine the stability of a steady state to perturbations in the biomasses and to invasion by new species. We report the results of applying our proposed framework to a small community of four E. coli mutants that compete for externally supplied glucose, as well as cooperate since the mutants are auxotrophic for metabolites exported by other mutants, and a more realistic model for a gut microbiome consisting of nine species.

Metabolic cross-feeding interactions modulate the dynamic community structure in microbial fuel cell under variable organic loading wastewaters

The efficiency of microbial fuel cells (MFCs) in industrial wastewater treatment is profoundly influenced by the microbial community, which can be disrupted by variable industrial operations. Although microbial guilds linked to MFC performance under specific conditions have been identified, comprehensive knowledge of the convergent community structure and pathways of adaptation is lacking. Here, we developed a microbe-microbe interaction genome-scale metabolic model (mmGEM) based on metabolic cross-feeding to study the adaptation of microbial communities in MFCs treating sulfide-containing wastewater from a canned-pineapple factory. The metabolic model encompassed three major microbial guilds: sulfate-reducing bacteria (SRB), methanogens (MET), and sulfide-oxidizing bacteria (SOB). Our findings revealed a shift from an SOB-dominant to MET-dominant community as organic loading rates (OLRs) increased, along with a decline in MFC performance. The mmGEM accurately predicted microbial relative abundance at low OLRs (L-OLRs) and adaptation to high OLRs (H-OLRs). The simulations revealed constraints on SOB growth under H-OLRs due to reduced sulfate-sulfide (S) cycling and acetate cross-feeding with SRB. More cross-fed metabolites from SRB were diverted to MET, facilitating their competitive dominance. Assessing cross-feeding dynamics under varying OLRs enabled the execution of practical scenario-based simulations to explore the potential impact of elevated acidity levels on SOB growth and MFC performance. This work highlights the role of metabolic cross-feeding in shaping microbial community structure in response to high OLRs. The insights gained will inform the development of effective strategies for implementing MFC technology in real-world industrial environments.

DyMMM-LEAPS: An ML-based framework for modulating evenness and stability in synthetic microbial communities

There have been a growing number of computational strategies to aid in the design of synthetic microbial consortia. A framework to identify regions in parametric space to maximize two essential properties, evenness and stability, is critical. In this study, we introduce DyMMM-LEAPS (dynamic multispecies metabolic modeling-locating evenness and stability in large parametric space), an extension of the DyMMM framework. Our method explores the large parametric space of genetic circuits in synthetic microbial communities to identify regions of evenness and stability. Due to the high computational costs of exhaustive sampling, we utilize adaptive sampling and surrogate modeling to reduce the number of simulations required to map the vast space. Our framework predicts engineering targets and computes their operating ranges to maximize the probability of the engineered community to have high evenness and stability. We demonstrate our approach by simulating five cocultures and one three-strain culture with different social interactions (cooperation, competition, and predation) employing quorum-sensing-based genetic circuits. In addition to guiding circuit tuning, our pipeline gives an opportunity for a detailed analysis of pockets of evenness and stability for the circuit under investigation, which can further help dissect the relationship between the two properties. DyMMM-LEAPS is easily customizable and can be expanded to a larger community with more complex interactions.

The optimization of microbial functions through rational environmental manipulations

Microorganisms play a central role in biotechnology and it is key that we develop strategies to engineer and optimize their functionality. To this end, most efforts have focused on introducing genetic manipulations in microorganisms which are then grown either in monoculture or in mixed-species consortia. An alternative strategy to optimize microbial processes is to rationally engineer the environment in which microbes grow. The microbial environment is multidimensional, including factors such as temperature, pH, salinity, nutrient composition, etc. These environmental factors all influence the growth and phenotypes of microorganisms and they generally "interact" with one another, combining their effects in complex, non-additive ways. In this piece, we overview the origins and consequences of these "interactions" between environmental factors and discuss how they have been built into statistical, bottom-up predictive models of microbial function to identify optimal environmental conditions for monocultures and microbial consortia. We also overview alternative "top-down" approaches, such as genetic algorithms, to finding optimal combinations of environmental factors. By providing a brief summary of the state of this field, we hope to stimulate further work on the rational manipulation and optimization of the microbial environment.

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

A new era of synthetic biology-microbial community design

Synthetic biology conceptualizes biological complexity as a network of biological parts, devices, and systems with predetermined functionalities and has had a revolutionary impact on fundamental and applied research. With the unprecedented ability to synthesize and transfer any DNA and RNA across organisms, the scope of synthetic biology is expanding and being recreated in previously unimaginable ways. The field has matured to a level where highly complex networks, such as artificial communities of synthetic organisms, can be constructed. In parallel, computational biology became an integral part of biological studies, with computational models aiding the unravelling of the escalating complexity and emerging properties of biological phenomena. However, there is still a vast untapped potential for the complete integration of modelling into the synthetic design process, presenting exciting opportunities for scientific advancements. Here, we first highlight the most recent advances in computer-aided design of microbial communities. Next, we propose that such a design can benefit from an organism-free modular modelling approach that places its emphasis on modules of organismal function towards the design of multispecies communities. We argue for a shift in perspective from single organism-centred approaches to emphasizing the functional contributions of organisms within the community. By assembling synthetic biological systems using modular computational models with mathematical descriptions of parts and circuits, we can tailor organisms to fulfil specific functional roles within the community. This approach aligns with synthetic biology strategies and presents exciting possibilities for the design of artificial communities. Graphical Abstract.

Synthetic microbial ecology: engineering habitats for modular consortia

Microbiomes, the complex networks of micro-organisms and the molecules through which they interact, play a crucial role in health and ecology. Over at least the past two decades, engineering biology has made significant progress, impacting the bio-based industry, health, and environmental sectors; but has only recently begun to explore the engineering of microbial ecosystems. The creation of synthetic microbial communities presents opportunities to help us understand the dynamics of wild ecosystems, learn how to manipulate and interact with existing microbiomes for therapeutic and other purposes, and to create entirely new microbial communities capable of undertaking tasks for industrial biology. Here, we describe how synthetic ecosystems can be constructed and controlled, focusing on how the available methods and interaction mechanisms facilitate the regulation of community composition and output. While experimental decisions are dictated by intended applications, the vast number of tools available suggests great opportunity for researchers to develop a diverse array of novel microbial ecosystems.

Revealing microbiota characteristics and predicting flavor-producing sub-communities in Nongxiangxing baijiu pit mud through metagenomic analysis and metabolic modeling

The microorganisms of the pit mud (PM) of Nongxiangxing baijiu (NXXB) have an important role in the synthesis of flavor substances, and they determine attributes and quality of baijiu. Herein, we utilize metagenomics and genome-scale metabolic models (GSMMs) to investigate the microbial composition, metabolic functions in PM microbiota, as well as to identify microorganisms and communities linked to flavor compounds. Metagenomic data revealed that the most prevalent assembly of bacteria and archaea was Proteiniphilum, Caproicibacterium, Petrimonas, Lactobacillus, Clostridium, Aminobacterium, Syntrophomonas, Methanobacterium, Methanoculleus, and Methanosarcina. The important enzymes ofPMwere in bothGH and GT familymetabolism. A total of 38 high-quality metagenome-assembled genomes (MAGs) were obtained, including those at the family level (n = 13), genus level (n = 17), and species level (n = 8). GSMMs of the 38 MAGs were then constructed. From the GSMMs, individual and community capabilities respectively were predicted to be able to produce 111 metabolites and 598 metabolites. Twenty-three predicted metabolites were consistent with the metabonomics detected flavors and served as targets. Twelve sub-community of were screened by cross-feeding of 38 GSMMs. Of them, Methanobacterium, Sphaerochaeta, Muricomes intestini, Methanobacteriaceae, Synergistaceae, and Caloramator were core microorganisms for targets in each sub-community. Overall, this study of metagenomic and target-community screening could help our understanding of the metabolite-microbiome association and further bioregulation of baijiu.

Emerging methods for genome-scale metabolic modeling of microbial communities

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels.

Rational construction of synthetic consortia: Key considerations and model-based methods for guiding the development of a novel biosynthesis platform

The rapid development of synthetic biology has significantly improved the capabilities of mono-culture systems in converting different substrates into various value-added bio-chemicals through metabolic engineering. However, overexpression of biosynthetic pathways in recombinant strains can impose a heavy metabolic burden on the host, resulting in imbalanced energy distribution and negatively affecting both cell growth and biosynthesis capacity. Synthetic consortia, consisting of two or more microbial species or strains with complementary functions, have emerged as a promising and efficient platform to alleviate the metabolic burden and increase product yield. However, research on synthetic consortia is still in its infancy, with numerous challenges regarding the design and construction of stable synthetic consortia. This review provides a comprehensive comparison of the advantages and disadvantages of mono-culture systems and synthetic consortia. Key considerations for engineering synthetic consortia based on recent advances are summarized, and simulation and computational tools for guiding the advancement of synthetic consortia are discussed. Moreover, further development of more efficient and cost-effective synthetic consortia with emerging technologies such as artificial intelligence and machine learning is highlighted.

Progress on network modeling and analysis of gut microecology: a review

The gut microecological network is a complex microbial community within the human body that plays a key role in linking dietary nutrition and host physiology. To understand the complex relationships among microbes and their functions within this community, network analysis has emerged as a powerful tool. By representing the interactions between microbes and their associated omics data as a network, we can gain a comprehensive understanding of the ecological mechanisms that drive the human gut microbiota. In addition, the network-based approach provides a more intuitive analysis of the gut microbiota, simplifying the study of its complex dynamics and interdependencies. This review provides a comprehensive overview of the methods used to construct and analyze networks in the context of gut microecological background. We discuss various types of network modeling approaches, including co-occurrence networks, causal networks, dynamic networks, and multi-omics networks, and describe the analytical techniques used to identify important network properties. We also highlight the challenges and limitations of network modeling in this area, such as data scarcity and heterogeneity, and provide future research directions to overcome these limitations. By exploring these network-based methods, researchers can gain valuable insights into the intricate relationships and functional roles of microbial communities within the gut, ultimately advancing our understanding of the gut microbiota's impact on human health.

Bacteriophage specificity is impacted by interactions between bacteria

Predators play a central role in shaping community structure, function, and stability. The degree to which bacteriophage predators (viruses that infect bacteria) evolve to be specialists with a single bacterial prey species versus generalists able to consume multiple types of prey has implications for their effect on microbial communities. The presence and abundance of multiple bacterial prey types can alter selection for phage generalists, but less is known about how interactions between prey shape predator specificity in microbial systems. Using a phenomenological mathematical model of phage and bacterial populations, we find that the dominant phage strategy depends on prey ecology. Given a fitness cost for generalism, generalist predators maintain an advantage when prey species compete, while specialists dominate when prey are obligately engaged in cross-feeding interactions. We test these predictions in a synthetic microbial community with interacting strains of Escherichia coli and Salmonella enterica by competing a generalist T5-like phage able to infect both prey against P22vir, an S. enterica-specific phage. Our experimental data conform to our modeling expectations when prey species are competing or obligately mutualistic, although our results suggest that the in vitro cost of generalism is caused by a combination of biological mechanisms not anticipated in our model. Our work demonstrates that interactions between bacteria play a role in shaping ecological selection on predator specificity in obligately lytic bacteriophages and emphasizes the diversity of ways in which fitness trade-offs can manifest.

IMPORTANCE

There is significant natural diversity in how many different types of bacteria a bacteriophage can infect, but the mechanisms driving this diversity are unclear. This study uses a combination of mathematical modeling and an in vitro system consisting of Escherichia coli, Salmonella enterica, a T5-like generalist phage, and the specialist phage P22vir to highlight the connection between bacteriophage specificity and interactions between their potential microbial prey. Mathematical modeling suggests that competing bacteria tend to favor generalist bacteriophage, while bacteria that benefit each other tend to favor specialist bacteriophage. Experimental results support this general finding. The experiments also show that the optimal phage strategy is impacted by phage degradation and bacterial physiology. These findings enhance our understanding of how complex microbial communities shape selection on bacteriophage specificity, which may improve our ability to use phage to manage antibiotic-resistant microbial infections.

Biophysical metabolic modeling of complex bacterial colony morphology

Microbial colony growth is shaped by the physics of biomass propagation and nutrient diffusion, and by the metabolic reactions that organisms activate as a function of the surrounding environment. While microbial colonies have been explored using minimal models of growth and motility, full integration of biomass propagation and metabolism is still lacking. Here, building upon our framework for Computation of Microbial Ecosystems in Time and Space (COMETS), we combine dynamic flux balance modeling of metabolism with collective biomass propagation and demographic fluctuations to provide nuanced simulations of E. coli colonies. Simulations produced realistic colony morphology, consistent with our experiments. They characterize the transition between smooth and furcated colonies and the decay of genetic diversity. Furthermore, we demonstrate that under certain conditions, biomass can accumulate along "metabolic rings" that are reminiscent of coffee-stain rings, but have a completely different origin. Our approach is a key step towards predictive microbial ecosystems modeling.

Life at the borderlands: microbiomes of interfaces critical to One Health

Microbiomes are foundational components of the environment that provide essential services relating to food security, carbon sequestration, human health, and the overall well-being of ecosystems. Microbiota exert their effects primarily through complex interactions at interfaces with their plant, animal, and human hosts, as well as within the soil environment. This review aims to explore the ecological, evolutionary, and molecular processes governing the establishment and function of microbiome-host relationships, specifically at interfaces critical to One Health-a transdisciplinary framework that recognizes that the health outcomes of people, animals, plants, and the environment are tightly interconnected. Within the context of One Health, the core principles underpinning microbiome assembly will be discussed in detail, including biofilm formation, microbial recruitment strategies, mechanisms of microbial attachment, community succession, and the effect these processes have on host function and health. Finally, this review will catalogue recent advances in microbiology and microbial ecology methods that can be used to profile microbial interfaces, with particular attention to multi-omic, advanced imaging, and modelling approaches. These technologies are essential for delineating the general and specific principles governing microbiome assembly and functions, mapping microbial interconnectivity across varying spatial and temporal scales, and for the establishment of predictive frameworks that will guide the development of targeted microbiome-interventions to deliver One Health outcomes.

Applications of synthetic yeast consortia for the production of native and non-native chemicals

The application of microbial consortia is a new approach in synthetic biology. Synthetic yeast consortia, simple or complex synthetic mixed cultures, have been used for the production of various metabolites. Cooperation between the members of a consortium and cross-feeding can be applied to create stable microbial communication. These consortia can: consume a variety of substrates, perform more complex functions, produce metabolites in high titer, rate, and yield (TRY), and show higher stability during industrial fermentations. Due to the new research context of synthetic consortia, few yeasts were used to build these consortia, including Saccharomyces cerevisiae, Pichia pastoris, and Yarrowia lipolytica. Here, application of the yeasts for design of synthetic microbial consortia and their advantages and bottlenecks for effective and robust production of valuable metabolites from bioresource, including: cellulose, xylose, glycerol and so on, have been reviewed. Key trends and challenges are also discussed for the future development of synthetic yeast consortia.

Harnessing root exudates for plant microbiome engineering and stress resistance in plants

A wide range of abiotic and biotic stresses adversely affect plant's growth and production. Under stress, one of the main responses of plants is the modulation of exudates excreted in the rhizosphere, which consequently leads to alterations in the resident microbiota. Thus, the exudates discharged into the rhizospheric environment play a preponderant role in the association and formation of plant-microbe interactions. In this review, we aimed to provide a synthesis of the latest and most pertinent literature on the diverse biochemical and structural compositions of plant root exudates. Also, this work investigates into their multifaceted role in microbial nutrition and intricate signaling processes within the rhizosphere, which includes quorum-sensing molecules. Specifically, it explores the contributions of low molecular weight compounds, such as carbohydrates, phenolics, organic acids, amino acids, and secondary metabolites, as well as the significance of high molecular weight compounds, including proteins and polysaccharides. It also discusses the state-of-the-art omics strategies that unveil the vital role of root exudates in plant-microbiome interactions, including defense against pathogens like nematodes and fungi. We propose multiple challenges and perspectives, including exploiting plant root exudates for host-mediated microbiome engineering. In this discourse, root exudates and their derived interactions with the rhizospheric microbiota should receive greater attention due to their positive influence on plant health and stress mitigation.

Predicting microbial interactions with approaches based on flux balance analysis: an evaluation

BACKGROUND

Given a genome-scale metabolic model (GEM) of a microorganism and criteria for optimization, flux balance analysis (FBA) predicts the optimal growth rate and its corresponding flux distribution for a specific medium. FBA has been extended to microbial consortia and thus can be used to predict interactions by comparing in-silico growth rates for co- and monocultures. Although FBA-based methods for microbial interaction prediction are becoming popular, a systematic evaluation of their accuracy has not yet been performed.

RESULTS

Here, we evaluate the accuracy of FBA-based predictions of human and mouse gut bacterial interactions using growth data from the literature. For this, we collected 26 GEMs from the semi-curated AGORA database as well as four previously published curated GEMs. We tested the accuracy of three tools (COMETS, Microbiome Modeling Toolbox and MICOM) by comparing growth rates predicted in mono- and co-culture to growth rates extracted from the literature and also investigated the impact of different tool settings and media. We found that except for curated GEMs, predicted growth rates and their ratios (i.e. interaction strengths) do not correlate with growth rates and interaction strengths obtained from in vitro data.

CONCLUSIONS

Prediction of growth rates with FBA using semi-curated GEMs is currently not sufficiently accurate to predict interaction strengths reliably.

Strategies for tailoring functional microbial synthetic communities

Natural ecosystems harbor a huge reservoir of taxonomically diverse microbes that are important for plant growth and health. The vast diversity of soil microorganisms and their complex interactions make it challenging to pinpoint the main players important for the life support functions microbes can provide to plants, including enhanced tolerance to (a)biotic stress factors. Designing simplified microbial synthetic communities (SynComs) helps reduce this complexity to unravel the molecular and chemical basis and interplay of specific microbiome functions. While SynComs have been successfully employed to dissect microbial interactions or reproduce microbiome-associated phenotypes, the assembly and reconstitution of these communities have often been based on generic abundance patterns or taxonomic identities and co-occurrences but have only rarely been informed by functional traits. Here, we review recent studies on designing functional SynComs to reveal common principles and discuss multidimensional approaches for community design. We propose a strategy for tailoring the design of functional SynComs based on integration of high-throughput experimental assays with microbial strains and computational genomic analyses of their functional capabilities.

Emergent antibiotic persistence in a spatially structured synthetic microbial mutualism

Antibiotic persistence (heterotolerance) allows a subpopulation of bacteria to survive antibiotic-induced killing and contributes to the evolution of antibiotic resistance. Although bacteria typically live in microbial communities with complex ecological interactions, little is known about how microbial ecology affects antibiotic persistence. Here, we demonstrated within a synthetic two-species microbial mutualism of Escherichia coli and Salmonella enterica that the combination of cross-feeding and community spatial structure can emergently cause high antibiotic persistence in bacteria by increasing the cell-to-cell heterogeneity. Tracking ampicillin-induced death for bacteria on agar surfaces, we found that E. coli forms up to 55 times more antibiotic persisters in the cross-feeding coculture than in monoculture. This high persistence could not be explained solely by the presence of S. enterica, the presence of cross-feeding, average nutrient starvation, or spontaneous resistant mutations. Time-series fluorescent microscopy revealed increased cell-to-cell variation in E. coli lag time in the mutualistic co-culture. Furthermore, we discovered that an E. coli cell can survive antibiotic killing if the nearby S. enterica cells on which it relies die first. In conclusion, we showed that the high antibiotic persistence phenotype can be an emergent phenomenon caused by a combination of cross-feeding and spatial structure. Our work highlights the importance of considering spatially structured interactions during antibiotic treatment and understanding microbial community resilience more broadly.

Metaphors as design tools for microbial consortia: An analysis of recent peer-reviewed literature

Single engineered microbial species cannot always conduct complex transformations, while complex, incompletely defined microbial consortia have heretofore been suited to a limited range of tasks. As biodesigners bridge this gap with intentionally designed microbial communities, they will, intentionally or otherwise, build communities that embody particular ideas about what microbial communities can and should be. Here, we suggest that metaphors-ideas about what microbial communities are like-are therefore important tools for designing synthetic consortia-based bioreactors. We identify a range of metaphors currently employed in peer-reviewed microbiome research articles, characterizing each through its potential structural implications and distinctive imagery. We present this metaphor catalogue in the interest of, first, making metaphors visible as design choices, second, enabling deliberate experimentation with them towards expanding the potential design space of the field, and third, encouraging reflection on the goals and values they embed.

Mutualism reduces the severity of gene disruptions in predictable ways across microbial communities

Predicting evolution in microbial communities is critical for problems from human health to global nutrient cycling. Understanding how species interactions impact the distribution of fitness effects for a focal population would enhance our ability to predict evolution. Specifically, does the type of ecological interaction, such as mutualism or competition, change the average effect of a mutation (i.e., the mean of the distribution of fitness effects)? Furthermore, how often does increasing community complexity alter the impact of species interactions on mutant fitness? To address these questions, we created a transposon mutant library in Salmonella enterica and measured the fitness of loss of function mutations in 3,550 genes when grown alone versus competitive co-culture or mutualistic co-culture with Escherichia coli and Methylorubrum extorquens. We found that mutualism reduces the average impact of mutations, while competition had no effect. Additionally, mutant fitness in the 3-species communities can be predicted by averaging the fitness in each 2-species community. Finally, we discovered that in the mutualism S. enterica obtained vitamins and more amino acids than previously known. Our results suggest that species interactions can predictably impact fitness effect distributions, in turn suggesting that evolution may ultimately be predictable in multi-species communities.

Xylose and shikimate transporters facilitates microbial consortium as a chassis for benzylisoquinoline alkaloid production

Plant-sourced aromatic amino acid (AAA) derivatives are a vast group of compounds with broad applications. Here, we present the development of a yeast consortium for efficient production of (S)-norcoclaurine, the key precursor for benzylisoquinoline alkaloid biosynthesis. A xylose transporter enables the concurrent mixed-sugar utilization in Scheffersomyces stipitis, which plays a crucial role in enhancing the flux entering the highly regulated shikimate pathway located upstream of AAA biosynthesis. Two quinate permeases isolated from Aspergillus niger facilitates shikimate translocation to the co-cultured Saccharomyces cerevisiae that converts shikimate to (S)-norcoclaurine, resulting in the maximal titer (11.5 mg/L), nearly 110-fold higher than the titer reported for an S. cerevisiae monoculture. Our findings magnify the potential of microbial consortium platforms for the economical de novo synthesis of complex compounds, where pathway modularization and compartmentalization in distinct specialty strains enable effective fine-tuning of long biosynthetic pathways and diminish intermediate buildup, thereby leading to increases in production.

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.

Inferring microbial interactions with their environment from genomic and metagenomic data

Microbial communities assemble through a complex set of interactions between microbes and their environment, and the resulting metabolic impact on the host ecosystem can be profound. Microbial activity is known to impact human health, plant growth, water quality, and soil carbon storage which has lead to the development of many approaches and products meant to manipulate the microbiome. In order to understand, predict, and improve microbial community engineering, genome-scale modeling techniques have been developed to translate genomic data into inferred microbial dynamics. However, these techniques rely heavily on simulation to draw conclusions which may vary with unknown parameters or initial conditions, rather than more robust qualitative analysis. To better understand microbial community dynamics using genome-scale modeling, we provide a tool to investigate the network of interactions between microbes and environmental metabolites over time. Using our previously developed algorithm for simulating microbial communities from genome-scale metabolic models (GSMs), we infer the set of microbe-metabolite interactions within a microbial community in a particular environment. Because these interactions depend on the available environmental metabolites, we refer to the networks that we infer as metabolically contextualized, and so name our tool MetConSIN: Metabolically Contextualized Species Interaction Networks.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

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'.

Mutualism reduces the severity of gene disruptions in predictable ways across microbial communities

Predicting evolution in microbial communities is critical for problems from human health to global nutrient cycling. Understanding how species interactions impact the distribution of fitness effects for a focal population would enhance our ability to predict evolution. Specifically, it would be useful to know if the type of ecological interaction, such as mutualism or competition, changes the average effect of a mutation (i.e., the mean of the distribution of fitness effects). Furthermore, how often does increasing community complexity alter the impact of species interactions on mutant fitness? To address these questions, we created a transposon mutant library in Salmonella enterica and measured the fitness of loss of function mutations in 3,550 genes when grown alone versus competitive co-culture or mutualistic co-culture with Escherichia coli and Methylorubrum extorquens. We found that mutualism reduces the average impact of mutations, while competition had no effect. Additionally, mutant fitness in the 3-species communities can be predicted by averaging the fitness in each 2-species community. Finally, the fitness effects of several knockouts in the mutualistic communities were surprising. We discovered that S. enterica is obtaining a different source of carbon and more vitamins and amino acids than we had expected. Our results suggest that species interactions can predictably impact fitness effect distributions, in turn suggesting that evolution may ultimately be predictable in multi-species communities.

Differences in Exudates Between Strains of Chlorella sorokiniana Affect the Interaction with the Microalga Growth-Promoting Bacteria Azospirillum brasilense : Differences in Exudates Between Strains of Chlorella sorokiniana Affect the Interaction with the Microalga Growth-Promoting Bacteria Azospirillum brasilense

The microalga Chlorella sorokiniana and the microalgae growth-promoting bacteria (MGPB) Azospirillum brasilense have a mutualistic interaction that can begin within the first hours of co-incubation; however, the metabolites participating in this initial interaction are not yet identified. Nuclear magnetic resonance (NMR) was used in the present study to characterize the metabolites exuded by two strains of C. sorokiniana (UTEX 2714 and UTEX 2805) and A. brasilense Cd when grown together in an oligotrophic medium. Lactate and myo-inositol were identified as carbon metabolites exuded by the two strains of C. sorokiniana; however, only the UTEX 2714 strain exuded glycerol as the main carbon compound. In turn, A. brasilense exuded uracil when grown on the exudates of either microalga, and both microalga strains were able to utilize uracil as a nitrogen source. Interestingly, although the total carbohydrate content was higher in exudates from C. sorokiniana UTEX 2805 than from C. sorokiniana UTEX 2714, the growth of A. brasilense was greater in the exudates from the UTEX 2714 strain. These results highlight the fact that in the exuded carbon compounds differ between strains of the same species of microalgae and suggest that the type, rather than the quantity, of carbon source is more important for sustaining the growth of the partner bacteria.

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.