Skip to main content

Causal association between gut microbiota and intrahepatic cholestasis of pregnancy: mendelian randomization study

Abstract

Background

Previous observational cohort studies have shown that the composition of the gut microbiota is related to the risk of intrahepatic cholestasis of pregnancy (ICP), although it is unclear if the association is causative. This study used Mendelian randomization (MR) to systematically examine whether the gut microbiota was causally linked to ICP.

Methods

We obtained the genome-wide association study (GWAS) summary statistics of gut microbiota and ICP from published GWASs. Maximum likelihood (ML), MR-Egger regression, weighted median, inverse variance weighted (IVW), and weighted model were used to investigate the causal association between gut microbiota and ICP. We further conducted a series of sensitivity analyses to confirm the robustness of the primary results of the MR analyses. Reverse MR analysis was performed on the bacterial taxa that were reported to be causally linked to ICP risk in forwarding MR analysis to evaluate the possibility of reverse causation.

Results

MR analysis revealed that phylum Tenericutes (OR: 1.670, 95%CI: 1.073–2.598, P = 0.023), class Bacteroidia (OR: 1.644, 95%CI: 1.031–2.622, P = 0.037), class Mollicutes (OR: 1.670, 95%CI: 1.073–2.598, P = 0.023), and order Bacteroidales (OR: 1.644, 95%CI: 1.031–2.622, P = 0.037), and were positively associated with the risk of ICP. And we identified that the relative abundance of genus Dialister (OR: 0.562, 95%CI: 0.323–0.977, P = 0.041), genus Erysipelatoclostridium (OR: 0.695, 95%CI: 0.490–0.987, P = 0.042), genus Eubacterium (brachy group) (OR: 0.661, 95%CI: 0.497–0.880, P = 0.005), genus Eubacterium (hallii group) (OR: 0.664, 95%CI: 0.451–0.977, P = 0.037), genus Holdemania (OR: 0.590, 95%CI: 0.414–0.840, P = 0.003), genus Ruminococcus (torques group) (OR: 0.448, 95%CI: 0.235–0.854, P = 0.015), and genus Veillonella (OR: 0.513, 95%CI: 0.294–0.893, P = 0.018) were related to a lower risk of ICP. Additional sensitivity analyses confirmed the robustness of the association between specific gut microbiota composition and ICP. No evidence of reverse causality from ICP to identified bacterial taxa was found in the findings of the reverse MR analyses.

Conclusions

Under MR assumptions, our findings propose new evidence of the relationship between gut microbiota and ICP risk. Our results show that the gut microbiota may be useful target of intervention for ICP.

Peer Review reports

Introduction

Intrahepatic cholestasis of pregnancy (ICP), typically presenting in the 2nd or 3rd trimester, is the most common liver disease specific to pregnancy [1]. ICP affects 0.3–5.6% of pregnant women, with marked differences by ethnicity [2]. ICP is characterized by the new-onset maternal pruritus and elevated serum bile acids concentration, and maternal symptoms and biochemical test abnormalities resolve after delivery [1, 3, 4]. Numerous studies revealed that the development of ICP is characterized by a decrease in bile flow through the liver and subsequent excretion, which would finally result in an intrahepatic accumulation of toxic bile acids [5]. ICP has impacts on both maternal and fetal health, with the effects on the fetus being the most concerning. Increased bile acid peak concentrations are associated with multiple adverse perinatal outcomes, including meconium-stained amniotic fluid, spontaneous preterm birth, stillbirth, fetal asphyxia, and neonatal unit admission [2, 6, 7]. The dangers of ICP call for more investigation into the causes of this condition.

Growing evidences have indicated the observational associations between gut microbiota and the metabolism of bile acids. Gut microbiota is involved in multi biological processes that contribute to the metabolism of bile acids. Sayin et al. demonstrated that gut microbiota played a vital role in the metabolism of deconjugated primary bile acids into secondary bile acids through a series of enzymatic reactions [8]. Gut microbiota could also regulate bile acid synthesis in the liver by alleviating farnesoid X receptor (FXR) inhibition in the ileum [8]. An animal study also discovered the probiotic could improve bile acid metabolism dysregulation in pregnant rats [9]. But the studies investigating the relationships between gut microbiota and ICP are largely of observational nature. Conventional observational studies have bias and confounding factors that make it difficult to draw firm conclusions about whether or not a particular bacterial taxon is significantly linked to the risk of ICP [10].

