Constraint-based stoichiometric modelling from single organisms to microbial communities

Deciphering and designing microbial communities by genome-scale metabolic modelling

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

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Genome-scale model of Rothia mucilaginosa predicts gene essentialities and reveals metabolic capabilities

Cystic fibrosis (CF), an inherited genetic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator gene, results in sticky and thick mucosal fluids. This environment facilitates the colonization of various microorganisms, some of which can cause acute and chronic lung infections, while others may positively impact the disease. Rothia mucilaginosa, an oral commensal, is relatively abundant in the lungs of CF patients. Recent studies have unveiled its anti-inflammatory properties using in vitro three-dimensional lung epithelial cell cultures and in vivo mouse models relevant to chronic lung diseases. Apart from this, R. mucilaginosa has been associated with severe infections. However, its metabolic capabilities and genotype-phenotype relationships remain largely unknown. To gain insights into its cellular metabolism and genetic content, we developed the first manually curated genome-scale metabolic model, iRM23NL. Through growth kinetics and high-throughput phenotypic microarray testings, we defined its complete catabolic phenome. Subsequently, we assessed the model's effectiveness in accurately predicting growth behaviors and utilizing multiple substrates. We used constraint-based modeling techniques to formulate novel hypotheses that could expedite the development of antimicrobial strategies. More specifically, we detected putative essential genes and assessed their effect on metabolism under varying nutritional conditions. These predictions could offer novel potential antimicrobial targets without laborious large-scale screening of knockouts and mutant transposon libraries. Overall, iRM23NL demonstrates a solid capability to predict cellular phenotypes and holds immense potential as a valuable resource for accurate predictions in advancing antimicrobial therapies. Moreover, it can guide metabolic engineering to tailor R. mucilaginosa's metabolism for desired performance.IMPORTANCECystic fibrosis (CF) is a genetic disorder characterized by thick mucosal secretions, leading to chronic lung infections. Rothia mucilaginosa is a common bacterium found in various parts of the human body, acting as a normal part of the flora. In people with weakened immune systems, it can become an opportunistic pathogen, while it is prevalent and active in CF airways. Recent studies have highlighted its anti-inflammatory properties in the lower pulmonary system, indicating the intricate relationship between microbes and human health. Herein, we have developed the first manually curated metabolic model of R. mucilaginosa. Our study examined the previously unknown relationships between the bacterium's genotype and phenotype and identified essential genes that impact the metabolism under various conditions. With this, we opt for paving the way for developing new strategies in antimicrobial therapy and metabolic engineering, leading to enhanced therapeutic outcomes in cystic fibrosis and related conditions.

Formation of a constructed microbial community in a nutrient-rich environment indicates bacterial interspecific competition

Understanding the organizational principles of microbial communities is essential for interpreting ecosystem stability. Previous studies have investigated the formation of bacterial communities under nutrient-poor conditions or obligate relationships to observe cooperative interactions among different species. How microorganisms form stabilized communities in nutrient-rich environments, without obligate metabolic interdependency for growth, is still not fully disclosed. In this study, three bacterial strains isolated from the Populus deltoides rhizosphere were co-cultured in complex medium, and their growth behavior was tracked. These strains co-exist in mixed culture over serial transfer for multiple growth-dilution cycles. Competition is proposed as an emergent interaction relationship among the three bacteria based on their significantly decreased growth levels. The effects of different initial inoculum ratios, up to three orders of magnitude, on community structure were investigated, and the final compositions of the mixed communities with various starting composition indicate that community structure is not dependent on the initial inoculum ratio. Furthermore, the competitive relationships within the community were not altered by different initial inoculum ratios. The community structure was simulated by generalized Lotka-Volterra and dynamic flux balance analysis to provide mechanistic predictions into emergence of community structure under a nutrient-rich environment. Metaproteomic analyses provide support for the metabolite exchanges predicted by computational modeling and for highly altered physiologies when microbes are grown in co-culture. These findings broaden our understanding of bacterial community dynamics and metabolic diversity in higher-order interactions and could be significant in the management of rhizospheric bacterial communities.

IMPORTANCE

Bacteria naturally co-exist in multispecies consortia, and the ability to engineer such systems can be useful in biotechnology. Despite this, few studies have been performed to understand how bacteria form a stable community and interact with each other under nutrient-rich conditions. In this study, we investigated the effects of initial inoculum ratios on bacterial community structure using a complex medium and found that the initial inoculum ratio has no significant impact on resultant community structure or on interaction patterns between community members. The microbial population profiles were simulated using computational tools in order to understand intermicrobial relationships and to identify potential metabolic exchanges that occur during stabilization of the bacterial community. Studying microbial community assembly processes is essential for understanding fundamental ecological principles in microbial ecosystems and can be critical in predicting microbial community structure and function.

