Genome-scale Metabolic Reconstructions Research Articles

Spatial transcriptome-guided multi-scale framework connects P. aeruginosa metabolic states to oxidative stress biofilm microenvironment.

With the generation of spatially resolved transcriptomics of microbial biofilms, computational tools can be used to integrate this data to elucidate the multi-scale mechanisms controlling heterogeneous biofilm metabolism. This work presents a Multi-scale model of Metabolism In Cellular Systems (MiMICS) which is a computational framework that couples a genome-scale metabolic network reconstruction (GENRE) with Hybrid Automata Library (HAL), an existing agent-based model and reaction-diffusion model platform. A key feature of MiMICS is the ability to incorporate multiple -omics-guided metabolic models, which can represent unique metabolic states that yield different metabolic parameter values passed to the extracellular models. We used MiMICS to simulate Pseudomonas aeruginosa regulation of denitrification and oxidative stress metabolism in hypoxic and nitric oxide (NO) biofilm microenvironments. Integration of P. aeruginosa PA14 biofilm spatial transcriptomic data into a P. aeruginosa PA14 GENRE generated four PA14 metabolic model states that were input into MiMICS. Characteristic of aerobic, denitrification, and oxidative stress metabolism, the four metabolic model states predicted different oxygen, nitrate, and NO exchange fluxes that were passed as inputs to update the agent's local metabolite concentrations in the extracellular reaction-diffusion model. Individual bacterial agents chose a PA14 metabolic model state based on a combination of stochastic rules, and agents sensing local oxygen and NO. Transcriptome-guided MiMICS predictions suggested microscale denitrification and oxidative stress metabolic heterogeneity emerged due to local variability in the NO biofilm microenvironment. MiMICS accurately predicted the biofilm's spatial relationships between denitrification, oxidative stress, and central carbon metabolism. As simulated cells responded to extracellular NO, MiMICS revealed dynamics of cell populations heterogeneously upregulating reactions in the denitrification pathway, which may function to maintain NO levels within non-toxic ranges. We demonstrated that MiMICS is a valuable computational tool to incorporate multiple -omics-guided metabolic models to mechanistically map heterogeneous microbial metabolic states to the biofilm microenvironment.

Metabolic modeling of Halomonas campaniensis improves polyhydroxybutyrate production under nitrogen limitation

Poly-hydroxybutyrate (PHB) is an environmentally friendly alternative for conventional fossil fuel-based plastics that is produced by various microorganisms. Large-scale PHB production is challenging due to the comparatively higher biomanufacturing costs. A PHB overproducer is the haloalkaliphilic bacterium Halomonas campaniensis, which has low nutritional requirements and can grow in cultures with high salt concentrations, rendering it resistant to contamination. Despite its virtues, the metabolic capabilities of H. campaniensis as well as the limitations hindering higher PHB production remain poorly studied. To address this limitation, we present HaloGEM, the first high-quality genome-scale metabolic network reconstruction, which encompasses 888 genes, 1528 reactions (1257 gene-associated), and 1274 metabolites. HaloGEM not only displays excellent agreement with previous growth data and experiments from this study, but it also revealed nitrogen as a limiting nutrient when growing aerobically under high salt concentrations using glucose as carbon source. Among different nitrogen source mixtures for optimal growth, HaloGEM predicted glutamate and arginine as a promising mixture producing increases of 54.2% and 153.4% in the biomass yield and PHB titer, respectively. Furthermore, the model was used to predict genetic interventions for increasing PHB yield, which were consistent with the rationale of previously reported strategies. Overall, the presented reconstruction advances our understanding of the metabolic capabilities of H. campaniensis for rationally engineering this next-generation industrial biotechnology platform.Key pointsA comprehensive genome-scale metabolic reconstruction of H. campaniensis was developed.Experiments and simulations predict N limitation in minimal media under aerobiosis.In silico media design increased experimental biomass yield and PHB titer.

Identifying metabolic adaptations characteristic of cardiotoxicity using paired transcriptomics and metabolomics data integrated with a computational model of heart metabolism.