Mendelian randomization (MR) is an approach integrating summary data of genome-wide association studies (GWAS) to estimate the causal link between risk factors and outcomes [11]. Genetic polymorphisms that are randomly assigned at conception are used in MR design as instrumental variables, which can fill in the gaps in the evidence by reducing confounding variables [12]. MR method must conform to three important assumptions. First, the instrumental variables must be associated with the risk factors of interest. Second, instrumental variables must be independent of confounders that influence risk factors and outcomes. Third, the instrumental variables can only influence the outcome through the risk factor [13, 14].

Thus, to determine the causal effect of gut microbiota on ICP risk, we carried out a bidirectional two sample MR design for the first time. We selected genetic variants significantly associated with specific intestinal flora as instrumental variables (IVs) to improve interference for a possible influence of gut microbiota on ICP.

Methods

Data sources

Summary-level data for human gut microbiota were collected from a GWAS meta-analysis published to date for gut microbiota composition performed by the MiBioGen consortium [15, 16]. To investigate how human genetics affect the gut microbiota, the study coordinated the 16 S rRNA gene sequencing profiles and genome-wide genotypes of 18,340 individuals from 24 cohorts, the majority of whom were of European ancestry [16]. In the study, we excluded the 15 bacterial taxa without specific species names (unknown family or genus). As a result, 196 bacterial taxa (119 genera, 32 families, 20 orders, 16 classes, and 9 phyla) were included in the current study for analysis. FinnGen collaboration provided the GWAS summary data for ICP, which included 940 cases and 122,639 controls of European ancestry. A brief description of the bidirectional MR design is displayed in Fig. 1.

Fig. 1
figure 1

Study design of the bidirectional MR study of the associations between gut microbiota and ICP

Instrumental variable selection

The selection criteria for selecting optimal IVs were as follows. First, to obtain more comprehensive results, we selected a set of single nucleotide polymorphisms (SNPs) strongly associated with specific bacterial taxa (p < 1.0 × 10–5) as IVs [17]. Second, SNPs associated with each bacterial taxa were clumped to retain only independent SNPs. The linkage disequilibrium (LD) threshold for clumping was set to r2 < 0.001, with a clumping distance of 10,000 kb. LD was calculated based on the 1000 Genomes project European samples data reference panel. Third, the proxy SNPs (r2 > 0.8) would be selected to substitute the variants of interest, which were absent in the outcome summary statistic. Fourth, palindromic SNPs are removed from the chosen IVs in this MR analysis. Fifth, the strength of IVs was evaluated by calculating the F-statistic. If the corresponding F-statistic was much greater than 10, it was considered that there was small possibility of weak IVs bias [18].

Statistical analysis

We applied fixed-effect inverse variance weighted (IVW) method for the primary MR analysis, with the random-effects modelling used in the presence of potential heterogeneity among selected SNPs [19]. The heterogeneity of IVs was measured using Cochran’s Q statistics. To derive a comprehensive evaluation of the impact of risk factors on outcomes, Maximum likelihood (ML) [20], Mendelian randomization Egger (MR Egger) regression [20], weighted median [21], and weighted model [21] were also conducted to infer the causality, although these methods have less statistical power than does the IVW test.

To evaluate the robustness of the detected causal effect of gut microbiota on ICP, we further carried out a series of sensitivity analyses, including the MR-Egger regression, Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) analysis, and leave-one-out sensitivity analysis. The intercept from the MR-Egger could provide a formal test of directional pleiotropy [22]. The MR-PRESSO analysis detects the pleiotropic biases and corrects horizontal pleiotropy by removing the outliers [23]. Each instrumental SNP was removed one at a time during the leave-one-out sensitivity analysis to determine whether any possible SNPs with strong influences existed. In addition, reverse MR analyses were carried out to infer whether there was a reverse causal link between identified bacteria in forward MR analysis and ICP.

All statistical analyses were performed in R (version 4.1.3) with TwosampleMR [24, 25] and MR-PRESSO packages [23].

Results

Causal effect of gut microbiota on ICP

According to the selection criteria of IVs, a total of 2108 SNPs were used as IVs for 195 intestinal floras. The details about the used genetic variants in the MR analysis were shown in the Table S1. The F-statistic of the selected SNPs were all larger than 10, indicating the research was not susceptible to weak instrumental variables bias.

