BacArena: Individual-based metabolic modeling of heterogeneous microbes in complex communities

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.

Examination of hydrogen cross-feeders using a colonic microbiota model

Background

Hydrogen cross-feeding microbes form a functionally important subset of the human colonic microbiota. The three major hydrogenotrophic functional groups of the colon: sulphate-reducing bacteria (SRB), methanogens and reductive acetogens, have been linked to wide ranging impacts on host physiology, health and wellbeing.

Results

An existing mathematical model for microbial community growth and metabolism was combined with models for each of the three hydrogenotrophic functional groups. The model was further developed for application to the colonic environment via inclusion of responsive pH, host metabolite absorption and the inclusion of host mucins. Predictions of the model, using two existing metabolic parameter sets, were compared to experimental faecal culture datasets. Model accuracy varied between experiments and measured variables and was most successful in predicting the growth of high relative abundance functional groups, such as the Bacteroides, and short chain fatty acid (SCFA) production. Two versions of the colonic model were developed: one representing the colon with sequential compartments and one utilising a continuous spatial representation. When applied to the colonic environment, the model predicted pH dynamics within the ranges measured in vivo and SCFA ratios comparable to those in the literature. The continuous version of the model simulated relative abundances of microbial functional groups comparable to measured values, but predictions were sensitive to the metabolic parameter values used for each functional group. Sulphate availability was found to strongly influence hydrogenotroph activity in the continuous version of the model, correlating positively with SRB and sulphide concentration and negatively with methanogen concentration, but had no effect in the compartmentalised model version.

Conclusions

Although the model predictions compared well to only some experimental measurements, the important features of the colon environment included make it a novel and useful contribution to modelling the colonic microbiota.

Synthetic gut microbiome: Advances and challenges

An exponential rise in studies regarding the association among human gut microbial communities, human health, and diseases is currently attracting the attention of researchers to focus on human gut microbiome research. However, even with the ever-growing number of studies on the human gut microbiome, translation into improved health is progressing slowly. This hampering is due to the complexities of the human gut microbiome, which is composed of >1,000 species of microorganisms, such as bacteria, archaea, viruses, and fungi. To overcome this complexity, it is necessary to reduce the gut microbiome, which can help simplify experimental variables to an extent, such that they can be deliberately manipulated and controlled. Reconstruction of synthetic or established gut microbial communities would make it easier to understand the structure, stability, and functional activities of the complex microbial community of the human gut. Here, we provide an overview of the developments and challenges of the synthetic human gut microbiome, and propose the incorporation of multi-omics and mathematical methods in a better synthetic gut ecosystem design, for easy translation of microbiome information to therapies.

Metabolic modeling predicts specific gut bacteria as key determinants for Candida albicans colonization levels

Candida albicans is a leading cause of life-threatening hospital-acquired infections and can lead to Candidemia with sepsis-like symptoms and high mortality rates. We reconstructed a genome-scale C. albicans metabolic model to investigate bacterial-fungal metabolic interactions in the gut as determinants of fungal abundance. We optimized the predictive capacity of our model using wild type and mutant C. albicans growth data and used it for in silico metabolic interaction predictions. Our analysis of more than 900 paired fungal–bacterial metabolic models predicted key gut bacterial species modulating C. albicans colonization levels. Among the studied microbes, Alistipes putredinis was predicted to negatively affect C. albicans levels. We confirmed these findings by metagenomic sequencing of stool samples from 24 human subjects and by fungal growth experiments in bacterial spent media. Furthermore, our pairwise simulations guided us to specific metabolites with promoting or inhibitory effect to the fungus when exposed in defined media under carbon and nitrogen limitation. Our study demonstrates that in silico metabolic prediction can lead to the identification of gut microbiome features that can significantly affect potentially harmful levels of C. albicans.

Genome-scale model reconstruction of C. albicans with 89% growth accuracy.

Mutualism and parasitism are the most common predicted C. albicans-gut bacteria interactions.

Metagenomic sequencing and in vitro assays reveal modulators of fungal growth.

Alistipes putredinis potentially prevents elevated C. albicans levels.

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.

Predicting Nash equilibria for microbial metabolic interactions

MOTIVATION

Microbial metabolic interactions impact ecosystems, human health and biotechnology profoundly. However, their determination remains elusive, invoking an urgent need for predictive models seamlessly integrating metabolism with evolutionary principles that shape community interactions.

RESULTS

Inspired by the evolutionary game theory, we formulated a bi-level optimization framework termed NECom for which any feasible solutions are Nash equilibria of microbial community metabolic models with/without an outer-level (community) objective function. Distinct from discrete matrix games, NECom models the continuous interdependent strategy space of metabolic fluxes. We showed that NECom successfully predicted several classical games in the context of metabolic interactions that were falsely or incompletely predicted by existing methods, including prisoner's dilemma, snowdrift and cooperation. The improved capability originates from the novel formulation to prevent 'forced altruism' hidden in previous static algorithms while allowing for sensing all potential metabolite exchanges to determine evolutionarily favorable interactions between members, a feature missing in dynamic methods. The results provided insights into why mutualism is favorable despite seemingly costly cross-feeding metabolites and demonstrated similarities and differences between games in the continuous metabolic flux space and matrix games. NECom was then applied to a reported algae-yeast co-culture system that shares typical cross-feeding features of lichen, a model system of mutualism. 488 growth conditions corresponding to 3,221 experimental data points were simulated. Without training any parameters using the data, NECom is more predictive of species' growth rates given uptake rates compared with flux balance analysis with an overall 63.5% and 81.7% reduction in root-mean-square error for the two species.

AVAILABILITY

Simulation code and data are available at https://github.com/Jingyi-Cai/NECom.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

Experimental and computational approaches to unravel microbial community assembly