Metabolic model predictions enable targeted microbiome manipulation through precision prebiotics

While numerous health-beneficial interactions between host and microbiota have been identified, there is still a lack of targeted approaches for modulating these interactions. Thus, we here identify precision prebiotics that specifically modulate the abundance of a microbiome member species of interest. In the first step, we show that defining precision prebiotics by compounds that are only taken up by the target species but no other species in a community is usually not possible due to overlapping metabolic niches. Subsequently, we use metabolic modeling to identify precision prebiotics for a two-member Caenorhabditis elegans microbiome community comprising the immune-protective target species Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71. We experimentally confirm four of the predicted precision prebiotics, L-serine, L-threonine, D-mannitol, and γ-aminobutyric acid, to specifically increase the abundance of MYb11. L-serine was further assessed in vivo, leading to an increase in MYb11 abundance also in the worm host. Overall, our findings demonstrate that metabolic modeling is an effective tool for the design of precision prebiotics as an important cornerstone for future microbiome-targeted therapies.IMPORTANCEWhile various mechanisms through which the microbiome influences disease processes in the host have been identified, there are still only few approaches that allow for targeted manipulation of microbiome composition as a first step toward microbiome-based therapies. Here, we propose the concept of precision prebiotics that allow to boost the abundance of already resident health-beneficial microbial species in a microbiome. We present a constraint-based modeling pipeline to predict precision prebiotics for a minimal microbial community in the worm Caenorhabditis elegans comprising the host-beneficial Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71 with the aim to boost the growth of MYb11. Experimentally testing four of the predicted precision prebiotics, we confirm that they are specifically able to increase the abundance of MYb11 in vitro and in vivo. These results demonstrate that constraint-based modeling could be an important tool for the development of targeted microbiome-based therapies against human diseases.

The proteome is a terminal electron acceptor

Microbial metabolism is impressively flexible, enabling growth even when available nutrients differ greatly from biomass in redox state. E. coli, for example, rearranges its physiology to grow on reduced and oxidized carbon sources through several forms of fermentation and respiration. To understand the limits on and evolutionary consequences of metabolic flexibility, we developed a mathematical model coupling redox chemistry with principles of cellular resource allocation. Our integrated model clarifies key phenomena, including demonstrating that autotrophs grow slower than heterotrophs because of constraints imposed by intracellular production of reduced carbon. Our model further indicates that growth is improved by adapting the redox state of biomass to nutrients, revealing an unexpected mode of evolution where proteins accumulate mutations benefiting organismal redox balance.

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

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

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Toward a modeling, optimization, and predictive control framework for fed-batch metabolic cybergenetics

Biotechnology offers many opportunities for the sustainable manufacturing of valuable products. The toolbox to optimize bioprocesses includes extracellular process elements such as the bioreactor design and mode of operation, medium formulation, culture conditions, feeding rates, and so on. However, these elements are frequently insufficient for achieving optimal process performance or precise product composition. One can use metabolic and genetic engineering methods for optimization at the intracellular level. Nevertheless, those are often of static nature, failing when applied to dynamic processes or if disturbances occur. Furthermore, many bioprocesses are optimized empirically and implemented with little-to-no feedback control to counteract disturbances. The concept of cybergenetics has opened new possibilities to optimize bioprocesses by enabling online modulation of the gene expression of metabolism-relevant proteins via external inputs (e.g., light intensity in optogenetics). Here, we fuse cybergenetics with model-based optimization and predictive control for optimizing dynamic bioprocesses. To do so, we propose to use dynamic constraint-based models that integrate the dynamics of metabolic reactions, resource allocation, and inducible gene expression. We formulate a model-based optimal control problem to find the optimal process inputs. Furthermore, we propose using model predictive control to address uncertainties via online feedback. We focus on fed-batch processes, where the substrate feeding rate is an additional optimization variable. As a simulation example, we show the optogenetic control of the ATPase enzyme complex for dynamic modulation of enforced ATP wasting to adjust product yield and productivity.

Searching for Principles of Microbial Ecology Across Levels of Biological Organization

Microbial communities play pivotal roles in ecosystems across different scales, from global elemental cycles to household food fermentations. These complex assemblies comprise hundreds or thousands of microbial species whose abundances vary over time and space. Unraveling the principles that guide their dynamics at different levels of biological organization, from individual species, their interactions, to complex microbial communities, is a major challenge. To what extent are these different levels of organization governed by separate principles, and how can we connect these levels to develop predictive models for the dynamics and function of microbial communities? Here, we will discuss recent advances that point towards principles of microbial communities, rooted in various disciplines from physics, biochemistry, and dynamical systems. By considering the marine carbon cycle as a concrete example, we demonstrate how the integration of levels of biological organization can offer deeper insights into the impact of increasing temperatures, such as those associated with climate change, on ecosystem-scale processes. We argue that by focusing on principles that transcend specific microbiomes, we can pave the way for a comprehensive understanding of microbial community dynamics and the development of predictive models for diverse ecosystems.