Improvements in the diagnosis and treatment of cancer have revealed long-term side effects of chemotherapeutics, particularly cardiotoxicity. Here, we present paired transcriptomics and metabolomics data characterizing in vitro cardiotoxicity to three compounds: 5-fluorouracil, acetaminophen, and doxorubicin. Standard gene enrichment and metabolomics approaches identify some commonly affected pathways and metabolites but are not able to readily identify metabolic adaptations in response to cardiotoxicity. The paired data was integrated with a genome-scale metabolic network reconstruction of the heart to identify shifted metabolic functions, unique metabolic reactions, and changes in flux in metabolic reactions in response to these compounds. Using this approach, we confirm previously seen changes in the p53 pathway by doxorubicin and RNA synthesis by 5-fluorouracil, we find evidence for an increase in phospholipid metabolism in response to acetaminophen, and we see a shift in central carbon metabolism suggesting an increase in metabolic demand after treatment with doxorubicin and 5-fluorouracil.

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.

The limitations of phenotype prediction in metabolism.

Phenotype prediction is at the center of many questions in biology. Prediction is often achieved by determining statistical associations between genetic and phenotypic variation, ignoring the exact processes that cause the phenotype. Here, we present a framework based on genome-scale metabolic reconstructions to reveal the mechanisms behind the associations. We calculated a polygenic score (PGS) that identifies a set of enzymes as predictors of growth, the phenotype. This set arises from the synergy of the functional mode of metabolism in a particular setting and its evolutionary history, and is suitable to infer the phenotype across a variety of conditions. We also find that there is optimal genetic variation for predictability and demonstrate how the linear PGS can still explain phenotypes generated by the underlying nonlinear biochemistry. Therefore, the explicit model interprets the black box statistical associations of the genotype-to-phenotype map and helps to discover what limits the prediction in metabolism.

An integrated systems biology approach reveals differences in formate metabolism in the genus Methanothermobacter

Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.

Unveiling metabolic pathways involved in the extreme desiccation tolerance of an Atacama cyanobacterium

Gloeocapsopsis dulcis strain AAB1 is an extremely xerotolerant cyanobacterium isolated from the Atacama Desert (i.e., the driest and oldest desert on Earth) that holds astrobiological significance due to its ability to biosynthesize compatible solutes at ultra-low water activities. We sequenced and assembled the G. dulcis genome de novo using a combination of long- and short-read sequencing, which resulted in high-quality consensus sequences of the chromosome and two plasmids. We leveraged the G. dulcis genome to generate a genome-scale metabolic model (iGd895) to simulate growth in silico. iGd895 represents, to our knowledge, the first genome-scale metabolic reconstruction developed for an extremely xerotolerant cyanobacterium. The model's predictive capability was assessed by comparing the in silico growth rate with in vitro growth rates of G. dulcis, in addition to the synthesis of trehalose. iGd895 allowed us to explore simulations of key metabolic processes such as essential pathways for water-stress tolerance, and significant alterations to reaction flux distribution and metabolic network reorganization resulting from water limitation. Our study provides insights into the potential metabolic strategies employed by G. dulcis, emphasizing the crucial roles of compatible solutes, metabolic water, energy conservation, and the precise regulation of reaction rates in their adaptation to water stress.

A genome-scale metabolic reconstruction of soybean and Bradyrhizobium diazoefficiens reveals the cost-benefit of nitrogen fixation.

Nitrogen-fixing symbioses allow legumes to thrive in nitrogen-poor soils at the cost of diverting some photoassimilate to their microsymbionts. Effort is being made to bioengineer nitrogen fixation into nonleguminous crops. This requires a quantitative understanding of its energetic costs and the links between metabolic variations and symbiotic efficiency. A whole-plant metabolic model for soybean (Glycine max) with its associated microsymbiont Bradyrhizobium diazoefficiens was developed and applied to predict the cost-benefit of nitrogen fixation with varying soil nitrogen availability. The model predicted a nitrogen-fixation cost of c. 4.13 g C g-1 N, which when implemented into a crop scale model, translated to a grain yield reduction of 27% compared with a non-nodulating plant receiving its nitrogen from the soil. Considering the lower nitrogen content of cereals, the yield cost to a hypothetical N-fixing cereal is predicted to be less than half that of soybean. Soybean growth was predicted to be c. 5% greater when the nodule nitrogen export products were amides versus ureides. This is the first metabolic reconstruction in a tropical crop species that simulates the entire plant and nodule metabolism. Going forward, this model will serve as a tool to investigate carbon use efficiency and key mechanisms within N-fixing symbiosis in a tropical species forming determinate nodules.

