Search

Facilitative interactions between microbial species are ubiquitous in various types of ecosystems on the Earth. Therefore, inferring how entangled webs of interspecific interactions shift through time in microbial ecosystems is an essential step for understanding ecological processes driving microbiome dynamics. By compiling shotgun metagenomic sequencing data of an experimental microbial community, we examined how the architectural features of facilitative interaction networks could change through time. A metabolic modeling approach for estimating dependence between microbial genomes (species) allowed us to infer the network structure of potential facilitative interactions at 13 time points through the 110-day monitoring of experimental microbiomes. We then found that positive feedback loops, which were theoretically predicted to promote cascade breakdown of ecological communities, existed within the inferred networks of metabolic interactions prior to the drastic community-compositional shift observed in the microbiome time-series. We further applied “directed-graph” analyses to pinpoint potential keystone species located at the “upper stream” positions of such feedback loops. These analyses on facilitative interactions will help us understand key mechanisms causing catastrophic shifts in microbial community structure. [ABSTRACT FROM AUTHOR]

The bacterial microbiota of plants is diverse, with 1000s of operational taxonomic units (OTUs) associated with any individual plant. In this work, we used phenotypic analysis, comparative genomics, and metabolic models to investigate the differences between 19 sequenced Pseudomonas fluorescens strains. These isolates represent a single OTU and were collected from the rhizosphere and endosphere of Populus deltoides. While no traits were exclusive to either endosphere or rhizosphere P. fluorescens isolates, multiple pathways relevant for plant-bacterial interactions are enriched in endosphere isolate genomes. Further, growth phenotypes such as phosphate solubilization, protease activity, denitrification and root growth promotion are biased toward endosphere isolates. Endosphere isolates have significantly more metabolic pathways for plant signaling compounds and an increased metabolic range that includes utilization of energy rich nucleotides and sugars, consistent with endosphere colonization. Rhizosphere P. fluorescens have fewer pathways representative of plant-bacterial interactions but show metabolic bias toward chemical substrates often found in root exudates. This work reveals the diverse functions that may contribute to colonization of the endosphere by bacteria and are enriched among closely related isolates.

Present in many industrial effluents and as common degradation product of organic matter, phenol is a widespread compound which may cause serious environmental problems, due to its toxicity to animals and humans. Degradation of phenol from the environment by mesophilic bacteria has been studied extensively over the past decades, but only little is known about phenol biodegradation at high temperatures or low pH. In this work we studied phenol degradation in the thermoacidophilic archaeon Saccharolobus solfataricus P2 (basonym: Sulfolobus solfataricus) under extreme conditions (80°C, pH 3.5). We combined metabolomics and transcriptomics together with metabolic modeling to elucidate the organism’s response to growth with phenol as sole carbon source. Although S. solfataricus is able to utilize phenol for biomass production, the carbon source induces profound stress reactions, including genome rearrangement as well as a strong intracellular accumulation of polyamines. Furthermore, computational modeling revealed a 40% higher oxygen demand for substrate oxidation, compared to growth on glucose. However, only 16.5% of oxygen is used for oxidation of phenol to catechol, resulting in a less efficient integration of carbon into the biomass. Finally, our data underlines the importance of the phenol meta-degradation pathway in S. solfataricus and enables us to predict enzyme candidates involved in the degradation processes downstream of 2-hydroxymucconic acid.

