- 1Guang’anmen Hospital, China Academy of Chinese Medical Sciences, Beijing, China
- 2Graduate School, Beijing University of Chinese Medicine, Beijing, China
- 3Peking Union Medical College Hospital, Chinese Academy of Medical Science, Beijing, China
Background: The association of gut microbiota (GM) and chronic kidney disease (CKD), and the relevancy of GM and chronic systemic inflammation in CKD, were revealed on the basis of researches on gut–kidney axis in previous studies. However, their causal relationships are still unclear.
Objective: To uncover the causal relationships between GM and CKD, as well as all known GM from eligible statistics and chronic systemic inflammation in CKD, we performed two-sample Mendelian randomization (MR) analysis.
Materials and methods: We acquired the latest and most comprehensive summary statistics of genome-wide association study (GWAS) from the published materials of GWAS involving GM, CKD, estimated glomerular filtration rate (eGFR), c-reactive protein (CRP) and urine albumin creatine ratio (UACR). Subsequently, two-sample MR analysis using the inverse-variance weighted (IVW) method was used to determine the causality of exposure and outcome. Based on it, additional analysis and sensitivity analysis verified the significant results, and the possibility of reverse causality was also assessed by reverse MR analysis during this study.
Results: At the locus-wide significance threshold, IVW method and additional analysis suggested that the protective factors for CKD included family Lachnospiraceae (P=0.049), genus Eubacterium eligens group (P=0.002), genus Intestinimonas (P=0.009), genus Streptococcu (P=0.003) and order Desulfovibrionales (P=0.001). Simultaneously, results showed that genus LachnospiraceaeUCG010 (P=0.029) was a risk factor for CKD. Higher abundance of genus Desulfovibrio (P=0.048) was correlated with higher eGFR; higher abundance of genus Parasutterella (P=0.018) was correlated with higher UACR; higher abundance of class Negativicutes (P=0.003), genus Eisenbergiella (P=0.021), order Selenomonadales (P=0.003) were correlated with higher CRP levels; higher abundance of class Mollicutes (0.024), family Prevotellaceae (P=0.030), phylum Tenericutes (P=0.024) were correlated with lower levels of CRP. No significant pleiotropy or heterogeneity was found in the results of sensitivity analysis, and no significant causality was found in reverse MR analysis.
Conclusion: This study highlighted associations within gut-kidney axis, and the causal relationships between GM and CKD, as well as GM and chronic systemic inflammation in CKD were also revealed. Meanwhile, we expanded specific causal gut microbiota through comprehensive searches. With further studies for causal gut microbiota, they may have the potential to be new biomarkers for targeted prevention of CKD and chronic systemic inflammation in CKD.
1 Introduction
Recently, a worldwide trend of increasing prevalence has appeared in chronic kidney disease (CKD) (1). The most important and commonly used CKD-classifying quantitative traits include estimated glomerular filtration rate (eGFR) and urinary albumin-to-creatinine ratio (UACR) (2, 3). Based on these traits, an epidemiological investigation revealed (4) that the prevalence of CKD is estimated to be as high as 13.4% (11.7-15.1%) worldwide. What’s more, CKD significantly increases the mortality of patients, and heavily burdens the financial and healthcare system (5). In addition, Patients with CKD are usually in a chronic inflammatory status, as manifested by the elevated levels of inflammatory markers in CKD patients, which are caused by multiple factors (6). To some extent, renal functions and severity of inflammation in CKD patients can be reflected by levels of c-reactive protein (CRP) or other inflammatory markers (7). Interestingly, levels of the same inflammatory markers are inconsistent within the same period of CKD patients. Besides explicit risk factors including diabetes mellitus (DM) (8), hypertension (HTN) (9), and oxidative st (10) that can influence the inflammation of CKD patients, more risk factors remain to be explored.
Targeting explicit risk factors, the only interventions available for improving renal function and chronic systemic inflammation in patients with CKD are lifestyle changes, pharmacotherapy, and optimization of dialysis conditions (11). However, these interventions are either costly or lack sufficient clinical evidence to support their therapeutic effects. Hence, it is essential to identify potential risk factors for CKD and chronic systemic inflammation in CKD, and develop new preventive and therapeutic measures for them.
Gut microbiota (GM) is a crucial regulator of human health (12). As a critical regulator of human health, it has gradually received attention recently for being considered one of the new potential risk factors of CKD and chronic systemic inflammation in CKD (13, 14). Meijers et al. (15), in 2011, proposed the concept known as gut-kidney axis, revealing that GM might interact with CKD and chronic systemic inflammation in CKD through it. Lau Wei Ling et al. (16) further indicated that GM of CKD patients appeared to be ecologically imbalanced, and the degree of GM imbalance directly influenced the severity of inflammation in the dialysis population. Some studies asserted that interventions on GM could reduce the levels of inflammatory factors and delay disease progression in CKD patients (17). Conversely, some studies suggested that altering types and abundance of GM could not improve renal function and chronic systemic inflammation in CKD patients (18). Nevertheless, their findings were limited by confounding bias, small sample size and reverse causality. Therefore, the association effects they showed were not equivalent to causality. The causal relationships between GM and CKD, as well as GM and chronic systemic inflammation in CKD are still unclear.
Mendelian randomization (MR), as a novel genetic statistical method, can be an effective alternative to traditional epidemiological study approaches. It uses genetic variants strongly associated with exposure factors as instrumental variables (IVs) to statistically evaluate the causality of exposure and outcome (19). The advantage of MR over traditional epidemiological study approaches is that its random assignment method is determined by the DNA genotype. Thus, the influence of external factors on the robustness of causality can be limited to the greatest extent possible (20). In this research, the latest genome-wide association study (GWAS) statistics were used in two-sample MR analysis, and four sets of causal relationships were analyzed at the genetic level: GM and CKD, GM and eGFR, GM and UACR, GM and CRP. Also, we determined the specific causal GM in these four sets. The results supplemented gaps in existing research and might provide new ideas for improving and enriching treatment measures.
2 Materials and methods
2.1 Study design
In this study, four sets of causal relationships were assessed by two-sample MR method: GM and CKD, GM and eGFR, GM and UACR, GM and CRP. The overall study flow is displayed in Figure 1. To ensure the reliability of MR results, three basic assumptions (21) must be conformed: (I) IVs significantly correlated with exposure are used in the analysis. (II) The IVs are independent of confounding factors affecting exposure and outcome. (III) The IVs are not horizontally pleiotropic, that is, the IVs have effects on outcome only through exposure.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1287698/fimmu-14-1287698-HTML/image_m/fimmu-14-1287698-g001.jpg)
Figure 1 The study design and major assumptions of the Two-sample MR study. GM, gut microbiota; LD, linkage disequilibrium; SNP, single nucleotide polymorphism; IV: instruments variable; CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; UACR, urine albumin creatine ratio; CRP, c-reactive protein; MR-PRESSO, Mendelian Randomization Pleiotropy RESidual Sum and Outlier.
2.2 Data sources of instruments variables
The largest GWAS meta-analysis of GM to date is from the MiBioGen consortium, containing 25 cohorts totaling 18,340 participants of European ancestry from 11 countries (22). From 211 bacterial taxa in total, this study ultimately identified 122,110 variant sites at 5 levels: phylum, class, order, family, and genus. To ensure the accuracy of the data, our study excluded 15 bacterial taxa of unknown family or genus, leaving 196.
Accurate measurement of renal function is difficult to achieve, so the use of biomarkers is necessary (23). Thus, we selected four sets of IVs to represent the renal function and chronic systemic inflammation. Among them, CKD is considered a chronic disease associated with impairment of renal function. eGFR is an important indicator of renal filtration function (24). UACR, a reflection of the degree of urinary protein, is the main clinical diagnostic criteria for CKD (25). CRP is one of the inflammatory markers measuring the severity of the inflammation status (26). The CKD, eGFR and UACR related summarized level GWAS data were collected from the CKDGen consortium’s meta-analysis of GWAS from participants of European ancestry. Wuttke et al. have reported CKDGen consortium in detail (27), so we did not elaborate here. With 480698 samples from 23 European ancestry cohorts, the GWAS meta-analysis of CKD comprised 41395 samples in the trial group and 439303 samples in the control group. GWAS meta-analysis of eGFR comprised 567460 European ancestry samples from 54 cohorts. The data of UCAR included 54 GWAS summary statistics of 564,257 participants. GWAS summary statistics related to CRP were obtained from The MRC IEU OpenGWAS data infrastructure (28). The dataset name is c-reactive protein (ID: ieu-b-4764) and it contains 8036590 single nucleotide polymorphism (SNP) in a sample size of 61308.
2.3 Selection of instruments variables
Quality checks of SNPs were performed to obtain eligible IVs: (I) GM related SNPs must reach a threshold (P < 5×10-8) with genome-wide significance. To obtain more comprehensive results, another set of SNPs reaching locus-wide significance level (P < 1×10-5) as IV was selected. (II) No linkage disequilibrium (LD) existed among GM-associated IVs, and a clumping process (r2 < 0.001, clumping distance = 10000kb) was performed on the screened SNPs to retain independent ones. (III) Moreover, for SNPs that were not available in GWASs of the outcome, we used the LD proxy search on the online platform (https://snipa.helmholtz-muenchen.de/snipa3/index.php/) to replace them with the proxy SNPs identified in high-LD (r2 > 0.8) or discard them if the proxies were not available (IV) The effects of GM-associated IVs on both exposure and outcome corresponded to the same alleles, so the palindromic SNPs were removed. (V) To avoid SNPs associated with potential risk factors for outcome, the PhenoScanner V2 website was used to retrieve these SNPs and exclude those associated with potential confounders or risk factors. (VI) To avoid bias caused by weak IVs, we calculated the intensity of IVs using F=R 2 (n-k-1)/k (1-R 2) (29, 30). R 2 represents the exposure variance explained by the selected SNPs, n is the sample size, and k represent the number of included instrumental variables. We excluded weak IVs with F < 10 (31).
2.4 Ethics statement
In this study, all the summary-level data were published available de-identified ones, which were authorized by the Ethical Standards Committee. No independent ethical approval was necessary during the research.
2.5 Mendelian randomization analysis
When no horizontal pleiotropy was available, the IVW method was used in this study as the primary method for inferring 4 sets of causal relationships: GM and CKD, GM and eGFR, GM and UACR, GM and CRP (32). To detect the presence of heterogeneity, we performed Cochran’ Q test. If there was significant heterogeneity (P < 0.05), a random-effects IVW model was adopted, otherwise a fixed-effects IVW model was applied (33). Moreover, to obtain more robust results under broader conditions, the weighted median (WM) approach and MR-Egger method were adopted to complement the IVW method. The additional methods need to satisfy the respective model assumptions: the WM approach assumes that at least half of the SNPs are free of pleiotropy (34). If the number of SNPs possessing pleiotropy exceeds 50%, the MR-Egger inference is still robust (35). In this study, the causality of exposure and outcome was considered to exist if the results of the main MR analysis reached a nominal significance (P < 0.05). The result would be regarded as significant and stable if it was supported by one or more additional methods simultaneously (36), and we would provide a focused discussion on such results.
To avoid the interference of pleiotropy to MR hypothesis, sensitivity analysis on the study results was performed: MR-Egger regression was used to estimate the potential horizontal pleiotropy of the included SNPs, and the results were considered to have horizontal pleiotropy if P < 0.05 (35). Considering the lower precision and statistical efficacy of MR-Egger regression, we used Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) to examine any deviations to reflect pleiotropy bias and to give the causal effect of excluding outliers (37). The results of the sensitivity analysis for all IVs are presented in the Supplementary Tables 1, 2. In addition, we conducted leave-one-out sensitivity analysis on significant results, in order to determine whether the significant causal association of the MR analysis was caused by a single IV (38).
2.6 Reverse-direction Mendelian randomization analysis
In MR analysis, additional reverse-direction MR analysis on stable and significant results was performed to test whether gene-predicted CKD, eGFR, UACR, and CRP would be causal to GM. The steps of reverse-direction MR analysis were the same as those of MR analysis.
3 Results
3.1 Selection of IVs related GM
Through a series of quality control measures, at the genome-wide significance level (P < 5×10-8), 12 SNPs, 11 SNPs, 12 SNPs and 12 SNPs were used for genetic prediction of CKD, eGFR, UACR and CRP, respectively. At the locus-wide significance level (P < 1×10-5), 2195 SNPs, 2166 SNPs, 2157 SNPs and 2129 SNPs were used for gene prediction of CKD, eGFR, UACR and CRP, respectively (Supplementary Tables 3, 4). And the F-statistic values of the SNPs were all met (Supplementary Tables 3, 4) (39). Based on the locus-wide significance level, we identified 16, 17, 4, and 10 bacterial taxa causally associated with CKD, eGFR, UACR, and CRP in the primary MR analysis, respectively, as detailed in Supplementary Table 5 and Figures 2, 3. After additional MR analysis and sensitivity analysis, only 6, 1, 1 and 6 bacterial taxa remained robust to the results of CKD, eGFR, UACR, and CRP, respectively (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1287698/fimmu-14-1287698-HTML/image_m/fimmu-14-1287698-g002.jpg)
Figure 2 Results of MR study and sensitivity analysis between GM and CKD, GM and eGFR, GM and UACR, GM and CRP (locus-wide significance, P<1×10-5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1287698/fimmu-14-1287698-HTML/image_m/fimmu-14-1287698-t001.jpg)
Table 1 Significant MR analysis results of causal links between GM and CKD, eGFR, UACR, CRP (P<1×10-5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1287698/fimmu-14-1287698-HTML/image_m/fimmu-14-1287698-g003.jpg)
Figure 3 MR results of GM taxa with a significant causal relationships to CKD, EGFR, UACR and CRP (locus-wide significance, P<1×10-5).
3.2 Locus-wide significance threshold P<1×10-5
3.2.1 CKD
Primary MR analysis showed that 16 bacterial taxa were associated with CKD risk (Supplementary Table 5; Figure 2). However, only 6 bacterial taxa remained robust in subsequent additional and sensitivity analysis (Table 1). Precisely, higher genetically predicted Lachnospiraceae, Eubacteriumeligens group, Intestinimonas, Streptococcus, and Desulfovibrionales were related to a lower risk of CKD [odds ratio (OR): 0.927, 95% confidence interval (CI): 0.79-0.98, P=0.049 for Lachnospiraceae; and OR: 0.830, 95% CI: 0.74-0.93, P=0.002 for Eubacteriumeligens group, and OR: 0.924, 95% CI: 0.87-0.98, P=0.009 for Intestinimonas, and OR: 0.892, 95% CI: 0.82-0.96, P=0.003 for Streptococcus, and OR: 0.873, 95% CI: 0.81-0.94, P=0.001 for Desulfovibrionales]. Conversely, higher genetic prediction of LachnospiraceaeUCG010 was related to a higher risk of CKD [OR: 1.096, 95% CI: 1.01-1.19, P=0.029].
3.2.2 eGFR
Seventeen bacterial taxa associated with eGFR were identified from primary MR analysis (Supplementary Table 5; Figure 2). After additional and sensitivity analysis, only the results for Desulfovibrio remained robust, as detailed in Table 1. Results from the IVW method showed that higher abundance of Desulfovibrio was related to higher eGFR [beta: 0.003, 95% CI: 2.37×10-5-6.48×10-3, P=0.048], and this result was also supported by the WM approach [beta: 0.005, 95% CI: 5.79×10-4-8.51×10-3, P=0.025], suggesting a protective effect of Desulfovibrio on eGFR decline.
3.2.3 UACR
Primary MR analysis confirmed that 4 bacterial taxa were associated with UACR (Supplementary Table 5; Figure 2). Sensitivity analysis and additional analysis only supported the results for Parasutterella (Table 1). IVW results showed that higher abundance of Parasutterella was a risk factor for proteinuria and was related to higher levels of UACR [beta: 0.021, 95% CI: 0.01-0.04, P=0.018].
3.2.4 CRP
Ten bacterial taxa were proved to be causally related to CRP in the primary MR analysis (Supplementary Table 5; Figure 2). However, sensitivity and additional analysis supported the results for only 6 bacterial taxa, as shown in Table 1. MR analysis of IVW showed that genetically predicted Mollicutes [beta: -0.079, 95% CI: -0.15–0.01, P=0.024], Prevotellaceae [beta: -0.059, 95% CI: -0.11–0.01, P=0.030] and Tenericutes [beta: -0.079, 95% CI: -0.15–0.01, P=0.024] were associated with lower CRP levels and were protective factors against inflammation. Conversely, genetically predicted Negativicutes [beta: 0.098, 95% CI: 0.03–0.16, P=0.003], Eisenbergiella [beta: 0.047, 95% CI: 0.01–0.09, P=0.021], and Selenomonadales [beta: 0.098, 95%CI: 0.03–0.16, P=0.003] were associated with higher CRP levels (Table 1).
As shown in Supplementary Table 6, the intercepts of the MR-Egger regressions all did not deviate from the null, suggesting no evidence for horizontal pleiotropy (P > 0.05 for all intercepts). We further validated the MR-Egger regression results using MR-PRESSO, and no evidence of outliers was found (Supplementary Table 5). Leave-one-out analysis found no outliers as in MR-PRESSO (Supplementary Figures 1A–N). Without heterogeneity and pleiotropy, outcomes of the above MR analysis were reliable and robust.
3.3 Genome-wide significance threshold P<5×10-8
At this stage, we analyzed gut microbiota at five levels (phylum, class, order, family, and genus) based on SNPs available in the gut microbiome GWAS summary data. The results showed no significant causal relationship with CKD, eGFR, UACR or CRP when MR analysis was conducted on each taxonomic level (Table 2). In sensitivity analysis, no intercepts of MR-Egger regression were diverged from the null value, suggesting the absence of evidence for horizontal pleiotropy (P > 0.05). MR-PRESSO results showed outliers in sensitivity analysis of UACR at the genus level, then we removed rs35866622, rs602075 and re-analyzed the results, the results did not change significantly (p=0.1095194, beta=0.0125, OR=1.0126). The results of MR analysis could not prove a significant causal relationship between GM and CKD, eGFR, UACR or CRP, because there were too few eligible IVs (Supplementary Table 2).
3.4 Reverse-direction Mendelian randomization analysis
Finally, reverse MR analysis was performed for the 14 bacterial taxa presenting significant results with CKD, eGFR, UACR and CRP. After additional analysis and sensitivity analysis, no significant and stable results were found. The specific information is presented in Supplementary Table 7.
4 Discussion
This research was not the first one exploring the cause-and-effect relationship between GM and CKD. However, the innovative points of this article remedied the shortcomings of the past studies: (I) GM data used here was the latest large GWAS dataset, which contains 196 bacterial taxa at 5 levels, consisting of 9 phyla, 16 classes, 20 orders, 32 families, 119 genera, in contrast, previous studies only estimated the effect of 8 microbiota genera on CKD (40). The bacterial taxa we studied in this research included but were not limited to those analyzed in previous studies. (II) The selection of IVs was more rigorous, and two thresholds were used to select independent SNPs. In contrast, previous studies did not select SNPs strongly associated with exposure and did not perform the aggregation process, which might seriously affect the reliability and robustness of MR analysis results. (III) Chronic systemic inflammation in CKD is a strong predictor of CKD prognosis, and UACR is sensitive and specific to the overall renal injury profile. Hence, this study revealed not only the causal relationship between GM and CKD, eGFR, but also GM and UACR, CRP. In conclusion, our findings expanded the bacterial taxa associated with CKD and revealed that GM played a regulatory role in CKD and its inflammation state.
The human intestine contains more than 1,000 bacterial taxa, 1014 in total, which are involved in the metabolic processes of substances in the body and regulate immune function. Several studies (41, 42) found that the composition and abundance of GM were different between healthy individuals and CKD patients. The diversity and different abundance of bacterial taxa may contribute to the development of CKD and CKD co-morbidities such as chronic systemic inflammation. Therefore, the causal relationships between GM and CKD, as well as GM and chronic systemic inflammation in CKD, were explored in this study. After additional analysis and sensitivity analysis, 14 bacterial taxa had a significant causal relationship with our outcome at the phylum, class, order, family, and genus levels (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1287698/fimmu-14-1287698-HTML/image_m/fimmu-14-1287698-g004.jpg)
Figure 4 The present MR study reveals that GM causally influences CKD, eGFR, UACR and CRP, which supports the existence of gut-kidney axis.
Family Lachnospiraceae isa critical member of GM. Previous studies demonstrated that Lachnospiraceae took an essential part in the pathogenesis of CKD (43, 44), however, there were conflicting findings among these studies. Lai et al. (45) found that the abundance of Lachnospiraceae was increased in mice with kidney damage and in patients with CKD, which might predict that it was a risk bacterial taxon for CKD. Conversely, other observational studies (44) showed that the abundance of Lachnospiraceae in stool samples from patients with CKD was reduced. Moreover, the abundance of Lachnospiraceae was negatively correlated with serum creatinine, urea nitrogen and other indicators of renal function, and our study supported these findings. The reason may be that Lachnospiraceae is the primary producer of short-chain fatty acids (SCFAs) (43). SCFAs act as a link between GM and within-host environments. They play essential roles in correcting pH value and ensuring the integrity of the intestinal epithelium (46). SCFAs also induce renin secretion and the development of renal cell subtypes to control blood pressure, thereby delaying the progression of CKD (47).
Furthermore, a study by Hu et al. (48) found an increased abundance of genus Intestinimonas in CKD patients, suggesting that it might be a risk factor for CKD. The mechanism may be the involvement of genus Intestinimonas in the production of tryptophan-derived uremic toxins, which damage the renal vasculature and the circulatory system (49). However, these findings are inconsistent with our conclusions. It is possible that gene-gene interactions and gene-environment interactions result in the same changes in GM but cause different results. In addition, these observational studies may draw biased conclusions because of their smaller sample size (12 rats) and shorter trial duration (2 weeks). In addition to the bacteria mentioned above that are causally related to CKD, in the present study, genus Eubacterium eligens group, genus Streptococcus, and order Desulfovibrionales, all showed a positive effect on CKD, whereas genus LachnospiraceaeUCG010 indicated a negative effect on CKD.
Two-sample MR analysis also found that higher abundance of genus Desulfovibrio was positively associated with higher levels of eGFR. An observational study (50) supporting our view confirmed that the abundance of order Desulfovibrionales was low in stool samples from patients presenting with eGFR <15 mL/min/1.73 m2, implying that order Desulfovibrionales was a beneficial bacterial taxon for eGFR. Genus Desulfovibrio is a bacterial genera pattern of order Desulfovibrionales, which may be relevant to the findings of this study. However, Zhao et al. (51) suggested that the abundance of genus Desulfovibrio was higher in patients with CKD whose eGFR was significantly reduced. Although the direction of association between genus Desulfovibrio and eGFR varied considerably in different studies, the current consensus concluded that genus Desulfovibrio could metabolize renoprotective SCFAs (14). It also produces more hydrogen sulfide, causing metabolic impairment and toxicity in the kidney (52, 53), which consequently affects eGFR.
In the MR analysis of GM and proteinuria, our results showed that genus Parasutterella was positively associated with higher UACR. This finding is also supported by previous studies: Feng et al. (54) found that GM disorders such as genus Parasutterella increased in abundance was positively correlated with increased proteinuria in CKD model rats. The mechanism may involve phenyl sulfate (PS) (17), an enteric-derived uremic toxin (55) that can increase proteinuria by damaging podocytes (56). The level of PS is found to be significantly associated with UACR in microalbuminuric patients. It is the only element that predicts the progression of UACR in microalbuminuric patients over the next two years (57). A few bacterial taxa such as Parasutterella and Citrobacter contain TPL, which synthesizes tyrosine into phenol, which is further metabolized to PS (58). Therefore, GM intervention will constitute a new approach towards the prevention and treatment of proteinuria.
Researchers have long worked on the relationship between GM and chronic systemic inflammation in CKD. One of the common comorbidities of CKD is chronic systemic inflammation. To some extent, chronic inflammation can reflect the renal function status of CKD patients, and it is predictive of CKD prognosis and related comorbidities (7, 59). Previous study (60) suggested that the amount and abundance of GM producing SCFAs were increased, whereas GM producing uremic toxins were increased in CKD patients. Toxins disrupted intestinal tight junctions, penetrated from the damaged intestinal mucosa into the blood, and produced systemic chronic inflammation (60). The regulation of imbalanced GM has become an increasingly popular topic in correcting chronic inflammation. It has been considered (61) that regulating GM abundance by supplementation with probiotics and prebiotics is a proven prevention and treatment method to modify the inflammatory status of CKD patients effectively. Our findings indirectly confirmed that differences in GM abundance do have an effect on inflammatory status.
In the present study, high abundance of family Prevotellaceae, class Mollicutes, and phylum Tenericutes were negatively correlated with CRP. Vaziri et al. (62) found that in stool samples from ESRD patients, family Prevotellaceae abundance decreased, leading to increased levels of inflammatory factors (such as CRP) in patients. This was consistent with some of our findings. Family Prevotellaceae is the main bacterial taxon producing SCFAs. SCFAs play an important role in maintaining intestinal epithelial integrity and activating FFAR2 to suppress inflammation (63). Li et al. (64) found that patients with CKD4-5 had improved inflammatory status after appropriate intake of SCFAs. Therefore, further study of GM such as family Prevotellaceae, class Mollicutes and phylum Tenericutes may provide theoretical support for the therapeutic strategy to improve chronic systemic inflammation in CKD by intervening GM. In addition to the mentioned GM above are associated with chronic systemic inflammation in CKD, our results also showed that high abundance of class Negativicutes, genus Eisenbergiella, and order Selenomonadales were positively associated with CRP levels.
The mechanisms of GM effects on chronic systemic inflammation in CKD have been proposed so far. The main mechanisms focus on the impact of intestinal endotoxins and GM metabolites (e.g., IS, PCS) (65) on inflammatory responses. This includes the induction of renal tubular oxidative stress and inflammatory responses by IS and PCS (65), and the induction of cytokine secretion by renal interstitial cells to promote inflammatory responses (66). However, no human or rodent model studies have yet explored class Negativicutes, genus Eisenbergiella, and order Selenomonadales, so researches for these bacteria taxa may become an approach to discovering novel markers for predicting chronic inflammation in CKD.
This study has certain advantages and limitations. First, the advantage is that MR analysis was used, which is a novel genetic statistical method to detect causality. It can avoid confounders and reverse causality. Limitations include (I) GM data were only classified above the genus level. Thus, the causality of CKD and chronic systemic inflammation in CKD could not be concluded at species, strain, or more specialized levels. (II) Since the GWAS data only included European ancestry, the conclusions of this research could not be extended to cover other races. (III) Considering the multi-stage statistics process and biological plausibility, stringent multiple testing calibration may overlook the potential GM with CKD and chronic systemic inflammation in CKD. Hence, we did not strictly follow the multi-corrected p-values to screen GM. (IV) Classes are subcategories of phyla. So the SNP data contained in phylum, class, order, family, and genus may have significant overlap, which may lead to the reproducibility of MR analysis results.
5 Conclusion
Four sets of causal relationships were confirmed in this work: GM and CKD, GM and eGFR, GM and UACR, GM and CRP; however, the reverse causality was not identified. In addition, the present study identified specific bacterial taxa associated with CKD, eGFR, UACR, and CRP, respectively, which might be novel biomarkers for further research of CKD and chronic systemic inflammation in CKD. The study of these bacterial taxa may contribute to the prevention and treatment of CKD and chronic systemic inflammation in CKD, and provide a theoretical basis for studying the mechanism governing gut-kidney axis.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Files, further inquiries can be directed to the corresponding author.
Author contributions
FR: Formal Analysis, Writing – original draft. QiuJ: Formal Analysis, Writing – review and editing. QiJ: Formal Analysis, Writing – review and editing. YQ: Data curation, Formal Analysis, Writing – review and editing. XR: Data curation, Writing – review and editing. TL: Writing – review and editing. YZ: Conceptualization, Funding acquisition, Investigation, Writing – review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was jointly supported by the National Natural Science Foundation of China (No. 82074393) and the Science and Technology Innovation Project of China Academy of Chinese Medical Sciences (No. CI2021A01210).
Acknowledgments
We are grateful to researchers who provided assistance with our study. We thank the MiBioGen consortium for publishing the GWAS summative data on GM. We appreciate CKDGen consortium for providing GWAS summary data of CKD, eGFR and UACR. We are thankful to IEU for providing GWAS data of CRP. We are grateful for expertise assistance in Figures 1, 4 from Figdraw (www.figdraw.com).
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/fimmu.2023.1287698/full#supplementary-material
Supplementary Figure 1 | (A) The result of sensitivity analysis of the causal effect of family Lachnospiraceae on CKD risk. (B) The result of sensitivity analysis of the causal effect of genus Eubacteriumeligens group on CKD risk. (C) The result of sensitivity analysis of the causal effect of genus Intestinimonas on CKD risk. (D) The result of sensitivity analysis of the causal effect of genus LachnospiraceaeUCG010 on CKD risk. (E) The result of sensitivity analysis of the causal effect of genus Streptococcus on CKD risk. (F) The result of sensitivity analysis of the causal effect of genus Desulfovibrionales on CKD risk. (G) The result of sensitivity analysis of the causal effect of order Desulfovibrio on eGFR. (H) The result of sensitivity analysis of the causal effect of class Mollicutes on CRP. (I) The result of sensitivity analysis of the causal effect of class Negativicutes on CRP. (J) The result of sensitivity analysis of the causal effect of family Prevotellaceae on CRP. (K) The result of sensitivity analysis of the causal effect of genus Eisenbergiella on CRP. (L) The result of sensitivity analysis of the causal effect of order Selenomonadales on CRP. (M) The result of sensitivity analysis of the causal effect of phylum Tenericutes onCRP. (N) The result of sensitivity analysis of the causal effect of genus Parasutterella on UACR.
References
1. Lameire NH, Levin A, Kellum JA, Cheung M, Jadoul M, Winkelmayer WC, et al. Harmonizing acute and chronic kidney disease definition and classification: report of a Kidney Disease: Improving Global Outcomes (KDIGO) Consensus Conference. Kidney Int (2021) 100:516–26. doi: 10.1016/j.kint.2021.06.028
2. Köttgen A, Pattaro C, Böger CA, Fuchsberger C, Olden M, Glazer NL, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet (2010) 42:376–84. doi: 10.1038/ng.568
3. Böger CA, Chen MH, Tin A, Olden M, Köttgen A, de Boer IH, et al. CUBN is a gene locus for albuminuria. J Am Soc Nephrol (2011) 22:555–70. doi: 10.1681/asn.2010060598
4. Lv JC, Zhang LX. Prevalence and disease burden of chronic kidney disease. Adv Exp Med Biol (2019) 1165:3–15. doi: 10.1007/978-981-13-8871-2_1
5. Glassock RJ, Warnock DG, Delanaye P. The global burden of chronic kidney disease: estimates, variability and pitfalls. Nat Rev Nephrol (2017) 13:104–14. doi: 10.1038/nrneph.2016.163
6. Akchurin OM, Kaskel F. Update on inflammation in chronic kidney disease. Blood Purif (2015) 39:84–92. doi: 10.1159/000368940
7. Bazeley J, Bieber B, Li Y, Morgenstern H, de Sequera P, Combe C, et al. C-reactive protein and prediction of 1-year mortality in prevalent hemodialysis patients. Clin J Am Soc Nephrol (2011) 6:2452–61. doi: 10.2215/cjn.00710111
8. Navaneethan SD, Zoungas S, Caramori ML, Chan JCN, Heerspink HJL, Hurst C, et al. Diabetes management in chronic kidney disease: synopsis of the 2020 KDIGO clinical practice guideline. Ann Intern Med (2021) 174:385–94. doi: 10.7326/m20-5938
9. Kim HW, Park JT, Joo YS, Kang SC, Lee JY, Lee S, et al. Systolic blood pressure and chronic kidney disease progression in patients with primary glomerular disease. J Nephrol (2021) 34:1057–67. doi: 10.1007/s40620-020-00930-x
10. Hsu CN, Tain YL. Developmental origins of kidney disease: why oxidative stress matters? Antioxidants (Basel) (2020) 10:33. doi: 10.3390/antiox10010033
11. Carrero JJ, Yilmaz MI, Lindholm B, Stenvinkel P. Cytokine dysregulation in chronic kidney disease: how can we treat it? Blood Purif (2008) 26:291–9. doi: 10.1159/000126926
12. Sabatino A, Regolisti G, Brusasco I, Cabassi A, Morabito S, Fiaccadori E. Alterations of intestinal barrier and microbiota in chronic kidney disease. Nephrol Dial Transplant (2015) 30:924–33. doi: 10.1093/ndt/gfu287
13. Yang T, Richards EM, Pepine CJ, Raizada MK. The gut microbiota and the brain-gut-kidney axis in hypertension and chronic kidney disease. Nat Rev Nephrol (2018) 14:442–56. doi: 10.1038/s41581-018-0018-2
14. Li F, Wang M, Wang J, Li R, Zhang Y. Alterations to the gut microbiota and their correlation with inflammatory factors in chronic kidney disease. Front Cell Infect Microbiol (2019) 9:206. doi: 10.3389/fcimb.2019.00206
15. Meijers BK, Evenepoel P. The gut-kidney axis: indoxyl sulfate, p-cresyl sulfate and CKD progression. Nephrol Dial Transplant (2011) 26:759–61. doi: 10.1093/ndt/gfq818
16. Lau WL, Kalantar-Zadeh K, Vaziri ND. The gut as a source of inflammation in chronic kidney disease. Nephron (2015) 130:92–8. doi: 10.1159/000381990
17. Kikuchi K, Saigusa D, Kanemitsu Y, Matsumoto Y, Thanai P, Suzuki N, et al. Gut microbiome-derived phenyl sulfate contributes to albuminuria in diabetic kidney disease. Nat Commun (2019) 10:1835. doi: 10.1038/s41467-019-09735-4
18. Tao S, Tao S, Cheng Y, Liu J, Ma L, Fu P. Effects of probiotic supplements on the progression of chronic kidney disease: A meta-analysis. Nephrol (Carlton) (2019) 24:1122–30. doi: 10.1111/nep.13549
19. Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32:1–22. doi: 10.1093/ije/dyg070
20. Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. Jama (2017) 318:1925–26. doi: 10.1001/jama.2017.17219
21. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet (2014) 23:R89–98. doi: 10.1093/hmg/ddu328
22. Kurilshikov A, Medina-Gomez C, Bacigalupe R, Radjabzadeh D, Wang J, Demirkan A, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet (2021) 53:156–65. doi: 10.1038/s41588-020-00763-1
23. Levey AS, Inker LA. GFR as the "Gold standard": estimated, measured, and true. Am J Kidney Dis (2016) 67:9–12. doi: 10.1053/j.ajkd.2015.09.014
24. Marini S, Georgakis MK, Chung J, Henry JQA, Dichgans M, Rosand J, et al. Genetic overlap and causal inferences between kidney function and cerebrovascular disease. Neurology (2020) 94:e2581–e91. doi: 10.1212/wnl.0000000000009642
25. Filippatos G, Anker SD, Agarwal R, Ruilope LM, Rossing P, Bakris GL, et al. Finerenone reduces risk of incident heart failure in patients with chronic kidney disease and type 2 diabetes: analyses from the FIGARO-DKD trial. Circulation (2022) 145:437–47. doi: 10.1161/circulationaha.121.057983
26. Carlsson AC, Juhlin CC, Larsson TE, Larsson A, Ingelsson E, Sundström J, et al. Soluble tumor necrosis factor receptor 1 (sTNFR1) is associated with increased total mortality due to cancer and cardiovascular causes - findings from two community based cohorts of elderly. Atherosclerosis (2014) 237:236–42. doi: 10.1016/j.atherosclerosis.2014.09.005
27. Wuttke M, Li Y, Li M, Sieber KB, Feitosa MF, Gorski M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet (2019) 51:957–72. doi: 10.1038/s41588-019-0407-x
28. Mitchell R, Elsworth BL, Raistrick CA, Paternoster L, Hemani G, Gaunt TR, et al. MRC IEU UK Biobank GWAS pipeline version 2. University of Bristol. (2019). doi: 10.5523/bris.pnoat8cxo0u52p6ynfaekeigi.
29. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med (2008) 27:1133–63. doi: 10.1002/sim.3034
30. Mensah-Kane J, Schmidt AF, Hingorani AD, Finan C, Chen Y, van Duijvenboden S, et al. No clinically relevant effect of heart rate increase and heart rate recovery during exercise on cardiovascular disease: A mendelian randomization analysis. Front Genet (2021) 12:569323. doi: 10.3389/fgene.2021.569323
31. Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol (2011) 40:755–64. doi: 10.1093/ije/dyr036
32. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol (2013) 37:658–65. doi: 10.1002/gepi.21758
33. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med (2015) 34:2926–40. doi: 10.1002/sim.6522
34. Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol (2017) 46:1985–98. doi: 10.1093/ije/dyx102
35. 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:512–25. doi: 10.1093/ije/dyv080
36. Ni JJ, Xu Q, Yan SS, Han BX, Zhang H, Wei XT, et al. Gut microbiota and psychiatric disorders: A two-sample mendelian randomization study. Front Microbiol (2021) 12:737197. doi: 10.3389/fmicb.2021.737197
37. 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:693–98. doi: 10.1038/s41588-018-0099-7
38. Xiang K, Wang P, Xu Z, Hu YQ, He YS, Chen Y, et al. Causal effects of gut microbiome on systemic lupus erythematosus: A two-sample mendelian randomization study. Front Immunol (2021) 12:667097. doi: 10.3389/fimmu.2021.667097
39. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol (2011) 40:740–52. doi: 10.1093/ije/dyq151
40. Mazidi M, Shekoohi N, Covic A, Mikhailidis DP, Banach M. Adverse Impact of Desulfovibrio spp. and Beneficial Role of Anaerostipes spp. on Renal Function: Insights from a Mendelian Randomization Analysis. Nutrients (2020) 12:2216. doi: 10.3390/nu12082216
41. Wang X, Yang S, Li S, Zhao L, Hao Y, Qin J, et al. Aberrant gut microbiota alters host metabolome and impacts renal failure in humans and rodents. Gut (2020) 69:2131–42. doi: 10.1136/gutjnl-2019-319766
42. Mahmoodpoor F, Rahbar Saadat Y, Barzegari A, Ardalan M, Zununi Vahed S. The impact of gut microbiota on kidney function and pathogenesis. BioMed Pharmacother (2017) 93:412–19. doi: 10.1016/j.biopha.2017.06.066
43. Zhu Y, He H, Tang Y, Peng Y, Hu P, Sun W, et al. Reno-protective effect of low protein diet supplemented with α-ketoacid through gut microbiota and fecal metabolism in 5/6 nephrectomized mice. Front Nutr (2022) 9:889131. doi: 10.3389/fnut.2022.889131
44. Guirong YE, Minjie Z, Lixin YU, Junsheng YE, Lin Y, Lisha S. [Gut microbiota in renal transplant recipients, patients with chronic kidney disease and healthy subjects]. Nan Fang Yi Ke Da Xue Xue Bao (2018) 38:1401–08. doi: 10.12122/j.issn.1673-4254.2018.12.01
45. Lai L, Li Y, Liu J, Luo L, Tang J, Xue J, et al. Bovine serum albumin aggravates macrophage M1 activation and kidney injury in heterozygous Klotho-deficient mice via the gut microbiota-immune axis. Int J Biol Sci (2021) 17:742–55. doi: 10.7150/ijbs.56424
46. Szeto CC, McIntyre CW, Li PK. Circulating bacterial fragments as cardiovascular risk factors in CKD. J Am Soc Nephrol (2018) 29:1601–08. doi: 10.1681/asn.2018010068
47. Andrade-Oliveira V, Amano MT, Correa-Costa M, Castoldi A, Felizardo RJ, de Almeida DC, et al. Gut bacteria products prevent AKI induced by ischemia-reperfusion. J Am Soc Nephrol (2015) 26:1877–88. doi: 10.1681/asn.2014030288
48. Hu X, Ouyang S, Xie Y, Gong Z, Du J. Characterizing the gut microbiota in patients with chronic kidney disease. Postgrad Med (2020) 132:495–505. doi: 10.1080/00325481.2020.1744335
49. Sallée M, Dou L, Cerini C, Poitevin S, Brunet P, Burtey S. The aryl hydrocarbon receptor-activating effect of uremic toxins from tryptophan metabolism: a new concept to understand cardiovascular complications of chronic kidney disease. Toxins (Basel) (2014) 6:934–49. doi: 10.3390/toxins6030934
50. Jiang S, Xie S, Lv D, Wang P, He H, Zhang T, et al. Alteration of the gut microbiota in Chinese population with chronic kidney disease. Sci Rep (2017) 7:2870. doi: 10.1038/s41598-017-02989-2
51. Zhao J, Ning X, Liu B, Dong R, Bai M, Sun S. Specific alterations in gut microbiota in patients with chronic kidney disease: an updated systematic review. Ren Fail (2021) 43:102–12. doi: 10.1080/0886022x.2020.1864404
52. Gibson GR. Physiology and ecology of the sulphate-reducing bacteria. J Appl Bacteriol (1990) 69:769–97. doi: 10.1111/j.1365-2672.1990.tb01575.x
53. Quigley EMM. Nutraceuticals as modulators of gut microbiota: Role in therapy. Br J Pharmacol (2020) 177:1351–62. doi: 10.1111/bph.14902
54. Feng YL, Cao G, Chen DQ, Vaziri ND, Chen L, Zhang J, et al. Microbiome-metabolomics reveals gut microbiota associated with glycine-conjugated metabolites and polyamine metabolism in chronic kidney disease. Cell Mol Life Sci (2019) 76:4961–78. doi: 10.1007/s00018-019-03155-9
55. Mishima E, Fukuda S, Mukawa C, Yuri A, Kanemitsu Y, Matsumoto Y, et al. Evaluation of the impact of gut microbiota on uremic solute accumulation by a CE-TOFMS-based metabolomics approach. Kidney Int (2017) 92:634–45. doi: 10.1016/j.kint.2017.02.011
56. Brinkkoetter PT, Ising C, Benzing T. The role of the podocyte in albumin filtration. Nat Rev Nephrol (2013) 9:328–36. doi: 10.1038/nrneph.2013.78
57. Hoshino J, Furuichi K, Yamanouchi M, Mise K, Sekine A, Kawada M, et al. A new pathological scoring system by the Japanese classification to predict renal outcome in diabetic nephropathy. PloS One (2018) 13:e0190923. doi: 10.1371/journal.pone.0190923
58. Watkins EB, Phillips RS. Inhibition of tyrosine phenol-lyase from Citrobacter freundii by 2-azatyrosine and 3-azatyrosine. Biochemistry (2001) 40:14862–8. doi: 10.1021/bi015707s
59. Carrero JJ, Stenvinkel P. Inflammation in end-stage renal disease–what have we learned in 10 years? Semin Dial (2010) 23:498–509. doi: 10.1111/j.1525-139X.2010.00784.x
60. Lau WL, Savoj J, Nakata MB, Vaziri ND. Altered microbiome in chronic kidney disease: systemic effects of gut-derived uremic toxins. Clin Sci (Lond) (2018) 132:509–22. doi: 10.1042/cs20171107
61. Tsai YL, Lin TL, Chang CJ, Wu TR, Lai WF, Lu CC, et al. Probiotics, prebiotics and amelioration of diseases. J BioMed Sci (2019) 26:3. doi: 10.1186/s12929-018-0493-6
62. Vaziri ND, Wong J, Pahl M, Piceno YM, Yuan J, DeSantis TZ, et al. Chronic kidney disease alters intestinal microbial flora. Kidney Int (2013) 83:308–15. doi: 10.1038/ki.2012.345
63. Hu Y, Chen D, Zheng P, Yu J, He J, Mao X, et al. The bidirectional interactions between resveratrol and gut microbiota: an insight into oxidative stress and inflammatory bowel disease therapy. BioMed Res Int (2019) 2019:5403761. doi: 10.1155/2019/5403761
64. Li YJ, Chen X, Kwan TK, Loh YW, Singer J, Liu Y, et al. Dietary Fiber Protects against Diabetic Nephropathy through Short-Chain Fatty Acid-Mediated Activation of G Protein-Coupled Receptors GPR43 and GPR109A. J Am Soc Nephrol (2020) 31:1267–81. doi: 10.1681/asn.2019101029
65. Frosali S, Pagliari D, Gambassi G, Landolfi R, Pandolfi F, Cianci R. How the intricate interaction among toll-like receptors, microbiota, and intestinal immunity can influence gastrointestinal pathology. J Immunol Res (2015) 2015:489821. doi: 10.1155/2015/489821
Keywords: two-sample Mendelian randomization, gut microbiota, chronic kidney disease, inflammation, causality two-sample mendelian randomization, causality
Citation: Ren F, Jin Q, Jin Q, Qian Y, Ren X, Liu T and Zhan Y (2023) Genetic evidence supporting the causal role of gut microbiota in chronic kidney disease and chronic systemic inflammation in CKD: a bilateral two-sample Mendelian randomization study. Front. Immunol. 14:1287698. doi: 10.3389/fimmu.2023.1287698
Received: 02 September 2023; Accepted: 20 October 2023;
Published: 02 November 2023.
Edited by:
Marijn M. Speeckaert, Ghent University Hospital, BelgiumReviewed by:
Dengpiao Xie, The Teaching Hospital of Chengdu University of Traditional Chinese Medicine, ChinaWondwossen Amogne Degu, Addis Ababa University, Ethiopia
Copyright © 2023 Ren, Jin, Jin, Qian, Ren, Liu and Zhan. 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: Yongli Zhan, zhanyongli2023@163.com
†These authors have contributed equally to this work