Dynamic flux balance analysis for synthetic microbial communities

Modeling bacterial interactions uncovers the importance of outliers in the coastal lignin-degrading consortium

Lignin, as the abundant carbon polymer, is essential for carbon cycle and biorefinery. Microorganisms interact to form communities for lignin biodegradation, yet it is a challenge to understand such complex interactions. Here, we develop a coastal lignin-degrading bacterial consortium (LD), through "top-down" enrichment. Sequencing and physiological analyses reveal that LD is dominated by the lignin degrader Pluralibacter gergoviae (>98%), with additional rare non-degraders. Interestingly, LD, cultured in lignin-MB medium, significantly enhances cell growth and lignin degradation as compared to P. gergoviae alone, implying a role of additional outliers. Using genome-scale metabolic models, metabolic profiling and culture experiments, modeling of inter-species interactions between P. gergoviae, Vibrio alginolyticus, Aeromonas hydrophila and Shewanella putrefaciens, unravels cross-feeding of amino acids, organic acids and alcohols between the degrader and non-degraders. Furthermore, the sub-population ratio is essential to enforce the synergy. Our study highlights the unrecognized role of outliers in lignin degradation.

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.

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Food-breastmilk combinations alter the colonic microbiome of weaning infants: an in silico study

The introduction of solid foods to infants, also known as weaning, is a critical point for the development of the complex microbial community inhabiting the human colon, impacting host physiology in infancy and later in life. This research investigated in silico the impact of food-breastmilk combinations on growth and metabolite production by colonic microbes of New Zealand weaning infants using the metagenome-scale metabolic model named Microbial Community. Eighty-nine foods were individually combined with breastmilk, and the 12 combinations with the strongest influence on the microbial production of short-chain fatty acids (SCFAs) and branched-chain fatty acids (BCFAs) were identified. Fiber-rich and polyphenol-rich foods, like pumpkin and blackcurrant, resulted in the greatest increase in predicted fluxes of total SCFAs and individual fluxes of propionate and acetate when combined, respectively, with breastmilk. Identified foods were further combined with other foods and breastmilk, resulting in 66 multiple food-breastmilk combinations. These combinations altered in silico the impact of individual foods on the microbial production of SCFAs and BCFAs, suggesting that the interaction between the dietary compounds composing a meal is the key factor influencing colonic microbes. Blackcurrant combined with other foods and breastmilk promoted the greatest increase in the production of acetate and total SCFAs, while pork combined with other foods and breastmilk decreased the production of total BCFAs.IMPORTANCELittle is known about the influence of complementary foods on the colonic microbiome of weaning infants. Traditional in vitro and in vivo microbiome methods are limited by their resource-consuming concerns. Modeling approaches represent a promising complementary tool to provide insights into the behavior of microbial communities. This study evaluated how foods combined with other foods and human milk affect the production of short-chain fatty acids and branched-chain fatty acids by colonic microbes of weaning infants using a rapid and inexpensive in silico approach. Foods and food combinations identified here are candidates for future experimental investigations, helping to fill a crucial knowledge gap in infant nutrition.

Microalgae growth-promoting bacteria for cultivation strategies: Recent updates and progress

Microalgae growth-promoting bacteria (MGPB), both actinobacteria and non-actinobacteria, have received considerable attention recently because of their potential to develop microalgae-bacteria co-culture strategies for improved efficiency and sustainability of the water-energy-environment nexus. Owing to their diverse metabolic pathways and ability to adapt to diverse conditions, microalgal-MGPB co-cultures could be promising biological systems under uncertain environmental and nutrient conditions. This review proposes the recent updates and progress on MGPB for microalgae cultivation through co-culture strategies. Firstly, potential MGPB strains for microalgae cultivation are introduced. Following, microalgal-MGPB interaction mechanisms and applications of their co-cultures for biomass production and wastewater treatment are reviewed. Moreover, state-of-the-art studies on synthetic biology and metabolic network analysis, along with the challenges and prospects of opting these approaches for microalgal-MGPB co-cultures are presented. It is anticipated that these strategies may significantly improve the sustainability of microalgal-MGPB co-cultures for wastewater treatment, biomass valorization, and bioproducts synthesis in a circular bioeconomy paradigm.

A long-term growth stable Halomonas sp. deleted with multiple transposases guided by its metabolic network model Halo-ecGEM

Microbial instability is a common problem during bio-production based on microbial hosts. Halomonas bluephagenesis has been developed as a chassis for next generation industrial biotechnology (NGIB) under open and unsterile conditions. However, the hidden genomic information and peculiar metabolism have significantly hampered its deep exploitation for cell-factory engineering. Based on the freshly completed genome sequence of H. bluephagenesis TD01, which reveals 1889 biological process-associated genes grouped into 84 GO-slim terms. An enzyme constrained genome-scale metabolic model Halo-ecGEM was constructed, which showed strong ability to simulate fed-batch fermentations. A visible salt-stress responsive landscape was achieved by combining GO-slim term enrichment and CVT-based omics profiling, demonstrating that cells deploy most of the protein resources by force to support the essential activity of translation and protein metabolism when exposed to salt stress. Under the guidance of Halo-ecGEM, eight transposases were deleted, leading to a significantly enhanced stability for its growth and bioproduction of various polyhydroxyalkanoates (PHA) including 3-hydroxybutyrate (3HB) homopolymer PHB, 3HB and 3-hydroxyvalerate (3HV) copolymer PHBV, as well as 3HB and 4-hydroxyvalerate (4HB) copolymer P34HB. This study sheds new light on the metabolic characteristics and stress-response landscape of H. bluephagenesis, achieving for the first time to construct a long-term growth stable chassis for industrial applications. For the first time, it was demonstrated that genome encoded transposons are the reason for microbial instability during growth in flasks and fermentors.

Boundary flux analysis: an emerging strategy for investigating metabolic pathway activity in large cohorts

