Metagenomic analysis of bile salt biotransformation in the human gut microbiome

doi:10.21203/rs.2.9769/v1

Download PDF

Research article

Metagenomic analysis of bile salt biotransformation in the human gut microbiome

https://doi.org/10.21203/rs.2.9769/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 29 Jun, 2019

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

Background: In the biochemical milieu of human colon, bile acids act as signaling mediators between the host and its gut microbiota. Biotransformation of primary to secondary bile acids have been known to be involved in the immune regulation of human physiology. Several 16S amplicon-based studies with inflammatory bowel disease (IBD) subjects were found to have an association with the level of fecal bile acids. However, a detailed investigation of all the bile acid biotransformation genes involved in the secondary bile acid metabolism has not been performed. Results: Here, we report a comprehensive analysis of the bile acid biotransformation genes and their distribution at the phyla level. Based on the analysis of shotgun metagenomes, we found that the IBD subjects harbored a significantly lower abundance of these genes compared to the healthy controls. Majority of these genes originated from Firmicutes in comparison to other phyla. From metabolomics data, we found that the IBD subjects were measured with a significantly low level of secondary bile acids and high levels of primary bile acids compared to that of the healthy controls. Conclusions: Our bioinformatics-driven approach of identifying bile acid biotransformation genes predicts the bile salt biotransformation potential in the gut microbiota of IBD subjects. The functional level of dysbiosis likely contributes to the variation in the bile acid pool. This study sets the stage to envisage potential solutions to modulate the gut microbiome with the objective to restore the bile acid pool in the gut.

Epigenetics & Genomics

Secondary bile acids

gut microbiota

IBD

bioinformatics

metagenomics

metabolomics

and bile salt biotransformation genes

Microbial modification of human-derived bile acids (BAs) is known as secondary BAs, which provide cross-talk between the human host and gut bacteria in the large intestine [1]. Primary BAs such as cholic acid (CA) and chenodeoxycholic acid (CDCA) that are produced in the liver are modified to deoxycholic acid (DCA) and lithocholic acid (LCA) in the colon respectively. These secondary BAs have been known to regulate several host metabolic processes via receptor signaling, including the farnesoid X receptor (FXR), the liver X receptor (LXR), the G-protein coupled receptor TGR5, and the vitamin D receptor present in the gut, the liver and in the periphery [2]. Concerning the anti-inflammatory and anti-microbial functions exerted by the secondary BAs, dysbiosis in the bile acid biotransformation potential of gut microbiota might influence the function of BAs in the host physiology.

Previous studies of small-sized cohorts or 16S amplicon based studies have identified significant differences in the fecal metabolites between the healthy and Inflammatory Bowel Disease (IBD) subjects [3–6]. From a growing body of evidence, two well-established IBD-associated taxonomic signatures have been observed: (i) phylum-level decrease in Firmicutes (especially of class Clostridia and family Lachnospiraceae), and (ii) phylum-level increase in Proteobacteria (especially of family Enterobacteriaceae) [7]. At metabolite level, mouse-models of inflammation were measured with a significant reduction of gut microbiota-derived bile acids [8–10]. This decrease has resulted in an altered pathophysiological response in the host via an altered-composition of the gut microbiome.

However, the extent of secondary BA metabolism by the gut microbiota remains less studied as the bacterial species that are specialized with bile acid metabolic genes have been characterized and studied in only few Clostridium species [11–14]. Furthermore, so far only a few of the enzymatic genes comprising the entire pathway of BA metabolism had been assessed from shotgun metagenomic sequences [15, 16]. None of the previous studies provided a complete picture considering all the relevant enzymatic genes, on the secondary bile acid metabolism, i.e., their prevalence in other bacterial species and abundances under different health states.

Henceforth, we sought to investigate secondary BA pathway-related protein and gene sequences based on sequence and structural domain conservation. Once they have been screened, estimation of these sequences from the shotgun metagenomes would allow predicting the bile acid metabolic potential in a microbial community. Based on this approach, the following research questions were formulated: (i) What is the distribution profile of bile salt biotransformation protein (BSBP) homologs at phyla level? (ii) How are these BSBP homologs prevalent at strain level? (iii) How conserved or divergent these homologous protein sequences are with respect to the phylum and enzymatic category? (iv) Does the abundance of bile salt biotransformation genes (BSBGs) vary among the fecal metagenomes of similar group and between different groups? and (v) Could parallel observations be deduced at metabolite level based on the findings from gene abundance estimation?

To provide a comprehensive view on gut microbial-mediated bile acid metabolism, the following analyses were performed: (i) Identification of BSBP homologs in microbial strains and their distribution at phyla level; (ii) Analysis of protein-sequence similarity network (PSSN) of identified homologs for each enzymatic function; (iii) Quantitative comparison of BSBG abundances between healthy controls and IBD subjects; (iv) Phylum-level quantitative comparison of BSBG abundances between healthy controls and IBD subjects, and (v) Quantitative comparison of fecal bile acid metabolites between healthy controls and IBD subjects.

Identification and distribution of bile salt biotransformation protein homologs

