#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

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


Authors: Ju-Sheng Zheng aff001;  Jian’an Luan aff001;  Eleni Sofianopoulou aff003;  Stephen J. Sharp aff001;  Felix R. Day aff001;  Fumiaki Imamura aff001;  Thomas E. Gundersen aff005;  Luca A. Lotta aff001;  Ivonne Sluijs aff006;  Isobel D. Stewart aff001;  Rupal L. Shah aff001;  Yvonne T. van der Schouw aff006;  Eleanor Wheeler aff001;  Eva Ardanaz aff007;  Heiner Boeing aff010;  Miren Dorronsoro aff011;  Christina C. Dahm aff012;  Niki Dimou aff013;  Douae El-Fatouhi aff014;  Paul W. Franks aff015;  Guy Fagherazzi aff014;  Sara Grioni aff017;  José María Huerta aff009;  Alicia K. Heath aff019;  Louise Hansen aff020;  Mazda Jenab aff013;  Paula Jakszyn aff021;  Rudolf Kaaks aff023;  Tilman Kühn aff023;  Kay-Tee Khaw aff024;  Nasser Laouali aff014;  Giovanna Masala aff025;  Peter M. Nilsson aff015;  Kim Overvad aff012;  Anja Olsen aff020;  Salvatore Panico aff027;  J. Ramón Quirós aff028;  Olov Rolandsson aff029;  Miguel Rodríguez-Barranco aff009;  Carlotta Sacerdote aff032;  Annemieke M. W. Spijkerman aff033;  Tammy Y. N. Tong aff034;  Rosario Tumino aff035;  Konstantinos K. Tsilidis aff019;  John Danesh aff003;  Elio Riboli aff019;  Adam S. Butterworth aff003;  Claudia Langenberg aff001;  Nita G. Forouhi aff001;  Nicholas J. Wareham aff001
Authors place of work: MRC Epidemiology Unit, University of Cambridge, Cambridge, United Kingdom aff001;  Westlake Laboratory of Life Sciences and Biomedicine, Key Laboratory of Growth Regulation and Translational Research of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, China aff002;  MRC/BHF Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge, United Kingdom aff003;  National Institute for Health Research Blood and Transplant Research Unit in Donor Health and Genomics, Department of Public Health and Primary Care, University of Cambridge, Cambridge, United Kingdom aff004;  VITAS, Oslo, Norway aff005;  Julius Center for Health Sciences and Primary Care, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands aff006;  Navarra Public Health Institute, Pamplona, Spain aff007;  Navarra Institute for Health Research (IdiSNA), Pamplona, Spain aff008;  CIBER Epidemiology and Public Health (CIBERESP), Madrid, Spain aff009;  Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Germany aff010;  Public Health Division of Gipuzkoa, San Sebastian, Spain aff011;  Department of Public Health, Aarhus University, Aarhus, Denmark aff012;  International Agency for Research on Cancer, Lyon, France aff013;  Center of Research in Epidemiology and Population Health, UMR 1018 Inserm, Institut Gustave Roussy, Paris South–Paris Saclay University, Villejuif, France aff014;  Department of Clinical Sciences, Lund University, Malmö, Sweden aff015;  Department of Population Health, Luxembourg Institute of Health, Strassen, Luxembourg aff016;  Epidemiology and Prevention Unit, Milan, Italy aff017;  Department of Epidemiology, Murcia Regional Health Council, Instituto Murciano de Investigación Biosanitaria Virgen de la Arrixaca, Murcia, Spain aff018;  Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, United Kingdom aff019;  Danish Cancer Society Research Center, Copenhagen, Denmark aff020;  Unit of Nutrition and Cancer, Cancer Epidemiology Research Program, Catalan Institute of Oncology–Institut d’Investigació Biomédica de Bellvitge, L’Hospitalet de Llobregat, Barcelona, Spain aff021;  Facultat Ciències Salut Blanquerna, Universitat Ramon Llull, Barcelona, Spain aff022;  Division of Cancer Epidemiology, German Cancer Research Center (DKFZ), Heidelberg, Germany aff023;  Department of Public Health and Primary Care, University of Cambridge, Cambridge, United Kingdom aff024;  Cancer Risk Factors and Life-Style Epidemiology Unit, Institute for Cancer Research, Prevention and Clinical Network (ISPRO), Florence, Italy aff025;  Department of Cardiology, Aalborg University Hospital, Aarhus, Denmark aff026;  Dipartimento di Medicina Clinica e Chirurgia, University of Naples Federico II, Naples, Italy aff027;  Public Health Directorate, Asturias, Spain aff028;  Family Medicine Division, Department of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden aff029;  Andalusian School of Public Health (EASP), Granada, Spain aff030;  Instituto de Investigación Biosanitaria de Granada, Universidad de Granada, Granada, Spain aff031;  Unit of Cancer Epidemiology, Città della Salute e della Scienza di Torino University Hospital–University of Turin and Center for Cancer Prevention (CPO), Torino, Italy aff032;  National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands aff033;  Cancer Epidemiology Unit, Nuffield Department of Population Health, University of Oxford, Oxford, United Kingdom aff034;  Azienda Sanitaria Provinciale, Ragusa, Italy aff035;  Department of Hygiene and Epidemiology, University of Ioannina School of Medicine, Ioannina, Greece aff036;  British Heart Foundation Cambridge Centre of Excellence, Division of Cardiovascular Medicine, Addenbrooke’s Hospital, Cambridge, United Kingdom aff037;  Department of Human Genetics, Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, United Kingdom aff038
Published in the journal: The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis. PLoS Med 17(10): e32767. doi:10.1371/journal.pmed.1003394
Category: Research Article
doi: https://doi.org/10.1371/journal.pmed.1003394