Inferring microbial interactions with their environment from genomic and metagenomic data

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

Unveiling abundance-dependent metabolic phenotypes of microbial communities

IMPORTANCE

In nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question. In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

Shifts in carbon partitioning by photosynthetic activity increase terpenoid synthesis in glandular trichomes

Several commercially important secondary metabolites are produced and accumulated in high amounts by glandular trichomes, giving the prospect of using them as metabolic cell factories. Due to extremely high metabolic fluxes through glandular trichomes, previous research focused on how such flows are achieved. The question regarding their bioenergetics became even more interesting with the discovery of photosynthetic activity in some glandular trichomes. Despite recent advances, how primary metabolism contributes to the high metabolic fluxes in glandular trichomes is still not fully elucidated. Using computational methods and available multi-omics data, we first developed a quantitative framework to investigate the possible role of photosynthetic energy supply in terpenoid production and next tested experimentally the simulation-driven hypothesis. With this work, we provide the first reconstruction of specialised metabolism in Type-VI photosynthetic glandular trichomes of Solanum lycopersicum. Our model predicted that increasing light intensities results in a shift of carbon partitioning from catabolic to anabolic reactions driven by the energy availability of the cell. Moreover, we show the benefit of shifting between isoprenoid pathways under different light regimes, leading to a production of different classes of terpenes. Our computational predictions were confirmed in vivo, demonstrating a significant increase in production of monoterpenoids while the sesquiterpenes remained unchanged under higher light intensities. The outcomes of this research provide quantitative measures to assess the beneficial role of chloroplast in glandular trichomes for enhanced production of secondary metabolites and can guide the design of new experiments that aim at modulating terpenoid production.

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

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

Metabolic model predictions enable targeted microbiome manipulation through precision prebiotics

The microbiome is increasingly receiving attention as an important modulator of host health and disease. However, while numerous mechanisms through which the microbiome influences its host have been identified, there is still a lack of approaches that allow to specifically modulate the abundance of individual microbes or microbial functions of interest. Moreover, current approaches for microbiome manipulation such as fecal transfers often entail a non-specific transfer of entire microbial communities with potentially unwanted side effects. To overcome this limitation, we here propose the concept of precision prebiotics that specifically modulate the abundance of a microbiome member species of interest. In a first step, we show that defining precision prebiotics by compounds that are only taken up by the target species but no other species in a community is usually not possible due to overlapping metabolic niches. Subsequently, we present a metabolic modeling network framework that allows us to define precision prebiotics for a two-member C. elegans microbiome model community comprising the immune-protective Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71. Thus, we predicted compounds that specifically boost the abundance of the host-beneficial MYb11, four of which were experimentally validated in vitro (L-serine, L-threonine, D-mannitol, and γ-aminobutyric acid). L-serine was further assessed in vivo, leading to an increase in MYb11 abundance also in the worm host. Overall, our findings demonstrate that constraint-based metabolic modeling is an effective tool for the design of precision prebiotics as an important cornerstone for future microbiome-targeted therapies.

Spatial self-organization of metabolism in microbial systems: A matter of enzymes and chemicals

Most bacteria live in dense, spatially structured communities such as biofilms. The high density allows cells to alter the local microenvironment, whereas the limited mobility can cause species to become spatially organized. Together, these factors can spatially organize metabolic processes within microbial communities so that cells in different locations perform different metabolic reactions. The overall metabolic activity of a community depends both on how metabolic reactions are arranged in space and on how they are coupled, i.e., how cells in different regions exchange metabolites. Here, we review mechanisms that lead to the spatial organization of metabolic processes in microbial systems. We discuss factors that determine the length scales over which metabolic activities are arranged in space and highlight how the spatial organization of metabolic processes affects the ecology and evolution of microbial communities. Finally, we define key open questions that we believe should be the main focus of future research.

Unlocking the magic in mycelium: Using synthetic biology to optimize filamentous fungi for biomanufacturing and sustainability