To identify BSBP homologs, BSBP sequences from Clostridium scindens were deployed as reference query sequences (Additional file 1: Table S1). The reasons for considering this species as reference are: (i) limited experimental data (in testing the bacterial growth sensitivity to BAs) or lack in availability of BSBP or BSBG reference dataset (like an ARG database for antibiotic-resistance genes) [17], and (ii) extensive characterization of enzymatic genes comprising the secondary bile acid metabolic pathway (as shown to be encoded by bai operon) in the species C. scindens [13, 14, 18]. Through sensitivity-specificity analysis and user-defined methods of criteria to select the optimal parameters (Additional file 1: Table S2 and S4), it was ensured that the selected sequences were true positive homologs and rejected sequences that were false positive and true negative homologs. Therefore, based on protein sequence and domain conservation, a total of 10,613 protein homologs were identified. These counts of protein sequences were distributed based on the phyla and enzyme-homology (Table 1).

Based on the taxonomic origin (i.e., at strain level) of protein sequences, all the identified protein homologs were inspected for the presence and absence of each functional enzyme category as presented in Additional file 2: Table S1. Microbial strains that were prevalent for more than half of the enzyme categories were shown in Fig 1. From the presence and absence heatmap (Fig 1), we observed that species from the same genus had a similar profile pattern of BSBP homologs. Moreover, in a majority of the species, the distribution pattern of these enzymatic sequences was not as similar to that in the reference bai operon from C. scindens (Additional file 2: Table S1). Microbial strains of Lactobacillus species were known to deconjugate bile salts [19, 20], and we observed the same with prevalence for BSH homologs. However, enteric pathogens like Shigella, Salmonella, and Klebsiella were prevalent for Bai-like proteins but not for BSH homolog. (Additional file 2: Table S1).

Protein-sequence similarity network of bile salt biotransformation protein homologs

Next, to evaluate the genetic diversity of identified homologs at the phylum level, protein sequences from each enzyme category were analyzed using sequence-based similarity network (Fig 2). For a multiple sequence alignment (for at least thousands of sequences), the low resolution of nodes is often encountered due to a low number of informative phylogenetic sites [21]. Hence, to overcome this challenge, we performed the following analysis. In a PSSN, the number of cluster formation would differ based on the threshold value, such as there would be fewer clusters for relaxed threshold value and more clusters for stringent threshold value. The optimal threshold value for each network shown in Fig 2 is provided in Additional file 1: Table S4. From a broad perspective, we observed that sequences from similar phylum were co-localized on the network, with minor exceptions where sequences from different phylum were clustered together.

However, on a detailed level of inspection, BaiCD/H showed a single cluster even at higher stringent edge thresholds in comparison to the optimal values of other networks, which denotes less divergence within these homologous sequences across diverse phylum. While for the other networks, there was a variation in the number of clusters from a range of relaxed to stringent threshold criteria (Fig 2 and Additional file 1: Table S4). For instance, a network for BaiG and BaiA/L (enzyme category with few protein connections or more clusters) identifies farther sequence divergence among the putative homologs. For each enzyme category, even though these protein homologs were identified through sequence conservation to a reference sequence, the clustering pattern reveals the biological degree of divergence or convergence among the identified homologs at their phyla-level.

Comparative metagenomic analysis of bile salt biotransformation gene abundances between healthy and IBD subjects

Through investigation of gene abundances (for the identified list of protein homologs, as described in the Methods section) in the fecal metagenomes, bile salt biotransformation potential was compared quantitatively in the gut microbiota of healthy and IBD subjects.

