The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis

Abstract

Background: Prior research suggested a differential association of 25-hydroxyvitamin D (25(OH)D) metabolites with type 2 diabetes (T2D), with total 25(OH)D and 25(OH)D(3)inversely associated with T2D, but the epimeric form (C3-epi-25(OH)D-3) positively associated with T2D. Whether or not these observational associations are causal remains uncertain. We aimed to examine the potential causality of these associations using Mendelian randomisation (MR) analysis. Methods: and findings We performed a meta-analysis of genome-wide association studies for total 25(OH)D (N =120,618), 25(OH)D-3(N =40,562), and C3-epi-25(OH)D-3(N =40,562) in participants of European descent (European Prospective Investigation into Cancer and Nutrition [EPIC]-InterAct study, EPIC-Norfolk study, EPIC-CVD study, Ely study, and the SUNLIGHT consortium). We identified genetic variants for MR analysis to investigate the causal association of the 25(OH)D metabolites with T2D (including 80,983 T2D cases and 842,909 non-cases). We also estimated the observational association of 25(OH)D metabolites with T2D by performing random effects meta-analysis of results from previous studies and results from the EPIC-InterAct study. We identified 10 genetic loci associated with total 25(OH)D, 7 loci associated with 25(OH)D(3)and 3 loci associated with C3-epi-25(OH)D-3. Based on the meta-analysis of observational studies, each 1-standard deviation (SD) higher level of 25(OH)D was associated with a 20% lower risk of T2D (relative risk [RR]: 0.80; 95% CI 0.77, 0.84;p <0.001), but a genetically predicted 1-SD increase in 25(OH)D was not significantly associated with T2D (odds ratio [OR]: 0.96; 95% CI 0.89, 1.03;p =0.23); this result was consistent across sensitivity analyses. In EPIC-InterAct, 25(OH)D-3(per 1-SD) was associated with a lower risk of T2D (RR: 0.81; 95% CI 0.77, 0.86;p <0.001), while C3-epi-25(OH)D-3(above versus below lower limit of quantification) was positively associated with T2D (RR: 1.12; 95% CI 1.03, 1.22;p =0.006), but neither 25(OH)D-3(OR: 0.97; 95% CI 0.93, 1.01;p =0.14) nor C3-epi-25(OH)D-3(OR: 0.98; 95% CI 0.93, 1.04;p =0.53) was causally associated with T2D risk in the MR analysis. Main limitations include the lack of a non-linear MR analysis and of the generalisability of the current findings from European populations to other populations of different ethnicities. Conclusions: Our study found discordant associations of biochemically measured and genetically predicted differences in blood 25(OH)D with T2D risk. The findings based on MR analysis in a large sample of European ancestry do not support a causal association of total 25(OH)D or 25(OH)D metabolites with T2D and argue against the use of vitamin D supplementation for the prevention of T2D. Author summaryWhy was this study done? There is ongoing uncertainty on whether the body's vitamin D status indicated by blood 25-hydroxyvitamin D (25(OH)D) is relevant to the prevention of type 2 diabetes. There are conflicting findings from observational studies and a limited number of randomised controlled trials. Prior research did not distinguish between different metabolites of total 25(OH)D, including 25(OH)D3 and C3-epi-25(OH)D3, an isomer of 25(OH)D3. It is not clear whether the associations of 25(OH)D metabolites with type 2 diabetes are causal, with conflicting findings from prior Mendelian randomisation studies of total 25(OH)D and no previous Mendelian randomisation studies appraising 25(OH)D metabolites. What did the researchers do and find? The current research compared observational estimates of the association between 25(OH)D metabolites and incident type 2 diabetes with Mendelian randomisation estimates based on genetic instruments. Using multiple data sources, we performed genome-wide association studies among 120,618 individuals for total 25(OH)D, and among 40,562 individuals for the other vitamin D metabolites. Among participants of European descent, 10 genetic loci were identified for total 25(OH)D, 7 loci for 25(OH)D3 and 3 loci for C3-epi-25(OH)D3. In meta-analysis of observational studies, we found that each 1-standard deviation higher level of total 25(OH)D was associated with 20% lower risk of type 2 diabetes. The result was similar for 25(OH)D3, but for C3-epi-25(OH)D3, a positive association with type 2 diabetes was found. With up to 80,983 type 2 diabetes cases and 842,909 controls, we assessed the association of genetically predicted differences in total 25(OH)D and its metabolites with type 2 diabetes. Neither genetically predicted higher total 25(OH)D level nor genetically predicted higher levels of 25(OH)D metabolites were significantly associated with type 2 diabetes. What do these findings mean? There were conflicting findings for a link with type 2 diabetes for the observational analysis of biochemically measured 25(OH)D metabolites versus the genetically predicted levels of these metabolites. The null findings based on Mendelian randomisation analysis indicate that blood levels of 25(OH)D or its metabolites are not likely to be causal factors for the development of type 2 diabetes. The current findings together with other evidence from randomised controlled trials do not support the use of vitamin D supplementation for the prevention of type 2 diabetes.