Summary

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)D3 inversely associated with T2D, but the epimeric form (C3-epi-25(OH)D3) 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)D3 (N = 40,562), and C3-epi-25(OH)D3 (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)D3 and 3 loci associated with C3-epi-25(OH)D3. 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)D3 (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)D3 (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)D3 (OR: 0.97; 95% CI 0.93, 1.01; p = 0.14) nor C3-epi-25(OH)D3 (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.

Keywords:

Genetics – Genetic loci – Genome-wide association studies – Metaanalysis – Metabolites – type 2 diabetes – vitamin D – Single nucleotide polymorphisms

Introduction

Adequate vitamin D status is crucial for maintaining bone homeostasis, and interest in its potential beneficial roles in cardiometabolic diseases, including type 2 diabetes (T2D), has increased recently. Prospective epidemiological studies have consistently reported an inverse association between circulating 25-hydroxyvitamin D (25(OH)D, a blood marker of vitamin D status) and T2D risk [1,2], but whether there is a causal relationship between vitamin D status and T2D remains uncertain. Older randomised controlled trials (RCTs) reported no beneficial effect of vitamin D supplementation on T2D risk [3,4], but they were limited by issues of sub-group or post hoc analyses, inadequate dose, or inability to separate the effect of vitamin D and calcium. The results of the VITAL study for the secondary endpoint of T2D are still awaited, with the primary endpoints being cancer and cardiovascular disease [5]. With T2D as a primary endpoint, the D2d RCT recently reported that vitamin D3 supplementation at a dose of 4,000 IU per day among people at high risk for T2D did not result in a lower risk of T2D than placebo [6]. This evidence argues against a benefit of vitamin D3 supplementation for the prevention of T2D, but there are still unresolved issues.

The D2d trial included participants with confirmed prediabetes, so it remains unclear whether vitamin D supplementation may be more effective at earlier stages in the natural history of T2D, before beta cell dysfunction. However, designing a trial to examine the effects of early supplementation poses several challenges including long follow-up, adherence, and cost. The mean baseline 25(OH)D level was 70 nmol/l in D2d, with the majority of participants vitamin D replete, and the study was not adequately powered to examine potential effects of supplementation among participants with lower levels of 25(OH)D. Moreover, studies of supplementation do not directly assess the causal relationship of long-term body vitamin D status as assessed by blood 25(OH)D metabolite levels, which is relevant biologically given that blood levels reflect both dietary source and endogenous synthesis. Genetic Mendelian randomisation (MR) analysis is a complementary approach to RCTs that, subject to several assumptions, can enable the estimation of causal associations using data from observational studies [7,8]. Given the lack of definitive evidence from RCTs and the challenges of implementing an ideal RCT, MR studies can be a useful supplementary tool.

Previous MR analyses of 25(OH)D and T2D risk reported conflicting results, with some MR studies suggesting no association [1,9,10] but others reporting inverse associations [11,12]. A prior constraint was the use of a genetic instrument that included a limited set of up to 4 single nucleotide polymorphisms (SNPs) in genes involved in vitamin D synthesis and metabolism. Specifically, synthesis pathways involve DHCR7 in 25(OH)D synthesis in the skin and CYP2R1 in hepatic 25-hydroxylation, while metabolism pathways involve GC (or DBP; encoding vitamin D binding protein) in 25(OH)D transport and CYP24A1 in 25(OH)D catabolism [13]. The inclusion of a greater number of genetic variants could generate a more powerful genetic instrumental variable for MR analysis. A further limitation to date has been the lack of a genetic instrument that can distinguish the 25(OH)D metabolites, including 25(OH)D3 and C3-epi-25(OH)D3, an isomer of 25(OH)D3 [14] that was positively associated with incident T2D [2]. There is no prior genome-wide association study (GWAS) of 25(OH)D3 to our knowledge, as past studies were restricted to analysing total 25(OH)D. Similarly, no GWAS to our knowledge has been performed for C3-epi-25(OH)D3.

Therefore, our objective was to assess the evidence for whether the association between total 25(OH)D and T2D is causal, using MR based on a genetic instrument derived from an updated meta-analysis of GWASs for total 25(OH)D. In addition, taking a similar approach, we also assessed the evidence to support any potentially causal associations with T2D for the vitamin D metabolites 25(OH)D3 and C3-epi-25(OH)D3 using a novel meta-analysis of GWASs for 25(OH)D3 and C3-epi-25(OH)D3.

Methods

Study design and populations

We adopted a multi-stage approach comprising 3 parts (defined in a prospective analysis plan; S1 Text): stage 1, GWASs of total 25(OH)D and vitamin D metabolites (25(OH)D3 and C3-epi-25(OH)D3); stage 2, MR analysis of T2D risk using GWAS-identified lead genetic variants; and stage 3, comparison of MR estimates with observational estimates on the association between vitamin D metabolites and T2D incidence.

At the GWAS stage, for total 25(OH)D, we performed a meta-analysis of GWASs including 120,618 participants of European origin from a number of cohorts as shown in Fig 1: the European Prospective Investigation into Cancer and Nutrition (EPIC)–InterAct study (n = 18,078) [15], EPIC-Norfolk study (n = 10,231) [16], EPIC-CVD study (n = 12,253) [17], Ely study (n = 690) [18], and publicly available GWAS summary statistics from the SUNLIGHT consortium (n = 79,366) [19]. We excluded participants duplicated among the EPIC-InterAct study, EPIC-Norfolk study, and EPIC-CVD study. We then used the UK Biobank dataset (n = 410,826) to replicate the findings from the above GWAS meta-analysis of total 25(OH)D. For meta-analysis of GWASs of both 25(OH)D3 and C3-epi-25(OH)D3, we included up to 40,562 participants of European origin from the EPIC-InterAct study (n = 18,078), EPIC-Norfolk study (n = 10,231), and EPIC-CVD study (n = 12,253).

Fig. 1. Design of the Mendelian randomisation study to estimate the causal association of 25-hydroxyvitamin D metabolites with type 2 diabetes.
Design of the Mendelian randomisation study to estimate the causal association of 25-hydroxyvitamin D metabolites with type 2 diabetes.
25(OH)D, 25-hydroxyvitamin D; GWAS, genome-wide association study.

At the MR stage, we estimated the association of GWAS-identified lead genetic variants with T2D risk, performing a meta-analysis including participants (80,983 T2D cases, and 842,909 non-cases) from the DIAMANTE study (55,005 T2D cases, 400,308 non-cases) [20], UK Biobank (24,758 T2D cases, 424,575 non-cases) [21], and EPIC-Norfolk study (additional cases not included in DIAMANTE study: 1,220 T2D cases, 18,026 non-cases) [16].

At stage 3, as a comparison to MR results, we estimated the observational association of total 25(OH)D, 25(OH)D3, and C3-epi-25(OH)D3 with T2D in the EPIC-InterAct study (8,331 T2D cases, 11,837 non-cases), adapted from the analysis in our prior publication [2]. For total 25(OH)D, we combined the results from the EPIC-InterAct study with those from a previous meta-analysis [1].

Ethics statement

Detailed description of each study is provided in S2 Text. All the studies included in the present analyses were approved by local ethical committees, and participants provided written informed consent. Specifically, EPIC-InterAct and EPIC-CVD were approved by local ethics committees in the participating countries and the institutional review board of the International Agency for Research on Cancer. EPIC-Norfolk was approved by the Norwich local ethics committee, and the Ely study was approved by the Ely local research ethics committee.

Circulating 25(OH)D measurement and genome-wide genotyping

In the EPIC-InterAct study, EPIC-Norfolk study, and EPIC-CVD study, plasma 25(OH)D metabolites were measured using liquid chromatography–tandem mass spectrometry (LC-MS/MS) at VITAS (Oslo, Norway; a reference lab in Europe with a Vitamin D External Quality Assessment Scheme certificate) [2]. In the Ely study, serum total 25(OH)D concentrations were measured by radioimmunoassay. The GWAS summary statistics from the SUNLIGHT consortium included data from 31 studies, within which radioimmunoassay, LC-MS/MS, enzyme-linked immunosorbent assay (ELISA), or chemiluminescence immunoassay was used for the measurement of the 25(OH)D metabolites [19]. In the UK Biobank, serum total 25(OH)D was measured using chemiluminescence immunoassay. Methods of genome-wide genotyping in each study are presented in S2 Text.

Statistical analysis

GWASs. For total 25(OH)D, we performed meta-analysis of 9 sets of GWASs (Table 1) including 4 in the EPIC-InterAct study (subcohort GWAS array [n = 3,844], subcohort core-exome array [n = 6,932], non-subcohort GWAS array [n = 3,188], and non-subcohort core-exome array [n = 4,114]), 1 in the EPIC-Norfolk study (n = 10,231), 2 in the EPIC-CVD study (subcohort [n = 887] and non-subcohort [n = 11,366]), 1 in the Ely study (n = 690), and 1 from the SUNLIGHT consortium (n = 79,366). For 25(OH)D3 and C3-epi-25(OH)D3 (as a binary variable: above versus below the lower limit of quantification [LLQ; 1 nmol/l]), the same datasets as total 25(OH)D were used except for the Ely study, where only total 25(OH)D was available. C3-epi-25(OH)D3 was treated as a binary variable as around half of values were below the LLQ (i.e., missing) in each individual cohort.

Tab. 1. Characteristics of the cohorts included in the genome-wide meta-analysis.
Characteristics of the cohorts included in the genome-wide meta-analysis.

For each participating cohort in the above GWAS discovery, standardised residuals of natural-log transformed 25(OH)D metabolites were calculated, adjusting for age, sex, BMI, season of blood collection, and study centre (where appropriate). Then, the GWAS was performed using linear regression with SNPTEST (v2.5.4) assuming an additive effect, adjusting for the first 10 genetic principal components of ancestry within each cohort. For the binary C3-epi-25(OH)D3 variable, the GWAS was performed using logistic regression (above versus below LLQ) with QUICKTEST (v6.5.2) assuming an additive effect, adjusting for age, sex, BMI, season of blood collection, study centre (where appropriate), and the first 10 genetic principal components of ancestry within each cohort.

For each of the 25(OH)D measures (total 25(OH)D, 25(OH)D3, and C3-epi-25(OH)D3), we performed fixed-effect inverse-variance-weighted (IVW) meta-analysis to combine our results with published GWAS summary statistics (from the SUNLIGHT consortium, only for total 25(OH)D) [19] using the software METAL [22]. The quality control thresholds were as follows: minor allele frequency ≥ 0.01, imputation info score ≥ 0.4, and p-value for Hardy–Weinberg equilibrium ≥ 1.0 × 10−6. Associated loci were identified using the conventional threshold for genome-wide statistical significance (p < 5 × 10−8). At each locus, the lead SNP was identified as the SNP with the lowest p-value within a 1 million base-pair window. To visualise the findings, we generated Manhattan plots and quantile–quantile plots for the genetic association using R package EasyStrata version 8.5, and regional association plots using LocusZoom software [23]. We used HaploReg v4.1 to explore annotations of the identified lead SNPs, including the nearest genes, eQTL, GRASP QTL, and previous GWAS hits [24].

We used UK Biobank to replicate our GWAS findings for total 25(OH)D. In the UK Biobank study, standardised residuals of natural-log transformed total 25(OH)D were calculated, adjusting for age, sex, BMI, season of blood collection, genotype chip, and aliquot number. GWAS analysis of total 25(OH)D was performed using a linear mixed model with BOLT-LMM, adjusting for the first 10 genetic principal components of ancestry.

Estimation of genetic correlations and variance explained

Genetic correlation of 25(OH)D variables with T2D and related glycaemic traits (fasting glucose, insulin, 2-hour glucose, homeostatic model assessment of insulin resistance [HOMA-IR], homeostatic model assessment of beta cell function [HOMA-B], and glycated haemoglobin [HbA1c]) was estimated with linkage disequilibrium score regression analysis using the meta-analysed GWAS summary statistics and publicly available datasets in the LD Hub platform [25,26]. Variance in 25(OH)D metabolites explained by the identified lead SNPs was estimated using linear regression models with individual-level data from the EPIC-Norfolk study and UK Biobank study. We also calculated the F-statistic in the EPIC-Norfolk study to evaluate the strength of the genetic instrument.

Observational analysis

For comparison to MR, we estimated the observational association between total 25(OH)D and T2D incidence by meta-analysing the results from a previous meta-analysis [1] together with results from the EPIC-InterAct study we previously published [2] (cases n = 15,611; non-cases n = 21,106). We also used the published effect estimate of T2D for 25(OH)D3 in the EPIC-InterAct study [2]. For the observational association between C3-epi-25(OH)D3 and T2D, we used results based on the EPIC-InterAct study, adapted from our prior analysis [2], to estimate the relative risk (RR) of T2D comparing those above versus below the LLQ. To minimise the possibility of residual confounding in the observational association between 25(OH)D and T2D by adiposity, we performed exploratory analyses examining the influence of additional adiposity-related covariates. Specifically, in the EPIC-InterAct study, we included a genetic risk score (GRS) for BMI (generated from 97 BMI-related genetic variants by summing up the number of risk alleles) [27], and quadratic terms for BMI and waist-to-hip ratio, to account for a potential non-linear association between adiposity and T2D.

MR analyses

We performed MR analysis to combine estimates of ‘SNP to 25(OH)D level’ and ‘SNP to T2D’ associations, to estimate the genetically predicted association of the 25(OH)D variable (either total 25(OH)D or 25(OH)D3) with T2D (odds ratio [OR] per 1-SD increase). Similarly, we estimated the association of genetically predicted higher C3-epi-25(OH)D3 level (above versus below the LLQ) with T2D risk. For each of the 25(OH)D variables, we used an IVW method, MR-Egger method, and weighted median method to pool the estimates from multiple SNPs [28,29]. We used MR-Egger regression to detect and adjust for potential unbalanced pleiotropy in the MR analysis, and we used the weighted median MR method to examine the robustness of the results and highlight the results if significant heterogeneity of the associations among different genetic variants was observed. We used the effect estimate from the 25(OH)D GWAS discovery cohort, not UK Biobank, for the above MR analyses, as there were a large number of overlapping UK Biobank samples in the T2D GWASs. In response to peer reviewer comments, we performed analyses using several other MR methods to test the robustness of our results (S1 Text), including MR-PRESSO [30], MR-RAPS [31], and MRMix [32], and we further performed a multivariable MR analysis [7] that jointly estimated the causal association of the highly related traits. In addition, we calculated the statistical power of the MR analysis, which suggested that our MR had 90% power to detect a 6% difference in T2D risk per 1-SD change in 25(OH)D (α = 0.05, assuming that the genetic instrument explains 3.95% of the variance of the exposure).

MR assumptions are violated if there is horizontal pleiotropy, i.e., MR is valid when the genetic instrument is associated only with the exposure, not other variables that may be potential confounders. To assess the plausibility of this assumption, we examined the association of the GWAS-identified lead SNPs with metabolic markers and lifestyle and demographic factors, using the PhenoScanner tool (v2) [33]. We performed secondary analyses stratifying the genetic variants into several groups according to their biological roles: for total 25(OH)D and 25(OH)D3, lead SNPs were divided into those in the vitamin D synthesis pathway (CYP2R1 and NADSYN1/DHCR7), in the vitamin D metabolism pathway (GC, CYP24A1), and the others (AMDHD1, SEC23A, PADI1, CRCT1, UGT1A5, SULT2A1); for C3-epi-25(OH)D3, lead SNPs were divided into 2 groups: the SDR9C7 group and the others (which overlapped with 25(OH)D genes).

We also performed MR analysis to examine the genetically predicted association of the 25(OH)D metabolites with glycaemic traits, including HOMA-IR, HOMA-B, HbA1c, fasting insulin, fasting glucose, and 2-hour glucose. R version 3.4.3 and Stata version 14.2 (StataCorp) were used for the above statistical analyses.

This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 STROBE Checklist).