Filamentous fungi drive carbon and nutrient cycling across our global ecosystems, through its interactions with growing and decaying flora and their constituent microbiomes. The remarkable metabolic diversity, secretion ability, and fiber-like mycelial structure that have evolved in filamentous fungi have been increasingly exploited in commercial operations. The industrial potential of mycelial fermentation ranges from the discovery and bioproduction of enzymes and bioactive compounds, the decarbonization of food and material production, to environmental remediation and enhanced agricultural production. Despite its fundamental impact in ecology and biotechnology, molds and mushrooms have not, to-date, significantly intersected with synthetic biology in ways comparable to other industrial cell factories (e.g. Escherichia coli,Saccharomyces cerevisiae, and Komagataella phaffii). In this review, we summarize a suite of synthetic biology and computational tools for the mining, engineering and optimization of filamentous fungi as a bioproduction chassis. A combination of methods across genetic engineering, mutagenesis, experimental evolution, and computational modeling can be used to address strain development bottlenecks in established and emerging industries. These include slow mycelium growth rate, low production yields, non-optimal growth in alternative feedstocks, and difficulties in downstream purification. In the scope of biomanufacturing, we then detail previous efforts in improving key bottlenecks by targeting protein processing and secretion pathways, hyphae morphogenesis, and transcriptional control. Bringing synthetic biology practices into the hidden world of molds and mushrooms will serve to expand the limited panel of host organisms that allow for commercially-feasible and environmentally-sustainable bioproduction of enzymes, chemicals, therapeutics, foods, and materials of the future.

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.

Flux balance analysis of metabolic networks for efficient engineering of microbial cell factories

Metabolic engineering principles have long been applied to explore the metabolic networks of complex microbial cell factories under a variety of environmental constraints for effective deployment of the microorganisms in the optimal production of biochemicals like biofuels, polymers, amino acids, recombinant proteins. One of the methodologies used for analyzing microbial metabolic networks is the Flux Balance Analysis (FBA), which employs applications of optimization techniques for forecasting biomass growth and metabolic flux distribution of industrially important products under specified environmental conditions. The in silico flux simulations are instrumental for designing the production-specific microbial cell factories. However, FBA has some inherent limitations. The present review emphasizes how the incorporation of additional kinetic, thermodynamic, expression and regulatory constraints and integration of omics data into the classical FBA platform improve the prediction accuracy of FBA. A programmed comparison of the simulated data with the experimental observations is presented for supporting the claim. The review further accounts for the successful implementation of classical FBA in biotechnological applications and identifies areas in which classical FBA fails to make correct predictions. The analysis of the predictive capabilities of the different FBA strategies presented here is expected to help researchers in finding new avenues in engineering highly efficient microbial metabolic pathways and identify the key metabolic bottlenecks during the process. Based on the appropriate metabolic network design, fermentation engineers will be able to effectively design the bioreactors and optimize large-scale biochemical production through suitable pathway modifications.

Towards the human nasal microbiome: Simulating D. pigrum and S. aureus

The human nose harbors various microbes that decisively influence the wellbeing and health of their host. Among the most threatening pathogens in this habitat is Staphylococcus aureus. Multiple epidemiological studies identify Dolosigranulum pigrum as a likely beneficial bacterium based on its positive association with health, including negative associations with S. aureus. Carefully curated GEMs are available for both bacterial species that reliably simulate their growth behavior in isolation. To unravel the mutual effects among bacteria, building community models for simulating co-culture growth is necessary. However, modeling microbial communities remains challenging. This article illustrates how applying the NCMW fosters our understanding of two microbes' joint growth conditions in the nasal habitat and their intricate interplay from a metabolic modeling perspective. The resulting community model combines the latest available curated GEMs of D. pigrum and S. aureus. This uses case illustrates how to incorporate genuine GEM of participating microorganisms and creates a basic community model mimicking the human nasal environment. Our analysis supports the role of negative microbe-microbe interactions involving D. pigrum examined experimentally in the lab. By this, we identify and characterize metabolic exchange factors involved in a specific interaction between D. pigrum and S. aureus as an in silico candidate factor for a deep insight into the associated species. This method may serve as a blueprint for developing more complex microbial interaction models. Its direct application suggests new ways to prevent disease-causing infections by inhibiting the growth of pathogens such as S. aureus through microbe-microbe interactions.

Integrating metagenomic binning with flux balance analysis to unravel syntrophies in anaerobic CO2 methanation

BACKGROUND

Carbon fixation through biological methanation has emerged as a promising technology to produce renewable energy in the context of the circular economy. The anaerobic digestion microbiome is the fundamental biological system operating biogas upgrading and is paramount in power-to-gas conversion. Carbon dioxide (CO₂) methanation is frequently performed by microbiota attached to solid supports generating biofilms. Despite the apparent simplicity of the microbial community involved in biogas upgrading, the dynamics behind most of the interspecies interaction remain obscure. To understand the role of the microbial species in CO₂ fixation, the biofilm generated during the biogas upgrading process has been selected as a case study. The present work investigates via genome-centric metagenomics, based on a hybrid Nanopore-Illumina approach the biofilm developed on the diffusion devices of four ex situ biogas upgrading reactors. Moreover, genome-guided metabolic reconstruction and flux balance analysis were used to propose a biological role for the dominant microbes.

