Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 09 October 2023
Sec. Intestinal Microbiome

Gut microbiota causally affects cholelithiasis: a two-sample Mendelian randomization study

Xin LiuXin LiuXingsi QiXingsi QiRongshuang HanRongshuang HanTao MaoTao MaoZibin Tian*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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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
www.frontiersin.org

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Birney, E. (2022). Mendelian randomization. Cold Spring Harb. Perspect. Med. 12, a041302. doi: 10.1101/cshperspect.a041302

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Holmes, M. V. (2019). Meta-analysis and Mendelian randomization: A review. Res. Synth. Methods 10, 486–496. doi: 10.1002/jrsm.1346

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Emdin, C. A., Khera, A. V., Kathiresan, S. (2017). Mendelian randomization. JAMA 318, 1925–1926. doi: 10.1001/jama.2017.17219

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukui, H. (2021). Leaky gut and gut-liver axis in liver cirrhosis: clinical studies update. Gut Liver 15, 666–676. doi: 10.5009/gnl20032

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Grigor’eva, I. N., Romanova, T. I. (2020). Gallstone disease and microbiome. Microorganisms 8, 835. doi: 10.3390/microorganisms8060835

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

PubMed Abstract | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zmora, N., Suez, J., Elinav, E. (2019). You are what you eat: diet, health and the gut microbiota. Nat. Rev. Gastroenterol. Hepatol. 16, 35–56. doi: 10.1038/s41575-018-0061-2

PubMed Abstract | CrossRef Full Text | Google Scholar

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, Mexico

Reviewed by:

Almagul Kushugulova, Nazarbayev University, Kazakhstan
Akihiko 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

Disclaimer: 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.