As illustrated in Table 1; Figs. 2 and 11 bacterial taxa, including phylum Tenericutes, class Bacteroidia, class Mollicutes, order Bacteroidales, genus Dialister, genus Erysipelatoclostridium, genus Eubacterium (brachy group), genus Eubacterium (hallii group), genus Holdemania, genus Ruminococcus (torques group), and genus Veillonella were associated with the risk of ICP in at least one MR method. Specifically, the results of IVW analysis showed that phylum Tenericutes (OR: 1.670, 95%CI: 1.073–2.598, P = 0.023), class Bacteroidia (OR: 1.644, 95%CI: 1.031–2.622, P = 0.037), class Mollicutes (OR: 1.670, 95%CI: 1.073–2.598, P = 0.023), order Bacteroidales (OR: 1.644, 95%CI: 1.031–2.622, P = 0.037), and were positively associated with the risk of ICP. The relative abundance of genus Dialister (OR: 0.562, 95%CI: 0.323–0.977, P = 0.041), genus Erysipelatoclostridium (OR: 0.695, 95%CI: 0.490–0.987, P = 0.042), genus Eubacterium (brachy group) (OR: 0.661, 95%CI: 0.497–0.880, P = 0.005), genus Eubacterium (hallii group) (OR: 0.664, 95%CI: 0.451–0.977, P = 0.037), genus Holdemania (OR: 0.590, 95%CI: 0.414–0.840, P = 0.003), genus Ruminococcus (torques group) (OR: 0.448, 95%CI: 0.235–0.854, P = 0.015), and genus Veillonella (OR: 0.513, 95%CI: 0.294–0.893, P = 0.018) were negatively related to the risk of ICP.

Table 1 MR estimates of causal effect of gut microbiota on ICP
Fig. 2
figure 2

Scatter plots for the causal effect of gut microbiota on ICP. (A) Phylum Tenericutes(B) Class Bacteroidia (C) Class Mollicutes(D) Order Bacteroidales(E) Genus Dialister(F) Genus Erysipelatoclostridium(G) Genus Eubacterium (brachy group) (H) Genus Eubacterium (hallii group) (I) Genus Holdemania(J) Genus Ruminococcus (torques group) (K) Genus Veillonella. Plot showing the effect sizes of the SNP effects on ICP (y-axes) and the SNP effects on bacterial traits (x-axes)

Sensitivity analysis

Cochran’s IVW Q test showed that there was no substantial heterogeneity among the chosen IVs (Table S2). There was no evidence of horizontal pleiotropy according to the results of the MR-Egger intercept analysis (Table S3). No pleiotropic SNP was found in the analysis of class Bacteroidia (P = 0.529), order Bacteroidales (P = 0.541), genus Dialister (P = 0.183), genus Erysipelatoclostridium (P = 0.804), genus Eubacterium (brachy group) (P = 0.523), genus Eubacterium (hallii group) (P = 0.715), genus Holdemania (P = 0.958), genus Ruminococcus (torques group) (P = 0.447), and genus Veillonella (P = 0.789) assessed by MRPRESSO analysis. We identified, rs10108398 as outlier for the association between phylum Tenericutes and ICP risk. Besides, the results of MR-PRESSO global test showed that there was significant horizontal pleiotropy between the IVs associated with class Mollicutes (P = 0.002), and we identified rs10108398 an outlier SNP. After removal of these pleiotropic SNPs, the outlier-corrected results showed that there is no evidence of horizontal pleiotropy of the remaining IVs (Table S3). In the leave-out sensitivity analysis, we discovered that none of the risk estimates between particular bacterial taxa and risk of ICP were caused by a single SNP (Fig. 3).

Fig. 3
figure 3

Leave-one-out plots for the causal effect of gut microbiota on ICP. (A) Phylum Tenericutes(B) Class Bacteroidia (C) Class Mollicutes(D) Order Bacteroidales(E) Genus Dialister(F) Genus Erysipelatoclostridium(G) Genus Eubacterium (brachy group) (H) Genus Eubacterium (hallii group) (I) Genus Holdemania(J) Genus Ruminococcus (torques group) (K) Genus Veillonella

Reverse MR analysis

We performed a reverse MR analysis to infer whether there was a causal link between ICP and relative abundance of 11 bacterial taxa. We selected the SNPs (P < 5 × 10− 8) significantly associated with the risk of ICP as IVs (Table S4). The results indicated there was no reverse causal link between ICP and identified bacterial features (Table S5). Q statistics of the IVW test demonstrated no significant heterogeneity across the included IVs (Table S6). No substantial horizontal pleiotropy was identified in the MR-Egger intercept analysis and MR-PRESSO global test (Table S7). In the leave-out sensitivity analysis, strongly affecting SNPs were identified in the IVs of genus Erysipelatoclostridium (rs17209837), genus Eubacterium (hallii_group) (rs17209837), and genus Veillonella (rs1260326). After excluding the one, substantially influential IV, the results were still valid (Fig. S1).