Microbial communities have a preponderant role in the life support processes of our common home planet Earth. These extremely diverse communities drive global biogeochemical cycles, and develop intimate relationships with most multicellular organisms, with a significant impact on their fitness. Our understanding of their composition and function has enjoyed a significant thrust during the last decade thanks to the rise of high-throughput sequencing technologies. Intriguingly, the diversity patterns observed in nature point to the possible existence of fundamental community assembly rules. Unfortunately, these rules are still poorly understood, despite the fact that their knowledge could spur a scientific, technological, and economic revolution, impacting, for instance, agricultural, environmental, and health-related practices. In this minireview, I recapitulate the most important wet lab techniques and computational approaches currently employed in the study of microbial community assembly, and briefly discuss various experimental designs. Most of these approaches and considerations are also relevant to the study of microbial microevolution, as it has been shown that it can occur in ecological relevant timescales. Moreover, I provide a succinct review of various recent studies, chosen based on the diversity of ecological concepts addressed, experimental designs, and choice of wet lab and computational techniques. This piece aims to serve as a primer to those new to the field, as well as a source of new ideas to the more experienced researchers.

Inferring composition and function of the human gut microbiome in time and space: A review of genome-scale metabolic modelling tools

The human gut hosts a complex community of microorganisms that directly influences gastrointestinal physiology, playing a central role in human health. Because of its importance, the metabolic interplay between the gut microbiome and host metabolism has gained special interest. While there has been great progress in the field driven by metagenomics and experimental studies, the mechanisms underpinning microbial composition and interactions in the microbiome remain poorly understood. Genome-scale metabolic models are mathematical structures capable of describing the metabolic potential of microbial cells. They are thus suitable tools for probing the metabolic properties of microbial communities. In this review, we discuss the most recent and relevant genome-scale metabolic modelling tools for inferring the composition, interactions, and ultimately, biological function of the constituent species of a microbial community with special emphasis in the gut microbiota. Particular attention is given to constraint-based metabolic modelling methods as well as hybrid agent-based methods for capturing the interactions and behavior of the community in time and space. Finally, we discuss the challenges hindering comprehensive modelling of complex microbial communities and its application for the in-silico design of microbial consortia with therapeutic functions.

Genome scale metabolic models and analysis for evaluating probiotic potentials

Probiotics are live beneficial microorganisms that can be consumed in the form of dairy and food products as well as dietary supplements to promote a healthy balance of gut bacteria in humans. Practically, the main challenge is to identify and select promising strains and formulate multi-strain probiotic blends with consistent efficacy which is highly dependent on individual dietary regimes, gut environments, and health conditions. Limitations of current in vivo and in vitro methods for testing probiotic strains can be overcome by in silico model guided systems biology approaches where genome scale metabolic models (GEMs) can be used to describe their cellular behaviors and metabolic states of probiotic strains under various gut environments. Here, we summarize currently available GEMs of microbial strains with probiotic potentials and propose a knowledge-based framework to evaluate metabolic capabilities on the basis of six probiotic criteria. They include metabolic characteristics, stability, safety, colonization, postbiotics, and interaction with the gut microbiome which can be assessed by in silico approaches. As such, the most suitable strains can be identified to design personalized multi-strain probiotics in the future.

Defining the nutritional input for genome-scale metabolic models: A roadmap

The reconstruction and application of genome-scale metabolic network models is a central topic in the field of systems biology with numerous applications in biotechnology, ecology, and medicine. However, there is no agreed upon standard for the definition of the nutritional environment for these models. The objective of this article is to provide a guideline and a clear paradigm on how to translate nutritional information into an in-silico representation of the chemical environment. Step-by-step procedures explain how to characterise and categorise the nutritional input and to successfully apply it to constraint-based metabolic models. In parallel, we illustrate the proposed procedure with a case study of the growth of Escherichia coli in a complex nutritional medium and show that an accurate representation of the medium is crucial for physiological predictions. The proposed framework will assist researchers to expand their existing metabolic models of their microbial systems of interest with detailed representations of the nutritional environment, which allows more accurate and reproducible predictions of microbial metabolic processes.

A guide to deciphering microbial interactions and metabolic fluxes in microbiome communities

Microbiomes occupy nearly all environments on Earth. These communities of interacting microorganisms are highly complex, dynamic biological systems that impact and reshape the molecular composition of their habitats by performing complex biochemical transformations. The structure and function of microbiomes are influenced by local environmental stimuli and spatiotemporal changes. In order to control the dynamics and ultimately the function of microbiomes, we need to develop a mechanistic and quantitative understanding of the ecological, molecular, and evolutionary driving forces that govern these systems. Here, we describe recent advances in developing computational and experimental approaches that can promote a more fundamental understanding of microbial communities through comprehensive model-based analysis of heterogeneous data types across multiple scales, from intracellular metabolism, to metabolite cross-feeding interactions, to the emergent macroscopic behaviors. Ultimately, harnessing the full potential of microbiomes for practical applications will require developing new predictive modeling approaches and better tools to manipulate microbiome interactions.

Following the community development of SIHUMIx - a new intestinal in vitro model for bioreactor use

Diverse intestinal microbiota is frequently used in in vitro bioreactor models to study the effects of diet, chemical contaminations, or medication. However, the reproducible cultivation of fecal microbiota is challenging and the resultant communities behave highly dynamic. To approach the issue of reproducibility in in vitro models, we established an intestinal microbiota model community of reduced complexity, SIHUMIx, as a valuable model for in vitro use. The development of the SIHUMIx community was monitored over time with methods covering the cellular and the molecular level. We used microbial flow cytometry, intact protein profiling and terminal restriction fragment length polymorphism analysis to assess community structure. In parallel, we analyzed the functional level by targeted analysis of short-chain fatty acids and untargeted metabolomics. The stability properties constancy, resistance, and resilience were approached both on the structural and functional level of the community. We show that the SIHUMIx community is highly reproducible and constant since day 5 of cultivation. Furthermore, SIHUMIx has the ability to resist and recover from a pulsed perturbation, with changes in community structure recovered earlier than functional changes. Since community structure and function changed divergently, both levels need to be monitored at the same time to gain a full overview of the community development. All five methods are highly suitable to follow the community dynamics of SIHUMIx and indicated stability on day five. This makes SIHUMIx a suitable in vitro model to investigate the effects of e.g. medical, chemical, or dietary interventions.

From bag-of-genes to bag-of-genomes: metabolic modelling of communities in the era of metagenome-assembled genomes

