Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genetic Variations in the Vitamin D Receptor Predict Type 2 Diabetes and Myocardial Infarction in a Community-Based Population: The Tromsø Study

  • Ieva Zostautiene ,

    ieva.zostautiene@unn.no

    Affiliations Tromsø Endocrine Research Group, Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Division of Internal Medicine, University Hospital of North Norway, Tromsø, Norway

    ⨯
  • Rolf Jorde,

    Affiliations Tromsø Endocrine Research Group, Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Division of Internal Medicine, University Hospital of North Norway, Tromsø, Norway

    ⨯
  • Henrik Schirmer,

    Affiliations Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Division of Cardiothoracic and Respiratory Medicine, University Hospital of North Norway, Tromsø, Norway

    ⨯
  • Ellisiv Bøgeberg Mathiesen,

    Affiliations Brain and Circulation Research Group, Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Department of Neurology and Neurophysiology, University Hospital of North Norway, Tromsø, Norway

    ⨯
  • Inger Njølstad,

    Affiliation Department of Community Medicine, UiT The Arctic University of Norway, Tromsø, Norway

    ⨯
  • Maja-Lisa Løchen,

    Affiliation Department of Community Medicine, UiT The Arctic University of Norway, Tromsø, Norway

    ⨯
  • Tom Wilsgaard,

    Affiliation Department of Community Medicine, UiT The Arctic University of Norway, Tromsø, Norway

    ⨯
  • Ragnar Martin Joakimsen,

    Affiliations Tromsø Endocrine Research Group, Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Division of Internal Medicine, University Hospital of North Norway, Tromsø, Norway

    ⨯
  • Elena Kamycheva

    Affiliations Tromsø Endocrine Research Group, Department of Clinical Medicine, UiT The Arctic University of Norway, Tromsø, Norway, Division of Internal Medicine, University Hospital of North Norway, Tromsø, Norway

    ⨯

Correction

19 Sep 2016: Martinaityte I, Jorde R, Schirmer H, Mathiesen EB, Njølstad I, et al. (2016) Correction: Genetic Variations in the Vitamin D Receptor Predict Type 2 Diabetes and Myocardial Infarction in a Community-Based Population: The Tromsø Study. PLOS ONE 11(9): e0163573. https://doi.org/10.1371/journal.pone.0163573 View correction

Abstract

Background

Though the associations between low serum 25-hydroxyvitamin D (25(OH)D) levels and health outcomes such as type 2 diabetes (T2D), myocardial infarction (MI), cancer, and mortality are well-studied, the effect of supplementation with vitamin D is uncertain. This may be related to genetic differences. Thus, rs7968585, a single nucleotide polymorphism (SNP) of the vitamin D receptor (VDR), has recently been reported as a predictor of composite health outcome. We therefore aimed to evaluate whether rs7968585 predicts separate clinical outcomes such as T2D, MI, cancer, and mortality in a community-based Norwegian population.

Methods and Findings

Measurements and DNA were obtained from the participants in the Tromsø Study in 1994–1995, registered with the outcomes of interest and a randomly selected control group. The impact of the rs7968585 genotypes was evaluated with Cox proportional hazards. A total of 8,461 subjects were included among whom 1,054 subjects were registered with T2D, 2,287 with MI, 3,166 with cancer, and 4,336 with death. Mean follow-up time from birth was 60.8 years for T2D and MI, 61.2 years for cancer, while mean follow-up time from examination date was 16.5 years for survival. Mean serum 25(OH)D levels did not differ across the rs7968585 genotypes. With the major homozygote genotype as reference, the minor homozygote subjects had hazard ratios of 1.25 (95% CI 1.05–1.49) for T2D and 1.14 (1.02–1.28) for MI (P = 0.011 and 0.023, respectively, without the Bonferroni correction). No significant interaction between serum 25(OH)D status and the rs7968585 genotype was found for any of the endpoints.

Conclusions

The VDR-related SNP rs7968585 minor allele is a significant and positive predictor for T2D and possibly for MI. Since the functional mechanism of this SNP is not yet understood, and the association with T2D is reported for the first time, confirmatory studies are needed.

Citation: Zostautiene I, Jorde R, Schirmer H, Mathiesen EB, Njølstad I, Løchen M-L, et al. (2015) Genetic Variations in the Vitamin D Receptor Predict Type 2 Diabetes and Myocardial Infarction in a Community-Based Population: The Tromsø Study. PLoS ONE 10(12): e0145359. https://doi.org/10.1371/journal.pone.0145359

Editor: Masaru Katoh, National Cancer Center, JAPAN

Received: July 14, 2015; Accepted: December 1, 2015; Published: December 23, 2015

Copyright: © 2015 Zostautiene et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: There are legal restrictions imposed on the data by the Norwegian Data Protection Authority that prevent these data from being made publicly available. All interested parties can request the data using the link: http://tromsoundersokelsen.uit.no/tromso.