RESULTS

The combined microbiome was composed of 59 species, with five being dominant (> 70% of total abundance); the metagenome-assembled genomes representing these species were refined to reach a high level of completeness. Genome-guided metabolic analysis appointed Firmicutes sp. GSMM966 as the main responsible for biofilm formation. Additionally, species interactions were investigated considering their co-occurrence in 134 samples, and in terms of metabolic exchanges through flux balance simulation in a simplified medium. Some of the most abundant species (e.g., Limnochordia sp. GSMM975) were widespread (~ 67% of tested experiments), while others (e.g., Methanothermobacter wolfeii GSMM957) had a scattered distribution. Genome-scale metabolic models of the microbial community were built with boundary conditions taken from the biochemical data and showed the presence of a flexible interaction network mainly based on hydrogen and carbon dioxide uptake and formate exchange.

CONCLUSIONS

Our work investigated the interplay between five dominant species within the biofilm and showed their importance in a large spectrum of anaerobic biogas reactor samples. Flux balance analysis provided a deeper insight into the potential syntrophic interaction between species, especially Limnochordia sp. GSMM975 and Methanothermobacter wolfeii GSMM957. Finally, it suggested species interactions to be based on formate and amino acids exchanges. Video Abstract.

Modeling-Guided Amendments Lead to Enhanced Biodegradation in Soil

Extensive use of agrochemicals is emerging as a serious environmental issue coming at the cost of the pollution of soil and water resources. Bioremediation techniques such as biostimulation are promising strategies used to remove pollutants from agricultural soils by supporting the indigenous microbial degraders. Though considered cost-effective and eco-friendly, the success rate of these strategies typically varies, and consequently, they are rarely integrated into commercial agricultural practices. In the current study, we applied metabolic-based community-modeling approaches for promoting realistic in terra solutions by simulation-based prioritization of alternative supplements as potential biostimulants, considering a collection of indigenous bacteria. Efficacy of biostimulants as enhancers of the indigenous degrader Paenarthrobacter was ranked through simulation and validated in pot experiments. A two-dimensional simulation matrix predicting the effect of different biostimulants on additional potential indigenous degraders (Pseudomonas, Clostridium, and Geobacter) was crossed with experimental observations. The overall ability of the models to predict the compounds that act as taxa-selective stimulants indicates that computational algorithms can guide the manipulation of the soil microbiome in situ and provides an additional step toward the educated design of biostimulation strategies. IMPORTANCE Providing the food requirements of a growing population comes at the cost of intensive use of agrochemicals, including pesticides. Native microbial soil communities are considered key players in the degradation of such exogenous substances. Manipulating microbial activity toward an optimized outcome in efficient biodegradation processes conveys a promise of maintaining intensive yet sustainable agriculture. Efficient strategies for harnessing the native microbiome require the development of approaches for processing big genomic data. Here, we pursued metabolic modeling for promoting realistic in terra solutions by simulation-based prioritization of alternative supplements as potential biostimulants, considering a collection of indigenous bacteria. Our genomic-based predictions point at strategies for optimizing biodegradation by the native community. Developing a systematic, data-guided understanding of metabolite-driven targeted enhancement of selected microorganisms lays the foundation for the design of ecologically sound methods for optimizing microbiome functioning.

Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework

Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this framework's success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that when genetic perturbations were implemented in the systems, the optimal constraint approach typically remained the same as for the wild-type regardless of the model topology or parameterization, indicating that just a single wild-type dataset could allow identification of the ideal constraint to enable model predictivity for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems.

Metagenome-Scale Metabolic Network Suggests Folate Produced by Bifidobacterium longum Might Contribute to High-Fiber-Diet-Induced Weight Loss in a Prader-Willi Syndrome Child

Gut-microbiota-targeted nutrition intervention has achieved success in the management of obesity, but its underlying mechanism still needs extended exploration. An obese Prader-Willi syndrome boy lost 25.8 kg after receiving a high-fiber dietary intervention for 105 days. The fecal microbiome sequencing data taken from the boy on intervention days 0, 15, 30, 45, 60, 75, and 105, along with clinical indexes, were used to construct a metagenome-scale metabolic network. Firstly, the abundances of the microbial strains were obtained by mapping the sequencing reads onto the assembly of gut organisms through use of reconstruction and analysis (AGORA) genomes. The nutritional components of the diet were obtained through the Virtual Metabolic Human database. Then, a community model was simulated using the Microbiome Modeling Toolbox. Finally, the significant Spearman correlations among the metabolites and the clinical indexes were screened and the strains that were producing these metabolites were identified. The high-fiber diet reduced the overall amount of metabolite secretions, but the secretions of folic acid derivatives by Bifidobacterium longum strains were increased and were significantly relevant to the observed weight loss. Reduced metabolites might also have directly contributed to the weight loss or indirectly contribute by enhancing leptin and decreasing adiponectin. Metagenome-scale metabolic network technology provides a cost-efficient solution for screening the functional microbial strains and metabolic pathways that are responding to nutrition therapy.