Metagenomic sequencing of complete microbial communities has greatly enhanced our understanding of the taxonomic composition of microbiotas. This has led to breakthrough developments in bioinformatic disciplines such as assembly, gene clustering, metagenomic binning of species genomes and the discovery of an incredible, so far undiscovered, taxonomic diversity. However, functional annotations and estimating metabolic processes from single species - or communities - is still challenging. Earlier approaches relied mostly on inferring the presence of key enzymes for metabolic pathways in the whole metagenome, ignoring the genomic context of such enzymes, resulting in the 'bag-of-genes' approach to estimate functional capacities of microbiotas. Here, we review recent developments in metagenomic bioinformatics, with a special focus on emerging technologies to simulate and estimate metabolic information, that can be derived from metagenomic assembled genomes. Genome-scale metabolic models can be used to model the emergent properties of microbial consortia and whole communities, and the progress in this area is reviewed. While this subfield of metagenomics is still in its infancy, it is becoming evident that there is a dire need for further bioinformatic tools to address the complex combinatorial problems in modelling the metabolism of large communities as a 'bag-of-genomes'.

Our natural "makeup" reveals more than it hides: Modeling the skin and its microbiome

Skin is our primary interface with the environment. A structurally and functionally complex organ that hosts a dynamic ecosystem of microbes, and synthesizes many compounds that affect our well-being and psychosocial interactions. It is a natural platform of signal exchange between internal organs, skin resident microbes, and the environment. These interactions have gained a great deal of attention due to the increased prevalence of atopic diseases, and the co-occurrence of multiple allergic diseases related to allergic sensitization in early life. Despite significant advances in experimentally characterizing the skin, its microbial ecology, and disease phenotypes, high-levels of variability in these characteristics even for the same clinical phenotype are observed. Addressing this variability and resolving the relevant biological processes requires a systems approach. This review presents some of our current understanding of the skin, skin-immune, skin-neuroendocrine, skin-microbiome interactions, and computer-based modeling approaches to simulate this ecosystem in the context of health and disease. The review highlights the need for a systems-based understanding of this sophisticated ecosystem. This article is categorized under: Infectious Diseases > Computational Models.

μBialSim: Constraint-Based Dynamic Simulation of Complex Microbiomes

Microbial communities are pervasive in the natural environment, associated with many hosts, and of increasing importance in biotechnological applications. The complexity of these microbial systems makes the underlying mechanisms driving their dynamics difficult to identify. While experimental meta-OMICS techniques are routinely applied to record the inventory and activity of microbiomes over time, it remains difficult to obtain quantitative predictions based on such data. Mechanistic, quantitative mathematical modeling approaches hold the promise to both provide predictive power and shed light on cause-effect relationships driving these dynamic systems. We introduce μbialSim (pronounced "microbial sim"), a dynamic Flux-Balance-Analysis-based (dFBA) numerical simulator which is able to predict the time course in terms of composition and activity of microbiomes containing 100s of species in batch or chemostat mode. Activity of individual species is simulated by using separate FBA models which have access to a common pool of compounds, allowing for metabolite exchange. A novel augmented forward Euler method ensures numerical accuracy by temporarily reducing the time step size when compound concentrations decrease rapidly due to high compound affinities and/or the presence of many consuming species. We present three exemplary applications of μbialSim: a batch culture of a hydrogenotrophic archaeon, a syntrophic methanogenic biculture, and a 773-species human gut microbiome which exhibits a complex and dynamic pattern of metabolite exchange. Focusing on metabolite exchange as the main interaction type, μbialSim allows for the mechanistic simulation of microbiomes at their natural complexity. Simulated trajectories can be used to contextualize experimental meta-OMICS data and to derive hypotheses on cause-effect relationships driving community dynamics based on scenario simulations. μbialSim is implemented in Matlab and relies on the COBRA Toolbox or CellNetAnalyzer for FBA calculations. The source code is available under the GNU General Public License v3.0 at https://git.ufz.de/UMBSysBio/microbialsim.

Competition for Hydrogen Prevents Coexistence of Human Gastrointestinal Hydrogenotrophs in Continuous Culture

Understanding the metabolic dynamics of the human gastrointestinal tract (GIT) microbiota is of growing importance as research continues to link the microbiome to host health status. Microbial strains that metabolize hydrogen have been associated with a variety of both positive and negative host nutritional and health outcomes, but limited data exists for their competition in the GIT. To enable greater insight into the behaviour of these microbes, a mathematical model was developed for the metabolism and growth of the three major hydrogenotrophic groups: sulphate-reducing bacteria (SRB), methanogens and reductive acetogens. In batch culture simulations with abundant sulphate and hydrogen, the SRB outcompeted the methanogen for hydrogen due to having a half-saturation constant 10⁶ times lower than that of the methanogen. The acetogen, with a high model threshold for hydrogen uptake of around 70 mM, was the least competitive. Under high lactate and zero sulphate conditions, hydrogen exchange between the SRB and the methanogen was the dominant interaction. The methanogen grew at 70% the rate of the SRB, with negligible acetogen growth. In continuous culture simulations, both the SRB and the methanogen were washed out at dilution rates above 0.15 h⁻¹ regardless of substrate availability, whereas the acetogen could survive under abundant hydrogen conditions. Specific combinations of conditions were required for survival of more than one hydrogenotroph in continuous culture, and survival of all three was not possible. The stringency of these requirements and the inability of the model to simulate survival of all three hydrogenotrophs in continuous culture demonstrates that factors outside of those modelled are vital to allow hydrogenotroph coexistence in the GIT.

ACBM: An Integrated Agent and Constraint Based Modeling Framework for Simulation of Microbial Communities

The development of new methods capable of more realistic modeling of microbial communities necessitates that their results be quantitatively comparable with experimental findings. In this research, a new integrated agent and constraint based modeling framework abbreviated ACBM has been proposed that integrates agent-based and constraint-based modeling approaches. ACBM models the cell population in three-dimensional space to predict spatial and temporal dynamics and metabolic interactions. When used to simulate the batch growth of C. beijerinckii and two-species communities of F. prausnitzii and B. adolescent., ACBM improved on predictions made by two previous models. Furthermore, when transcriptomic data were integrated with a metabolic model of E. coli to consider intracellular constraints in the metabolism, ACBM accurately predicted growth rate, half-rate constant, and concentration of biomass, glucose, and acidic products over time. The results also show that the framework was able to predict the metabolism changes in the early stationary compared to the log phase. Finally, ACBM was implemented to estimate starved cells under heterogeneous feeding and it was concluded that a percentage of cells are always subject to starvation in a bioreactor with high volume.

