The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes.

For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical 'Rosetta Stone' to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org/biochem and KBase.

Distinct microbes, metabolites, and ecologies define the microbiome in deficient and proficient mismatch repair colorectal cancers.

BACKGROUND: Links between colorectal cancer (CRC) and the gut microbiome have been established, but the specific microbial species and their role in carcinogenesis remain an active area of inquiry. Our understanding would be enhanced by better accounting for tumor subtype, microbial community interactions, metabolism, and ecology. METHODS: We collected paired colon tumor and normal-adjacent tissue and mucosa samples from 83 individuals who underwent partial or total colectomies for CRC. Mismatch repair (MMR) status was determined in each tumor sample and classified as either deficient MMR (dMMR) or proficient MMR (pMMR) tumor subtypes. Samples underwent 16S rRNA gene sequencing and a subset of samples from 50 individuals were submitted for targeted metabolomic analysis to quantify amino acids and short-chain fatty acids. A PERMANOVA was used to identify the biological variables that explained variance within the microbial communities. dMMR and pMMR microbial communities were then analyzed separately using a generalized linear mixed effects model that accounted for MMR status, sample location, intra-subject variability, and read depth. Genome-scale metabolic models were then used to generate microbial interaction networks for dMMR and pMMR microbial communities. We assessed global network properties as well as the metabolic influence of each microbe within the dMMR and pMMR networks. RESULTS: We demonstrate distinct roles for microbes in dMMR and pMMR CRC. Bacteroides fragilis and sulfidogenic Fusobacterium nucleatum were significantly enriched in dMMR CRC, but not pMMR CRC. These findings were further supported by metabolic modeling and metabolomics indicating suppression of B. fragilis in pMMR CRC and increased production of amino acid proxies for hydrogen sulfide in dMMR CRC. CONCLUSIONS: Integrating tumor biology and microbial ecology highlighted distinct microbial, metabolic, and ecological properties unique to dMMR and pMMR CRC. This approach could critically improve our ability to define, predict, prevent, and treat colorectal cancers.

Inferring Metabolic Mechanisms of Interaction within a Defined Gut Microbiota.

The diversity and number of species present within microbial communities create the potential for a multitude of interspecies metabolic interactions. Here, we develop, apply, and experimentally test a framework for inferring metabolic mechanisms associated with interspecies interactions. We perform pairwise growth and metabolome profiling of co-cultures of strains from a model mouse microbiota. We then apply our framework to dissect emergent metabolic behaviors that occur in co-culture. Based on one of the inferences from this framework, we identify and interrogate an amino acid cross-feeding interaction and validate that the proposed interaction leads to a growth benefit in vitro. Our results reveal the type and extent of emergent metabolic behavior in microbial communities composed of gut microbes. We focus on growth-modulating interactions, but the framework can be applied to interspecies interactions that modulate any phenotype of interest within microbial communities.

Synthesis of multi-omic data and community metabolic models reveals insights into the role of hydrogen sulfide in colon cancer.

Multi-omic data and genome-scale microbial metabolic models have allowed us to examine microbial communities, community function, and interactions in ways that were not available to us historically. Now, one of our biggest challenges is determining how to integrate data and maximize data potential. Our study demonstrates one way in which to test a hypothesis by combining multi-omic data and community metabolic models. Specifically, we assess hydrogen sulfide production in colorectal cancer based on stool, mucosa, and tissue samples collected on and off the tumor site within the same individuals. 16S rRNA microbial community and abundance data were used to select and inform the metabolic models. We then used MICOM, an open source platform, to track the metabolic flux of hydrogen sulfide through a defined microbial community that either represented on-tumor or off-tumor sample communities. We also performed targeted and untargeted metabolomics, and used the former to quantitatively evaluate our model predictions. A deeper look at the models identified several unexpected but feasible reactions, microbes, and microbial interactions involved in hydrogen sulfide production for which our 16S and metabolomic data could not account. These results will guide future in vitro, in vivo, and in silico tests to establish why hydrogen sulfide production is increased in tumor tissue.

ProbAnnoWeb and ProbAnnoPy: probabilistic annotation and gap-filling of metabolic reconstructions.

Summary: Gap-filling is a necessary step to produce quality genome-scale metabolic reconstructions capable of flux-balance simulation. Most available gap-filling tools use an organism-agnostic approach, where reactions are selected from a database to fill gaps without consideration of the target organism. Conversely, our likelihood based gap-filling with probabilistic annotations selects candidate reactions based on a likelihood score derived specifically from the target organism's genome. Here, we present two new implementations of probabilistic annotation and likelihood based gap-filling: a web service called ProbAnnoWeb, and a standalone python package called ProbAnnoPy. Availability and implementation: Our tools are available as a web service with no installation needed (ProbAnnoWeb) at probannoweb.systemsbiology.net, and as a local python package implementation (ProbAnnoPy) at github.com/PriceLab/probannopy. Contact: evangelos.simeonidis@systemsbiology.org or nathan.price@systemsbiology.org. Supplementary information: Supplementary data are available at Bioinformatics online.

Mackinac: a bridge between ModelSEED and COBRApy to generate and analyze genome-scale metabolic models.

SUMMARY: Reconstructing and analyzing a large number of genome-scale metabolic models is a fundamental part of the integrated study of microbial communities; however, two of the most widely used frameworks for building and analyzing models use different metabolic network representations. Here we describe Mackinac, a Python package that combines ModelSEED's ability to automatically reconstruct metabolic models with COBRApy's advanced analysis capabilities to bridge the differences between the two frameworks and facilitate the study of the metabolic potential of microorganisms. AVAILABILITY AND IMPLEMENTATION: This package works with Python 2.7, 3.4, and 3.5 on MacOS, Linux and Windows. The source code is available from https://github.com/mmundy42/mackinac . CONTACT: mundy.michael@mayo.edu or soares.maria@mayo.edu.