Present in many industrial effluents and as common degradation product of organic matter, phenol is a widespread compound which may cause serious environmental problems, due to its toxicity to animals and humans. Degradation of phenol from the environment by mesophilic bacteria has been studied extensively over the past decades, but only little is known about phenol biodegradation at high temperatures or low pH. In this work we studied phenol degradation in the thermoacidophilic archaeon Saccharolobus solfataricus P2 (basonym: Sulfolobus solfataricus) under extreme conditions (80°C, pH 3.5). We combined metabolomics and transcriptomics together with metabolic modeling to elucidate the organism's response to growth with phenol as sole carbon source. Although S. solfataricus is able to utilize phenol for biomass production, the carbon source induces profound stress reactions, including genome rearrangement as well as a strong intracellular accumulation of polyamines. Furthermore, computational modeling revealed a 40% higher oxygen demand for substrate oxidation, compared to growth on glucose. However, only 16.5% of oxygen is used for oxidation of phenol to catechol, resulting in a less efficient integration of carbon into the biomass. Finally, our data underlines the importance of the phenol meta -degradation pathway in S. solfataricus and enables us to predict enzyme candidates involved in the degradation processes downstream of 2-hydroxymucconic acid. [ABSTRACT FROM AUTHOR]

Salt tolerance in the γ-proteobacterium Halomonas elongata is linked to its ability to produce the compatible solute ectoine. The metabolism of ectoine production is of great interest since it can shed light on the biochemical basis of halotolerance as well as pave the way for the improvement of the biotechnological production of such compatible solute. Ectoine belongs to the biosynthetic family of aspartate-derived amino-acids. Aspartate is formed from oxaloacetate, thereby connecting ectoine production to the anaplerotic reactions that refill carbon into the tricarboxylic acid cycle (TCA cycle). This places a high demand on these reactions and creates the need to regulate them not only in response to growth but also in response to extracellular salt concentration. In this work, we combine modeling and experiments to analyze how these different needs shape the anaplerotic reactions in H. elongata. First, the stoichiometric and thermodynamic factors that condition the flux distributions are analyzed, then the optimal patterns of operation for oxaloacetate production are calculated. Finally, the phenotype of two deletion mutants lacking potentially relevant anaplerotic enzymes: phosphoenolpyruvate carboxylase (Ppc) and oxaloacetate decarboxylase (Oad) are experimentally characterized. The results show that the anaplerotic reactions in H. elongata are indeed subject to evolutionary pressures that differ from those faced by other gram-negative bacteria. Ectoine producing halophiles must meet a higher metabolic demand for oxaloacetate and the reliance of many marine bacteria on the Entner-Doudoroff pathway compromises the anaplerotic efficiency of Ppc, which is usually one of the main enzymes fulfilling this role. The anaplerotic flux in H. elongata is contributed not only by Ppc but also by Oad, an enzyme that has not yet been shown to play this role in vivo. Ppc is necessary for H. elongata to grow normally at low salt concentrations but it is not required to achieve near maximal growth rates as long as there is a steep sodium gradient. On the other hand, the lack of Oad presents serious difficulties to grow at high salt concentrations. This points to a shared role of these two enzymes in guaranteeing the supply of oxaloacetate for biosynthetic reactions.

Metabolic adaptation to the host cell is important for obligate intracellular pathogens such as Chlamydia trachomatis (Ct). Here we infer the flux differences for Ct from proteome and qRT-PCR data by comprehensive pathway modeling. We compare the comparatively inert infectious elementary body (EB) and the active replicative reticulate body (RB) systematically using a genome-scale metabolic model with 321 metabolites and 277 reactions. This did yield 84 extreme pathways based on a published proteomics dataset at three different time points of infection. Validation of predictions was done by quantitative RT-PCR of enzyme mRNA expression at three time points. Ct’s major active pathways are glycolysis, gluconeogenesis, glycerol-phospholipid (GPL) biosynthesis (support from host acetyl-CoA) and pentose phosphate pathway (PPP), while its incomplete TCA and fatty acid biosynthesis are less active. The modeled metabolic pathways are much more active in RB than in EB. Our in silico model suggests that EB and RB utilize folate to generate NAD(P)H using independent pathways. The only low metabolic flux inferred for EB involves mainly carbohydrate metabolism. RB utilizes energy -rich compounds to generate ATP in nucleic acid metabolism. Validation data for the modeling include proteomics experiments (model basis) as well as qRT-PCR confirmation of selected metabolic enzyme mRNA expression differences. The metabolic modeling is made fully available here. Its detailed insights and models on Ct metabolic adaptations during infection are a useful modeling basis for future studies.