Results

Total 25(OH)D and T2D

The average 25(OH)D concentrations (SD) of the cohorts ranged from 38.0 (17.2) nmol/l in EPIC-InterAct to 58.6 (23.9) nmol/l in the Ely study (Table 1). For total 25(OH)D level (Fig 2; Table 2), in the meta-analysis of 9 sets of GWASs among 120,618 participants, we identified 10 genetic loci. Of these 10 loci, we confirmed 6 previously discovered loci at GC (rs3755967, p = 2.48 × 10−465), CYP2R1 (rs116970203, p = 1.19 × 10−64), NADSYN1/DHCR7 (rs12785878, p = 5.60 × 10−87), AMDHD1 (rs3213737, p = 2.05 × 10−19), SEC23A (rs8018720, p = 1.46 × 10−10), and CYP24A1 (rs17216707, p = 1.61 × 10−29). Additionally, we identified 4 novel loci at PADI1 (rs11203339, p = 4.64 × 10−8), CRCT1 (rs7529325, p = 2.09 × 10−9), UGT1A5 (rs17862870, p = 5.57 × 10−9), and SULT2A1 (rs9304669, p = 4.53 × 10−8) (Table 2; S1S3 Figs). The association of total 25(OH)D SNPs with each 25(OH)D metabolite is shown in S1 Table. In the UK Biobank, we replicated associations of all the above 10 SNPs, with p-values < 10−33 (S2 Table).