Many biological phenotypes are rooted in metabolic pathway activity rather than the concentrations of individual metabolites. Despite this, most metabolomics studies only capture steady-state metabolism - not metabolic flux. Although sophisticated metabolic flux analysis strategies have been developed, these methods are technically challenging and difficult to implement in large-cohort studies. Recently, a new boundary flux analysis (BFA) approach has emerged that captures large-scale metabolic flux phenotypes by quantifying changes in metabolite levels in the media of cultured cells. This approach is advantageous because it is relatively easy to implement yet captures complex metabolic flux phenotypes. We describe the opportunities and challenges of BFA and illustrate how it can be harnessed to investigate a wide transect of biological phenomena.

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.

Stochastic biological system-of-systems modelling for iPSC culture

Large-scale manufacturing of induced pluripotent stem cells (iPSCs) is essential for cell therapies and regenerative medicines. Yet, iPSCs form large cell aggregates in suspension bioreactors, resulting in insufficient nutrient supply and extra metabolic waste build-up for the cells located at the core. Since subtle changes in micro-environment can lead to a heterogeneous cell population, a novel Biological System-of-Systems (Bio-SoS) framework is proposed to model cell-to-cell interactions, spatial and metabolic heterogeneity, and cell response to micro-environmental variation. Building on stochastic metabolic reaction network, aggregation kinetics, and reaction-diffusion mechanisms, the Bio-SoS model characterizes causal interdependencies at individual cell, aggregate, and cell population levels. It has a modular design that enables data integration and improves predictions for different monolayer and aggregate culture processes. In addition, a variance decomposition analysis is derived to quantify the impact of factors (i.e., aggregate size) on cell product health and quality heterogeneity.

Microbial interactions from a new perspective: reinforcement learning reveals new insights into microbiome evolution

MOTIVATION

Microbes are essential part of all ecosystems, influencing material flow and shaping their surroundings. Metabolic modeling has been a useful tool and provided tremendous insights into microbial community metabolism. However, current methods based on flux balance analysis (FBA) usually fail to predict metabolic and regulatory strategies that lead to long-term survival and stability especially in heterogenous communities.

RESULTS

Here, we introduce a novel reinforcement learning algorithm, Self-Playing Microbes in Dynamic FBA, which treats microbial metabolism as a decision-making process, allowing individual microbial agents to evolve by learning and adapting metabolic strategies for enhanced long-term fitness. This algorithm predicts what microbial flux regulation policies will stabilize in the dynamic ecosystem of interest in the presence of other microbes with minimal reliance on predefined strategies. Throughout this article, we present several scenarios wherein our algorithm outperforms existing methods in reproducing outcomes, and we explore the biological significance of these predictions.

AVAILABILITY AND IMPLEMENTATION

The source code for this article is available at: https://github.com/chan-csu/SPAM-DFBA.

'Social' versus 'asocial' cells-dynamic competition flux balance analysis

In multicellular organisms cells compete for resources or growth factors. If any one cell type wins, the co-existence of diverse cell types disappears. Existing dynamic Flux Balance Analysis (dFBA) does not accommodate changes in cell density caused by competition. Therefore we here develop 'dynamic competition Flux Balance Analysis' (dcFBA). With total biomass synthesis as objective, lower-growth-yield cells were outcompeted even when cells synthesized mutually required nutrients. Signal transduction between cells established co-existence, which suggests that such 'socialness' is required for multicellularity. Whilst mutants with increased specific growth rate did not outgrow the other cell types, loss of social characteristics did enable a mutant to outgrow the other cells. We discuss that 'asocialness' rather than enhanced growth rates, i.e., a reduced sensitivity to regulatory factors rather than enhanced growth rates, may characterize cancer cells and organisms causing ecological blooms. Therapies reinforcing cross-regulation may therefore be more effective than those targeting replication rates.

Metabolic engineering of human gut microbiome: Recent developments and future perspectives

Many studies have demonstrated that the gut microbiota is associated with human health and disease. Manipulation of the gut microbiota, e.g. supplementation of probiotics, has been suggested to be feasible, but subject to limited therapeutic efficacy. To develop efficient microbiota-targeted diagnostic and therapeutic strategies, metabolic engineering has been applied to construct genetically modified probiotics and synthetic microbial consortia. This review mainly discusses commonly adopted strategies for metabolic engineering in the human gut microbiome, including the use of in silico, in vitro, or in vivo approaches for iterative design and construction of engineered probiotics or microbial consortia. Especially, we highlight how genome-scale metabolic models can be applied to advance our understanding of the gut microbiota. Also, we review the recent applications of metabolic engineering in gut microbiome studies as well as discuss important challenges and opportunities.

Elucidating dynamic anaerobe metabolism with HRMAS 13C NMR and genome-scale modeling

Anaerobic microbial metabolism drives critical functions within global ecosystems, host-microbiota interactions, and industrial applications, yet remains ill-defined. Here we advance a versatile approach to elaborate cellular metabolism in obligate anaerobes using the pathogen Clostridioides difficile, an amino acid and carbohydrate-fermenting Clostridia. High-resolution magic angle spinning nuclear magnetic resonance (NMR) spectroscopy of C. difficile, grown with fermentable 13C substrates, informed dynamic flux balance analysis (dFBA) of the pathogen's genome-scale metabolism. Analyses identified dynamic recruitment of oxidative and supporting reductive pathways, with integration of high-flux amino acid and glycolytic metabolism at alanine's biosynthesis to support efficient energy generation, nitrogen handling and biomass generation. Model predictions informed an approach leveraging the sensitivity of 13C NMR spectroscopy to simultaneously track cellular carbon and nitrogen flow from [U-13C]glucose and [15N]leucine, confirming the formation of [13C,15N]alanine. Findings identify metabolic strategies used by C. difficile to support its rapid colonization and expansion in gut ecosystems.

Synthetic microbial communities (SynComs) of the human gut: design, assembly, and applications