Funding: The authors are indebted to the Norwegian Insitute of Public Health for their participation in collection of data in the fourth survey of the Tromsø Study (http://www.fhi.no/artikler/?id=28291), The Northern Norway Regional Health Authority (http://www.helse-nord.no/about-helse-nord/category5770.html; IZ grant No. SFP1215-14), and The Research Council of Norway (https://www.forskningsradet.no/prosjektbanken/#!/project/213787/no; RJ grant 213787). The funding institutions had no role in study design, data analysis and interpretation, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Vitamin D is a biologically active substance important not only for maintenance of calcium homeostasis but also for several nonskeletal metabolic pathways [1]. Despite the association between low serum 25-hydroxyvitamin D (25(OH)D) levels and increased risk of type 2 diabetes (T2D), cardiovascular diseases, cancer, and mortality [1–5], intervention studies with vitamin D thus far have been inconclusive [6]. Possibly, this may be due to genetic differences in vitamin D metabolism [7].

The vitamin D receptor (VDR) has been identified in most of the body’s tissues, including pancreatic islets, myocardium, and fibroblasts [8, 9]. The activated VDR regulates approximately 3% of all genes [10, 11] and may protect against diabetes by regulating insulin secretion and resistance and by reducing inflammatory damage to the pancreatic islets [12, 13]. Furthermore, the VDR regulates calcium homeostasis, may act antiatherosclerotic, may suppress the renin system, and may protect against myocardial hypertrophy and heart failure as well as coronary artery disease (CAD) [1, 14, 15]. In cancer pathogenesis, the activated VDR is suggested to act antiproliferatively and to regulate cell differentiation and apoptosis [16, 17]. Thus, the role of the VDR and its individual variations, instead of the serum concentration of vitamin D metabolites per se might be important in the pathogenesis of diabetes, myocardial infarction (MI), cancer, and mortality.

Many recent studies have focused on VDR polymorphisms as risk factors for these health outcomes, and some VDR single nucleotide polymorphisms (SNPs) are reported to be associated with risk of T2D [18], several forms of cancer including breast, prostate cancer, and malignant melanoma [19, 20], cancer-related mortality [21] and CAD [22–25], but not overall mortality [26].

Recently, the association between the common VDR gene variant rs7968585 and risk of a composite health outcome including hip fracture, MI, cancer, and mortality has been reported in subjects with low serum 25(OH)D levels [27]. However, in a subsequent study no associations between rs739837, a SNP in high linkage disequilibrium (LD) with rs7968585, and cardiometabolic outcomes were found [28].

The Tromsø Study is a large community-based study in northern Norway with repeated surveys every 6–7 years [29]. In the fourth survey performed in 1994–1995, lifestyle information and blood samples were collected from more than 27,000 subjects [30], and detailed endpoint registers were created. We were therefore able to evaluate the reported association between the rs7968585 genotypes, serum 25(OH)D and separate clinical outcomes such as T2D, MI, cancer, and mortality. Furthermore, since we also had genotype data on five other VDR SNPs (FokI (rs2228570/rs10735810), BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232) and Cdx2 (rs11568820)), we used the same approach to evaluate their association with the above mentioned outcomes.

Materials and Methods

The Tromsø Study

The Tromsø study initiated in 1974 is a longitudinal, community-based multipurpose Norwegian study focusing on lifestyle-related diseases [29]. It is performed in the Tromsø municipality, which is situated at 69° North with a current population of about 72,000 subjects. In the fourth survey in 1994–1995, all individuals aged 25 years or older were invited to participate; 77% participated, comprising a total of 27,158 subjects [30] of whom 26,956 were included in our analysis. A group of subjects were invited for a second visit for more extensive clinical examination and measurements, and 7,965 subjects, or 78% of the invited, participated [31].

Selection of the study subjects

The selection of subjects for the genetic studies has been described in detail previously [32]. Briefly, the subjects were primarily selected for evaluation of genetic polymorphisms and clinical endpoints registered in the Tromsø Study, including T2D, MI, cancer, and death. Since these endpoints were of potential interest regarding genetic polymorphisms, and limited funding did not allow genetic analyses of the entire Tromsø Study cohort, a case—cohort design was used with randomly selected controls from the entire cohort who participated in the fourth survey [33]. Accordingly, genotyping was performed in all subjects included in the endpoint registers as well as the randomly selected controls and was successful for 11,752 subjects. Since this control group was randomly selected from the entire cohort, the group included numerous subjects with one or more endpoints. When analyzing a specific endpoint, subjects in the control group with that specific endpoint were moved to the case group. Therefore, the size of the control groups varied depending on the endpoint in question.

Definition of T2D, MI, cancer, and death endpoints

The identification of subjects with endpoints has also previously been described in detail [34, 35]. In short, subjects with T2D were identified with self-reported diabetes questionnaires, elevated HbA1c (≥6.5%), and a search in medical records and national registries. Differentiation between T1D and T2D was based on clinical judgment and blood tests (C-peptide, anti-glutamic acid decarboxylase (anti-GAD)) when available. Cases of unknown types of diabetes were excluded from the analyses.