The Computational Diet: A Review of Computational Methods Across Diet, Microbiome, and Health

Food and human health are inextricably linked. As such, revolutionary impacts on health have been derived from advances in the production and distribution of food relating to food safety and fortification with micronutrients. During the past two decades, it has become apparent that the human microbiome has the potential to modulate health, including in ways that may be related to diet and the composition of specific foods. Despite the excitement and potential surrounding this area, the complexity of the gut microbiome, the chemical composition of food, and their interplay in situ remains a daunting task to fully understand. However, recent advances in high-throughput sequencing, metabolomics profiling, compositional analysis of food, and the emergence of electronic health records provide new sources of data that can contribute to addressing this challenge. Computational science will play an essential role in this effort as it will provide the foundation to integrate these data layers and derive insights capable of revealing and understanding the complex interactions between diet, gut microbiome, and health. Here, we review the current knowledge on diet-health-gut microbiota, relevant data sources, bioinformatics tools, machine learning capabilities, as well as the intellectual property and legislative regulatory landscape. We provide guidance on employing machine learning and data analytics, identify gaps in current methods, and describe new scenarios to be unlocked in the next few years in the context of current knowledge.

The engineering of spatially linked microbial consortia - potential and perspectives

Traditional biotechnological applications of microorganisms employ mono-cultivation or co-cultivation in well-mixed vessels disregarding the potential of spatially organized cultures. Metabolic specialization and guided species interactions facilitated through spatial isolation would enable consortia of microbes to accomplish more complex functions than currently possible, for bioproduction as well as biodegradation processes. Here, we review concepts of spatially linked microbial consortia in which spatial arrangement is optimized to increase control and facilitate new species combinations. We highlight that genome-scale metabolic network models can inform the design and tuning of synthetic microbial consortia and suggest that a standardized assembly of such systems allows the combination of 'incompatibles', potentially leading to countless novel applications.

The microbial exometabolome: ecological resource and architect of microbial communities

All microorganisms release many metabolites, collectively known as the exometabolome. The resultant multi-way cross-feeding of metabolites among microorganisms distributes resources, thereby increasing total biomass of the microbial community, and promotes the recruitment and persistence of phylogenetically and functionally diverse taxa in microbial communities. Metabolite transfer can also select for evolutionary diversification, yielding multiple closely related but functionally distinct strains. Depending on starting conditions, the evolved strains may be auxotrophs requiring metabolic outputs from producer cells or, alternatively, display loss of complementary reactions in metabolic pathways, with increased metabolic efficiency. Metabolite cross-feeding is widespread in many microbial communities associated with animals and plants, including the animal gut microbiome, and these metabolic interactions can yield products valuable to the host. However, metabolite exchange between pairs of intracellular microbial taxa that share the same host cell or organ can be very limited compared to pairs of free-living microorganisms, perhaps as a consequence of host controls over the metabolic function of intracellular microorganisms. Priorities for future research include the development of tools for improved quantification of metabolite exchange in complex communities and greater integration of the roles of metabolic cross-feeding and other ecological processes, including priority effects and antagonistic interactions, in shaping microbial communities. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Computational Modeling of the Human Microbiome

The impact of microorganisms on human health has long been acknowledged and studied, but recent advances in research methodologies have enabled a new systems-level perspective on the collections of microorganisms associated with humans, the human microbiome. Large-scale collaborative efforts such as the NIH Human Microbiome Project have sought to kick-start research on the human microbiome by providing foundational information on microbial composition based upon specific sites across the human body. Here, we focus on the four main anatomical sites of the human microbiome: gut, oral, skin, and vaginal, and provide information on site-specific background, experimental data, and computational modeling. Each of the site-specific microbiomes has unique organisms and phenomena associated with them; there are also high-level commonalities. By providing an overview of different human microbiome sites, we hope to provide a perspective where detailed, site-specific research is needed to understand causal phenomena that impact human health, but there is equally a need for more generalized methodology improvements that would benefit all human microbiome research.

The functional repertoire contained within the native microbiota of the model nematode Caenorhabditis elegans

The microbiota is generally assumed to have a substantial influence on the biology of multicellular organisms. The exact functional contributions of the microbes are often unclear and cannot be inferred easily from 16S rRNA genotyping, which is commonly used for taxonomic characterization of bacterial associates. In order to bridge this knowledge gap, we here analyzed the metabolic competences of the native microbiota of the model nematode Caenorhabditis elegans. We integrated whole-genome sequences of 77 bacterial microbiota members with metabolic modeling and experimental characterization of bacterial physiology. We found that, as a community, the microbiota can synthesize all essential nutrients for C. elegans. Both metabolic models and experimental analyses revealed that nutrient context can influence how bacteria interact within the microbiota. We identified key bacterial traits that are likely to influence the microbe's ability to colonize C. elegans (i.e., the ability of bacteria for pyruvate fermentation to acetoin) and affect nematode fitness (i.e., bacterial competence for hydroxyproline degradation). Considering that the microbiota is usually neglected in C. elegans research, the resource presented here will help our understanding of this nematode's biology in a more natural context. Our integrative approach moreover provides a novel, general framework to characterize microbiota-mediated functions.

Defining and Evaluating Microbial Contributions to Metabolite Variation in Microbiome-Metabolome Association Studies