Salt tolerance in the γ-proteobacterium Halomonas elongata is linked to its ability to produce the compatible solute ectoine. The metabolism of ectoine production is of great interest since it can shed light on the biochemical basis of halotolerance as well as pave the way for the improvement of the biotechnological production of such compatible solute. Ectoine belongs to the biosynthetic family of aspartate-derived amino-acids. Aspartate is formed from oxaloacetate, thereby connecting ectoine production to the anaplerotic reactions that refill carbon into the tricarboxylic acid cycle (TCA cycle). This places a high demand on these reactions and creates the need to regulate them not only in response to growth but also in response to extracellular salt concentration. In this work, we combine modeling and experiments to analyze how these different needs shape the anaplerotic reactions in H. elongata. First, the stoichiometric and thermodynamic factors that condition the flux distributions are analyzed, then the optimal patterns of operation for oxaloacetate production are calculated. Finally, the phenotype of two deletion mutants lacking potentially relevant anaplerotic enzymes: phosphoenolpyruvate carboxylase (Ppc) and oxaloacetate decarboxylase (Oad) are experimentally characterized. The results show that the anaplerotic reactions in H. elongata are indeed subject to evolutionary pressures that differ from those faced by other gram-negative bacteria. Ectoine producing halophiles must meet a higher metabolic demand for oxaloacetate and the reliance of many marine bacteria on the Entner-Doudoroff pathway compromises the anaplerotic efficiency of Ppc, which is usually one of the main enzymes fulfilling this role. The anaplerotic flux in H. elongata is contributed not only by Ppc but also by Oad, an enzyme that has not yet been shown to play this role in vivo. Ppc is necessary for H. elongata to grow normally at low salt concentrations but it is not required to achieve near maximal growth rates as long as there is a steep sodium gradient. On the other hand, the lack of Oad presents serious difficulties to grow at high salt concentrations. This points to a shared role of these two enzymes in guaranteeing the supply of oxaloacetate for biosynthetic reactions. [ABSTRACT FROM AUTHOR]

Metabolic adaptation to the host cell is important for obligate intracellular pathogens such as Chlamydia trachomatis (Ct). Here we infer the flux differences for Ct from proteome and qRT-PCR data by comprehensive pathway modeling. We compare the comparatively inert infectious elementary body (EB) and the active replicative reticulate body (RB) systematically using a genome-scale metabolic model with 321 metabolites and 277 reactions. This did yield 84 extreme pathways based on a published proteomics dataset at three different time points of infection. Validation of predictions was done by quantitative RT-PCR of enzyme mRNA expression at three time points. Ct 's major active pathways are glycolysis, gluconeogenesis, glycerol-phospholipid (GPL) biosynthesis (support from host acetyl-CoA) and pentose phosphate pathway (PPP), while its incomplete TCA and fatty acid biosynthesis are less active. The modeled metabolic pathways are much more active in RB than in EB. Our in silico model suggests that EB and RB utilize folate to generate NAD(P)H using independent pathways. The only low metabolic flux inferred for EB involves mainly carbohydrate metabolism. RB utilizes energy -rich compounds to generate ATP in nucleic acid metabolism. Validation data for the modeling include proteomics experiments (model basis) as well as qRT-PCR confirmation of selected metabolic enzyme mRNA expression differences. The metabolic modeling is made fully available here. Its detailed insights and models on Ct metabolic adaptations during infection are a useful modeling basis for future studies. [ABSTRACT FROM AUTHOR]