COVRECON: automated integration of genome- and metabolome-scale network reconstruction and data-driven inverse modeling of metabolic interaction networks.

One central goal of systems biology is to infer biochemical regulations from large-scale OMICS data. Many aspects of cellular physiology and organismal phenotypes can be understood as results of metabolic interaction network dynamics. Previously, we have proposed a convenient mathematical method, which addresses this problem using metabolomics data for the inverse calculation of biochemical Jacobian matrices revealing regulatory checkpoints of biochemical regulations. The proposed algorithms for this inference are limited by two issues: they rely on structural network information that needs to be assembled manually, and they are numerically unstable due to ill-conditioned regression problems for large-scale metabolic networks. To address these problems, we developed a novel regression loss-based inverse Jacobian algorithm, combining metabolomics COVariance and genome-scale metabolic RECONstruction, which allows for a fully automated, algorithmic implementation of the COVRECON workflow. It consists of two parts: (i) Sim-Network and (ii) inverse differential Jacobian evaluation. Sim-Network automatically generates an organism-specific enzyme and reaction dataset from Bigg and KEGG databases, which is then used to reconstruct the Jacobian's structure for a specific metabolomics dataset. Instead of directly solving a regression problem as in the previous workflow, the new inverse differential Jacobian is based on a substantially more robust approach and rates the biochemical interactions according to their relevance from large-scale metabolomics data. The approach is illustrated by in silico stochastic analysis with differently sized metabolic networks from the BioModels database and applied to a real-world example. The characteristics of the COVRECON implementation are that (i) it automatically reconstructs a data-driven superpathway model; (ii) more general network structures can be investigated, and (iii) the new inverse algorithm improves stability, decreases computation time, and extends to large-scale models. The code is available in the website https://bitbucket.org/mosys-univie/covrecon.

Transcriptome-guided metabolic network analysis reveals rearrangements of carbon flux distribution in Neisseria gonorrhoeae during neutrophil co-culture.

The ability of bacterial pathogens to metabolically adapt to the environmental conditions of their hosts is critical to both colonization and invasive disease. Infection with Neisseria gonorrhoeae (the gonococcus, Gc) is characterized by the influx of neutrophils [polymorphonuclear leukocytes (PMNs)], which fail to clear the bacteria and make antimicrobial products that can exacerbate tissue damage. The inability of the human host to clear Gc infection is particularly concerning in light of the emergence of strains that are resistant to all clinically recommended antibiotics. Bacterial metabolism represents a promising target for the development of new therapeutics against Gc. Here, we generated a curated genome-scale metabolic network reconstruction (GENRE) of Gc strain FA1090. This GENRE links genetic information to metabolic phenotypes and predicts Gc biomass synthesis and energy consumption. We validated this model with published data and in new results reported here. Contextualization of this model using the transcriptional profile of Gc exposed to PMNs revealed substantial rearrangements of Gc central metabolism and induction of Gc nutrient acquisition strategies for alternate carbon source use. These features enhanced the growth of Gc in the presence of neutrophils. From these results, we conclude that the metabolic interplay between Gc and PMNs helps define infection outcomes. The use of transcriptional profiling and metabolic modeling to reveal new mechanisms by which Gc persists in the presence of PMNs uncovers unique aspects of metabolism in this fastidious bacterium, which could be targeted to block infection and thereby reduce the burden of gonorrhea in the human population. IMPORTANCE The World Health Organization designated Gc as a high-priority pathogen for research and development of new antimicrobials. Bacterial metabolism is a promising target for new antimicrobials, as metabolic enzymes are widely conserved among bacterial strains and are critical for nutrient acquisition and survival within the human host. Here we used genome-scale metabolic modeling to characterize the core metabolic pathways of this fastidious bacterium and to uncover the pathways used by Gc during culture with primary human immune cells. These analyses revealed that Gc relies on different metabolic pathways during co-culture with human neutrophils than in rich media. Conditionally essential genes emerging from these analyses were validated experimentally. These results show that metabolic adaptation in the context of innate immunity is important to Gc pathogenesis. Identifying the metabolic pathways used by Gc during infection can highlight new therapeutic targets for drug-resistant gonorrhea.