The human gut harbors native microbial communities, forming a highly complex ecosystem. Synthetic microbial communities (SynComs) of the human gut are an assembly of microorganisms isolated from human mucosa or fecal samples. In recent decades, the ever-expanding culturing capacity and affordable sequencing, together with advanced computational modeling, started a ''golden age'' for harnessing the beneficial potential of SynComs to fight gastrointestinal disorders, such as infections and chronic inflammatory bowel diseases. As simplified and completely defined microbiota, SynComs offer a promising reductionist approach to understanding the multispecies and multikingdom interactions in the microbe-host-immune axis. However, there are still many challenges to overcome before we can precisely construct SynComs of designed function and efficacy that allow the translation of scientific findings to patients' treatments. Here, we discussed the strategies used to design, assemble, and test a SynCom, and address the significant challenges, which are of microbiological, engineering, and translational nature, that stand in the way of using SynComs as live bacterial therapeutics.

Metabolic models of human gut microbiota: Advances and challenges

The human gut is a complex ecosystem consisting of hundreds of microbial species interacting with each other and with the human host. Mathematical models of the gut microbiome integrate our knowledge of this system and help to formulate hypotheses to explain observations. The generalized Lotka-Volterra model has been widely used for this purpose, but it does not describe interaction mechanisms and thus does not account for metabolic flexibility. Recently, models that explicitly describe gut microbial metabolite production and consumption have become popular. These models have been used to investigate the factors that shape gut microbial composition and to link specific gut microorganisms to changes in metabolite concentrations found in diseases. Here, we review how such models are built and what we have learned so far from their application to human gut microbiome data. In addition, we discuss current challenges of these models and how these can be addressed in the future.

Alpha-Hemolysin from Staphylococcus aureus Obstructs Yeast-Hyphae Switching and Diminishes Pathogenicity in Candida albicans

The use of antibiotics can disrupt the body's natural balance and increase the susteptibility of patients towards fungal infections. Candida albicans is a dimorphic opportunistic fungal pathogen with niches similar to those of bacteria. Our aim was to study the interaction between this pathogen and bacteria to facilitate the control of C. albicans infection. Alpha-hemolysin (Hla), a protein secreted from Staphylococcus aureus, causes cell wall damage and impedes the yeast-hyphae transition in C. albicans. Mechanistically, Hla stimulation triggered the formation of reactive oxygen species that damaged the cell wall and mitochondria of C. albicans. The cell cycle was arrested in the G0/G1 phase, CDC42 was downregulated, and Ywp1 was upregulated, disrupting yeast hyphae switching. Subsequently, hyphae development was inhibited. In mouse models, C. albicans pretreated with Hla reduced the C. albicans burden in skin and vaginal mucosal infections, suggesting that S. aureus Hla can inhibit hyphal development and reduce the pathogenicity of candidiasis in vivo.

Optimization and Scale-Up of Fermentation Processes Driven by Models

In the era of sustainable development, the use of cell factories to produce various compounds by fermentation has attracted extensive attention; however, industrial fermentation requires not only efficient production strains, but also suitable extracellular conditions and medium components, as well as scaling-up. In this regard, the use of biological models has received much attention, and this review will provide guidance for the rapid selection of biological models. This paper first introduces two mechanistic modeling methods, kinetic modeling and constraint-based modeling (CBM), and generalizes their applications in practice. Next, we review data-driven modeling based on machine learning (ML), and highlight the application scope of different learning algorithms. The combined use of ML and CBM for constructing hybrid models is further discussed. At the end, we also discuss the recent strategies for predicting bioreactor scale-up and culture behavior through a combination of biological models and computational fluid dynamics (CFD) models.

Bottom-up synthetic ecology study of microbial consortia to enhance lignocellulose bioconversion

Lignocellulose is the most abundant organic carbon polymer on the earth. Its decomposition and conversion greatly impact the global carbon cycle. Furthermore, it provides feedstock for sustainable fuel and other value-added products. However, it continues to be underutilized, due to its highly recalcitrant and heterogeneric structure. Microorganisms, which have evolved versatile pathways to convert lignocellulose, undoubtedly are at the heart of lignocellulose conversion. Numerous studies that have reported successful metabolic engineering of individual strains to improve biological lignin valorization. Meanwhile, the bottleneck of single strain modification is becoming increasingly urgent in the conversion of complex substrates. Alternatively, increased attention has been paid to microbial consortia, as they show advantages over pure cultures, e.g., high efficiency and robustness. Here, we first review recent developments in microbial communities for lignocellulose bioconversion. Furthermore, the emerging area of synthetic ecology, which is an integration of synthetic biology, ecology, and computational biology, provides an opportunity for the bottom-up construction of microbial consortia. Then, we review different modes of microbial interaction and their molecular mechanisms, and discuss considerations of how to employ these interactions to construct synthetic consortia via synthetic ecology, as well as highlight emerging trends in engineering microbial communities for lignocellulose bioconversion.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Advances in flux balance analysis by integrating machine learning and mechanism-based models

The availability of multi-omics data sets and genome-scale metabolic models for various organisms provide a platform for modeling and analyzing genotype-to-phenotype relationships. Flux balance analysis is the main tool for predicting flux distributions in genome-scale metabolic models and various data-integrative approaches enable modeling context-specific network behavior. Due to its linear nature, this optimization framework is readily scalable to multi-tissue or -organ and even multi-organism models. However, both data and model size can hamper a straightforward biological interpretation of the estimated fluxes. Moreover, flux balance analysis simulates metabolism at steady-state and thus, in its most basic form, does not consider kinetics or regulatory events. The integration of flux balance analysis with complementary data analysis and modeling techniques offers the potential to overcome these challenges. In particular machine learning approaches have emerged as the tool of choice for data reduction and selection of most important variables in big data sets. Kinetic models and formal languages can be used to simulate dynamic behavior. This review article provides an overview of integrative studies that combine flux balance analysis with machine learning approaches, kinetic models, such as physiology-based pharmacokinetic models, and formal graphical modeling languages, such as Petri nets. We discuss the mathematical aspects and biological applications of these integrated approaches and outline challenges and future perspectives.

Modelling microbial communities: Harnessing consortia for biotechnological applications