Fig. 2. Genome-wide association of circulating 25-hydroxyvitamin D metabolites: Manhattan plot and quantile–quantile plot of all genetic variants from the meta-analysis.
Genome-wide association of circulating 25-hydroxyvitamin D metabolites: Manhattan plot and quantile–quantile plot of all genetic variants from the meta-analysis.
Manhattan plot and quantile–quantile plot for each of total 25-hydroxyvitamin D (A), 25-hydroxyvitamin D3 (B), and C3-epi-25-hydroxyvitamin D3 (C). Manhattan plot (left): SNPs are plotted on the x-axis according to their positions on each chromosome. The red line indicates the threshold for genome-wide significance (p < 5 × 10−8). Blue points represent SNPs in the ±100-kb region of a top hit. Loci are annotated with the gene names closest to the SNPs with lowest p-values (red indicates novel loci identified in the present study). Quantile–quantile plot (right): The x-axis shows the expected −log10 p-values, and the y-axis shows the observed −log10 p-values. Each SNP is plotted as a red dot, and the black line indicates the null hypothesis of no true association.
Tab. 2. Genetic loci identified in the genome-wide analyses for circulating 25(OH)D metabolites.
Genetic loci identified in the genome-wide analyses for circulating 25(OH)D metabolites.

Total 25(OH)D was not genetically correlated with T2D or glycaemic traits including HOMA-IR, HOMA-B, HbA1c, fasting insulin, fasting glucose, and 2-hour glucose (S3 Table). In the EPIC-Norfolk study, the variance explained by the 10 lead SNPs for total 25(OH)D was 3.95% (with F-statistic = 43.1), which was slightly higher than that explained (3.66%) by the 6 previously known loci (GC, CYP2R1, NADSYN1, CYP24A1, AMDHD1, and SEC23A). In the UK Biobank study, the variance explained by the 10 lead SNPs and the 6 previously known loci was 3.28% and 3.08%, respectively.