In the METAHIT-Spanish cohort, on comparing the normalized gene abundance of total BSBGs between healthy (n = 14) and IBD (n = 25) subjects, we found a significant difference between the abundance values of healthy and IBD subjects. The results in Fig 3a show that the mean of the normalized abundance of BSBGs was lower in IBD subjects (i.e., 7.77e-05) than that of healthy controls (i.e., 1.08e-04) (Mann–Whitney–Wilcoxon test, p < 0.05), suggesting lower potency of BA biotransformation in the IBD subjects than that of the healthy controls. Our results here were consistent with findings from similar analyses performed in Labbe et al. [22], where samples from the MetaHIT cohort were also analyzed. However, further statistical comparisons between the subtypes of IBD subjects (i.e., Ulcerative colitis (UC) and Crohn's disease (CD)) could not be performed due to an insufficient number of CD samples.

In the iHMP-American cohort, on comparing the normalized gene abundance of total BSBGs in the healthy (n = 18) and IBD (n = 65) subjects, we could not find any significant difference between the abundance values of healthy and IBD subjects. However, on all possible combination of statistical comparisons between the healthy, UC and CD subjects, the results in Fig 3b show that the mean of normalized abundance of BSBGs was lower in the CD subjects (i.e., 3.68e-05) than that of the healthy controls (i.e., 4.34e-05) (Mann–Whitney–Wilcoxon test, p < 0.05), suggesting lower potency of BA biotransformation in the CD subjects than that of the healthy controls. No significant differences were observed between the following comparisons (i) healthy and UC subjects, and (ii) UC and CD subjects.

Comparative metagenomic analysis of bile salt biotransformation gene abundances among healthy subjects from different countries

Similar to above, through investigation of gene abundances (for the identified list of protein homologs) in the fecal metagenomes, we compared the bile acid biotransformation potential in the gut microbiota of healthy subjects from three different countries, i.e., USA (n = 18), Denmark (n = 86) and Spain (n = 14). Additional file 1: Figure S2 revealed variation in the normalized abundances of the total BSBGs. From all possible combination of pairwise comparisons among different countries, we observed that the significant differences among healthy individuals of different countries. While individuals from the USA showed a significant difference from that of Denmark and Spain (Mann Whitney Wilcoxon, p < 0.01), no significant difference was found between the healthy individuals of Denmark and Spain. These findings hint at the possibility of diet as an important factor that could influence the composition of the gut microbiota, thereby leading to the differences in the abundance of BSBGs [23–25]. This analysis also highlights the importance of considering demographically matched samples as a prerequisite for comparative functional studies.

Phyla-level comparative metagenomic analysis of bile salt biotransformation gene abundances between healthy and IBD subjects

To evaluate the contribution of phyla accounting for BSBG abundances, each mapped gene was retrieved for its corresponding phylum, and the normalized abundances of phyla-level BSBGs were computed accordingly for each group in both the cohort (Fig 4). The results show that the genes from Firmicutes phylum accounted for the most considerable abundance for the BSBG in comparison to other phyla (Kruskal-Wallis test, p < 0.01). This observation, i.e., high abundance of genes from the Firmicutes phylum, is consistent with the existing literature, which states that species from Firmicutes phylum dominate in BA metabolism [16]. However, lower abundance of Firmicutes phylum-specific BSBGs in IBD subjects compared to that of healthy controls further confirms one of the typical IBD signature (i.e., marked by a decreased abundance of species from Firmicutes) (Mann-Whitney-Wilcoxon test, p < 0.05) [26–28]. Furthermore, no significant differences were observed in the level of BSBGs for the phyla: Actinobacteria, Bacteroidetes and Proteobacteria between the healthy and IBD subjects.

Comparative fecal metabolomics analysis of bile acids between healthy and IBD subjects

To verify if the IBD individuals have reduced levels of secondary BAs corresponding to their predicted low abundance of BSBGs than that of healthy controls, metabolomics data from the iHMP cohort were analyzed. Fecal bile acid metabolites were measured in LC-MS (C18 negative ion mode analysis) [29, 30]. To compare the bile acid profile between healthy and IBD subjects, quantification of BAs was expressed in proportion, where the level of each BA was calculated to the total level of BAs. Fig 5 shows that the proportion of conjugated primary BAs (i.e., cholate and chenodeoxycholate) and secondary BAs (i.e., deoxycholate and lithocholate) were found to be at increased and decreased levels in IBD subjects compared to that of healthy controls respectively (Mann-Whitney-Wilcoxon test, p < 0.01). The decrease in the level of secondary BAs compared to the increased level of primary BAs suggests a decrease in the potency of BA biotransformation, which is consistent with our comparative metagenomic analysis. These results were in line with a similar analysis performed on this cohort [29, 31]. Furthermore, the proportions of other different forms of conjugated BA metabolites were found at higher levels in IBD subjects than that of healthy controls (Additional file 1: Figure S3). These observations were congruent to that found in Duboc et al. [32] where low levels of secondary BAs were found in IBD subjects, who were in the active phase and clinical remission state.

In this study, we performed bioinformatic analysis on secondary bile acid metabolism and explored their abundances in the whole-genome metagenomic sequences of healthy and IBD subjects. The overall aim of this work was to gain an understanding of the bile salt biotransformation potential of the human gut microbiota. This provided insights into how the change in the abundance of BSBGs altered the level of BA metabolic profiles in a dysbiotic state.

Based on the analyses that have been performed, we highlight the following results resolving back to our formulated objectives (i) BSBP homologs were prevalent across abundant phyla, such as Bacteroidetes, Firmicutes, Actinobacteria, and Proteobacteria. However, phyla-level differential abundant analysis revealed the dominant role of Firmicutes for BSBGs; (ii) The prevalence of BSBP profile varied from a low to a high scale of spectrum; (iii) PSSN revealed the degree of divergence for different enzymatic categories; (iv) BSBG-specific quantitative estimation of shotgun metagenomes suggests lower potency of bile salt biotransformation in the gut microbiota of IBD subjects; (v) Evaluation at metabolomics level revealed parallel reduction of gut microbial-derived bile acids in IBD subjects compared to healthy controls.

Homology-based detection of BSBP homologs suggests three patterns as discussed: firstly, the presence of a partial set of homologs in comparison to the whole set of enzymatic categories might suggest a similar situation as found in C. hiradonis and C. sordellii. For instance, C. hiradonis carried baiBCDEA2FGHJ while C. sordellii carried only baiCDA2HE [33]. Alternatively, it could be a possibility of similar function with similar substrates to that of bile acid metabolism. Secondly, experimental approaches to elucidate the requirement for minimum bai gene set could clear the specific distribution of these genes across gut species. Thirdly, it is to be established which combination of Bai and BSH protein contribute to the phenotype of resistance or sensitivity to BAs, because, only a subset of species was found to carry BSH homologs, while other species lacked in BSH homolog. Extracellular enzymatic activities similar to BSH could perhaps provide functional compensation for species lacking the enzyme [34, 35].

Besides evaluation of BSBGs in IBD subjects, similar dysbiosis has been observed in cirrhotic patients, where a low-level input of primary BAs to the gut leads to a reduction of BA metabolizers in the large intestine. Accordingly, reduced levels of secondary BAs in the feces and increased levels of primary BAs were found in the serum and feces of cirrhotic patients [36–38]. As both the cases converge to a similar metabolic phenotype, further investigation at advanced level could avoid misclassification of these subjects. Contrarily, high physiological levels of secondary BA have been associated with colorectal cancer patients [39]. Hence, any deviation from the physiological range of secondary BAs could contribute to a disease phenotype.

The biological relevance of these secondary and tertiary BAs, such as DCA, LCA, ursodeoxycholic acid (UCDA) have known to be potent agonists for several BA receptors and exert anti-inflammatory effects in the colon. For instance, they bind to TGR5 and inhibit the production of cytokines once encountered with lipopolysaccharides from gram-negative bacteria [40]. LCA could also regulate adaptive immune responses via VDR signaling [41] whereas FXR activation maintains the production of BAs and their homeostasis [42]. These properties of secondary BAs suggest their beneficial role in reducing inflammation in IBD [43].

Considering the relation between BA metabolism and IBD, one potential application using probiotics has alleviated disease remission of relapsing UC patients [44, 45]. To tolerate bile acid stress in the gut, the ability to hydrolyze bile salts is often considered for probiotic strain selection [20, 46]. Therefore, to improve the gut health of IBD subjects, had they been supplemented with BSH-active strain formulated in probiotics along with their routine anti-inflammatory drug, synergistic effects of this combination could help in the conversion of bile salts to secondary BAs. This approach has improved the clinical features of relapsing UC patients suggesting that interventions targeting BA metabolism may have therapeutic implications [44, 45].

Taken together, our approach of targeted identification and quantification of bile salt biotransformation has enabled us to realize the impact of gut microbial mediated metabolic effects on the host system. This functional dysbiosis was observed in the gut microbiota of IBD individuals, who were in early-onset IBD or clinical remission stage. We also observed that the reduction of Firmicutes-specific bile acid-related genes is in line with the compositional reduction of Firmicutes (especially species with bile metabolic potential). From a broad view of the omics lens, the metagenomic predictions and the metabolomics evaluations were consistent in both approaches. However, verification of identified homologs from few of the cultivable species, through a reductionist approach, could further verify the biological significance of sequence conservation. Nonetheless, gut microbial composition and their function could potentially tip a compensated state of immune balance in favor of chronic disease in IBD hosts.

IBD cohort description

Shotgun metagenomes were obtained from European Nucleotide Archive at EMBL-EBI under the accession number PRJEB2054 and PRJNA389280.

From PRJEB2054, 14, 4, and 21 samples of Healthy (H), Crohn’s disease (CD), and Ulcerative colitis (UC) respectively were analyzed for differential BSBG abundance analysis [47]. Both subtypes of IBD subjects were in clinical remission for at least 3 months and had stable maintenance therapy with mesalazine or azathioprine. Individuals who were undergoing antibiotic treatment for at least 4 weeks before fecal sample collection were excluded [48].

From PRJNA389280, 18, 44, and 21 samples of Healthy (H), Crohn's disease (CD) and Ulcerative colitis (UC) respectively were analyzed for differential BSBG abundance and BA metabolite analysis. These samples had paired or coordinated metagenomics and metabolomics data available. None of these IBD subjects had a history of terminal ileal resection. Detailed information for the sample selection criteria on IBD subjects is available in the original article [29].

Bioinformatic identification of bile salt biotransformation protein homologs

To identify BSBP homologs from the UniProt database [49], experimentally verified BSBPs from Clostridium scindens (strain JCM 10418 / VPI 12708) and other species were deployed as query sequences (Additional file 1: Table S1 and Figure S1). The selection process of BSBP homologs was performed in two sections. First, screening the candidate BSBPs using BLASTp [50] under optimal threshold parameters (Additional file 1: Table S2). To determine the threshold for BLASTp search, protein sequences with known BSB function and non-BSB function were used as test data. For instance, BSB proteins from Eggerthella lenta was used as positive control [51, 52], while protein sequences from Helicobacter, Prevotella, and Porphyromonas were used as negative control [53–58]. The accuracy of BLASTp was estimated at different similarity thresholds with defined e-value and coverage (Additional file 1: Table S2).

Using the BSH-specific HMM and PVA-specific HMM [59, 60], protein sequences of bile salt hydrolase (BSH) category were identified to make a distinct classification from penicillin-v-acylase (PVA) enzyme sequences, as both shared identical domain CBAH superfamily. Sequences that had higher e-value in the BSH cluster compared to that in the PVA cluster were selected for further analysis.

Secondly, the candidate BSBP homologs were searched against the Pfam 31.0 database to identify their functional domains using the default GA (gathering cut-offs) [61]. Sequences with identical domain organizations were retained for the following analysis (Additional file 1: Table S3). Additionally, homologous protein sequences that were filtered across the Pfam domain analysis were checked for their functional annotation using the eggNOG algorithm [62] with DIAMOND as mapping mode and other default parameters. Sequences that mapped to similar annotations to the reference protein were considered for the final reference database. The taxonomic lineage of the identified gene sequences was retrieved from the UniProt website (dated on October 2018) [49].

Construction of protein sequence similarity network (PSSN)

An all-by-all protein BLAST was performed for all the homologous sequences to be analyzed in the dataset. Pairs of sequences with significant similarity were selected through the application of a series of pairwise protein BLAST e-value thresholds, where their distribution was used to define the optimal e-value cut-off. Subsequently, significant alignments (with e-value less than the optimal e-value and minimum of 30% pairwise similarity) were used to construct the PSSN (Additional file 1: Table S4). Each node of the network indicates a protein sequence, and the edge represents the relatedness between two nodes sharing significant similarity as defined by the edge metric (i.e., e-value and percentage identity). Networks were visualized with the ggnetwork package in R [63, 64].

Functional metagenomic mapping

For quality assessment, publicly available fastq reads were subjected to a quality check using FastQC [65] with a minimum Phred score of 30. Each pair of fastq files was repaired to remove the singletons (if needed) using the script available at (https://github.com/BioInfoTools/BBMap/blob/master/sh/repair.sh). Bile salt biotransformation gene abundances were calculated using FMAP [66] with the default parameters. Gene sequences that were binned into the functional category of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway ko00121 were determined by mapping of UniRef to KEGG Orthologies (KOs) database [67, 68]. Then, the numbers of mapped reads were normalized to the total number of paired reads in the metagenomes.

Statistical analysis

All the statistical analysis and graphics were performed using R version 3.5.0 [69]. In each cohort, samples were country-matched to avoid cross country-effects. Normality of the given dataset was assessed through the Shapiro–Wilk test. Differential abundance analysis between healthy and IBD subjects was tested using Mann-Whitney-Wilcoxon tests (for comparison between two groups), and Kruskal Wallis tests (for comparison among three groups). Multivariable linear regression was performed for age, BMI and gender on the gene abundance to identify confounding variable (p < 0.01) by unadjusted analysis.

IBD: Inflammatory Bowel Disease

CD: Crohn’ Disease

UC: Ulcerative Colitis

BSBG: Bile Salt Biotransformation Gene

BSBP: Bile Salt Biotransformation Protein

PSSN: Protein Sequence Similarity Network

KEGG: Kyoto Encyclopedia of Genes and Genomes

KOs: KEGG Orthologies

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and materials

All the relevant identifiers and other related information employed in this work have been provided in the additional files. Analyzed metagenomics datasets were obtained from https://www.ebi.ac.uk/ena/data/view/PRJEB2054 and https://www.ebi.ac.uk/ena/data/view/PRJNA389280.

Competing interests

The authors declare that they have no competing interests.

Funding

This study was supported by a grant from the European Commission FP7 project METACARDIS with the grant agreement HEALTH-F4-2012-305312. The funding body did not have any role in the study design, data collection, analysis and interpretation, the writing of the manuscript, or any influence on the content of the manuscript.

Authors’ contributions

PD, BJ, and JN designed the study. PD and SM created the reference dataset. PD performed, analyzed the metagenomics and metabolomics data, and wrote the draft. BJ and JN supervised the study. All authors read and approved the final manuscript.

Acknowledgments

We thank our research engineer Pierre-Etienne Cholley for his computational support. The computations were performed on resources at Chalmers Centre for Computational Science and Engineering (C3SE) provided by the Swedish National Infrastructure for Computing (SNIC).

1. Begley M, Gahan CGM, Hill C. The interaction between bacteria and bile. FEMS Microbiol Rev. 2005;29:625–51. doi:10.1016/j.femsre.2004.09.003.

2. De Aguiar Vallim TQ, Tarling EJ, Edwards PA. Pleiotropic roles of bile acids in metabolism. Cell Metab. 2013;17:657–69. doi:10.1016/j.cmet.2013.03.013.

3. Bjerrum JT, Wang Y, Hao F, Coskun M, Ludwig C, Günther U, et al. Metabonomics of human fecal extracts characterize ulcerative colitis, Crohn’s disease and healthy individuals. Metabolomics. 2015;11:122–33. doi:10.1007/s11306-014-0677-3.

4. Le Gall G, Noor SO, Ridgway K, Scovell L, Jamieson C, Johnson IT, et al. Metabolomics of Fecal Extracts Detects Altered Metabolic Activity of Gut Microbiota in Ulcerative Colitis and Irritable Bowel Syndrome. J Proteome Res. 2011;10:4208–18. doi:10.1021/pr2003598.

5. Jacobs JP, Goudarzi M, Singh N, Tong M, McHardy IH, Ruegger P, et al. A Disease-Associated Microbial and Metabolomics State in Relatives of Pediatric Inflammatory Bowel Disease Patients. Cell Mol Gastroenterol Hepatol. 2016;2:750–66. doi:10.1016/J.JCMGH.2016.06.004.

6. Jansson J, Willing B, Lucio M, Fekete A, Dicksved J, Halfvarson J, et al. Metabolomics Reveals Metabolic Biomarkers of Crohn’s Disease. PLoS One. 2009;4:e6386. doi:10.1371/journal.pone.0006386.

7. Huttenhower C, Kostic AD, Xavier RJ. Inflammatory Bowel Disease as a Model for Translating the Microbiome. Immunity. 2014;40:843–54. doi:10.1016/J.IMMUNI.2014.05.013.

8. Wahlström A, Sayin SI, Marschall H-U, Bäckhed F. Intestinal Crosstalk between Bile Acids and Microbiota and Its Impact on Host Metabolism. Cell Metab. 2016;24:41–50. doi:10.1016/J.CMET.2016.05.005.

9. Wahlström A, Kovatcheva-Datchary P, Ståhlman M, Bäckhed F, Marschall H-U. Crosstalk between Bile Acids and Gut Microbiota and Its Impact on Farnesoid X Receptor Signalling. Dig Dis. 2017;35:246–50. doi:10.1159/000450982.

10. Marcobal A, Kashyap PC, Nelson TA, Aronov PA, Donia MS, Spormann A, et al. A metabolomic view of how the human gut microbiota impacts the host metabolome using humanized and gnotobiotic mice. ISME J. 2013;7:1933–43. doi:10.1038/ismej.2013.89.

11. Mallonee DH, White WB, Hylemon PB. Cloning and sequencing of a bile acid-inducible operon from Eubacterium sp. strain VPI 12708. J Bacteriol. 1990;172:7011–9. http://www.ncbi.nlm.nih.gov/pubmed/2254270. Accessed 6 May 2018.

12. Mallonee DH, Hylemon PB. Sequencing and Expression of a Gene Encoding a Bile Acid Transporter from Eubacterium sp. Strain VPI 12708. J Bacteriol. 1996;178:7053–8. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC178615/pdf/1787053.pdf. Accessed 17 Apr 2018.

13. Kang D-J, Ridlon JM, Moore DR, Barnes S, Hylemon PB, Hylemon PB. Clostridium scindens baiCD and baiH genes encode stereo-specific 7alpha/7beta-hydroxy-3-oxo-delta4-cholenoic acid oxidoreductases. Biochim Biophys Acta. 2008;1781:16–25. doi:10.1016/j.bbalip.2007.10.008.

14. Ridlon JM, Harris SC, Bhowmik S, Kang D-J, Hylemon PB. Consequences of bile salt biotransformations by intestinal bacteria. Gut Microbes. 2016;7:22–39. doi:10.1080/19490976.2015.1127483.

15. Gothe F, Beigel F, Rust C, Hajji M, Koletzko S, Freudenberg F. Bile acid malabsorption assessed by 7 alpha-hydroxy-4-cholesten-3-one in pediatric inflammatory bowel disease: Correlation to clinical and laboratory findings. J Crohn’s Colitis. 2014;8:1072–8. doi:10.1016/j.crohns.2014.02.027.

16. Jones B V., Begley M, Hill C, Gahan CGM, Marchesi JR. Functional and comparative metagenomic analysis of bile salt hydrolase activity in the human gut microbiome. Proc Natl Acad Sci U S A. 2008;105:13580–5. doi:10.1073/pnas.0804437105.

17. Liu B, Pop M. ARDB—antibiotic resistance genes database. Nucleic Acids Res. 2008;37 suppl_1:D443–7.

18. Ridlon JM, Hylemon PB. Identification and characterization of two bile acid coenzyme A transferases from Clostridium scindens, a bile acid 7α-dehydroxylating intestinal bacterium. J Lipid Res. 2012;53:66–76. doi:10.1194/jlr.M020313.

19. Gilliland2 SE, Speck ML. Deconjugation of Bile Acids by Intestinal Lactobacilli1. Appl Environ Microbiol. 1977;33:15–8. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC170567/pdf/aem00012-0033.pdf. Accessed 3 Jun 2018.

20. Ruiz L, Margolles A, Sánchez B. Bile resistance mechanisms in Lactobacillus and Bifidobacterium. Front Microbiol. 2013;4:396. doi:10.3389/fmicb.2013.00396.

21. Frickey T, Lupas A. CLANS: a Java application for visualizing protein families based on pairwise similarity. Bioinforma Appl NOTE. 2004;20:3702–4. doi:10.1093/bioinformatics/bth444.

22. Labbé A, Ganopolsky JG, Martoni CJ, Prakash S, Jones ML. Bacterial Bile Metabolising Gene Abundance in Crohn’s, Ulcerative Colitis and Type 2 Diabetes Metagenomes. PLoS One. 2014;9:e115175. doi:10.1371/journal.pone.0115175.

23. Nagao M, Ohhira S, Kishi H, Komatsu W, Kobashi G, Uchiyama K. Lipid and Bile Acid Dysmetabolism in Crohn’s Disease. J Immunol Res. 2018;2018:1–6.

24. Amre DK, D’souza S, Morgan K, Seidman G, Lambrette P, Grimard G, et al. Imbalances in dietary consumption of fatty acids, vegetables, and fruits are associated with risk for Crohn’s disease in children. Am J Gastroenterol. 2007;102:2016.

25. Watanabe K, Igarashi M, Li X, Nakatani A, Miyamoto J, Inaba Y, et al. Dietary soybean protein ameliorates high-fat diet-induced obesity by modifying the gut microbiota-dependent biotransformation of bile acids. PLoS One. 2018;13:e0202083. doi:10.1371/journal.pone.0202083.

26. Becker C, Neurath MF, Wirtz S. The Intestinal Microbiota in Inflammatory Bowel Disease. ILAR J. 2015;56:192–204. doi:10.1093/ilar/ilv030.

27. Matsuoka K, Kanai T. The gut microbiota and inflammatory bowel disease. Semin Immunopathol. 2015;37:47–55. doi:10.1007/s00281-014-0454-4.

28. Hold GL, Smith M, Grange C, Watt ER, El-Omar EM, Mukhopadhya I. Role of the gut microbiota in inflammatory bowel disease pathogenesis: what have we learnt in the past 10 years? World J Gastroenterol. 2014;20:1192–210. doi:10.3748/wjg.v20.i5.1192.

29. Schirmer M, Franzosa EA, Lloyd-Price J, McIver LJ, Schwager R, Poon TW, et al. Dynamics of metatranscription in the inflammatory bowel disease gut microbiome. Nat Microbiol. 2018;3:337–46. doi:10.1038/s41564-017-0089-z.

30. The Integrative HMP Research Network Consortium. 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–89.

31. Franzosa EA, Sirota-Madi A, Avila-Pacheco J, Fornelos N, Haiser HJ, Reinker S, et al. Gut microbiome structure and metabolic activity in inflammatory bowel disease. Nat Microbiol. 2018;:1. doi:10.1038/s41564-018-0306-4.

32. Duboc H, Rajca S, Rainteau D, Benarous D, Maubert M-A, Quervain E, et al. Connecting dysbiosis, bile-acid dysmetabolism and gut inflammation in inflammatory bowel diseases. Gut. 2013;62:531–9. doi:10.1136/gutjnl-2012-302578.

33. Long SL, Gahan CGM, Joyce SA. Interactions between gut bacteria and bile in health and disease. Mol Aspects Med. 2017;56:54–65. doi:10.1016/J.MAM.2017.06.002.

34. Lundeen SG, Savage DC. Characterization of an extracellular factor that stimulates bile salt hydrolase activity in Lactobacillus sp. strain 100-100. FEMS Microbiol Lett. 1992;94:121–6. doi:10.1016/0378-1097(92)90594-E.

35. Travers M-A, Sow C, Zirah S, Deregnaucourt C, Chaouch S, Queiroz RML, et al. Deconjugated Bile Salts Produced by Extracellular Bile-Salt Hydrolase-Like Activities from the Probiotic Lactobacillus johnsonii La1 Inhibit Giardia duodenalis In vitro Growth. Front Microbiol. 2016;7:1453. doi:10.3389/fmicb.2016.01453.

36. Kakiyama G, Pandak WM, Gillevet PM, Hylemon PB, Heuman DM, Daita K, et al. Modulation of the fecal bile acid profile by gut microbiota in cirrhosis. J Hepatol. 2013;58:949–55. doi:10.1016/j.jhep.2013.01.003.

37. Bajaj JS, Heuman DM, Hylemon PB, Sanyal AJ, White MB, Monteith P, et al. Altered profile of human gut microbiome is associated with cirrhosis and its complications. J Hepatol. 2014;60:940–7. doi:10.1016/j.jhep.2013.12.019.

38. Ridlon JM, Alves JM, Hylemon PB, Bajaj JS. Cirrhosis, bile acids and gut microbiota: unraveling a complex relationship. Gut Microbes. 2013;4:382–7. doi:10.4161/gmic.25723.

39. Ajouz H, Mukherji D, Shamseddine A. Secondary bile acids: an underrecognized cause of colon cancer. World J Surg Oncol. 2014;12:164. doi:10.1186/1477-7819-12-164.

40. Keitel V, Donner M, Winandy S, Kubitz R, Häussinger D. Expression and function of the bile acid receptor TGR5 in Kupffer cells. Biochem Biophys Res Commun. 2008;372:78–84. doi:10.1016/J.BBRC.2008.04.171.

41. Pols TWH, Puchner T, Korkmaz HI, Vos M, Soeters MR, de Vries CJM. Lithocholic acid controls adaptive immune responses by inhibition of Th1 activation through the Vitamin D receptor. PLoS One. 2017;12:e0176715. doi:10.1371/journal.pone.0176715.

42. Martinot E, Sèdes L, Baptissart M, Lobaccaro J-M, Caira F, Beaudoin C, et al. Bile acids and their receptors. Mol Aspects Med. 2017;56:2–9.

43. Ward JBJ, Lajczak NK, Kelly OB, O’dwyer AM, Giddam AK, Gabhann JN, et al. Ursodeoxycholic acid and lithocholic acid exert anti-inflammatory actions in the colon. Am J Physiol Gastrointest Liver Physiol. 2017;312:550–8. doi:10.1152/ajpgi.00256.2016.

44. Tursi A, Brandimarte G, Papa A, Giglio A, Elisei W, Giorgetti GM, et al. Treatment of Relapsing Mild-to-Moderate Ulcerative Colitis With the Probiotic VSL#3 as Adjunctive to a Standard Pharmaceutical Treatment: A Double-Blind, Randomized, Placebo-Controlled Study. Am J Gastroenterol. 2010;105:2218–27. doi:10.1038/ajg.2010.218.

45. Ghouri YA, Richards DM, Rahimi EF, Krill JT, Jelinek KA, DuPont AW. Systematic review of randomized controlled trials of probiotics, prebiotics, and synbiotics in inflammatory bowel disease. Clin Exp Gastroenterol. 2014;7:473–87. doi:10.2147/CEG.S27530.

46. Begley M, Hill C, Gahan CGM. Bile salt hydrolase activity in probiotics. Appl Environ Microbiol. 2006;72:1729–38. doi:10.1128/AEM.72.3.1729-1738.2006.

47. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65. doi:10.1038/nature08821.

48. Arumugam M et al. Enterotypes of the human gut microbiome. Nature. 2011;473:174–80. doi:10.1038/nature09944.

49. Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, et al. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32 suppl_1:D115–9.

50. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10. doi:10.1016/S0022-2836(05)80360-2.

51. Hirano S, Masuda N. Transformation of bile acids by Eubacterium lentum. Appl Environ Microbiol. 1981;42:912–5. http://www.ncbi.nlm.nih.gov/pubmed/6947718. Accessed 9 Jun 2018.

52. Harris SC, Devendran S, Méndez- García C, Mythen SM, Wright CL, Fields CJ, et al. Bile acid oxidation by Eggerthella lenta strains C592 and DSM 2243 T. Gut Microbes. 2018;:1–17.

53. Yokota A, Fukiya S, Islam KBMS, Ooka T, Ogura Y, Hayashi T, et al. Is bile acid a determinant of the gut microbiota on a high-fat diet? Gut Microbes. 2012;3:455–9. doi:10.4161/gmic.21216.

54. Han SW, Evans DG, el-Zaatari FA, Go MF, Graham DY. The interaction of pH, bile, and Helicobacter pylori may explain duodenal ulcer. Am J Gastroenterol. 1996;91:1135–7. http://www.ncbi.nlm.nih.gov/pubmed/8651159. Accessed 6 Jun 2018.

55. Itoh M, Wada K, Tan S, Kitano Y, Kai J, Makino I. Antibacterial action of bile acids against Helicobacter pylori and changes in its ultrastructural morphology: effect of unconjugated dihydroxy bile acid. J Gastroenterol. 1999;34:571–6. http://www.ncbi.nlm.nih.gov/pubmed/10535483. Accessed 6 Jun 2018.

56. Finegold SM. Anaerobic Gram-Negative Bacilli. University of Texas Medical Branch at Galveston; 1996.

57. Shah HN, Collins MD, Olsen I, Paster BJ, Dewhirsp FE. Reclassification of Bacteroides Zevii (Holdeman, Cato, and Moore) in the Genus Polphyromonas, as Porphyromonas Zevii comb. nov. Int J Syst Bacteriol. 1995;45:586–8. http://www.microbiologyresearch.org/docserver/fulltext/ijsem/45/3/ijs-45-3-586.pdf?expires=1528287987&id=id&accname=guest&checksum=99A85DF43380B23AC251826BA8B82EEE. Accessed 6 Jun 2018.

58. Gillespie SH, Hawkey PM. Principles and practice of clinical bacteriology. John Wiley & Sons; 2006.

59. Lambert JM, Siezen RJ, de Vos WM, Kleerebezem M. Improved annotation of conjugated bile acid hydrolase superfamily members in Gram-positive bacteria. Microbiology. 2008;154:2492–500.

60. Panigrahi P, Sule M, Sharma R, Ramasamy S, Suresh CG. An improved method for specificity annotation shows a distinct evolutionary divergence among the microbial enzymes of the cholylglycine hydrolase family. Microbiology. 2014;160 Pt_6:1162–74.

61. Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44:D279–85. doi:10.1093/nar/gkv1344.

62. Jensen LJ, Julien P, Kuhn M, von Mering C, Muller J, Doerks T, et al. eggNOG: automated construction and annotation of orthologous groups of genes. Nucleic Acids Res. 2008;36 Database issue:D250-4. doi:10.1093/nar/gkm796.

63. Wickham H. ggplot2: elegant graphics for data analysis. Cham: Springer International Publishing; 2016.

64. Tyner S, Briatte F, Hofmann H. Network Visualization with ggplot2. R J. 2017.

65. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010.

66. Kim J, Kim MS, Koh AY, Xie Y, Zhan X. FMAP: Functional Mapping and Analysis Pipeline for metagenomics and metatranscriptomics studies. BMC Bioinformatics. 2016;17:420. doi:10.1186/s12859-016-1278-0.

67. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2011;40:D109–14.

68. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32 suppl_1:D277–80.

69. Team RC. R: A language and environment for statistical computing. 2016.

Table 1. Distribution of the total number of bile biotransformation protein homologs grouped by respective phylum.

Phylum	BaiA	BaiB	BaiCD	BaiE	BaiF	BaiG	BaiH	BaiI	BaiJ	BaiK	BaiL	Bsh	HSD
Actinobacteria	86	317	25	538	146	475	67	0	15	68	44	36	182
Bacteroidetes	158	1	2	1	0	1	2	0	3	0	23	4	264
Firmicutes	385	13	435	20	96	244	424	2	130	92	576	889	409
Fusobacteria	5	0	2	0	0	0	2	0	44	0	0	0	11
Proteobacteria	114	152	92	104	460	73	85	0	31	534	84	2	598
Verrucomicrobia	1	0	0	1	0	0	0	0	0	0	2	0	2

Download PDF

Journal Publication

published 29 Jun, 2019

Read the published version in BMC Genomics →

Submission checks completed at journal
22 May, 2019
Editorial decision: Major revision
21 May, 2019
Review #1 received at journal
20 May, 2019
Reviewer #1 agreed at journal
10 May, 2019
Editor invited by journal
07 May, 2019
Editor assigned by journal
07 May, 2019
Reviewers invited by journal
07 May, 2019
First submitted to journal
29 Mar, 2019

You are reading this older preprint version

Read the latest preprint version →

Metagenomic analysis of bile salt biotransformation in the human gut microbiome

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Results

Identification and distribution of bile salt biotransformation protein homologs

Protein-sequence similarity network of bile salt biotransformation protein homologs

Comparative metagenomic analysis of bile salt biotransformation gene abundances between healthy and IBD subjects

Comparative metagenomic analysis of bile salt biotransformation gene abundances among healthy subjects from different countries

Phyla-level comparative metagenomic analysis of bile salt biotransformation gene abundances between healthy and IBD subjects

Comparative fecal metabolomics analysis of bile acids between healthy and IBD subjects

Discussion

Conclusions

Methods

IBD cohort description

Bioinformatic identification of bile salt biotransformation protein homologs

Construction of protein sequence similarity network (PSSN)

Functional metagenomic mapping

Statistical analysis

Abbreviations

Declarations

Ethics approval and consent to participate

Consent for publication

Availability of data and materials

Competing interests

Funding

Authors’ contributions

Acknowledgments

References

Table 1

Supplementary Files

Status:

Journal Publication

Version 1