Microbes propagate and thrive in complex communities, and there are many benefits to studying and engineering microbial communities instead of single strains. Microbial communities are being increasingly leveraged in biotechnological applications, as they present significant advantages such as the division of labour and improved substrate utilisation. Nevertheless, they also present some interesting challenges to surmount for the design of efficient biotechnological processes. In this review, we discuss key principles of microbial interactions, followed by a deep dive into genome-scale metabolic models, focussing on a vast repertoire of constraint-based modelling methods that enable us to characterise and understand the metabolic capabilities of microbial communities. Complementary approaches to model microbial communities, such as those based on graph theory, are also briefly discussed. Taken together, these methods provide rich insights into the interactions between microbes and how they influence microbial community productivity. We finally overview approaches that allow us to generate and test numerous synthetic community compositions, followed by tools and methodologies that can predict effective genetic interventions to further improve the productivity of communities. With impending advancements in high-throughput omics of microbial communities, the stage is set for the rapid expansion of microbial community engineering, with a significant impact on biotechnological processes.

An evolutionary algorithm for designing microbial communities via environmental modification

Despite a growing understanding of how environmental composition affects microbial communities, it remains difficult to apply this knowledge to the rational design of synthetic multispecies consortia. This is because natural microbial communities can harbour thousands of different organisms and environmental substrates, making up a vast combinatorial space that precludes exhaustive experimental testing and computational prediction. Here, we present a method based on the combination of machine learning and metabolic modelling that selects optimal environmental compositions to produce target community phenotypes. In this framework, dynamic flux balance analysis is used to model the growth of a community in candidate environments. A genetic algorithm is then used to evaluate the behaviour of the community relative to a target phenotype, and subsequently adjust the environment to allow the organisms to approach this target. We apply this iterative process to thousands of in silico communities of varying sizes, showing how it can rapidly identify environments that yield desired taxonomic compositions and patterns of metabolic exchange. Moreover, this combination of approaches produces testable predictions for the assembly of experimental microbial communities with specific properties and can facilitate rational environmental design processes for complex microbiomes.

Microfluidic and mathematical modeling of aquatic microbial communities

Aquatic microbial communities contribute fundamentally to biogeochemical transformations in natural ecosystems, and disruption of these communities can lead to ecological disasters such as harmful algal blooms. Microbial communities are highly dynamic, and their composition and function are tightly controlled by the biophysical (e.g., light, fluid flow, and temperature) and biochemical (e.g., chemical gradients and cell concentration) parameters of the surrounding environment. Due to the large number of environmental factors involved, a systematic understanding of the microbial community-environment interactions is lacking. In this article, we show that microfluidic platforms present a unique opportunity to recreate well-defined environmental factors in a laboratory setting in a high throughput way, enabling quantitative studies of microbial communities that are amenable to theoretical modeling. The focus of this article is on aquatic microbial communities, but the microfluidic and mathematical models discussed here can be readily applied to investigate other microbiomes.

Understanding the host-microbe interactions using metabolic modeling

The human gut harbors an enormous number of symbiotic microbes, which is vital for human health. However, interactions within the complex microbiota community and between the microbiota and its host are challenging to elucidate, limiting development in the treatment for a variety of diseases associated with microbiota dysbiosis. Using in silico simulation methods based on flux balance analysis, those interactions can be better investigated. Flux balance analysis uses an annotated genome-scale reconstruction of a metabolic network to determine the distribution of metabolic fluxes that represent the complete metabolism of a bacterium in a certain metabolic environment such as the gut. Simulation of a set of bacterial species in a shared metabolic environment can enable the study of the effect of numerous perturbations, such as dietary changes or addition of a probiotic species in a personalized manner. This review aims to introduce to experimental biologists the possible applications of flux balance analysis in the host-microbiota interaction field and discusses its potential use to improve human health. Video abstract.

Metabolic modelling approaches for describing and engineering microbial communities

Microbes do not live in isolation but in microbial communities. The relevance of microbial communities is increasing due to growing awareness of their influence on a huge number of environmental, health and industrial processes. Hence, being able to control and engineer the output of both natural and synthetic communities would be of great interest. However, most of the available methods and biotechnological applications involving microorganisms, both in vivo and in silico, have been developed in the context of isolated microbes. In vivo microbial consortia development is extremely difficult and costly because it implies replicating suitable environments in the wet-lab. Computational approaches are thus a good, cost-effective alternative to study microbial communities, mainly via descriptive modelling, but also via engineering modelling. In this review we provide a detailed compilation of examples of engineered microbial communities and a comprehensive, historical revision of available computational metabolic modelling methods to better understand, and rationally engineer wild and synthetic microbial communities.

An insight to flux-balance analysis for biochemical networks

Systems biology is one of the integrated ways to study biological systems and is more favourable than the earlier used approaches. It includes metabolic pathway analysis, modelling, and regulatory as well as signal transduction for getting insights into cellular behaviour. Among the various techniques of modelling, simulation, analysis of networks and pathways, flux-based analysis (FBA) has been recognised because of its extensibility as well as simplicity. It is widely accepted because it is not like a mechanistic simulation which depends on accurate kinetic data. The study of fluxes through the network is informative and can give insights even in the absence of kinetic data. FBA is one of the widely used tools to study biochemical networks and needs information of reaction stoichiometry, growth requirements, specific measurement parameters of the biological system, in particular the reconstruction of the metabolic network for the genome-scale, many of which have already been built previously. It defines the boundaries of flux distributions which are possible and achievable with a defined set of genes. This review article gives an insight into FBA, from the extension of flux balancing to mathematical representation followed by a discussion about the formulation of flux-balance analysis problems, defining constraints for the stoichiometry of the pathways and the tools that can be used in FBA such as FASIMA, COBRA toolbox, and OptFlux. It also includes broader areas in terms of applications which can be covered by FBA as well as the queries which can be addressed through FBA.

Dynamic co-culture metabolic models reveal the fermentation dynamics, metabolic capacities and interplays of cheese starter cultures