The observational analyses suggested that higher levels of 25(OH)D were associated with lower risk of T2D (per 1-SD RR: 0.80; 95% CI 0.77, 0.84; p < 0.001) (Fig 3). Further adjustment for BMI GRS (weighted or unweighted) and quadratic terms for BMI and waist-to-hip ratio had very little effect on the results (S4 Table).

Fig. 3. Association of 25(OH)D metabolites with type 2 diabetes from observational and MR analyses.
Association of 25(OH)D metabolites with type 2 diabetes from observational and MR analyses.
Estimates (95% CIs) were scaled to represent RRs from observational analyses or ORs from MR per 1–standard deviation difference in each 25(OH)D metabolite, except for the binary C3-epi-25(OH)D3 variable (above versus below the lower limit of quantification). For the MR sensitivity analysis for total 25(OH)D, the 4 prior known genes are GC, CYP2R1, NADSYN1/DHCR7, and CYP24A1; the 25(OH)D synthesis pathway genes are CYP2R1 and NADSYN1/DHCR7; the 25(OH)D metabolism genes are GC and CYP24A1; and the 6 recently identified genes are PADI1, CRCT1, UGT1A5, AMDHD1, SEC23A, and SULT2A1. For each analysis, 3 MR methods were used, including the IVW method, MR-Egger method, and weighted median method. For MR analyses based on 2 SNPs, the MR-Egger and weighted median methods could not be applied. For the observational estimate, the total number of type 2 diabetes cases and non-cases was 15,611 and 21,106, respectively, for total 25(OH)D, and 8,331 and 11,837, respectively, for 25(OH)D3 (or C3-epi-25(OH)D3). For the MR estimate of all 3 vitamin D variables, the total number of cases and non-cases was 80,983 and 842,909, respectively. 25(OH)D, 25-hydroxyvitamin D; CI, confidence interval; IVW, inverse-variance-weighted; MR, Mendelian randomisation; OR, odds ratio; RR, relative risk; SNP, single nucleotide polymorphism.

In the MR analysis (Fig 3), we did not find evidence that genetically predicted higher total 25(OH)D was associated with T2D, with per 1-SD ORs of 0.96 (95% CI 0.89, 1.03; p = 0.23), 1.00 (95% CI 0.91, 1.10; p = 0.99), and 0.98 (95% CI 0.90, 1.06; p = 0.58) for the IVW, MR-Egger, and weighted median methods, respectively. Individual 25(OH)D genetic variants were not significantly associated with metabolic markers or lifestyle or demographic factors (S4 Fig). Testing the MR-Egger intercept (beta = −0.005; 95% CI −0.011, 0.001; p = 0.15) did not yield evidence of directional pleiotropy. MR-PRESSO, MR-RAPS, MRMix, and multivariable MR produced similar null results (S5 Fig). MR analyses suggested that total 25(OH)D was not significantly associated with any glycaemic trait (S5 Table).

Individual 25(OH)D metabolites and T2D

For 25(OH)D3 (Fig 2; Table 2), we identified 7 novel genetic loci, of which one was a unique locus at SHQ1 (rs13084927, p = 1.94 × 10−8), while the other 6 were already known to affect total 25(OH)D: GC (rs4588, p = 6.55 × 10−266), CYP2R1 (rs116970203, p = 3.22 × 10−60), NADSYN1/DHCR7 (rs28364617, p = 4.08 × 10−59), AMDHD1 (rs3819817, p = 3.59 × 10−16), SULT2A1 (rs9304669, p = 1.27 × 10−8), and CYP24A1 (rs17216707, p = 1.09 × 10−15) (S6S8 Figs).

For C3-epi-25(OH)D3, we identified 2 loci that overlapped with total 25(OH)D at GC (rs4588, p = 2.48 × 10−32) and NADSYN1/DHCR7 (rs28364617, p = 1.95 × 10−15), and 1 unique locus at SDR9C7 (rs11172066, p = 8.08 × 10−15) (Table 2; S9S11 Figs).

Neither 25(OH)D3 nor C3-epi-25(OH)D3 was genetically correlated with T2D or glycaemic traits including HOMA-IR, HOMA-B, HbA1c, fasting insulin, fasting glucose, and 2-hour glucose (S4 Table). In the EPIC-Norfolk study, the variance explained by the lead SNPs was 4.58% for 25(OH)D3 (with F-statistic = 71.2) and 0.41% for C3-epi-25(OH)D3 (binary variable) (with F-statistic = 15.2).

In the non-genetic observational analyses, circulating 25(OH)D3 was inversely associated with T2D (per 1-SD RR: 0.81; 95% CI 0.77, 0.86; p < 0.001), and C3-epi-25(OH)D3 was positively associated with T2D (above versus below LLQ RR: 1.12; 95% CI 1.03, 1.22; p = 0.006) (Fig 3). Further adjustment for BMI GRS (weighted or unweighted) and quadratic terms for BMI and waist-to-hip ratio did not change the estimate for 25(OH)D3, while slightly strengthening the positive association for C3-epi-25(OH)D3 (S5 Table).

There was no evidence of pleiotropic associations of the genetic variants with metabolic markers or lifestyle or demographic factors (S12 and S13 Figs). Genetically predicted increases in 25(OH)D3 or high levels of C3-epi-25(OH)D3 (binary) were not associated with T2D (Fig 3). The findings did not change substantially in the analyses using other MR methods (S5 Fig) and in secondary analyses stratifying SNPs into different groups (S14 Fig). Testing the MR-Egger intercept did not yield evidence of directional pleiotropy for 25(OH)D3 (beta = −0.008; 95% CI −0.018, 0.002; p = 0.16) or C3-epi-25(OH)D3 (beta = −0.021; 95% CI −0.073, 0.03; p = 0.65). None of the genetically predicted variations of vitamin D metabolites were significantly associated with the glycaemic outcomes (S5 Table).

Discussion