The abundance of methane in shale gas and of other gases such as carbon monoxide, hydrogen, and carbon dioxide as chemical process byproducts has motivated the use of gas fermentation for bioproduction. Recent advances in metabolic engineering and synthetic biology allow for engineering of microbes metabolizing a variety of chemicals including gaseous feeds into a number of biorenewables and transportation liquid fuels. New computational tools enable the systematic exploration of all feasible conversion alternatives. Here we computationally assessed all thermodynamically feasible ways of co-utilizing CH4, CO, and CO2 using ferric as terminal electron acceptor for the production of all key precursor metabolites. We identified the thermodynamically feasible co-utilization ratio ranges of CH4, CO, and CO2 toward production of the target metabolite(s) as a function of ferric uptake. A revised version of the iMAC868 genome-scale metabolic model of Methanosarcina acetivorans was chosen to assess co-utilization of CH4, CO, and CO2 and their conversion into selected target products using the optStoic pathway design tool. This revised version contains the latest information on electron flow mechanisms by the methanogen while supplied with methane as the sole carbon source. The interplay between different gas co-utilization ratios and the energetics of reverse methanogenesis were also analyzed using the same metabolic model.

Bacterial endosymbionts and their insect hosts establish an intimate metabolic relationship. Bacteria offer a variety of essential nutrients to their hosts, whereas insect cells provide the necessary sources of matter and energy to their tiny metabolic allies. These nutritional complementations sustain themselves on a diversity of metabolite exchanges between the cell host and the reduced yet highly specialized bacterial metabolism—which, for instance, overproduces a small set of essential amino acids and vitamins. A well-known case of metabolic complementation is provided by the cedar aphid Cinara cedri that harbors two co-primary endosymbionts, Buchnera aphidicola BCc and Ca. Serratia symbiotica SCc, and in which some metabolic pathways are partitioned between different partners. Here we present a genome-scale metabolic network (GEM) for the bacterial consortium from the cedar aphid iBSCc. The analysis of this GEM allows us the confirmation of cases of metabolic complementation previously described by genome analysis (i.e., tryptophan and biotin biosynthesis) and the redefinition of an event of metabolic pathway sharing between the two endosymbionts, namely the biosynthesis of tetrahydrofolate. In silico knock-out experiments with iBSCc showed that the consortium metabolism is a highly integrated yet fragile network. We also have explored the evolutionary pathways leading to the emergence of metabolic complementation between reduced metabolisms starting from individual, complete networks. Our results suggest that, during the establishment of metabolic complementation in endosymbionts, adaptive evolution is significant in the case of tryptophan biosynthesis, whereas vitamin production pathways seem to adopt suboptimal solutions.

The abundance of methane in shale gas and of other gases such as carbon monoxide, hydrogen, and carbon dioxide as chemical process byproducts has motivated the use of gas fermentation for bioproduction. Recent advances in metabolic engineering and synthetic biology allow for engineering of microbes metabolizing a variety of chemicals including gaseous feeds into a number of biorenewables and transportation liquid fuels. New computational tools enable the systematic exploration of all feasible conversion alternatives. Here we computationally assessed all thermodynamically feasible ways of co-utilizing CH4, CO, and CO2 using ferric as terminal electron acceptor for the production of all key precursor metabolites. We identified the thermodynamically feasible co-utilization ratio ranges of CH4, CO, and CO2 toward production of the target metabolite(s) as a function of ferric uptake. A revised version of the iMAC868 genome-scale metabolic model of Methanosarcina acetivorans was chosen to assess co-utilization of CH4, CO, and CO2 and their conversion into selected target products using the optStoic pathway design tool. This revised version contains the latest information on electron flow mechanisms by the methanogen while supplied with methane as the sole carbon source. The interplay between different gas co-utilization ratios and the energetics of reverse methanogenesis were also analyzed using the same metabolic model. [ABSTRACT FROM AUTHOR]