In this study, we have investigated the cheese starter culture as a microbial community through a question: can the metabolic behaviour of a co-culture be explained by the characterized individual organism that constituted the co-culture? To address this question, the dairy-origin lactic acid bacteria Lactococcus lactis subsp. cremoris, Lactococcus lactis subsp. lactis, Streptococcus thermophilus and Leuconostoc mesenteroides, commonly used in cheese starter cultures, were grown in pure and four different co-cultures. We used a dynamic metabolic modelling approach based on the integration of the genome-scale metabolic networks of the involved organisms to simulate the co-cultures. The strain-specific kinetic parameters of dynamic models were estimated using the pure culture experiments and they were subsequently applied to co-culture models. Biomass, carbon source, lactic acid and most of the amino acid concentration profiles simulated by the co-culture models fit closely to the experimental results and the co-culture models explained the mechanisms behind the dynamic microbial abundance. We then applied the co-culture models to estimate further information on the co-cultures that could not be obtained by the experimental method used. This includes estimation of the profile of various metabolites in the co-culture medium such as flavour compounds produced and the individual organism level metabolic exchange flux profiles, which revealed the potential metabolic interactions between organisms in the co-cultures.

Minimizing the number of optimizations for efficient community dynamic flux balance analysis

Dynamic flux balance analysis uses a quasi-steady state assumption to calculate an organism's metabolic activity at each time-step of a dynamic simulation, using the well-known technique of flux balance analysis. For microbial communities, this calculation is especially costly and involves solving a linear constrained optimization problem for each member of the community at each time step. However, this is unnecessary and inefficient, as prior solutions can be used to inform future time steps. Here, we show that a basis for the space of internal fluxes can be chosen for each microbe in a community and this basis can be used to simulate forward by solving a relatively inexpensive system of linear equations at most time steps. We can use this solution as long as the resulting metabolic activity remains within the optimization problem's constraints (i.e. the solution to the linear system of equations remains a feasible to the linear program). As the solution becomes infeasible, it first becomes a feasible but degenerate solution to the optimization problem, and we can solve a different but related optimization problem to choose an appropriate basis to continue forward simulation. We demonstrate the efficiency and robustness of our method by comparing with currently used methods on a four species community, and show that our method requires at least 91% fewer optimizations to be solved. For reproducibility, we prototyped the method using Python. Source code is available at https://github.com/jdbrunner/surfin_fba.

Bridging substrate intake kinetics and bacterial growth phenotypes with flux balance analysis incorporating proteome allocation

Empirical kinetic models such as the Monod equation have been widely applied to relate the cell growth with substrate availability. The Monod equation shares a similar form with the mechanistically-based Michaelis-Menten kinetics for enzymatic processes, which has provoked long-standing and un-concluded conjectures on their relationship. In this work, we integrated proteome allocation principles into a Flux Balance Analysis (FBA) model of Escherichia coli, which quantitatively revealed potential mechanisms that underpin the phenomenological Monod parameters: the maximum specific growth rate could be dictated by the abundance of growth-controlling proteome and growth-pertinent proteome cost; more importantly, the Monod constant (K_s) was shown to relate to the Michaelis constant for substrate transport (K_m,g), with the link being dependent on the cell's metabolic strategy. Besides, the proposed model was able to predict glucose uptake rate at given external glucose concentration through the size of available proteome resource for substrate transport and its enzymatic cost, while growth rate and acetate overflow were accurately simulated for two E. coli strains. Bridging the enzymatic kinetics of substrate intake and overall growth phenotypes, this work offers a mechanistic interpretation to the empirical Monod law, and demonstrates the potential of coupling local and global cellular constrains in predictive modelling.

Metabolic pathway analysis for in silico design of efficient autotrophic production of advanced biofuels

Herein, autotrophic metabolism of Cupriavidus necator H16 growing on CO2, H2 and O2 gas mixture was analyzed by metabolic pathway analysis tools, specifically elementary mode analysis (EMA) and flux balance analysis (FBA). As case studies, recombinant strains of C. necator H16 for the production of short-chain (isobutanol) and long-chain (hexadecanol) alcohols were constructed and examined by a combined tools of EMA and FBA to comprehensively identify the cell’s metabolic flux profiles and its phenotypic spaces for the autotrophic production of recombinant products. The effect of genetic perturbations via gene deletion and overexpression on phenotypic space of the organism was simulated to improve strain performance for efficient bioconversion of CO2 to products at high yield and high productivity. EMA identified multiple gene deletion together with controlling gas input composition to limit phenotypic space and push metabolic fluxes towards high product yield, while FBA identified target gene overexpression to debottleneck rate-limiting fluxes, hence pulling more fluxes to enhance production rate of the products. A combination of gene deletion and overexpression resulted in designed mutant strains with a predicted yield of 0.21–0.42 g/g for isobutanol and 0.20–0.34 g/g for hexadecanol from CO2. The in silico-designed mutants were also predicted to show high productivity of up to 38.4 mmol/cell-h for isobutanol and 9.1 mmol/cell-h for hexadecanol under autotrophic cultivation. The metabolic modeling and analysis presented in this study could potentially serve as a valuable guidance for future metabolic engineering of C. necator H16 for an efficient CO2-to-biofuels conversion.

Nature's recyclers: anaerobic microbial communities drive crude biomass deconstruction

Microbial communities within anaerobic ecosystems have evolved to degrade and recycle carbon throughout the earth. A number of strains have been isolated from anaerobic microbial communities, which are rich in carbohydrate active enzymes (CAZymes) to liberate fermentable sugars from crude plant biomass (lignocellulose). However, natural anaerobic communities host a wealth of microbial diversity that has yet to be harnessed for biotechnological applications to hydrolyze crude biomass into sugars and value-added products. This review highlights recent advances in 'omics' techniques to sequence anaerobic microbial genomes, decipher microbial membership, and characterize CAZyme diversity in anaerobic microbiomes. With a focus on the herbivore rumen, we further discuss methods to discover new CAZymes, including those found within multi-enzyme fungal cellulosomes. Emerging techniques to characterize the interwoven metabolism and spatial interactions between anaerobes are also reviewed, which will prove critical to developing a predictive understanding of anaerobic communities to guide in microbiome engineering.

