Systems modeling approaches for microbial community studies: from metagenomics to inference of the community structure

A framework for integrating genomics, microbial traits, and ecosystem biogeochemistry

Microbes drive the biogeochemical cycles of earth systems, yet the long-standing goal of linking emerging genomic information, microbial traits, mechanistic ecosystem models, and projections under climate change has remained elusive despite a wealth of emerging genomic information. Here we developed a general genome-to-ecosystem (G2E) framework for integrating genome-inferred microbial kinetic traits into mechanistic models of terrestrial ecosystems and applied it at a well-studied Arctic wetland by benchmarking predictions against observed greenhouse gas emissions. We found variation in genome-inferred microbial kinetic traits resulted in large differences in simulated annual methane emissions, quantitatively demonstrating that the genomically observable variations in microbial capacity are consequential for ecosystem functioning. Applying microbial community-aggregated traits via genome relative-abundance-weighting gave better methane emissions predictions (i.e., up to 54% decrease in bias) compared to ignoring the observed abundances, highlighting the value of combined trait inferences and abundances. This work provides an example of integrating microbial functional trait-based genomics, mechanistic and pragmatic trait parameterizations of diverse microbial metabolisms, and mechanistic ecosystem modeling. The generalizable G2E framework will enable the use of abundant microbial metagenomics data to improve predictions of microbial interactions in many complex systems, including oceanic microbiomes.

Modeling Microbial Communities: Perspective and Challenges

Microbial communities are immensely important due to their widespread presence and profound impact on various facets of life. Understanding these complex systems necessitates mathematical modeling, a powerful tool for simulating and predicting microbial community behavior. This review offers a critical analysis of metabolic modeling and highlights key areas that would greatly benefit from broader discussion and collaboration. Moreover, we explore the challenges and opportunities linked to the intricate nature of these communities, spanning data generation, modeling, and validation. We are confident that ongoing advancements in modeling techniques, such as machine learning, coupled with interdisciplinary collaborations, will unlock the full potential of microbial communities across diverse applications.

Effects of thermophilic and acidophilic microbial consortia on maize wet-milling steeping

To understand the ecology of species and promote biotechnology through beneficial strain selection for improving starch yield in maize wet-milling steeping, bacterial diversity and community structure during the counter-current steeping process in a commercial steeping system were characterized and investigated. The microbial diversity in the steeping liquor, which consisted of 16 phyla, 131 families, and 290 genera, was more abundant compared to those present on the surface of unsteeped maize. As the counter-current steeping progressed, exposing newer maize to the older steepwater, Lactobacillus dominated, replacing Rahnella, Pseudomonas, Pantoea, and Serratia. The thermophilic and acidophilic microbial consortia were enriched through adaptive evolution engineering and employed to improve starch yield. Several steeping strategies were evaluated, including water alone, SO2 alone, mono-culture of B. coagulans, microbial consortia, and a combination of consortium and SO2. Combining the microbial consortium with SO2 significantly increased the starch yield to, about 66.4 ± 0.5%, a 22% and 46% increase over SO2 alone and the consortium alone, respectively. Scanning electron microscope (SEM) of steeped maize structure indicated that the combination of consortium and SO2 disrupted the protein matrix and widened gaps between starch granules in maize endosperm. This released proteins into the steepwater and left starch granules in the aleurone layer. The steeping strategy of using thermophilic and acidophilic microbial consortium as additives shows potential application as an environmentally friendly alternative to conventional maize steeping procedures.

Microbially mediated climate feedbacks from wetland ecosystems

Wetlands are crucial nodes in the carbon cycle, emitting approximately 20% of global CH4 while also sequestering 20%-30% of all soil carbon. Both greenhouse gas fluxes and carbon storage are driven by microbial communities in wetland soils. However, these key players are often overlooked or overly simplified in current global climate models. Here, we first integrate microbial metabolisms with biological, chemical, and physical processes occurring at scales from individual microbial cells to ecosystems. This conceptual scale-bridging framework guides the development of feedback loops describing how wetland-specific climate impacts (i.e., sea level rise in estuarine wetlands, droughts and floods in inland wetlands) will affect future climate trajectories. These feedback loops highlight knowledge gaps that need to be addressed to develop predictive models of future climates capturing microbial contributions. We propose a roadmap connecting environmental scientific disciplines to address these knowledge gaps and improve the representation of microbial processes in climate models. Together, this paves the way to understand how microbially mediated climate feedbacks from wetlands will impact future climate change.

Upgrading dilute ethanol to odd-chain carboxylic acids by a synthetic co-culture of Anaerotignum neopropionicum and Clostridium kluyveri

BACKGROUND

Dilute ethanol streams generated during fermentation of biomass or syngas can be used as feedstocks for the production of higher value products. In this study, we describe a novel synthetic microbial co-culture that can effectively upgrade dilute ethanol streams to odd-chain carboxylic acids (OCCAs), specifically valerate and heptanoate. The co-culture consists of two strict anaerobic microorganisms: Anaerotignum neopropionicum, a propionigenic bacterium that ferments ethanol, and Clostridium kluyveri, well-known for its chain-elongating metabolism. In this co-culture, A. neopropionicum grows on ethanol and CO₂ producing propionate and acetate, which are then utilised by C. kluyveri for chain elongation with ethanol as the electron donor.

RESULTS