Subjects with MI were identified by searching the related diagnoses in all available registries and records using the World Health Organization (WHO) Multinational MONItoring of Trends and Determinants in CArdiovascular Disease/MOnica Risk, Genetics, Archiving and Monograph (MONICA/MORGAM) diagnostic criteria [36]. Silent MIs defined with electrocardiogram only were not included as cases because of difficulties in determining the exact date of the event. The T2D and MI endpoints were included until the end of 2011 but completely updated until the end of 2010.

Cancer events were obtained from the Cancer Registry of Norway, with endpoints registered until the end of 2010. The whole cancer group was also divided into four subgroups according to location: breast, lung, colorectum, and prostate.

Information about death was retrieved from the Causes of Death Registry, and information about participants who moved out of the Tromsø area was obtained from the Norwegian Registry of Vital Statistics, updated until January 2013.

Measurements

Baseline information was collected by physical examination, nonfasting blood samples, and self-administered questionnaires in 1994–1995. Blood pressure, weight, height, total cholesterol, and HbA1c were measured/analyzed as previously described [37]. Information on current smoking and current physical activity was collected from these questionnaires. Physical activity was defined as the presence or absence of light or hard physical activity during leisure time. Body mass index (BMI) was calculated as weight divided by height squared.

Sera, collected during the second visit of the Tromsø Study in 1994–1995, were stored at –70°C, thawed in 2008, and analyzed for serum 25(OH)D using an automated clinical chemistry analyzer (Modular E170, Roche).

DNA was prepared with the manual isolation method from whole blood samples collected during the first visit. Genotyping was performed using the KASP (KBioScience Allele-Specific Polymorphism) SNP genotyping system as previously described [32].

Statistical analyses

Distribution of the variables serum 25(OH)D, systolic blood pressure, BMI, total cholesterol, and HbA1c were evaluated with skewness, kurtosis, and visual inspection of histogram and found normal except for HbA1c, which, however, could not be normalized with log transformation. Trends across the VDR SNP genotypes were evaluated with linear regression for continuous variables, using age, gender, and for 25(OH)D also season (months, using dummy variables), as covariates. For categorical variables, trends across genotypes were evaluated with the chi-square test with linear-by-linear association, and for HbA1c, the Kruskal—Wallis test was used for testing differences between the genotypes. Since serum 25(OH)D levels in smokers are overestimated when analyzed with the Modular E170 (Roche) [38], smokers and nonsmokers were analyzed separately unless otherwise specified. The Student’s t-test, Mann-Whitney U test, or Pearson’s chi-square test was used to compare variables in endpoint groups and their controls.

The genotype frequencies were evaluated for Hardy-Weinberg equilibrium using the chi-square test [39].

The hazard ratio (HR) of the VDR SNP genotypes for T2D, MI, cancer, and mortality was estimated with Cox regression with the major homozygote genotype used as the reference. Observation time was from birth until the event of interest for all endpoints except for death, where the observation time was from 1994.

In Model 1, for all endpoints, an adjustment for age and gender was made. In Model 2 for T2D and MI additional adjustment for data obtained in 1994/1995 on BMI, systolic blood pressure, total cholesterol, current smoking, and physical activity was performed. For MI, T2D was additionally included in Model 3. Since the effects of genotypes on cancer and mortality risk were not significant in any model regardless of covariates included, only Model 1 results were presented.

Additionally, serum 25(OH)D levels below or above the 20th percentile as well as serum 25(OH)D used as a continuous variable were evaluated for association with T2D, MI, cancer, and death with age and gender as covariates. Similarly, interaction between serum 25(OH)D below or above the 20th percentile as well as 25(OH)D as a continuous variable and the VDR SNP genotypes were also tested with Cox regression with age and gender as covariates. In the analysis of a specific endpoint other endpoints were not included as covariates unless related to rs7968585 genotypes which was the case only for MI and T2D.

In the 20th percentile analyses, smokers and nonsmokers were analyzed together using month-specific 20th percentile cutoffs calculated separately for smokers and nonsmokers. When 25(OH)D was used as a continuous variable, smokers and nonsmokers were analyzed separately. These analyses with 25(OH)D were performed on incident cases only (observation time from Tromsø 4, 1994–1995).

The data are shown as mean ± standard deviation (SD). All tests are presented two-sided, and P-value <0.05 was considered statistically significant. The data were analyzed with IBM SPSS Statistics 22 (SPSS Inc., Chicago, IL, USA). The data are presented without correction for multiple analyses unless specified in the text if the Bonferroni correction was applied (P-value multiplied by a factor of 4 (four main endpoints)).

Ethics

The study was approved by the Regional Committee for Medical and Health Research Ethics (REK Nord) (reference 2010/2913-4). Only participants with valid written consent were included.

Results