High-Quality Genome-Scale Reconstruction of Corynebacterium glutamicum ATCC 13032

Corynebacterium glutamicum belongs to the microbes of enormous biotechnological relevance. In particular, its strain ATCC 13032 is a widely used producer of L-amino acids at an industrial scale. Its apparent robustness also turns it into a favorable platform host for a wide range of further compounds, mainly because of emerging bio-based economies. A deep understanding of the biochemical processes in C. glutamicum is essential for a sustainable enhancement of the microbe's productivity. Computational systems biology has the potential to provide a valuable basis for driving metabolic engineering and biotechnological advances, such as increased yields of healthy producer strains based on genome-scale metabolic models (GEMs). Advanced reconstruction pipelines are now available that facilitate the reconstruction of GEMs and support their manual curation. This article presents iCGB21FR, an updated and unified GEM of C. glutamicum ATCC 13032 with high quality regarding comprehensiveness and data standards, built with the latest modeling techniques and advanced reconstruction pipelines. It comprises 1042 metabolites, 1539 reactions, and 805 genes with detailed annotations and database cross-references. The model validation took place using different media and resulted in realistic growth rate predictions under aerobic and anaerobic conditions. The new GEM produces all canonical amino acids, and its phenotypic predictions are consistent with laboratory data. The in silico model proved fruitful in adding knowledge to the metabolism of C. glutamicum: iCGB21FR still produces L-glutamate with the knock-out of the enzyme pyruvate carboxylase, despite the common belief to be relevant for the amino acid's production. We conclude that integrating high standards into the reconstruction of GEMs facilitates replicating validated knowledge, closing knowledge gaps, and making it a useful basis for metabolic engineering. The model is freely available from BioModels Database under identifier MODEL2102050001.

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

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

OrtSuite: from genomes to prediction of microbial interactions within targeted ecosystem processes

OrtSuite predicts synergistic species interactions using the genomic potential of microbial communities

The high complexity found in microbial communities makes the identification of microbial interactions challenging. To address this challenge, we present OrtSuite, a flexible workflow to predict putative microbial interactions based on genomic content of microbial communities and targeted to specific ecosystem processes. The pipeline is composed of three user-friendly bash commands. OrtSuite combines ortholog clustering with genome annotation strategies limited to user-defined sets of functions allowing for hypothesis-driven data analysis such as assessing microbial interactions in specific ecosystems. OrtSuite matched, on average, 96% of experimentally verified KEGG orthologs involved in benzoate degradation in a known group of benzoate degraders. We evaluated the identification of putative synergistic species interactions using the sequenced genomes of an independent study that had previously proposed potential species interactions in benzoate degradation. OrtSuite is an easy-to-use workflow that allows for rapid functional annotation based on a user-curated database and can easily be extended to ecosystem processes where connections between genes and reactions are known. OrtSuite is an open-source software available at https://github.com/mdsufz/OrtSuite.

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.

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

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

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

Integrating Systems and Synthetic Biology to Understand and Engineer Microbiomes

Microbiomes are complex and ubiquitous networks of microorganisms whose seemingly limitless chemical transformations could be harnessed to benefit agriculture, medicine, and biotechnology. The spatial and temporal changes in microbiome composition and function are influenced by a multitude of molecular and ecological factors. This complexity yields both versatility and challenges in designing synthetic microbiomes and perturbing natural microbiomes in controlled, predictable ways. In this review, we describe factors that give rise to emergent spatial and temporal microbiome properties and the meta-omics and computational modeling tools that can be used to understand microbiomes at the cellular and system levels. We also describe strategies for designing and engineering microbiomes to enhance or build novel functions. Throughout the review, we discuss key knowledge and technology gaps for elucidating the networks and deciphering key control points for microbiome engineering, and highlight examples where multiple omics and modeling approaches can be integrated to address these gaps.

Modelling microbial communities: Harnessing consortia for biotechnological applications

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

Optimization of a modeling platform to predict oncogenes from genome-scale metabolic networks of non-small-cell lung cancers