A co-culture of A. neopropionicum and C. kluyveri was established in serum bottles with 50 mM ethanol, leading to the production of valerate (5.4 ± 0.1 mM) as main product of ethanol-driven chain elongation. In a continuous bioreactor supplied with 3.1 g ethanol L^-1 d^-1, the co-culture exhibited high ethanol conversion (96.6%) and produced 25% (mol/mol) valerate, with a steady-state concentration of 8.5 mM and a rate of 5.7 mmol L^-1 d^-1. In addition, up to 6.5 mM heptanoate was produced at a rate of 2.9 mmol L^-1 d^-1. Batch experiments were also conducted to study the individual growth of the two strains on ethanol. A. neopropionicum showed the highest growth rate when cultured with 50 mM ethanol (μ_max = 0.103 ± 0.003 h^-1) and tolerated ethanol concentrations of up to 300 mM. Cultivation experiments with C. kluyveri showed that propionate and acetate were used simultaneously for chain elongation. However, growth on propionate alone (50 mM and 100 mM) led to a 1.8-fold reduction in growth rate compared to growth on acetate. Our results also revealed sub-optimal substrate use by C. kluyveri during odd-chain elongation, where excessive ethanol was oxidised to acetate.

CONCLUSIONS

This study highlights the potential of synthetic co-cultivation in chain elongation processes to target the production of OCCAs. Furthermore, our findings shed light on to the metabolism of odd-chain elongation by C. kluyveri.

Systematizing Microbial Bioplastic Production for Developing Sustainable Bioeconomy: Metabolic Nexus Modeling, Economic and Environmental Technologies Assessment

The excessive usage of non-renewable resources to produce plastic commodities has incongruously influenced the environment's health. Especially in the times of COVID-19, the need for plastic-based health products has increased predominantly. Given the rise in global warming and greenhouse gas emissions, the lifecycle of plastic has been established to contribute to it significantly. Bioplastics such as polyhydroxy alkanoates, polylactic acid, etc. derived from renewable energy origin have been a magnificent alternative to conventional plastics and reconnoitered exclusively for combating the environmental footprint of petrochemical plastic. However, the economically reasonable and environmentally friendly procedure of microbial bioplastic production has been a hard nut to crack due to less scouted and inefficient process optimization and downstream processing methodologies. Thereby, meticulous employment of computational tools such as genome-scale metabolic modeling and flux balance analysis has been practiced in recent times to understand the effect of genomic and environmental perturbations on the phenotype of the microorganism. In-silico results not only aid us in determining the biorefinery abilities of the model microorganism but also curb our reliance on equipment, raw materials, and capital investment for optimizing the best conditions. Additionally, to accomplish sustainable large-scale production of microbial bioplastic in a circular bioeconomy, extraction, and refinement of bioplastic needs to be investigated extensively by practicing techno-economic analysis and life cycle assessment. This review put forth state-of-the-art know-how on the proficiency of these computational techniques in laying the foundation of an efficient bioplastic manufacturing blueprint, chiefly focusing on microbial polyhydroxy alkanoates (PHA) production and its efficacy in outplacing fossil based plastic products.

A Review of Basic Bioinformatic Techniques for Microbial Community Analysis in an Anaerobic Digester

Biogas production involves various types of intricate microbial populations in an anaerobic digester (AD). To understand the anaerobic digestion system better, a broad-based study must be conducted on the microbial population. Deep understanding of the complete metagenomics including microbial structure, functional gene form, similarity/differences, and relationships between metabolic pathways and product formation, could aid in optimization and enhancement of AD processes. With advancements in technologies for metagenomic sequencing, for example, next generation sequencing and high-throughput sequencing, have revolutionized the study of microbial dynamics in anaerobic digestion. This review includes a brief introduction to the basic process of metagenomics research and includes a detailed summary of the various bioinformatics approaches, viz., total investigation of data obtained from microbial communities using bioinformatics methods to expose metagenomics characterization. This includes (1) methods of DNA isolation and sequencing, (2) investigation of anaerobic microbial communities using bioinformatics techniques, (3) application of the analysis of anaerobic microbial community and biogas production, and (4) restriction and prediction of bioinformatics analysis on microbial metagenomics. The review has been concluded, giving a summarized insight into bioinformatic tools and also promoting the future prospects of integrating humungous data with artificial intelligence and neural network software.

A Two-Compartment Fermentation System to Quantify Strain-Specific Interactions in Microbial Co-Cultures

To fulfil the growing interest in investigating microbial interactions in co-cultures, a novel two-compartment bioreactor system was developed, characterised, and implemented. The system allowed for the exchange of amino acids and peptides via a polyethersulfone membrane that retained biomass. Further system characterisation revealed a Bodenstein number of 18, which hints at backmixing. Together with other physical settings, the existence of unwanted inner-compartment substrate gradients could be ruled out. Furthermore, the study of Damkoehler numbers indicated that a proper metabolite supply between compartments was enabled. Implementing the two-compartment system (2cs) for growing Streptococcus thermophilus and Lactobacillus delbrueckii subs. bulgaricus, which are microorganisms commonly used in yogurt starter cultures, revealed only a small variance between the one-compartment and two-compartment approaches. The 2cs enabled the quantification of the strain-specific production and consumption rates of amino acids in an interacting S. thermophilus-L. bulgaricus co-culture. Therefore, comparisons between mono- and co-culture performance could be achieved. Both species produce and release amino acids. Only alanine was produced de novo from glucose through potential transaminase activity by L. bulgaricus and consumed by S. thermophilus. Arginine availability in peptides was limited to S. thermophilus' growth, indicating active biosynthesis and dependency on the proteolytic activity of L. bulgaricus. The application of the 2cs not only opens the door for the quantification of exchange fluxes between microbes but also enables continuous production modes, for example, for targeted evolution studies.

Biology and medicine in the landscape of quantum advantages

Quantum computing holds substantial potential for applications in biology and medicine, spanning from the simulation of biomolecules to machine learning methods for subtyping cancers on the basis of clinical features. This potential is encapsulated by the concept of a quantum advantage, which is contingent on a reduction in the consumption of a computational resource, such as time, space or data. Here, we distill the concept of a quantum advantage into a simple framework to aid researchers in biology and medicine pursuing the development of quantum applications. We then apply this framework to a wide variety of computational problems relevant to these domains in an effort to (i) assess the potential of practical advantages in specific application areas and (ii) identify gaps that may be addressed with novel quantum approaches. In doing so, we provide an extensive survey of the intersection of biology and medicine with the current landscape of quantum algorithms and their potential advantages. While we endeavour to identify specific computational problems that may admit practical advantages throughout this work, the rapid pace of change in the fields of quantum computing, classical algorithms and biological research implies that this intersection will remain highly dynamic for the foreseeable future.