Bacterial endosymbionts and their insect hosts establish an intimate metabolic relationship. Bacteria offer a variety of essential nutrients to their hosts, whereas insect cells provide the necessary sources of matter and energy to their tiny metabolic allies. These nutritional complementations sustain themselves on a diversity of metabolite exchanges between the cell host and the reduced yet highly specialized bacterial metabolism--which, for instance, overproduces a small set of essential amino acids and vitamins. A well-known case of metabolic complementation is provided by the cedar aphid Cinara cedri that harbors two co-primary endosymbionts, Buchnera aphidicola BCc and Ca. Serratia symbiotica SCc, and in which some metabolic pathways are partitioned between different partners. Here we present a genome-scale metabolic network (GEM) for the bacterial consortium from the cedar aphid iBSCc. The analysis of this GEM allows us the confirmation of cases of metabolic complementation previously described by genome analysis (i.e., tryptophan and biotin biosynthesis) and the redefinition of an event of metabolic pathway sharing between the two endosymbionts, namely the biosynthesis of tetrahydrofolate. In silico knock-out experiments with iBSCc showed that the consortium metabolism is a highly integrated yet fragile network. We also have explored the evolutionary pathways leading to the emergence of metabolic complementation between reduced metabolisms starting from individual, complete networks. Our results suggest that, during the establishment of metabolic complementation in endosymbionts, adaptive evolution is significant in the case of tryptophan biosynthesis, whereas vitamin production pathways seem to adopt suboptimal solutions. [ABSTRACT FROM AUTHOR]

In this study we have applied an integrated system biology approach to characterize the metabolic landscape of Streptomyces ambofaciens and to identify a list of potential metabolic engineering targets for the overproduction of the secondary metabolites in this microorganism. We focused on an often overlooked growth period (i.e., post-first rapid growth phase) and, by integrating constraint-based metabolic modeling with time resolved RNA-seq data, we depicted the main effects of changes in gene expression on the overall metabolic reprogramming occurring in S. ambofaciens. Moreover, through metabolic modeling, we unraveled a set of candidate overexpression gene targets hypothetically leading to spiramycin overproduction. Model predictions were experimentally validated by genetic manipulation of the recently described ethylmalonyl-CoA metabolic node, providing evidence that spiramycin productivity may be increased by enhancing the carbon flow through this pathway. The goal was achieved by over-expressing the ccr paralog srm4 in an ad hoc engineered plasmid. This work embeds the first metabolic reconstruction of S. ambofaciens and the successful experimental validation of model predictions and demonstrates the validity and the importance of in silico modeling tools for the overproduction of molecules with a biotechnological interest. Finally, the proposed metabolic reconstruction, which includes manually refined pathways for several secondary metabolites with antimicrobial activity, represents a solid platform for the future exploitation of S. ambofaciens biotechnological potential.

In this study we have applied an integrated system biology approach to characterize the metabolic landscape of Streptomyces ambofaciens and to identify a list of potential metabolic engineering targets for the overproduction of the secondary metabolites in this microorganism. We focused on an often overlooked growth period (i.e., post-first rapid growth phase) and by integrating constraint-based metabolic modeling with time resolved RNA-seq data, we depicted the main effects of changes in gene expression on the overall metabolic reprogramming occurring in S. Ambofaciens. Moreover, through metabolic modeling, we unraveled a set of candidate overexpression gene targets hypothetically leading to spiramycin overproduction. Model predictions were experimentally validated by genetic manipulation of the recently described ethylmalonyl-CoA metabolic node, providing evidence that spiramycin productivity may be increased by enhancing the carbon flow through this pathway. The goal was achieved by overexpressing the ccr paralog srm4 in an ad hoc engineered plasmid. This work embeds the first metabolic reconstruction of S. ambofaciens and the successful experimental validation of model predictions and demonstrates the validity and the importance of in silico modeling tools for the overproduction of molecules with a biotechnological interest. Finally, the proposed metabolic reconstruction, which includes manually refined pathways for several secondary metabolites with antimicrobial activity, represents a solid platform for the future exploitation of S. Ambofaciens biotechnological potential. [ABSTRACT FROM AUTHOR]