The baseline characteristics of the entire Tromsø Study and the endpoint groups relevant for those successfully genotyped for rs7968585 are shown in Table 1. Serum 25(OH)D values both for smokers and nonsmokers were lower in the T2D and mortality groups than in the respective control groups. There were significantly more men and smokers in all case groups compared to the respective control groups. The rs7968585 genotype frequencies were not found to deviate.

thumbnail
Download:
Table 1. Baseline characteristics in 1994–1995 in the entire Tromsø Study population and in the endpoint groups related to the rs7968585 analyses. The Tromsø Study.

https://doi.org/10.1371/journal.pone.0145359.t001

Age, sex, systolic blood pressure, BMI, total cholesterol, HbA1c, physical activity, and serum 25(OH)D values in the 11,752 subjects successfully genotyped for rs7968585 are presented in Table 2 with no significant trends observed across the genotypes.

thumbnail
Download:
Table 2. Baseline characteristics of the 11,752 subjects in 1994 according to rs7968585 genotype. The Tromsø Study.

https://doi.org/10.1371/journal.pone.0145359.t002

Among the 11,752 successfully genotyped subjects for rs7968585, 4,198 had been randomly selected for the control cohort, 1,054 subjects were registered with T2D, 2,287 subjects had confirmed MI, 3,166 had cancer, among whom there were 431 breast, 385 lung, 501 colorectal, and 406 prostate cancer cases, and 4,336 deaths. Mean follow-up time from birth was 60.8 years for T2D and MI, 61.2 years for cancer, while the mean follow-up time from examination date was 16.5 years for death. The distribution of major endpoints within the rs7968585 genotype groups is presented in Table 3, with significant trend across genotypes for T2D and MI (P<0.05).

thumbnail
Download:
Table 3. Distribution of endpoints in rs7968585 genotype groups. The Tromsø Study.

https://doi.org/10.1371/journal.pone.0145359.t003

The HRs for all endpoints in regard to rs7968585 genotype are shown in Tables 4 and 5. Considering major homozygote (T:T) as a reference, subjects with the minor homozygote (C:C) had a significantly (P = 0.044 after Bonferroni correction) 25% increased risk of developing T2D in Models 1 and 2. Similarly for MI, subjects with the minor homozygotes had a 14% and 13% increased risk in Models 1 and 2, respectively, but the P-value was >0.05 after the Bonferroni correction for multiple testing. Furthermore, when T2D was included as a covariate, the increased risk of MI in subjects with the minor homozygote genotype was reduced to 11% (P = 0.080). No significant impact of the rs7968585 genotypes on the risk of total cancer (as well as for the four subtypes breast, lung, colorectal, and prostate, data not shown) and mortality was observed.

thumbnail
Download:
Table 4. Hazard ratios (HR) for the rs7968585 genotypes regarding T2D, MI, cancer and mortality analyzed with the Cox regression, adjusted for age and gender. The Tromsø Study.

https://doi.org/10.1371/journal.pone.0145359.t004

thumbnail
Download:
Table 5. Hazard ratios (HR) for the rs7968585 genotypes regarding T2D and MI analyzed with the Cox regression, adjusted for age, gender and additional risk factors. The Tromsø Study.

https://doi.org/10.1371/journal.pone.0145359.t005

There were ~8,500 subjects successfully genotyped for FokI, BsmI, TaqI, ApaI and Cdx2 groups within the control cohort and the endpoint groups. However, we did not observe significant trends across their genotypes for serum 25(OH)D status. Furthermore, none of the FokI, BsmI, TaqI, ApaI and Cdx2 genotypes showed significantly increased risk of any of the four endpoints in Cox regression analyses after Bonferroni correction (data not shown).

Serum 25(OH)D status was significantly associated with T2D, MI, and death but not cancer. When the cohort was divided according to the serum 25(OH)D 20th percentile, the subjects in the low serum 25(OH)D group had an 73% increased risk of T2D (95% confidence interval (CI) 1.40–2.14, P<0.01), a 20% increased risk of MI (95% CI 1.02–1.41, P<0.05), and a 21% increased risk of death (95% CI 1.09–1.34, P<0.01). Same trend was observed analyzing serum 25(OH)D as a continuous variable: HR per SD decrease was 1.45 for T2D in smokers, 1.26 in nonsmokers. For MI, HR per SD decrease was 1.14 in smokers, not significant in nonsmokers. For mortality, HR per SD was 1.10 in smokers and 1.07 in nonsmokers.

The 20th percentile for serum 25(OH)D was 33.9–46.6 nmol/L in nonsmokers and 50.4–64.2 nmol/L in smokers, differing according to the month. The number of subjects with serum 25(OH)D below the 20th percentile (month and smoking specific) was 115, 182, 188, and 471 for the T2D, MI, cancer, and mortality case groups, respectively. However, there were no significant interactions between serum 25(OH)D status (above/below 20th percentile, as well as when serum 25(OH)D was used as a continuous variable) and rs7968585 genotype regarding any of the endpoints.