Correlation-based analysis of paired microbiome-metabolome data sets is becoming a widespread research approach, aiming to comprehensively identify microbial drivers of metabolic variation. To date, however, the limitations of this approach and other microbiome-metabolome analysis methods have not been comprehensively evaluated. To address this challenge, we have introduced a mathematical framework to quantify the contribution of each taxon to metabolite variation based on uptake and secretion fluxes. We additionally used a multispecies metabolic model to simulate simplified gut communities, generating idealized microbiome-metabolome data sets. We then compared observed taxon-metabolite correlations in these data sets to calculated ground truth taxonomic contribution values. We found that in simulations of both a representative simple 10-species community and complex human gut microbiota, correlation-based analysis poorly identified key contributors, with an extremely low predictive value despite the idealized setting. We further demonstrate that the predictive value of correlation analysis is strongly influenced by both metabolite and taxon properties, as well as by exogenous environmental variation. We finally discuss the practical implications of our findings for interpreting microbiome-metabolome studies.IMPORTANCE Identifying the key microbial taxa responsible for metabolic differences between microbiomes is an important step toward understanding and manipulating microbiome metabolism. To achieve this goal, researchers commonly conduct microbiome-metabolome association studies, comprehensively measuring both the composition of species and the concentration of metabolites across a set of microbial community samples and then testing for correlations between microbes and metabolites. Here, we evaluated the utility of this general approach by first developing a rigorous mathematical definition of the contribution of each microbial taxon to metabolite variation and then examining these contributions in simulated data sets of microbial community metabolism. We found that standard correlation-based analysis of our simulated microbiome-metabolome data sets can identify true contributions with very low predictive value and that its performance depends strongly on specific properties of both metabolites and microbes, as well as on those of the surrounding environment. Combined, our findings can guide future interpretation and validation of microbiome-metabolome studies.

The Simplified Human Intestinal Microbiota (SIHUMIx) Shows High Structural and Functional Resistance against Changing Transit Times in In Vitro Bioreactors

Many functions in host-microbiota interactions are potentially influenced by intestinal transit times, but little is known about the effects of altered transition times on the composition and functionality of gut microbiota. To analyze these effects, we cultivated the model community SIHUMIx in bioreactors in order to determine the effects of varying transit times (TT) on the community structure and function. After five days of continuous cultivation, we investigated the influence of different medium TT of 12 h, 24 h, and 48 h. For profiling the microbial community, we applied flow cytometric fingerprinting and revealed changes in the community structure of SIHUMIx during the change of TT, which were not associated with changes in species abundances. For pinpointing metabolic alterations, we applied metaproteomics and metabolomics and found, along with shortening the TT, a slight decrease in glycan biosynthesis, carbohydrate, and amino acid metabolism and, furthermore, a reduction in butyrate, methyl butyrate, isobutyrate, valerate, and isovalerate concentrations. Specifically, B. thetaiotaomicron was identified to be affected in terms of butyrate metabolism. However, communities could recover to the original state afterward. This study shows that SIHUMIx showed high structural stability when TT changed-even four-fold. Resistance values remained high, which suggests that TTs did not interfere with the structure of the community to a certain degree.

Developing a Microbial Consortium for Enhanced Metabolite Production from Simulated Food Waste

Food waste disposal and transportation of commodity chemicals to the point-of-need are substantial challenges in military environments. Here, we propose addressing these challenges via the design of a microbial consortium for the fermentation of food waste to hydrogen. First, we simulated the exchange metabolic fluxes of monocultures and pairwise co-cultures using genome-scale metabolic models on a food waste proxy. We identified that one of the top hydrogen producing co-cultures comprised Clostridium beijerinckii NCIMB 8052 and Yokenella regensburgei ATCC 43003. A consortium of these two strains produced a similar amount of hydrogen gas and increased butyrate compared to the C. beijerinckii monoculture, when grown on an artificial garbage slurry. Increased butyrate production in the consortium can be attributed to cross-feeding of lactate produced by Y. regensburgei. Moreover, exogenous lactate promotes the growth of C. beijerinckii with or without a limited amount of glucose. Increasing the scale of the consortium fermentation proved challenging, as two distinct attempts to scale-up the enhanced butyrate production resulted in different metabolic profiles than observed in smaller scale fermentations. Though the genome-scale metabolic model simulations provided a useful starting point for the design of microbial consortia to generate value-added products from waste materials, further model refinements based on experimental results are required for more robust predictions.

Systems biology perspective for studying the gut microbiota in human physiology and liver diseases

The advancement in high-throughput sequencing technologies and systems biology approaches have revolutionized our understanding of biological systems and opened a new path to investigate unacknowledged biological phenomena. In parallel, the field of human microbiome research has greatly evolved and the relative contribution of the gut microbiome to health and disease have been systematically explored. This review provides an overview of the network-based and translational systems biology-based studies focusing on the function and composition of gut microbiota. We also discussed the association between the gut microbiome and the overall human physiology, as well as hepatic diseases and other metabolic disorders.

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.

FLYCOP: metabolic modeling-based analysis and engineering microbial communities

Motivation

Synthetic microbial communities begin to be considered as promising multicellular biocatalysts having a large potential to replace engineered single strains in biotechnology applications, in pharmaceutical, chemical and living architecture sectors. In contrast to single strain engineering, the effective and high-throughput analysis and engineering of microbial consortia face the lack of knowledge, tools and well-defined workflows. This manuscript contributes to fill this important gap with a framework, called FLYCOP (FLexible sYnthetic Consortium OPtimization), which contributes to microbial consortia modeling and engineering, while improving the knowledge about how these communities work. FLYCOP selects the best consortium configuration to optimize a given goal, among multiple and diverse configurations, in a flexible way, taking temporal changes in metabolite concentrations into account.

Results

In contrast to previous systems optimizing microbial consortia, FLYCOP has novel characteristics to face up to new problems, to represent additional features and to analyze events influencing the consortia behavior. In this manuscript, FLYCOP optimizes a Synechococcus elongatus-Pseudomonas putida consortium to produce the maximum amount of bio-plastic (PHA, polyhydroxyalkanoate), and highlights the influence of metabolites exchange dynamics in a four auxotrophic Escherichia coli consortium with parallel growth. FLYCOP can also provide an explanation about biological evolution driving evolutionary engineering endeavors by describing why and how heterogeneous populations emerge from monoclonal ones.

Availability and implementation

Code reproducing the study cases described in this manuscript are available on-line: https://github.com/beatrizgj/FLYCOP.

Supplementary information

Supplementary data are available at Bioinformatics online.

Scalable and exhaustive screening of metabolic functions carried out by microbial consortia