Facilitative interactions between microbial species are ubiquitous in various types of ecosystems on the Earth. Therefore, inferring how entangled webs of interspecific interactions shift through time in microbial ecosystems is an essential step for understanding ecological processes driving microbiome dynamics. By compiling shotgun metagenomic sequencing data of an experimental microbial community, we examined how the architectural features of facilitative interaction networks could change through time. A metabolic modeling approach for estimating dependence between microbial genomes (species) allowed us to infer the network structure of potential facilitative interactions at 13 time points through the 110-day monitoring of experimental microbiomes. We then found that positive feedback loops, which were theoretically predicted to promote cascade breakdown of ecological communities, existed within the inferred networks of metabolic interactions prior to the drastic community-compositional shift observed in the microbiome time-series. We further applied "directed-graph" analyses to pinpoint potential keystone species located at the "upper stream" positions of such feedback loops. These analyses on facilitative interactions will help us understand key mechanisms causing catastrophic shifts in microbial community structure., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Fujita, Ushio, Suzuki, Abe, Yamamichi, Okazaki, Canarini, Hayashi, Fukushima, Fukuda, Kiers and Toju.)

Bacterial communities may display metabolic complementation, in which different members of the association partially contribute to the same biosynthetic pathway. In this way, the end product of the pathway is synthesized by the community as a whole. However, the emergence and the benefits of such complementation are poorly understood. Herein, we present a simple model to analyze the metabolic interactions among bacteria, including the host in the case of endosymbiotic bacteria. The model considers two cell populations, with both cell types encoding for the same linear biosynthetic pathway. We have found that, for metabolic complementation to emerge as an optimal strategy, both product inhibition and large permeabilities are needed. In the light of these results, we then consider the patterns found in the case of tryptophan biosynthesis in the endosymbiont consortium hosted by the aphid Cinara cedri. Using in-silico computed physicochemical properties of metabolites of this and other biosynthetic pathways, we verified that the splitting point of the pathway corresponds to the most permeable intermediate. [ABSTRACT FROM AUTHOR]

The ability to build accurate protein families is a fundamental operation in bioinformatics that influences comparative analyses, genome annotation and metabolic modeling. For several years we have been maintaining protein families for all microbial genomes in the PATRIC database (Pathosystems Resource Integration Center, patricbrc.org) in order to drive many of the comparative analysis tools that are available through the PATRIC website. However, due to the burgeoning number of genomes, traditional approaches for generating protein families are becoming prohibitive. In this report, we describe a new approach for generating protein families, which we call PATtyFams. This method uses the k-mer-based function assignments available through RAST (Rapid Annotation using Subsystem Technology) to rapidly guide family formation, and then differentiates the function-based groups into families using a Markov Cluster algorithm (MCL). This new approach for generating protein families is rapid, scalable and has properties that are consistent with alignment-based methods.

The bacterial microbiota of plants is diverse, with 1,000s of operational taxonomic units (OTUs) associated with any individual plant. In this work we investigate the differences between 19 sequenced Pseudomonas fluorescens strains, isolated from Populus deltoides rhizosphere and endosphere and which represent a single OTU, using phenotypic analysis, comparative genomics, and metabolic models. While no traits were exclusive to either endosphere or rhizosphere P. fluorescens isolates, multiple pathways relevant for plant-bacterial interactions are enriched in endosphere isolate genomes. Further, growth phenotypes such as phosphate solubilization, protease activity, denitrification and root growth promotion are biased towards endosphere isolates. Endosphere isolates have significantly more metabolic pathways for plant signaling compounds and an increased metabolic range that includes utilization of energy rich nucleotides and sugars, consistent with endosphere colonization. Rhizosphere P. fluorescens have fewer pathways representative of plant-bacterial interactions but show metabolic bias towards chemical substrates often found in root exudates. This work reveals the diverse functions that may contribute to colonization of the endosphere by bacteria and are enriched among closely related isolates.

