Next Article in Journal
Knockout of Arabidopsis thaliana VEP1, Encoding a PRISE (Progesterone 5β-Reductase/Iridoid Synthase-Like Enzyme), Leads to Metabolic Changes in Response to Exogenous Methyl Vinyl Ketone (MVK)
Next Article in Special Issue
A Checklist for Reproducible Computational Analysis in Clinical Metabolomics Research
Previous Article in Journal
Modelling Metabolic Shifts during Cardiomyocyte Differentiation, Iron Deficiency and Transferrin Rescue Using Human Pluripotent Stem Cells
Previous Article in Special Issue
MSCAT: A Machine Learning Assisted Catalog of Metabolomics Software Tools
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 1
10.3390/metabo12010010
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TrpNet: Understanding Tryptophan Metabolism across Gut Microbiome

by
Yao Lu
1,
Jasmine Chong
2,
Shiqian Shen
3,
Joey-Bahige Chammas
4,5,
Lorraine Chalifour
4,5 and
Jianguo Xia
1,2,*
1
Department of Microbiology and Immunology, McGill University, Montreal, QC H3A 2T5, Canada
2
Institute of Parasitology, McGill University, Montreal, QC H3A 2T5, Canada
3
MGH Center for Translational Pain Research, Massachusetts General Hospital, Harvard Medical School, Boston, MA 02215, USA
4
Department of Medicine, McGill University, Montreal, QC H3A 2T5, Canada
5
Lady Davis Institute for Medical Research, Montreal, QC H3T 1E2, Canada
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(1), 10; https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12010010
Submission received: 23 November 2021 / Revised: 20 December 2021 / Accepted: 20 December 2021 / Published: 23 December 2021
(This article belongs to the Special Issue Advances in Computational Metabolomics)
Browse Figures
Versions Notes

Abstract

:
Crosstalk between the gut microbiome and the host plays an important role in animal development and health. Small compounds are key mediators in this host–gut microbiome dialogue. For instance, tryptophan metabolites, generated by biotransformation of tryptophan through complex host–microbiome co-metabolism can trigger immune, metabolic, and neuronal effects at local and distant sites. However, the origin of tryptophan metabolites and the underlying tryptophan metabolic pathway(s) are not well characterized in the current literature. A large number of the microbial contributors of tryptophan metabolism remain unknown, and there is a growing interest in predicting tryptophan metabolites for a given microbiome. Here, we introduce TrpNet, a comprehensive database and analytics platform dedicated to tryptophan metabolism within the context of host (human and mouse) and gut microbiome interactions. TrpNet contains data on tryptophan metabolism involving 130 reactions, 108 metabolites and 91 enzymes across 1246 human gut bacterial species and 88 mouse gut bacterial species. Users can browse, search, and highlight the tryptophan metabolic pathway, as well as predict tryptophan metabolites on the basis of a given taxonomy profile using a Bayesian logistic regression model. We validated our approach using two gut microbiome metabolomics studies and demonstrated that TrpNet was able to better predict alterations in in indole derivatives compared to other established methods.

Graphical Abstract

1. Introduction