Species abundance correlations carry limited information about microbial network interactions

Unraveling the network of interactions in ecological communities is a daunting task. Common methods to infer interspecific interactions from cross-sectional data are based on co-occurrence measures. For instance, interactions in the human microbiome are often inferred from correlations between the abundances of bacterial phylogenetic groups across subjects. We tested whether such correlation-based methods are indeed reliable for inferring interaction networks. For this purpose, we simulated bacterial communities by means of the generalized Lotka-Volterra model, with variation in model parameters representing variability among hosts. Our results show that correlations can be indicative for presence of bacterial interactions, but only when measurement noise is low relative to the variation in interaction strengths between hosts. Indication of interaction was affected by type of interaction network, process noise and sampling under non-equilibrium conditions. The sign of a correlation mostly coincided with the nature of the strongest pairwise interaction, but this is not necessarily the case. For instance, under rare conditions of identical interaction strength, we found that competitive and exploitative interactions can result in positive as well as negative correlations. Thus, cross-sectional abundance data carry limited information on specific interaction types. Correlations in abundance may hint at interactions but require independent validation.

The bacteria in and on our body (the human microbiome) largely determine how our body functions, and whether we stay healthy or get sick. These bacteria do not live on their own, but interact among each other and with their human host. Finding out which bacteria interact with each other is cumbersome, but patterns of joint occurrence between species might provide a clue to their ecological dependencies. We investigated whether correlations in species abundance can be used for the purpose of ecological network reconstruction. We simulated different bacterial communities with known interactions according to a theoretical population model. After having collected virtual samples from our simulated data, we performed a correlation analysis and then compared the correlation network with our known interaction network. We found that correlations can be informative for underlying interactions, but ecological conclusions should be drawn carefully. An obvious limitation of correlation analysis is that direction of interaction cannot be recovered from co-occurrence data, making correlations insensitive for detection of asymmetric interactions. In addition, we found that competitive and exploitative interactions can induce positive as well as negative correlations. We recommend careful interpretation and validation when inferring networks from cross-sectional abundance data.

NetCom: A Network-Based Tool for Predicting Metabolic Activities of Microbial Communities Based on Interpretation of Metagenomics Data

The study of microbial activity can be viewed as a triangle with three sides: environment (dominant resources in a specific habitat), community (species dictating a repertoire of metabolic conversions) and function (production and/or utilization of resources and compounds). Advances in metagenomics enable a high-resolution description of complex microbial communities in their natural environments and support a systematic study of environment-community-function associations. NetCom is a web-tool for predicting metabolic activities of microbial communities based on network-based interpretation of assembled and annotated metagenomics data. The algorithm takes as an input, lists of differentially abundant enzymatic reactions and generates the following outputs: (i) pathway associations of differently abundant enzymes; (ii) prediction of environmental resources that are unique to each treatment, and their pathway associations; (iii) prediction of compounds that are produced by the microbial community, and pathway association of compounds that are treatment-specific; (iv) network visualization of enzymes, environmental resources and produced compounds, that are treatment specific (2 and 3D). The tool is demonstrated on metagenomic data from rhizosphere and bulk soil samples. By predicting root-specific activities, we illustrate the relevance of our framework for forecasting the impact of soil amendments on the corresponding microbial communities. NetCom is available online.

Using Functional Annotations to Study Pairwise Interactions in Urinary Tract Infection Communities

The behaviour of microbial communities depends on environmental factors and on the interactions of the community members. This is also the case for urinary tract infection (UTI) microbial communities. Here, we devise a computational approach that uses indices of complementarity and competition based on metabolic gene annotation to rapidly predict putative interactions between pair of organisms with the aim to explain pairwise growth effects. We apply our method to 66 genomes selected from online databases, which belong to 6 genera representing members of UTI communities. This resulted in a selection of metabolic pathways with high correlation for each pairwise combination between a complementarity index and the experimentally derived growth data. Our results indicated that Enteroccus spp. were most complemented in its metabolism by the other members of the UTI community. This suggests that the growth of Enteroccus spp. can potentially be enhanced by complementary metabolites produced by other community members. We tested a few putative predicted interactions by experimental supplementation of the relevant predicted metabolites. As predicted by our method, folic acid supplementation led to the increase in the population density of UTI Enterococcus isolates. Overall, we believe our method is a rapid initial in silico screening for the prediction of metabolic interactions in microbial communities.

Mining Synergistic Microbial Interactions: A Roadmap on How to Integrate Multi-Omics Data

Mining interspecies interactions remain a challenge due to the complex nature of microbial communities and the need for computational power to handle big data. Our meta-analysis indicates that genetic potential alone does not resolve all issues involving mining of microbial interactions. Nevertheless, it can be used as the starting point to infer synergistic interspecies interactions and to limit the search space (i.e., number of species and metabolic reactions) to a manageable size. A reduced search space decreases the number of additional experiments necessary to validate the inferred putative interactions. As validation experiments, we examine how multi-omics and state of the art imaging techniques may further improve our understanding of species interactions’ role in ecosystem processes. Finally, we analyze pros and cons from the current methods to infer microbial interactions from genetic potential and propose a new theoretical framework based on: (i) genomic information of key members of a community; (ii) information of ecosystem processes involved with a specific hypothesis or research question; (iii) the ability to identify putative species’ contributions to ecosystem processes of interest; and, (iv) validation of putative microbial interactions through integration of other data sources.

Formation, characterization and modeling of emergent synthetic microbial communities

