Modeling the differences in biochemical capabilities of pseudomonas species by flux balance analysis: how good are genome-scale metabolic networks at predicting the differences?

Bacteria and Methanogens in the Human Microbiome: a Review of Syntrophic Interactions

Methanogens are microorganisms belonging to the Archaea domain and represent the primary source of biotic methane. Methanogens encode a series of enzymes which can convert secondary substrates into methane following three major methanogenesis pathways. Initially recognized as environmental microorganisms, methanogens have more recently been acknowledged as host-associated microorganisms after their detection and initial isolation in ruminants in the 1950s. Methanogens have also been co-detected with bacteria in various pathological situations, bringing their role as pathogens into question. Here, we review reported associations between methanogens and bacteria in physiological and pathological situations in order to understand the metabolic interactions explaining these associations. To do so, we describe the origin of the metabolites used for methanogenesis and highlight the central role of methanogens in the syntrophic process during carbon cycling. We then focus on the metabolic abilities of co-detected bacterial species described in the literature and infer from their genomes the probable mechanisms of their association with methanogens. The syntrophic interactions between bacteria and methanogens are paramount to gut homeostasis. Therefore, any dysbiosis affecting methanogens might impact human health. Thus, the monitoring of methanogens may be used as a bio-indicator of dysbiosis. Moreover, new therapeutic approaches can be developed based on their administration as probiotics. We thus insist on the importance of investigating methanogens in clinical microbiology.

Genome-Scale Metabolic Network Models of Bacillus Species Suggest that Model Improvement is Necessary for Biotechnological Applications

Background

A genome-scale metabolic network model (GEM) is a mathematical representation of an organism’s metabolism. Today, GEMs are popular tools for computationally simulating the biotechnological processes and for predicting biochemical properties of (engineered) strains.

Objectives

In the present study, we have evaluated the predictive power of two GEMs, namely iBsu1103 (for Bacillus subtilis 168) and iMZ1055 (for Bacillus megaterium WSH002).

Materials and Methods

For comparing the predictive power of Bacillus subtilis and Bacillus megaterium GEMs, experimental data were obtained from previous wet-lab studies included in PubMed. By using these data, we set the environmental, stoichiometric and thermodynamic constraints on the models, and FBA is performed to predict the biomass production rate, and the values of other fluxes. For simulating experimental conditions in this study, COBRA toolbox was used.

Results

By using the wealth of data in the literature, we evaluated the accuracy of in silico simulations of these GEMs. Our results suggest that there are some errors in these two models which make them unreliable for predicting the biochemical capabilities of these species. The inconsistencies between experimental and computational data are even greater where B. subtilis and B. megaterium do not have similar phenotypes.

Conclusions

Our analysis suggests that literature-based improvement of genome-scale metabolic network models of the two Bacillus species is essential if these models are to be successfully applied in biotechnology and metabolic engineering.

Aldehyde dehydrogenase diversity in bacteria of the Pseudomonas genus

Aldehyde dehydrogenases (ALDHs) comprise one of the most ancient protein superfamilies widely distributed in the three domains of life. Their members have been extensively studied in animals and plants, sorted out in different ALDH protein families and their participation in a broad variety of metabolic pathways has been documented. Paradoxically, no systematic studies comprising ALDHs from bacteria have been performed so far. Among bacteria, the genus Pseudomonas occupies numerous ecological niches, and is one of the most complex bacterial genera with the largest number of known species. For these reasons, we selected Pseudomonas as a paradigm to analyze the diversity of ALDHs in bacteria. With this aim, complete Pseudomonas genome sequences and annotations were retrieved from NCBI's RefSeq genome database. The 258 Pseudomonas strains belong to 46 different species, along with 23 with no species designation. The genomes of these Pseudomonas contain from 3,315 to 6,825 annotated protein coding genes. A total of 6,510 ALDH sequences were found in the selected Pseudomonas, with a median of 24 ALDH-coding genes per strain (by comparison humans possess only 19 different ALDH loci). Pseudomonas saudiphocaensis possesses the lowest number of aldh genes (9), while Pseudomonas pseudoalcaligenes KF707 NBRC110670 possesses the maximum number of aldh genes (49). The ALDHs found in Pseudomonas can be sorted out into 42 protein families, with a predominance of 14 families, which contained 76% of all ALDHs found. In this regard, it is important to note that many Pseudomonas genomes have multiple aldh genes coding for proteins belonging to the same family. Given that all strains contained members of families ALDH4, ALDH5, ALDH6, ALDH14, ALDH18 and ALDH27, we consider these families to be part of the core Pseudomonas genome.

CiliateGEM: an open-project and a tool for predictions of ciliate metabolic variations and experimental condition design

BACKGROUND

