|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: KSV DJB MM OR JV BDM EAM NJW TM TL WHO CC TJW MRJ JCW ADH EH. Performed the experiments: KSV DJB CL ET SP LTH JDC KHL ZD RL ARW KM LZ LMY JD MCK KJ VS BDM EAS DKH TL PK CC RJFL JCW ADH EH. Analyzed the data: KSV DJB CL ET SP LTH JDC KHL ZD RL ARW KM LZ LMY JD MCK KJ VS BDM EAS DKH TL PK CC RJFL JCW ADH EH. Contributed reagents/materials/analysis tools: KSV DJB SP LTH KM LMY MM JD MK NA OR JV SBK HM EI LB LL VS BDM KH EAS DKH TL PK CC WM RJFL CP MRJ JCW ADH EH. Wrote the first draft of the manuscript: KSV DJB EH. Contributed to the writing of the manuscript: KSV DJB EH. ICMJE criteria for authorship read and met: KSV DJB CL ET SP LTH JDC ZD RL DKH ARW KM LV LZ LMY MM JD MK MEK KJ NA OR JV KKL LF HM EI LB LL ML VS HC MD BDM KHH AP ALH EAS ET AJ NJW CO TMF SBK TDS JBR TL WHO PK CC WM CP RJFL TJW MRJ JCW ADH EH. Agree with manuscript results and conclusions: KSV DJB CL ET SP LTH JDC ZD RL DKH ARW KM LV LZ LMY MM JD MK MEK KJ NA OR JV KKL LF HM EI LB LL ML VS HC MD BDM KHH AP ALH EAS ET AJ NJW CO TMF SBK TDS JBR TL WHO PK CC WM CP RJFL TJW MRJ JCW ADH EH. Obtained the data: KSV DJB SP LTH ZD RL KM LV MM MK MCK KJ NA OR JV SBK LF HM EI LB LL ML VS HC MD BDM AP AH ET NJW CO TMF DKH JBR TL WHO PK CC WM CP ADH MRJ EH. Provided the administrative, technical, or material support: KSV DJB KM MK OR JV SBK LF HM EI LL VS TMF DKH TL CP EH. Supervised the study: JD OR JV SBK HC MD BDM EAS AJ ET TS TMF JBR TL TJW MRJ JCW ADH EH.
Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.
We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects.
Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n=123,864). Each 1 kg/m2 higher BMI was associated with 1.15% lower 25(OH)D (p=6.52×10−27). The BMI allele score was associated both with BMI (p=6.30×10−62) and 25(OH)D (−0.06% [95% CI −0.10 to −0.02], p=0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10−57 for both scores) but not with BMI (synthesis score, p=0.88; metabolism score, p=0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: −4.2 [95% CI −7.1 to −1.3], p=0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores).
On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.
Please see later in the article for the Editors' Summary
The prevalence of obesity has increased in the last two decades and it is presently the most common and costly nutritional problem –. In the United States, one-third of the population is affected by obesity, according to the National Health and Nutrition Examination Survey . Despite a known genetic contribution, the increase in obesity prevalence has been largely attributed to lifestyle changes, which means that it is amenable to modification through public health and other interventions .
Vitamin D deficiency is another increasingly prevalent public health concern in developed countries –, and there is evidence that vitamin D metabolism, storage, and action both influence and are influenced by adiposity. Observational studies have reported an increased risk of vitamin D deficiency in those who are obese; however, the underlying explanations and direction of causality are unclear . Active vitamin D (1,25-dihydroxyvitamin D) may influence the mobilisation of free fatty acids from the adipose tissue . In vitro experiments in rats have also shown that large doses of vitamin D2 lead to increases in energy expenditure due to uncoupling of oxidative phosphorylation in adipose tissues . However, randomized controlled trials (RCTs) testing the effect of vitamin D supplementation on weight loss in obese or overweight individuals have provided inconsistent findings –. It has also been suggested that obesity could result from an excessive adaptive winter response, and that the decline in vitamin D skin synthesis due to reduced sunlight exposure contributes to the tendency to increase fat mass during the colder periods of the year ,. However, vitamin D is stored in the adipose tissue and, hence, perhaps the most likely explanation for the association is that the larger storage capacity for vitamin D in obese individuals leads to lower circulating 25-hydroxyvitamin D [25(OH)D] concentrations, a marker for nutritional status .
In the Mendelian randomization (MR) approach, causality is inferred from associations between genetic variants that mimic the influence of a modifiable environmental exposure and the outcome of interest . If lower vitamin D intake/status is causally related to obesity, a genetic variant associated with lower 25(OH)D concentrations should be associated with higher body mass index (BMI) (in proportion to the effect on 25(OH)D). Conversely, if obesity leads to lower vitamin D status, then genetic variants associated with higher BMI should be related to lower 25(OH)D concentrations. The genetic associations, unlike the directly observed associations for vitamin D intake/status, should be less prone to confounding by lifestyle and socio-economic factors and be free from reverse causation as genotypes are invariant and assigned at random before conception . The use of multiple SNPs to index the intermediate exposure of interest increases power and reduces the risk of alternative biological pathways (pleiotropy) affecting the observed associations between the genotype and the outcome ,.
In the present study, we investigated the relationship between BMI, a commonly used measure for monitoring the prevalence of obesity at the population level, and vitamin D status and we inferred causality by using genetic variants as instruments in bi-directional MR analyses. Meta-analysis included data from 21 studies comprising up to 42,024 individuals.
All participants provided written, informed consent, and ethical permission was granted by the local research ethics committees for all participating studies.
The collaboration investigating the association of vitamin D and the risk of cardiovascular disease and related traits (D-CarDia) consists of European ancestry cohorts from the United Kingdom (UK), United States (US), Canada, Finland, Germany, and Sweden. This study comprised a meta-analysis of directly genotyped and imputed SNPs from 21 cohorts totalling 42,024 individuals (Table 1). An expanded description of the participating studies is provided in the Text S2.
To replicate our findings on the association between the vitamin D-related SNPs and allele scores with BMI, we used the data from the genome-wide meta-analyses on BMI conducted as part of the Genetic Investigation of Anthropometric Traits (GIANT) consortium . The GIANT meta-analyses consisted of 46 studies with up to 123,865 adults of European ancestry, including the 1958 British Birth Cohort, Framingham Heart study, Nurses' Health Study, Twins UK, UK Blood Services Common Control Collection, the Amish Family Osteoporosis Study, Health2000 GENMETS sub-sample, and Northern Finland Birth Cohort 1966, which were also part of the D-CarDia collaboration.
We selected 12 established BMI-related SNPs (fat mass and obesity-associated, [FTO]- rs9939609, melanocortin 4 receptor [MC4R]- rs17782313, transmembrane protein 18 [TMEM18]- rs2867125, SH2B adaptor protein 1 [SH2B1]- rs7498665, brain-derived neurotrophic factor [BDNF]- rs4074134, potassium channel tetramerisation domain containing 15 [KCTD15]- rs29941, ets variant 5 [ETV5]- rs7647305, SEC16 homolog B [SEC16B]- rs10913469, Fas apoptotic inhibitory molecule 2 [FAIM2]- rs7138803, neuronal growth regulator 1 [NEGR1]- rs3101336, mitochondrial carrier 2 [MTCH2]- rs10838738, and glucosamine-6-phosphate deaminase 2 [GNPDA2]- rs10938397) for our analysis based on the study by Li et al.  and previously published genome-wide association studies for obesity-related traits ,,. The four vitamin D-related SNPs (DHCR7- rs12785878, CYP2R1- rs10741657, GC- rs2282679, and CYP24A1- rs6013897) were chosen on the basis of the recent genome-wide association study on 25(OH)D . The studies that did not have genotyped data analysed imputed or proxy SNPs (r2=1) as available (with a call threshold of 0.9 for the SNPs imputed with Impute; for those imputed with MACH, a call threshold of 0.8 was used) . The genetic data for most studies were obtained from genome-wide association platforms, but for some studies, variants were genotyped de novo (MRC Ely, the Canadian Multicentre Osteoporosis Study, the Hertfordshire cohort study) or obtained through metabochip custom array (MRC Ely). Five studies did not have all the BMI-related SNPs (Framingham Heart Study [one missing SNP], Hertfordshire cohort study [three missing SNPs], InCHIANTI [two missing SNPs], PIVUS [two missing SNPs], and ULSAM [three missing SNPs]) and were still included in the BMI allele score analysis. Table S1 shows the minor allele frequencies for the BMI and vitamin D SNPs that were included in the analysis. A detailed description of the genotyping methods is provided in Text S2.
Analyses in each study were performed according to a standardized analysis plan. When used as outcome variables, 25(OH)D and BMI were natural log transformed to be more closely approximated by normal distributions. If multiplied by 100, coefficients from linear regression models with ln transformed outcomes can be interpreted as the percentage difference in the outcome . Models with BMI as an outcome were adjusted for age, sex, geographical site, and/or principal components from population stratification analysis (depending on data available); models with 25(OH)D as the outcome were additionally adjusted for month of blood sample collection (as a categorical variable) to account for seasonal variation and laboratory batch, where relevant. To assess the BMI relationship with 25(OH)D and vice versa, each study ran linear regression models adjusting for the covariates listed for each outcome, and the models were repeated stratifying by sex.
For the BMI SNPs, the effect allele was the BMI raising allele as established by Speliotes et al. . We created a weighted score in each study , by multiplying each SNP (coded as 0–2) by a weight based on its effect size with BMI in the meta-analysis by Speliotes et al. . The weighted BMI allele score was rescaled over the sum of weights for the available SNPs in each study to facilitate interpretation . For the vitamin D SNPs, the effect allele was the 25(OH)D lowering allele as established by the SUNLIGHT Consortium . As external weights were not available and the use of internal weights could bias the instrumental variable (IV) results , we performed an unweighted allele score analysis for the vitamin D SNPs. Vitamin D SNPs were used to form two separate allele scores : a “synthesis” allele score, created by summing the risk alleles in DHCR7 and CYP2R1, and a “metabolism” allele score, created by summing the risk alleles in GC and CYP24A1 (Figure S1). Synthesis allele score was not created for the LURIC study (one missing SNP) and both synthesis and metabolism allele scores were not created for the MRC Ely study (two missing SNPs). The synthesis allele score included the SNPs that contribute directly to the production of 25(OH)D, and hence, for which the association with the outcome can be readily estimated based on the magnitude of the association between the score and 25(OH)D . All analyses were done separately for the “metabolism” SNPs that are involved in the clearance or transport of 25(OH)D (with possible influences on bioavailability ) as the quantification of the association with the outcome based on the observed SNP-25(OH)D association is more difficult . We also evaluated the joint contribution of synthesis and metabolism scores on BMI by including both vitamin D scores as separate variables in a multiple regression model. To examine the strength of the allele scores as instruments, the F-statistic was approximated from the proportion of variation in the respective phenotype (R2) explained by the allele score, [F-stat=(R2×(n−2))/(1−R2)] .
To confirm our findings on the association between the vitamin D-related SNPs and allele scores with BMI in a larger sample, we used the summary statistics for the four vitamin D-related SNPs from the GIANT consortium. These SNPs were combined into synthesis and metabolism allele scores using an approximation method as previously described . The individual SNP association with BMI is then weighted according to its predefined effect size and meta-analysed using the inverse-variance method with the other SNPs in the score . The formal MR analyses to estimate the possible causal effect of BMI on 25(OH)D (and vice versa) were done using the IV ratio method ,. To estimate the IV ratio for the BMI effect on 25(OH)D, the meta-analysed association of the BMI allele score with 25(OH)D was divided by the association of BMI allele score with BMI. The variance for the IV ratio was estimated using a Taylor expansion . The corresponding calculation was done to establish the 25(OH)D effect on BMI, with the IV ratio method applied separately for the two vitamin D allele scores. The joint contribution of the two vitamin D scores on BMI was assessed by multivariate meta-analysis , which incorporated the covariance matrix as estimated by study specific analyses.
In the presence of heterogeneity of association between the studies, random effects meta-analyses  were run, otherwise fixed effects models were used. Univariate meta-regression models were run to assess differences in the observed associations by study level factors of sex, average BMI (BMI≤25 kg/m2 versus >25 kg/m2), the average age of participants (≤40, 41–60, and ≥61 y old), continent (North America versus Europe), and vitamin D assay (radio-immunoassay, enzyme-linked radio-immunoassay, and mass spectrometry). Power calculations for IV regression were performed by simulation  on the basis of associations observed between the phenotypes and their genetic proxies. For comparability across instruments/outcomes, power was determined for 0.02 log unit increase/decrease by decile, approximately corresponding to the association observed between BMI and 25(OH)D. To evaluate the ability to detect weaker effects on BMI using the synthesis and metabolism scores, power was also calculated for a 50% weaker effect (0.01 log unit increase/decrease). All meta-analyses and power calculations were performed at the Institute of Child Health (University College London, London) using STATA version 12 .
In the meta-analyses of 21 studies, each unit (kg/m2) increase in BMI was associated with 1.15% (95% CI 0.94%–1.36%, p=6.52×10−27) lower concentrations of 25(OH)D after adjusting for age, sex, laboratory batch, month of measurement, and principal components. The inverse association between BMI and 25(OH)D was stronger among the studies from North America than those from Europe (−1.58% [−1.81% to −1.36%], p=1.01×10−43 versus −0.91% [−1.18% to −0.64%], p=4.55×10−11; pmeta-regression=0.004) and for women than men (−1.43% [−1.65% to −1.22%], p=1.13×10−38 versus −0.75% [−1.00% to −0.50%], p=3.89×10−9; pmeta-regression=4.10×10−4) while no variation was seen by average age (pmeta-regression=0.78) or BMI (pmeta-regression=0.48) (Figure 1A and 1B).
The BMI allele score created from the 12 BMI-related SNPs showed a positive dose-response association with BMI (per unit increase 0.14% [0.12%–0.16%], p=6.30×10−62), and both vitamin D allele scores showed the expected strong associations with 25(OH)D (per allele in synthesis score: −3.47% [−3.90% to −3.05%], p=8.07×10−57; metabolism allele score: −5.38% [−5.84% to −4.93%], p=1.07×10−118) (Figures 2, S2, and S3). The BMI allele score was also associated with 25(OH)D concentrations (per unit increase −0.06%, [−0.10% to −0.02%], p=0.004) (Figure 3), while no association with BMI was seen for either the vitamin D synthesis or metabolism allele scores (per allele in synthesis score: 0.01% [−0.17% to 0.20%], p=0.88, metabolism allele score: 0.17% [−0.02% to 0.35%], p=0.08]) (Figure 4A and 4B). Analyses of joint effects by synthesis and metabolism scores provided no evidence for an association between 25(OH)D and BMI (per allele in synthesis score −0.03% [−0.23% to 0.16%] and metabolism score 0.17% [−0.04% to 0.37%], joint contribution p=0.26).
In the analyses to establish the direction and causality of BMI–25(OH)D association by the use of the IV ratio, BMI was associated with 25(OH)D: each 10% increase in BMI lead to a 4.2% decrease in 25(OH)D concentrations (−7.1% to −1.3%; p=0.005). However, the IV ratio analyses provided little evidence for a causal effect of 25(OH)D on BMI (p≥0.08 for both). We have summarised the coefficients for the MR analyses in Table 2.
The lack of association of the vitamin D allele scores with BMI was further confirmed using the GIANT consortium including 123,864 individuals in 46 studies : neither the synthesis nor the metabolism allele score showed any evidence for an association with BMI (p≥0.57 for both) (Table 3).
The BMI SNPs and the vitamin D SNPs were all individually associated with BMI and 25(OH)D, respectively (Figures S4 and S5). The exception was KCTD15 SNP, which despite previous evidence for an association , was not associated with BMI in our meta-analyses. Across the studies, the 12 BMI SNPs combined as the BMI allele score explained 0.97% of the variation in BMI (F-statistic=316; n=32,391). The synthesis allele score explained 0.64% (F-statistic=230; n=35,873) and the metabolism allele score 1.26% (F statistic=489; n=38,191) of the variation in 25(OH)D. There was no evidence for variation in the BMI allele score–BMI association by continent (p meta-regression=0.15) or BMI (p meta-regression=0.83). However, the BMI allele score–BMI association was slightly weaker in studies with older compared to younger participants (−0.03% [−0.05% to −0.002%], p meta-regression=0.03). The vitamin D allele score–25(OH)D association did not vary by age, BMI, continent, or assay (p meta-regression≥0.09 for all comparisons).
Of the 12 individual BMI SNPs, the SNP for FTO was the only one that showed evidence of a univariate association with 25(OH)D (p=0.050) (Figure S6). None of the four 25(OH)D SNPs were individually associated with BMI (p≥0.10) (Figure S7). The lack of association of the four vitamin D SNPs with BMI was further confirmed using the summary data from the GIANT consortium (p>0.30 for all the SNPs) (Table 3).
The association between BMI allele score and 25(OH)D did not vary by study level factors, including age (pmeta-regression=0.40), BMI (pmeta-regression=0.18), continent of study (pmeta-regression=0.78), or vitamin D assay (pmeta-regression=0.23). Similarly, there was no evidence for variation in the vitamin D allele score–BMI association by age (pmeta-regression≥0.25 for both scores), or continent (pmeta-regression≥0.50 for both scores). There was also no strong evidence for variation in the vitamin D allele score–BMI association by average BMI of the study (≤25 kg/m2 versus ≥25 kg/m2), although for the synthesis score the meta-regression coefficient was of borderline significance (pmeta-regression=0.053, Figure S8; pmeta-regression=0.78 for metabolism score).
Illustrative power calculations are presented in Figure S9. In theory, we had greater power to detect an association between 25(OH)D and BMI using the metabolism score as an instrument, compared with an equal sized association between BMI and 25(OH)D using the BMI risk score. However, if the size of the association between 25(OH)D and BMI was only half that seen between BMI and 25(OH)D, our study would not have been adequately powered even with the inclusion of the GIANT results.
Obesity, and perhaps vitamin D deficiency, are among the most important modifiable risk factors for a number of chronic diseases. Obesity and vitamin D status are known to be associated but the direction of the association and whether it is causal has been uncertain. We have presented genetic evidence that higher BMI leads to lower vitamin D status. Conversely, our analyses provided no evidence for a causal role of vitamin D in the development of obesity, although our study was not powered to detect very small effects. These results suggest that although increases in vitamin D status are not likely to help with weight regulation, increased risk of vitamin D deficiency could contribute to the adverse health effects associated with obesity.
The association between obesity and vitamin D status was remarkably consistent across the different populations included in our meta-analyses, being apparent both in men and in women, and in the young and older cohorts alike. Interestingly, the association between obesity and 25(OH)D concentrations appeared stronger for populations in North America compared to Europe, possibly reflecting differences in the distribution of BMI across the continents. Recent intervention studies have shown that obese individuals need higher vitamin D dosages than lean individuals to achieve the same 25(OH)D concentrations ,. Given that North America has one of the highest rates of obesity in the world , our study highlights the importance of considering obesity as a risk factor for vitamin D deficiency with implications on the dosage requirements and possible targeting of relevant health promotion strategies.
The lack of any suggestion for an association between the vitamin D SNPs and BMI in the GIANT consortium (n=123,864) alongside our own large meta-analyses provides a strong case against linear increases in 25(OH)D having a substantive influence on BMI. This conclusion is in accordance with a recent study on Chinese women (n=7,000), which also failed to observe evidence for an association with BMI for genetic variants in the vitamin D pathway . Although a recent RCT (n=77) suggested greater loss in fat mass for women receiving vitamin D , previous trials have failed to show any evidence for an effect despite larger treatment groups (n=200–445), use of higher vitamin D dosages, and equal duration of treatment (12 mo) ,. Dilution related to the greater volume of distribution has been recently proposed as the most likely explanation for the lower 25(OH)D concentrations in obese individuals . In that study, no evidence was found for reduced bioavailability through increased sequestration of vitamin D in the adipose tissue, which had previously been suggested to contribute to the low 25(OH)D concentrations in obesity . In contrast, intact parathyroid hormone (iPTH) levels , which stimulate the 1-α-hydroxylase (CYP27B1) enzyme that converts 25(OH)D to 1,25-dihydroxyvitamin D (the active hormonal form), have been found to be elevated in obesity , which could to some extent also contribute to the lower 25(OH)D concentrations in obese individuals. It is also possible that differences in lifestyle could contribute to lower 25(OH)D concentrations in obese compared to normal weight individuals, although the association between obesity and low 25(OH)D concentrations has been found to only modestly attenuate after adjustment for vitamin D-related lifestyle and dietary factors .
The main strengths of this study are the large sample size and the individual level population-based data from North America and Europe. We used a bi-directional MR approach to investigate the causal directions between obesity and vitamin D deficiency, observing evidence for reductions in 25(OH)D by BMI but not vice versa. However, based on the biological pathways proposed, a possible effect of 25(OH)D on BMI could be expected to be weaker than the effect of BMI on 25(OH)D. Despite including data from the large GIANT consortium to narrow the range of effects compatible with the data, we are unable to exclude very small effects. Furthermore, while the MR approach enables the approximation of life-long differences in average concentrations, with genetic markers it is not possible to examine the influences arising from the extremes of non-linear distributions . Consequently, we cannot discount a possible effect of severe vitamin D deficiency on BMI due to evidence of non-linearity seen in some studies . In contrast, associations between BMI and 25(OH)D within levels in the obesity range were consistently linear in studies included in our analyses (unpublished data), hence the observed association between higher BMI and lower 25(OH)D is likely to be informative in the context of obesity.
One of the methodological challenges of the MR approach relates to the large sample size requirement, arising from the availability of relatively weak instruments for most exposures ,. This aspect of the MR approach is also reflected in our study, notably in the relatively small amount of variation explained by all the instruments used. We used the IV ratio method on meta-analyzed coefficients since all studies were not able to share individual level participant data. This method assumes linear relationships and may have less power to detect an effect than other IV methods . However, as shown by the clear outcome of these analyses, we were able to overcome these issues by combining several cohorts with comparable information, allowing us to achieve the large numbers required (maximum n=42,024) . To confirm the lack of association between vitamin D-related genetic variations and BMI, we were able to expand the analyses by using data from the large GIANT meta-analyses (n=123,864) . However, this cannot be considered an independent replication, as eight of the studies that were part of the D-CarDia Collaboration were also included in GIANT. The F-statistic is used to measure the strength of an instrument, and an instrument that has a value greater than 10 is considered strong enough to use in IV analyses . In our analyses, the F-statistic was greater than 200 for all instruments used due to our large sample size.
Combining large population-based studies from North America and Europe could lead to confounding by population stratification; however, we adjusted for geographical variation/principal components in all analyses, which appeared adequate, as there was no evidence for heterogeneity by continent for the allele score meta-analyses. An important benefit of the MR approach is that it helps to overcome problems of confounding and reverse causality, which limit the ability to draw causal inferences in non-genetic observational studies ,. However, it could be argued that as the biological function for some of the BMI SNPs is yet to be established, there could be alternative biological pathways explaining their association with BMI. Using multiple SNPs to index BMI, we were able to minimise the risk of pleiotropic effects, as the effects of alternative pathways reflected by individual SNPs would be expected to be strongly diluted when combined in a multi marker score ,.
In conclusion, we demonstrated that the association between BMI and lower 25(OH)D concentrations in Caucasian populations from North America and Europe can be seen across different age groups and in both men and women. We also show that higher BMI leads to lower vitamin D status, providing evidence for the role of obesity as a causal risk factor for the development of vitamin D deficiency. Together with the suggested increases in vitamin D requirements in obese individuals ,, our study highlights the importance of monitoring and treating vitamin D deficiency as a means of alleviating the adverse influences of excess adiposity on health. Our findings suggest that population level interventions to reduce obesity would be expected to lead to a reduction in the prevalence of vitamin D deficiency.
Vitamin D pathway showing the “synthesis” and “metabolism” indicators.
Meta-analysis of BMI allele score association with BMI in collaborating studies ( n =32,391).
Meta-analysis of synthesis (A) ( n =35,873) and metabolism (B) ( n =38,191) allele score associations with 25(OH)D in collaborating studies.
Association of the 12 BMI-related SNPs with BMI.
Association of the four vitamin D SNPs with 25(OH)D.
Association of the 12 BMI-related SNPs with 25(OH)D.
Association of the four vitamin D SNPs with BMI.
Meta-analysis of synthesis allele score association with BMI stratified by mean BMI of collaborating studies ( n =36,553).
Power to detect an association between BMI and 25(OH)D using genetic proxies in instrumental variable regression. The solid lines represent the power to detect the same effect size for a BMI association with 25(OH)D using BMI risk score (black line), 25(OH)D association with BMI using metabolism score (mid grey line), and 25(OH)D association with BMI using synthesis score (dark grey line). The dash lines represent the power to detect an effect size half that of the same coloured solid lines for the 25(OH)D association with BMI.
Minor allele frequency (MAF) for the BMI and vitamin D-related SNPs.
STROBE statement—checklist of items.
Study details of the 21 participating cohorts and the funding information.
We thank Alana Cavadino (UCL ICH) for technical assistance in formatting the figures.
Members of the Genetic Investigation of Anthropometric Traits (GIANT) Consortium include:
Elizabeth K Speliotes, Cristen J Willer, Sonja I Berndt, Keri L Monda, Gudmar Thorleifsson, Anne U Jackson, Hana Lango Allen, Cecilia M Lindgren, Jian'an Luan, Reedik Mägi, Joshua C Randall, Sailaja Vedantam, Thomas W Winkler, Lu Qi, Tsegaselassie Workalemahu, Iris M Heid, Valgerdur Steinthorsdottir, Heather M Stringham, Michael N Weedon, Eleanor Wheeler, Andrew R Wood, Teresa Ferreira, Robert J Weyant, Ayellet V Segrè, Karol Estrada, Liming Liang, James Nemesh, Ju-Hyun Park, Stefan Gustafsson, Tuomas O Kilpeläinen, Jian Yang, Nabila Bouatia-Naji, Tõnu Esko, Mary F Feitosa, Zoltán Kutalik, Massimo Mangino, Soumya Raychaudhuri, Andre Scherag, Albert Vernon Smith, Ryan Welch, Jing Hua Zhao, Katja K Aben, Devin M Absher, Najaf Amin, Anna L Dixon, Eva Fisher, Nicole L Glazer, Michael E Goddard, Nancy L Heard-Costa, Volker Hoesel, Jouke-Jan Hottenga, Åsa Johansson, Toby Johnson, Shamika Ketkar, Claudia Lamina, Shengxu Li, Miriam F Moffatt, Richard H Myers, Narisu Narisu, John R B Perry, Marjolein J Peters, Michael Preuss, Samuli Ripatti, Fernando Rivadeneira, Camilla Sandholt, Laura J Scott, Nicholas J Timpson, Jonathan P Tyrer, Sophie van Wingerden, Richard M Watanabe, Charles C White, Fredrik Wiklund, Christina Barlassina, Daniel I Chasman, Matthew N Cooper, John-Olov Jansson, Robert W Lawrence, Niina Pellikka, Inga Prokopenko, Jianxin Shi, Elisabeth Thiering, Helene Alavere, Maria T S Alibrandi, Peter Almgren, Alice M Arnold, Thor Aspelund, Larry D Atwood, Beverley Balkau, Anthony J Balmforth, Amanda J Bennett, Yoav Ben-Shlomo, Richard N Bergman, Sven Bergmann, Heike Biebermann, Alexandra I F Blakemore, Tanja Boes, Lori L Bonnycastle, Stefan R Bornstein, Morris J Brown, Thomas A Buchanan, Fabio Busonero, Harry Campbell, Francesco P Cappuccio, Christine Cavalcanti-Proença, Yii-Der Ida Chen, Chih-Mei Chen, Peter S Chines, Robert Clarke, Lachlan Coin, John Connell, Ian N M Day, Martin den Heijer, Jubao Duan, Shah Ebrahim, Paul Elliott, Roberto Elosua, Gudny Eiriksdottir, Michael R Erdos, Johan G Eriksson, Maurizio F Facheris, Stephan B Felix, Pamela Fischer-Posovszky, Aaron R Folsom, Nele Friedrich, Nelson B Freimer, Mao Fu, Stefan Gaget, Pablo V Gejman, Eco J C Geus, Christian Gieger, Anette P Gjesing, Anuj Goel, Philippe Goyette, Harald Grallert, Jürgen Gräßler, Danielle M Greenawalt, Christopher J Groves, Vilmundur Gudnason, Candace Guiducci, Anna-Liisa Hartikainen, Neelam Hassanali, Alistair S Hall, Aki S Havulinna, Caroline Hayward, Andrew C Heath, Christian Hengstenberg, Andrew A Hicks, Anke Hinney, Albert Hofman, Georg Homuth, Jennie Hui, Wilmar Igl, Carlos Iribarren, Bo Isomaa, Kevin B Jacobs, Ivonne Jarick, Elizabeth Jewell, Ulrich John, Torben Jørgensen, Pekka Jousilahti, Antti Jula, Marika Kaakinen, Eero Kajantie, Lee M Kaplan, Sekar Kathiresan, Johannes Kettunen, Leena Kinnunen, Joshua W Knowles, Ivana Kolcic, Inke R König, Seppo Koskinen, Peter Kovacs, Johanna Kuusisto, Peter Kraft, Kirsti Kvaløy, Jaana Laitinen, Olivier Lantieri, Chiara Lanzani, Lenore J Launer, Cecile Lecoeur, Terho Lehtimäki, Guillaume Lettre, Jianjun Liu, Marja-Liisa Lokki, Mattias Lorentzon, Robert N Luben, Barbara Ludwig, Paolo Manunta, Diana Marek, Michel Marre, Nicholas G Martin, Wendy L McArdle, Anne McCarthy, Barbara McKnight, Thomas Meitinger, Olle Melander, David Meyre, Kristian Midthjell, Grant W Montgomery, Mario A Morken, Andrew P Morris, Rosanda Mulic, Julius S Ngwa, Mari Nelis, Matt J Neville, Dale R Nyholt, Christopher J O'Donnell, Stephen O'Rahilly, Ken K Ong, Ben Oostra, Guillaume Paré, Alex N Parker, Markus Perola, Irene Pichler, Kirsi H Pietiläinen, Carl G P Platou, Ozren Polasek, Anneli Pouta, Suzanne Rafelt, Olli Raitakari, Nigel W Rayner, Martin Ridderstråle, Winfried Rief, Aimo Ruokonen, Neil R Robertson, Peter Rzehak, Veikko Salomaa, Alan R Sanders, Manjinder S Sandhu, Serena Sanna, Jouko Saramies, Markku J Savolainen, Susann Scherag, Sabine Schipf, Stefan Schreiber, Heribert Schunkert, Kaisa Silander, Juha Sinisalo, David S Siscovick, Jan H Smit, Nicole Soranzo, Ulla Sovio, Jonathan Stephens, Ida Surakka, Amy J Swift, Mari-Liis Tammesoo, Jean-Claude Tardif, Maris Teder-Laving, Tanya M Teslovich, John R Thompson, Brian Thomson, Anke Tönjes, Tiinamaija Tuomi, Joyce B J van Meurs, Gert-Jan van Ommen, Vincent Vatin, Jorma Viikari, Sophie Visvikis-Siest, Veronique Vitart, Carla I G Vogel, Benjamin F Voight, Lindsay L Waite, Henri Wallaschofski, G Bragi Walters, Elisabeth Widen, Susanna Wiegand, Sarah H Wild, Gonneke Willemsen, Daniel R Witte, Jacqueline C Witteman, Jianfeng Xu, Qunyuan Zhang, Lina Zgaga, Andreas Ziegler, Paavo Zitting, John P Beilby, I Sadaf Farooqi, Johannes Hebebrand, Heikki V Huikuri,Alan L James, Mika Kähönen, Douglas F Levinson, Fabio Macciardi, Markku S Nieminen, Claes Ohlsson, Lyle J Palmer, Paul M Ridker, Michael Stumvoll, Jacques S Beckmann, Heiner Boeing, Eric Boerwinkle, Dorret I Boomsma, Mark J Caulfield, Stephen J Chanock, Francis S Collins, L Adrienne Cupples, George Davey Smith, Jeanette Erdmann, Philippe Froguel, Henrik Grönberg, Ulf Gyllensten, Per Hall, Torben Hansen, Tamara B Harris, Andrew T Hattersley, Richard B Hayes, Joachim Heinrich, Frank B Hu, Kristian Hveem, Thomas Illig, Marjo-Riitta Jarvelin, Jaakko Kaprio, Fredrik Karpe, Kay-Tee Khaw, Lambertus A Kiemeney, Heiko Krude, Markku Laakso, Debbie A Lawlor, Andres Metspalu, Patricia B Munroe, Willem H Ouwehand, Oluf Pedersen, Brenda W Penninx, Annette Peters, Peter P Pramstaller, Thomas Quertermous, Thomas Reinehr, Aila Rissanen, Igor Rudan, Nilesh J Samani, Peter E H Schwarz, Alan R Shuldiner, Timothy D Spector, Jaakko Tuomilehto, Manuela Uda, André Uitterlinden, Timo T Valle, Martin Wabitsch, Gérard Waeber, Nicholas J Wareham,Hugh Watkins, James F Wilson, Alan F Wright, M Carola Zillikens, Nilanjan Chatterjee, Steven A McCarroll, Shaun Purcell, Eric E Schadt, Peter M Visscher, Themistocles L Assimes, Ingrid B Borecki, Panos Deloukas, Caroline S Fox, Leif C Groop, Talin Haritunians, David J Hunter, Robert C Kaplan, Karen L Mohlke, Jeffrey R O'Connell, Leena Peltonen, David Schlessinger, David P Strachan, Cornelia M van Duijn, H-Erich Wichmann, Timothy M Frayling, Unnur Thorsteinsdottir, Gonçalo R Abecasis, Inês Barroso, Michael Boehnke, Kari Stefansson, Kari E North, Mark I McCarthy, Joel N Hirschhorn, Erik Ingelsson and Ruth J F Loos.
The authors thank the British Heart Foundation (grant PG/09/023) and the UK Medical Research Council (MRC; grant G0601653) for funding this work. ADH is a British Heart Foundation Senior Research Fellow (Award FS05/125). EH is a Department of Health (UK) Public Health Career Scientist. This work was undertaken at the Centre for Paediatric Epidemiology and Biostatistics, which benefits from funding support from the MRC in its capacity as the MRC Centre of Epidemiology for Child Health. Research at the University College London Institute of Child Health and Great Ormond Street Hospital for Children NHS Trust benefits from R&D funding received from the NHS Executive. No funding bodies had any role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.