Discussion

To the best of our knowledge, this study is the first to use published available GWAS summary statistics to undertake a bidirectional two-sample MR analysis to investigate the causal associations between gut microbiota and ICP risk. We found that ICP risk was associated with a higher abundance of the phylum Tenericutes, class Bacteroidia, class Mollicutes, and order Bacteroidale, and. By contrast, we observed that higher levels of the genus Dialister, genus Erysipelatoclostridium, genus Eubacterium (brachy group), genus Eubacterium (hallii group), genus Holdemania, genus Ruminococcus (torques group), and genus Veillonella were related to a lower risk of ICP.

Although the precise biological mechanism underlying ICP’s occurrence is still unknown, earlier observational studies have suggested that gut microbiota played a significant part in the disease’s etiology. A case-control study performed by Li et al. showed that the butyrate-producing bacteria including genus Eubacterium (hallii group), were depleted in ICP patients, which was in line with our results [26]. According to Zhan et al., the severe ICP group displayed considerable gut microbiota dysbiosis and had increased microbial gene functions for propanoate metabolism and the manufacture of unsaturated fatty acids [27]. However, there was no consensus on whether there was a causal effect of gut microbiota on the occurrence of ICP. In this study, a two-sample MR analysis successfully identified that phylum Tenericutes, class Bacteroidia, and order Bacteroidale were positively associated with ICP risk. The phylum Bacteroidetes possess the enzyme bile salt hydrolase which could deconjugate bile acids, reduce ileal bile acid uptake, decrease the induction of FXR, and reduce the level of fibroblast growth factor 19/15 (FGF19/15), thereby increasing synthesis of bile acid in the liver [28, 29]. Tang et al. found that a notable difference in microbial profiles between ICP cases and healthy controls, and the microbiomes of patients with ICP were primarily characterized by Bacteriodes fragilis, which was consistent with our findings [30]. The phylum Tenericutes are positively associated with pro-inflammatory factors IL-6, TNF-α, and IL-17 A, which were significantly increased in ICP patients [31, 32]. Furthermore, we proved some bacterial taxa were negatively associated with ICP risk. Genus Dialister [33], genus Eryipelatoclostridium [34], genus Eubacterium [35], genus Ruminococcus [36] produce butyrate, which plays an important role in energy homeostasis, maintenance of the gut barrier functions, immunomodulation, and anti-inflammation [37]. The genus Eubacterium could also regulate the expression of several enzymes during the metabolism of bile acids, including, 7α-hydroxylase (Cyp7a1), oxysterol 7α-hydroxylase (Cyp7b1), and sterol 27-hydroxylase (Cyp27a1), and to determine the amount of bile acids produced [38]. We also discovered that class Mollicutes, genus Holdemania, and genus Veillonella were associated with ICP. However, there have not been many prior investigations on how these bacteria affect ICP etiology. A theoretical foundation for the use of probiotics to treat ICP should be provided by investigations on the molecular mechanisms underlying these gut microbiota compositions in ICP.

This research has several strengths. First, MR analysis is unlikely to be influenced by the interference of confounding factors and reverse causality compared with traditional observational study design. In addition, we collected the genetic variants from the largest available GWAS meta-analysis for human gut microbiota composition, ensuring the strength of genetic instruments in this MR analysis. To ensure that our results are not the product of pleiotropic effects, we also employed the MR-Egger intercept analysis and MR-PRESSO method to examine for any potential horizontal pleiotropy on the observed causal correlations.

Still, our analysis has several potential limitations. First, the present study included only participants of European ancestry, which minimize the risk of confounding due to population admixture but may limit the generalizability of our findings to different populations. Second, the two-sample MR design only tested the linear effect of the relative abundance of gut microbiota composition on ICP risk in the general population. Third, The GWAS summary statistics for gut microbiota were not restricted to the female population [16]. Although the genetic variants located on the sex chromosomes were excluded, as well as sex was adjusted in the analysis [16], the potential bias due to sex could not be totally excluded. Fourth, the precise biological mechanisms of the influence of specific intestinal flora on the pathogenesis of ICP were still unclear. To better understand the impact of gut microbiota on ICP, additional clinical and functional research was required.

Conclusions