Our updated meta-analysis of GWASs for total 25(OH)D levels using data from 120,618 European-descent participants identified 4 novel genetic loci, in addition to 6 loci previously described [19]. The 4 novel loci were replicated in the UK Biobank study. With GWAS-identified lead SNPs as a genetic instrument, our MR analysis did not find evidence that total 25(OH)D was causally associated with T2D. To the best of our knowledge, the present study is the first GWAS discovery for 25(OH)D3, the major metabolite of circulating total 25(OH)D, identifying 7 genetic loci, with 1 unique locus and 6 loci overlapping with those for total 25(OH)D. Similarly, for C3-epi-25(OH)D3, we identified 1 unique locus associated with C3-epi-25(OH)D3. In the MR analysis, we did not find evidence supporting a causal association between T2D and 25(OH)D3 or C3-epi-25(OH)D3.

Our current finding of lack of evidence for a causal association of 25(OH)D with T2D in MR analysis is at variance with the strong inverse observational association between 25(OH)D and T2D [1,34]. This discrepancy was substantial, with non-overlapping confidence intervals (RR 0.80 [95% CI 0.77, 0.84] and 0.96 [95% CI 0.89, 1.03], respectively, for the observational and MR findings), and was not explained by detailed adjustment for the confounding effect of adiposity in the observational analysis. The reasons for the difference between the observational and MR findings remain unclear, but it may be that adjustment for adiposity, diet, and physical activity using a single imprecise measure of these variables at baseline only partially reduced the confounding effects of these variables. However, the current lack of evidence for a causal association between 25(OH)D levels and T2D is consistent with recent RCT evidence (D2d trial) of a lack of benefit from vitamin D supplementation for the prevention of T2D [6]. The D2d trial was conducted among individuals with prediabetes and thus at high risk for developing T2D [6]; prior RCTs among postmenopausal women or elderly people also had null findings [3,4], and the awaited results from the VITAL trial will be further informative.

Several prior MR analyses found inconsistent results [1,912]. In a MR analysis among 28,144 T2D cases and 76,344 controls, we previously used 4 prior GWAS-discovered genetic variants (at GC, CYP2R1, DHCR7, and CYP24A1) as genetic instruments, and found a null association between 25(OH)D and T2D (OR 1.01 [95% CI 0.75, 1.36] per 1-SD reduction in 25(OH)D) [1]. When we restricted our analysis to the 2 synthesis-related SNPs (at DHCR7 and CYP2R1), we still had null findings [1], but this differed from other studies [11,12]. Afzal et al. examined variants in the 2 synthesis-related genes DHCR7 and CYP2R1 among 96,423 white Danish adults, reporting that per 20-nmol/l genetically determined reduction in plasma 25(OH)D there was an allelic effect for DHCR7, with an OR of 1.51 (95% CI 0.98, 2.33), but not for CYP2R1 (OR 1.02 [95% CI 0.75–1.37]) [11]. The authors acknowledged that their results were weak and generated a hypothesis for a possible causal inverse association with endogenously synthesised 25(OH)D. Lu et al. detected a significant causal protective effect on T2D risk using 2 synthesis SNPs related to DHCR7 and CYP2R1 only when including Chinese and European populations in a meta-analysis [12]. The possibility of uncorrected population stratification by different ethnic groups [35] cannot be excluded, though the authors made attempts to reduce this possibility by analysing area-specific estimates and combining them using inverse-variance weighting. Notably, a recent Chinese study did not find a causal association between 25(OH)D and T2D [10].

Described for the first time, to the best of our knowledge, we found a genome-wide significant signal for 25(OH)D3 at SHQ1, which encodes H/ACA ribonucleoprotein assembly factor and has functions in the processing of ribosomal RNAs, modification of spliceosomal small nuclear RNAs, and stabilisation of telomerase [36]. In addition, we identified a novel locus at SDR9C7 for C3-epi-25(OH)D3. SDR9C7 encodes short chain dehydrogenase, whose relationship to the epimerisation of 25(OH)D3 to C3-epi-25(OH)D3 or other metabolites is unknown. Since the biological role of the epimerase is yet to be discovered [37], the present GWAS results may provide important insights for future investigation of mechanisms.

Although 25(OH)D3 is usually the major component of total 25(OH)D [2], the causal association of 25(OH)D3, on its own, with T2D risk has not previously been evaluated in MR analysis. Prior research mainly focused on total 25(OH)D as a biomarker of vitamin D status, which included a combination of 25(OH)D2 and 25(OH)D3, and, less frequently, C3-epi-25(OH)D3, the epimeric form of 25(OH)D3, depending on the assay methods. Blood 25(OH)D3 is derived both from diet and from biosynthesis in the skin upon exposure to sunlight, and 25(OH)D2 is mainly derived from diet, while C3-epi-25(OH)D3 is a metabolite of 25(OH)D3 via a C3 epimerisation process [37]. Traditional high-performance liquid chromatography, liquid chromatography–mass spectrometry, ELISA, and chemiluminescent immunoassay methods do not distinguish C3-epi-25(OH)D3 from 25(OH)D3 and thus include it within the definition of 25(OH)D3 [38]. Our recent study [2] suggested that blood 25(OH)D3 had an inverse association with T2D risk, while C3-epi-25(OH)D3 was positively associated. Therefore, these current findings are an important extension to previous research, by including MR analysis for T2D for total 25(OH)D as well as for 25(OH)D3 and C3-epi-25(OH)D3 separately, as the results for total 25(OH)D may be potentially confounded by other vitamin D metabolites.

There are several strengths of this study. To the best of our knowledge, this is the largest GWAS meta-analysis to date of total 25(OH)D (N = 120,618), and used the most comprehensive genetic instrument for the MR estimate for 25(OH)D and T2D risk (based on 10 SNPs versus the 4 or 2 used previously [1,11,12]). We replicated our total 25(OH)D GWAS results in the UK Biobank study, but were unable to evaluate genetic associations for 25(OH)D metabolites since these were not available in UK Biobank. In our MR analysis, we included 80,983 T2D cases, a larger sample than in the largest previous MR study (58,312 T2D cases) [12]. We performed a novel GWAS on individual vitamin D metabolites—25(OH)D3 and C3-epi-25(OH)D3—and used results from these in the MR analysis.