Cancer cell dysregulations result in the abnormal regulation of cellular metabolic pathways. By simulating this metabolic reprogramming using constraint-based modeling approaches, oncogenes can be predicted, and this knowledge can be used in prognosis and treatment. We introduced a trilevel optimization problem describing metabolic reprogramming for inferring oncogenes. First, this study used RNA-Seq expression data of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) samples and their healthy counterparts to reconstruct tissue-specific genome-scale metabolic models and subsequently build the flux distribution pattern that provided a measure for the oncogene inference optimization problem for determining tumorigenesis. The platform detected 45 genes for LUAD and 84 genes for LUSC that lead to tumorigenesis. A high level of differentially expressed genes was not an essential factor for determining tumorigenesis. The platform indicated that pyruvate kinase (PKM), a well-known oncogene with a low level of differential gene expression in LUAD and LUSC, had the highest fitness among the predicted oncogenes based on computation. By contrast, pyruvate kinase L/R (PKLR), an isozyme of PKM, had a high level of differential gene expression in both cancers. Phosphatidylserine synthase 1 (PTDSS1), an oncogene in LUAD, was inferred to have a low level of differential gene expression, and overexpression could significantly reduce survival probability. According to the factor analysis, PTDSS1 characteristics were close to those of the template, but they were unobvious in LUSC. Angiotensin-converting enzyme 2 (ACE2) has recently garnered widespread interest as the SARS-CoV-2 virus receptor. Moreover, we determined that ACE2 is an oncogene of LUSC but not of LUAD. The platform developed in this study can identify oncogenes with low levels of differential expression and be used to identify potential therapeutic targets for cancer treatment.

MiMeNet: Exploring microbiome-metabolome relationships using neural networks

The advance in microbiome and metabolome studies has generated rich omics data revealing the involvement of the microbial community in host disease pathogenesis through interactions with their host at a metabolic level. However, the computational tools to uncover these relationships are just emerging. Here, we present MiMeNet, a neural network framework for modeling microbe-metabolite relationships. Using ten iterations of 10-fold cross-validation on three paired microbiome-metabolome datasets, we show that MiMeNet more accurately predicts metabolite abundances (mean Spearman correlation coefficients increase from 0.108 to 0.309, 0.276 to 0.457, and -0.272 to 0.264) and identifies more well-predicted metabolites (increase in the number of well-predicted metabolites from 198 to 366, 104 to 143, and 4 to 29) compared to state-of-art linear models for individual metabolite predictions. Additionally, we demonstrate that MiMeNet can group microbes and metabolites with similar interaction patterns and functions to illuminate the underlying structure of the microbe-metabolite interaction network, which could potentially shed light on uncharacterized metabolites through “Guilt by Association”. Our results demonstrated that MiMeNet is a powerful tool to provide insights into the causes of metabolic dysregulation in disease, facilitating future hypothesis generation at the interface of the microbiome and metabolomics.

The microbiome has shown to functionally interact with its host or environment at a metabolic level, however the exact nature of these interactions is not well understood. In addition, metabolic dysregulation caused by the microbiome is believed to contribute to the development of diseases such as inflammatory bowel disease, diabetes mellitus, and obesity. In this manuscript, we introduce a computational framework to integrate microbiome and metabolome data to uncover microbe-metabolite interactions in a data-driven manner. Our model uses neural networks to predict metabolite abundances from microbe abundances. The trained models are then used to derive microbe-metabolite feature scores, which are used for clustering microbes and metabolites into functional modules. These module-based interactions are useful in generating biological insights and facilitating hypothesis generation for the investigation of their roles in various metabolic diseases. The software of our model is made freely available to interested researchers.

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.

CODY enables quantitatively spatiotemporal predictions on in vivo gut microbial variability induced by diet intervention

Quantitatively understanding and predicting spatiotemporal dynamics of microbiota is imperative for development of tailored microbiome-directed therapeutics treatments. However, the complexity of microbial variations, due to interactions with the host, other microbes, and environmental factors, makes it challenging to identify how microbiota colonize in the human gut. Here, we describe a novel multiscale framework for COmputing the DYnamics of the gut microbiota (CODY), which enables the quantification of spatiotemporal-specific variations of gut microbiome abundance profiles, without a prior knowledge of microbiome interactions. Importantly, the predictive power of CODY is demonstrated using cross-sectional data from two longitudinal metagenomics studies—the microbiota development during early infancy and during short-term diet intervention of obese adults.