The gut microbiome is a community of metabolically active microorganisms inhabiting all niches along the intestines that coevolves with its host. Growing evidence has shown that the gut microbiome plays a critical role in animal development and health [1]. Disruptions in microbiome composition, termed dysbiosis, are implicated in various diseases including gastrointestinal diseases [2], infectious diseases [3], metabolic diseases [4,5], and neurological disorders [6]. Dysbiosis leads to a shift in the production of various microbial metabolites which then influence the physiology and immune status of the host [7]. Among these bioactive metabolites, short-chain fatty acids (SCFAs, produced by bacteria from fermenting dietary fibers), secondary bile acids (originated in liver and transformed by gut microbiome), and tryptophan-derived metabolites are most well known.
Microbes can degrade tryptophan to a range of indoles including indolelactate (ILA), indoleacetic acid (IAA), indolealdehyde (IAld), indoleacrylic acid (IA), and indolepropionate (IPA). These can activate the aryl hydrocarbon receptor (AhR), a transcription factor widely expressed by cells in the immune system, regulate intestinal homeostasis [8], initiate an immune response [9], and control oxidative stress defense [10]. AhR activation is associated with multiple diseases such as inflammatory bowel disease (IBD) [11], type 2 diabetes [12], and central nervous system (CNS)-related disorders [13].
To understand the health impact of tryptophan metabolites in host–microbiome interactions, it is essential to have detailed knowledge of tryptophan metabolism. Researchers can find related reactions and enzymes in several databases including Kyoto Encyclopedia of Genes and Genomes (KEGG) [14], BioCyc [15], Small Molecule Pathway Database (SMPDB) [16], and WikiPathways [17]. Among them, SMPDB and WikiPathways are mainly concerned with host tryptophan metabolism, while KEGG and BioCyc provide very limited information on tryptophan degradation across the gut microbiome. In particular, the origins of tryptophan metabolism are assigned at the enzyme level, and it is difficult to obtain the reactions and metabolites of a single species or from the microbial community. Moreover, not all tryptophan metabolites can be found in these databases. For instance, IPA, an important neuroprotective antioxidant produced by the human gut microbe represented by Clostridium sporogenes, is currently missing in KEGG.
Predicting metagenomic functions on the basis of microbiome composition has attracted great attention in recent years. Current tools, such as PICRUSt [18] and Tax4Fun [19], focus on enzyme- and pathway-level predictions and cannot be directly used to understand the biological effects driven by metabolites. In addition, the bias and missing information in their underlying databases pose inherent limitations on more focused analyses, such as on microbial tryptophan metabolism. To elucidate the microbial contributors for tryptophan metabolites, it is necessary to have a dedicated resource covering all the existing reactions in tryptophan degradation.
The rapid advancements in omics technologies and bioinformatics have joined forces to systematically decipher the function of gut microbiome. For instance, metagenomics and metabolomics led to the discovery of thousands of microbe-derived small molecules, as well as the genes associated with their productions [20]. Meanwhile, genome-scale metabolic models (GEMs) aim to accurately capture an organism’s metabolism by integrating information obtained from genome annotation, biochemical reaction, manual curation, and literature review. High-quality GEMs are now available for numerous microorganisms and serve as good resources for large-scale investigations of microbial tryptophan metabolism [21]. Several recent studies have employed GEMs to infer the relationships between phenotypic differentiation and metabolic capacities [22], to predict drug targets [23], and to understand tryptophan-metabolizing microbes involved in murine diarrhea [24]. The tryptophan metabolism pathways can be extracted from the GEMs to help understand the tryptophan metabolism distribution across the gut microbiome.
Here, we introduce TrpNet (https://www.trpnet.ca) which systematically collects the metabolites, reactions, and enzymes involved in tryptophan metabolism with specific attention to their microbial producers. TrpNet currently describes the tryptophan degradation pathways across ~5000 bacterial species including most known gut microbial species. Users can easily navigate TrpNet to find the origins of tryptophan metabolites or to generate tryptophan metabolism networks for selected microbes. Lastly, TrpNet allows researchers to predict tryptophan metabolites from any given taxonomic profiles on the basis of a Bayesian logistics regression model.

2. Results

The overall workflow for the development of TrpNet is shown in Figure 1. The details of key steps are described in the subsequent sections.

2.1. Literature Search and Intestinal Tryptophan Metabolism

We define the tryptophan metabolism pathway as a set of reactions that transfer tryptophan to an end product without further searchable reactions or breakdown of tryptophan until the energy metabolism. To enumerate the metabolic reactions and metabolites regarding tryptophan biotransformation in microbes and their mammalian hosts, we manually searched and compared over 300 biochemical and metabolomic research papers on tryptophan metabolites, 37 reviews published in the last 5 years, and 14 public databases.
As shown in Figure 2, tryptophan metabolism generates 29 bioactive metabolites via three major pathways—indole pathway, serotonin pathway, and kynurenine pathway. The indole pathway, converting tryptophan into indole derivatives, including AhR ligands, predominates in gut microbes, while the serotonin and kynurenine pathways predominate in mammalian hosts. However, the origins of some tryptophan metabolites, whether microbe-derived or host-derived, are inconsistent across previously published reviews. In addition, most reports focused on the indirect roles of the gut microbiome in modulating kynurenine and serotonin production through non-tryptophan metabolites such as butyrate, an important SCFAs derived from gut microbes [25]. The direct production of kynurenines and serotonin by gut microbes is not well characterized. For instance, serotonin was termed as a host-limited metabolite [26], yet our investigation showed that several species such as Lactococcus lactis, Lactobacillus plantarum, and Klebsiella pneumoniae produced serotonin in a similar way to their mammalian host via aromatic amino-acid decarboxylase (AAAD) [27,28,29]. Another important neurotransmitter, tryptamine, was traditionally regarded as a microbial metabolite produced by Clostridium, Ruminococcus, Blautia, and Lactobacillus through tryptophan decarboxylases [13]. However, it is also reported to be produced by brain cells in certain cases and may play specific roles in the mammalian brain [30]. In addition, some gut microbes can degrade tryptophan in a different way through the kynurenine pathway. For instance, gut species Burkholderia cepacia was reported to convert tryptophan to 2-amino-3-carboxymuconate semialdehyde, which was further enzymatically degraded to pyruvate and acetate via the intermediates 2-aminomuconate and 4-oxalocrotonate rather than the known mammalian pathway which transforms 2-aminomuconate to 2-ketoadipate and, ultimately, glutaryl-coenzyme [31]. Compared with the most recent reviews [26,32,33,34] on the bioactive tryptophan metabolites, we updated the origin of all collected tryptophan metabolites including three inconsistent annotations of the kynurenines according to the current literature searches. The results were further cross-validated and enhanced with the information obtained from mining the GEMs, as described below.

2.2. Curation of Genome-Scale Metabolic Models

GEMs are knowledge-based stoichiometric-balanced metabolic networks containing the entire set of metabolic reactions, genes, and metabolites in the target organism [35]. Current developments in systems biology allow for the large-scale reconstruction of GEMs for numerous microorganisms. For instance, AGORA is a set of semiautomatically generated GEMs for 818 gut bacteria [36], and EMBL_GEMs is another large collection (5584 bacteria) for all reference and representative bacterial genomes of NCBI RefSeq [37] using CarveMe [38]. The reconstruction tools for both AGORA (assembly of gut organisms through reconstruction and analysis) and EMBL_GEMs were evaluated outstanding among the general tools, especially in gap-filling the network [39].
A total of 6402 GEMs covering 41 phyla were collected from AGORA and EMBL_GEMs. Most GEMs are at a strain level except for 73 at the species level and 333 models belonging to same strains shared between the two datasets. GEMs were manually annotated according to literature searches [40,41,42,43,44], of which 2114 models were labeled as the human gut microbe covering 1380 species of 30 phyla, and 177 were part of the mouse gut microbiome from 98 species of 10 phyla. The reactions, metabolites, and enzymes involved in microbial tryptophan metabolism were extracted from GEMs. These include nearly 5000 species belonging to 39 phyla and involve tryptophan metabolism covering 1246 species in the human gut and 88 species in the mouse gut (Figure 3). The results were corroborated by literature searches to reconcile inconsistencies between AGORA and EMBL_GEMs, as well as to make the GEM data as complete as possible.

2.3. Development of a Database for Tryptophan Metabolism and Functional Prediction

Following two major procedures described above, the final tryptophan metabolism pathway contains the entries for 130 reactions and 108 metabolites (excluding currency compounds such as water, hydrogen, oxygen, etc.) linking to 91 enzymes and more than 5000 GEMs. We developed a user-friendly web-based database and visual analytics tool - TrpNet (https://www.trpnet.ca/, accessed on 23 November 2021) to share this resource with the community. Users can browse, search, and filter reactions, metabolites, or microbes involved in tryptophan metabolism and visualize more detailed information and summary tables in multiple formats. Whenever possible, different entries are hyperlinked to PubMed, KEGG [14], BioCyc [15], and ModelSEED [21].
A main motivation of developing TrpNet is to help understand the relationship between the gut microbiome composition and the capacity for tryptophan metabolism. We designed the interface and functions to allow users to easily obtain the distribution of tryptophan metabolite production at different taxonomy levels. Figure 4 shows the pairwise distance between the phylogenetic tree from the dominant genus in the host gut and the corresponding metabolic clusters, according to the presence or absence of tryptophan metabolite production. It can be observed that phylogenetically close species may differ in their capacities in metabolite production. These data will help to resolve some inconsistencies between microbiome and metabolome divergence and the coexistence of specific species [22,45]. For instance, Bacteroides were found to be relatively conservative while Lactobacillus fluctuated in tryptophan metabolite production depending on whether they produced indole derivatives. Human- and mouse-specific gut microbes differed in the production of several AhR ligands such as IA, IAA, ILA, and IPA. This may help explain the different affinities of human AhR and mouse AhR in selecting exogenous ligands as reported in several studies [46,47] and shown in Figure S1.
Several computational tools such as PICRUSt2 [18] and Tax4Fun2 [19] are available for predicting functional profiles from 16S rRNA gene sequence data. Their performances are inherently limited by the known annotated enzyme groups which may not represent the metabolite generation. Specifically, public databases used in current tools are not tailored for tryptophan metabolism, and this may lead to bias due to incomplete information. TrpNet provides a more complete tryptophan metabolism according to literature curation and GEMs that describe metabolism at strain level with the potential to predict unknown enzymatic reactions. Here, we explored whether we could better predict the microbial tryptophan metabolism using the TrpNet database.
One constraint is that 16S rRNA data cannot reach the resolution of strain level but usually identify the microbiome at the genus level. To address this issue, we used a logistic regression model to estimate the tryptophan metabolite production potential of the interested genus depending on the metabolite distribution collected by TrpNet. This approach was used in previous studies [48,49,50] to model microbiome compositional data and to identify informative microbiome features. To acquire more accurate models for our prediction, we fit Bayesian logistic regression models for each tryptophan metabolite according to their distributions across the taxonomy levels. In this model, the human/mouse gut origin was included as a nonrandom covariate as tryptophan metabolite production differs by the niche. Table 1 and Table 2 show the estimated odds ratios for the prevalence genus in producing bioactive indoles generated from mouse model and human model, respectively. The models were firstly validated by randomly split TrpNet database, whereby 80% was used for training and 20% was used to evaluate the model performance. We found that genus levels provided relatively reliable results for different metabolites in general. Figure 5 shows the ROC curves of the prediction models comparing different taxonomic levels in predicting IAA production. Please note that the performance measures are likely to be inflated as the same database was used for calculating the parameters of the regression models.
A network visualization page was implemented to allow users to search metabolites of interest in the network or to customize the tryptophan metabolism network according to a user-specified list of microbes (Figure 6). The result can be highlighted against the whole network or downloaded as a table. Another key feature of TrpNet lies in the annotation for the origin beyond the enzyme level. Reactions and metabolites were individually checked against the literature to label them as host-derived or microbial-derived, to help decipher the host–microbe interactions and co-metabolism.

2.4. Case Studies

2.4.1. Myocardial Infarct (MI) Case Study

Disturbed tryptophan metabolism is known to alter the host inflammation status and affect many diseases including heart diseases such as myocardial infarction (MI) with an increased ratio of KYN/TRP [51]. To understand gut microbiome and host MI status with tryptophan metabolism, we collected 16 cecal samples from 16 mice (8 with MI and 8 control) day 3 post MI. Each sample was processed for 16S rRNA bacterial sequencing and untargeted metabolomics based on LC–MS and MS/MS. As it has been reported that females and males have differences in the risk of MI, we included data from male mice in our case study to exclude any additional effects of sex [52].
DADA2 [53] was used to assign taxonomy to amplicon sequence variants (ASVs). After filtering the 712 low-quality features, the remaining 304 ASVs were attributed to 69 genera dominated by Lachnospiraceae spp. And Ruminococcaceae spp. For metabolomics data, XCMS [54] and metID [55] were used for spectrum processing and peak annotation. A total of 24 microbial tryptophan metabolites were detected in LC–MS/MS, of which nine metabolites were significantly different including IAA, IAM, IalD, and serotonin. Statistical analyses of microbiome data were performed using MicrobiomeAnalyst [56]. Principal component analysis (PCA) evaluation showed that male mice without an MI differed from male mice post MI in microbiome composition. This was caused by a lack of Proteobacteria and Verrucomicrobia, which are active tryptophan metabolites producers in the no MI group (Figure S2).
Prediction models built on the TrpNet database were used to predict tryptophan metabolite production as a function of the genus-level data. Figure 7 shows the prediction result from the gut microbiome, as well as the comparison with metabolomics data and related enzymes predicted by PICRUSt2. Tryptophan degradation of the MI group differed significantly from their counterparts without MI, which may be explained by their diverse gut microbiome composition. According to our prediction, the MI group is more likely to produce greater amounts of tryptophan metabolites, including AhR ligands such as indole, IAM, and IAA, supporting the metabolomics data. Previous evidence has suggested that AhR activity is a critical modulator in the development and pathogenesis of the cardiovascular system [57]. AhR knockout mice were reported to be more susceptible to cardiac hypertrophy, vascular remodeling and systemic hypertension [58]. However, AhR activation can also contribute to the formation and promotion of atherosclerosis through inducing vascular inflammation [59]. Further studies are necessary in order to elucidate the effects on MI progression triggered by microbial AhR ligands from tryptophan metabolism.
Our prediction also shows increased activation of the kynurenine pathway post MI. This suggests that gut microbes may directly contribute to the increased KYN/TRP ratio, leading to the decreased level of beneficial serotonins and accumulation of neurotoxic KYN metabolites during the disease process. In parallel, we performed analyses using PICRUSt2. Only enzyme EC 4.1.99.1 relating to indole production was identified as significantly increased in the MI group, similar to the prediction by TrpNet and metabolomics data. Thus, TrpNet can serve as a better resource for exploring intestinal tryptophan metabolism.

2.4.2. IBD Case Study

Previous studies have demonstrated the key role of the gut microbiome in IBD, and some highlighted the potential link to gut tryptophan metabolism [60]. The 16S rRNA and metabolomics data were collected from 26 participants between age 6 and 19 randomly selected from the Inflammatory Bowel Disease Multi-omics Database (http://ibdmdb.org, accessed on 17 October 2021) [61]. For each data type, 20 samples from pediatric Crohn’s disease (CD) patients and 20 from pediatric healthy controls were also included.
From the metabolomics data annotation, nine tryptophan-derived metabolites were observed among which seven could be produced by the microbiome. IPA was significantly decreased in the CD group. Regarding the ASV sequencing data, 147 ASVs were annotated to 44 genera after filtering out the low-abundance features. However, there were no significant differences observed regarding microbiome composition between the CD and control group (Figure S3). TrpNet was then used for tryptophan metabolite prediction for each sample using the established model. Figure 8 shows the predicted distribution of tryptophan metabolites in comparison with metabolomics data and EC identified by PICRUSt2. Our prediction found the alteration of IPA validated by metabolomics and the obligatory intermediate IA in producing IPA from tryptophan. In contrast, PICRUSt2 did not contain the information for IPA, and the enzyme for IA was not significantly different between CD patients and healthy controls. Indole derivatives were predicted by TrpNet to be more abundant in healthy people than in CD patients, which is consistent with a previous report showing a reduction in AhR ligands by the microbiota in IBD patients. Most metabolites did not show significant differences between the two groups, probably due to the disperse microbiome structure. Interestingly, although indolepyruvate, which can improve intestinal epithelial barrier function during challenges with inflammatory stimuli, was not annotated by the metabolomics data, our prediction shows its decrease in CD patients, replicating previous results [60]. Despite serotonin being found increased in the CD group by metabolomics analysis, which was possibly due to the decreased expression of SERT in the ileum and colon [62], there was no significant difference according to our prediction. Similarly, the kynurenines were not different between the two groups using all the methods. Consequently, we can envision that gut microbes may affect IBD processing through tryptophan-derived AhR ligands such as IPA and IPY.

3. Discussion

Tryptophan metabolism plays a central role in host physiologic and pathologic processes. The balance among microbial tryptophan metabolism, supplementation, and microbial modulation exerts a great impact on local gastrointestinal and circulating tryptophan availability for its host and ultimately contributes to host health and disease. Hence, it is important to fully characterize tryptophan metabolism within a host or within its resident gut microbes. TrpNet, a first step toward addressing this gap, includes a collection of all currently known reactions and metabolites relating to tryptophan according to comprehensive literature reviews and large-scale data mining across >5000 GEMs. However, despite our intensive curation efforts, several reactions and metabolites are still left without related literature reports. For example, no reaction details are currently available for several tryptophan metabolites such as for Iald, an important AhR ligand.
One of the major challenges in microbiome studies is to determine the causal role that the gut microbiome composition plays in specific phenotypes. This is difficult due to the complexity of host–microbe interactions and microbe–microbe interactions. TrpNet can help decipher this co-metabolism by providing the detailed tryptophan metabolism within specific microbial species according to GEMs and literature annotation. Many current studies are based on 16S rRNA sequencing, making it essential to improve functional prediction and maximize the information gained from these relatively low-resolution taxonomic profiles. Here, we took an initial trial to predict tryptophan metabolism from genus-level bacterial identification using a logit regression model based on the TrpNet database. It should be noted that this prediction is limited by the current knowledge of the tryptophan metabolism, as well as algorithm for GEM construction or function prediction. Optimized methods are needed to improve the annotation of microbe to metabolite levels for mechanical and therapeutical insights. For instance, an increased KYN/Trp ratio has been reported as a potential biomarker for inflammation status, and supplementation of gut species that can naturally produce AhR ligands such as Lactobacillus spp. could help recover the AhR signaling. This microbe-based therapeutic approach was successfully applied in a mouse model of colitis [11]. As the gut microbiome can also modulate tryptophan metabolism indirectly by producing other small molecules such as bile acids, it will be useful to gather the information of microbes involved in these relevant processes to further improve TrpNet.

4. Materials and Methods

4.1. Literature Review

Review papers were searched from PubMed, Web of Science, and bioRxiv (www.biorxiv.org/, accessed on 8 October 2021) using the search term “tryptophan metabolism AND gut microbiome” since 2017. Those studies providing a global view of and tryptophan metabolism and focusing on the host–microbe interaction were included. Furthermore, for each tryptophan metabolite, research paper surveys were conducted to determine its origin. These papers showed showing at least one of the sources of genetic, enzymatic, or metabolic evidence in certain microbial species were prioritized.

4.2. GEMs Collection

A total of 818 GEMs in AGORA were collected from the Virtual Metabolic Human (VMH) database that can be accessed via the website (http://vmh.life, accessed on 3 September 2021), and EMBL GEMs were download from the EMBL BioModels website (https://www.ebi.ac.uk/biomodels, accessed on 3 September 2021). SBML files were parsed using R studio (version 4.1.1). GEMs were first annotated to human and/or mouse gut microbes on the basis of several large-scale studies and public gut microbiome databases. The models without records were then manually searched in PubMed to annotate their habitat.

4.3. TrpNet Implementation

The web-based database was developed on the basis of the JavaServer Faces (JSF) technology using the PrimeFaces framework (v11). The network visualization was implemented using D3 (version 5.0).

4.4. Sample Collection for MI Case Study

The murine experiments were approved by the Lady Davis Facility Animal Care Committee and followed the guidelines described by the Canadian Council on Animal Care. Retired breeder male mice were purchased from Charles River, St. Constant, PQ, Canada. Mice were housed in single cages on irradiated corn cob bedding in a vented rack, fed an irradiated Harlan Teklad Global 2018 diet which contains no animal protein and acidified tap water, and acclimated to the facility for 1 month before use. Surgery to create an MI was performed by the Surgery Core of the Lady Davis Institute [63,64]. Samples of cecal contents were collected day 3 post MI from a total of eight male mice, as well as from eight male mice which did not experience MI surgery. DNA for 16S rRNA sequencing was isolated using a Qiagen QIAamp PowerFecal DNA kit according to the manufacturer’s instructions. DNA samples were quantified, purity was determined, and samples sent to the McGill Genome Center. There, the bacterial V4 region was PCR-amplified from bases 515F to 806R, sequenced using a MiSeq Reagent Kit v3 (600-cycle), and run on an Illumina MiSeq. Data were processed and returned as amplified sequence variants (ASVs). R package DADA2 v 1.20.0 [53] was used to determine the abundance and gut bacterial species assignment. Gut metabolomics data from the same cecal contents were processed using an Orbitrap Q-Exactive LC–MS system in both positive and negative mode using a C18 column. MS/MS spectra were collected using data-independent acquisition (DIA). Raw LC–MS spectra were processed by MetaboAnalyst v5.0 [65] to generate a peak list table. About 150 MS1 peaks were found to be from potential tryptophan metabolites. According to MS/MS data, 24 tryptophan metabolites were identified using the metID package [52] (Table S1).

4.5. Sample Collection for IBD Case Study

The dataset of pediatric IBD stool samples was downloaded from the Integrative Human Microbiome Project Consortium (iHMP) [66]. For evaluation purposes, we randomly selected individuals between age 6 and 19 for disease (diagnosed with Crohn’s disease) and control groups. The information of the sample is listed in Table S2, and the original data can be found at https://ibdmdb.org/ (accessed on 17 October 2021). The tryptophan metabolites were extracted on the basis of annotation information provided by the authors (Table S3).

4.6. Logistic Regression Model for Predicting Metabolite Profiles

The logistic regression model was used to infer the metabolic profile from known taxonomy compositions. This method is from the generalized linear model family and can learn probabilistic models to predict the outcome of a binary variable from one or more response categorical or continuous variables. In our case, we aimed to predict tryptophan metabolite production using taxonomy profile and host type. The algorithm involved four key steps as described below.
  • Different taxonomy levels and their combinations were evaluated for their predictive values. Models were ranked by Akaike information criterion (AIC). The genus level combined with the host type was selected as the best predictor;
  • The models were further optimized by Bayesian logistic regression coupled with a fast Pareto smoothed leave-one-out cross-validation for the penalized likelihood estimation [67]. These models capture the metabolite production potential (PM, G) for the underlying metabolite (M) of interest in every genus (G) for a given host type;
  • The predicted probability (PM,G) was multiplied by the genus abundance table obtained from 16S rRNA sequencing data to compute the accumulated production potential for each metabolite of interest for each sample;
  • The results of all samples were normalized by total sum scaling to be comparable with metabolomics data.

5. Conclusions

Understanding molecular dialogues between the gut microbiome and the host is critical for developing microbiome-based diagnostic and therapeutic approaches. In this manuscript, we focused on improving our knowledge on tryptophan metabolism by integrating information from >5000 GEMs, 14 databases, and >300 literature reports. Through its user-friendly interface and interactive visualization, TrpNet provides the most up-to-date information for researchers to study tryptophan metabolism within the context of host and microbiome interactions. According to this information, we further developed an algorithm for predicting the microbial tryptophan metabolism from the 16S rRNA abundance profiles. Our two case studies demonstrated that our approach gives more accurate results compared to other established methods. We hope that TrpNet will be a useful resource that allows researchers to better understand the gut microbial tryptophan metabolism in the context of the gut microbiome for translational applications.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo12010010/s1. Figure S1. Observed distribution of tryptophan metabolites production across the predominant genus in host gut. Microbial tryptophan metabolite production differs according to their host niche, Figure S2. Comparation of microbiome composition between MI and NO.MI group by (a) PCA plot and (b) bar plot, Figure S3. Comparation of microbiome composition between CD patients and health control by (a) PCA plot and (b) stack bar plot, Table S1. Metabolite annotation of MI case study, Table S2. IBD sample information, Table S3. Metabolite annotation of IBD case study.

Author Contributions

Conceptualization, J.X.; methodology, Y.L. and J.X.; software, Y.L. and J.X.; data curation, Y.L., J.C. and J.-B.C.; writing—original draft preparation, Y.L.; writing—review and editing, J.X., S.S., L.C. and J.C.; supervision, J.X.; funding acquisition J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant and McGill Interdisciplinary Initiative in Infection and Immunity (MI4) Grant (J.X.), NIH R35 GM128692 and RF1AG070141 (S.S.), and the Heart and Stroke Foundation of Canada (L.C.). Y. Lu was supported by a PhD Scholarship from the China Scholarship Council.

Institutional Review Board Statement

The animal study protocol was approved by the Animal Care Committee of McGill University (protocol code: 2001-3634, 4 January 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The IBD data is available from http://ibdmdb.org (accessed on 23 November 2021). The MI data is available upon request.

Acknowledgments

The authors acknowledge Surgery Core of the Lady Davis Institute for the performance of the surgery to create the MI model. The McGill Genome Center and Ken Dewar are acknowledged for the amplification and sequencing of the 16S rRNA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dethlefsen, L.; McFall-Ngai, M.; Relman, D.A. An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature 2007, 449, 811–818. [Google Scholar] [CrossRef]
  2. Chassaing, B.; Darfeuille-Michaud, A. The commensal microbiota and enteropathogens in the pathogenesis of inflammatory bowel diseases. Gastroenterology 2011, 140, 1720–1728. [Google Scholar] [CrossRef] [PubMed]
  3. Khan, N.; Mendonca, L.; Dhariwal, A.; Fontes, G.; Menzies, D.; Xia, J.; Divangahi, M.; King, I.L. Intestinal dysbiosis compromises alveolar macrophage immunity to Mycobacterium tuberculosis. Mucosal. Immunol. 2019, 12, 772–783. [Google Scholar] [CrossRef]
  4. Qin, J.; Li, Y.; Cai, Z.; Li, S.; Zhu, J.; Zhang, F.; Liang, S.; Zhang, W.; Guan, Y.; Shen, D.; et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 2012, 490, 55–60. [Google Scholar] [CrossRef]
  5. Ley, R.E.; Turnbaugh, P.J.; Klein, S.; Gordon, J.I. Microbial ecology: Human gut microbes associated with obesity. Nature 2006, 444, 1022–1023. [Google Scholar] [CrossRef]
  6. Sampson, T.R.; Mazmanian, S.K. Control of brain development, function, and behavior by the microbiome. Cell Host Microbe 2015, 17, 565–576. [Google Scholar] [CrossRef] [Green Version]
  7. Blacher, E.; Levy, M.; Tatirovsky, E.; Elinav, E. Microbiome-Modulated Metabolites at the Interface of Host Immunity. J. Immunol. 2017, 198, 572–580. [Google Scholar] [CrossRef] [Green Version]
  8. Alexeev, E.E.; Lanis, J.M.; Kao, D.J.; Campbell, E.L.; Kelly, C.J.; Battista, K.D.; Gerich, M.E.; Jenkins, B.R.; Walk, S.T.; Kominsky, D.J.; et al. Microbiota-Derived Indole Metabolites Promote Human and Murine Intestinal Homeostasis through Regulation of Interleukin-10 Receptor. Am. J. Pathol. 2018, 188, 1183–1194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Stockinger, B.; Hirota, K.; Duarte, J.; Veldhoen, M. External influences on the immune system via activation of the aryl hydrocarbon receptor. Semin. Immunol. 2011, 23, 99–105. [Google Scholar] [CrossRef]
  10. Jaichander, P.; Selvarajan, K.; Garelnabi, M.; Parthasarathy, S. Induction of paraoxonase 1 and apolipoprotein A-I gene expression by aspirin. J. Lipid Res. 2008, 49, 2142–2148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Lamas, B.; Richard, M.L.; Leducq, V.; Pham, H.P.; Michel, M.L.; Da Costa, G.; Bridonneau, C.; Jegou, S.; Hoffmann, T.W.; Natividad, J.M.; et al. CARD9 impacts colitis by altering gut microbiota metabolism of tryptophan into aryl hydrocarbon receptor ligands. Nat. Med. 2016, 22, 598–605. [Google Scholar] [CrossRef]
  12. Qi, Q.; Li, J.; Yu, B.; Moon, J.Y.; Chai, J.C.; Merino, J.; Hu, J.; Ruiz-Canela, M.; Rebholz, C.; Wang, Z.; et al. Host and gut microbial tryptophan metabolism and type 2 diabetes: An integrative analysis of host genetics, diet, gut microbiome and circulating metabolites in cohort studies. Gut 2021. [Google Scholar] [CrossRef]
  13. Roth, W.; Zadeh, K.; Vekariya, R.; Ge, Y.; Mohamadzadeh, M. Tryptophan Metabolism and Gut-Brain Homeostasis. Int. J. Mol. Sci. 2021, 22, 2973. [Google Scholar] [CrossRef] [PubMed]
  14. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Karp, P.D.; Billington, R.; Caspi, R.; Fulcher, C.A.; Latendresse, M.; Kothari, A.; Keseler, I.M.; Krummenacker, M.; Midford, P.E.; Ong, Q.; et al. The BioCyc collection of microbial genomes and metabolic pathways. Brief. Bioinform. 2019, 20, 1085–1093. [Google Scholar] [CrossRef] [PubMed]
  16. Jewison, T.; Su, Y.; Disfany, F.M.; Liang, Y.; Knox, C.; Maciejewski, A.; Poelzer, J.; Huynh, J.; Zhou, Y.; Arndt, D.; et al. SMPDB 2.0: Big improvements to the Small Molecule Pathway Database. Nucleic Acids Res 2014, 42, D478–D484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Martens, M.; Ammar, A.; Riutta, A.; Waagmeester, A.; Slenter, D.N.; Hanspers, K.; Miller, R.A.; Digles, D.; Lopes, E.N.; Ehrhart, F.; et al. WikiPathways: Connecting communities. Nucleic Acids Res 2021, 49, D613–D621. [Google Scholar] [CrossRef]
  18. Douglas, G.M.; Maffei, V.J.; Zaneveld, J.R.; Yurgel, S.N.; Brown, J.R.; Taylor, C.M.; Huttenhower, C.; Langille, M.G.I. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 2020, 38, 685–688. [Google Scholar] [CrossRef]
  19. Wemheuer, F.; Taylor, J.A.; Daniel, R.; Johnston, E.; Meinicke, P.; Thomas, T.; Wemheuer, B. Tax4Fun2: Prediction of habitat-specific functional profiles and functional redundancy based on 16S rRNA gene sequences. Environ. Microbiome 2020, 15, 11. [Google Scholar] [CrossRef]
  20. Han, S.; Van Treuren, W.; Fischer, C.R.; Merrill, B.D.; DeFelice, B.C.; Sanchez, J.M.; Higginbottom, S.K.; Guthrie, L.; Fall, L.A.; Dodd, D.; et al. A metabolomics pipeline for the mechanistic interrogation of the gut microbiome. Nature 2021, 595, 415–420. [Google Scholar] [CrossRef]
  21. Seaver, S.M.D.; Liu, F.; Zhang, Q.; Jeffryes, J.; Faria, J.P.; Edirisinghe, J.N.; Mundy, M.; Chia, N.; Noor, E.; Beber, M.E.; et al. The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes. Nucleic Acids Res. 2020, 49, D575–D588. [Google Scholar] [CrossRef]
  22. Bauer, E.; Laczny, C.C.; Magnusdottir, S.; Wilmes, P.; Thiele, I. Phenotypic differentiation of gastrointestinal microbes is reflected in their encoded metabolic repertoires. Microbiome 2015, 3, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Abdel-Haleem, A.M.; Hefzi, H.; Mineta, K.; Gao, X.; Gojobori, T.; Palsson, B.O.; Lewis, N.E.; Jamshidi, N. Functional interrogation of Plasmodium genus metabolism identifies species- and stage-specific differences in nutrient essentiality and drug targeting. PLoS Comput. Biol. 2018, 14, e1005895. [Google Scholar] [CrossRef]
  24. Zhang, C.; Shao, H.; Li, D.; Xiao, N.; Tan, Z. Role of tryptophan-metabolizing microbiota in mice diarrhea caused by Folium sennae extracts. BMC Microbiol. 2020, 20, 185. [Google Scholar] [CrossRef]
  25. Kennedy, P.J.; Cryan, J.F.; Dinan, T.G.; Clarke, G. Kynurenine pathway metabolism and the microbiota-gut-brain axis. Neuropharmacology 2017, 112, 399–412. [Google Scholar] [CrossRef]
  26. Agus, A.; Planchais, J.; Sokol, H. Gut Microbiota Regulation of Tryptophan Metabolism in Health and Disease. Cell Host Microbe 2018, 23, 716–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. O’Mahony, S.M.; Clarke, G.; Borre, Y.E.; Dinan, T.G.; Cryan, J.F. Serotonin, tryptophan metabolism and the brain-gut-microbiome axis. Behav. Brain Res. 2015, 277, 32–48. [Google Scholar] [CrossRef] [PubMed]
  28. Özoğul, F.; Kuley, E.; Özoğul, Y.; Özoğul, İ. The Function of Lactic Acid Bacteria on Biogenic Amines Production by Food-Borne Pathogens in Arginine Decarboxylase Broth. Food Sci. Technol. Res. 2012, 18, 795–804. [Google Scholar] [CrossRef] [Green Version]
  29. Özoğul, F. Production of biogenic amines by Morganella morganii, Klebsiella pneumoniae and Hafnia alvei using a rapid HPLC method. Eur. Food Res. Technol. 2004, 219, 465–469. [Google Scholar] [CrossRef]
  30. Kitahama, K.; Ikemoto, K.; Jouvet, A.; Nagatsu, I.; Sakamoto, N.; Pearson, J. Aromatic L-amino acid decarboxylase- and tyrosine hydroxylase-immunohistochemistry in the adult human hypothalamus. J. Chem. Neuroanat. 1998, 16, 43–55. [Google Scholar] [CrossRef]
  31. Colabroy, K.L.; Begley, T.P. Tryptophan catabolism: Identification and characterization of a new degradative pathway. J. Bacteriol. 2005, 187, 7866–7869. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, Y.; Hou, Y.; Wang, G.; Zheng, X.; Hao, H. Gut Microbial Metabolites of Aromatic Amino Acids as Signals in Host-Microbe Interplay. Trends Endocrinol. Metab. 2020, 31, 818–834. [Google Scholar] [CrossRef]
  33. Roager, H.M.; Licht, T.R. Microbial tryptophan catabolites in health and disease. Nat. Commun. 2018, 9, 3294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Krautkramer, K.A.; Fan, J.; Bäckhed, F. Gut microbial metabolites as multi-kingdom intermediates. Nat. Rev. Microbiol. 2021, 19, 77–94. [Google Scholar] [CrossRef] [PubMed]
  35. Nielsen, J. Systems Biology of Metabolism. Annu. Rev. Biochem. 2017, 86, 245–275. [Google Scholar] [CrossRef]
  36. Magnúsdóttir, S.; Heinken, A.; Kutt, L.; Ravcheev, D.A.; Bauer, E.; Noronha, A.; Greenhalgh, K.; Jäger, C.; Baginska, J.; Wilmes, P.; et al. Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat. Biotechnol. 2017, 35, 81–89. [Google Scholar] [CrossRef] [PubMed]
  37. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI Reference Sequence (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, 33, D501–D504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Machado, D.; Andrejev, S.; Tramontano, M.; Patil, K.R. Fast automated reconstruction of genome-scale metabolic models for microbial species and communities. Nucleic Acids Res. 2018, 46, 7542–7553. [Google Scholar] [CrossRef]
  39. Mendoza, S.N.; Olivier, B.G.; Molenaar, D.; Teusink, B. A systematic assessment of current genome-scale metabolic reconstruction tools. Genome Biol. 2019, 20, 158. [Google Scholar] [CrossRef] [Green Version]
  40. Forster, S.C.; Kumar, N.; Anonye, B.O.; Almeida, A.; Viciani, E.; Stares, M.D.; Dunn, M.; Mkandawire, T.T.; Zhu, A.; Shao, Y.; et al. A human gut bacterial genome and culture collection for improved metagenomic analyses. Nat. Biotechnol. 2019, 37, 186–192. [Google Scholar] [CrossRef]
  41. Zou, Y.; Xue, W.; Luo, G.; Deng, Z.; Qin, P.; Guo, R.; Sun, H.; Xia, Y.; Liang, S.; Dai, Y.; et al. 1520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat. Biotechnol. 2019, 37, 179–185. [Google Scholar] [CrossRef] [PubMed]
  42. Arumugam, M.; Raes, J.; Pelletier, E.; Le Paslier, D.; Yamada, T.; Mende, D.R.; Fernandes, G.R.; Tap, J.; Bruls, T.; Batto, J.-M.; et al. Enterotypes of the human gut microbiome. Nature 2011, 473, 174–180. [Google Scholar] [CrossRef]
  43. Lagkouvardos, I.; Pukall, R.; Abt, B.; Foesel, B.U.; Meier-Kolthoff, J.P.; Kumar, N.; Bresciani, A.; Martínez, I.; Just, S.; Ziegler, C.; et al. The Mouse Intestinal Bacterial Collection (miBC) provides host-specific insight into cultured diversity and functional potential of the gut microbiota. Nat. Microbiol. 2016, 1, 16131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Liu, C.; Zhou, N.; Du, M.-X.; Sun, Y.-T.; Wang, K.; Wang, Y.-J.; Li, D.-H.; Yu, H.-Y.; Song, Y.; Bai, B.-B.; et al. The Mouse Gut Microbial Biobank expands the coverage of cultured bacteria. Nat. Commun. 2020, 11, 79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wilmes, P.; Bowen, B.P.; Thomas, B.C.; Mueller, R.S.; Denef, V.J.; VerBerkmoes, N.C.; Hettich, R.L.; Northen, T.R.; Banfield, J.F. Metabolome-proteome differentiation coupled to microbial divergence. MBio 2010, 1, e00246-10. [Google Scholar] [CrossRef] [Green Version]
  46. Hubbard, T.D.; Murray, I.A.; Bisson, W.H.; Lahoti, T.S.; Gowda, K.; Amin, S.G.; Patterson, A.D.; Perdew, G.H. Adaptation of the human aryl hydrocarbon receptor to sense microbiota-derived indoles. Sci. Rep. 2015, 5, 12689. [Google Scholar] [CrossRef]
  47. Ramadoss, P.; Perdew, G.H. Use of 2-azido-3-[125I]iodo-7,8-dibromodibenzo-p-dioxin as a probe to determine the relative ligand affinity of human versus mouse aryl hydrocarbon receptor in cultured cells. Mol. Pharmacol. 2004, 66, 129–136. [Google Scholar] [CrossRef]
  48. Dong, M.; Li, L.; Chen, M.; Kusalik, A.; Xu, W. Predictive analysis methods for human microbiome data with application to Parkinson’s disease. PLoS ONE 2020, 15, e0237779. [Google Scholar] [CrossRef]
  49. Marcos-Zambrano, L.J.; Karaduzovic-Hadziabdic, K.; Loncar Turukalo, T.; Przymus, P.; Trajkovik, V.; Aasmets, O.; Berland, M.; Gruca, A.; Hasic, J.; Hron, K.; et al. Applications of Machine Learning in Human Microbiome Studies: A Review on Feature Selection, Biomarker Identification, Disease Prediction and Treatment. Front. Microbiol. 2021, 12, 634511. [Google Scholar] [CrossRef]
  50. Xia, F.; Chen, J.; Fung, W.K.; Li, H. A logistic normal multinomial regression model for microbiome compositional data analysis. Biometrics 2013, 69, 1053–1063. [Google Scholar] [CrossRef]
  51. Mangge, H.; Stelzer, I.; Reininghaus, E.Z.; Weghuber, D.; Postolache, T.T.; Fuchs, D. Disturbed tryptophan metabolism in cardiovascular disease. Curr. Med. Chem. 2014, 21, 1931–1937. [Google Scholar] [CrossRef]
  52. Millett, E.R.C.; Peters, S.A.E.; Woodward, M. Sex differences in risk factors for myocardial infarction: Cohort study of UK Biobank participants. BMJ 2018, 363, k4247. [Google Scholar] [CrossRef] [Green Version]
  53. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  54. Smith, C.A.; Want, E.J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 2006, 78, 779–787. [Google Scholar] [CrossRef]
  55. Shen, X.; Wu, S.; Liang, L.; Chen, S.; Contrepois, K.; Zhu, Z.J.; Snyder, M. metID: A R package for automatable compound annotation for LC-MS-based data. Bioinformatics 2021. [Google Scholar] [CrossRef]
  56. Chong, J.; Liu, P.; Zhou, G.; Xia, J. Using MicrobiomeAnalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat. Protoc. 2020, 15, 799–821. [Google Scholar] [CrossRef] [PubMed]
  57. Kerley-Hamilton, J.S.; Trask, H.W.; Ridley, C.J.; Dufour, E.; Lesseur, C.; Ringelberg, C.S.; Moodie, K.L.; Shipman, S.L.; Korc, M.; Gui, J.; et al. Inherent and benzo[a]pyrene-induced differential aryl hydrocarbon receptor signaling greatly affects life span, atherosclerosis, cardiac gene expression, and body and heart growth in mice. Toxicol. Sci. 2012, 126, 391–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Huang, S.; Shui, X.; He, Y.; Xue, Y.; Li, J.; Li, G.; Lei, W.; Chen, C. AhR expression and polymorphisms are associated with risk of coronary arterial disease in Chinese population. Sci. Rep. 2015, 5, 8022. [Google Scholar] [CrossRef] [Green Version]
  59. Wu, D.; Nishimura, N.; Kuo, V.; Fiehn, O.; Shahbaz, S.; Van Winkle, L.; Matsumura, F.; Vogel, C.F. Activation of aryl hydrocarbon receptor induces vascular inflammation and promotes atherosclerosis in apolipoprotein E-/- mice. Arterioscler. Thromb. Vasc. Biol. 2011, 31, 1260–1267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Scott, S.A.; Fu, J.; Chang, P.V. Microbial tryptophan metabolites regulate gut barrier function via the aryl hydrocarbon receptor. Proc. Natl. Acad. Sci. USA 2020, 117, 19376–19387. [Google Scholar] [CrossRef] [PubMed]
  61. Lloyd-Price, J.; Arze, C.; Ananthakrishnan, A.N.; Schirmer, M.; Avila-Pacheco, J.; Poon, T.W.; Andrews, E.; Ajami, N.J.; Bonham, K.S.; Brislawn, C.J.; et al. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases. Nature 2019, 569, 655–662. [Google Scholar] [CrossRef]
  62. Manzella, C.R.; Jayawardena, D.; Pagani, W.; Li, Y.; Alrefai, W.A.; Bauer, J.; Jung, B.; Weber, C.R.; Gill, R.K. Serum Serotonin Differentiates Between Disease Activity States in Crohn’s Patients. Inflamm. Bowel Dis. 2020, 26, 1607–1618. [Google Scholar] [CrossRef] [PubMed]
  63. Patel, B.B.; Kasneci, A.; Bolt, A.M.; Di Lalla, V.; Di Iorio, M.R.; Raad, M.; Mann, K.K.; Chalifour, L.E. Chronic Exposure to Bisphenol a Reduces Successful Cardiac Remodeling After an Experimental Myocardial Infarction in Male C57bl/6n Mice. Toxicol. Sci. 2015, 146, 101–115. [Google Scholar] [CrossRef] [Green Version]
  64. Kasneci, A.; Lee, J.S.; Yun, T.J.; Shang, J.; Lampen, S.; Gomolin, T.; Cheong, C.C.; Chalifour, L.E. From the Cover: Lifelong Exposure of C57bl/6n Male Mice to Bisphenol A or Bisphenol S Reduces Recovery from a Myocardial Infarction. Toxicol. Sci. 2017, 159, 189–202. [Google Scholar] [CrossRef]
  65. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.E.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef] [PubMed]
  66. Integrative, H.M.P.R.N.C. The Integrative Human Microbiome Project: Dynamic analysis of microbiome-host omics profiles during periods of human health and disease. Cell Host Microbe 2014, 16, 276–289. [Google Scholar] [CrossRef] [Green Version]
  67. Vehtari, A.; Gelman, A.; Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat Comput. 2017, 27, 1413–1432. [Google Scholar] [CrossRef] [Green Version]
Metabolites 12 00010 g001 550
Figure 1. The development of TrpNet.
Figure 1. The development of TrpNet.
Metabolites 12 00010 g001
Metabolites 12 00010 g002 550
Figure 2. Host and microbial tryptophan metabolism pathways.
Figure 2. Host and microbial tryptophan metabolism pathways.
Metabolites 12 00010 g002
Metabolites 12 00010 g003 550
Figure 3. Distribution of tryptophan metabolite production across the human gut microbiome. Each branch indicates a microbe strain colored on the basis of their phyla (3HAA: 3-hydroxyanthranilate; FKYN: formylkynurenine; HKYN: 3-hydroxykynurenine; IA: 3-indoleacrylate; IAA: indole-3-acetate; IAAlD: indole-3-acetaldehyde; Iald: indole-3-carboxaldehyde; IAM: indole-3-acetamide; IE: indole-3-ethanol; ILA: indolelactate; IPA: indolepropionate; IPY: indolepyruvate; KYN: l-kynurenine; QA: quinolinate).
Figure 3. Distribution of tryptophan metabolite production across the human gut microbiome. Each branch indicates a microbe strain colored on the basis of their phyla (3HAA: 3-hydroxyanthranilate; FKYN: formylkynurenine; HKYN: 3-hydroxykynurenine; IA: 3-indoleacrylate; IAA: indole-3-acetate; IAAlD: indole-3-acetaldehyde; Iald: indole-3-carboxaldehyde; IAM: indole-3-acetamide; IE: indole-3-ethanol; ILA: indolelactate; IPA: indolepropionate; IPY: indolepyruvate; KYN: l-kynurenine; QA: quinolinate).
Metabolites 12 00010 g003
Metabolites 12 00010 g004 550
Figure 4. Tanglegram between the dendrograms of phylogenetic and metabolic distance. The phylogenetic dendrogram generated by hierarchical clustering with complete linkage of the taxonomy rank of maximum likelihood tree. The dendrogram of metabolic distance was calculated on the basis of the presence or absence of tryptophan metabolite production. Lines are colored by genus and connect the same microbes.
Figure 4. Tanglegram between the dendrograms of phylogenetic and metabolic distance. The phylogenetic dendrogram generated by hierarchical clustering with complete linkage of the taxonomy rank of maximum likelihood tree. The dendrogram of metabolic distance was calculated on the basis of the presence or absence of tryptophan metabolite production. Lines are colored by genus and connect the same microbes.
Metabolites 12 00010 g004
Metabolites 12 00010 g005 550
Figure 5. ROC plot for top-ranked IAA models based on prediction of stimulated data. The curves are colored by different models which use different predictors as listed for predicting IAA production.
Figure 5. ROC plot for top-ranked IAA models based on prediction of stimulated data. The curves are colored by different models which use different predictors as listed for predicting IAA production.
Metabolites 12 00010 g005
Metabolites 12 00010 g006 550
Figure 6. A screenshot of TrpNet showing the overall tryptophan metabolic network.
Figure 6. A screenshot of TrpNet showing the overall tryptophan metabolic network.
Metabolites 12 00010 g006
Metabolites 12 00010 g007 550
Figure 7. (a) Tryptophan metabolite prediction from microbiome data; (b) TrpNet prediction compared with metabolomics data and PICRUSt prediction.
Figure 7. (a) Tryptophan metabolite prediction from microbiome data; (b) TrpNet prediction compared with metabolomics data and PICRUSt prediction.
Metabolites 12 00010 g007
Metabolites 12 00010 g008 550
Figure 8. (a) Tryptophan metabolite prediction from microbiome data; (b) TrpNet prediction compared with metabolomics data and PICRUSt2 prediction.
Figure 8. (a) Tryptophan metabolite prediction from microbiome data; (b) TrpNet prediction compared with metabolomics data and PICRUSt2 prediction.
Metabolites 12 00010 g008
Table
Table 1. Odds ratio of dominate genera in mouse gut for bioactive indole generation (red: p-value < 0.001; orange: p-value < 0.01, yellow: p-value < 0.05, blue p-value < 0.1 based on Wald test).
Table 1. Odds ratio of dominate genera in mouse gut for bioactive indole generation (red: p-value < 0.001; orange: p-value < 0.01, yellow: p-value < 0.05, blue p-value < 0.1 based on Wald test).
PredictorsIAIndoleIAAlDIAMIAAILAIPATryptamine
Bacteroides0.8786310.41181.52560.56212.942469.00480.85150.1484
Bifidobacterium0.87120.04210.48790.60811.0597103.14760.83938.5582
Clostridium413.06812.25261.6328106.63080.840189.0063638.31643.473
Desulfovibrio0.97381.52261.24510.813914.92150.89310.96760.3856
Enterococcus0.92252.16670.88610.70170.03920.64460.90171.0872
Escherichia0.936241.12310.0870.79973.9424226.0360.91619.0402
Eubacterium0.97833.71210.6220.861516.72530.87920.97230.4817
Lactobacillus0.85360.03241.97940.50871.76536.94240.82270.4338
Mouse.gut2.59422.29311.14240.9371.83196.92112.91190.4471
Parabacteroides0.96680.18410.93440.83841.260513.82910.95690.4585
Prevotella0.91451.68550.70290.59470.4730.72860.89690.1589
Ruminococcus0.97180.74490.40290.80357.28950.88610.96510.3699
Streptococcus0.89590.62880.72870.55330.5270.67560.874919.0245
Table
Table 2. Odds ratio of dominate genera in human gut for bioactive indole generation (red: p-value < 0.001; orange: p-value < 0.01, yellow: p-value < 0.05, blue p-value < 0.1 based on Wald test).
Table 2. Odds ratio of dominate genera in human gut for bioactive indole generation (red: p-value < 0.001; orange: p-value < 0.01, yellow: p-value < 0.05, blue p-value < 0.1 based on Wald test).
PredictorsIAIndoleIAAlDIAMIAAILAIPATryptamine
Bacteroides0.88551595.58321.45950.52652.648979.56180.85750.1214
Bifidobacterium0.90250.03710.76580.56740.7158175.30570.87815.5164
Clostridium414.02541.86061.236691.49661.626881.5643606.56032.8663
Desulfovibrio0.96831.550.60040.7947.25920.85120.95930.3552
Enterococcus0.94782.5611.26730.70370.03220.78580.93330.9879
Escherichia0.9639318.8560.081.0.76774.0763607.72440.95314.6206
Eubacterium0.9812.02081.24660.855910.53650.90350.97530.463
Human.gut21.32041.44130.7778342.74062.168520.485337.3321876.3277
Lactobacillus0.87570.12562.44920.50551.534352.76020.84620.412
Parabacteroides0.9760.19641.30150.82711.68128.29690.96870.4121
Prevotella0.94793.52732.00260.70351.2940.78790.93320.2527
Ruminococcus0.96630.60210.34770.77844.98190.84870.95620.3397
Streptococcus0.88710.4510.74350.53040.34670.63570.859618.3925
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, Y.; Chong, J.; Shen, S.; Chammas, J.-B.; Chalifour, L.; Xia, J. TrpNet: Understanding Tryptophan Metabolism across Gut Microbiome. Metabolites 2022, 12, 10. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12010010

AMA Style

Lu Y, Chong J, Shen S, Chammas J-B, Chalifour L, Xia J. TrpNet: Understanding Tryptophan Metabolism across Gut Microbiome. Metabolites. 2022; 12(1):10. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12010010

Chicago/Turabian Style

Lu, Yao, Jasmine Chong, Shiqian Shen, Joey-Bahige Chammas, Lorraine Chalifour, and Jianguo Xia. 2022. "TrpNet: Understanding Tryptophan Metabolism across Gut Microbiome" Metabolites 12, no. 1: 10. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12010010

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop