Genome-scale Network Reconstructions Research Articles

Genome-scale metabolic modeling of the human milk oligosaccharide utilization by Bifidobacterium longum subsp. infantis.

Bifidobacterium longum subsp. infantis is a representative and dominant species in the infant gut and is considered a beneficial microbe. This organism displays multiple adaptations to thrive in the infant gut, regarded as a model for human milk oligosaccharides (HMOs) utilization. These carbohydrates are abundant in breast milk and include different molecules based on lactose. They contain fucose, sialic acid, and N-acetylglucosamine. Bifidobacterium metabolism is complex, and a systems view of relevant metabolic pathways and exchange metabolites during HMO consumption is missing. To address this limitation, a refined genome-scale network reconstruction of this bacterium is presented using a previous reconstruction of B. infantis ATCC 15967 as a template. The latter was expanded based on an extensive revision of genome annotations, current literature, and transcriptomic data integration. The metabolic reconstruction (iLR578) accounted for 578 genes, 1,047 reactions, and 924 metabolites. Starting from this reconstruction, we built context-specific genome-scale metabolic models using RNA-seq data from cultures growing in lactose and three HMOs. The models revealed notable differences in HMO metabolism depending on the functional characteristics of the substrates. Particularly, fucosyl-lactose showed a divergent metabolism due to a fucose moiety. High yields of lactate and acetate were predicted under growth rate maximization in all conditions, whereas formate, ethanol, and 1,2-propanediol were substantially lower. Similar results were also obtained under near-optimal growth on each substrate when varying the empirically observed acetate-to-lactate production ratio. Model predictions displayed reasonable agreement between central carbon metabolism fluxes and expression data across all conditions. Flux coupling analysis revealed additional connections between succinate exchange and arginine and sulfate metabolism and a strong coupling between central carbon reactions and adenine metabolism. More importantly, specific networks of coupled reactions under each carbon source were derived and analyzed. Overall, the presented network reconstruction constitutes a valuable platform for probing the metabolism of this prominent infant gut bifidobacteria.IMPORTANCEThis work presents a detailed reconstruction of the metabolism of Bifidobacterium longum subsp. infantis, a prominent member of the infant gut microbiome, providing a systems view of its metabolism of human milk oligosaccharides.

A computational model of Pseudomonas syringae metabolism unveils a role for branched-chain amino acids in Arabidopsis leaf colonization.

Bacterial pathogens adapt their metabolism to the plant environment to successfully colonize their hosts. In our efforts to uncover the metabolic pathways that contribute to the colonization of Arabidopsis thaliana leaves by Pseudomonas syringae pv tomato DC3000 (Pst DC3000), we created iPst19, an ensemble of 100 genome-scale network reconstructions of Pst DC3000 metabolism. We developed a novel approach for gene essentiality screens, leveraging the predictive power of iPst19 to identify core and ancillary condition-specific essential genes. Constraining the metabolic flux of iPst19 with Pst DC3000 gene expression data obtained from naïve-infected or pre-immunized-infected plants, revealed changes in bacterial metabolism imposed by plant immunity. Machine learning analysis revealed that among other amino acids, branched-chain amino acids (BCAAs) metabolism significantly contributed to the overall metabolic status of each gene-expression-contextualized iPst19 simulation. These predictions were tested and confirmed experimentally. Pst DC3000 growth and gene expression analysis showed that BCAAs suppress virulence gene expression in vitro without affecting bacterial growth. In planta, however, an excess of BCAAs suppress the expression of virulence genes at the early stages of infection and significantly impair the colonization of Arabidopsis leaves. Our findings suggesting that BCAAs catabolism is necessary to express virulence and colonize the host. Overall, this study provides valuable insights into how plant immunity impacts Pst DC3000 metabolism, and how bacterial metabolism impacts the expression of virulence.

Genome-scale metabolic network reconstructions of diverse Escherichia strains reveal strain-specific adaptations.

