Genome-wide Association Study for Vitamin D Levels Reveals 69 Independent Loci

We aimed to increase our understanding of the genetic determinants of vitamin D levels by undertaking a large-scale genome-wide association study (GWAS) of serum 25 hydroxyvitamin D (25OHD). To do so, we used imputed genotypes from 401,460 white British UK Biobank participants with available 25OHD levels, retaining single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) > 0.1% and imputation quality score > 0.3. We performed a linear mixed model GWAS on standardized log-transformed 25OHD, adjusting for age, sex, season of measurement, and vitamin D supplementation. These results were combined with those from a previous GWAS including 42,274 Europeans. In silico functional follow-up of the GWAS results was undertaken to identify enrichment in gene sets, pathways, and expression in tissues, and to investigate the partitioned heritability of 25OHD and its shared heritability with other traits. Using this approach, the SNP heritability of 25OHD was estimated to 16.1%. 138 conditionally independent SNPs were detected (p value < 6.6 × 10^-9) among which 53 had MAF < 5%. Single variant association signals mapped to 69 distinct loci, among which 63 were previously unreported. We identified enrichment in hepatic and lipid metabolism gene pathways and enriched expression of the 25OHD genes in liver, skin, and gastrointestinal tissues. We observed partially shared heritability between 25OHD and socio-economic traits, a feature which may be mediated through time spent outdoors. Therefore, through a large 25OHD GWAS, we identified 63 loci that underline the contribution of genes outside the vitamin D canonical metabolic pathway to the genetic architecture of 25OHD.