Gimibacter soli gen. nov. sp. nov., isolated from mangrove soil and insight into its ecological distribution and metabolic potential.

A Gram-stain-negative, facultatively anaerobic, motile and rod-shaped bacterium, designated as 6D33T, was isolated from mangrove soil. Growth was found to occur at 15-32 °C (optimum, 28 °C), at pH 6-9 (optimum, pH 7) and in 0-3 % NaCl (optimum, 1 %, w/v). The results of 16S rRNA gene-based analysis showed that strain 6D33T belonged to the family Temperatibacteraceae, sharing 93.1-94.4 % identity with its close neighbours within the genus Kordiimonas. The phylogenomic results indicated that strain 6D33T formed an independent branch distinct from type strains of the genus Kordiimonas. The overall genome relatedness indices of digital DNA-DNA hybridization, average nucleotide identity and amino acid identity values showed that strain 6D33T represents a novel species of a novel genus. The results of chemotaxonomic characterization indicated that the major cellular fatty acids of strain 6D33T were summed feature 9 (C16 : 0 10-methyl and/or iso-C17 : 1 ω9c), summed feature 3 (C16 : 1 ω6c and/or C16 : 1 ω7c) and iso-C15 : 0; the polar lipids comprised diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, an unidentified aminolipid and three unidentified lipids; the only respiratory quinone was ubiquinone-10. The genomic size and DNA G+C contents were 3.59 Mbp and 60.84 mol%, respectively. The 16S rRNA gene sequence reads abundance profiles revealed that the rare taxon is prevalent in marine environments, especially in sediments. Genome-scale metabolic reconstruction of strain 6D33T revealed a heterotrophic lifestyle and many pathways responsible for the degradation of aromatic compounds, suggesting application potential in aromatic hydrocarbon removal. Based on its genotypic and phenotypic characteristics, strain 6D33T is concluded to represent a novel species of the novel genus in the family Temperatibacteraceae, for which the name Gimibacter soli gen. nov. sp. nov. is proposed. The type strain of the type species is 6D33T (=GDMCC 1.1959T=KCTC 82335T).

Genome-Scale Metabolic Reconstruction, Non-Targeted LC-QTOF-MS Based Metabolomics Data, and Evaluation of Anticancer Activity of Cannabis sativa Leaf Extracts.

Over the past decades, Colombia has suffered complex social problems related to illicit crops, including forced displacement, violence, and environmental damage, among other consequences for vulnerable populations. Considerable effort has been made in the regulation of illicit crops, predominantly Cannabis sativa, leading to advances such as the legalization of medical cannabis and its derivatives, the improvement of crops, and leaving an open window to the development of scientific knowledge to explore alternative uses. It is estimated that C. sativa can produce approximately 750 specialized secondary metabolites. Some of the most relevant due to their anticancer properties, besides cannabinoids, are monoterpenes, sesquiterpenoids, triterpenoids, essential oils, flavonoids, and phenolic compounds. However, despite the increase in scientific research on the subject, it is necessary to study the primary and secondary metabolism of the plant and to identify key pathways that explore its great metabolic potential. For this purpose, a genome-scale metabolic reconstruction of C. sativa is described and contextualized using LC-QTOF-MS metabolic data obtained from the leaf extract from plants grown in the region of Pesca-Boyaca, Colombia under greenhouse conditions at the Clever Leaves facility. A compartmentalized model with 2101 reactions and 1314 metabolites highlights pathways associated with fatty acid biosynthesis, steroids, and amino acids, along with the metabolism of purine, pyrimidine, glucose, starch, and sucrose. Key metabolites were identified through metabolomic data, such as neurine, cannabisativine, cannflavin A, palmitoleic acid, cannabinoids, geranylhydroquinone, and steroids. They were analyzed and integrated into the reconstruction, and their potential applications are discussed. Cytotoxicity assays revealed high anticancer activity against gastric adenocarcinoma (AGS), melanoma cells (A375), and lung carcinoma cells (A549), combined with negligible impact against healthy human skin cells.