Microbial variations in the human gut are harbored in temporal and spatial heterogeneity, and quantitative prediction of spatiotemporal dynamic changes in the gut microbiota is imperative for development of tailored microbiome-directed therapeutics treatments, e.g. precision nutrition. Given the high-degree complexity of microbial variations, subject to the dynamic interactions among host, microbial, and environmental factors, identifying how microbiota colonize in the gut represents an important challenge. Here we present COmputing the DYnamics of microbiota (CODY), a multiscale framework that integrates species-level modeling of microbial dynamics and ecosystem-level interactions into a mathematical model that characterizes spatial-specific in vivo microbial residence in the colon as impacted by host physiology. The framework quantifies spatiotemporal resolution of microbial variations on species-level abundance profiles across site-specific colon regions and in feces, independent of a priori knowledge. We demonstrated the effectiveness of CODY using cross-sectional data from two longitudinal metagenomics studies—the microbiota development during early infancy and during short-term diet intervention of obese adults. For each cohort, CODY correctly predicts the microbial variations in response to diet intervention, as validated by available metagenomics and metabolomics data. Model simulations provide insight into the biogeographical heterogeneity among lumen, mucus, and feces, which provides insight into how host physical forces and spatial structure are shaping microbial structure and functionality.

Understanding the host-microbe interactions using metabolic modeling

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

Metabolic modelling approaches for describing and engineering microbial communities

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

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.

VFFVA: dynamic load balancing enables large-scale flux variability analysis

BACKGROUND

Genome-scale metabolic models are increasingly employed to predict the phenotype of various biological systems pertaining to healthcare and bioengineering. To characterize the full metabolic spectrum of such systems, Fast Flux Variability Analysis (FFVA) is commonly used in parallel with static load balancing. This approach assigns to each core an equal number of biochemical reactions without consideration of their solution complexity.

RESULTS

Here, we present Very Fast Flux Variability Analysis (VFFVA) as a parallel implementation that dynamically balances the computation load between the cores in runtime which guarantees equal convergence time between them. VFFVA allowed to gain a threefold speedup factor with coupled models and up to 100 with ill-conditioned models along with a 14-fold decrease in memory usage.

CONCLUSIONS

VFFVA exploits the parallel capabilities of modern machines to enable biological insights through optimizing systems biology modeling. VFFVA is available in C, MATLAB, and Python at https://github.com/marouenbg/VFFVA .

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

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

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

Oncogene inference optimization using constraint-based modelling incorporated with protein expression in normal and tumour tissues

Cancer cells are known to exhibit unusual metabolic activity, and yet few metabolic cancer driver genes are known. Genetic alterations and epigenetic modifications of cancer cells result in the abnormal regulation of cellular metabolic pathways that are different when compared with normal cells. Such a metabolic reprogramming can be simulated using constraint-based modelling approaches towards predicting oncogenes. We introduced the tri-level optimization problem to use the metabolic reprogramming towards inferring oncogenes. The algorithm incorporated Recon 2.2 network with the Human Protein Atlas to reconstruct genome-scale metabolic network models of the tissue-specific cells at normal and cancer states, respectively. Such reconstructed models were applied to build the templates of the metabolic reprogramming between normal and cancer cell metabolism. The inference optimization problem was formulated to use the templates as a measure towards predicting oncogenes. The nested hybrid differential evolution algorithm was applied to solve the problem to overcome solving difficulty for transferring the inner optimization problem into the single one. Head and neck squamous cells were applied as a case study to evaluate the algorithm. We detected 13 of the top-ranked one-hit dysregulations and 17 of the top-ranked two-hit oncogenes with high similarity ratios to the templates. According to the literature survey, most inferred oncogenes are consistent with the observation in various tissues. Furthermore, the inferred oncogenes were highly connected with the TP53/AKT/IGF/MTOR signalling pathway through PTEN, which is one of the most frequently detected tumour suppressor genes in human cancer.

Genome-Scale Metabolic Reconstruction of Acetobacter pasteurianus 386B, a Candidate Functional Starter Culture for Cocoa Bean Fermentation

Acetobacter pasteurianus 386B is a candidate functional starter culture for the cocoa bean fermentation process. To allow in silico simulations of its related metabolism in response to different environmental conditions, a genome-scale metabolic model for A. pasteurianus 386B was reconstructed. This is the first genome-scale metabolic model reconstruction for a member of the genus Acetobacter. The metabolic network reconstruction process was based on extensive genome re-annotation and comparative genomics analyses. The information content related to the functional annotation of metabolic enzymes and transporters was placed in a metabolic context by exploring and curating a Pathway/Genome Database of A. pasteurianus 386B using the Pathway Tools software. Metabolic reactions and curated gene-protein-reaction associations were bundled into a genome-scale metabolic model of A. pasteurianus 386B, named iAp386B454, containing 454 genes, 322 reactions, and 296 metabolites embedded in two cellular compartments. The reconstructed model was validated by performing growth experiments in a defined medium, which revealed that lactic acid as the sole carbon source could sustain growth of this strain. Further, the reconstruction of the A. pasteurianus 386B genome-scale metabolic model revealed knowledge gaps concerning the metabolism of this strain, especially related to the biosynthesis of its cell envelope and the presence or absence of metabolite transporters.

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.

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.

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.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.