Microbial communities colonize plant tissues and contribute to host function. How these communities form and how individual members contribute to shaping the microbial community are not well understood. Synthetic microbial communities, where defined individual isolates are combined, can serve as valuable model systems for uncovering the organizational principles of communities. Using genome-defined organisms, systematic analysis by computationally-based network reconstruction can lead to mechanistic insights and the metabolic interactions between species. In this study, 10 bacterial strains isolated from the Populus deltoides rhizosphere were combined and passaged in two different media environments to form stable microbial communities. The membership and relative abundances of the strains stabilized after around 5 growth cycles and resulted in just a few dominant strains that depended on the medium. To unravel the underlying metabolic interactions, flux balance analysis was used to model microbial growth and identify potential metabolic exchanges involved in shaping the microbial communities. These analyses were complemented by growth curves of the individual isolates, pairwise interaction screens, and metaproteomics of the community. A fast growth rate is identified as one factor that can provide an advantage for maintaining presence in the community. Final community selection can also depend on selective antagonistic relationships and metabolic exchanges. Revealing the mechanisms of interaction among plant-associated microorganisms provides insights into strategies for engineering microbial communities that can potentially increase plant growth and disease resistance. Further, deciphering the membership and metabolic potentials of a bacterial community will enable the design of synthetic communities with desired biological functions.

One versus Many: Polymicrobial Communities and the Cystic Fibrosis Airway

Culture-independent studies have revealed that chronic lung infections in persons with cystic fibrosis (pwCF) are rarely limited to one microbial species. Interactions among bacterial members of these polymicrobial communities in the airways of pwCF have been reported to modulate clinically relevant phenotypes. Furthermore, it is clear that a single polymicrobial community in the context of CF airway infections cannot explain the diversity of clinical outcomes. While large 16S rRNA gene-based studies have allowed us to gain insight into the microbial composition and predicted functional capacities of communities found in the CF lung, here we argue that in silico approaches can help build clinically relevant in vitro models of polymicrobial communities that can in turn be used to experimentally test and validate computationally generated hypotheses. Furthermore, we posit that combining computational and experimental approaches will enhance our understanding of mechanisms that drive microbial community function and identify new therapeutics to target polymicrobial infections.

Modeling a co-culture of Clostridium autoethanogenum and Clostridium kluyveri to increase syngas conversion to medium-chain fatty-acids

•

We model a co-culture of C. autoethanogenum and C. kluyveri for syngas fermentation.

•

Biomass species ratio affects ethanol and acetate profiles in the co-culture.

•

The model predicts that addition of succinate increases caproate production.

•

Genetic interventions  in C. autoethanogenum could increase caproate production.

Microbial fermentation of synthesis gas (syngas) is becoming more attractive for sustainable production of commodity chemicals. To date, syngas fermentation focuses mainly on the use of Clostridium species for the production of small organic molecules such as ethanol and acetate. The co-cultivation of syngas-fermenting microorganisms with chain-elongating bacteria can expand the range of possible products, allowing, for instance, the production of medium-chain fatty acids (MCFA) and alcohols from syngas. To explore these possibilities, we report herein a genome-scale, constraint-based metabolic model to describe growth of a co-culture of Clostridium autoethanogenum and Clostridium kluyveri on syngas for the production of valuable compounds. Community flux balance analysis was used to gain insight into the metabolism of the two strains and their interactions, and to reveal potential strategies enabling production of butyrate and hexanoate.

The model suggests that one strategy to optimize the production of medium-chain fatty-acids from syngas would be the addition of succinate. According to the prediction, addition of succinate would increase the pool of crotonyl-CoA and the ethanol/acetate uptake ratio in C. kluyveri, resulting in a flux of up to 60% of electrons into hexanoate. Another potential way to further optimize butyrate and hexanoate production would be an increase of C. autoethanogenum ethanol production. Blocking either acetaldehyde dehydrogenase or formate dehydrogenase (ferredoxin) activity or formate transport, in the C. autoethanogenum metabolic model could potentially lead to an up to 150% increase in ethanol production.

Modeling microbial cross-feeding at intermediate scale portrays community dynamics and species coexistence

Social interaction between microbes can be described at many levels of details: from the biochemistry of cell-cell interactions to the ecological dynamics of populations. Choosing an appropriate level to model microbial communities without losing generality remains a challenge. Here we show that modeling cross-feeding interactions at an intermediate level between genome-scale metabolic models of individual species and consumer-resource models of ecosystems is suitable to experimental data. We applied our modeling framework to three published examples of multi-strain Escherichia coli communities with increasing complexity: uni-, bi-, and multi-directional cross-feeding of either substitutable metabolic byproducts or essential nutrients. The intermediate-scale model accurately fit empirical data and quantified metabolic exchange rates that are hard to measure experimentally, even for a complex community of 14 amino acid auxotrophies. By studying the conditions of species coexistence, the ecological outcomes of cross-feeding interactions, and each community’s robustness to perturbations, we extracted new quantitative insights from these three published experimental datasets. Our analysis provides a foundation to quantify cross-feeding interactions from experimental data, and highlights the importance of metabolic exchanges in the dynamics and stability of microbial communities.

The behavior of microbial communities such as the human microbiome is hard to predict by its species composition alone. Our efforts to engineer microbiomes—for example to improve human health—would benefit from mathematical models that accurately describe how microbes exchange metabolites with each other and how their environment shapes these exchanges. But what is an appropriate level of details for those models? We propose an intermediate level to model metabolic exchanges between microbes. We show that these models can accurately describe population dynamics in three laboratory communities and predicts their stability in response to perturbations such as changes in the nutrients available in the medium that they grow on. Our work suggests that a highly detailed metabolic network model is unnecessary for extracting ecological insights from experimental data and improves mathematical models so that one day we may be able to predict the behavior of real-world communities such as the human microbiome.

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