A limitation of this study is the lack of generalisability of our results from European populations to other populations of different ethnicities. Another limitation is that we combined studies with different study designs to maximise sample size and power for both MR analysis and observational analysis, even though definitions of endpoints (e.g., contributing to the degree of misclassification of T2D status) and exposure measurements (e.g., different assays for 25(OH)D measurement) were not standardised across these studies. In addition, we tested a linear causal association using the MR approach, but we cannot rule out a potential weak causal association among people with vitamin D deficiency or insufficiency. Furthermore, although we did not find convincing evidence for an association of these variants with a variety of metabolic markers and lifestyle and demographic factors, we cannot rule out the possibility that horizontal pleiotropy exists for the genetic variants used in the MR analysis. Finally, use of a binary exposure (C3-epi-25(OH)D3) in MR analysis can potentially lead to violation of the exclusion restriction assumption: The genetic variant can influence the outcome via the continuous risk factor even if the binary exposure does not change [39].

The public health implication of the present study, together with prior evidence of null findings from RCTs of vitamin D supplementation, is that focusing on population-wide supplementation to raise blood vitamin D levels is not likely to be an effective strategy for the prevention of T2D in European populations. However, our current limitations also indicate there are unresolved issues with important implications for further research. Next steps to address these issues include, but are not limited to, investigating the reasons for the discrepancy between the observational and MR findings, which likely reflect residual confounding as discussed above; including populations from different ethnic groups with varying distributions of 25(OH)D levels; and further increasing sample size and hence statistical power for conducting MR analysis of non-linear associations. The last of these would help to address whether there may be effects specifically in those who have hypovitaminosis D or are vitamin D deficient.

In summary, the present findings using MR do not support a causal association of 25(OH)D and related metabolites with T2D, despite the strong association seen in observational studies. The totality of the available evidence from our study and RCTs to date does not justify the use of vitamin D supplementation for the prevention of T2D.

Supporting information

S1 STROBE Checklist [docx]
STROBE checklist.

S1 Fig [tif]
Regional plots of the genome-wide significant loci for total 25-hydroxyvitamin D.

S2 Fig [tif]
QQ plots by GWAS performed for total 25-hydroxyvitamin D.

S3 Fig [tif]
Association of novel identified genetic variants with total 25-hydroxyvitamin D in each GWAS.

S4 Fig [tif]
Association of total 25-hydroxyvitamin D genetic variants with metabolic markers and lifestyle and demographic factors.

S5 Fig [oh]
Forest plots for the Mendelian randomisation analysis of 25-hydroxyvitamin D metabolites and type 2 diabetes.

S6 Fig [tif]
Regional plots of the genome-wide significant loci for 25-hydroxyvitamin D.

S7 Fig [tif]
QQ plots by GWAS performed for total 25-hydroxyvitamin D.

S8 Fig [tif]
Association of novel identified genetic variants with 25-hydroxyvitamin D in each GWAS.

S9 Fig [tif]
Regional plots of the genome-wide significant loci for C3-epi-25-hydroxyvitamin D.

S10 Fig [tif]
QQ plots by GWAS performed for C3-epi-25-hydroxyvitamin D.

S11 Fig [tif]
Association of novel identified genetic variants with C3-epi-25-hydroxyvitamin D in each GWAS.

S12 Fig [tif]
Association of 25-hydroxyvitamin D genetic variants with metabolic markers and lifestyle and demographic factors.

S13 Fig [tif]
Association of C3-epi-25-hydroxyvitamin D genetic variants with metabolic markers and lifestyle and demographic factors.

S14 Fig [mr]
Mendelian randomisation analysis of 25-hydroxyvitamin D metabolites and type 2 diabetes with stratification of single nucleotide polymorphisms.

S1 Table [docx]
Association of GWAS-identified genetic variants with each of the 25-hydroxyvitamin D metabolites.

S2 Table [docx]
Replication of the GWAS-identified variants for total 25-hydroxyvitamin D in the UK Biobank study.

S3 Table [docx]
Genetic correlation of 25-hydroxyvitamin D metabolites with type 2 diabetes and glycaemic traits.

S4 Table [docx]
Observational estimates of the association of 25-hydroxyvitamin D metabolites with incident type 2 diabetes in the EPIC-InterAct study.

S5 Table [docx]
Mendelian randomisation analysis of 25-hydroxyvitamin D metabolites and glycaemic traits.

S1 Text [docx]
Study protocol.

S2 Text [docx]
Description of the participating studies in the present genome-wide association study for vitamin D metabolites.


Zdroje

1. Ye Z, Sharp SJ, Burgess S, Scott RA, Imamura F, InterAct Consortium, et al. Association between circulating 25-hydroxyvitamin D and incident type 2 diabetes: a Mendelian randomisation study. Lancet Diabetes Endocrinol. 2015;3:35–42. doi: 10.1016/S2213-8587(14)70184-6 25281353

2. Zheng J-S, Imamura F, Sharp SJ, van der Schouw YT, Sluijs I, Gundersen TE, et al. Association of plasma vitamin D metabolites with incident type 2 diabetes: EPIC-InterAct case-cohort study. J Clin Endocrinol Metab. 2019;104:1293–303. doi: 10.1210/jc.2018-01522 30418614

3. Avenell A, Cook JA, MacLennan GS, McPherson GC. Vitamin D supplementation and type 2 diabetes: a substudy of a randomised placebo-controlled trial in older people (RECORD trial, ISRCTN 51647438). Age Ageing. 2009;38(5):606–9. doi: 10.1093/ageing/afp109 19617604

4. De Boer IH, Tinker LF, Connelly S, Curb JD, Howard BV, Kestenbaum B, et al. Calcium plus vitamin D supplementation and the risk of incident diabetes in the women’s health initiative. Diabetes Care. 2008;31:701–7. doi: 10.2337/dc07-1829 18235052

5. Manson JE, Cook NR, Lee I-M, Christen W, Bassuk SS, Mora S, et al. Vitamin D supplements and prevention of cancer and cardiovascular disease. N Engl J Med. 2019;380:33–44. doi: 10.1056/NEJMoa1809944 30415629

6. Pittas AG, Dawson-Hughes B, Sheehan P, Ware JH, Knowler WC, Aroda VR, et al. Vitamin D supplementation and prevention of type 2 diabetes. N Engl J Med. 2019;380:23–32. doi: 10.1056/NEJMoa1811403 30415637

7. Burgess S, Thompson SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181:251–60. doi: 10.1093/aje/kwu283 25632051

8. Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA. 2017;318:1925. doi: 10.1001/jama.2017.17219 29164242

