Skip to main content

Mendelian randomization study of urolithiasis: exploration of risk factors using human blood metabolites

Abstract

Background

Urolithiasis is a highly prevalent global disease closely associated with metabolic factors; however, the causal relationship between blood metabolites and urolithiasis remains poorly understood.

Method

In our study, we employed a bi-directional two-sample Mendelian randomization (MR) analysis to investigate the causal associations between urolithiasis and metabolites. The random-effects inverse-variance weighted (IVW) estimation method was utilized as the primary approach, complemented by several other estimators including MR-Egger, weighted median, colocalization and MR-PRESSO. Furthermore, the study included replication and meta-analysis. Finally, we conducted metabolic pathway analysis to elucidate potential metabolic pathways.

Results

After conducting multiple tests for correction, glycerol might contribute to the urolithiasis and dehydroisoandrosterone sulfate (DHEA-S) might inhibit this process. Furthermore, several blood metabolites had shown potential associations with a causal relationship. Among the protective metabolites were lipids (dehydroisoandrosterone sulfate and 1-stearoylglycerol (1-monostearin)), amino acids (isobutyrylcarnitine and 2-aminobutyrate), a keto acid (acetoacetate) and a carbohydrate (mannose). The risk metabolites included lipids (1-palmitoylglycerophosphoethanolamine, glycerol and cortisone), a carbohydrate (erythronate), a peptide (pro-hydroxy-pro) and a fatty acid (eicosenoate). In reverse MR analysis, urolithiasis demonstrated a statistically significant causal relationship with butyrylcarnitine, 3-methyl-2-oxobutyrate, scyllo-inositol, leucylleucine and leucylalanine. However, it was worth noting that none of the blood metabolites exhibited statistical significance after multiple corrections. Additionally, we identified one metabolic pathway associated with urolithiasis.

Conclusion

The results we obtained demonstrate the causal relevance between two metabolites and urolithiasis, as well as identify one metabolic pathway potentially associated with its development. Given the high prevalence of urolithiasis, further investigations are encouraged to elucidate the mechanisms of these metabolites and explore novel therapeutic strategies.

Peer Review reports

Introduction

Urolithiasis is common in urologic diseases, with a continually increasing prevalence and incidence [1, 2]. Approximately 10% of the global population experiences kidney stone occurrence at least once during their lifetime, with a recurrence rate of 2% among affected individuals [3]. The metabolic risk factors associated with urinary stones are receiving increasing attention as our understanding of the etiology of urinary stones continues to deepen [4]. In order to evaluate the metabolic disorder associated with stone formation, some studies have applied metabolomics to urolithiasis. Through preliminary metabolic analysis of urine in renal stone patients, Duan et al. [5]. have found that four metabolic pathways, namely acetic acid and dicarboxylic acid metabolism, glycine, serine, and threonine metabolism, phenylalanine metabolism, and citric acid cycle, are closely associated with urolithiasis. Lately, Zhang et al. [6]. applied metabolomics technologies to discover the role of succinate in combating stone formation. In Agudelo et al.’s study [7], the significant enrichment of metabolites in the stone-forming group suggested that lithogenic metabolites in the urinary tract might be a crucial driver of stone formation. These studies are indeed promising to contribute to the targeted exploration of certain metabolites or metabolic pathways to identify biomarkers for urolithiasis. However, their analysis was limited to investigating the association between urolithiasis and urine metabolomics. In comparison to urinary metabolomics, blood metabolites offer the advantages of being easily obtainable in large quantities with good stability. Moreover, they provide a wealth of information [8], making them a promising option for early disease detection [9]. However, the detailed pathophysiological mechanisms of blood metabolites in urolithiasis have not yet been elucidated. Therefore, to clarify the causal relationship of blood metabolites in the pathogenesis of urolithiasis, a comprehensive and complete analysis is urgently needed. The best methods to investigate causality are randomized controlled trials (RCTs) due to their ability to mitigate reverse causality and residual confounding through randomization. However, the lengthy duration and high cost of RCTs pose significant challenges. Under this background, Mendelian randomization (MR) is a new way which can examine the causality between exposure and outcome, using single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) [10]. Additionally, the results of traditional observational studies are biased in the estimation of causal effects due to reverse causality and residual confounders whereas MR studies are generally unlikely impacted by confounders because the genotypes assignment from parents to offspring is random. MR studies are generally unlikely impacted by confounders because the genotypes assignment from parents to offspring is random [11]. Here, we use MR methods to analyze to determine the potential causal associations of metabolites with the risk of urolithiasis.