Optimizing Strategies for Bio-Based Ethanol Production Using Genome-Scale Metabolic Modeling of the Hyperthermophilic Archaeon, Pyrococcus furiosus.

A genome-scale metabolic model, encompassing a total of 623 genes, 727 reactions, and 865 metabolites, was developed for Pyrococcus furiosus, an archaeon that grows optimally at 100°C by carbohydrate and peptide fermentation. The model uses subsystem-based genome annotation, along with extensive manual curation of 237 gene-reaction associations including those involved in central carbon metabolism, amino acid metabolism, and energy metabolism. The redox and energy balance of P. furiosus was investigated through random sampling of flux distributions in the model during growth on disaccharides. The core energy balance of the model was shown to depend on high acetate production and the coupling of a sodium-dependent ATP synthase and membrane-bound hydrogenase, which generates a sodium gradient in a ferredoxin-dependent manner, aligning with existing understanding of P. furiosus metabolism. The model was utilized to inform genetic engineering designs that favor the production of ethanol over acetate by implementing an NADPH and CO-dependent energy economy. The P. furiosus model is a powerful tool for understanding the relationship between generation of end products and redox/energy balance at a systems-level that will aid in the design of optimal engineering strategies for production of bio-based chemicals and fuels. IMPORTANCE The bio-based production of organic chemicals provides a sustainable alternative to fossil-based production in the face of today's climate challenges. In this work, we present a genome-scale metabolic reconstruction of Pyrococcus furiosus, a well-established platform organism that has been engineered to produce a variety of chemicals and fuels. The metabolic model was used to design optimal engineering strategies to produce ethanol. The redox and energy balance of P. furiosus was examined in detail, which provided useful insights that will guide future engineering designs.

Reconstructor: a COBRApy compatible tool for automated genome-scale metabolic network reconstruction with parsimonious flux-based gap-filling.

Genome-scale metabolic network reconstructions (GENREs) are valuable for understanding cellular metabolism in silico. Several tools exist for automatic GENRE generation. However, these tools frequently (i) do not readily integrate with some of the widely-used suites of packaged methods available for network analysis, (ii) lack effective network curation tools, (iii) are not sufficiently user-friendly, and (iv) often produce low-quality draft reconstructions. Here, we present Reconstructor, a user-friendly, COBRApy-compatible tool that produces high-quality draft reconstructions with reaction and metabolite naming conventions that are consistent with the ModelSEED biochemistry database and includes a gap-filling technique based on the principles of parsimony. Reconstructor can generate SBML GENREs from three input types: annotated protein .fasta sequences (Type 1 input), a BLASTp output (Type 2), or an existing SBML GENRE that can be further gap-filled (Type 3). While Reconstructor can be used to create GENREs of any species, we demonstrate the utility of Reconstructor with bacterial reconstructions. We demonstrate how Reconstructor readily generates high-quality GENRES that capture strain, species, and higher taxonomic differences in functional metabolism of bacteria and are useful for further biological discovery. The Reconstructor Python package is freely available for download. Complete installation and usage instructions and benchmarking data are available at http://github.com/emmamglass/reconstructor.

Systems biology elucidates the distinctive metabolic niche filled by the human gut microbe Eggerthella lenta.

Human gut bacteria perform diverse metabolic functions with consequences for host health. The prevalent and disease-linked Actinobacterium Eggerthella lenta performs several unusual chemical transformations, but it does not metabolize sugars and its core growth strategy remains unclear. To obtain a comprehensive view of the metabolic network of E. lenta, we generated several complementary resources: defined culture media, metabolomics profiles of strain isolates, and a curated genome-scale metabolic reconstruction. Stable isotope-resolved metabolomics revealed that E. lenta uses acetate as a key carbon source while catabolizing arginine to generate ATP, traits which could be recapitulated in silico by our updated metabolic model. We compared these in vitro findings with metabolite shifts observed in E. lenta-colonized gnotobiotic mice, identifying shared signatures across environments and highlighting catabolism of the host signaling metabolite agmatine as an alternative energy pathway. Together, our results elucidate a distinctive metabolic niche filled by E. lenta in the gut ecosystem. Our culture media formulations, atlas of metabolomics data, and genome-scale metabolic reconstructions form a freely available collection of resources to support further study of the biology of this prevalent gut bacterium.