The study of cell metabolism is becoming central in several fields such as biotechnology, evolution/adaptation and human disease investigations. Here we present CiliateGEM, the first metabolic network reconstruction draft of the freshwater ciliate Tetrahymena thermophila. We also provide the tools and resources to simulate different growth conditions and to predict metabolic variations. CiliateGEM can be extended to other ciliates in order to set up a meta-model, i.e. a metabolic network reconstruction valid for all ciliates. Ciliates are complex unicellular eukaryotes of presumably monophyletic origin, with a phylogenetic position that is equal from plants and animals. These cells represent a new concept of unicellular system with a high degree of species, population biodiversity and cell complexity. Ciliates perform in a single cell all the functions of a pluricellular organism, including locomotion, feeding, digestion, and sexual processes.

RESULTS

After generating the model, we performed an in-silico simulation with the presence and absence of glucose. The lack of this nutrient caused a 32.1% reduction rate in biomass synthesis. Despite the glucose starvation, the growth did not stop due to the use of alternative carbon sources such as amino acids.

CONCLUSIONS

The future models obtained from CiliateGEM may represent a new approach to describe the metabolism of ciliates. This tool will be a useful resource for the ciliate research community in order to extend these species as model organisms in different research fields. An improved understanding of ciliate metabolism could be relevant to elucidate the basis of biological phenomena like genotype-phenotype relationships, population genetics, and cilia-related disease mechanisms.

A Genome-Scale Metabolic Reconstruction of Phytophthora infestans With the Integration of Transcriptional Data Reveals the Key Metabolic Patterns Involved in the Interaction of Its Host

Phytophthora infestans, the causal agent of late blight disease, affects potatoes and tomatoes worldwide. This plant pathogen has a hemibiotrophic lifestyle, having an initial biotrophic infection phase during which the pathogen spreads within the host tissue, followed by a necrotrophic phase in which host cell death is induced. Although increasing information is available on the molecular mechanisms, underlying the distinct phases of the hemibiotrophic lifestyle, studies that consider the entire metabolic processes in the pathogen while undergoing the biotrophic, transition to necrotrophic, and necrotrophic phases have not been conducted. In this study, the genome-scale metabolic reconstruction of P. infestans was achieved. Subsequently, transcriptional data (microarrays, RNA-seq) was integrated into the metabolic reconstruction to obtain context-specific (metabolic) models (CSMs) of the infection process, using constraint-based reconstruction and analysis. The goal was to identify specific metabolic markers for distinct stages of the pathogen's life cycle. Results indicate that the overall metabolism show significant changes during infection. The most significant changes in metabolism were observed at the latest time points of infection. Metabolic activity associated with purine, pyrimidine, fatty acid, fructose and mannose, arginine, glycine, serine, and threonine amino acids appeared to be the most important metabolisms of the pathogen during the course of the infection, showing high number of reactions associated with them and expression switches at important stages of the life cycle. This study provides a framework for future throughput studies of the metabolic changes during the hemibiotrophic life cycle of this important plant pathogen.

In silico engineering of Pseudomonas metabolism reveals new biomarkers for increased biosurfactant production

BACKGROUND

Rhamnolipids, biosurfactants with a wide range of biomedical applications, are amphiphilic molecules produced on the surfaces of or excreted extracellularly by bacteria including Pseudomonas aeruginosa. However, Pseudomonas putida is a non-pathogenic model organism with greater metabolic versatility and potential for industrial applications.

METHODS

We investigate in silico the metabolic capabilities of P. putida for rhamnolipids biosynthesis using statistical, metabolic and synthetic engineering approaches after introducing key genes (RhlA and RhlB) from P. aeruginosa into a genome-scale model of P. putida. This pipeline combines machine learning methods with multi-omic modelling, and drives the engineered P. putida model toward an optimal production and export of rhamnolipids out of the membrane.

RESULTS

We identify a substantial increase in synthesis of rhamnolipids by the engineered model compared to the control model. We apply statistical and machine learning techniques on the metabolic reaction rates to identify distinct features on the structure of the variables and individual components driving the variation of growth and rhamnolipids production. We finally provide a computational framework for integrating multi-omics data and identifying latent pathways and genes for the production of rhamnolipids in P. putida.

CONCLUSIONS

We anticipate that our results will provide a versatile methodology for integrating multi-omics data for topological and functional analysis of P. putida toward maximization of biosurfactant production.

Efficient Reconstruction of Predictive Consensus Metabolic Network Models

Understanding cellular function requires accurate, comprehensive representations of metabolism. Genome-scale, constraint-based metabolic models (GSMs) provide such representations, but their usability is often hampered by inconsistencies at various levels, in particular for concurrent models. COMMGEN, our tool for COnsensus Metabolic Model GENeration, automatically identifies inconsistencies between concurrent models and semi-automatically resolves them, thereby contributing to consolidate knowledge of metabolic function. Tests of COMMGEN for four organisms showed that automatically generated consensus models were predictive and that they substantially increased coherence of knowledge representation. COMMGEN ought to be particularly useful for complex scenarios in which manual curation does not scale, such as for eukaryotic organisms, microbial communities, and host-pathogen interactions.

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures-Systems biology-based interpretation using genome-scale metabolic flux balance model and multivariate data analysis

Genome-scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745-753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome-scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163-1173, 2016.