Gut microbiota causally affects cholelithiasis: a two-sample Mendelian randomization study
Xin Liu 
Xingsi Qi
Zibin Tian- Department of Gastroenterology, The Affiliated Hospital of Qingdao University, Qingdao, China
Background: The gut microbiota is closely linked to cholesterol metabolism-related diseases such as obesity and cardiovascular diseases. However, whether gut microbiota plays a causal role in cholelithiasis remains unclear.
Aims: This study explored the causal relationship between gut microbiota and cholelithiasis. We hypothesize that the gut microbiota influences cholelithiasis development.
Methods: A two-sample Mendelian randomization method was combined with STRING analysis to test this hypothesis. Summary data on gut microbiota and cholelithiasis were obtained from the MiBioGen (n=13,266) and FinnGen R8 consortia (n=334,367), respectively.
Results: Clostridium senegalense, Coprococcus3, and Lentisphaerae increased the risk of cholelithiasis and expressed more bile salt hydrolases. In contrast, Holdemania, Lachnospiraceae UCG010, and Ruminococcaceae NK4A214 weakly expressed bile salt hydrolases and were implied to have a protective effect against cholelithiasis by Mendelian randomization analysis.
Conclusion: Gut microbiota causally influences cholelithiasis and may be related to bile salt hydrolases. This work improves our understanding of cholelithiasis causality to facilitate the development of treatment strategies.
Introduction
Cholelithiasis (also known as gallstones) is defined as a solid clot in the gallbladder or biliary system (Lammert et al., 2016). Approximately 90% of cholelithiasis occurrences are cholesterol gallstones, while the incidence of other stone types (including black and brown pigment stones) is below 10% (Sun et al., 2022). Approximately 10–20% of the global population has gallstones, and over 20% of cases develop gallstone diseases, such as acute cholecystitis, acute cholangitis, and obstructive jaundice (Cortés et al., 2020). Gallstone disease is one of the most expensive gastrointestinal conditions from a societal perspective (Grigor’eva and Romanova, 2020). Epidemiological studies identified numerous risk factors for cholesterol stones, which include type 2 diabetes, physical inactivity, and over-nutrition (Di Ciaula et al., 2019). This finding can be attributed to the risk factors that lead to excess cholesterol or disruption of cholesterol homeostasis (Rudling et al., 2019).
The gut microbiota has a non-negligible impact on metabolic disorders, including insulin resistance (Gomes et al., 2018), obesity (Abenavoli et al., 2019), and hyperlipidemia (Michels et al., 2022). These disorders are known risk factors for increased hepatic cholesterol synthesis, gallstone formation, and symptomatic gallstones (Di Ciaula et al., 2019). However, the role of intestinal microbiota in gallstone development remains unclear. Elevated levels of cholesterol and bilirubin in the bile and decreased levels of bile salts cause cholesterol gallstones (Sun et al., 2022). Decreased bile salt levels are observed in liver disease and in conditions such as Crohn’s disease or in individuals that have undergone colectomy or intestinal resection, where the enterohepatic circulation of bile salts is impaired (Molina-Molina et al., 2018). These findings led us to hypothesize that the gut microbiota and intestine may play key roles in influencing gallstone formation in the host. The diversity and taxonomy of gut microbiota are associated with bile acid levels in gallstone disease, and increased concentrations of taurodeoxycholic acid and taurocholic acid are associated with the presence of conditionally pathogenic bacteria (Petrov et al., 2020). However, intestinal flora is often influenced by confounding factors, including age, sex, environment, alcohol consumption, diet, and lifestyle (Zmora et al., 2019). Eliminating these confounders in observational studies can be challenging, thereby limiting research on the causal role of gut microbiota in gallstones. Fortunately, Mendelian randomization (MR) analysis can be employed to explore the causal role of intestinal microbiota in the etiology of human diseases independent of confounding factors (O’Donnell et al., 2023).
MR borrows economically inspired statistical techniques to enable researchers to examine the causal factors affecting human diseases (Birney, 2022). Numerous risk factors related to diseases have not established causation owing to the limitations of observational studies that cannot avoid the influence of confounders. Thus, MR has become one of the most effective methods for addressing issues in human biology and epidemiology, including the relationship between intestinal microbiota and disease (Georgakis and Gill, 2021). The correlation between genetic variation and outcome is independent of confounders in MR analysis (Bowden and Holmes, 2019). For instance, MR analysis has revealed that Bifidobacterium is causally linked to preeclampsia-eclampsia (Li et al., 2022), and the fecal abundance of Oscillibacter and Alistipes is causally associated with decreased triglyceride levels (Liu et al., 2022).
This study used summary-level statistics of the genome-wide association study (GWAS) from the MiBioGen and FinnGen consortia to perform a two-sample MR analysis to investigate the causal relationship between intestinal microbiota and cholelithiasis.
Materials and methods
The two-sample MR study relied on three assumptions to draw conclusions regarding causation. Genetic variants strongly predict microbiome exposure independent of confounding factors and outcomes. Genetic variants influence outcomes through exposure (Emdin et al., 2017).
Data sources
Single nucleotide polymorphisms (SNPs) related to intestinal microbiota were obtained from the GWAS dataset of the International Consortium MiBioGen and were used as instrumental variables (IVs). This study included 34,024 individuals from 18 cohorts predominantly of European ancestry (including those from the United States, Canada, Israel, South Korea, Germany, Denmark, The Netherlands, Belgium, Sweden, Finland, and the United Kingdom). The dataset provided genotyping data that coordinated 16S ribosomal RNA gene sequencing to examine the relationship between genetic variants and intestinal microbiota by profiling taxonomic classification. A total of 211 taxa were included in the analysis. We selected SNPs that showed significant correlations with genus at a suggestive of genome-wide significance thresholds (P < 1× 10−5, F >10) as potential IVs from MiBioGen (Li et al., 2022).
Data outcome
GWAS summary statistics for cholelithiasis were acquired from the FinnGen Consortium R8 release data. The phenotype “cholelithiasis” was adopted in our research. The GWAS consisted of 32,894 cases and 301,383 controls. The mean age of the patients was 52.11 years old (female: 48.62, male: 59.36). The principal components (sex and age), and the genotyping batch were corrected during the analysis. SNPs were identified with genome-wide significant correlations with cholelithiasis (P < 5 × 10−8).
Instrumental variable filtering
First, SNPs showing significant correlations with the genus were selected at suggestive genome-wide significance thresholds from the MiBioGen as potential IVs (P < 1× 10−5, F >10). A chain imbalance threshold (R2 < 0.001) and linkage disequilibrium threshold of 10,000 kbp was then applied to ensure independence among the selected SNPs. Single nucleotide polymorphisms that were unclear, duplicated, or palindromic were removed to ensure consistent SNP orientation in the exposure and results. Next, SNPs should have P >1×10-5 to ensure that SNPs are independent of the outcome. Finally, the PhenoScanner online tool was used to check whether the SNPs affected the outcome, followed by their manual exclusion. The MR-Egger intercept analysis was employed to test for horizontal pleiotropy (P> 0.05) and leave-one-out analysis was conducted to assess whether each IV affected the overall estimates of the remaining IVs.
MR analysis
The most frequently used MR methods were used to analyze the causal relationship between intestinal microbiota and cholelithiasis. They were inverse variance weighted (IVW), MR-Egger, weighted median, weighted mode, and simple mode.
Reverse MR analysis
MR inverse analysis was performed to explore whether cholelithiasis had a causal effect on the identified microbiomes. During this analysis, cholelithiasis and bacteria was considered as an exposure and an outcome, respectively, SNPs with notably correlated with cholelithiasis were taken as IVs.
All data analysis was used R software (version 4.2.1) to conduct. We applied the IVW, weighted median, and MR-Egger regression methods using the R packages of “TwoSampleMR” (version 0.5.6). And “MRPRESSO” package was used to perform the MR-PRESSO analysis.
SRTING analysis
The STRING database collects and integrates protein-protein interaction (PPI), performs enrichment analysis and highlights proteins in the PPI network (Szklarczyk et al., 2023). We tried to use the STRING to predict the function of BSH in the bacteria (https://string-db.org), and conduct ontology (GO)and KEGG enrichment analysis. Statistical significance was defined as p<0.05.
Results
Single nucleotide polymorphism selection
A total of 2557 SNPs related to 211 taxa were identified as IVs of the gut microbiota. A series of quality control steps was performed resulting in the selection of 72 SNPs associated with six genera and one phylum (based on the IVW P-value < 0.05). Analysis of these 72 SNPs using PhenoScanner showed that no SNPs were related to confounding factors. Heterogeneity between the two samples was tested using Cochran’s Q statistics and no evidence of heterogeneity (p > 0.05) or horizontal pleiotropy of the IVs was observed (MRPRESSO-global, p > 0.05; MR-Egger intercept, p > 0.05) (Table S1). The MR analysis results are shown in Figure 1.
Figure 1 Scatter plots of Mendelian randomization (MR) analysis. The scatter plot of Clostridiumsenegalense (A), Coprococcus3 (B), Holdemania (C), Lachnospiraceae UCG010 (D), Ruminococcaceae 617 NK4A214 (E), and Lentisphaerae (F).
Detailed MR results
Five genera (Clostridiumsenegalense, Coprococcus3, Holdemania, Lachnospiraceae UCG010, and Ruminococcaceae NK4A214) and one phylum (Lentisphaerae) of the microbiome were suggested to have causal relationships with cholelithiasis. IVW analysis suggested that Clostridiumsenegalense (OR=1.244, 95% CI:1.091–1.419), Coprococcus3 (OR=1.236, 95% CI: 1.047–1.458), and the Lentisphaerae phylum (OR=1.081, 95% CI: 1.004–1.164) increased the risk of cholelithiasis. Meanwhile, Holdemania (OR=0.911, 95% CI: 0.841–0.987), Lachnospiraceae UCG010 (OR=0.835, 95% CI:0.740–0.942), and Ruminococcaceae NK4A214 (OR=0.860, 95% CI:0.768–0.963) had a protective effect on cholelithiasis (Table 1).
Table 1 The MR analysis of causal effects between gut microbiota and cholelithiasis.
Sensitivity analysis
Sensitivity analysis is necessary to assess the effectiveness of IVWs. The MR-Egger method was used to assess horizontal multi-effectiveness. The MR-Egger intercept and MR-PRESSO global tests indicated a low likelihood of horizontal multiplicity (Table S1, p > 0.05). All I2 values in the heterogeneity tests were <50%, and all p-values were >0.05. This finding indicated that our findings were probably not influenced by heterogeneity bias. The MR-Egger intercept and the MR-PRESSO global test showed no significant horizontal pleiotropy. This result indicated that the outliers did not significantly affect the results. Simultaneously, the consistency of the robust IVW (adjusted for the effect of outliers) and MR-PRESSO (adjusted for the effect of horizontal multiplicity) supported the lack of significant outlier effects on the results (Table S1). A leave-one-out analysis was also conducted to assess whether each IV affects the overall estimates of the remaining IVs. These results suggested that none of the single IV treatments affected the results (Figure 2).
Figure 2 Leave-one-out sensitivity analysis. Leave-one-out sensitivity analysis of Clostridiumsenegalense (A), Coprococcus3 (B), Holdemania (C), Lachnospiraceae UCG010 (D), Ruminococcaceae NK4A214 (E), and Lentisphaerae (F).
The causal effects of gut microbiota and cholelithiasis via inverse MR analysis
The reverse causal effects were examined using cholelithiasis as the exposure and intestinal microbiota as the outcome. Forty-three SNPs associated with cholelithiasis were selected as IVs (P < 5 × 10-8). Cholelithiasis was causally associated with Actinobacteria, Lachnospiraceae, Clostridium innocuum, Eggerthella, Eubacterium brachy, Intestinimonas, Paraprevotella, and Mollicutes RF9 (Table 2). None of the bacteria exhibited a bidirectional causal relationship with cholelithiasis and gut bacteria.
Table 2 Inverse MR analysis the causal effects of gut microbiota and cholelithiasis.
Microbiota causally linked to cholelithiasis may be associated with bile salt hydrolase
Bile acids (BAs) are divided into primary and secondary categories. Primary BAs are excreted into the intestine and converted into secondary BAs by intestinal microorganisms (Sinha et al., 2020). The initial step involves hydrolysis of the amino acid fraction by BSH during secondary BA metabolism (Smirnova et al., 2022). BSH (also known as choloylglycine hydrolase) is present in the intestinal microbiome to maintain BA balance. Imbalances in BAs are associated with gallstones, gallbladder disease, obesity, and diabetes (Cai et al., 2022a). BSHs are highly conserved in all major gut microbial phyla (including Bacteroidetes, Firmicutes, and Actinobacteria); however, they are bacterially different owing to their preferential activity toward glycine- or taurine-conjugated BAs. The Human Microbiome Project reported that 26.03% bacterial strains encode BSHs (Song et al., 2019). Thus, gut microbiota-related BSH directly determines the synthesis of secondary bile acids, which involve in regulating cholesterol metabolism. We hypothesized that BSH may be a link between intestinal microbiota and gallstone causation. Thus, we selected BSHs of the microbiome (that are associated with cholelithiasis by MR analysis) in the Human Microbiome Project database and National Center for Biotechnology Information database for protein–protein interaction (PPI) analysis. Clostridium and Coprococcus were predicted as a risk factor for cholelithiasis by MR and they expressed more BSHs than the protective bacteria (Holdemania, Lachnospiraceae UCG010, and Ruminococcaceae NK4A214) (Figures 3A, B, D, E). Similarly, inverse MR analysis predicted cholelithiasis as a risk factor for Intestinimonas, which expresses BSH (Figures 3C, F). We also used STRING to predicted the function of Clostridium, Coprococcus and Intestinimonas. Enrichment analysis showed that the bile acid catabolic process and secondary bile acid biosynthesis were enriched in Clostridium and Coprococcus according to Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis. Intestinimonas was only enriched for secondary bile acid biosynthesis (Table S2). No enrichment of secondary bile acid biosynthesis in microbiota was observed that did not contain the BSH protein. These results suggest that BSH may serve as a link between gut microbiota and cholelithiasis (Figure 4).
Figure 3 The number of bile salt hydrolases (BSHs) in microbiota. (A) The number of BSHs in microbiota genera base on HMP database. (B) The number of BSHs in microbiota genera based on Mendelian randomization (MR) analysis and reverse MR analysis (C). The number of BSH- related proteins in (D) Coprococcus, (E) Clostridium, (F) Intestinimonas, red represents BSH, blue and green dots represent BSH-related proteins.
Figure 4 The framework of two-sample Mendelian randomization analysis results. Gut microbiota causally influences cholelithiasis, while cholelithiasis impacts the composition of gut microbiota. GCA, glycocholic acid; TCA, taurocholic acid. GCDCA, Glycochenodeoxycholic acid; TCDCA, Taurochnodeoxycholic acid; CA, Cholic acid; CDCA, Chenodeoxycholic acid; DCA, Deoxycholic acid; LCA, Lithocholic acid; UDCA, Ursodeoxycholic acid; BSH, bile salt hydrolase.
Microbiota linked to serum total cholesterol may be associated with BSH
Lower total cholesterol levels in serum may be an independent risk factor for cholelithiasis (Chen et al., 2022). It has been proved that one of the effects of BSH on host is cholesterol lowering (Ertürkmen et al., 2023), and Terrisporobacter was associated with higher total cholesterol levels (Guo et al., 2023). We detected the expression of HBS in Terrisporobacter, the red node represented BSH, and the functional enrichments was predicted by STRING, no enrichments were found (PPI enrichment P value: 0.457) and the function was described protein lipoylation and CoA hydrolase activity rather than bile acid or fatty acid metabolism (Figure 5). the result indicated that BSH of Terrisporobacter may be little function in bile acid or fatty acid metabolism.
Figure 5 The analysis of Terrisporobacter othiniensis by protein-protein interaction (PPI). (A) The number of BSHs in Terrisporobacter othiniensis, red represents BSH. (B) Functional prediction of Terrisporobacter othiniensis by STRING. BSH, bile salt hydrolase; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein-protein interaction.
Discussion
This study used the MiBioGen database and cholelithiasis data from FinnGen to investigate the causal relationship between gut microbes and cholelithiasis using MR analysis. Clostridium senegalense, Coprococcus3, and Lentisphaerae increase the risk of cholelithiasis. In contrast, Holdemania, Lachnospiraceae UCG010, and Ruminococcaceae NK4A214 showed protective effects. The causal relationship identified BSH as a potential link between bile salt metabolism and the gut microbiota.
The intestinal microbiota is a metabolic organ that produces numerous metabolites (Connell et al., 2022) (including BAs and indole derivatives) (Agus et al., 2018) that play crucial roles in regulating host metabolism (Cai et al., 2022b). Cholesterol oxidized by liver enzymes results in the production of BAs that are further metabolized by the intestinal microbiota (Collins et al., 2023). The gut microbiota regulates key enzymes, such as cholesterol-7α-hydroxylase (CYP7A1), involved in BA synthesis (Hartmann et al., 2018; Fukui, 2021), and BA synthesis is tightly controlled by negative feedback inhibition through the farnesoid X receptor (FXR) (Jia et al., 2018). Bile acid deconjugation is primarily mediated by bacteria with BSH activity (Smirnova et al., 2022). Thus, intestinal microbiota plays a significant role in the key enzymatic processes involved in cholesterol synthesis in the liver, secondary BA production, and enterohepatic circulation of BAs.
The balance between cholesterol and bile salts is a critical factor for gallstones (Ye et al., 2021). Previous studies have revealed direct associations between taurocholic acid, taurochenodeoxycholic acid, and alpha diversity of the microbiota, together with positive associations with the genera Chitinophagaceae, Microbacterium, Lutibacterium, and Prevotella intermedia (Petrov et al., 2020). Patients with gallstones exhibit an increased richness of 7a-dehydroxylating microbiota and decreased levels of Fimicutes and diversity of gut microbiota (Wang et al., 2020). Additionally, the abundance of Lactobacillus strains significantly reduced in lithogenic diet-induced gallstones through the mediation of FXR signaling (Ye et al., 2022). However, the causal correlation between the intestinal flora and cholelithiasis remains unclear. Our study predicted that presence of abundance of Clostridium senegalense and Coprococcus3 have a causal relationship with cholelithiasis.
Many observational studies report a correlation between Clostridium and gallstone disease. For instance, patients with gallstones have elevated levels of Clostridium in the stool (Grigor’eva and Romanova, 2020), and Clostridium species were isolated from gallbladder stones (Liu et al., 2000), consistent with the results of our study. We found a causal association between Clostridium senegalense and cholelithiasis. Clostridium-encoded protein analysis revealed that Clostridium had two genes of BSH consistent with a previous study (Song et al., 2019). The presence of greater amounts of bacteria with BSH leads to increased bile salt deconjugation, resulting in elevated biliary deoxycholate levels, positive regulation of hepatic cholesterol secretion, and cholesterol crystallization (Lammert et al., 2016). Furthermore, increased hydrolysis can lead to steatosis, and increased secondary bile acid levels are associated with colorectal cancer (Horáčková et al., 2018). Recent research showing that theabrownin can reduce hypercholesterolemia by inhibiting the intestinal microbiota related to BSH activity indirectly supports our hypothesis. The underlying mechanism may be that low BSH activity increased the concentration of conjugated BAs, which could inhibit the FXR-FGF15 signaling pathway, thereby led to decreased hepatic cholesterol synthesis (Huang et al., 2019). However, the effect of BSH on the microbiota in cholelithiasis causation requires verification and further exploration.
Coprococcus was also predicted to be a risk factor for gallstone disease, and encoded BSH proteins which was similar to the previous study (Song et al., 2019). Moreover, the relative abundance of BSH in the microbiota is significantly associated with mortality from diabetes and cardiovascular disease (Song et al., 2019). This finding further emphasizes the significance of BSH in cholelithiasis. Further research is required to understand how BSH in Clostridium and Coprococcus affect the formation of cholesterol stones. Recent studies show that Desulfovibrionales are enriched in cholelithiasis patients, increase the synthesis of secondary BAs and intestinal cholesterol uptake, stimulate biliary secretion, and affect FXR and CYP7A expression [15]. Further research is required to determine whether BSHs in Clostridium and Coprococcus promote cholesterol stone formation by increasing the secondary BA synthesis.
The present study revealed the protective effects of Lachnospiraceae UCG010 and Ruminococcaceae NK4A214 against cholelithiasis. These strains did not express BSH. This result suggests that BSH deficiency might be a protective factor against gallstone formation; however, this hypothesis requires verification. Microbial conversion of cholesterol to coprostanol is another mechanism that reduces cholesterol and decreases the formation of cholesterol stones (Kenny et al., 2020). Lactobacillus curvatus KFP419 strain reduce cholesterol levels by increasing the conversion of cholesterol to coprostanol (Park et al., 2018). Lachnospiraceae and Ruminococcaceae are reportedly associated with high coprostanol levels (Antharam et al., 2016). An elevation of sterol metabolites coincides with increases in the Lachnospiraceae family in vitro, implying that this family promotes sterol metabolism (Blanco-Morales et al., 2020). These findings suggest that the protective effects of Lachnospiraceae and Ruminococcaceae against cholelithiasis may be related to cholesterol conversion; however, the underlying mechanisms require further exploration.
To the best of our knowledge, this is the first MR investigation revealing a link between gut microbiota and cholelithiasis. Robust genetic instruments were employed with a test of horizontal pleiotropy and two-sample heterogeneity to obtain reliable results from the MR analysis. Additionally, a leave-one-out analysis was performed to examine potential biases introduced by individual SNPs. Five sets of genetic instruments were used for MR analysis. In addition, the concept of BSH was introduced as an innovative approach to investigate the underlying mechanisms of the causal association between microbiota and cholelithiasis.
Nevertheless, this study has limitations. First, most patients with cholelithiasis in the analysis were of European descent, whereas the gut flora database encompasses other populations. Second, the potential association between microbiota and cholelithiasis by BSH was only superficially speculated. Further research is required to elucidate the underlying mechanisms and causal relationship. Third, subgroup analyses such as distinguishing between symptomatic and asymptomatic cholelithiasis was not possible because summary data for cholelithiasis was used in our analysis. Fourth, our exploration was limited to the genus level owing to the lowest taxonomic level available in the gut microbiota dataset, thus impeding a more detailed investigation at the species level. Additionally, some bacteria were predicted to be present at the phylum level; we could not analyze their BSH proteins. Finally, the sample size of the exposure group was relatively small; therefore, the reverse MR analysis results could not completely exclude the possibility of reverse causality.
Conclusion
Our investigations demonstrated a causal relationship between cholelithiasis and Clostridium senegalense and Coprococcus3, whereas Holdemania, Lachnospiraceae UCG010, and Ruminococcaceae NK4A214 had a protective effect. The causal relationship between the gut microbiota and cholelithiasis may be mediated by BSH. In addition, reverse MR analysis supported a causal relationship between cholelithiasis and the intestinal microbiota. Moreover, these findings suggest that cholelithiasis may influence the gut flora. Further validation and mechanistic studies are required.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
XL: Data curation, Formal Analysis, Investigation, Methodology, Resources, Visualization, Writing – original draft, Writing – review & editing. XQ: Data curation, Formal Analysis, Methodology, Writing – review & editing. RH: Data curation, Validation, Writing – original draft. TM: Supervision, Writing – review & editing. ZT: Funding acquisition, Investigation, Project administration, Writing – review & editing.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 81970461).
Acknowledgments
We thank all the participants, including the patients, the healthy control population, and the investigators that were devoted to data collation, analysis, summary, and public accessibility.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1253447/full#supplementary-material
References
Abenavoli, L., Scarpellini, E., Colica, C., Boccuto, L., Salehi, B., Sharifi-Rad, J., et al. (2019). Gut microbiota and obesity: A role for probiotics. Nutrients 11. doi: 10.3390/nu11112690
Agus, A., Planchais, J., Sokol, H. (2018). Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe 23, 716–724. doi: 10.1016/j.chom.2018.05.003
Antharam, V. C., McEwen, D. C., Garrett, T. J., Dossey, A. T., Li, E. C., Kozlov, A. N., et al. (2016). An integrated metabolomic and microbiome analysis identified specific gut microbiota associated with fecal cholesterol and coprostanol in clostridium difficile infection. PloS One 11, e0148824. doi: 10.1371/journal.pone.0148824
Birney, E. (2022). Mendelian randomization. Cold Spring Harb. Perspect. Med. 12, a041302. doi: 10.1101/cshperspect.a041302
Blanco-Morales, V., Garcia-Llatas, G., Yebra, M. J., Sentandreu, V., Lagarda, M. J., Alegría, A. (2020). Impact of a plant sterol- and galactooligosaccharide-enriched beverage on colonic metabolism and gut microbiota composition using an in vitro dynamic model. J. Agric. Food Chem. 68, 1884–1895. doi: 10.1021/acs.jafc.9b04796
Bowden, J., Holmes, M. V. (2019). Meta-analysis and Mendelian randomization: A review. Res. Synth. Methods 10, 486–496. doi: 10.1002/jrsm.1346
Cai, J., Rimal, B., Jiang, C., Chiang, J. Y. L., Patterson, A. D. (2022a). Bile acid metabolism and signaling, the microbiota, and metabolic disease. Pharmacol. Ther. 237, 108238. doi: 10.1016/j.pharmthera.2022.108238
Cai, J., Sun, L., Gonzalez, F. J. (2022b). Gut microbiota-derived bile acids in intestinal immunity, inflammation, and tumorigenesis. Cell Host Microbe 30, 289–300. doi: 10.1016/j.chom.2022.02.004
Chen, L., Yang, H., Li, H., He, C., Yang, L., Lv, G. (2022). Insights into modifiable risk factors of cholelithiasis: A Mendelian randomization study. Hepatol. Baltim. Md 75, 785–796. doi: 10.1002/hep.32183
Collins, S. L., Stine, J. G., Bisanz, J. E., Okafor, C. D., Patterson, A. D. (2023). Bile acids and the gut microbiota: metabolic interactions and impacts on disease. Nat. Rev. Microbiol. 21, 236–247. doi: 10.1038/s41579-022-00805-x
Connell, E., Le Gall, G., Pontifex, M. G., Sami, S., Cryan, J. F., Clarke, G., et al. (2022). Microbial-derived metabolites as a risk factor of age-related cognitive decline and dementia. Mol. Neurodegener. 17, 43. doi: 10.1186/s13024-022-00548-6
Cortés, V. A., Barrera, F., Nervi, F. (2020). Pathophysiological connections between gallstone disease, insulin resistance, and obesity. Obes. Rev. Off. J. Int. Assoc. Study Obes. 21, e12983. doi: 10.1111/obr.12983
Di Ciaula, A., Garruti, G., Frühbeck, G., De Angelis, M., de Bari, O., Wang, D. Q.-H., et al. (2019). The role of diet in the pathogenesis of cholesterol gallstones. Curr. Med. Chem. 26, 3620–3638. doi: 10.2174/0929867324666170530080636
Emdin, C. A., Khera, A. V., Kathiresan, S. (2017). Mendelian randomization. JAMA 318, 1925–1926. doi: 10.1001/jama.2017.17219
Ertürkmen, P., Fırıncıoğulları, B., Öner, Z. (2023). The Expression Levels of Genes Responsible for the Enzymatic Activity of Bile Salt Hydrolase (BSH) and the Relationship of Cholesterol Assimilation in L. plantarum and L. paracasei. Curr. Microbiol. 80, 205. doi: 10.1007/s00284-023-03311-2
Fukui, H. (2021). Leaky gut and gut-liver axis in liver cirrhosis: clinical studies update. Gut Liver 15, 666–676. doi: 10.5009/gnl20032
Georgakis, M. K., Gill, D. (2021). Mendelian randomization studies in stroke: exploration of risk factors and drug targets with human genetic data. Stroke 52, 2992–3003. doi: 10.1161/STROKEAHA.120.032617
Gomes, A. C., Hoffmann, C., Mota, J. F. (2018). The human gut microbiota: Metabolism and perspective in obesity. Gut Microbes 9, 308–325. doi: 10.1080/19490976.2018.1465157
Grigor’eva, I. N., Romanova, T. I. (2020). Gallstone disease and microbiome. Microorganisms 8, 835. doi: 10.3390/microorganisms8060835
Guo, G., Wu, Y., Liu, Y., Wang, Z., Xu, G., Wang, X., et al. (2023). Exploring the causal effects of the gut microbiome on serum lipid levels: A two-sample Mendelian randomization analysis. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1113334
Hartmann, P., Hochrath, K., Horvath, A., Chen, P., Seebauer, C. T., Llorente, C., et al. (2018). Modulation of the intestinal bile acid/farnesoid X receptor/fibroblast growth factor 15 axis improves alcoholic liver disease in mice. Hepatol. Baltim. Md 67, 2150–2166. doi: 10.1002/hep.29676
Horáčková, Š., Plocková, M., Demnerová, K. (2018). Importance of microbial defence systems to bile salts and mechanisms of serum cholesterol reduction. Biotechnol. Adv. 36, 682–690. doi: 10.1016/j.bioteChadv.2017.12.005
Huang, F., Zheng, X., Ma, X., Jiang, R., Zhou, W., Zhou, S., et al. (2019). Theabrownin from Pu-erh tea attenuates hypercholesterolemia via modulation of gut microbiota and bile acid metabolism. Nat. Commun. 10, 4971. doi: 10.1038/s41467-019-12896-x
Jia, W., Xie, G., Jia, W. (2018). Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis. Nat. Rev. Gastroenterol. Hepatol. 15, 111–128. doi: 10.1038/nrgastro.2017.119
Kenny, D. J., Plichta, D. R., Shungin, D., Koppel, N., Hall, A. B., Fu, B., et al. (2020). Cholesterol metabolism by uncultured human gut bacteria influences host cholesterol level. Cell Host Microbe 28, 245–257.e6. doi: 10.1016/j.chom.2020.05.013
Lammert, F., Gurusamy, K., Ko, C. W., Miquel, J.-F., Méndez-Sánchez, N., Portincasa, P., et al. (2016). Gallstones. Nat. Rev. Dis. Primer 2, 16024. doi: 10.1038/nrdp.2016.24
Li, P., Wang, H., Guo, L., Gou, X., Chen, G., Lin, D., et al. (2022). Association between gut microbiota and preeclampsia-eclampsia: a two-sample Mendelian randomization study. BMC Med. 20, 443. doi: 10.1186/s12916-022-02657-x
Liu, X., Tong, X., Zou, Y., Lin, X., Zhao, H., Tian, L., et al. (2022). Mendelian randomization analyses support causal relationships between blood metabolites and the gut microbiome. Nat. Genet. 54, 52–61. doi: 10.1038/s41588-021-00968-y
Liu, Y., Lam, K., Tang, Y. J., Gumerlock, P. H., Lee, D. K., Kim, M. H., et al. (2000). Anaerobic bacteria and intrahepatic stones: detections of Clostridium sp. and Bacteroides fragilis. Chin. Med. J. (Engl.) 113, 858–861.
Michels, N., Zouiouich, S., Vanderbauwhede, B., Vanacker, J., Indave Ruiz, B. I., Huybrechts, I. (2022). Human microbiome and metabolic health: An overview of systematic reviews. Obes. Rev. Off. J. Int. Assoc. Study Obes. 23, e13409. doi: 10.1111/obr.13409
Molina-Molina, E., Lunardi Baccetto, R., Wang, D. Q.-H., de Bari, O., Krawczyk, M., Portincasa, P. (2018). Exercising the hepatobiliary-gut axis. The impact of physical activity performance. Eur. J. Clin. Invest. 48, e12958. doi: 10.1111/eci.12958
O’Donnell, J. A., Zheng, T., Meric, G., Marques, F. Z. (2023). The gut microbiome and hypertension. Nat. Rev. Nephrol. 19, 153–167. doi: 10.1038/s41581-022-00654-0
Park, M. Y., Kim, J., Kim, S., Whang, K.-Y. (2018). Lactobacillus curvatus KFP419 and Leuconostoc mesenteroides subsp. mesenteroides KDK411 Isolated from Kimchi Ameliorate Hypercholesterolemia in Rats. J. Med. Food 21, 647–653. doi: 10.1089/jmf.2017.4125
Petrov, V. A., Fernández-Peralbo, M. A., Derks, R., Knyazeva, E. M., Merzlikin, N. V., Sazonov, A. E., et al. (2020). Biliary microbiota and bile acid composition in cholelithiasis. BioMed. Res. Int. 2020, 1242364. doi: 10.1155/2020/1242364
Rudling, M., Laskar, A., Straniero, S. (2019). Gallbladder bile supersaturated with cholesterol in gallstone patients preferentially develops from shortage of bile acids. J. Lipid Res. 60, 498–505. doi: 10.1194/jlr.S091199
Sinha, S. R., Haileselassie, Y., Nguyen, L. P., Tropini, C., Wang, M., Becker, L. S., et al. (2020). Dysbiosis-induced secondary bile acid deficiency promotes intestinal inflammation. Cell Host Microbe 27, 659–670.e5. doi: 10.1016/j.chom.2020.01.021
Smirnova, E., Muthiah, M. D., Narayan, N., Siddiqui, M. S., Puri, P., Luketic, V. A., et al. (2022). Metabolic reprogramming of the intestinal microbiome with functional bile acid changes underlie the development of NAFLD. Hepatol. Baltim. Md 76, 1811–1824. doi: 10.1002/hep.32568
Song, Z., Cai, Y., Lao, X., Wang, X., Lin, X., Cui, Y., et al. (2019). Taxonomic profiling and populational patterns of bacterial bile salt hydrolase (BSH) genes based on worldwide human gut microbiome. Microbiome 7, 9. doi: 10.1186/s40168-019-0628-3
Sun, H., Warren, J., Yip, J., Ji, Y., Hao, S., Han, W., et al. (2022). Factors influencing gallstone formation: A review of the literature. Biomolecules 12, 550. doi: 10.3390/biom12040550
Szklarczyk, D., Kirsch, R., Koutrouli, M., Nastou, K., Mehryary, F., Hachilif, R., et al. (2023). The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 51, D638–D646. doi: 10.1093/nar/gkac1000
Wang, Q., Hao, C., Yao, W., Zhu, D., Lu, H., Li, L., et al. (2020). Intestinal flora imbalance affects bile acid metabolism and is associated with gallstone formation. BMC Gastroenterol. 20, 59. doi: 10.1186/s12876-020-01195-1
Ye, X., Huang, D., Dong, Z., Wang, X., Ning, M., Xia, J., et al. (2022). FXR signaling-mediated bile acid metabolism is critical for alleviation of cholesterol gallstones by lactobacillus strains. Microbiol. Spectr. 10, e0051822. doi: 10.1128/spectrum.00518-22
Ye, X., Shen, S., Xu, Z., Zhuang, Q., Xu, J., Wang, J., et al. (2021). Sodium butyrate alleviates cholesterol gallstones by regulating bile acid metabolism. Eur. J. Pharmacol. 908, 174341. doi: 10.1016/j.ejphar.2021.174341
Keywords: gallstones, causal relationship, bile salt hydrolase, gut microbiota, Mendelian randomization
Citation: Liu X, Qi X, Han R, Mao T and Tian Z (2023) Gut microbiota causally affects cholelithiasis: a two-sample Mendelian randomization study. Front. Cell. Infect. Microbiol. 13:1253447. doi: 10.3389/fcimb.2023.1253447
Received: 07 July 2023; Accepted: 25 September 2023;
Published: 09 October 2023.
Edited by:
Yolanda López-Vidal, National Autonomous University of Mexico, MexicoReviewed by:
Almagul Kushugulova, Nazarbayev University, KazakhstanAkihiko Oka, Shimane University, Japan
Copyright © 2023 Liu, Qi, Han, Mao and Tian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zibin Tian, tianzb@qdu.edu.cn