MMinte: an application for predicting metabolic interactions among the microbial species in a community.

BACKGROUND: The explosive growth of microbiome research has yielded great quantities of data. These data provide us with many answers, but raise just as many questions. 16S rDNA-the backbone of microbiome analyses-allows us to assess α-diversity, ß-diversity, and microbe-microbe associations, which characterize the overall properties of an ecosystem. However, we are still unable to use 16S rDNA data to directly assess the microbe-microbe and microbe-environment interactions that determine the broader ecology of that system. Thus, properties such as competition, cooperation, and nutrient conditions remain insufficiently analyzed. Here, we apply predictive community metabolic models of microbes identified with 16S rDNA data to probe the ecology of microbial communities. RESULTS: We developed a methodology for the large-scale assessment of microbial metabolic interactions (MMinte) from 16S rDNA data. MMinte assesses the relative growth rates of interacting pairs of organisms within a community metabolic network and whether that interaction has a positive or negative effect. Moreover, MMinte's simulations take into account the nutritional environment, which plays a strong role in determining the metabolism of individual microbes. We present two case studies that demonstrate the utility of this software. In the first, we show how diet influences the nature of the microbe-microbe interactions. In the second, we use MMinte's modular feature set to better understand how the growth of Desulfovibrio piger is affected by, and affects the growth of, other members in a simplified gut community under metabolic conditions suggested to be determinant for their dynamics. CONCLUSION: By applying metabolic models to commonly available sequence data, MMinte grants the user insight into the metabolic relationships between microbes, highlighting important features that may relate to ecological stability, susceptibility, and cross-feeding. These relationships are at the foundation of a wide range of ecological questions that impact our ability to understand problems such as microbially-derived toxicity in colon cancer.

Effect of individualized social activities on quality of life among older adults with mild to moderate cognitive impairment in a geriatric psychiatry facility.

OBJECTIVES: The study examined the effect of an individualized social activities intervention (ISAI) on quality of life among older adults with mild to moderate cognitive impairment in a geriatric psychiatry facility. METHOD: This randomized control trial consisted of 52 older adults (M = 70.63, SD = 5.62) with mild to moderate cognitive impairment in a geriatric inpatient psychiatry facility. A 2 (group condition) × 2 (time of measurement) design was used to compare the control (treatment-as-usual) and intervention (treatment-as-usual plus ISAI) conditions at pre- and post-treatment. ISAI consisted of 30- to 60-minute sessions for up to 15 consecutive days. The Dementia Quality of Life instrument and Neurobehavioral Rating Scale-Revised were used to examine quality of life and behavioral and psychological symptoms of dementia at pre- and post-treatment. RESULTS: Intent-to-treat analyses indicated a significant time × group condition interaction on quality of life, with this effect remaining when only completer data were included. There was no evidence of a significant treatment effect on behavioral and psychological symptoms of dementia. CONCLUSION: Findings suggest that individualized social activities are a promising treatment for cognitively impaired geriatric inpatients.

Predictors in geriatric psychiatry hospital length of stay.

This paper examined predictors of length of stay in a freestanding geriatric psychiatry hospital. Data on patient and treatment characteristics of geriatric inpatients (N = 1,593) were extracted from an archival administrative tracking database from Mary Starke Geriatric Harper Center. Five independent variables (length of time between last discharge and most recent admission, number of previous admissions, number of assaults, co-morbid medical condition, and admitting psychiatric diagnosis) were entered into a hierarchical regression model as potential predictors of length of stay in a geriatric psychiatry hospital. Number of assaults committed by the patient was the only significant predictor of length of stay, such that patients that had a greater number of assaults were more likely to have longer lengths of stay than those with fewer assaults. These findings highlight the importance of identifying patients at risk for assaultive behavior and developing effective interventions for aggression in geriatric psychiatry hospitals.

Likelihood-based gene annotations for gap filling and quality assessment in genome-scale metabolic models.

Genome-scale metabolic models provide a powerful means to harness information from genomes to deepen biological insights. With exponentially increasing sequencing capacity, there is an enormous need for automated reconstruction techniques that can provide more accurate models in a short time frame. Current methods for automated metabolic network reconstruction rely on gene and reaction annotations to build draft metabolic networks and algorithms to fill gaps in these networks. However, automated reconstruction is hampered by database inconsistencies, incorrect annotations, and gap filling largely without considering genomic information. Here we develop an approach for applying genomic information to predict alternative functions for genes and estimate their likelihoods from sequence homology. We show that computed likelihood values were significantly higher for annotations found in manually curated metabolic networks than those that were not. We then apply these alternative functional predictions to estimate reaction likelihoods, which are used in a new gap filling approach called likelihood-based gap filling to predict more genomically consistent solutions. To validate the likelihood-based gap filling approach, we applied it to models where essential pathways were removed, finding that likelihood-based gap filling identified more biologically relevant solutions than parsimony-based gap filling approaches. We also demonstrate that models gap filled using likelihood-based gap filling provide greater coverage and genomic consistency with metabolic gene functions compared to parsimony-based approaches. Interestingly, despite these findings, we found that likelihoods did not significantly affect consistency of gap filled models with Biolog and knockout lethality data. This indicates that the phenotype data alone cannot necessarily be used to discriminate between alternative solutions for gap filling and therefore, that the use of other information is necessary to obtain a more accurate network. All described workflows are implemented as part of the DOE Systems Biology Knowledgebase (KBase) and are publicly available via API or command-line web interface.