Resources and tools for the high-throughput, multi-omic study of intestinal microbiota

The human gut microbiome impacts several aspects of human health and disease, including digestion, drug metabolism and the propensity to develop various inflammatory, autoimmune and metabolic diseases. Many of the molecular processes that play a role in the activity and dynamics of the microbiota go beyond species and genic composition and thus, their understanding requires advanced bioinformatics support. This article aims to provide an up-to-date view of the resources and software tools that are being developed and used in human gut microbiome research, in particular data integration and systems-level analysis efforts. These efforts demonstrate the power of standardized and reproducible computational workflows for integrating and analysing varied omics data and gaining deeper insights into microbe community structure and function as well as host-microbe interactions.

Modelling microbial communities using biochemical resource allocation analysis

To understand the functioning and dynamics of microbial communities is a fundamental challenge in current biology. To tackle this challenge, the construction of computational models of interacting microbes is an indispensable tool. There is, however, a large chasm between ecologically motivated descriptions of microbial growth used in many current ecosystems simulations, and detailed metabolic pathway and genome-based descriptions developed in the context of systems and synthetic biology. Here, we seek to demonstrate how resource allocation models of microbial growth offer the potential to advance ecosystem simulations and their parametrization. In particular, recent work on quantitative resource allocation allow us to formulate mechanistic models of microbial growth that are physiologically meaningful while remaining computationally tractable. These models go beyond Michaelis-Menten and Monod-type growth models, and are capable of accounting for emergent properties that underlie the remarkable plasticity of microbial growth. We outline the utility and advantages of using biochemical resource allocation models by considering a coarse-grained model of cyanobacterial growth and demonstrate how the model allows us to address specific questions of relevance for the simulation of marine microbial ecosystems, including the physiological acclimation of protein expression to different environments, the description of co-limitation by several nutrients and the differential use of alternative nutrient sources, as well as the description of metabolic diversity based on our increasing knowledge about quantitative cell physiology.

Metabolic modelling of chronic wound microbiota predicts mutualistic interactions that drive community composition

AIMS

To identify putative mutualistic interactions driving community composition in polymicrobial chronic wound infections using metabolic modelling.

METHODS AND RESULTS

We developed a 12 species metabolic model that covered 74% of 16S rDNA pyrosequencing reads of dominant genera from 2963 chronic wound patients. The community model was used to predict species abundances averaged across this large patient population. We found that substantially improved predictions were obtained when the model was constrained with genera prevalence data and predicted abundances were averaged over 5000 ensemble simulations with community participants randomly determined according to the experimentally determined prevalences. Staphylococcus and Pseudomonas were predicted to exhibit a strong mutualistic relationship that resulted in community growth rate and diversity simultaneously increasing, suggesting that these two common chronic wound pathogens establish dominance by cooperating with less harmful commensal species. In communities lacking one or both dominant pathogens, other mutualistic relationship including Staphylococcus/Acinetobacter, Pseudomonas/Serratia and Streptococcus/Enterococcus were predicted consistent with published experimental data.

CONCLUSIONS

Mutualistic interactions were predicted to be driven by crossfeeding of organic acids, alcohols and amino acids that could potentially be disrupted to slow chronic wound disease progression.

SIGNIFICANCE AND IMPACT OF THE STUDY

Approximately 2% of the US population suffers from nonhealing chronic wounds infected by a combination of commensal and pathogenic bacteria. These polymicrobial infections are often resilient to antibiotic treatment due to the nutrient-rich wound environment and species interactions that promote community stability and robustness. The simulation results from this study were used to identify putative mutualistic interactions between bacteria that could be targeted to enhance treatment efficacy.

Synthetic microbial consortia for biosynthesis and biodegradation: promises and challenges

Functional differentiation and metabolite exchange enable microbial consortia to perform complex metabolic tasks and efficiently cycle the nutrients. Inspired by the cooperative relationships in environmental microbial consortia, synthetic microbial consortia have great promise for studying the microbial interactions in nature and more importantly for various engineering applications. However, challenges coexist with promises, and the potential of consortium-based technologies is far from being fully harnessed. Thorough understanding of the underlying molecular mechanisms of microbial interactions is greatly needed for the rational design and optimization of defined consortia. These knowledge gaps could be potentially filled with the assistance of the ongoing revolution in systems biology and synthetic biology tools. As current fundamental and technical obstacles down the road being removed, we would expect new avenues with synthetic microbial consortia playing important roles in biological and environmental engineering processes such as bioproduction of desired chemicals and fuels, as well as biodegradation of persistent contaminants.

Augmenting Biogas Process Modeling by Resolving Intracellular Metabolic Activity

The process of anaerobic digestion in which waste biomass is transformed to methane by complex microbial communities has been modeled for more than 16 years by parametric gray box approaches that simplify process biology and do not resolve intracellular microbial activity. Information on such activity, however, has become available in unprecedented detail by recent experimental advances in metatranscriptomics and metaproteomics. The inclusion of such data could lead to more powerful process models of anaerobic digestion that more faithfully represent the activity of microbial communities. We augmented the Anaerobic Digestion Model No. 1 (ADM1) as the standard kinetic model of anaerobic digestion by coupling it to Flux-Balance-Analysis (FBA) models of methanogenic species. Steady-state results of coupled models are comparable to standard ADM1 simulations if the energy demand for non-growth associated maintenance (NGAM) is chosen adequately. When changing a constant feed of maize silage from continuous to pulsed feeding, the final average methane production remains very similar for both standard and coupled models, while both the initial response of the methanogenic population at the onset of pulsed feeding as well as its dynamics between pulses deviates considerably. In contrast to ADM1, the coupled models deliver predictions of up to 1,000s of intracellular metabolic fluxes per species, describing intracellular metabolic pathway activity in much higher detail. Furthermore, yield coefficients which need to be specified in ADM1 are no longer required as they are implicitly encoded in the topology of the species’ metabolic network. We show the feasibility of augmenting ADM1, an ordinary differential equation-based model for simulating biogas production, by FBA models implementing individual steps of anaerobic digestion. While cellular maintenance is introduced as a new parameter, the total number of parameters is reduced as yield coefficients no longer need to be specified. The coupled models provide detailed predictions on intracellular activity of microbial species which are compatible with experimental data on enzyme synthesis activity or abundance as obtained by metatranscriptomics or metaproteomics. By providing predictions of intracellular fluxes of individual community members, the presented approach advances the simulation of microbial community driven processes and provides a direct link to validation by state-of-the-art experimental techniques.