Motivation

The selection of species exhibiting metabolic behaviors of interest is a challenging step when switching from the investigation of a large microbiota to the study of functions effectiveness. Approaches based on a compartmentalized framework are not scalable. The output of scalable approaches based on a non-compartmentalized modeling may be so large that it has neither been explored nor handled so far.

Results

We present the Miscoto tool to facilitate the selection of a community optimizing a desired function in a microbiome by reporting several possibilities which can be then sorted according to biological criteria. Communities are exhaustively identified using logical programming and by combining the non-compartmentalized and the compartmentalized frameworks. The benchmarking of 4.9 million metabolic functions associated with the Human Microbiome Project, shows that Miscoto is suited to screen and classify metabolic producibility in terms of feasibility, functional redundancy and cooperation processes involved. As an illustration of a host-microbial system, screening the Recon 2.2 human metabolism highlights the role of different consortia within a family of 773 intestinal bacteria.

Availability and implementation

Miscoto source code, instructions for use and examples are available at: https://github.com/cfrioux/miscoto.

Consumer Safety Considerations of Skin and Oral Microbiome Perturbation

Microbiomes associated with human skin and the oral cavity are uniquely exposed to personal care regimes. Changes in the composition and activities of the microbial communities in these environments can be utilized to promote consumer health benefits, for example, by reducing the numbers, composition, or activities of microbes implicated in conditions such as acne, axillary odor, dandruff, and oral diseases. It is, however, important to ensure that innovative approaches for microbiome manipulation do not unsafely disrupt the microbiome or compromise health, and where major changes in the composition or activities of the microbiome may occur, these require evaluation to ensure that critical biological functions are unaffected. This article is based on a 2-day workshop held at SEAC Unilever, Sharnbrook, United Kingdom, involving 31 specialists in microbial risk assessment, skin and oral microbiome research, microbial ecology, bioinformatics, mathematical modeling, and immunology. The first day focused on understanding the potential implications of skin and oral microbiome perturbation, while approaches to characterize those perturbations were discussed during the second day. This article discusses the factors that the panel recommends be considered for personal care products that target the microbiomes of the skin and the oral cavity.

Modelling approaches for studying the microbiome

Advances in metagenome sequencing of the human microbiome have provided a plethora of new insights and revealed a close association of this complex ecosystem with a range of human diseases. However, there is little knowledge about how the different members of the microbial community interact with each other and with the host, and we lack basic mechanistic understanding of these interactions related to health and disease. Mathematical modelling has been demonstrated to be highly advantageous for gaining insights into the dynamics and interactions of complex systems and in recent years, several modelling approaches have been proposed to enhance our understanding of the microbiome. Here, we review the latest developments and current approaches, and highlight how different modelling strategies have been applied to unravel the highly dynamic nature of the human microbiome. Furthermore, we discuss present limitations of different modelling strategies and provide a perspective of how modelling can advance understanding and offer new treatment routes to impact human health.

Predicting the Longitudinally and Radially Varying Gut Microbiota Composition Using Multi-Scale Microbial Metabolic Modeling

Background: The gut microbiota is a heterogeneous group of microbes that is spatially distributed along various sections of the intestines and across the mucosa and lumen in each section. Understanding the dynamics between the spatially differential microbial populations and the driving forces for the observed spatial organization will provide valuable insights into important questions such as the nature of colonization of the infant gut and different types of inflammatory bowel disease localized in different regions of the intestines. However, in most studies, the microbiota is sampled only at a single site (often feces) or from a particular anatomical site of the intestines. Differential oxygen availability is putatively a key factor shaping the spatial organization. Results: To test this hypothesis, we constructed a community genome-scale metabolic model consisting of representative organisms for the major phyla present in the human gut microbiome. By solving step-wise optimization problems embedded in a dynamic framework to predict community metabolism and integrate the mucosally-adherent with the luminal microbiome between consecutive sections along the intestines, we were able to capture (i) the essential features of the spatially differential composition of obligate anaerobes vs. facultative anaerobes and aerobes determined experimentally, and (ii) the accumulation of microbial biomass in the lumen. Sensitivity analysis suggests that the spatial organization depends primarily on the oxygen-per-microbe availability in each region. Oxygen availability is reduced relative to the ~100-fold increase in mucosal microbial density along the intestines, causing the switch between aerobes and anaerobes. Conclusion: The proposed integrated dynamic framework is able to predict spatially differential gut microbiota composition using microbial genome-scale metabolic models and test hypotheses regarding the dynamics of the gut microbiota. It can potentially become a valuable tool for exploring therapeutic strategies for site-specific perturbation of the gut microbiota and the associated metabolic activities.

Modeling metabolic networks of individual bacterial agents in heterogeneous and dynamic soil habitats (IndiMeSH)

Natural soil is characterized as a complex habitat with patchy hydrated islands and spatially variable nutrients that is in a constant state of change due to wetting-drying dynamics. Soil microbial activity is often concentrated in sparsely distributed hotspots that contribute disproportionally to macroscopic biogeochemical nutrient cycling and greenhouse gas emissions. The mechanistic representation of such dynamic hotspots requires new modeling approaches capable of representing the interplay between dynamic local conditions and the versatile microbial metabolic adaptations. We have developed IndiMeSH (Individual-based Metabolic network model for Soil Habitats) as a spatially explicit model for the physical and chemical microenvironments of soil, combined with an individual-based representation of bacterial motility and growth using adaptive metabolic networks. The model uses angular pore networks and a physically based description of the aqueous phase as a backbone for nutrient diffusion and bacterial dispersal combined with dynamic flux balance analysis to calculate growth rates depending on local nutrient conditions. To maximize computational efficiency, reduced scale metabolic networks are used for the simulation scenarios and evaluated strategically to the genome scale model. IndiMeSH was compared to a well-established population-based spatiotemporal metabolic network model (COMETS) and to experimental data of bacterial spatial organization in pore networks mimicking soil aggregates. IndiMeSH was then used to strategically study dynamic response of a bacterial community to abrupt environmental perturbations and the influence of habitat geometry and hydration conditions. Results illustrate that IndiMeSH is capable of representing trophic interactions among bacterial species, predicting the spatial organization and segregation of bacterial populations due to oxygen and carbon gradients, and provides insights into dynamic community responses as a consequence of environmental changes. The modular design of IndiMeSH and its implementation are adaptable, allowing it to represent a wide variety of experimental and in silico microbial systems.

