Previous observational studies have revealed the association between serum uric acid and 25-hydroxyvitamin D. However, the causality and the direction of the associations remain unknown. Thus, we performed a two-sample bidirectional Mendelian Randomization (MR) analysis to investigate the causal association between uric acid and 25-hydroxyvitamin D and to determine the direction of the association.
Based on the summary-level GWAS data from large genome-wide association studies, several steps were taken in our analysis to select eligible single-nucleotide polymorphisms (SNPs), which were strongly related to exposure as the instrumental variables. We used different analytical methods, such as inverse-variance weighting (IVW) method, weighted median, MR-Egger regression, and weighted mode method, to make our result more robust and reliable. The IVW method was used as the primary analysis. The Cochran’s Q test, MR-Egger intercept test, MR-PRESSO method, and “leave-one-out” sensitivity analysis was performed to evaluate the heterogeneities, horizontal pleiotropy, and robustness of the results. MR analyses were also conducted using genetic risk scores (GRS) as instrumental variables in both directions by using the same summary-level GWAS data.
Our two-sample MR analysis suggested a causal association of genetically predicted uric acid on 25-hydroxyvitamin D [IVW method: β(SE), −0.0352(0.0149);
Serum uric acid may causally affect the 25-hydroxyvitamin D levels, whereas the causal role of 25-hydroxyvitamin D on uric acid was not supported in our MR analysis. Our findings suggest that increased levels of uric acid should prompt investigation for vitamin D deficiency.