Metabolic Robustness to Growth Temperature of a Cold- Adapted Marine Bacterium.

Microbial communities experience continuous environmental changes, with temperature fluctuations being the most impacting. This is particularly important considering the ongoing global warming but also in the "simpler" context of seasonal variability of sea-surface temperature. Understanding how microorganisms react at the cellular level can improve our understanding of their possible adaptations to a changing environment. In this work, we investigated the mechanisms through which metabolic homeostasis is maintained in a cold-adapted marine bacterium during growth at temperatures that differ widely (15 and 0°C). We have quantified its intracellular and extracellular central metabolomes together with changes occurring at the transcriptomic level in the same growth conditions. This information was then used to contextualize a genome-scale metabolic reconstruction, and to provide a systemic understanding of cellular adaptation to growth at 2 different temperatures. Our findings indicate a strong metabolic robustness at the level of the main central metabolites, counteracted by a relatively deep transcriptomic reprogramming that includes changes in gene expression of hundreds of metabolic genes. We interpret this as a transcriptomic buffering of cellular metabolism, able to produce overlapping metabolic phenotypes, despite the wide temperature gap. Moreover, we show that metabolic adaptation seems to be mostly played at the level of few key intermediates (e.g., phosphoenolpyruvate) and in the cross talk between the main central metabolic pathways. Overall, our findings reveal a complex interplay at gene expression level that contributes to the robustness/resilience of core metabolism, also promoting the leveraging of state-of-the-art multi-disciplinary approaches to fully comprehend molecular adaptations to environmental fluctuations. IMPORTANCE This manuscript addresses a central and broad interest topic in environmental microbiology, i.e. the effect of growth temperature on microbial cell physiology. We investigated if and how metabolic homeostasis is maintained in a cold-adapted bacterium during growth at temperatures that differ widely and that match measured changes on the field. Our integrative approach revealed an extraordinary robustness of the central metabolome to growth temperature. However, this was counteracted by deep changes at the transcriptional level, and especially in the metabolic part of the transcriptome. This conflictual scenario was interpreted as a transcriptomic buffering of cellular metabolism, and was investigated using genome-scale metabolic modeling. Overall, our findings reveal a complex interplay at gene expression level that contributes to the robustness/resilience of core metabolism, also promoting the use of state-of-the-art multi-disciplinary approaches to fully comprehend molecular adaptations to environmental fluctuations.

Computational Insight into Intraspecies Distinctions in Pseudoalteromonas distincta: Carotenoid-like Synthesis Traits and Genomic Heterogeneity

Advances in the computational annotation of genomes and the predictive potential of current metabolic models, based on more than thousands of experimental phenotypes, allow them to be applied to identify the diversity of metabolic pathways at the level of ecophysiology differentiation within taxa and to predict phenotypes, secondary metabolites, host-associated interactions, survivability, and biochemical productivity under proposed environmental conditions. The significantly distinctive phenotypes of members of the marine bacterial species Pseudoalteromonas distincta and an inability to use common molecular markers make their identification within the genus Pseudoalteromonas and prediction of their biotechnology potential impossible without genome-scale analysis and metabolic reconstruction. A new strain, KMM 6257, of a carotenoid-like phenotype, isolated from a deep-habituating starfish, emended the description of P. distincta, particularly in the temperature growth range from 4 to 37 °C. The taxonomic status of all available closely related species was elucidated by phylogenomics. P. distincta possesses putative methylerythritol phosphate pathway II and 4,4'-diapolycopenedioate biosynthesis, related to C30 carotenoids, and their functional analogues, aryl polyene biosynthetic gene clusters (BGC). However, the yellow-orange pigmentation phenotypes in some strains coincide with the presence of a hybrid BGC encoding for aryl polyene esterified with resorcinol. The alginate degradation and glycosylated immunosuppressant production, similar to brasilicardin, streptorubin, and nucleocidines, are the common predicted features. Starch, agar, carrageenan, xylose, lignin-derived compound degradation, polysaccharide, folate, and cobalamin biosynthesis are all strain-specific.

Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine.