Soil bacterial communities are key players in global biogeochemical cycles and drive other soil regulatory and provisional ecosystem functions. Despite the relatively high bacterial abundance found in fertile soil, bacteria occupy only a small fraction of the soil surfaces and often form hotspots with disproportionate contributions to observed biogenic fluxes. As soil opacity and complexity limit detailed observations of such hotspots in situ, we have developed a modeling platform, IndiMeSH (Individual-based Metabolic network model for Soil Habitats), to enable systematic study of dense multispecies bacterial communities within a structured habitat resembling (but not limited) to soil. The model is capable of representing multispecies trophic interactions and spatial self-organization in response to nutrient gradients, as confirmed in comparison with published results. IndiMeSH offers new opportunities for quantifying bacterial hotspot formation and dynamics and observe their resilience and response to perturbations in hydration and nutrient conditions.

Meta-Omics- and Metabolic Modeling-Assisted Deciphering of Human Microbiota Metabolism

The human microbiota is a complex community of commensal, symbiotic, and pathogenic microbes that play a crucial role in maintaining the homeostasis of human health. Such a homeostasis is maintained through the collective functioning of enzymatic genes responsible for the production of metabolites, enabling the interaction and signaling within microbiota as well as between microbes and the human host. Understanding microbial genes, their associated chemistries and functions would be valuable for engineering systemic metabolic pathways within the microbiota to manage human health and diseases. Given that there are many unknown gene metabolic functions and interactions, increasing efforts have been made to gain insights into the underlying functions of microbiota metabolism. This can be achieved through culture-independent metagenomic approaches and metabolic modeling to simulate the microenvironment of human microbiota. In this article, the recent advances in metagenome mining and functional profiling for the discovery of the genetic and biochemical links in human microbiota metabolism as well as metabolic modeling for simulation and prediction of metabolic fluxes in the human microbiota are reviewed. This review provides useful insights into the understanding, reconstruction, and modulation of the human microbiota guided by the knowledge acquired from the basic understanding of the human microbiota metabolism.

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.

Metastasis of pancreatic cancer: An uninflamed liver micromilieu controls cell growth and cancer stem cell properties by oxidative phosphorylation in pancreatic ductal epithelial cells

Pancreatic ductal adenocarcinoma (PDAC) is commonly diagnosed when liver metastases already emerged. We recently demonstrated that hepatic stromal cells determine the dormancy status along with cancer stem cell (CSC) properties of pancreatic ductal epithelial cells (PDECs) during metastasis. This study investigated the influence of the hepatic microenvironment - and its inflammatory status - on metabolic alterations and how these impact cell growth and CSC-characteristics of PDECs. Coculture with hepatic stellate cells (HSCs), simulating a physiological liver stroma, but not with hepatic myofibroblasts (HMFs) representing liver inflammation promoted expression of Succinate Dehydrogenase subunit B (SDHB) and an oxidative metabolism along with a quiescent phenotype in PDECs. SiRNA-mediated SDHB knockdown increased cell growth and CSC-properties. Moreover, liver micrometastases of tumor bearing KPC mice strongly expressed SDHB while expression of the CSC-marker Nestin was exclusively found in macrometastases. Consistently, RNA-sequencing and in silico modeling revealed significantly altered metabolic fluxes and enhanced SDH activity predominantly in premalignant PDECs in the presence of HSC compared to HMF. Overall, these data emphasize that the hepatic microenvironment determines the metabolism of disseminated PDECs thereby controlling cell growth and CSC-properties during liver metastasis.

Streptomyces Volatile Compounds Influence Exploration and Microbial Community Dynamics by Altering Iron Availability

Microbial growth and community interactions are influenced by a multitude of factors. A new mode of Streptomyces growth—exploration—is promoted by interactions with the yeast Saccharomyces cerevisiae and requires the emission of trimethylamine (TMA), a pH-raising volatile compound. We show here that TMA emission also profoundly alters the environment around exploring cultures. It specifically reduces iron availability, and this in turn adversely affects the viability of surrounding microbes. Paradoxically, Streptomyces bacteria thrive in these iron-depleted niches, both rewiring their gene expression and metabolism to facilitate iron uptake and increasing their exploration rate. Growth in close proximity to other microbes adept at iron uptake also enhances exploration. Collectively, the data from this work reveal a new role for bacterial volatile compounds in modulating nutrient availability and microbial community behavior. The results further expand the repertoire of interspecies interactions and nutrient cues that impact Streptomyces exploration and provide new mechanistic insight into this unique mode of bacterial growth.

Bacteria and fungi produce a wide array of volatile organic compounds (VOCs), and these can act as chemical cues or as competitive tools. Recent work has shown that the VOC trimethylamine (TMA) can promote a new form of Streptomyces growth, termed “exploration.” Here, we report that TMA also serves to alter nutrient availability in the area surrounding exploring cultures: TMA dramatically increases the environmental pH and, in doing so, reduces iron availability. This, in turn, compromises the growth of other soil bacteria and fungi. In response to this low-iron environment, Streptomyces venezuelae secretes a suite of differentially modified siderophores and upregulates genes associated with siderophore uptake. Further reducing iron levels by limiting siderophore uptake or growing cultures in the presence of iron chelators enhanced exploration. Exploration was also increased when S. venezuelae was grown in association with the related low-iron- and TMA-tolerant Amycolatopsis bacteria, due to competition for available iron. We are only beginning to appreciate the role of VOCs in natural communities. This work reveals a new role for VOCs in modulating iron levels in the environment and implies a critical role for VOCs in modulating the behavior of microbes and the makeup of their communities. It further adds a new dimension to our understanding of the interspecies interactions that influence Streptomyces exploration and highlights the importance of iron in exploration modulation.

Metabolic Modeling of Human Gut Microbiota on a Genome Scale: An Overview

