Levin GP, Robinson-Cohen C, de Boer IH, Houston DK, Lohman K, Liu Y, et al. Genetic variants and associations of 25-hydroxyvitamin D concentrations with major clinical outcomes. JAMA. 2012;308(18):doi:10.1001/jama.2012.17304.
eTable 1. Detailed methods for measurement of 25(OH)D concentrations and genotyping in the Cardiovascular Health Study, Health ABC, InCHIANTI, and ULSAM cohorts
eTable 2. Incident events in the Cardiovascular Health Study, Health ABC, InCHIANTI, and ULSAM cohorts
eTable 3. Detailed results in the Cardiovascular Health Study, Health ABC, InCHIANTI, and ULSAM cohorts for interactions between low vitamin D and VDR SNPs rs7968585 and rs2239179 on risk of the composite outcome
eFigure 1. Flow chart defining the Cardiovascular Health Study discovery cohort
eFigure 2. Penalized spline describing the association between season-specific 25-hydroxyvitamin D Z score and risk of the composite outcome, adjusted for age and sex, in the Cardiovascular Health Study
eFigure 3. Estimates of interaction between low vitamin D and the number of minor alleles on risk of the composite outcome, by study and overall, for the two VDR variants: rs7968585 and rs2239179
eFigure 4. Kaplan-Meier plots describing the association of low vitamin D with risk of the composite outcome in the Cardiovascular Health Study, stratified according to the number of copies of the minor allele for the 2 VDR variants: rs7968585 and rs2239179
Gregory P. Levin, Cassianne Robinson-Cohen, Ian H. de Boer, Denise K. Houston, Kurt Lohman, Yongmei Liu, Stephen B. Kritchevsky, Jane A. Cauley, Toshiko Tanaka, Luigi Ferrucci, Stefania Bandinelli, Kushang V. Patel, Emil Hagström, Karl Michaëlsson, Håkan Melhus, Thomas Wang, Myles Wolf, Bruce M. Psaty, David Siscovick, Bryan Kestenbaum. Genetic Variants and Associations of 25-Hydroxyvitamin D Concentrations With Major Clinical Outcomes. JAMA. 2012;308(18):1898–1905. doi:10.1001/jama.2012.17304
Author Affiliations: Departments of Biostatistics (Dr Levin), Epidemiology (Drs Robinson-Cohen, de Boer, Psaty, Siscovick, and Kestenbaum), Medicine (Drs de Boer, Psaty, Siscovick, and Kestenbaum), and Anasthesiology and Pain Medicine (Dr Patel), Division of Nephrology and Kidney Research Institute (Drs de Boer and Kestenbaum), Center for Pain Research on Impact, Measurement, and Effectiveness (Dr Patel), Cardiovascular Health Research Unit (Drs Psaty and Siscovick), Group Health Research Institute and Cooperative (Dr Psaty), University of Washington, Seattle; Departments of Internal Medicine (Drs Houston and Kritchevsky), Biostatistical Sciences (Mr Lohman), and Epidemiology and Prevention (Dr Liu), Division of Public Health Sciences (Mr Lohman and Dr Liu), School of Medicine (Drs Houston and Kritchevsky and Mr Lohman), Wake Forest University, Winston-Salem, North Carolina; Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania (Dr Cauley); Clinical Research Branch, National Institute on Aging, Baltimore, Maryland (Drs Tanaka and Ferrucci); Geriatric Unit, Azienda Sanitaria Firenze, Florence, Italy (Dr Bandinelli); Departments of Medical Sciences (Drs Hagström and Melhus) and Surgical Sciences (Dr Michaëlsson), Uppsala University, Uppsala, Sweden; Division of Cardiology, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts (Dr Wang); and Miller School of Medicine, University of Miami, Miami, Florida (Dr Wolf).
Context Lower serum 25-hydroxyvitamin D concentrations are associated with greater risks of many chronic diseases across large, prospective community-based studies. Substrate 25-hydroxyvitamin D must be converted to 1,25-dihydroxyvitamin D for full biological activity, and complex metabolic pathways suggest that interindividual variability in vitamin D metabolism may alter the clinical consequences of measured serum 25-hydroxyvitamin D.
Objective To investigate whether common variation within genes encoding the vitamin D–binding protein, megalin, cubilin, CYP27B1, CYP24A1, and the vitamin D receptor (VDR) modify associations of low 25-hydroxyvitamin D with major clinical outcomes.
Design, Setting, and Participants Examination of 141 single-nucleotide polymorphisms in a discovery cohort of 1514 white participants (who were recruited from 4 US regions) from the community-based Cardiovascular Health Study. Participants had serum 25-hydroxyvitamin D measurements in 1992-1993 and were followed up for a median of 11 years (through 2006). Replication meta-analyses were conducted across the independent, community-based US Health, Aging, and Body Composition (n = 922; follow-up: 1998-1999 through 2005), Italian Invecchiare in Chianti (n = 835; follow-up: 1998-2000 through 2006), and Swedish Uppsala Longitudinal Study of Adult Men (n = 970; follow-up: 1991-1995 through 2008) cohort studies.
Main Outcome Measure Composite outcome of incident hip facture, myocardial infarction, cancer, and mortality over long-term follow-up.
Results Interactions between 5 single-nucleotide polymorphisms and low 25-hydroxyvitamin D concentration were identified in the discovery phase and 1 involving a variant in the VDR gene replicated in independent meta-analysis. Among Cardiovascular Health Study participants, low 25-hydroxyvitamin D concentration was associated with hazard ratios for risk of the composite outcome of 1.40 (95% CI, 1.12-1.74) for those who had 1 minor allele at rs7968585 and 1.82 (95% CI, 1.31-2.54) for those with 2 minor alleles at rs7968585. In contrast, there was no evidence of an association (estimated hazard ratio, 0.93 [95% CI, 0.70-1.24]) among participants who had 0 minor alleles at this single-nucleotide polymorphism.
Conclusion Known associations of low 25-hydroxyvitamin D with major health outcomes may vary according to common genetic differences in the vitamin D receptor.
Vitamin D status is defined by the circulating concentration of 25-hydroxyvitamin D.1,2 In large, prospective cohort studies, lower serum 25-hydroxyvitamin D concentrations are associated with greater risks of hip fracture, myocardial infarction (MI), cancer, and death.3- 10 In experimental models, disruption of the vitamin D–endocrine axis stimulates inflammatory cytokines, activates the renin-angiotensin system, and impairs skeletal mineralization.11- 13 The totality of these findings suggests that low 25-hydroxyvitamin D concentration may be a modifiable risk factor for many chronic diseases, which has motivated ongoing clinical trials to test whether vitamin D supplementation can reduce the risk of disease development.14 Laboratory testing for serum 25-hydroxyvitamin D and empirical therapy with vitamin D supplements have increased dramatically worldwide, with substantial associated costs.
Substrate 25-hydroxyvitamin D must be converted to 1,25-dihydroxyvitamin D, which is the potent hormonal form of vitamin D, for full biological activity. Conversion requires transportation in the blood by the vitamin D–binding protein, internalization via cell surface proteins megalin and cubilin, and metabolism by the 1-α hydroxylase enzyme. Activated 1,25-dihydroxyvitamin D then binds to the vitamin D receptor and regulates expression of a diverse array of vitamin D responsive genes. Elimination of 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D is primarily catalyzed by the 24-α hydroxylase enzyme.1,2 These complex metabolic pathways suggest that interindividual variability in vitamin D metabolism may alter the clinical consequences of measured serum 25-hydroxyvitamin D.
We hypothesized that known serum 25-hydroxyvitamin D disease relationships would differ according to common variation in 25-hydroxyvitamin D metabolism genes. We investigated whether single-nucleotide polymorphisms (SNPs) within genes encoding proteins that reside downstream from 25-hydroxyvitamin D modified associations of low serum 25-hydroxyvitamin D concentration with a composite outcome of incident hip fracture, MI, cancer, and mortality over long-term follow-up.
The Cardiovascular Health Study (CHS) is a cohort study of cardiovascular disease among 5888 ambulatory adults aged 65 years or older living in 1 of 4 US communities.15 We measured serum 25-hydroxyvitamin D concentrations during 1992-1993 in 2312 CHS participants without prevalent cardiovascular disease.16 The Health, Aging, and Body Composition (Health ABC) study is a cohort study of changes in body composition among 3075 community-dwelling adults aged 70 to 79 years.17 Serum 25-hydroxyvitamin D concentrations were measured in 2998 participants during 1998-1999. The Invecchiare in Chianti (InCHIANTI) study is a population-based cohort study of 1453 primarily older persons living in the Chianti region of Tuscany, Italy, with 25-hydroxyvitamin D concentration ascertainment during 1998-2000.18 The Uppsala Longitudinal Study of Adult Men (ULSAM) is a cohort study of 2322 initially 50-year-old Swedish men aimed at identifying metabolic risk factors for cardiovascular disease.19 Concentrations of 25-hydroxyvitamin D were measured during 1991-1995 in 1221 men. We excluded participants from these 4 study cohorts with self-reported nonwhite race (to reduce potential confounding effects of population stratification), as well as those who had prevalent disease (hip fracture, MI, or cardiovascular disease, depending on the cohort, or cancer), unsuccessful genotyping, or failed 25-hydroxyvitamin D ascertainment. All participants provided informed consent, and institutional review boards reviewed and approved the procedures at all sites.
Circulating serum 25-hydroxyvitamin D concentrations were measured using mass spectrometry (CHS and ULSAM) or radioimmunoassay (DiaSorin RIA20; Health ABC and InCHIANTI). All interassay coefficients of variation were less than 10.2% (eTable 1). We evaluated 25-hydroxyvitamin D as a dichotomous variable because we and others3,4,6,8,10,21 have observed threshold associations of 25-hydroxyvitamin D with disease risks. We used season-specific cut points because of the known seasonal variability in 25-hydroxyvitamin D concentrations,22- 25 and its associated impact on modeling.26 We defined low vitamin D concentration as the lowest season-specific quintile to ensure a similar definition across the study cohorts, which used different methods to measure serum 25-hydroxyvitamin D.
Genotyping was performed using Illumina platforms and software (eTable 1). Genome-wide association data for the CHS was drawn from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium,27 and SNPs were excluded for a call rate of less than 97% or a Hardy-Weinberg equilibrium P value of less than 10−5. We identified the following 6 genes based on their known role in metabolic pathways downstream from 25-hydroxyvitamin D: vitamin D–binding protein (GC; NCBI Entrez Gene NG_012837.2), megalin (LRP2; NCBI Entrez Gene NG_012634.1), cubilin (CUBN; NCBI Entrez Gene NG_008967.1), 1-α hydroxylase (CYP27B1; NCBI Entrez Gene NG_007076.1), 24-α hydroxylase (CYP24A1; NCBI Entrez Gene NG_008334.1), and the vitamin D receptor (VDR; NCBI Entrez Gene NG_008731.1). We studied 141 SNPs located within 20 kb of these genes that were successfully genotyped and had a minor allele frequency in the CHS of at least 5%.
The study outcome was the time from 25-hydroxyvitamin D measurement to the first occurrence of incident hip fracture, incident MI, incident cancer (excluding non–melanoma skin cancer), or death from any cause. This composite outcome was selected prior to the analyses to capture previously described associations of 25-hydroxyvitamin D with clinically important disease outcomes and to maximize statistical power to detect potential interactions. In the CHS (follow-up until 2005 or 2006), an events committee adjudicated MI cases.28 Hip fracture was defined by the International Classification of Diseases, Ninth Revision codes, and incident cancer cases were identified by linking the CHS records with population-based cancer registries. InCHIANTI investigators ascertained outcome information by examining and interviewing participants every 3 years (through 2006), and mortality data was obtained by the Tuscany Regional Health Authority. In the Health ABC study (follow-up until December 31, 2005), participants or their proxies reported any hospitalizations, outpatient cancer, fracture, or angioplasty events. Incident hip fracture, MI, and cancer events were adjudicated using medical records and other supporting documents. Incident events in the ULSAM (through December 31, 2008) were identified using the Swedish National Patient Register and the Swedish Cause of Death Register.
All survival analyses used Cox proportional hazards models with robust standard errors, adjusting for age and sex. Additional variables such as physical activity, hypertension, diabetes, and body mass index were not included in the models because the candidate genes (selected specifically on the basis of their role in the downstream metabolism of 25-hydroxyvitamin D) are unlikely to alter the impact of these potential confounding factors on the composite disease outcome through mechanisms outside of the vitamin D pathway. We graphically displayed a penalized spline29 to describe the functional form of the association between season-specific 25-hydroxyvitamin D concentration and incident events in the CHS.
We estimated the multiplicative interaction between low vitamin D concentration and each candidate SNP on risk of the composite outcome, and computed interaction P values using Wald tests. In these models, an exponential coefficient provides an estimate of the relative difference in the association of low 25-hydroxyvitamin D with risk of the composite outcome corresponding to the presence of 1 additional copy of the minor allele (ie, an estimated hazard ratio ratio [HRR]). To account for multiple statistical testing in the discovery phase, false discovery rate q values were calculated.30 We used a false discovery rate threshold of .25, meaning that we expected up to 25% of our declared discoveries to be false. It is common to use a threshold such as .20 or .25 in candidate gene studies, which use prior knowledge to select candidate genes and which conduct replication in independent cohorts.31,32
Genetic variants with interaction q values below the .25 threshold in the CHS discovery cohort were tested in the Health ABC, InCHIANTI, and ULSAM cohorts for replication. We calculated summary HRRs and 95% confidence intervals using Mantel-Haenszel methods and a fixed-effects meta-analysis model. Given that estimates from the discovery phase are likely biased away from the null,33 we used the pooled results from the independent cohorts alone to determine successful replication at the .05 significance level. We performed tests for cohort-level heterogeneity of the interaction estimates using the Woolf method.34 We carried out additional analyses in the CHS cohort to explore the nature of any replicated interactions. Within strata defined by each genetic variant, we computed Kaplan-Meier estimates and unadjusted incidence rates of the composite outcome according to 25-hydroxyvitamin D status, as well as adjusted HRs comparing participants with normal and low 25-hydroxyvitamin D concentrations. We jointly modeled interactions between the 2 VDR SNPs and 25-hydroxyvitamin D concentration to examine whether both variants independently modified the association of low 25-hydroxyvitamin D with risk of the composite outcome.
We report 95% confidence intervals and 2-sided P values. No adjustments for multiple comparisons in the replication phase were conducted. Statistical analyses were completed using R version 2.11.0 (R Project for Statistical Computing), Stata versions 10.1 and 11.2 (StataCorp), and SAS version 9.1 (SAS Institute Inc).
In the CHS cohort, we excluded 341 participants who had prevalent hip fracture or cancer, 309 participants of nonwhite race, and 148 participants who were unable to be genotyped, resulting in a discovery study population of 1514 (eFigure 1). We excluded Health ABC participants who had prevalent cardiovascular disease (n = 771), cancer (n = 414), or hip fracture (n = 24), missing (n = 102) or high 25-hydroxyvitamin D concentrations (>100 ng/mL; n = 2), nonwhite race (n = 690), or unsuccessful genotyping (n = 73), creating a cohort of 922 participants. In the InCHIANTI study, we excluded participants who did not return for future visits or have outcomes adjudicated (n = 277), and those with missing 25-hydroxyvitamin D concentrations (n = 66), unsuccessful genotyping (n = 114), or prevalent hip fracture, MI, or cancer (n = 161), leaving 835 participants. We excluded ULSAM participants with missing kidney function measurement (n = 65), and those with a history of MI (n = 94), hip fracture (n = 14), or cancer (n = 78), creating a replication cohort of 970 men.
The CHS discovery cohort was 70% female, with a mean age of 74 years and mean serum 25-hydroxyvitamin D concentration of 26.7 ng/mL. Participants who had low 25-hydroxyvitamin D concentrations were more likely to be female, had a higher mean body mass index, and were more likely to have diabetes or hypertension vs those who had normal 25-hydroxyvitamin D concentrations. Demographic and comorbidity profiles of participants from the Health ABC, InCHIANTI, and ULSAM replication cohorts were generally similar to those of the CHS discovery cohort, with a few notable exceptions. The Health ABC population had a considerably higher prevalence of male participants and vitamin D supplement use, the InCHIANTI participants were younger and more likely to be male, and the ULSAM cohort was composed entirely of men. Mean 25-hydroxyvitamin D concentration was lower among the InCHIANTI participants, and season-specific thresholds defining low vitamin D concentration varied by cohort (Table 1).
In the CHS cohort, 948 participants (63%) experienced an event contained within the composite outcome over a median (maximum) follow-up period of 11 (15) years. Consistent with previous studies, there was evidence of a threshold association between serum 25-hydroxyvitamin D concentration and risk of the composite outcome (eFigure 2). With adjustment for age and sex, low 25-hydroxyvitamin D concentration was associated with an HR of 1.32 (95% CI, 1.13-1.54; P = .001).
Of the 141 SNPs tested for discovery in the CHS cohort, 6 met the false-discovery rate threshold of 0.25. Two are located (±20 kb) within VDR, 2 in CUBN, and 2 in CYP27B1. The 2 identified CYP27B1 variants are in near-perfect linkage disequilibrium35; we carried only rs703842 forward to the replication phase because direct information on rs2069502 was not ascertained in the ULSAM. Each additional copy of the minor allele in the 5 remaining identified SNPs was associated with a 30% to 40% difference in the HR (P value range: .002-.008). The magnitudes and directions of each interaction were consistent across the individual events that comprised the composite outcome (Table 2). An additional 5 SNPs (all in VDR or CUBN) had P values of less than .05, but q values of greater than .25. This included the well-known VDR variant BsmI (rs1544410; HRR, 0.76, P = .02).
In the Health ABC, ULSAM, and InCHIANTI cohorts, a total of 1216 participants experienced the composite outcome. Estimates of association between low 25-hydroxyvitamin D concentration and risk of the composite outcome were in the same direction in all 4 cohorts, but statistical significance was attained only in the Health ABC cohort (eTable 2). Of the 5 variants identified in the CHS cohort, only the VDR SNP rs7968585 reached statistical significance for modifying the association between low 25-hydroxyvitamin D concentration and risk of the composite outcome in independent replication meta-analyses. An additional VDR SNP (rs2239179) attained statistical significance in a meta-analysis that included CHS results (Table 3). Each additional copy of the minor allele at rs7968585 was associated with an HRR of 1.22 (95% CI, 1.09-1.36) and each additional copy of the minor allele at rs2239179 was associated with an HRR of 0.85 (95% CI, 0.76-0.97). The HRRs for these 2 interactions were reasonably consistent across the independent cohorts (eFigure 3), and there was no statistical evidence of heterogeneity. The meta-analysis interactions also demonstrated similar magnitudes and directions across the individual events that comprised the composite outcome (eTable 3). These 2 VDR SNPs are in moderate linkage disequilibrium (R2 = 0.55).35
The nature of the observed interactions can be seen by comparing estimates of disease-free survival, as well as incidence rates of the composite outcome, across groups defined by each variant and by 25-hydroxyvitamin D concentration status (CHS data; Table 4 and eFigure 4). Without accounting for 25-hydroxyvitamin D concentration, there was no evidence in the CHS cohort of main effect associations between either rs7968585 (HR per additional minor allele, 1.00 [95% CI, 0.91-1.06]; P = .93) or rs2239179 (HR, 0.93 [95% CI, 0.85-1.02]; P = .13) and risk of the composite outcome. In an analysis jointly modeling both SNPs and their interactions with 25-hydroxyvitamin D concentration, there was insufficient statistical evidence to suggest that both rs7968585 (HRR, 1.26 [95% CI, 0.95-1.68]; P = .11) and rs2239179 (HRR, 0.86 [95% CI, 0.64-1.16]; P = .08) independently modified the association of low vitamin D concentration with risk of the composite outcome.
We found a SNP within the VDR gene that significantly modified associations of low serum 25-hydroxyvitamin D concentration with major health outcomes of hip fracture, MI, cancer, and death over long-term follow-up. Findings were observed within a large community-based study of older adults in the United States and were consistent in magnitude and direction across individual disease outcomes, and replicated in a meta-analysis of 3 large independent cohorts. An additional VDR SNP significantly modified the low 25-hydroxyvitamin D–disease association in a meta-analysis that included results from the discovery and replication cohorts. The discovered SNPs, which are common in European populations, identified subsets of individuals for whom associations between low 25-hydroxyvitamin D concentration and disease outcomes were either strongly positive vs null. These results suggest that individuals with specific 25-hydroxyvitamin D metabolism genotypes may be particularly susceptible to, or protected from, the potential adverse health effects of low vitamin D.
The 2 identified VDR variants rs7968585 and rs2239179 are common SNPs, with minor allele frequencies of 0.48 and 0.42, respectively. The SNPs have moderate linkage (R2 = 0.55) and exploratory analyses provided insufficient evidence (P = .08) of independent interactions between each of these variants and low 25-hydroxyvitamin D concentration on risk of the composite outcome. These results make it unclear whether the interactions between the 2 observed SNPs and 25-hydroxyvitamin D concentration represent unique signals.
The rs7968585 and rs2239179 SNPs are not located in a coding region or GT-AG splice site.36 However, the 2 variants are in linkage disequilibrium with a cluster of common SNPs (BsmI, Taq I, and ApaI) at the 3′ end of the n VDR gene that are associated with clinical outcomes in other studies. The rs7968585 SNP is intergenic, located 3.2 kb downstream from VDR, and in moderate to high linkage disequilibrium with ApaI (rs7975232, R2 = 0.87), Taq I (rs731236, R2 = 0.65), and BsmI (rs1544410, R2 = 0.65).35 One recent study in postmenopausal women with osteoporosis found Taq I to be associated with differential response to calcitriol and calcium therapy.37 Candidate gene studies and meta-analyses link the ApaI and BsmI variants with the development of several types of cancer.38- 41 Specific ApaI, Taq I, and BsmI genotypes may also be related to decreased bone mineral density and elevated fracture risk.42- 44 The ApaI and Taq I variants were not genotyped in the CHS and thus were not candidates in the discovery phase. The SNP BsmI was genotyped in the CHS, and its interaction with low 25-hydroxyvitamin D concentration was of borderline statistical significance (HRR, 0.76; P = .02, q = 0.29). The other genetic variant identified in this study, rs2239179, is an intronic VDR SNP that may be related to the risk of renal cancer.45
There is biological plausibility that genetic variation within the vitamin D receptor could alter associations of 25-hydroxyvitamin D concentrations with disease outcomes. The VDR, a member of the steroid-receptor gene superfamily, plays a central role in mediating vitamin D signaling. The biologically active 1,25-dihydroxyvitamin D ligand binds to the VDR, which then forms a heterodimer with the retinoid X receptor. The VDR then enters the nucleus to bind vitamin D response elements in multiple target genes through a conserved DNA binding sequence.46 Genetic variation within the VDR could alter the response to 25-hydroxyvitamin D stores. For example, greater VDR affinity for 1,25-dihydroxyvitamin D or greater VDR activity for a given amount of 25-hydroxyvitamin D could provide protection in situations of low 25-hydroxyvitamin D substrate, such as residence in communities at high latitudes. These conjectures remain speculative because the biological effect of the discovered genetic variants on VDR function remains unknown.
Strengths of our study include the use of large, community-based cohort studies that include long-term follow-up, excellent retention, and independent adjudication of major health events to identify potentially important interactions. The use of a composite outcome captures many of the suspected biological effects of low vitamin D across a wide range of health processes. We performed a replication phase in 3 independent cohorts that included long-term follow-up for the same clinical disease events as the discovery cohort to confirm findings and obtain more reliable and precise estimates of the identified interactions.
Our study has several important limitations. First, there was some heterogeneity of the estimated interactions across the study cohorts, which could be explained by differences in the demographic, comorbidity, or 25-hydroxyvitamin D distributions at baseline, or by contrasting incidence rates of the composite outcome. Second, genetic information in our study was derived from SNPs measured in whole genome chip platforms, which incompletely characterize the target genes and cannot precisely identify potential causal variants. In addition, we lack functional data to explore why the identified VDR variants may modify 25-hydroxyvitamin D associations. With a single discovery cohort and candidate gene approach, we could have missed other important interactions between genotype and 25-hydroxyvitamin D concentration. Our assessment of low 25-hydroxyvitamin D concentrations is based on a single measurement that may not best reflect typical circulating concentrations over a longer period. Furthermore, the demographic group included in our study was primarily older white adults, and the findings may not be generalizable to other demographic or ethnic groups.
In summary, results of this candidate gene study indicate that known associations of low serum vitamin D concentration with clinical outcomes may vary according to genetic differences in the vitamin D receptor. These findings represent a first step toward identifying what may be clinically relevant effects of 25-hydroxyvitamin D metabolism genes and may contribute to a better understanding of the biological impact of genetic variation within the vitamin D receptor. Further studies are needed to confirm these observed associations and to enhance knowledge of how variation in vitamin D metabolism genes may stratify individuals as to their susceptibility to vitamin D deficiency. Evaluating the identified interactions in randomized clinical trials of vitamin D supplementation, when available, would help to assess the validity of our results and pave the way toward identifying individual patients who may benefit most from vitamin D interventions.
Corresponding Author: Gregory P. Levin, PhD, Department of Biostatistics, University of Washington, Box 357232, Seattle, WA 98195 (firstname.lastname@example.org).
Author Contributions: Dr Levin had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Levin, Robinson-Cohen, de Boer, Liu, Ferrucci, Psaty, Siscovick, Kestenbaum.
Acquisition of data: Levin, Houston, Liu, Kritchesvsky, Cauley, Bandinelli, Hagström, Michaëlsson, Melhus, Wang, Psaty, Siscovick, Kestenbaum.
Analysis and interpretation of data: Levin, Robinson-Cohen, de Boer, Houston, Lohman, Liu, Tanaka, Patel, Hagström, Michaëlsson, Wolf, Siscovick, Kestenbaum.
Drafting of the manuscript: Levin, Robinson-Cohen, Hagström, Kestenbaum.
Critical revision of the manuscript for important intellectual content: Levin, de Boer, Houston, Lohman, Liu, Kritchesvsky, Cauley, Tanaka, Ferrucci, Bandinelli, Patel, Hagström, Michaëlsson, Melhus, Wang, Wolf, Psaty, Siscovick, Kestenbaum.
Statistical analysis: Levin, Robinson-Cohen, Houston, Lohman, Liu, Tanaka, Hagström.
Obtained funding: Houston, Liu, Kritchesvsky, Michaëlsson, Melhus, Psaty, Siscovick, Kestenbaum.
Administrative, technical, or material support: Levin, Robinson-Cohen, Liu, Cauley, Bandinelli, Hagström, Michaëlsson, Wolf, Psaty.
Study supervision: de Boer, Kritchesvsky, Ferrucci, Kestenbaum.
Conflict of Interest Disclosures: The authors have completed and submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Dr de Boer reported receiving research grant funding from Abbott Laboratories. Dr Houston reported pending institutional grant support from the National Institute on Aging, National Institutes of Health; and receiving payment for lectures from the American Society for Nutrition, the Department of Veteran Affairs, and Abbott Nutrition Health Institute. Dr Wang reported serving as a consultant to Diasorin; receiving an investigator-initiated grant from Diasorin; and receiving payment for lectures from Diasorin. Dr Wolf reported serving as a consultant to Abbott, Genzyme, Luitpold, Mitsubishi, Cytochroma, Astellas, and Kai; receiving institutional grants from the National Institutes of Health, Shire, and Amgen; receiving payment for lectures from Abbott, Genzyme, and Shire; and having a pending patent. Dr Psaty reported serving on a data and safety monitoring board for a clinical trial of a device funded by Zoll ZifeCor; and serving on a steering committee for the Yale Open Data Access Project funded by Medtronic. Dr Kestenbaum reported receiving institutional grant support from Amgen. No other authors reported disclosures.
Funding/Support: The Cardiovascular Health Study was supported by contracts HHSN268201200036C, N01-HC-35129, N01-HC-45133, N01-HC-75150, N01-HC-85079 through N01-HC-85086, N01 HC-15103, N01 HC-55222, U01 HL080295, HL087652, and HL096875 from the National Heart, Lung, and Blood Institute; by funding from the National Institute of Neurological Disorders and Stroke; and by grant R01AG027002 from the National Institute on Aging. Handling of DNA and genotyping was supported in part by National Center for Research Resources grant M01RR00069 awarded to the Cedars-Sinai General Clinical Research Center Genotyping Core, and the National Institute of Diabetes and Digestive and Kidney Diseases grant DK063491 awarded to the Southern California Diabetes Endocrinology Research Center. Serum vitamin D measurements in the Cardiovascular Health Study were funded by National Institutes of Health grant R01 HL084443 (Dr Kestenbaum). The Health, Aging, and Body Composition study was supported by National Institute on Aging contracts N01AG62101, N01AG62103, and N01AG62106. Assessment of 25-hydroxyvitamin D levels was funded by National Institute on Aging grant R01-AG029364. The genome-wide association study was funded by National Institute on Aging grant R01-AG032098, and genotyping services were provided by the Center for Inherited Disease Research, which is fully funded through contract HHSN268200782096C from the National Institutes of Health to Johns Hopkins University. The baseline Invecchiare in Chianti study (1998-2000) was supported as a targeted project (ICS110.1/RF97.71) by the Italian Ministry of Health, in part by the US National Institute on Aging contracts 263 MD 9164, 263 MD 821336, N.1-AG-1-1, N.1-AG-1-2111, and N01-AG-5-0002, and in part by the Intramural Research Program of the National Institute on Aging, National Institutes of Health. The Uppsala Longitudinal Study of Adult Men was supported by the Swedish Research Council (2006), Swedish Heart-Lung Foundation, Thuréus Foundation, and Uppsala University.
Role of the Sponsors: None of the funding sources had any role in the design of the study; in the analysis and interpretation of the data; or in the preparation of the manuscript. The National Institutes of Health was involved in the original design and conduct of the Cardiovascular Health Study and the Health, Aging, and Body Composition study, and in the data collection methods and final approval of the manuscript.
Additional Contributions: We thank the other investigators, the staff, and the participants of the Cardiovascular Health Study, the Health, Aging, and Body Composition study, the Invecchiare in Chianti study, and the Uppsala Longitudinal Study of Adult Men for their valuable contributions.