Methods

Data sets

Human blood metabolome genome-wide association study (GWAS) dataset was obtained from IEU OpenGWAS database (https://gwas.mrcieu.ac.uk/) [12, 13]. It’s worth noting that this study identified 2.1 million SNPs for 486 metabolites from 7,824 European which was conducted by Shin et al [14]. In addition, we also analyzed GWAS dataset from 24,925 European which was conducted by Kettumen et al [15].

Urolithiasis GWAS summary statistics data released FinnGen consortium (URL: https://r8.finngen.fi/pheno/N14_CALCUKIDUR), which included 4,969 cases and 213,445 controls [16]. The summary of data sources was presented in Supplementary Table S1.

Mendelian randomization analysis

The MR flowchart is shown in Fig. 1. In order to reduce the deviation caused by genetic variables, IVs should rely on three essential assumptions, which are elucidated in Fig. 1: (1) the SNPs should be closely associated with metabolites; (2) the SNPs are independent with any confounders; (3) the SNPs should affect the risk of urolithiasis only via metabolites and not through any other pathway [17, 18]. The Inverse Variance Weighted (IVW) method was widely acknowledged as a more quick, convenient and common approach for analysis, thus we adopted IVW as the primary analytical method. In cases where heterogeneity existed among causal estimates of different variants (as demonstrated in this article), the random effects model became more appropriate [19]. Furthermore, two additional methods were used to complement the analysis presented in this study. Specifically, MR egger [20] was employed for the identification and adjustment of pleiotropy effects, while weighted median [21] analysis was utilized to mitigate potential biases in strong hypotheses that hold true for all instrumental variables (IVs) in IVW. Firstly, we relaxed standards and chose IVs with the significance threshold (p < 1 × 10− 5) given that the scarcity of SNPs reaching genome-wide significance [22], and we used the clumping method (r2 < 0.001 and clump distance < 10,000 kb) to exclude SNPs based on European lineage reference data from the 1000 Genome Project, refraining from biased results originated from strong linkage disequilibrium (LD) [23, 24]. Another important point is that during the harmonizing process, palindromic SNPs were excluded to ensure the effects of SNPs on exposure accorded with the same allele as the effects of SNPs on the outcome.

Fig. 1
figure 1

Overview of the present MR study design. SNP, single-nucleotide polymorphism; IV, instrumental variable; MR, Mendelian randomization

Compared with two sample MR analysis, bidirectional MR analysis can solve the potential problem of causal entanglement. By conducting two-sample MR analysis from both directions to ascertain the direction of causal relationships, we were able to mitigate confusion arising from reverse causality and achieve a more comprehensive understanding of causal pathways. Therefore, we also performed reverse MR analysis on urolithiasis to assess its potential impact on blood metabolites. During this stage, we applied a P-value threshold of P < 5 × 10− 8, which aligns with the approach employed in forward MR analysis.

Sensitivity analysis

Sensitivity analysis includes Cochran’s Q test and the MR-Egger test in order to assess the significance of our results. Cochran’s Q statistic was applied to estimate the heterogeneity among SNPs associated with each metabolite [25]. We used MR-Egger regression and MR-Presso tests to evaluate whether genetic instruments had made pleiotropic effects on the outcome [20]. In addition, we excluded IVs with F statistics < 10, the F-statistic was defined as the ratio of the model’s mean square to that of the error: \(\:\text{F}=\frac{{R}^{2}(n-1-k)}{(1-{R}^{2})k}\:,\) in accordance with its academic and professional significance [26]. We implemented all MR analyses in R (version 4.2.1) using R package TwoSampleMR [13] to detect the causal effects of different blood metabolites on the risk of urolithiasis. The statistical significance was considered when the P-value was less than 0.05. In addition, the OR value was calculated based on the results obtained from IVW, and if it exceeded 1 and the P-value was less than 0.05, it indicated a significant risk factor for urolithiasis; conversely, if it was below 1 and the P-value was less than 0.05, it suggested a protective factor against urolithiasis.

Multiple-testing correction

The FDR (False Discovery Rate) was utilized for the correction of all P-values. The significance threshold, denoted as q, was adjusted to be less than 0.05. When a correlation between urolithiasis and blood metabolites was observed with P-values below 0.05 and q-values greater than or equal to 0.05, it suggests a potential association.

Power calculation

We employed a specialized online tool (https://shiny.cnsgenomics.com/mRnd/) [27], which utilizes asymptotic theory to estimate power values for detecting causal effects derived from IVs, to assess the statistical power of MR. We conducted power calculations at a type I error rate of 0.05, considering factors such as OR obtained from MR analyses utilizing the IVW approach, R2 of IVs and the proportion of cases of urolithiasis GWAS.

Replication and meta-analysis

In order to strengthen the robustness of our results, the replication and meta-analysis of MR were expanded by integrating additional GWAS datasets. The GWAS datasets were accessible through the IEU OpenGWAS database (https://gwas.mrcieu.ac.uk/). The GWAS datasets information could be found in Supplementary Table S1. The IVW method was the main methods for the replication and meta-analysis was utilized for combining the outcomes from these GWAS datasets.

Colocalization analysis

In order to determine whether the associations of the identified blood metabolites with urolithiasis were influenced by a shared causal variant, we utilized the R package coloc (v5.2.3) to complete this step, which employed Bayesian colocalization analysis. This analysis assessed five corresponding posterior probabilities of its following hypotheses, including H0 (no correlation with either trait); H1 (solely associated with Trait 1); H2 (solely associated with Trait 2); H3 (two traits are associated but with different causal variations) and H4 (two traits are associated and share a causal variation) [28]. H4/(H3 + H4) reveals the probability of colocalization given the presence of a causal variant for urolithiasis [29].

Metabolic pathway analysis

Using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/) to investigate the association of metabolic pathways with urolithiasis, the Kyoto Encyclopedia of Genes and Genomes (KEGG) [30] database was analyzed.

Results

A total of 296 SNPs (P < 1 × 10− 5) associated with 12 traits were identified for human blood metabolites. None of the F- statistics were less than 10, indicating a significant correlation between SNPs and metabolites (Supplementary Table S2).

MR analysis results of human blood metabolite

Firstly, we identified 12 human blood metabolites that were significantly associated with urolithiasis (P < 0.05), including lipids, amino acids, fatty acids, carbohydrates, peptide and keto acid (Fig. 2; Table 1). Several significant findings regarding clinically relevant blood metabolites were found in our study. Specifically, cortisone (OR: 2.18 (95%CI: 1.05–4.52), P = 0.035), glycerol (OR:1.38 (95%CI: 1.18–1.62), P < 0.001) were identified as risk factors for urolithiasis, whereas 2-aminobutyrate (OR: 0.47 (95%CI: 0.23–0.95), P = 0.035), dehydroisoandrosterone sulfate (DHEA-S) (OR: 0.57 (95%CI: 0.40–0.82), P = 0.002), mannose (OR: 0.35 (95%CI: 0.16–0.76), P = 0.008) and acetoacetate (OR:0.83 (95%CI:0.70–0.99), P = 0.040) were identified as protective factors for urolithiasis (Fig. 2; Table 1). After applying FDR correction, only dehydroisoandrosterone sulfonate (DHEA-S) (q = 0.045) and glycerol (q = 0.001) exhibited statistically significant differences. Pro-hydroxy-pro, Erythronate and Cortisone had a high statistical power of 1.00 (Supplementary Table S3). Furthermore, we conducted metabolite pathway analysis on all metabolites discovered by the IVW method (P < 0.05). As shown in Table 2, one metabolic pathway was significantly causal with urolithiasis.

Table 1 MR analysis for the association between blood metabolites and urolithiasis
Table 2 Sensitivity analysis of the causal association between blood metabolites and urolithiasis
Fig. 2
figure 2

Forest plot for evaluating the causal relationship between urolithiasis and blood metabolites based on the values obtained from the IVW method

Sensitivity analysis results

Except for mannose (P = 0.044) and erythronate (P = 0.011), no significant heterogeneity of IVs was showed according to the results of Cochran’s MR Egger Q test and Cochran’s IVW Q test (Table 3). Through the MR-PRESSO test, the IV rs6860069 (Rssobs = 0.0131, P < 0.015) was identified as outlier and removed from the next analysis. In addition, according to the results of the MR-Egger intercept tests and the MR-Presso tests (Table 3), it suggested that erythronate has horizontal pleiotropy. However, the subsequent distortion test P-value was greater than 0.05, which indicated that it did not affect our results. Therefore, our process was consistent with MR assumptions.

Table 3 Significant metabolic pathways in the pathogenesis of urolithiasis

Reverse MR analyses results

Furthermore, employing reverse MR analysis revealed a modest association among urolithiasis and butyrylcarnitine (β = 0.04, P = 0.029), 3-methyl-2-oxobutyrate (β = -0.02, P = 0.043), scyllo-inositol (β = 0.04, P = 0.024), leucylleucine (β = 0.05, P = 0.007), x-14,304—leucylalanine (β = 0.05, P = 0.028); however, after implementing multiple corrections, none of the individual findings reached the threshold for statistical significance (Table 4).

Table 4 Reverse MR analysis for the association between blood metabolites and urolithiasis

Replication and meta-analysis

In order to strengthen the robustness of our results, the replication and meta-analysis of MR were expanded by integrating additional GWAS datasets for the positively identified metabolites after multiple-testing correction. The findings revealed that glycerol demonstrated similar trends of causal associations with urolithiasis in other GWAS datasets (Supplementary Figures S1-S2). However, the results of the replication and meta-analysis were not significant, possibly due to the broader heterogeneity across GWAS datasets (Heterogeneity: P < 0.01) and including the heterogeneity statistics.

Colocalization analyses

Two metabolites with FDR significant MR associations with urolithiasis were performed colocalization analysis, and the results were presented in the Supplementary Table S4. The results of the colocalization analysis suggested that the connections between urolithiasis and the two established metabolites were not linked to shared causal variant sites. The regional associations identified in the colocalization results were illustrated on Supplementary Figures S3–S4.

Conclusion

By applying mendelian randomization and after applying FDR correction, only dehydroisoandrosterone sulfonate (DHEA-S) and glycerol exhibited statistically significant differences. 2-aminobutyrate, Isobutyrylcarnitine, mannose, acetoacetate and 1-stearoylglycerol (1-monostearin) were detected to possess suggestive protective effects against urolithiasis. On the contrary, a number of metabolites, incorporating eicosenoate, pro-hydroxy-pro, erythronate, 1-palmitoylglycerophosphoethanolamine and cortisone had suggestive negative effects on urolithiasis. Given that the primary purpose of this study was to explore and discover as many potential significant metabolites as possible, we posited that metabolites with p-values less than 0.05 and q-values greater than 0.05 also warranted consideration. Nevertheless, it should be noted that these potential findings have not been subjected to FDR correction, thus necessitating cautious validation in larger sample sizes in future studies.

In our study, we observed that DHEA-S exerted inhibitory effects on stone formation, highlighting its potential as a promising preventive agent in the field of urolithiasis. In an animal experiment, raising castrated rats with DHEA may directly modulate the hepatic enzyme activities of GRHPR and AGXT which subsequently regulate the endogenous oxalate production in the liver [31]. On the basis of it, Fuster et al [32]. conducted a cross-sectional analysis aiming to reveal the relationship between urinary sex hormones and excretion of urinary components in kidney stone formers. Of note, their result shown that DHEA had an inverse association with urinary oxalate excretion and supported our finding. Additionally, In Franca Serafini-Cessi et al.’s study, they described N-Glycans, which are rich in mannose, were capable to resist urological diseases [33]. S Proietti et al. assessed a D-mannose-containing product possessing protective effects against infection-related urinary stones [34]. Their experimental conclusions fully support the analysis results obtained in our MR analysis regarding the protective effect of mannose.

In parallel with these protective metabolites, it was of great interest to focus on the risk factors derived from our result. The most representative of these was cortisone, a glucocorticoid, which captured our attention. Several observational studies on hormone concentration in urine samples had indicated that glucocorticoids can mediate a negative impact on the excretion of inorganic salts and uric acid, even at normal physiological levels [35, 36]. Likewise, as an inactive precursor of glucocorticoids such as hydrocortisone, cortisone was a significant risk factor for stone formation in our analysis.

Apart from the metabolites mentioned, we found a significant relationship between steroid hormone biosynthesis and the formation of urolithiasis. In Wen et al.’s study [37], steroid biosynthesis was found to be altered in patients with urolithiasis. In addition, it was manifested that steroid derivatives are also potentially relevant to urolithiasis, which coincides with the protective and lithogenic effects of DHEA and cortisone, as concluded in our study, respectively. These metabolites and the altered metabolic pathway may collectively suggest the formation or compensatory onset of urolithiasis and may possess profound value as future biomarkers or therapeutic target sites for urolithiasis.

However, our study also had some limitations. Firstly, to address the issue of limited availability of SNPs for the exposure of interest at a genome-wide level, we had set a more relaxed threshold which was also commonly employed in other studies. Although relaxing the threshold might increase the likelihood of horizontal pleiotropy occurring, our findings confirmed the absence of any additional level of pleiotropy. In addition, the F-statistic value of selected SNPs all exceeded 10, indicating that our IVs were robust enough. Secondly, we did not prove our study in other populations such as Asians which may affect the generalizability of the results. Moreover, it is significant that our study needs to be confirmed through careful basic research.

In general, the results we obtained show us the causal relevance between two metabolites and urolithiasis, and we also ascertained one metabolic pathway that may be related to the development of urolithiasis. Facing the high prevalence of urolithiasis, further investigations are encouraged to clarify the mechanisms of these metabolites and explore new therapeutic strategies.

Data availability

All of the data analyzed in current research are available here: the FinnGen consortium (https://www.finngen.fi/en) and the IEU OpenGWAS database (https://gwas.mrcieu.ac.uk/).

Abbreviations

MR:

Mendelian randomization

IVW:

Inverse-variance weighted

RCTs:

Randomized controlled trials

SNPs:

Single nucleotide polymorphisms

GWAS:

Genome-wide association study

IVs:

Instrumental variables

LD:

Linkage disequilibrium

KEGG:

Kyoto Encyclopedia of Genes and Genomes

SMPDB:

Small Molecular Pathways Database

OR:

Odds ratio

References

  1. López M, Hoppe B. History, epidemiology and regional diversities of urolithiasis. Pediatr Nephrol. 2010;25(1):49–59. https://doi.org/10.1007/s00467-008-0960-5.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Khan A. Prevalence, pathophysiological mechanisms and factors affecting urolithiasis. Int Urol Nephrol. 2018;50(5):799–806. https://doi.org/10.1007/s11255-018-1849-2.

    Article  CAS  PubMed  Google Scholar 

  3. Wagner CA. Etiopathogenic factors of urolithiasis. Arch Esp Urol. 2021;74(1):16–23.

    PubMed  Google Scholar 

  4. Scales CD Jr., Tasian GE, Schwaderer AL, Goldfarb DS, Star RA, Kirkali Z. Urinary Stone Disease: advancing knowledge, patient care, and Population Health. Clin J Am Soc Nephrology: CJASN. 2016;11(7):1305–12. https://doi.org/10.2215/cjn.13251215.

    Article  Google Scholar 

  5. Duan X, Zhang T, Ou L, Kong Z, Wu W, Zeng G. (1)H NMR-based metabolomic study of metabolic profiling for the urine of kidney stone patients. Urolithiasis. 2020;48(1):27–35. https://doi.org/10.1007/s00240-019-01132-2.

    Article  CAS  PubMed  Google Scholar 

  6. Zhang XZ, Lei XX, Jiang YL, Zhao LM, Zou CY, Bai YJ, et al. Application of metabolomics in urolithiasis: the discovery and usage of succinate. Signal Transduct Target Therapy. 2023;8(1):41. https://doi.org/10.1038/s41392-023-01311-z.

    Article  CAS  Google Scholar 

  7. Agudelo J, Fedrigon D, Faris A, Wilkins L, Monga M, Miller AW. Delineating the role of the urinary metabolome in the lithogenesis of calcium-based kidney stones. Urology. 2022;167:49–55. https://doi.org/10.1016/j.urology.2022.06.004.

    Article  PubMed  Google Scholar 

  8. Zhang A, Sun H, Wang X. Serum metabolomics as a novel diagnostic approach for disease: a systematic review. Anal Bioanal Chem. 2012;404(4):1239–45. https://doi.org/10.1007/s00216-012-6117-1.

    Article  CAS  PubMed  Google Scholar 

  9. Amantonico A, Urban PL, Zenobi R. Analytical techniques for single-cell metabolomics: state of the art and trends. Anal Bioanal Chem. 2010;398(6):2493–504. https://doi.org/10.1007/s00216-010-3850-1.

    Article  CAS  PubMed  Google Scholar 

  10. Smith GD, Ebrahim S. Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22. https://doi.org/10.1093/ije/dyg070.

    Article  PubMed  Google Scholar 

  11. TS BS. Mendelian randomization: methods for causal inference using genetic variants. CRC; 2021.

  12. Elsworth B, Lyon M, Alexander T, Liu Y, Matthews P, Hallett J, et al. The MRC IEU OpenGWAS data infrastructure. bioRxiv. 2020. https://doi.org/10.1101/2020.08.10.244293:2020.08.10.244293.

    Article  Google Scholar 

  13. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife. 2018. https://doi.org/10.7554/eLife.34408.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Shin SY, Fauman EB, Petersen AK, Krumsiek J, Santos R, Huang J, et al. An atlas of genetic influences on human blood metabolites. Nat Genet. 2014;46(6):543–50. https://doi.org/10.1038/ng.2982.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kettunen J, Demirkan A, Würtz P, Draisma HHM, Haller T, Rawal R, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7(1):11122. https://doi.org/10.1038/ncomms11122.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613(7944):508–18. https://doi.org/10.1038/s41586-022-05473-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Yuan S, Liu J, Larsson SC. Smoking, alcohol and coffee consumption and pregnancy loss: a mendelian randomization investigation. Fertil Steril. 2021;116(4):1061–7. https://doi.org/10.1016/j.fertnstert.2021.05.103.

    Article  PubMed  Google Scholar 

  18. Burgess S, Thompson SG. Multivariable mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181(4):251–60. https://doi.org/10.1093/aje/kwu283.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiol (Cambridge Mass). 2017;28(1):30–42. https://doi.org/10.1097/ede.0000000000000559.

    Article  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. https://doi.org/10.1093/ije/dyv080.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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(6):1985–98. https://doi.org/10.1093/ije/dyx102.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Zhang T, Cao Y, Zhao J, Yao J, Liu G. Assessing the causal effect of genetically predicted metabolites and metabolic pathways on stroke. J Translational Med. 2023;21(1):822. https://doi.org/10.1186/s12967-023-04677-4.

    Article  CAS  Google Scholar 

  23. Shen HH, Zhang YY, Wang XY, Wang CJ, Wang Y, Ye JF, et al. Potential Causal Association between plasma metabolites, immunophenotypes, and Female Reproductive disorders: a two-sample mendelian randomization analysis. Biomolecules. 2024;14(1). https://doi.org/10.3390/biom14010116.

  24. Roze D. Causes and consequences of linkage disequilibrium among transposable elements within eukaryotic genomes. Genetics. 2023;224(2). https://doi.org/10.1093/genetics/iyad058.

  25. Egger M, Smith GD, Phillips AN. Meta-analysis: principles and procedures. BMJ (Clinical Res ed). 1997;315(7121):1533–7. https://doi.org/10.1136/bmj.315.7121.1533.

    Article  CAS  Google Scholar 

  26. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for mendelian randomization studies using multiple genetic variants. Int J Epidemiol. 2011;40(3):740–52. https://doi.org/10.1093/ije/dyq151.

    Article  PubMed  Google Scholar 

  27. Brion MJ, Shakhbazov K, Visscher PM. Calculating statistical power in mendelian randomization studies. Int J Epidemiol. 2013;42(5):1497–501. https://doi.org/10.1093/ije/dyt179.

    Article  PubMed  Google Scholar 

  28. Foley CN, Staley JR, Breen PG, Sun BB, Kirk PDW, Burgess S, et al. A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits. Nat Commun. 2021;12(1):764. https://doi.org/10.1038/s41467-020-20885-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhao SS, Yiu ZZN, Barton A, Bowes J. Association of lipid-lowering drugs with risk of Psoriasis: a mendelian randomization study. JAMA Dermatology. 2023;159(3):275–80. https://doi.org/10.1001/jamadermatol.2022.6051.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. Nucleic Acids Res. 2012;40(Database issue):D109–14. https://doi.org/10.1093/nar/gkr988. KEGG for integration and interpretation of large-scale molecular data sets.

  31. Nishijima S, Sugaya K, Hokama S, Morozumi M, Ogawa Y. Effect of dehydroepiandrosterone on oxalate metabolism in rats. Front Bioscience: J Virtual Libr. 2004;9:1360–4. https://doi.org/10.2741/1340.

    Article  CAS  Google Scholar 

  32. Fuster DG, Morard GA, Schneider L, Mattmann C, Lüthi D, Vogt B, et al. Association of urinary sex steroid hormones with urinary calcium, oxalate and citrate excretion in kidney stone formers. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association -. Eur Ren Association. 2022;37(2):335–48. https://doi.org/10.1093/ndt/gfaa360.

    Article  CAS  Google Scholar 

  33. Serafini-Cessi F, Monti A, Cavallone D. N-Glycans carried by Tamm-Horsfall glycoprotein have a crucial role in the defense against urinary tract diseases. Glycoconj J. 2005;22:7–9. https://doi.org/10.1007/s10719-005-2142-z.

    Article  CAS  Google Scholar 

  34. Proietti S, Giannantoni A, Luciani LG, Sortino G, Graziotti P, Giusti G. Cystoman® and calculi: a good alternative to standard therapies in preventing stone recurrence. Urolithiasis. 2014;42(4):285–90. https://doi.org/10.1007/s00240-014-0675-y.

    Article  CAS  PubMed  Google Scholar 

  35. Shi L, Berkemeyer S, Buyken AE, Maser-Gluth C, Remer T. Glucocorticoids and body fat associated with renal uric acid and oxalate, but not calcium excretion, in healthy children. Metab Clin Exp. 2010;59(1):134–9. https://doi.org/10.1016/j.metabol.2009.06.027.

    Article  CAS  PubMed  Google Scholar 

  36. Hua Y, Esche J, Hartmann MF, Maser-Gluth C, Wudy SA, Remer T. Cortisol and 11 beta-hydroxysteroid dehydrogenase type 2 as potential determinants of renal citrate excretion in healthy children. Endocrine. 2020;67(2):442–8. https://doi.org/10.1007/s12020-019-02151-0.

    Article  CAS  PubMed  Google Scholar 

  37. Wen J, Cao Y, Li Y, Zhu F, Yuan M, Xu J, et al. Metabolomics analysis of the serum from children with urolithiasis using UPLC-MS. Clin Transl Sci. 2021;14(4):1327–37. https://doi.org/10.1111/cts.12984.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We express our gratitude to the FinnGen consortium and IEU OpenGWAS database for generously sharing GWAS summary statistics data, which greatly facilitated our analysis.

Funding

This work was supported by the National Natural Science Foundation of China (82070724) and National Natural Science Foundation of China (82370768).

Author information

Authors and Affiliations

Authors

Contributions

HZY designed the study. HDK collected and analyzed the data. DAQ and PJS prepared the manuscript and provided the figures. YR and GDF reviewed the manuscript. HBB revised the manuscript and HZY provided constructive suggestions. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Bingbing Hou or Zongyao Hao.

Ethics declarations

Ethics approval and consent to participate

Not applicable. Our study used available summary statistics.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

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

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, D., Pan, J., Deng, A. et al. Mendelian randomization study of urolithiasis: exploration of risk factors using human blood metabolites. BMC Urol 24, 182 (2024). https://doi.org/10.1186/s12894-024-01568-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12894-024-01568-8

Keywords