Bottom-up approaches to systems biology rely on constructing a mechanistic basis for the biochemical and genetic processes that underlie cellular functions. Genome-scale network reconstructions of metabolism are built from all known metabolic reactions and metabolic genes in a target organism. A network reconstruction can be converted into a mathematical format and thus lend itself to mathematical analysis. Genome-scale models (GEMs) of metabolism enable a systems approach to characterize the pan and core metabolic capabilities of the Escherichia genus. In this work, GEMs were constructed for 222 representative strains of Escherichia across HC1100 levels spanning the known Escherichia phylogeny. The models were used to study Escherichia metabolic diversity and speciation on a large scale. The results show that unique strain-specific metabolic capabilities correspond to different species and nutrient niches. This work is a first step towards a curated reconstruction of pan-Escherichia metabolism.This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

Pan-Genome-Scale Network Reconstruction: Harnessing Phylogenomics Increases the Quantity and Quality of Metabolic Models.

A genome-scale network reconstruction (GENRE) is a knowledgebase for an organism and has various applications. Available genome sequences have risen in recent years, but the number of curated GENREs has not kept pace. Existing yeast GENREs contain significant commission and omission errors. Current practices limit the quantity and quality of GENREs. An open and transparent phylogenomic-driven framework is outlined to address these issues. The method is demonstrated with 33 yeasts and fungi in Dikarya. A pan-fungal metabolic network called FYRMENT (Fungal and Yeast Metabolic Network) (https://github.com/LMSE/FYRMENT) is created, and annotated with ortholog groups from AYbRAH (https://github.com/LMSE/AYbRAH). Metabolic models for lower-level taxons are compiled. The fungal pan-GENRE contains 1553 orthologs, 2759 reactions, 2251 metabolites. The GENREs have higher genomic and metabolic coverage than existing yeast and fungal GENREs created with other methods. Metabolic simulations show the maximum amino acid yields from glucose differs between yeast lineages, indicating metabolic networks have evolved. Curating genomes and reactions at higher taxonomic-levels increases the quantity and quality of GENREs than conventional approaches. This approach can scale to other branches in the tree of life.

In-depth characterization of genome-scale network reconstructions for the in vitro synthesis in cell-free systems.

Cell-free systems containing multiple enzymes are becoming an increasingly interesting tool for one-pot syntheses of biochemical compounds. To extensively explore the enormous wealth of enzymes in the biological space, we present methods for assembling and curing data from databases to apply them for the prediction of pathway candidates for directed enzymatic synthesis. We use Kyoto Encyclopedia of Genes and Genomes to establish single organism models and a pan-organism model that is combining the available data from all organisms listed there. We introduce a filtering scheme to remove data that are not suitable, for example, generic metabolites and general reactions. In addition, a valid stoichiometry of reactions is required for acceptance. The networks created are analyzed by graph theoretical methods to identify a set of metabolites that are potentially reachable from a defined set of starting metabolites. Thus, metabolites not connected to such starting metabolites cannot be produced unless new starting metabolites or reactions are introduced. The network models also comprise stoichiometric and thermodynamic data that allow the definition of constraints to identify potential pathways. The resulting data can be directly applied using existing or future pathway finding tools.

A simplified metabolic network reconstruction to promote understanding and development of flux balance analysis tools

GEnome-scale Network REconstructions (GENREs) mathematically describe metabolic reactions of an organism or a specific cell type. GENREs can be used with a number of constraint-based reconstruction and analysis (COBRA) methods to make computational predictions on how a system changes in different environments. We created a simplified GENRE (referred to as iSIM) that captures central energy metabolism with nine metabolic reactions to illustrate the use of and promote the understanding of GENREs and constraint-based methods. We demonstrate the simulation of single and double gene deletions, flux variability analysis (FVA), and test a number of metabolic tasks with the GENRE. Code to perform these analyses is provided in Python, R, and MATLAB. Finally, with iSIM as a guide, we demonstrate how inaccuracies in GENREs can limit their use in the interrogation of energy metabolism.

Metabolic network-based predictions of toxicant-induced metabolite changes in the laboratory rat

In order to provide timely treatment for organ damage initiated by therapeutic drugs or exposure to environmental toxicants, we first need to identify markers that provide an early diagnosis of potential adverse effects before permanent damage occurs. Specifically, the liver, as a primary organ prone to toxicants-induced injuries, lacks diagnostic markers that are specific and sensitive to the early onset of injury. Here, to identify plasma metabolites as markers of early toxicant-induced injury, we used a constraint-based modeling approach with a genome-scale network reconstruction of rat liver metabolism to incorporate perturbations of gene expression induced by acetaminophen, a known hepatotoxicant. A comparison of the model results against the global metabolic profiling data revealed that our approach satisfactorily predicted altered plasma metabolite levels as early as 5 h after exposure to 2 g/kg of acetaminophen, and that 10 h after treatment the predictions significantly improved when we integrated measured central carbon fluxes. Our approach is solely driven by gene expression and physiological boundary conditions, and does not rely on any toxicant-specific model component. As such, it provides a mechanistic model that serves as a first step in identifying a list of putative plasma metabolites that could change due to toxicant-induced perturbations.

Recon3D enables a three-dimensional view of gene variation in human metabolism.

Genome-scale network reconstructions have helped uncover the molecular basis of metabolism. Here we present Recon3D, a computational resource that includes three-dimensional (3D) metabolite and protein structure data and enables integrated analyses of metabolic functions in humans. We use Recon3D to functionally characterize mutations associated with disease, and identify metabolic response signatures that are caused by exposure to certain drugs. Recon3D represents the most comprehensive human metabolic network model to date, accounting for 3,288 open reading frames (representing 17% of functionally annotated human genes), 13,543 metabolic reactions involving 4,140 unique metabolites, and 12,890 protein structures. These data provide a unique resource for investigating molecular mechanisms of human metabolism. Recon3D is available at http://vmh.life.

Ssbio: a Python framework for structural systems biology.

Working with protein structures at the genome-scale has been challenging in a variety of ways. Here, we present ssbio, a Python package that provides a framework to easily work with structural information in the context of genome-scale network reconstructions, which can contain thousands of individual proteins. The ssbio package provides an automated pipeline to construct high quality genome-scale models with protein structures (GEM-PROs), wrappers to popular third-party programs to compute associated protein properties, and methods to visualize and annotate structures directly in Jupyter notebooks, thus lowering the barrier of linking 3D structural data with established systems workflows. ssbio is implemented in Python and available to download under the MIT license at http://github.com/SBRG/ssbio. Documentation and Jupyter notebook tutorials are available at http://ssbio.readthedocs.io/en/latest/. Interactive notebooks can be launched using Binder at https://mybinder.org/v2/gh/SBRG/ssbio/master?filepath=Binder.ipynb. Supplementary data are available at Bioinformatics online.

Reconciled rat and human metabolic networks for comparative toxicogenomics and biomarker predictions

The laboratory rat has been used as a surrogate to study human biology for more than a century. Here we present the first genome-scale network reconstruction of Rattus norvegicus metabolism, iRno, and a significantly improved reconstruction of human metabolism, iHsa. These curated models comprehensively capture metabolic features known to distinguish rats from humans including vitamin C and bile acid synthesis pathways. After reconciling network differences between iRno and iHsa, we integrate toxicogenomics data from rat and human hepatocytes, to generate biomarker predictions in response to 76 drugs. We validate comparative predictions for xanthine derivatives with new experimental data and literature-based evidence delineating metabolite biomarkers unique to humans. Our results provide mechanistic insights into species-specific metabolism and facilitate the selection of biomarkers consistent with rat and human biology. These models can serve as powerful computational platforms for contextualizing experimental data and making functional predictions for clinical and basic science applications.

A mixed-integer linear programming approach to the reduction of genome-scale metabolic networks

BackgroundConstraint-based analysis has become a widely used method to study metabolic networks. While some of the associated algorithms can be applied to genome-scale network reconstructions with several thousands of reactions, others are limited to small or medium-sized models. In 2015, Erdrich et al. introduced a method called NetworkReducer, which reduces large metabolic networks to smaller subnetworks, while preserving a set of biological requirements that can be specified by the user. Already in 2001, Burgard et al. developed a mixed-integer linear programming (MILP) approach for computing minimal reaction sets under a given growth requirement.ResultsHere we present an MILP approach for computing minimum subnetworks with the given properties. The minimality (with respect to the number of active reactions) is not guaranteed by NetworkReducer, while the method by Burgard et al. does not allow specifying the different biological requirements. Our procedure is about 5-10 times faster than NetworkReducer and can enumerate all minimum subnetworks in case there exist several ones. This allows identifying common reactions that are present in all subnetworks, and reactions appearing in alternative pathways.ConclusionsApplying complex analysis methods to genome-scale metabolic networks is often not possible in practice. Thus it may become necessary to reduce the size of the network while keeping important functionalities. We propose a MILP solution to this problem. Compared to previous work, our approach is more efficient and allows computing not only one, but even all minimum subnetworks satisfying the required properties.

Thermodynamics and H2 Transfer in a Methanogenic, Syntrophic Community.

Microorganisms in nature do not exist in isolation but rather interact with other species in their environment. Some microbes interact via syntrophic associations, in which the metabolic by-products of one species serve as nutrients for another. These associations sustain a variety of natural communities, including those involved in methanogenesis. In anaerobic syntrophic communities, energy is transferred from one species to another, either through direct contact and exchange of electrons, or through small molecule diffusion. Thermodynamics plays an important role in governing these interactions, as the oxidation reactions carried out by the first community member are only possible because degradation products are consumed by the second community member. This work presents the development and analysis of genome-scale network reconstructions of the bacterium Syntrophobacter fumaroxidans and the methanogenic archaeon Methanospirillum hungatei. The models were used to verify proposed mechanisms of ATP production within each species. We then identified additional constraints and the cellular objective function required to match experimental observations. The thermodynamic S. fumaroxidans model could not explain why S. fumaroxidans does not produce H2 in monoculture, indicating that current methods might not adequately estimate the thermodynamics, or that other cellular processes (e.g., regulation) play a role. We also developed a thermodynamic coculture model of the association between the organisms. The coculture model correctly predicted the exchange of both H2 and formate between the two species and suggested conditions under which H2 and formate produced by S. fumaroxidans would be fully consumed by M. hungatei.

Computing autocatalytic sets to unravel inconsistencies in metabolic network reconstructions.

Genome-scale metabolic network reconstructions have been established as a powerful tool for the prediction of cellular phenotypes and metabolic capabilities of organisms. In recent years, the number of network reconstructions has been constantly increasing, mostly because of the availability of novel (semi-)automated procedures, which enabled the reconstruction of metabolic models based on individual genomes and their annotation. The resulting models are widely used in numerous applications. However, the accuracy and predictive power of network reconstructions are commonly limited by inherent inconsistencies and gaps. Here we present a novel method to validate metabolic network reconstructions based on the concept of autocatalytic sets. Autocatalytic sets correspond to collections of metabolites that, besides enzymes and a growth medium, are required to produce all biomass components in a metabolic model. These autocatalytic sets are well-conserved across all domains of life, and their identification in specific genome-scale reconstructions allows us to draw conclusions about potential inconsistencies in these models. The method is capable of detecting inconsistencies, which are neglected by other gap-finding methods. We tested our method on the Model SEED, which is the largest repository for automatically generated genome-scale network reconstructions. In this way, we were able to identify a significant number of missing pathways in several of these reconstructions. Hence, the method we report represents a powerful tool to identify inconsistencies in large-scale metabolic networks. The method is available as source code on http://users.minet.uni-jena.de/∼m3kach/ASBIG/ASBIG.zip. christoph.kaleta@uni-jena.de Supplementary data are available at Bioinformatics online.

Optimizing genome-scale network reconstructions.

Metabolic reconstructions remain limited in their scope and content, and improvements in biochemical knowledge and collaborative research are required.

Quantitative Assessment of Thermodynamic Constraints on the Solution Space of Genome-Scale Metabolic Models

Constraint-based methods provide powerful computational techniques to allow understanding and prediction of cellular behavior. These methods rely on physiochemical constraints to eliminate infeasible behaviors from the space of available behaviors. One such constraint is thermodynamic feasibility, the requirement that intracellular flux distributions obey the laws of thermodynamics. The past decade has seen several constraint-based methods that interpret this constraint in different ways, including those that are limited to small networks, rely on predefined reaction directions, and/or neglect the relationship between reaction free energies and metabolite concentrations. In this work, we utilize one such approach, thermodynamics-based metabolic flux analysis (TMFA), to make genome-scale, quantitative predictions about metabolite concentrations and reaction free energies in the absence of prior knowledge of reaction directions, while accounting for uncertainties in thermodynamic estimates. We applied TMFA to a genome-scale network reconstruction of Escherichia coli and examined the effect of thermodynamic constraints on the flux space. We also assessed the predictive performance of TMFA against gene essentiality and quantitative metabolomics data, under both aerobic and anaerobic, and optimal and suboptimal growth conditions. Based on these results, we propose that TMFA is a useful tool for validating phenotypes and generating hypotheses, and that additional types of data and constraints can improve predictions of metabolite concentrations.

An Integrative Bioinformatics Framework for Genome-scale Multiple Level Network Reconstruction of Rice

Summary Understanding how metabolic reactions translate the genome of an organism into its phenotype is a grand challenge in biology. Genome-wide association studies (GWAS) statistically connect genotypes to phenotypes, without any recourse to known molecular interactions, whereas a molecular mechanistic description ties gene function to phenotype through gene regulatory networks (GRNs), protein-protein interactions (PPIs) and molecular pathways. Integration of different regulatory information levels of an organism is expected to provide a good way for mapping genotypes to phenotypes. However, the lack of curated metabolic model of rice is blocking the exploration of genome-scale multi-level network reconstruction. Here, we have merged GRNs, PPIs and genome-scale metabolic networks (GSMNs) approaches into a single framework for rice via omics’ regulatory information reconstruction and integration. Firstly, we reconstructed a genome-scale metabolic model, containing 4,462 function genes, 2,986 metabolites involved in 3,316 reactions, and compartmentalized into ten subcellular locations. Furthermore, 90,358 pairs of protein-protein interactions, 662,936 pairs of gene regulations and 1,763 microRNA-target interactions were integrated into the metabolic model. Eventually, a database was developped for systematically storing and retrieving the genome-scale multi-level network of rice. This provides a reference for understanding genotype-phenotype relationship of rice, and for analysis of its molecular regulatory network.

Mapping condition-dependent regulation of metabolism in yeast through genome-scale modeling

BackgroundThe genome-scale metabolic model of Saccharomyces cerevisiae, first presented in 2003, was the first genome-scale network reconstruction for a eukaryotic organism. Since then continuous efforts have been made in order to improve and expand the yeast metabolic network.ResultsHere we present iTO977, a comprehensive genome-scale metabolic model that contains more reactions, metabolites and genes than previous models. The model was constructed based on two earlier reconstructions, namely iIN800 and the consensus network, and then improved and expanded using gap-filling methods and by introducing new reactions and pathways based on studies of the literature and databases. The model was shown to perform well both for growth simulations in different media and gene essentiality analysis for single and double knock-outs. Further, the model was used as a scaffold for integrating transcriptomics, and flux data from four different conditions in order to identify transcriptionally controlled reactions, i.e. reactions that change both in flux and transcription between the compared conditions.ConclusionWe present a new yeast model that represents a comprehensive up-to-date collection of knowledge on yeast metabolism. The model was used for simulating the yeast metabolism under four different growth conditions and experimental data from these four conditions was integrated to the model. The model together with experimental data is a useful tool to identify condition-dependent changes of metabolism between different environmental conditions.

An integrative bioinformatics framework for genome-scale multiple level network reconstruction of rice.

Understanding how metabolic reactions translate the genome of an organism into its phenotype is a grand challenge in biology. Genome-wide association studies (GWAS) statistically connect genotypes to phenotypes, without any recourse to known molecular interactions, whereas a molecular mechanistic description ties gene function to phenotype through gene regulatory networks (GRNs), protein-protein interactions (PPIs) and molecular pathways. Integration of different regulatory information levels of an organism is expected to provide a good way for mapping genotypes to phenotypes. However, the lack of curated metabolic model of rice is blocking the exploration of genome-scale multi-level network reconstruction. Here, we have merged GRNs, PPIs and genome-scale metabolic networks (GSMNs) approaches into a single framework for rice via omics’ regulatory information reconstruction and integration. Firstly, we reconstructed a genome-scale metabolic model, containing 4,462 function genes, 2,986 metabolites involved in 3,316 reactions, and compartmentalized into ten subcellular locations. Furthermore, 90,358 pairs of protein-protein interactions, 662,936 pairs of gene regulations and 1,763 microRNA-target interactions were integrated into the metabolic model. Eventually, a database was developped for systematically storing and retrieving the genome-scale multi-level network of rice. This provides a reference for understanding genotype-phenotype relationship of rice, and for analysis of its molecular regulatory network.

Genome-scale bacterial transcriptional regulatory networks: reconstruction and integrated analysis with metabolic models

Advances in sequencing technology are resulting in the rapid emergence of large numbers of complete genome sequences. High-throughput annotation and metabolic modeling of these genomes is now a reality. The high-throughput reconstruction and analysis of genome-scale transcriptional regulatory networks represent the next frontier in microbial bioinformatics. The fruition of this next frontier will depend on the integration of numerous data sources relating to mechanisms, components and behavior of the transcriptional regulatory machinery, as well as the integration of the regulatory machinery into genome-scale cellular models. Here, we review existing repositories for different types of transcriptional regulatory data, including expression data, transcription factor data and binding site locations and we explore how these data are being used for the reconstruction of new regulatory networks. From template network-based methods to de novo reverse engineering from expression data, we discuss how regulatory networks can be reconstructed and integrated with metabolic models to improve model predictions and performance. We also explore the impact these integrated models can have in simulating phenotypes, optimizing the production of compounds of interest or paving the way to a whole-cell model.

Integrating Cellular Metabolism into a Multiscale Whole-Body Model

Cellular metabolism continuously processes an enormous range of external compounds into endogenous metabolites and is as such a key element in human physiology. The multifaceted physiological role of the metabolic network fulfilling the catalytic conversions can only be fully understood from a whole-body perspective where the causal interplay of the metabolic states of individual cells, the surrounding tissue and the whole organism are simultaneously considered. We here present an approach relying on dynamic flux balance analysis that allows the integration of metabolic networks at the cellular scale into standardized physiologically-based pharmacokinetic models at the whole-body level. To evaluate our approach we integrated a genome-scale network reconstruction of a human hepatocyte into the liver tissue of a physiologically-based pharmacokinetic model of a human adult. The resulting multiscale model was used to investigate hyperuricemia therapy, ammonia detoxification and paracetamol-induced toxication at a systems level. The specific models simultaneously integrate multiple layers of biological organization and offer mechanistic insights into pathology and medication. The approach presented may in future support a mechanistic understanding in diagnostics and drug development.