9. Buijsse B, Boeing H, Hirche F, Weikert C, Schulze MB, Gottschald M, et al. Plasma 25-hydroxyvitamin D and its genetic determinants in relation to incident type 2 diabetes: a prospective case-cohort study. Eur J Epidemiol. 2013;28:743–52. doi: 10.1007/s10654-013-9844-5 24002339

10. Wang N, Wang C, Chen X, Wan H, Chen Y, Chen C, et al. Vitamin D, prediabetes and type 2 diabetes: bidirectional Mendelian randomization analysis. Eur J Nutr. 2019;59:1379–88. doi: 10.1007/s00394-019-01990-x 31076857

11. Afzal S, Brøndum-Jacobsen P, Bojesen SE, Nordestgaard BG. Vitamin D concentration, obesity, and risk of diabetes: a Mendelian randomisation study. Lancet Diabetes Endocrinol. 2014;2:298–306. doi: 10.1016/S2213-8587(13)70200-6 24703048

12. Lu L, Bennett DA, Millwood IY, Parish S, McCarthy MI, Mahajan A, et al. Association of vitamin D with risk of type 2 diabetes: a Mendelian randomisation study in European and Chinese adults. PLOS Med. 2018;15:e1002566. doi: 10.1371/journal.pmed.1002566 29718904

13. Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet. 2010;376:180–8. doi: 10.1016/S0140-6736(10)60588-0 20541252

14. Bailey D, Veljkovic K, Yazdanpanah M, Adeli K. Analytical measurement and clinical relevance of vitamin D3C3-epimer. Clin Biochem. 2013;46(3):190–6. doi: 10.1016/j.clinbiochem.2012.10.037 23153571

15. Consortium InterAct, Langenberg C, Sharp S, Forouhi NG, Franks PW, Schulze MB, et al. Design and cohort description of the InterAct Project: an examination of the interaction of genetic and lifestyle factors on the incidence of type 2 diabetes in the EPIC Study. Diabetologia. 2011;54:2272–82. doi: 10.1007/s00125-011-2182-9 21717116

16. Day N, Oakes S, Luben R, Khaw KT, Bingham S, Welch A, et al. EPIC-Norfolk: study design and characteristics of the cohort. European Prospective Investigation of Cancer. Br J Cancer. 1999;80(Suppl 1):95–103.

17. Danesh J, Saracci R, Berglund G, Feskens E, Overvad K, Panico S, et al. EPIC-Heart: the cardiovascular component of a prospective study of nutritional, lifestyle and biological factors in 520,000 middle-aged participants from 10 European countries. Eur J Epidemiol. 2007;22:129–41. doi: 10.1007/s10654-006-9096-8 17295097

18. Forouhi NG, Luan J, Hennings S, Wareham NJ. Incidence of type 2 diabetes in England and its association with baseline impaired fasting glucose: the Ely study 1990–2000. Diabet Med. 2007;24:200–7. doi: 10.1111/j.1464-5491.2007.02068.x 17257284

19. Jiang X, O’Reilly PF, Aschard H, Hsu Y-H, Richards JB, Dupuis J, et al. Genome-wide association study in 79,366 European-ancestry individuals informs the genetic architecture of 25-hydroxyvitamin D levels. Nat Commun. 2018;9:260. doi: 10.1038/s41467-017-02662-2 29343764

20. Mahajan A, Taliun D, Thurner M, Robertson NR, Torres JM, Rayner NW, et al. Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat Genet. 2018;50:1505–13. doi: 10.1038/s41588-018-0241-6 30297969

21. Collins R. What makes UK Biobank special? Lancet. 2012;379:1173–4. doi: 10.1016/S0140-6736(12)60404-8 22463865

22. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–1. doi: 10.1093/bioinformatics/btq340 20616382

23. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26:2336–7. doi: 10.1093/bioinformatics/btq419 20634204

24. Ward LD, Kellis M. HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016;44:D877–81. doi: 10.1093/nar/gkv1340 26657631

25. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47:1236–41. doi: 10.1038/ng.3406 26414676

26. Zheng J, Erzurumluoglu AM, Elsworth BL, Kemp JP, Howe L, Haycock PC, et al. LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics. 2017;33:272–9. doi: 10.1093/bioinformatics/btw613 27663502

27. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518:197–206 doi: 10.1038/nature14177 25673413

28. 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:304–14. doi: 10.1002/gepi.21965 27061298

29. 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 24114802

30. 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–8. doi: 10.1038/s41588-018-0099-7 29686387

31. Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Stat. 2020;3:1742–69. doi: 10.1214/19-AOS1866

32. Qi G, Chatterjee N. Mendelian randomization analysis using mixture models for robust and efficient estimation of causal effects. Nat Commun. 2019;10:1941. doi: 10.1038/s41467-019-09432-2 31028273

33. Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. 2016;32:3207–9. doi: 10.1093/bioinformatics/btw373 27318201

34. Song Y, Wang L, Pittas AG, Del Gobbo LC, Zhang C, Manson JE, et al. Blood 25-hydroxy vitamin D levels and incident type 2 diabetes: a meta-analysis of prospective studies. Diabetes Care. 2013;36:1422–8. doi: 10.2337/dc12-0962 23613602

35. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601. doi: 10.1136/bmj.k601 30002074

36. Grozdanov PN, Roy S, Kittur N, Meier UT. SHQ1 is required prior to NAF1 for assembly of H/ACA small nucleolar and telomerase RNPs. RNA. 2009;15:1188–97. doi: 10.1261/rna.1532109 19383767

37. Bikle DD. Vitamin D metabolism, mechanism of action, and clinical applications. Chem Biol. 2014;21(3):319–29. doi: 10.1016/j.chembiol.2013.12.016 24529992

38. Le Goff C, Cavalier E, Souberbielle JC, González-Antuña A, Delvin E. Measurement of circulating 25-hydroxyvitamin D: a historical review. Pract Lab Med. 2015;2:1–14. doi: 10.1016/j.plabm.2015.04.001 28932799

39. Burgess S, Labrecque JA. Mendelian randomization with a binary exposure variable: interpretation and presentation of causal estimates. Eur J Epidemiol. 2018;33:947–52. doi: 10.1007/s10654-018-0424-6 30039250


Článek vyšel v časopise

PLOS Medicine


2020 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#