Systematic assessment of secondary bile acid metabolism in gut microbes reveals distinct metabolic capabilities in inflammatory bowel disease

BACKGROUND

The human gut microbiome performs important functions in human health and disease. A classic example for host-gut microbial co-metabolism is host biosynthesis of primary bile acids and their subsequent deconjugation and transformation by the gut microbiome. To understand these system-level host-microbe interactions, a mechanistic, multi-scale computational systems biology approach that integrates the different types of omic data is needed. Here, we use a systematic workflow to computationally model bile acid metabolism in gut microbes and microbial communities.

RESULTS

Therefore, we first performed a comparative genomic analysis of bile acid deconjugation and biotransformation pathways in 693 human gut microbial genomes and expanded 232 curated genome-scale microbial metabolic reconstructions with the corresponding reactions (available at https://vmh.life ). We then predicted the bile acid biotransformation potential of each microbe and in combination with other microbes. We found that each microbe could produce maximally six of the 13 secondary bile acids in silico, while microbial pairs could produce up to 12 bile acids, suggesting bile acid biotransformation being a microbial community task. To investigate the metabolic potential of a given microbiome, publicly available metagenomics data from healthy Western individuals, as well as inflammatory bowel disease patients and healthy controls, were mapped onto the genomes of the reconstructed strains. We constructed for each individual a large-scale personalized microbial community model that takes into account strain-level abundances. Using flux balance analysis, we found considerable variation in the potential to deconjugate and transform primary bile acids between the gut microbiomes of healthy individuals. Moreover, the microbiomes of pediatric inflammatory bowel disease patients were significantly depleted in their bile acid production potential compared with that of controls. The contributions of each strain to overall bile acid production potential across individuals were found to be distinct between inflammatory bowel disease patients and controls. Finally, bottlenecks limiting secondary bile acid production potential were identified in each microbiome model.

CONCLUSIONS

This large-scale modeling approach provides a novel way of analyzing metagenomics data to accelerate our understanding of the metabolic interactions between the host and gut microbiomes in health and diseases states. Our models and tools are freely available to the scientific community.

Synthetic Biology Tools to Engineer Microbial Communities for Biotechnology

Microbial consortia have been used in biotechnology processes, including fermentation, waste treatment, and agriculture, for millennia. Today, synthetic biologists are increasingly engineering microbial consortia for diverse applications, including the bioproduction of medicines, biofuels, and biomaterials from inexpensive carbon sources. An improved understanding of natural microbial ecosystems, and the development of new tools to construct synthetic consortia and program their behaviors, will vastly expand the functions that can be performed by communities of interacting microorganisms. Here, we review recent advancements in synthetic biology tools and approaches to engineer synthetic microbial consortia, discuss ongoing and emerging efforts to apply consortia for various biotechnological applications, and suggest future applications.

Model Microbial Consortia as Tools for Understanding Complex Microbial Communities

A major biological challenge in the postgenomic era has been untangling the composition and functions of microbes that inhabit complex communities or microbiomes. Multi-omics and modern bioinformatics have provided the tools to assay molecules across different cellular and community scales; however, mechanistic knowledge over microbial interactions often remains elusive. This is due to the immense diversity and the essentially undiminished volume of not-yet-cultured microbes. Simplified model communities hold some promise in enabling researchers to manage complexity so that they can mechanistically understand the emergent properties of microbial community interactions. In this review, we surveyed several approaches that have effectively used tractable model consortia to elucidate the complex behavior of microbial communities. We go further to provide some perspectives on the limitations and new opportunities with these approaches and highlight where these efforts are likely to lead as advances are made in molecular ecology and systems biology.

Machine Learning Reveals Missing Edges and Putative Interaction Mechanisms in Microbial Ecosystem Networks

Microbes affect each other's growth in multiple, often elusive, ways. The ensuing interdependencies form complex networks, believed to reflect taxonomic composition as well as community-level functional properties and dynamics. The elucidation of these networks is often pursued by measuring pairwise interactions in coculture experiments. However, the combinatorial complexity precludes an exhaustive experimental analysis of pairwise interactions, even for moderately sized microbial communities. Here, we used a machine learning random forest approach to address this challenge. In particular, we show how partial knowledge of a microbial interaction network, combined with trait-level representations of individual microbial species, can provide accurate inference of missing edges in the network and putative mechanisms underlying the interactions. We applied our algorithm to three case studies: an experimentally mapped network of interactions between auxotrophic Escherichia coli strains, a community of soil microbes, and a large in silico network of metabolic interdependencies between 100 human gut-associated bacteria. For this last case, 5% of the network was sufficient to predict the remaining 95% with 80% accuracy, and the mechanistic hypotheses produced by the algorithm accurately reflected known metabolic exchanges. Our approach, broadly applicable to any microbial or other ecological network, may drive the discovery of new interactions and new molecular mechanisms, both for therapeutic interventions involving natural communities and for the rational design of synthetic consortia. IMPORTANCE Different organisms in a microbial community may drastically affect each other's growth phenotypes, significantly affecting the community dynamics, with important implications for human and environmental health. Novel culturing methods and the decreasing costs of sequencing will gradually enable high-throughput measurements of pairwise interactions in systematic coculturing studies. However, a thorough characterization of all interactions that occur within a microbial community is greatly limited both by the combinatorial complexity of possible assortments and by the limited biological insight that interaction measurements typically provide without laborious specific follow-ups. Here, we show how a simple and flexible formal representation of microbial pairs can be used for the classification of interactions via machine learning. The approach we propose predicts with high accuracy the outcome of yet-to-be performed experiments and generates testable hypotheses about the mechanisms of specific interactions.

Metabolic Modeling of Microbial Community Interactions for Health, Environmental and Biotechnological Applications

In nature, microbes do not exist in isolation but co-exist in a variety of ecological and biological environments and on various host organisms. Due to their close proximity, these microbes interact among themselves, and also with the hosts in both positive and negative manners. Moreover, these interactions may modulate dynamically upon external stimulus as well as internal community changes. This demands systematic techniques such as mathematical modeling to understand the intrinsic community behavior. Here, we reviewed various approaches for metabolic modeling of microbial communities. If detailed species-specific information is available, segregated models of individual organisms can be constructed and connected via metabolite exchanges; otherwise, the community may be represented as a lumped network of metabolic reactions. The constructed models can then be simulated to help fill knowledge gaps, and generate testable hypotheses for designing new experiments. More importantly, such community models have been developed to study microbial interactions in various niches such as host microbiome, biogeochemical and bioremediation, waste water treatment and synthetic consortia. As such, the metabolic modeling efforts have allowed us to gain new insights into the natural and synthetic microbial communities, and design interventions to achieve specific goals. Finally, potential directions for future development in metabolic modeling of microbial communities were also discussed.

Genome-Scale, Constraint-Based Modeling of Nitrogen Oxide Fluxes during Coculture of Nitrosomonas europaea and Nitrobacter winogradskyi

Modern agriculture is sustained by application of inorganic nitrogen (N) fertilizer in the form of ammonium (NH₄⁺). Up to 60% of NH₄⁺-based fertilizer can be lost through leaching of nitrifier-derived nitrate (NO₃⁻), and through the emission of N oxide gases (i.e., nitric oxide [NO], N dioxide [NO₂], and nitrous oxide [N₂O] gases), the latter being a potent greenhouse gas. Our approach to modeling of nitrification suggests that both biotic and abiotic mechanisms function as important sources and sinks of N oxides during microaerobic conditions and that previous models might have underestimated gross NO production during nitrification.

Nitrification, the aerobic oxidation of ammonia to nitrate via nitrite, emits nitrogen (N) oxide gases (NO, NO₂, and N₂O), which are potentially hazardous compounds that contribute to global warming. To better understand the dynamics of nitrification-derived N oxide production, we conducted culturing experiments and used an integrative genome-scale, constraint-based approach to model N oxide gas sources and sinks during complete nitrification in an aerobic coculture of two model nitrifying bacteria, the ammonia-oxidizing bacterium Nitrosomonas europaea and the nitrite-oxidizing bacterium Nitrobacter winogradskyi. The model includes biotic genome-scale metabolic models (iFC578 and iFC579) for each nitrifier and abiotic N oxide reactions. Modeling suggested both biotic and abiotic reactions are important sources and sinks of N oxides, particularly under microaerobic conditions predicted to occur in coculture. In particular, integrative modeling suggested that previous models might have underestimated gross NO production during nitrification due to not taking into account its rapid oxidation in both aqueous and gas phases. The integrative model may be found at https://github.com/chaplenf/microBiome-v2.1.

IMPORTANCE Modern agriculture is sustained by application of inorganic nitrogen (N) fertilizer in the form of ammonium (NH₄⁺). Up to 60% of NH₄⁺-based fertilizer can be lost through leaching of nitrifier-derived nitrate (NO₃⁻), and through the emission of N oxide gases (i.e., nitric oxide [NO], N dioxide [NO₂], and nitrous oxide [N₂O] gases), the latter being a potent greenhouse gas. Our approach to modeling of nitrification suggests that both biotic and abiotic mechanisms function as important sources and sinks of N oxides during microaerobic conditions and that previous models might have underestimated gross NO production during nitrification.

Microbiota dysbiosis in inflammatory bowel diseases: in silico investigation of the oxygen hypothesis

BACKGROUND

Inflammatory bowel diseases (IBD), which include ulcerative colitis and Crohn's disease, cause chronic inflammation of the digestive tract in approximately 1.6 million Americans. A signature of IBD is dysbiosis of the gut microbiota marked by a significant reduction of obligate anaerobes and a sharp increase in facultative anaerobes. Numerous experimental studies have shown that IBD is strongly correlated with a decrease of Faecalibacterium prausnitzii and an increase of Escherichia coli. One hypothesis is that chronic inflammation induces increased oxygen levels in the gut, which in turn causes an imbalance between obligate and facultative anaerobes.

RESULTS

To computationally investigate the oxygen hypothesis, we developed a multispecies biofilm model based on genome-scale metabolic reconstructions of F. prausnitzii, E. coli and the common gut anaerobe Bacteroides thetaiotaomicron. Application of low bulk oxygen concentrations at the biofilm boundary reproduced experimentally observed behavior characterized by a sharp decrease of F. prausnitzii and a large increase of E. coli, demonstrating that dysbiosis consistent with IBD disease progression could be qualitatively predicted solely based on metabolic differences between the species. A diet with balanced carbohydrate and protein content was predicted to represent a metabolic "sweet spot" that increased the oxygen range over which F. prausnitzii could remain competitive and IBD could be sublimated. Host-microbiota feedback incorporated via a simple linear feedback between the average F. prausnitzii concentration and the bulk oxygen concentration did not substantially change the range of oxygen concentrations where dysbiosis was predicted, but the transition from normal species abundances to severe dysbiosis was much more dramatic and occurred over a much longer timescale. Similar predictions were obtained with sustained antibiotic treatment replacing a sustained oxygen perturbation, demonstrating how IBD might progress over several years with few noticeable effects and then suddenly produce severe disease symptoms.

CONCLUSIONS

The multispecies biofilm metabolic model predicted that oxygen concentrations of ∼1 micromolar within the gut could cause microbiota dysbiosis consistent with those observed experimentally for inflammatory bowel diseases. Our model predictions could be tested directly through the development of an appropriate in vitro system of the three species community and testing of microbiota-host interactions in gnotobiotic mice.

Constraint-based stoichiometric modelling from single organisms to microbial communities

Microbial communities are ubiquitously found in Nature and have direct implications for the environment, human health and biotechnology. The species composition and overall function of microbial communities are largely shaped by metabolic interactions such as competition for resources and cross-feeding. Although considerable scientific progress has been made towards mapping and modelling species-level metabolism, elucidating the metabolic exchanges between microorganisms and steering the community dynamics remain an enormous scientific challenge. In view of the complexity, computational models of microbial communities are essential to obtain systems-level understanding of ecosystem functioning. This review discusses the applications and limitations of constraint-based stoichiometric modelling tools, and in particular flux balance analysis (FBA). We explain this approach from first principles and identify the challenges one faces when extending it to communities, and discuss the approaches used in the field in view of these challenges. We distinguish between steady-state and dynamic FBA approaches extended to communities. We conclude that much progress has been made, but many of the challenges are still open.

Parametric studies of metabolic cooperativity in Escherichia coli colonies: Strain and geometric confinement effects

Characterizing the complex spatial and temporal interactions among cells in a biological system (i.e. bacterial colony, microbiome, tissue, etc.) remains a challenge. Metabolic cooperativity in these systems can arise due to the subtle interplay between microenvironmental conditions and the cells' regulatory machinery, often involving cascades of intra- and extracellular signalling molecules. In the simplest of cases, as demonstrated in a recent study of the model organism Escherichia coli, metabolic cross-feeding can arise in monoclonal colonies of bacteria driven merely by spatial heterogeneity in the availability of growth substrates; namely, acetate, glucose and oxygen. Another recent study demonstrated that even closely related E. coli strains evolved different glucose utilization and acetate production capabilities, hinting at the possibility of subtle differences in metabolic cooperativity and the resulting growth behavior of these organisms. Taking a first step towards understanding the complex spatio-temporal interactions within microbial populations, we performed a parametric study of E. coli growth on an agar substrate and probed the dependence of colony behavior on: 1) strain-specific metabolic characteristics, and 2) the geometry of the underlying substrate. To do so, we employed a recently developed multiscale technique named 3D dynamic flux balance analysis which couples reaction-diffusion simulations with iterative steady-state metabolic modeling. Key measures examined include colony growth rate and shape (height vs. width), metabolite production/consumption and concentration profiles, and the emergence of metabolic cooperativity and the fractions of cell phenotypes. Five closely related strains of E. coli, which exhibit large variation in glucose consumption and organic acid production potential, were studied. The onset of metabolic cooperativity was found to vary substantially between these five strains by up to 10 hours and the relative fraction of acetate utilizing cells within the colonies varied by a factor of two. Additionally, growth with six different geometries designed to mimic those that might be found in a laboratory, a microfluidic device, and inside a living organism were considered. Geometries were found to have complex, often nonlinear effects on colony growth and cross-feeding with "hard" features resulting in larger effect than "soft" features. These results demonstrate that strain-specific features and spatial constraints imposed by the growth substrate can have significant effects even for microbial populations as simple as isogenic E. coli colonies.

Compartmentalized metabolic network reconstruction of microbial communities to determine the effect of agricultural intervention on soils

Soil microbial communities are responsible for a wide range of ecological processes and have an important economic impact in agriculture. Determining the metabolic processes performed by microbial communities is crucial for understanding and managing ecosystem properties. Metagenomic approaches allow the elucidation of the main metabolic processes that determine the performance of microbial communities under different environmental conditions and perturbations. Here we present the first compartmentalized metabolic reconstruction at a metagenomics scale of a microbial ecosystem. This systematic approach conceives a meta-organism without boundaries between individual organisms and allows the in silico evaluation of the effect of agricultural intervention on soils at a metagenomics level. To characterize the microbial ecosystems, topological properties, taxonomic and metabolic profiles, as well as a Flux Balance Analysis (FBA) were considered. Furthermore, topological and optimization algorithms were implemented to carry out the curation of the models, to ensure the continuity of the fluxes between the metabolic pathways, and to confirm the metabolite exchange between subcellular compartments. The proposed models provide specific information about ecosystems that are generally overlooked in non-compartmentalized or non-curated networks, like the influence of transport reactions in the metabolic processes, especially the important effect on mitochondrial processes, as well as provide more accurate results of the fluxes used to optimize the metabolic processes within the microbial community.

Engineering microbial consortia for controllable outputs

Much research has been invested into engineering microorganisms to perform desired biotransformations; nonetheless, these efforts frequently fall short of expected results due to the unforeseen effects of biofeedback regulation and functional incompatibility. In nature, metabolic function is compartmentalized into diverse organisms assembled into robust consortia, in which the division of labor is thought to lead to increased community efficiency and productivity. Here we consider whether and how consortia can be designed to perform bioprocesses of interest beyond the metabolic flexibility limitations of a single organism. Advances in post-genomic analysis of microbial consortia and application of high-resolution global measurements now offer the promise of systems-level understanding of how microbial consortia adapt to changes in environmental variables and inputs of carbon and energy. We argue that, when combined with appropriate modeling frameworks, systems-level knowledge can markedly improve our ability to predict the fate and functioning of consortia. Here we articulate our collective perspective on the current and future state of microbial community engineering and control while placing specific emphasis on ecological principles that promote control over community function and emergent properties.

Genome-scale modelling of microbial metabolism with temporal and spatial resolution

Most natural microbial systems have evolved to function in environments with temporal and spatial variations. A major limitation to understanding such complex systems is the lack of mathematical modelling frameworks that connect the genomes of individual species and temporal and spatial variations in the environment to system behaviour. The goal of this review is to introduce the emerging field of spatiotemporal metabolic modelling based on genome-scale reconstructions of microbial metabolism. The extension of flux balance analysis (FBA) to account for both temporal and spatial variations in the environment is termed spatiotemporal FBA (SFBA). Following a brief overview of FBA and its established dynamic extension, the SFBA problem is introduced and recent progress is described. Three case studies are reviewed to illustrate the current state-of-the-art and possible future research directions are outlined. The author posits that SFBA is the next frontier for microbial metabolic modelling and a rapid increase in methods development and system applications is anticipated.