Systems Biology - A Guide for Understanding and Developing Improved Strains of Lactic Acid Bacteria

Lactic Acid Bacteria (LAB) are extensively employed in the production of various fermented foods, due to their safe status, ability to affect texture and flavor and finally due to the beneficial effect they have on shelf-life. More recently, LAB have also gained interest as production hosts for various useful compounds, particularly compounds with sensitive applications, such as food ingredients and therapeutics. As for all industrial microorganisms, it is important to have a good understanding of the physiology and metabolism of LAB in order to fully exploit their potential, and for this purpose, many systems biology approaches are available. Systems metabolic engineering, an approach that combines optimization of metabolic enzymes/pathways at the systems level, synthetic biology as well as in silico model simulation, has been used to build microbial cell factories for production of biofuels, food ingredients and biochemicals. When developing LAB for use in foods, genetic engineering is in general not an accepted approach. An alternative is to screen mutant libraries for candidates with desirable traits using high-throughput screening technologies or to use adaptive laboratory evolution to select for mutants with special properties. In both cases, by using omics data and data-driven technologies to scrutinize these, it is possible to find the underlying cause for the desired attributes of such mutants. This review aims to describe how systems biology tools can be used for obtaining both engineered as well as non-engineered LAB with novel and desired properties.

Metabolic Modeling of Cystic Fibrosis Airway Communities Predicts Mechanisms of Pathogen Dominance

Cystic fibrosis (CF) is a genetic disease in which chronic airway infections and lung inflammation result in respiratory failure. CF airway infections are usually caused by bacterial communities that are difficult to eradicate with available antibiotics. Using species abundance data for clinically stable adult CF patients assimilated from three published studies, we developed a metabolic model of CF airway communities to better understand the interactions between bacterial species and between the bacterial community and the lung environment. Our model predicted that clinically observed CF pathogens could establish dominance over other community members across a range of lung nutrient conditions. Heterogeneity of species abundances across 75 patient samples could be predicted by assuming that sample-to-sample heterogeneity was attributable to random variations in the CF nutrient environment. Our model predictions provide new insights into the metabolic determinants of pathogen dominance in the CF lung and could facilitate the development of improved treatment strategies.

Cystic fibrosis (CF) is a fatal genetic disease characterized by chronic lung infections due to aberrant mucus production and the inability to clear invading pathogens. The traditional view that CF infections are caused by a single pathogen has been replaced by the realization that the CF lung usually is colonized by a complex community of bacteria, fungi, and viruses. To help unravel the complex interplay between the CF lung environment and the infecting microbial community, we developed a community metabolic model comprised of the 17 most abundant bacterial taxa, which account for >95% of reads across samples, from three published studies in which 75 sputum samples from 46 adult CF patients were analyzed by 16S rRNA gene sequencing. The community model was able to correctly predict high abundances of the “rare” pathogens Enterobacteriaceae, Burkholderia, and Achromobacter in three patients whose polymicrobial infections were dominated by these pathogens. With these three pathogens removed, the model correctly predicted that the remaining 43 patients would be dominated by Pseudomonas and/or Streptococcus. This dominance was predicted to be driven by relatively high monoculture growth rates of Pseudomonas and Streptococcus as well as their ability to efficiently consume amino acids, organic acids, and alcohols secreted by other community members. Sample-by-sample heterogeneity of community composition could be qualitatively captured through random variation of the simulated metabolic environment, suggesting that experimental studies directly linking CF lung metabolomics and 16S sequencing could provide important insights into disease progression and treatment efficacy.

IMPORTANCE Cystic fibrosis (CF) is a genetic disease in which chronic airway infections and lung inflammation result in respiratory failure. CF airway infections are usually caused by bacterial communities that are difficult to eradicate with available antibiotics. Using species abundance data for clinically stable adult CF patients assimilated from three published studies, we developed a metabolic model of CF airway communities to better understand the interactions between bacterial species and between the bacterial community and the lung environment. Our model predicted that clinically observed CF pathogens could establish dominance over other community members across a range of lung nutrient conditions. Heterogeneity of species abundances across 75 patient samples could be predicted by assuming that sample-to-sample heterogeneity was attributable to random variations in the CF nutrient environment. Our model predictions provide new insights into the metabolic determinants of pathogen dominance in the CF lung and could facilitate the development of improved treatment strategies.

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.

Constraint-based modeling in microbial food biotechnology

Genome-scale metabolic network reconstruction offers a means to leverage the value of the exponentially growing genomics data and integrate it with other biological knowledge in a structured format. Constraint-based modeling (CBM) enables both the qualitative and quantitative analyses of the reconstructed networks. The rapid advancements in these areas can benefit both the industrial production of microbial food cultures and their application in food processing. CBM provides several avenues for improving our mechanistic understanding of physiology and genotype–phenotype relationships. This is essential for the rational improvement of industrial strains, which can further be facilitated through various model-guided strain design approaches. CBM of microbial communities offers a valuable tool for the rational design of defined food cultures, where it can catalyze hypothesis generation and provide unintuitive rationales for the development of enhanced community phenotypes and, consequently, novel or improved food products. In the industrial-scale production of microorganisms for food cultures, CBM may enable a knowledge-driven bioprocess optimization by rationally identifying strategies for growth and stability improvement. Through these applications, we believe that CBM can become a powerful tool for guiding the areas of strain development, culture development and process optimization in the production of food cultures. Nevertheless, in order to make the correct choice of the modeling framework for a particular application and to interpret model predictions in a biologically meaningful manner, one should be aware of the current limitations of CBM.