Conflict and cooperation between bacterial species drive the composition and function of microbial communities. Stability of these emergent properties will be influenced by the degree to which species’ interactions are robust to genetic perturbations. We use genome-scale metabolic modeling to computationally analyze the impact of genetic changes when Escherichia coli and Salmonella enterica compete, or cooperate. We systematically knocked out in silico each reaction in the metabolic network of E. coli to construct all 2,583 mutant stoichiometric models. Then, using a recently developed multi-scale computational framework, we simulated the growth of each mutant E. coli in the presence of S. enterica. The type of interaction between species was set by modulating the initial metabolites present in the environment. We found that the community was most robust to genetic perturbations when the organisms were cooperating. Species ratios were more stable in the cooperative community, and community biomass had equal variance in the two contexts. Additionally, the number of mutations that have a substantial effect is lower when the species cooperate than when they are competing. In contrast, when mutations were added to the S. enterica network the system was more robust when the bacteria were competing. These results highlight the utility of connecting metabolic mechanisms and studies of ecological stability. Cooperation and conflict alter the connection between genetic changes and properties that emerge at higher levels of biological organization.

Numerous methods for classifying gene activity states based on gene expression data have been proposed for use in downstream applications, such as incorporating transcriptomics data into metabolic models in order to improve resulting flux predictions. These methods often attempt to classify gene activity for each gene in each experimental condition as belonging to one of two states: active (the gene product is part of an active cellular mechanism) or inactive (the cellular mechanism is not active). These existing methods of classifying gene activity states suffer from multiple limitations, including enforcing unrealistic constraints on the overall proportions of active and inactive genes, failing to leverage a priori knowledge of gene co-regulation, failing to account for differences between genes, and failing to provide statistically meaningful confidence estimates. We propose a flexible Bayesian approach to classifying gene activity states based on a Gaussian mixture model. The model integrates genome-wide transcriptomics data from multiple conditions and information about gene co-regulation to provide activity state confidence estimates for each gene in each condition. We compare the performance of our novel method to existing methods on both simulated data and real data from 907 E. coli gene expression arrays, as well as a comparison with experimentally measured flux values in 29 conditions, demonstrating that our method provides more consistent and accurate results than existing methods across a variety of metrics. [ABSTRACT FROM AUTHOR]

Microbiological studies are increasingly relying on in silico methods to perform exploration and rapid analysis of genomic data, and functional genomics studies are supplemented by the new perspectives that genome-scale metabolic models offer. A mathematical model consisting of a microbe's entire metabolic map can be rapidly determined from whole-genome sequencing and annotating the genomic material encoded in its DNA. Flux-balance analysis (FBA), a linear programming technique that uses metabolic models to predict the phenotypic responses imposed by environmental elements and factors, is the leading method to simulate and manipulate cellular growth in silico. However, the process of creating an accurate model to use in FBA consists of a series of steps involving a multitude of connections between bioinformatics databases, enzyme resources, and metabolic pathways. We present the methodology and procedure to obtain a metabolic model using PyFBA, an extensible Pythonbased open-source software package aimed to provide a platform where functional annotations are used to build metabolic models (http://linsalrob.github.io/PyFBA). Backed by the Model SEED biochemistry database, PyFBA contains methods to reconstruct a microbe's metabolic map, run FBA upon different media conditions, and gap-fill its metabolism. The extensibility of PyFBA facilitates novel techniques in creating accurate genome-scale metabolic models. [ABSTRACT FROM AUTHOR]

The ability to build accurate protein families is a fundamental operation in bioinformatics that influences comparative analyses, genome annotation, and metabolic modeling. For several years we have been maintaining protein families for all microbial genomes in the PATRIC database (Pathosystems Resource Integration Center, patricbrc.org) in order to drive many of the comparative analysis tools that are available through the PATRIC website. However, due to the burgeoning number of genomes, traditional approaches for generating protein families are becoming prohibitive. In this report, we describe a new approach for generating protein families, which we call PATtyFams. This method uses the k-mer-based function assignments available through RAST (Rapid Annotation using Subsystem Technology) to rapidly guide family formation, and then differentiates the function-based groups into families using a Markov Cluster algorithm (MCL). This new approach for generating protein families is rapid, scalable and has properties that are consistent with alignment-based methods. [ABSTRACT FROM AUTHOR]