In particular, for T2D, the positive association with low serum 25(OH)D (below the 20th percentile compared to those above) was seen for subjects in all three rs7968585 genotypes with an increased risk of 123% for major homozygotes (95% CI 1.54–3.22, P<0.01) 46% for heterozygotes (95% CI 1.30–2.19, P<0.05), and 79% for minor homozygotes (95% CI 1.10–2.92, P<0.05).

Discussion

In our large, community-based study of 8,461 subjects, we found that the minor allele (C) at rs7968585 is a significant risk factor for T2D, possibly also for MI, but not for cancer or death. To our knowledge, the association between rs7968585 and T2D is reported for the first time; meanwhile, the relation between this genotype and MI was probable, based on previous observation by Levin et al., where the risk of a composite health outcome including MI was found to be modified by the minor allele of the rs7968585 genotype in subjects with low serum 25(OH)D levels [27].

The 25% increased risk of T2D in the rs7968585 minor homozygote subjects contrasts the results from the 1958 British Birth cohort study where the VDR SNP rs739837 (which is in high LD with rs7968585 (r2 = 0.87)) was not associated with T2D [28]. In addition to the possible effect of the slight genetic difference between the two SNPs, one possible explanation could be the lack of power in the 1958 British Birth cohort study that included only 5,160 subjects. In addition, the subjects were considerably younger than in our study, and therefore, probably had fewer cases with T2D included. In support of our findings, there was a positive, although not statistically significant, association between minor homozygote and T2D.

In addition to T2D, we also observed a 14% increased risk of MI in the minor homozygote subjects. This increased risk was significant only before correction for multiple testing, but as the pathogeneses of the endpoints probably are biologically separate, the necessity for the correction might be questionable. However, when T2D was included as a confounder in the Cox regression for the risk of MI, the association between rs7968585 and MI was decreased, indicating that the genetic influence on MI was partly mediated by its effects on glucose metabolism. The association between rs7968585 and MI was also reported by Levin et al. from the Cardiovascular Health Study Discovery Cohort that included 1,514 subjects with 214 cases with MI [27]. Similar to our study, Levin et al. found a 20% increased risk for each additional minor allele, but this did not reach statistical significance. However, when viewing the two studies combined, a true association appears to be likely.

However, we could not replicate the positive association between the minor homozygote rs7968585 and cancer as well as death reported in the Cardiovascular Health Study Discovery Cohort [27]. We found only a slight nonpositive association with cancer but found no relevant association for death. Accordingly, more studies are needed to establish associations between rs7968585 and these clinical endpoints.

If such associations for rs7968585 are established, the plausible biological explanation would be through an effect on the VDR [40]. Thus, rs7968585 is located in a noncoding region close to the VDR gene and is in high to moderate LD with the most studied VDR SNPs Apal (r2 = 0.87), Taql (r2 = 0.65) and Bsml (r2 = 0.65) [41, 42]. The observed associations between rs7968585 and T2D as well as MI could therefore reflect covariation with these or other SNPs as well as effects on the VDR via still not proved coding of the VDR or the noncoding RNAs [43].

Thus, the nuclear receptor VDR is expressed in pancreatic cells, is involved in regulating insulin resistance, secretion, and inflammation of pancreatic islets [8, 12, 13], and may accordingly play a role in pathogenesis of T2D. The VDR is also expressed in cardiomyocytes, vessels, and other tissues involved in CAD pathogenesis, and when activated seems to protect for CAD [1, 14, 15].

Finally, the VDR seems to modify cancer-related processes by regulating cell proliferation, differentiation, and apoptosis [16, 17]. In line with this, the commonly studied VDR SNP FokI (in low LD with rs7968585) appears to be associated with T2D, at least in Asian populations [13, 18] and cancer [44, 45]; and the SNPs ApaI, BsmI, as well as TaqI appear to be associated with CAD [22–24] as well as with cancer [44, 45] and cancer-related mortality [21]; and some Cdx2 (in low LD with rs7968585) haplotypes appear to be associated with increased risk of cancer [46]. However, in our study, none of the FokI, BsmI, TaqI, ApaI, Cdx2 SNPs showed any significant associations with the risk for these outcomes, which might be due to inter-population genetic variances.

Another major objective of our study was to evaluate possible interactions between rs7968585 and the serum 25(OH)D levels as described by Levin et al. regarding their composite endpoint [27]. This interaction is plausible given that the rs7968585 minor allele makes the VDR less active. If so, for subjects with minor alleles, the presumed positive effects by the VDR would need higher serum 25(OH)D levels to be activated [40]. However, in spite of the expected significant associations between serum 25(OH)D and subsequent T2D, MI and death [1, 6], no interactions between rs7968585 and serum 25(OH)D levels regarding the clinical endpoints were found, which might be due to generally sufficient serum 25(OH)D levels in our population (48.5–52.7 nmol/L in nonsmokers) in contrast to the study of Levin et al. [27]. Thus, in our study the negative effects of the minor allele were seen regardless of serum 25(OH)D levels, and a similar conclusion was reached by Vimaleswaran et al. regarding serum 25(OH)D and T2D in their British Birth cohort study [28]. One probable explanation could be that rs7968585 acts through modulating the response to vitamin D, independently from vitamin D status as suggested by Barry et al. [47]. However, given the attractiveness of such an interaction between serum 25(OH)D and the VDR, this should be tested not only for rs7968585 but also in all studies in which both genetic as well as serum 25(OH)D data are available.