Review and perspective on mathematical modeling of microbial ecosystems

Understanding microbial ecosystems means unlocking the path toward a deeper knowledge of the fundamental mechanisms of life. Engineered microbial communities are also extremely relevant to tackling some of today's grand societal challenges. Advanced meta-omics experimental techniques provide crucial insights into microbial communities, but have been so far mostly used for descriptive, exploratory approaches to answer the initial 'who is there?'

QUESTION

An ecosystem is a complex network of dynamic spatio-temporal interactions among organisms as well as between organisms and the environment. Mathematical models with their abstraction capability are essential to capture the underlying phenomena and connect the different scales at which these systems act. Differential equation models and constraint-based stoichiometric models are deterministic approaches that can successfully provide a macroscopic description of the outcome from microscopic behaviors. In this mini-review, we present classical and recent applications of these modeling methods and illustrate the potential of their integration. Indeed, approaches that can capture multiple scales are needed in order to understand emergent patterns in ecosystems and their dynamics regulated by different spatio-temporal phenomena. We finally discuss promising examples of methods proposing the integration of differential equations with constraint-based stoichiometric models and argue that more work is needed in this direction.

Modeling Microbial Communities: A Call for Collaboration between Experimentalists and Theorists

With our growing understanding of the impact of microbial communities, understanding how such communities function has become a priority. The influence of microbial communities is widespread. Human-associated microbiota impacts health, environmental microbes determine ecosystem sustainability, and microbe-driven industrial processes are expanding. This broad range of applications has led to a wide range of approaches to analyze and describe microbial communities. In particular, theoretical work based on mathematical modeling has been a steady source of inspiration for explaining and predicting microbial community processes. Here, we survey some of the modeling approaches used in different contexts. We promote classifying different approaches using a unified platform, and encourage cataloging the findings in a database. We believe that the synergy emerging from a coherent collection facilitates a better understanding of important processes that determine microbial community functions. We emphasize the importance of close collaboration between theoreticians and experimentalists in formulating, classifying, and improving models of microbial communities.

Analysis of Microbial Functions in the Rhizosphere Using a Metabolic-Network Based Framework for Metagenomics Interpretation

Advances in metagenomics enable high resolution description of complex bacterial communities in their natural environments. Consequently, conceptual approaches for community level functional analysis are in high need. Here, we introduce a framework for a metagenomics-based analysis of community functions. Environment-specific gene catalogs, derived from metagenomes, are processed into metabolic-network representation. By applying established ecological conventions, network-edges (metabolic functions) are assigned with taxonomic annotations according to the dominance level of specific groups. Once a function-taxonomy link is established, prediction of the impact of dominant taxa on the overall community performances is assessed by simulating removal or addition of edges (taxa associated functions). This approach is demonstrated on metagenomic data describing the microbial communities from the root environment of two crop plants – wheat and cucumber. Predictions for environment-dependent effects revealed differences between treatments (root vs. soil), corresponding to documented observations. Metabolism of specific plant exudates (e.g., organic acids, flavonoids) was linked with distinct taxonomic groups in simulated root, but not soil, environments. These dependencies point to the impact of these metabolite families as determinants of community structure. Simulations of the activity of pairwise combinations of taxonomic groups (order level) predicted the possible production of complementary metabolites. Complementation profiles allow formulating a possible metabolic role for observed co-occurrence patterns. For example, production of tryptophan-associated metabolites through complementary interactions is unique to the tryptophan-deficient cucumber root environment. Our approach enables formulation of testable predictions for species contribution to community activity and exploration of the functional outcome of structural shifts in complex bacterial communities. Understanding community-level metabolism is an essential step toward the manipulation and optimization of microbial function. Here, we introduce an analysis framework addressing three key challenges of such data: producing quantified links between taxonomy and function; contextualizing discrete functions into communal networks; and simulating environmental impact on community performances. New technologies will soon provide a high-coverage description of biotic and a-biotic aspects of complex microbial communities such as these found in gut and soil. This framework was designed to allow the integration of high-throughput metabolomic and metagenomic data toward tackling the intricate associations between community structure, community function, and metabolic inputs.

A Diverse Community To Study Communities: Integration of Experiments and Mathematical Models To Study Microbial Consortia

The last few years have seen the advancement of high-throughput experimental techniques that have produced an extraordinary amount of data. Bioinformatics and statistical analyses have become instrumental to interpreting the information coming from, e.g., sequencing data and often motivate further targeted experiments. The broad discipline of “computational biology” extends far beyond the well-established field of bioinformatics, but it is our impression that more theoretical methods such as the use of mathematical models are not yet as well integrated into the research studying microbial interactions. The empirical complexity of microbial communities presents challenges that are difficult to address with in vivo/in vitro approaches alone, and with microbiology developing from a qualitative to a quantitative science, we see stronger opportunities arising for interdisciplinary projects integrating theoretical approaches with experiments. Indeed, the addition of in silico experiments, i.e., computational simulations, has a discovery potential that is, unfortunately, still largely underutilized and unrecognized by the scientific community. This minireview provides an overview of mathematical models of natural ecosystems and emphasizes that one critical point in the development of a theoretical description of a microbial community is the choice of problem scale. Since this choice is mostly dictated by the biological question to be addressed, in order to employ theoretical models fully and successfully it is vital to implement an interdisciplinary view at the conceptual stages of the experimental design.

Perspectives and Challenges in Microbial Communities Metabolic Modeling