The human microbiome influences the efficacy and safety of a wide variety of commonly prescribed drugs. Designing precision medicine approaches that incorporate microbial metabolism would require strain- and molecule-resolved, scalable computational modeling. Here, we extend our previous resource of genome-scale metabolic reconstructions of human gut microorganisms with a greatly expanded version. AGORA2 (assembly of gut organisms through reconstruction and analysis, version 2) accounts for 7,302 strains, includes strain-resolved drug degradation and biotransformation capabilities for 98 drugs, and was extensively curated based on comparative genomics and literature searches. The microbial reconstructions performed very well against three independently assembled experimental datasets with an accuracy of 0.72 to 0.84, surpassing other reconstruction resources and predicted known microbial drug transformations with an accuracy of 0.81. We demonstrate that AGORA2 enables personalized, strain-resolved modeling by predicting the drug conversion potential of the gut microbiomes from 616 patients with colorectal cancer and controls, which greatly varied between individuals and correlated with age, sex, body mass index and disease stages. AGORA2 serves as a knowledge base for the human microbiome and paves the way to personalized, predictive analysis of host-microbiome metabolic interactions.

Identifying metabolic shifts in Crohn's disease using 'omics-driven contextualized computational metabolic network models

Crohn's disease (CD) is a chronic inflammatory disease of the gastrointestinal tract. A clear gap in our existing CD diagnostics and current disease management approaches is the lack of highly specific biomarkers that can be used to streamline or personalize disease management. Comprehensive profiling of metabolites holds promise; however, these high-dimensional profiles need to be reduced to have relevance in the context of CD. Machine learning approaches are optimally suited to bridge this gap in knowledge by contextualizing the metabolic alterations in CD using genome-scale metabolic network reconstructions. Our work presents a framework for studying altered metabolic reactions between patients with CD and controls using publicly available transcriptomic data and existing gene-driven metabolic network reconstructions. Additionally, we apply the same methods to patient-derived ileal enteroids to explore the utility of using this experimental in vitro platform for studying CD. Furthermore, we have piloted an untargeted metabolomics approach as a proof-of-concept validation strategy in human ileal mucosal tissue. These findings suggest that in silico metabolic modeling can potentially identify pathways of clinical relevance in CD, paving the way for the future discovery of novel diagnostic biomarkers and therapeutic targets.

Genome-scale metabolic modeling reveals increased reliance on valine catabolism in clinical isolates of Klebsiella pneumoniae

Infections due to carbapenem-resistant Enterobacteriaceae have recently emerged as one of the most urgent threats to hospitalized patients within the United States and Europe. By far the most common etiological agent of these infections is Klebsiella pneumoniae, frequently manifesting in hospital-acquired pneumonia with a mortality rate of ~50% even with antimicrobial intervention. We performed transcriptomic analysis of data collected previously from in vitro characterization of both laboratory and clinical isolates which revealed shifts in expression of multiple master metabolic regulators across isolate types. Metabolism has been previously shown to be an effective target for antibacterial therapy, and genome-scale metabolic network reconstructions (GENREs) have provided a powerful means to accelerate identification of potential targets in silico. Combining these techniques with the transcriptome meta-analysis, we generated context-specific models of metabolism utilizing a well-curated GENRE of K. pneumoniae (iYL1228) to identify novel therapeutic targets. Functional metabolic analyses revealed that both composition and metabolic activity of clinical isolate-associated context-specific models significantly differs from laboratory isolate-associated models of the bacterium. Additionally, we identified increased catabolism of L-valine in clinical isolate-specific growth simulations. These findings warrant future studies for potential efficacy of valine transaminase inhibition as a target against K. pneumoniae infection.