Our study has several limitations. Firstly, after the Bonferroni correction, the association between rs7968585 and the endpoints was significant only for T2D. This may reflect that for this type of study even a cohort of more than 8,000 subjects may be too small. Secondly, not all subjects had serum 25(OH)D measurements, which substantially reduced the power for the interaction part of the study. Thirdly, we relied on only one serum 25(OH)D measurement, and even though there is a high degree of tracking for serum 25(OH)D [48], several consecutive measurements would have been a great advantage. Finally, we do not know the consequences of the rs7968585 SNP for the function of the VDR, and it is therefore difficult to fit this SNP into a biological context.

Our study also has several strengths. The endpoint registers were of high quality and community-based in a large north Norwegian population with long-term follow-up. We also had the opportunity to take advantage of a previous publications regarding rs7968585 and therefore focused on this SNP only, which increased the statistical power.

In conclusion, we found a significant association between the VDR SNP rs7968585 and risk of T2D and possibly also MI in a general population, whereas interactions between rs7968585 and serum 25(OH)D levels regarding clinical endpoints were not found. Confirmatory studies are needed.

Acknowledgments

We are indebted to the Norwegian Institute of Public Health for their participation in collection of data in the fourth survey of the Tromsø Study, The Northern Norway Regional Health Authority, and The Research Council of Norway.

Author Contributions

Conceived and designed the experiments: IZ RJ. Performed the experiments: HS EBM IN M-LL TW RMJ. Analyzed the data: IZ RJ TW EK. Contributed reagents/materials/analysis tools: IZ RJ HS TW RMJ. Wrote the paper: IZ. Revised the manuscript critically for important intellectual content: RJ HS EBM IN M-LL TW RMJ EK.