There is growing interest in the metabolic interplay between the gut microbiome and host metabolism. Taxonomic and functional profiling of the gut microbiome by next-generation sequencing (NGS) has unveiled substantial richness and diversity. However, the mechanisms underlying interactions between diet, gut microbiome and host metabolism are still poorly understood. Genome-scale metabolic modeling (GSMM) is an emerging approach that has been increasingly applied to infer diet⁻microbiome, microbe⁻microbe and host⁻microbe interactions under physiological conditions. GSMM can, for example, be applied to estimate the metabolic capabilities of microbes in the gut. Here, we discuss how meta-omics datasets such as shotgun metagenomics, can be processed and integrated to develop large-scale, condition-specific, personalized microbiota models in healthy and disease states. Furthermore, we summarize various tools and resources available for metagenomic data processing and GSMM, highlighting the experimental approaches needed to validate the model predictions.

Understanding and Engineering Distributed Biochemical Pathways in Microbial Communities

Microbiomes impact nearly every environment on Earth by modulating the molecular composition of the environment. Temporally changing environmental stimuli and spatial organization are major variables shaping the structure and function of microbiomes. The web of interactions among members of these communities and between the organisms and the environment dictates microbiome functions. Microbial interactions are major drivers of microbiomes and are modulated by spatiotemporal parameters. A mechanistic and quantitative understanding of ecological, molecular, and environmental forces shaping microbiomes could inform strategies to control microbiome dynamics and functions. Major challenges for harnessing the potential of microbiomes for diverse applications include the development of predictive modeling frameworks and tools for precise manipulation of microbiome behaviors.

Suboptimal community growth mediated through metabolite crossfeeding promotes species diversity in the gut microbiota

The gut microbiota represent a highly complex ecosystem comprised of approximately 1000 species that forms a mutualistic relationship with the human host. A critical attribute of the microbiota is high species diversity, which provides system robustness through overlapping and redundant metabolic capabilities. The gradual loss of bacterial diversity has been associated with a broad array of gut pathologies and diseases including malnutrition, obesity, diabetes and inflammatory bowel disease. We formulated an in silico community model of the gut microbiota by combining genome-scale metabolic reconstructions of 28 representative species to explore the relationship between species diversity and community growth. While the individual species offered a broad range of metabolic capabilities, communities optimized for maximal growth on simulated Western and high-fiber diets had low diversities and imbalances in short-chain fatty acid (SCFA) synthesis characterized by acetate overproduction. Community flux variability analysis performed with the 28-species model and a reduced 20-species model suggested that enhanced species diversity and more balanced SCFA production were achievable at suboptimal growth rates. We developed a simple method for constraining species abundances to sample the growth-diversity tradeoff and used the 20-species model to show that tradeoff curves for Western and high-fiber diets resembled Pareto-optimal surfaces. Compared to maximal growth solutions, suboptimal growth solutions were characterized by higher species diversity, more balanced SCFA synthesis and lower exchange rates of crossfed metabolites between more species. We hypothesized that modulation of crossfeeding relationships through host-microbiota interactions could be an important means for maintaining species diversity and suggest that community metabolic modeling approaches that allow multiobjective optimization of growth and diversity are needed for more realistic simulation of complex communities.

Human gut microbe co-cultures have greater potential than monocultures for food waste remediation to commodity chemicals

Food waste represents an underutilized resource for commodity chemical generation. Constituents of the human gut microbiota that are already adapted to a food waste stream could be repurposed for useful chemical production. Industrial fermentations utilizing these microbes maintain organisms in isolation; however, microbial consortia offer an attractive alternative to monocultures in that metabolic interactions may result in more efficient processes with higher yields. Here we computationally assess the ability of co-cultures vs. monocultures to anaerobically convert a Western diet to commodity chemicals. The combination of genome-scale metabolic models with flux-balance analysis predicts that every organism analyzed can benefit from interactions with another microbe, as evidenced by increased biomass fluxes in co-culture vs. monoculture. Furthermore, microbe combinations result in emergent or increased commodity chemical production including butanol, methane, formaldehyde, propionate, hydrogen gas, and urea. These overproducing co-cultures are enriched for mutualistic and commensal interactions. Using Clostridium beijerinckii co-cultures as representative examples, models predict cross-fed metabolites will simultaneously modify multiple internal pathways, evident by different internal metabolic network structures. Differences in degree and betweenness centrality of hub precursor metabolites were correlated to C. beijerinckii metabolic outputs, and thus demonstrate the potential of co-cultures to differentially direct metabolisms to useful products.

Inferring Metabolic Mechanisms of Interaction within a Defined Gut Microbiota

The diversity and number of species present within microbial communities create the potential for a multitude of interspecies metabolic interactions. Here, we develop, apply, and experimentally test a framework for inferring metabolic mechanisms associated with interspecies interactions. We perform pairwise growth and metabolome profiling of co-cultures of strains from a model mouse microbiota. We then apply our framework to dissect emergent metabolic behaviors that occur in co-culture. Based on one of the inferences from this framework, we identify and interrogate an amino acid cross-feeding interaction and validate that the proposed interaction leads to a growth benefit in vitro. Our results reveal the type and extent of emergent metabolic behavior in microbial communities composed of gut microbes. We focus on growth-modulating interactions, but the framework can be applied to interspecies interactions that modulate any phenotype of interest within microbial communities.

From metagenomic data to personalized in silico microbiotas: predicting dietary supplements for Crohn's disease

Crohn's disease (CD) is associated with an ecological imbalance of the intestinal microbiota, consisting of hundreds of species. The underlying complexity as well as individual differences between patients contributes to the difficulty to define a standardized treatment. Computational modeling can systematically investigate metabolic interactions between gut microbes to unravel mechanistic insights. In this study, we integrated metagenomic data of CD patients and healthy controls with genome-scale metabolic models into personalized in silico microbiotas. We predicted short chain fatty acid (SFCA) levels for patients and controls, which were overall congruent with experimental findings. As an emergent property, low concentrations of SCFA were predicted for CD patients and the SCFA signatures were unique to each patient. Consequently, we suggest personalized dietary treatments that could improve each patient's SCFA levels. The underlying modeling approach could aid clinical practice to find dietary treatment and guide recovery by rationally proposing food aliments.