Lignocellulosic hydrolysate (LCH) inhibitors are a large class of bioactive molecules that arise from pretreatment, hydrolysis, and fermentation of plant biomass. These diverse compounds reduce lignocellulosic biofuel yields by inhibiting cellular processes and diverting energy into cellular responses. LCH inhibitors present one of the most significant challenges to efficient biofuel production by microbes. Development of new strains that lessen the effects of LCH inhibitors is an economically favorable strategy relative to expensive detoxification methods that also can reduce sugar content in deconstructed biomass. Systems biology analyses and metabolic modeling combined with directed evolution and synthetic biology are successful strategies for biocatalyst development, and methods that leverage state-of-the-art tools are needed to overcome inhibitors more completely. This perspective considers the energetic costs of LCH inhibitors and technologies that can be used to overcome their drain on conversion efficiency. We suggest academic and commercial research groups could benefit by sharing data on LCH inhibitors and implementing translational biofuel research.

Conflict and cooperation between bacterial species drive the composition and function of microbial communities. Stability of these emergent properties will be influenced by the degree to which species' interactions are robust to genetic perturbations. We use genome-scale metabolic modeling to computationally analyze the impact of genetic changes when Escherichia coli and Salmonella enterica compete, or cooperate. We systematically knocked out in silico each reaction in the metabolic network of E. coli to construct all 2583 mutant stoichiometric models. Then, using a recently developed multi-scale computational framework, we simulated the growth of each mutant E. coli in the presence of S. enterica. The type of interaction between species was set by modulating the initial metabolites present in the environment. We found that the community was most robust to genetic perturbations when the organisms were cooperating. Species ratios were more stable in the cooperative community, and community biomass had equal variance in the two contexts. Additionally, the number of mutations that have a substantial effect is lower when the species cooperate than when they are competing. In contrast, when mutations were added to the S. enterica network the system was more robust when the bacteria were competing. These results highlight the utility of connecting metabolic mechanisms and studies of ecological stability. Cooperation and conflict alter the connection between genetic changes and properties that emerge at higher levels of biological organization. [ABSTRACT FROM AUTHOR]

Bacterial communities may display metabolic complementation, in which different members of the association partially contribute to the same biosynthetic pathway. In this way, the end product of the pathway is synthesized by the community as a whole. However, the emergence and the benefits of such complementation are poorly understood. Herein, we present a simple model to analyze the metabolic interactions among bacteria, including the host in the case of endosymbiotic bacteria. The model considers two cell populations, with both cell types encoding for the same linear biosynthetic pathway. We have found that, for metabolic complementation to emerge as an optimal strategy, both product inhibition and large permeabilities are needed. In the light of these results, we then consider the patterns found in the case of tryptophan biosynthesis in the endosymbiont consortium hosted by the aphid Cinara cedri . Using in-silico computed physicochemical properties of metabolites of this and other biosynthetic pathways, we verified that the splitting point of the pathway corresponds to the most permeable intermediate.

Lignocellulosic hydrolysate (LCH) inhibitors are a large class of bioactive molecules that arise from pretreatment, hydrolysis, and fermentation of plant biomass. These diverse compounds reduce lignocellulosic biofuel yields by inhibiting cellular processes and diverting energy into cellular responses. LCH inhibitors present one of the most significant challenges to efficient biofuel production by microbes. Development of new strains that lessen the effects of LCH inhibitors is an economically favorable strategy relative to expensive detoxification methods that also can reduce sugar content in deconstructed biomass. Systems biology analyses and metabolic modeling combined with directed evolution and synthetic biology are successful strategies for biocatalyst development, and methods that leverage state-of-the-art tools are needed to overcome inhibitors more completely. This perspective considers the energetic costs of LCH inhibitors and technologies that can be used to overcome their drain on conversion efficiency. We suggest academic and commercial research groups could benefit by sharing data on LCH inhibitors and implementing "translational biofuel research."