References

  1. 1. Pludowski P, Holick MF, Pilz S, Wagner CL, Hollis BW, Grant WB, et al. Vitamin D effects on musculoskeletal health, immunity, autoimmunity, cardiovascular disease, cancer, fertility, pregnancy, dementia and mortality-a review of recent evidence. Autoimmunity reviews. 2013;12(10):976–89. pmid:23542507
  2. 2. Anderson JL, May HT, Horne BD, Bair TL, Hall NL, Carlquist JF, et al. Relation of vitamin D deficiency to cardiovascular risk factors, disease status, and incident events in a general healthcare population. The American journal of cardiology. 2010;106(7):963–8. pmid:20854958
  3. 3. Hutchinson MS, Figenschau Y, Njolstad I, Schirmer H, Jorde R. Serum 25-hydroxyvitamin D levels are inversely associated with glycated haemoglobin (HbA(1c)). The Tromso Study. Scandinavian journal of clinical and laboratory investigation. 2011;71(5):399–406. pmid:21517716
  4. 4. Schottker B, Haug U, Schomburg L, Kohrle J, Perna L, Muller H, et al. Strong associations of 25-hydroxyvitamin D concentrations with all-cause, cardiovascular, cancer, and respiratory disease mortality in a large cohort study. The American journal of clinical nutrition. 2013;97(4):782–93. pmid:23446902
  5. 5. Hutchinson MS, Grimnes G, Joakimsen RM, Figenschau Y, Jorde R. Low serum 25-hydroxyvitamin D levels are associated with increased all-cause mortality risk in a general population: the Tromso study. European journal of endocrinology / European Federation of Endocrine Societies. 2010;162(5):935–42. pmid:20185562
  6. 6. Theodoratou E, Tzoulaki I, Zgaga L, Ioannidis JP. Vitamin D and multiple health outcomes: umbrella review of systematic reviews and meta-analyses of observational studies and randomised trials. Bmj. 2014;348:g2035. pmid:24690624
  7. 7. Saccone D, Asani F, Bornman L. Regulation of the vitamin D receptor gene by environment, genetics and epigenetics. Gene. 2015;561(2):171–80. pmid:25682935
  8. 8. Chuang JC, Cha JY, Garmey JC, Mirmira RG, Repa JJ. Research resource: nuclear hormone receptor expression in the endocrine pancreas. Molecular endocrinology. 2008;22(10):2353–63. pmid:18669644
  9. 9. Chen S, Glenn DJ, Ni W, Grigsby CL, Olsen K, Nishimoto M, et al. Expression of the vitamin d receptor is increased in the hypertrophic heart. Hypertension. 2008;52(6):1106–12. pmid:18936343
  10. 10. Hossein-nezhad A, Spira A, Holick MF. Influence of vitamin D status and vitamin D3 supplementation on genome wide expression of white blood cells: a randomized double-blind clinical trial. PloS one. 2013;8(3):e58725. pmid:23527013
  11. 11. 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(9736):180–8. pmid:20541252
  12. 12. Pilz S, Kienreich K, Rutters F, de Jongh R, van Ballegooijen AJ, Grubler M, et al. Role of vitamin D in the development of insulin resistance and type 2 diabetes. Current diabetes reports. 2013;13(2):261–70. pmid:23264189
  13. 13. Van Belle TL, Gysemans C, Mathieu C. Vitamin D and diabetes: the odd couple. Trends in endocrinology and metabolism: TEM. 2013;24(11):561–8. pmid:23972673
  14. 14. Gangula PR, Dong YL, Al-Hendy A, Richard-Davis G, Montgomery-Rice V, Haddad G, et al. Protective cardiovascular and renal actions of vitamin D and estrogen. Frontiers in bioscience. 2013;5:134–48.
  15. 15. Gardner DG, Chen S, Glenn DJ. Vitamin D and the heart. American journal of physiology Regulatory, integrative and comparative physiology. 2013;305(9):R969–77. pmid:24026071
  16. 16. Deeb KK, Trump DL, Johnson CS. Vitamin D signalling pathways in cancer: potential for anticancer therapeutics. Nature reviews Cancer. 2007;7(9):684–700. pmid:17721433
  17. 17. Hansen CM, Binderup L, Hamberg KJ, Carlberg C. Vitamin D and cancer: effects of 1,25(OH)2D3 and its analogs on growth control and tumorigenesis. Frontiers in bioscience: a journal and virtual library. 2001;6:D820–48.
  18. 18. Li L, Wu B, Liu JY, Yang LB. Vitamin D receptor gene polymorphisms and type 2 diabetes: a meta-analysis. Archives of medical research. 2013;44(3):235–41. pmid:23506721
  19. 19. Ahn J, Albanes D, Berndt SI, Peters U, Chatterjee N, Freedman ND, et al. Vitamin D-related genes, serum vitamin D concentrations and prostate cancer risk. Carcinogenesis. 2009;30(5):769–76. pmid:19255064
  20. 20. Kostner K, Denzer N, Muller CS, Klein R, Tilgen W, Reichrath J. The relevance of vitamin D receptor (VDR) gene polymorphisms for cancer: a review of the literature. Anticancer research. 2009;29(9):3511–36. pmid:19667145
  21. 21. Buttigliero C, Monagheddu C, Petroni P, Saini A, Dogliotti L, Ciccone G, et al. Prognostic role of vitamin d status and efficacy of vitamin D supplementation in cancer patients: a systematic review. The oncologist. 2011;16(9):1215–27. pmid:21835895
  22. 22. Ortlepp JR, Krantz C, Kimmel M, von Korff A, Vesper K, Schmitz F, et al. Additive effects of the chemokine receptor 2, vitamin D receptor, interleukin-6 polymorphisms and cardiovascular risk factors on the prevalence of myocardial infarction in patients below 65 years. International journal of cardiology. 2005;105(1):90–5. pmid:16207551
  23. 23. Van Schooten FJ, Hirvonen A, Maas LM, De Mol BA, Kleinjans JC, Bell DA, et al. Putative susceptibility markers of coronary artery disease: association between VDR genotype, smoking, and aromatic DNA adduct levels in human right atrial tissue. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 1998;12(13):1409–17.
  24. 24. Ferrarezi DA, Bellili-Munoz N, Dubois-Laforgue D, Cheurfa N, Lamri A, Reis AF, et al. Allelic variations of the vitamin D receptor (VDR) gene are associated with increased risk of coronary artery disease in type 2 diabetics: the DIABHYCAR prospective study. Diabetes & metabolism. 2013;39(3):263–70.
  25. 25. Vaidya A, Sun B, Forman JP, Hopkins PN, Brown NJ, Kolatkar NS, et al. The Fok1 vitamin D receptor gene polymorphism is associated with plasma renin activity in Caucasians. Clinical endocrinology. 2011;74(6):783–90. pmid:21521263
  26. 26. Perna L, Hoffmeister M, Schottker B, Arndt V, Haug U, Holleczek B, et al. Vitamin D receptor polymorphism and colorectal cancer-specific and all-cause mortality. Cancer epidemiology. 2013;37(6):905–7. pmid:24075799
  27. 27. 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):1898–905. pmid:23150009
  28. 28. Vimaleswaran KS, Power C, Hypponen E. Interaction between vitamin D receptor gene polymorphisms and 25-hydroxyvitamin D concentrations on metabolic and cardiovascular disease outcomes. Diabetes & metabolism. 2014;40(5):386–9.
  29. 29. Thelle DS, Foorde OH, Try K, Lehmann EH. The Tromsoo heart study. Methods and main results of the cross-sectional study. Acta medica Scandinavica. 1976;200(1–2):107–18. pmid:785953
  30. 30. Jacobsen BK, Eggen AE, Mathiesen EB, Wilsgaard T, Njolstad I. Cohort profile: the Tromso Study. International journal of epidemiology. 2012;41(4):961–7. pmid:21422063
  31. 31. Berntsen GK, Fonnebo V, Tollan A, Sogaard AJ, Magnus JH. Forearm bone mineral density by age in 7,620 men and women: the Tromso study, a population-based study. American journal of epidemiology. 2001;153(5):465–73. pmid:11226978
  32. 32. Jorde R, Schirmer H, Wilsgaard T, Joakimsen RM, Mathiesen EB, Njolstad I, et al. Polymorphisms related to the serum 25-hydroxyvitamin D level and risk of myocardial infarction, diabetes, cancer and mortality. The Tromso Study. PloS one. 2012;7(5):e37295. pmid:22649517
  33. 33. Kulathinal S, Karvanen J, Saarela O, Kuulasmaa K. Case-cohort design in practice—experiences from the MORGAM Project. Epidemiologic perspectives & innovations: EP+I. 2007;4:15.
  34. 34. Joseph J, Svartberg J, Njolstad I, Schirmer H. Incidence of and risk factors for type-2 diabetes in a general population: the Tromso Study. Scandinavian journal of public health. 2010;38(7):768–75. pmid:20696770
  35. 35. Mannsverk J, Wilsgaard T, Njolstad I, Hopstock LA, Lochen ML, Mathiesen EB, et al. Age and gender differences in incidence and case fatality trends for myocardial infarction: a 30-year follow-up. The Tromso Study. European journal of preventive cardiology. 2012;19(5):927–34.
  36. 36. Evans A, Salomaa V, Kulathinal S, Asplund K, Cambien F, Ferrario M, et al. MORGAM (an international pooling of cardiovascular cohorts). International journal of epidemiology. 2005;34(1):21–7. pmid:15561751
  37. 37. Borch KH, Braekkan SK, Mathiesen EB, Njolstad I, Wilsgaard T, Stormer J, et al. Abdominal obesity is essential for the risk of venous thromboembolism in the metabolic syndrome: the Tromso study. Journal of thrombosis and haemostasis: JTH. 2009;7(5):739–45. pmid:19036065
  38. 38. Grimnes G, Almaas B, Eggen AE, Emaus N, Figenschau Y, Hopstock LA, et al. Effect of smoking on the serum levels of 25-hydroxyvitamin D depends on the assay employed. European journal of endocrinology / European Federation of Endocrine Societies. 2010;163(2):339–48. pmid:20479012
  39. 39. Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. American journal of epidemiology. 2009;169(4):505–14. pmid:19126586
  40. 40. Cheskis B, Freedman LP. Ligand modulates the conversion of DNA-bound vitamin D3 receptor (VDR) homodimers into VDR-retinoid X receptor heterodimers. Molecular and cellular biology. 1994;14(5):3329–38. pmid:8164684
  41. 41. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, de Bakker PI. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24(24):2938–9. pmid:18974171
  42. 42. Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, et al. The NCBI BioSystems database. Nucleic acids research. 2010;38(Database issue):D492–6. pmid:19854944
  43. 43. Pescador N, Perez-Barba M, Ibarra JM, Corbaton A, Martinez-Larrad MT, Serrano-Rios M. Serum circulating microRNA profiling for identification of potential type 2 diabetes and obesity biomarkers. PloS one. 2013;8(10):e77251. pmid:24204780
  44. 44. Gandini S, Gnagnarella P, Serrano D, Pasquali E, Raimondi S. Vitamin D receptor polymorphisms and cancer. Advances in experimental medicine and biology. 2014;810:69–105. pmid:25207361
  45. 45. Xu Y, He B, Pan Y, Deng Q, Sun H, Li R, et al. Systematic review and meta-analysis on vitamin D receptor polymorphisms and cancer risk. Tumour biology: the journal of the International Society for Oncodevelopmental Biology and Medicine. 2014;35(5):4153–69.
  46. 46. Slattery ML, Herrick J, Wolff RK, Caan BJ, Potter JD, Sweeney C. CDX2 VDR polymorphism and colorectal cancer. Cancer epidemiology, biomarkers & prevention: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2007;16(12):2752–5.
  47. 47. Barry EL, Rees JR, Peacock JL, Mott LA, Amos CI, Bostick RM, et al. Genetic variants in CYP2R1, CYP24A1, and VDR modify the efficacy of vitamin D3 supplementation for increasing serum 25-hydroxyvitamin D levels in a randomized controlled trial. The Journal of clinical endocrinology and metabolism. 2014;99(10):E2133–7. pmid:25070320
  48. 48. Jorde R, Sneve M, Hutchinson M, Emaus N, Figenschau Y, Grimnes G. Tracking of serum 25-hydroxyvitamin D levels during 14 years in a population-based study and during 12 months in an intervention study. American journal of epidemiology. 2010;171(8):903–8. pmid:20219763