Bacteria have evolved to efficiently interact each other, forming complex entities known as microbial communities. These "super-organisms" play a central role in maintaining the health of their eukaryotic hosts and in the cycling of elements like carbon and nitrogen. However, despite their crucial importance, the mechanisms that influence the functioning of microbial communities and their relationship with environmental perturbations are obscure. The study of microbial communities was boosted by tremendous advances in sequencing technologies, and in particular by the possibility to determine genomic sequences of bacteria directly from environmental samples. Indeed, with the advent of metagenomics, it has become possible to investigate, on a previously unparalleled scale, the taxonomical composition and the functional genetic elements present in a specific community. Notwithstanding, the metagenomic approach per se suffers some limitations, among which the impossibility of modeling molecular-level (e.g., metabolic) interactions occurring between community members, as well as their effects on the overall stability of the entire system. The family of constraint-based methods, such as flux balance analysis, has been fruitfully used to translate genome sequences in predictive, genome-scale modeling platforms. Although these techniques have been initially developed for analyzing single, well-known model organisms, their recent improvements allowed engaging in multi-organism in silico analyses characterized by a considerable predictive capability. In the face of these advances, here we focus on providing an overview of the possibilities and challenges related to the modeling of metabolic interactions within a bacterial community, discussing the feasibility and the perspectives of this kind of analysis in the (near) future.

Model-based quantification of metabolic interactions from dynamic microbial-community data

An important challenge in microbial ecology is to infer metabolic-exchange fluxes between growing microbial species from community-level data, concerning species abundances and metabolite concentrations. Here we apply a model-based approach to integrate such experimental data and thereby infer metabolic-exchange fluxes. We designed a synthetic anaerobic co-culture of Clostridium acetobutylicum and Wolinella succinogenes that interact via interspecies hydrogen transfer and applied different environmental conditions for which we expected the metabolic-exchange rates to change. We used stoichiometric models of the metabolism of the two microorganisms that represents our current physiological understanding and found that this understanding - the model - is sufficient to infer the identity and magnitude of the metabolic-exchange fluxes and it suggested unexpected interactions. Where the model could not fit all experimental data, it indicates specific requirement for further physiological studies. We show that the nitrogen source influences the rate of interspecies hydrogen transfer in the co-culture. Additionally, the model can predict the intracellular fluxes and optimal metabolic exchange rates, which can point to engineering strategies. This study therefore offers a realistic illustration of the strengths and weaknesses of model-based integration of heterogenous data that makes inference of metabolic-exchange fluxes possible from community-level experimental data.

A multi-objective constraint-based approach for modeling genome-scale microbial ecosystems

Interplay within microbial communities impacts ecosystems on several scales, and elucidation of the consequent effects is a difficult task in ecology. In particular, the integration of genome-scale data within quantitative models of microbial ecosystems remains elusive. This study advocates the use of constraint-based modeling to build predictive models from recent high-resolution -omics datasets. Following recent studies that have demonstrated the accuracy of constraint-based models (CBMs) for simulating single-strain metabolic networks, we sought to study microbial ecosystems as a combination of single-strain metabolic networks that exchange nutrients. This study presents two multi-objective extensions of CBMs for modeling communities: multi-objective flux balance analysis (MO-FBA) and multi-objective flux variability analysis (MO-FVA). Both methods were applied to a hot spring mat model ecosystem. As a result, multiple trade-offs between nutrients and growth rates, as well as thermodynamically favorable relative abundances at community level, were emphasized. We expect this approach to be used for integrating genomic information in microbial ecosystems. Following models will provide insights about behaviors (including diversity) that take place at the ecosystem scale.

High definition for systems biology of microbial communities: metagenomics gets genome-centric and strain-resolved

The systems biology of microbial communities, organismal communities inhabiting all ecological niches on earth, has in recent years been strongly facilitated by the rapid development of experimental, sequencing and data analysis methods. Novel experimental approaches and binning methods in metagenomics render the semi-automatic reconstructions of near-complete genomes of uncultivable bacteria possible, while advances in high-resolution amplicon analysis allow for efficient and less biased taxonomic community characterization. This will also facilitate predictive modeling approaches, hitherto limited by the low resolution of metagenomic data. In this review, we pinpoint the most promising current developments in metagenomics. They facilitate microbial systems biology towards a systemic understanding of mechanisms in microbial communities with scopes of application in many areas of our daily life.

Systems Biology of Microbial Exopolysaccharides Production

Exopolysaccharides (EPSs) produced by diverse group of microbial systems are rapidly emerging as new and industrially important biomaterials. Due to their unique and complex chemical structures and many interesting physicochemical and rheological properties with novel functionality, the microbial EPSs find wide range of commercial applications in various fields of the economy such as food, feed, packaging, chemical, textile, cosmetics and pharmaceutical industry, agriculture, and medicine. EPSs are mainly associated with high-value applications, and they have received considerable research attention over recent decades with their biocompatibility, biodegradability, and both environmental and human compatibility. However, only a few microbial EPSs have achieved to be used commercially due to their high production costs. The emerging need to overcome economic hurdles and the increasing significance of microbial EPSs in industrial and medical biotechnology call for the elucidation of the interrelations between metabolic pathways and EPS biosynthesis mechanism in order to control and hence enhance its microbial productivity. Moreover, a better understanding of biosynthesis mechanism is a significant issue for improvement of product quality and properties and also for the design of novel strains. Therefore, a systems-based approach constitutes an important step toward understanding the interplay between metabolism and EPS biosynthesis and further enhances its metabolic performance for industrial application. In this review, primarily the microbial EPSs, their biosynthesis mechanism, and important factors for their production will be discussed. After this brief introduction, recent literature on the application of omics technologies and systems biology tools for the improvement of production yields will be critically evaluated. Special focus will be given to EPSs with high market value such as xanthan, levan, pullulan, and dextran.