In summary, we comprehensively assess the potential causal relationship between gut microbiota and ICP. This two-sample MR study provide precise evidence that the relative abundance of several gut microbiota was causally associated with ICP. To better understand how probiotics affect ICP and its precise biological mechanism, additional research that incorporates the findings of well-designed randomized controlled trials is necessary. Additionally, although reverse MR analysis did not support the causal effect of ICP on gut microbiota, it cannot be ruled out that ICP may influence intestinal microorganisms, which needs to be verified by further research.

Data Availability

The datasets analyzed in this current study can be accessed through the following links: microbiota, https://mibiogen.gcc.rug.nl/; intrahepatic cholestasis of pregnancy, https://r5.finngen.fi/pheno/O15_ICP.

Abbreviations

GWAS:

Genome-wide association study

MR:

Mendelian randomization

SNP:

Single nucleotide polymorphism

IV:

Instrumental variable

IVW:

Inverse variance weighted

ML:

Maximum likelihood

MR-PRESSO:

Mendelian randomization pleiotropy residual sum and outlier

LD:

Linkage disequilibrium

ICP:

Intrahepatic cholestasis of pregnancy

OR:

Odds ration

CI:

Confidence interval

SCFA:

Short-chain fatty acid

References

  1. Williamson C, Geenes V. Intrahepatic cholestasis of pregnancy. Obstet Gynecol. 2014;124(1):120–33.

    Article  PubMed  Google Scholar 

  2. Ovadia C, Sajous J, Seed PT, et al. Ursodeoxycholic acid in intrahepatic cholestasis of pregnancy: a systematic review and individual participant data meta-analysis. Lancet Gastroenterol Hepatol. 2021;6(7):547–58.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Axelsen SM, Kampmann U, Koefoed AS, McIntyre D, Ovesen PG, Fuglsang J. Intrahepatic cholestasis of pregnancy: Association with glycaemic control in gestational diabetes. Diabet Med. 2021;38(8):e14574.

    Article  CAS  PubMed  Google Scholar 

  4. Roediger R, Fleckenstein J. Intrahepatic cholestasis of pregnancy: natural history and current management. Semin Liver Dis. 2021;41(1):103–8.

    Article  PubMed  Google Scholar 

  5. Shao Y, Chen S, Li H, Tang Q, Xu D. Maternal bile acid profile and subtype analysis of intrahepatic cholestasis of pregnancy. Orphanet J Rare Dis. 2021;16(1):259.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Chappell LC, Bell JL, Smith A, et al. Ursodeoxycholic acid versus placebo in women with intrahepatic cholestasis of pregnancy (PITCHES): a randomised controlled trial. Lancet. 2019;394(10201):849–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ovadia C, Seed PT, Sklavounos A, et al. Association of adverse perinatal outcomes of intrahepatic cholestasis of pregnancy with biochemical markers: results of aggregate and individual patient data meta-analyses. Lancet. 2019;393(10174):899–909.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Sayin SI, Wahlström A, Felin J, et al. Gut microbiota regulates bile acid metabolism by reducing the levels of tauro-beta-muricholic acid, a naturally occurring FXR antagonist. Cell Metab. 2013;17(2):225–35.

    Article  CAS  PubMed  Google Scholar 

  9. Geenes V, Lövgren-Sandblom A, Benthin L, et al. The reversed feto-maternal bile acid gradient in intrahepatic cholestasis of pregnancy is corrected by ursodeoxycholic acid. PLoS ONE. 2014;9(1):e83828.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Hartwig FP, Borges MC, Horta BL, Bowden J, Davey Smith G. Inflammatory biomarkers and risk of Schizophrenia: a 2-Sample mendelian randomization study. JAMA Psychiatry. 2017;74(12):1226–33.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Xiang K, Wang P, Xu Z, et al. Causal Effects of Gut Microbiome on systemic lupus erythematosus: a two-sample mendelian randomization study. Front Immunol. 2021;12:667097.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Chan II, Kwok MK, Schooling CM. Blood pressure and risk of cancer: a mendelian randomization study. BMC Cancer. 2021;21(1):1338.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hernán MA, Robins JM. Instruments for causal inference: an epidemiologist’s dream. Epidemiology. 2006;17(4):360–72.

    Article  PubMed  Google Scholar 

  14. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, EPIC- InterAct Consortium. Using published data in mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Wang J, Kurilshikov A, Radjabzadeh D, et al. Meta-analysis of human genome-microbiome association studies: the MiBioGen consortium initiative. Microbiome. 2018;6(1):101.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Kurilshikov A, Medina-Gomez C, Bacigalupe R, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet. 2021;53(2):156–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Sanna S, van Zuydam NR, Mahajan A, et al. Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet. 2019;51(4):600–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Burgess S, Thompson SG, CRP CHD Genetics Collaboration. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.

    Article  PubMed  Google Scholar 

  19. Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in mendelian randomization: comparison of allele score and summarized data methods. Stat Med. 2016;35(11):1880–906.

    Article  PubMed  Google Scholar 

  20. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid Instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11):e1007081.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Hemani G, Zheng J, Elsworth B et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7.

  26. Li GH, Huang SJ, Li X, Liu XS, Du QL. Response of gut microbiota to serum metabolome changes in intrahepatic cholestasis of pregnant patients. World J Gastroenterol. 2020;26(46):7338–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhan Q, Qi X, Weng R, et al. Alterations of the human gut microbiota in intrahepatic cholestasis of pregnancy. Front Cell Infect Microbiol. 2021;11:635680.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ovadia C, Perdones-Montero A, Spagou K, et al. Enhanced microbial bile acid deconjugation and impaired ileal uptake in pregnancy repress intestinal regulation of bile acid synthesis. Hepatology. 2019;70(1):276–93.

    Article  CAS  PubMed  Google Scholar 

  29. Cai J, Sun L, Gonzalez FJ. Gut microbiota-derived bile acids in intestinal immunity, inflammation, and tumorigenesis. Cell Host Microbe. 2022;30(3):289–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tang B, Tang L, Li S, et al. Gut microbiota alters host bile acid metabolism to contribute to intrahepatic cholestasis of pregnancy. Nat Commun. 2023;14(1):1305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Larson SP, Kovilam O, Agrawal DK. Immunological basis in the pathogenesis of intrahepatic cholestasis of pregnancy. Expert Rev Clin Immunol. 2016;12(1):39–48.

    Article  CAS  PubMed  Google Scholar 

  32. Li K, Zhang L, Xue J, et al. Dietary inulin alleviates diverse stages of type 2 diabetes mellitus via anti-inflammation and modulating gut microbiota in db/db mice. Food Funct. 2019;10(4):1915–27.

    Article  CAS  PubMed  Google Scholar 

  33. Pan L, Ai X, Fu T, et al. In vitro fermentation of hyaluronan by human gut microbiota: changes in microbiota community and potential degradation mechanism. Carbohydr Polym. 2021;269:118313.

    Article  CAS  PubMed  Google Scholar 

  34. Beisner J, Gonzalez-Granda A, Basrai M, Damms-Machado A, Bischoff SC. Fructose-Induced Intestinal Microbiota Shift following two types of short-term high-fructose Dietary Phases. Nutrients. 2020;12(11).

  35. Mukherjee A, Lordan C, Ross RP, Cotter PD. Gut microbes from the phylogenetically diverse genus Eubacterium and their various contributions to gut health. Gut Microbes. 2020;12(1):1802866.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Shi J, Yang Y, Xu W et al. Sex-Specific Associations between Gut Microbiome and non-alcoholic fatty liver disease among urban chinese adults. Microorganisms. 2021;9(10).

  37. Rivière A, Selak M, Lantin D, Leroy F, De Vuyst L. Bifidobacteria and Butyrate-Producing Colon bacteria: importance and strategies for their stimulation in the human gut. Front Microbiol. 2016;7:979.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ren S, Zhou Y, Xuan R. Research progress in the role of gut microbiota and its metabolites in intrahepatic cholestasis of pregnancy. Expert Rev Gastroenterol Hepatol. 2021;15(12):1361–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors want to acknowledge the participants and investigators of the FinnGen and MiBioGen consortium.

Funding

This work was supported by the Natural Science Foundation of Liaoning Province. (Grant no.2022-MS-13)

Author information

Authors and Affiliations

Authors

Contributions

Chuang Li and Na Li interpreted and analyzed the data, and drafted the manuscript. Shaowei Yin and Caixia Liu concepted and designed the study and revised the manuscript. All authors contributed to and approved the final manuscript.

Corresponding author

Correspondence to Shaowei Yin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

All authors have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, C., Li, N., Liu, C. et al. Causal association between gut microbiota and intrahepatic cholestasis of pregnancy: mendelian randomization study. BMC Pregnancy Childbirth 23, 568 (2023). https://doi.org/10.1186/s12884-023-05889-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12884-023-05889-8

Keywords