Skip to main content

Real-life use of vitamin D3-fortified bread and milk during a winter season: the effects of CYP2R1 and GC genes on 25-hydroxyvitamin D concentrations in Danish families, the VitmaD study

Abstract

Common genetic variants rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC and a combined genetic risk score (GRS) of these four variants influence late summer 25-hydroxyvitamin D (25(OH)D) concentrations. The objectives were to identify those who are most at risk of developing low vitamin D status during winter and to assess whether vitamin D3-fortified bread and milk will increase 25(OH)D concentrations in those with genetically determined low 25(OH)D concentrations at late summer. We used data from the VitmaD study. Participants were allocated to either vitamin D3-fortified bread and milk or non-fortified bread and milk during winter. In the fortification group, CYP2R1 (rs10741657) and GC (rs4588 and rs842999) were statistically significantly associated with winter 25(OH)D concentrations and CYP2R1 (rs10766197) was borderline significant. There was a negative linear trend between 25(OH)D concentrations and carriage of 0–8 risk alleles (p < 0.0001). No association was found for the control group (p = 0.1428). There was a significant positive linear relationship between different quintiles of total vitamin D intake and the increase in 25(OH)D concentrations among carriers of 0–2 (p = 0.0012), 3 (p = 0.0001), 4 (p = 0.0118) or 5 (p = 0.0029) risk alleles, but not among carriers of 6–8 risk alleles (p = 0.1051). Carriers of a high GRS were more prone to be vitamin D deficient compared to carriers of a low GRS. Furthermore, rs4588-AA carriers have a low but very stable 25(OH)D concentration, and interestingly, also low PTH level.

Introduction

In northern latitudes (>40°N), low vitamin D status in humans, measured as 25-hydroxyvitamin D (25(OH)D) concentrations, is common during winter months. This is because vitamin D cannot be synthesized in the skin due to the lack of solar ultraviolet B radiation (UVB) and because the average dietary intake of vitamin D is insufficient (Thuesen et al. 2012). Moreover, twin- and family-based studies indicate that genetic factors may influence 25(OH)D concentrations appreciably (Engelman et al. 2008; Shea et al. 2009; Karohl et al. 2010). Two genome-wide association studies (GWAS) and several candidate gene studies have shown single-nucleotide polymorphisms (SNPs) to influence 25(OH)D concentrations (Engelman et al. 2008; Sinotte et al. 2009; Bogh et al. 2010; Ahn et al. 2010; Bu et al. 2010; Zhang et al. 2012; Monticielo et al. 2012; Engelman et al. 2013; Zhang et al. 2013; Nissen et al. 2014). These SNPs are located in the group-specific component also known as Gc globulin (GC) and in or near genes involved in vitamin D synthesis, activation or degradation. These findings indicate that 25(OH)D concentrations do not only depend on vitamin D intake and sun exposure, but also on genetic factors. Thus, genetic factors may help to identify individuals at risk of low vitamin D status.

We have previously found genetic variants in CYP2R1 (rs10741657 and rs10766197) and GC (rs4588 and rs842999) genes to predict late summer 25(OH)D concentrations in Danish families in a study of 25 SNPs in vitamin D metabolism (Nissen et al. 2014). The main focus of this study is therefore on the influence of rs10741657 and rs10766197 in CYP2R1, and rs842999 and rs4588 in GC on 25(OH)D concentrations in participants allocated to either vitamin D3-fortified bread and milk or non-fortified bread and milk during winter.

CYP2R1, a member of the cytochrome P450 family of enzymes, is the primary enzyme that hydroxylates vitamin D to 25(OH)D in the liver. Genetic variants of the CYP2R1 gene are strongly associated with 25(OH)D concentration (Wjst et al. 2006; Ramos-lopez and Brück 2007; Bu et al. 2010; Zhang et al. 2012, 2013; Nissen et al. 2014) and reached a high score in two GWAS (Ahn et al. 2010; Wang et al. 2010). Furthermore, Ahn et al. (Ahn et al. 2010) observed heterogeneity between different cohorts in the GWAS and the association of 25(OH)D concentration with CYP2R1. A missense mutation in CYP2R1 in exon 2 (L99P) is known to lead to vitamin D deficiency (Cheng et al. 2004).

Genetic variants in the GC gene reached the highest score in two GWAS (Ahn et al. 2010; Wang et al. 2010), and several candidate gene studies have found association with 25(OH)D concentrations (Lauridsen et al. 2005; Kurylowicz et al. 2006; Abbas et al. 2008; Engelman et al. 2008; Sinotte et al. 2009; Fu et al. 2009; Gozdzik et al. 2011; Lu et al. 2012; Nissen et al. 2014).

The GC gene encodes the vitamin D-binding protein (DBP), which is the primary vitamin D carrier protein. DBP binds with high affinity 85–90 % of circulating 25(OH)D, albumin binds with low affinity 10–15 % of circulating 25(OH)D and less than 1 % of 25(OH)D is in the free form (Bikle et al. 1986). The main function of DBP is to stabilize and prolong the half-life of 25(OH)D and other vitamin D metabolites (Speeckaert et al. 2006). DBP has several other important biological functions including fatty acid transportation, extracellular actin scavenging, leucocyte C5a-mediated chemotaxis, macrophage activation and stimulation of osteoclasts (Pekkinen et al. 2014).

The most studied GC SNPs are rs4588 and rs7041 that give rise to three common DBP isoforms, GC1F (rs7041-T, rs4588-C), GC1S (rs7041-G, rs4588-C) and GC2 (rs7041-T, rs4588-A), which differ by amino acid composition and glycosylation (Gozdzik et al. 2011). Vitamin D status differed significantly depending on rs4588 (or rs2282679, r 2 > 0.99) and/or rs7041 genotypes, where the A-allele of rs4588 and/or the T-allele of rs7041 were consistently associated with lower 25(OH)D concentrations (Lauridsen et al. 2005; Kurylowicz et al. 2006; Abbas et al. 2008; Engelman et al. 2008; Sinotte et al. 2009; Fu et al. 2009; Gozdzik et al. 2011; Lu et al. 2012). In Caucasian, rs4588 and rs7041 are in almost complete linkage disequilibrium (LD) (Haploview software version 4.2). There is biological support that the affinity to both 25(OH)D and 1,25(OH)2D is higher for the rs4588 C-allele isoform than for the A-allele isoform (Arnaud and Constans 1993). Based on glycosylation patterns, it is suggested that the GC2 phenotype is fast metabolizer. Kawakami et al. (1981) observed that the metabolic rate indeed was higher in GC2-2 individuals than in GC1-1 individuals. In addition, the GC2 genotype, which is associated with lower 25(OH)D concentrations, is also associated with low mean DBP concentration (Lauridsen et al. 2005). The GC2 and GC1S isoforms are more frequent in people with light skin whereas the GC1F isoform is more frequent in people with dark skin (Kamboh and Ferrell 1986).

Measurement of 25(OH)D concentration in blood is currently the best biological marker of vitamin D status and reflects total vitamin D exposure—from diet, supplements and cutaneous synthesis. Severe vitamin D deficiency (<12 nmol/L) is a medical condition associated with osteomalacia in adults and rickets in children. Vitamin D deficiency can lead to osteoporosis due to increased bone resorption caused by increased serum concentrations of parathyroid hormone (PTH) (Holick 2007). Moreover, vitamin D deficiency is associated with muscle weakness, falls and osteoporotic fractures (Lips and van Schoor 2011). Maintaining a sufficient vitamin D status (>50 nmol/L) is important, not only for bone health, but also because vitamin D deficiency may be associated with various non-skeletal health outcomes (Borradale and Kimlin 2009). Thus, a sufficient vitamin D status may have a disease risk-reduction potential (Grant 2011). Moreover, a U-shaped association exists between 25(OH)D concentrations and risk of cardiovascular disease, certain cancers and overall mortality (Ross et al. 2011).

There is an on-going international discussion regarding which cut-off values should define sufficient 25(OH)D concentrations. There is a general agreement that a 25(OH)D concentration of at least 50 nmol/L is sufficient (Ross et al. 2011; Nordic Council of Ministers 2014). Concurrently, some experts argue that a 25(OH)D concentration >75 nmol/L is required to achieve sufficient vitamin D status and non-skeletal benefits (Holick and Chen 2008; Zhang and Naughton 2010).

It is not easy to determine which doses of vitamin D are required to achieve sufficient 25(OH)D concentrations. The Institute of Medicine (IOM) recently reported that a recommended dietary allowance (RDA) of 15 μg/day for individuals aged 1–70 years will cover the requirement for 97.5 % of the population in the USA and Canada, corresponding to 25(OH)D concentrations of at least 50 nmol/L (Ross et al. 2011). Recently, the recommended intakes (RI) for vitamin D in the Nordic countries were increased from 7.5 to 10 μg/day for individuals aged 2–60 years. This will cover the requirement for 95 % of the Nordic population (Nordic Council of Ministers 2004; Nordic Council of Ministers 2014). Both IOM and Nordic nutrition recommendations (NNRs) 2012 based their RDA and RI on the relationship between 25(OH)D concentrations and bone health.

It is a public health concern that vitamin D intakes in most populations are lower than the RDA or RI (Andersen et al. 2005; Madsen et al. 2013; Nordic Council of Ministers 2014). Food fortification is an effective way to increase vitamin D intake in the general population (O’Mahony et al. 2011), thus ensuring that the general vitamin D intake aligns with the recommendations. During wintertime, a dietary intake of 10 μg/day is needed to maintain 25(OH)D concentrations around 50 nmol/L for the majority of the population in the Nordic countries. For people with little or no sun-exposure, an intake of 20 μg/day of vitamin D is recommended (Nordic Council of Ministers 2014). In Denmark, the mean dietary vitamin D intake is between 2.0 and 2.9 μg/day and does not meet the recommendations for the majority of the population (Tetens et al. 2011). Thus, during wintertime in Denmark, 50–90 % of the population will develop deficient vitamin D status between 30 and 50 nmol/L (Andersen et al. 2005; Thuesen et al. 2012; Madsen et al. 2013).

The main objective of this study was to assess the effect of real-life use of vitamin D3-fortified bread and milk on 25(OH)D concentrations in relation to common genetic variants in CYP2R1 (rs10741657 and rs10766197) and GC (rs4588 and rs842999) in ethnic Danish families with dependent children during a 6-month winter period and furthermore to assess whether vitamin D supplementation will increase 25(OH)D concentration in those with genetically determined low 25(OH)D concentrations. A secondary objective was to evaluate the amount of vitamin D needed to maintain a sufficient 25(OH)D concentrations >50 nmol/L.

Participants and methods

Study design

The present study used data from the VitmaD intervention conducted in Gladsaxe Municipality in Denmark (latitude 56°N). The study design and methods are described in detail elsewhere (Madsen et al. 2013). Briefly, a double-blinded, randomized placebo-controlled intervention trial with apparently healthy ethnically Danish children and adults recruited as families was randomly allocated to either vitamin D3-fortified bread and milk or non-fortified placebo bread and milk during a 6-month winter period (September 2010 to April 2011) without sunlight exposure. The aim of the study design was to investigate a realistic D3-fortification strategy in real-life settings. Participants were instructed to replace their usual consumption of bread and milk with the products provided and in all other aspects, to live a normal life without changing any habits. The study was conducted according to the guideline in the Declaration of Helsinki, and the protocol was approved by the Danish ethics committee (H-4-2010-020) and registered in ClinicalTrials.gov (NCT01184716).

Study population

A total of 201 Danish families with dependent children (n = 782), 4–60 years of age, randomly drawn from the Danish Civil Registration System, participated in the study. Inclusion criteria were age between 4 and 60 years and a permanent address in the Gladsaxe Municipality in Denmark. Exclusion criteria were pregnancy, disease or medication influencing vitamin D metabolism, including dietary supplements with >10 or >5 µg vitamin D/day for children or adults, respectively. All the adult participants and guardians of the children gave written informed consent.

Vitamin D intakes

The participants’ vitamin D intakes were obtained from a self-administered web-based questionnaire based on a semi-quantitative food frequency questionnaire (Andersen et al. 2005) at baseline and at the end of the study. Dietary vitamin D intake was calculated based on the self-reported consumption frequencies and dietary contents of vitamin D (National Food Institute, Technical University of Denmark). Vitamin D intake from dietary supplements was calculated as self-reported frequency of use multiplied with the self-reported vitamin D content of the supplements. The contribution of vitamin D from intakes of vitamin D3-fortified bread and milk was calculated based on the self-reported consumption frequencies, amount and the measured vitamin D contents in the fortified products (5.2 ± 0.3 μg/100 g in wheat bread, 4.3 ± 0.3 μg/100 g in rye bread and 0.38 μg/100 mL in milk) (Madsen et al. 2013). The fortification strategy was to increase vitamin D intake to 7.5 μg/day as recommended in the Nordic nutrition recommendations (NNRs) until September 2013 (Nordic Council of Ministers 2004). Total vitamin D intake was estimated as the sum of dietary vitamin D, usage of multivitamin and vitamin D supplementation and furthermore intake of vitamin D3-fortified bread and milk for the fortification group.

Biochemical analyses

Non-fasting venous blood samples were drawn, and serum and plasma were stored at −80 °C until analysis at Clinical Biochemical Department, Holbæk Hospital, Denmark. Measurements of serum 25(OH)D concentrations relied on the determination of both 25(OH)D2 and 25(OH)D3 and were conducted by isotope dilution liquid chromatography tandem mass spectrometry (LC–MS/MS). As primary calibrator, the standard reference material, vitamin D, in humans (SRM 972) from the National Institute of Standards and Technology was used. The analytic quality of 25(OH)D assay was assured by Vitamin D External Quality Assessment Scheme certification, and the mean bias was −3.2 %. The inter-assay CVs for 25(OH)D2 were 7.6 and 4.6 % at 43 and 150 nmol/L, respectively, and for 25(OH)D3 2.2 and 2.8 % at 30 and 180 nmol/L, respectively, (Madsen et al. 2013). In Denmark, 25(OH)D concentrations <25 nmol/L are defined as vitamin D deficient, between 25 and 50 nmol/L as vitamin D insufficient and >50 nmol/L as vitamin D sufficient for the majority of the population (National Board of Health 2010). 25(OH)D concentrations can be divided by 2.496 to convert from nmol/L to ng/ml.

Plasma PTH levels (CV: 3.4 %) was measured by using immunology analyser Cobas e601 (Roche Diagnostics), and total calcium (CV 3.4 %) was measured by using a chemistry analyser Cobas c501 (Roche Diagnostics).

SNP selection and genotyping

In a previous study (Nissen et al. 2014), we genotyped 25 SNPs in seven vitamin D-related genes (CYP2R1, CYP24A1, CYP27B1, C10orf88, DHCR7/NADSYN1, GC and VDR) selected based on the reports from two GWAS and several candidate gene studies. We found a strong association between common SNPs in CYP2R1 and GC genes and baseline 25(OH)D concentrations in the presently studied 201 healthy Danish families with dependent children. We found that four SNPs, rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC, predicted baseline 25(OH)D concentrations. None of the four SNPs were in LD with each other: rs10741657 and rs10766196 (Pearson’s r = 0.60), rs10741567 and rs842999 (Pearson’s r = 0.03), rs10741657 and rs4588 (Pearson’s r = 0.10), rs10766197 and rs842999 (Pearson’s r = 0.09), rs10766197 and rs4588 (Pearson’s r = 0.05) and rs842999 and rs4588 (Pearson’s r = 0.0.31) were in LD. For the tri-allelic rs842999, there was a dose-dependent relationship between 25(OH)D concentrations and carriers of none, one or two copies of the G-allele and genotypes are presented as GG, GX and XX, where X represents C- or A-alleles.

DNA was purified from buffy coats as described by Miller et al. (1988). SNPs were genotyped using a Sequenom® platform (San Diego, California) and the iPLEX Gold reaction. The SNPs and the primers used are listed in Supplementary Table 1. Each PCR reaction contained 10 ng genomic DNA, 0.5 U HotStart Taq (Qiagen), 1.25 × Enzyme Buffer (Qiagen), 3.5 mM MgCl2, 1 mM of each deoxynucleotide. The primers were added to a final concentration of 500 nM each. The PCRs were performed at the following cycling parameters: 15 min preheat to 94 °C, 45 cycles (20 s 94 °C, 30 s 56 °C, 1 min 72 °C) followed by 3 min 72 °C and stored at −20 °C. The PCR products were treated with shrimp alkaline phosphatase, dephosphorylate unincorporated dNTPs and extension with molecular weight-modified nucleotides were performed in concordance to the manufacturer’s recommendations. The PCRs were cleaned with resin and dispend on SpectroCHIP® bioarrays. The SpectroCHIP® bioarrays were placed in a MALDI-TOF mass spectrometer, and the results were analysed by MassARRAY Type 4.0 (Sequenom) (Nissen et al. 2014).

Of the 782 recruited children and adults, DNA was obtained from 769 participants (98.3 %). A total of 762 (99.1 %) were successfully genotyped. For quality control, 344 duplicated samples (44 %) were randomly placed throughout each of the 384-well plates and the reproducibility was 100 %. No deviation from Hardy–Weinberg equilibrium was observed for the adult population (χ 2 testing, p > 0.05).

Statistical analysis

All statistical analyses were carried out using SAS Enterprise Guide 4.3 (SAS Institute, Inc., Cary. USA). Linear mixed models with family as a random factor were applied in all analyses to account for the non-independency of the participants. Before analysis, 25(OH)D concentrations and PTH levels were log-transformed to approximate a normal distribution and all means are presented as geometric means, unless otherwise specified. A nominal p value of 0.05 was considered statistically significant.

The following categorical variables were used: age (4–11, 12–17, 18–40, 41–60 years), sex (male, female), BMI (underweight, normal weight, overweight, obese) according to standards for children (Cole et al. 2000) and the WHO International standards for adults (World Health Organization 2000) measured at baseline, went on ski and sun vacation during the study period (yes, no), solarium use at least once a week (yes, no) and total calcium at baseline and at the end of the study. The continuous variables are log 25(OH)D concentrations and log PTH levels at baseline and at the end of the study, total vitamin D intake from diet, multivitamins and vitamin D supplements (μg/day).

A genetic risk score (GRS) was calculated as the sum of number of risk alleles. The GRS (range 0–8) was calculated as the sum of number of G-alleles of rs10741657, A-alleles of rs10766197, A-alleles of rs4588 and C/A-alleles of rs842999. A linear mixed model, adjusted for family and confounding variables, was fitted to the log 25(OH)D concentration with GRS as an explanatory factor. The adjusted mean concentration of 25(OH)D was calculated for each GRS. All the analyses were performed for control and fortification group and separately for adults and children.

Furthermore, each GRS category was stratified by quintile of total vitamin D intake (Q1: 0–2.9 μg/day; Q2: 3–7.4 μg/day; Q3: 7.5–9.9 μg/day; Q4: 10.0–14.9 μg/day; and Q5: >15.0 μg/day). Total vitamin D intake was estimated as the sum of dietary vitamin D, use of multivitamin and vitamin D supplementation and, for the fortification group, intake of vitamin D3-fortified bread and milk. The final concentration of 25(OH)D was estimated for each GRS by intake groups adjusted for family and confounding variables.

The prevalence (%) of participants with sufficient (>50 nmol/L) 25(OH)D concentrations was estimated for each GRS by intake groups adjusted for family and confounding variables.

Results

Of the 782 recruited children and adults, 762 participants had complete questionnaire data, genotypes and 25(OH)D concentrations measured at baseline. At the end of the study, a total of 756 participants (control group n = 384 and fortification group n = 384) had complete questionnaire data, genotypes and 25(OH)D concentrations measured. Characteristics of the study population are listed in Table 1, as previously described in detail elsewhere (Madsen et al. 2013; Nissen et al. 2014). At baseline, participants in the control group had significantly higher total calcium levels (p = 0.0438) compared to participants in the fortification group, as previously reported (Madsen et al. 2013). Furthermore, there was a statistically significant difference between the control and fortification group for the use of solarium (p = 0.0059), and ski and sun vacation (p = 0.0006) during the intervention period as previously reported (Madsen et al. 2013).

Table 1 Basic characteristics of the study population (n = 762)

In a previous study (Nissen et al. 2014), we found that at baseline, CYP2R1 (rs10741657 and rs10766197) and GC (rs4588 and rs842999) were strongly associated with 25(OH)D concentrations among all participants (Table 2). At the end of the study, no associations between SNPs rs10741657 and rs10766197 in CYP2R1 or rs4588 and rs842999 in GC and 25(OH)D concentrations were found for the control group. For the fortification group, rs10741657 in CYP2R1 and rs4588 and rs842999 in GC were statistically significantly associated with 25(OH)D concentrations. The association with CYP2R1 (rs10766197) was borderline significant (p = 0.0599). At the end of the study, total vitamin D intake (p < 0.0001) and 25(OH)D concentrations (p < 0.0001) were, as expected, significantly higher in the fortified group compared to the control group as previously reported (Madsen et al. 2013).

Table 2 The association between SNPs in CYP2R1 and GC gene, and 25(OH)D concentrations at baseline (all) and at the end of the study (control/fortification)

There was no difference in PTH levels when stratified by rs10741657, rs10766197, rs4588 and rs842999 for all the participants at baseline (Fig. 1). As anticipated, PTH levels were significantly higher in the control group compared to the fortification group (p = 0.0199) at the end of the study (Table 1). Furthermore, there was a significant difference in PTH levels for rs4588 in both the fortification group (p = 0.0064) and control group (p = 0.0132) at the end of the study. Carriers of the rs4588-AA genotype had significantly lower PTH levels compared to carriers of either the rs4588-CA or rs4588-CC genotype.

Fig. 1
figure 1

Association of rs10741657, rs10766197, rs4588 and rs842999 with PTH levels at baseline for all the participants and stratified by fortification and group control at the end of the study. Results are presented as unadjusted and adjusted geometric means. At baseline, the following variables were adjusted for age, sex, BMI, vacation and baseline total calcium, and at end of the study, the following variables were adjusted for age, sex, BMI, vacation, baseline 25(OH)D concentration, baseline PTH levels and end total calcium. Adjusted p values are given for each genotype. The numbers in the columns present the total numbers of participants carrying this genotype. Error bars indicate 95 % confidence interval. A statistically significant difference in PTH levels was observed for rs4588 in both the fortification and control group at the end of the study

The prevalence of participants with 25(OH)D concentration <30 nmol/L and <50 nmol/L was estimated for each genotype of rs10741657, rs10766197, rs4588 and rs842999 for all the participants at baseline and separately for the control and fortification group at the end of the study (Fig. 2a, b).

Fig. 2
figure 2

The prevalence (%) of <30 nmol/L a and <50 nmol/L b 25(OH)D concentrations in carriers of different genotypes of rs10741657, rs10766197, rs4588 and rs842999 at baseline for all the participants and at the end of the study stratified by control and fortification group. Cut-off value of 25(OH)D <50 nmol/L defines the requirement for optimal bone health for the majority of the population, and cut-off value <30 nmol/L defines the 25(OH)D concentration at which adverse effects on bone health may be expected (Ross et al. 2011)

At baseline, there was no difference in the prevalence of participants having 25(OH)D concentrations <30 nmol/L stratifying by genotype rs10741657, rs10766197, rs4588 and rs842999 (p = 0.2269, 0.1715, 0.6953 and 0.5111), respectively. In contrast, there was significant difference in the prevalence of participants having 25(OH)D concentrations <50 nmol/L for rs10741657, rs4588 and rs842999 (p = 0.0004, <0.0001 and 0.0435), respectively, and rs10766197 was borderline significantly associated (p = 0.0743).

At the end of the study, for the fortification group, a significant difference in the prevalence of participants having 25(OH)D concentrations was found for rs4588 (p = 0.0023 < 30 nmol and for <50 nmol/L p = 0.0002) and rs842999 (p = 0.0029 < 50 nmol/L). No difference in prevalence was observed for rs10741657 and rs10766197 (p = 0.5830 and 0.2348 for <30 nmol/L and for <50 nmol/L p = 0.5466 and 0.6652), respectively. Furthermore, no difference in prevalence was found for rs842999 (p = 0.1194 for <30 nmol/L).

For the control group, only rs842999 <30 nmol/L was significant (p = 0.0455). No significant difference was observed for rs10741657, rs10766197 and rs4588 (p = 0.8694, 0.6130 and 0.2651 < 30 nmol/L and p = 0.5645, 0.4948 and 0.2641 for <50 nmol/L), respectively. Furthermore, rs842999 <50 nmol/L was also found to be non-significant (p = 0.3402). In general, the lowest prevalence of vitamin D deficiency <30 and <50 nmol/L was observed at baseline (p = 0.0001 and 0.0001), respectively. Participants in the control group presented more often with vitamin D deficiency <30 and <50 nmol/L compared to the fortification group (p = 0.0001 for <30 nmol/L and for <50 nmol/L p = 0.0001),

At the end of the study, to determine the combined contributions of rs10741657, rs10766197, rs4588 and rs842999, a GRS was calculated individually for the control and fortification group and separately for all, adults and children (Fig. 3a–c). Participants carrying seven or eight (7–8) risk alleles were combined due to small sample size. The coefficients for rs10741657, rs10766197, rs4588 and rs842999 were very similar in a mixed regression model including all SNPs, and therefore, it was not necessary to weight the different risk alleles by the correlation coefficient. A linear mixed model with family as a random factor, adjusted for age, sex, BMI, total vitamin D intake, and ski and sun vacation showed that for the control group, there was no difference in 25(OH)D concentrations for carriers of 0 to 7–8 risk alleles (p = 0.1428, 0.2881 and 0.7667) for all, adults and children, respectively. For the fortification group, there was a negative linear trend between 25(OH)D concentrations and carriers of 0 to 7–8 risk alleles for all, adults and children (p < 0.0001, 0.0025 and 0.0023, respectively). Overall, there was a mean difference in 25(OH)D concentrations of 28.2, 28.6 and 31.9 nmol/L between carriers of no risk alleles and carriers of all 7–8 risk alleles in all, adults and children, respectively. Overall, the same GRS pattern was observed for adults and children.

Fig. 3
figure 3

Estimated mean 25(OH)D concentrations at the end of the study for each genetic risk score category stratified by control and fortification group, separately for all (a), adults (b) and children (c). Individuals carrying 7 or 8 (7–8) risk alleles were combined due to small sample size. Genetic risk score (range 0 to 7–8) was calculated as the sum of number of G-alleles of rs10741657, A-alleles of rs10766197, A-alleles of rs4588 and C/A-alleles of rs842999. The numbers in the columns present the total numbers of participants carrying the risk score. Error bars indicate 95 % confidence interval

We estimated the effect of total vitamin D intake for each category GRS (range 0–8), for the combined contributions of rs10741657, rs10766197, rs4588 and rs842999 (Fig. 4). Each participant was stratified by quintile of total vitamin D intake. Total vitamin D intake was estimated as the sum of dietary vitamin D, use of multivitamin and vitamin D supplements and, for the fortification group, self-reported intake of vitamin D3-fortified bread and milk. Quintile stratification for total vitamin D intake was based on different RDA or RI: <3 μg/day (no supplementation), <7.5 μg/day (old NNRs 2004), <10 μg/day (present NNRs 2012), <15 μg/day (IOM) or >15 μg/day. The following quintile stratification cut-off values were used: quintile 1: 0–2.9 μg/day; quintile 2: 3–7.4 μg/day; quintile 3: 7.5–9.9 μg/day; quintile 4: 10.0–14.9 μg/day and quintile 5: >15.0 μg/day. The control and fortification groups were combined in the linear mixed model. Individuals carrying 0, 1 or 2 (0–2) risk alleles or individuals carrying 6, 7 or 8 (6–8) risk alleles were combined due to small sample sizes after quintile stratification by total vitamin D intake. A total of 25.1, 22.4, 23.4, 15.6 and 13.6 % of the adult participants carried 0–2, 3, 4, 5 or 6–8 risk alleles, respectively. The majority of the participants in the control group had low total vitamin D intake and were therefore primarily located in the first two quintiles. In general, there was a statistically significant, positive linear relationship between total vitamin D intake and 25(OH)D concentrations among carriers of 0–2, 3, 4 or 5 risk alleles, (p = 0.0012, 0.0001, 0.0118 and 0.0029, respectively). For individuals carrying 6–8 risk alleles, there was no statistically significant relationship between total vitamin D intake and 25(OH)D concentrations (p = 0.1051).

Fig. 4
figure 4

Mean 25(OH)D concentrations at the end of the study for each genetic risk score category stratified by total vitamin D intakes for the study population. Total vitamin D intake was estimated as the sum of dietary vitamin D, usage of multivitamin and vitamin D supplementation and, for the fortification group, intake of vitamin D3-fortified bread and milk. The following quintile stratification was used: quintile 1: 0–2.9 μg/day; quintile 2: 3–7.4 μg/day; quintile 3: 7.5–9.9 μg/day; quintile 4: 10.0–14.9 μg/day; and quintile 5: >15.0 μg/day. Genetic risk score (range 0–8) was calculated as the sum of number of G-alleles of rs10741657, A-alleles of rs10766197, A-alleles of rs4588 and C/A-alleles of rs842999. Individuals carrying 0, 1 or 2 (0–2) risk alleles and individuals carrying 6, 7 or 8 (6–8) risk alleles were combined due to small sample size after quintile stratification by total vitamin D intake. The numbers in the columns present the total numbers of participants carrying this risk score. Error bars indicate 95 % confidence interval

A the end of the winter season in Denmark, a total vitamin D intake of <3 μg/day was not sufficient for 95 % of the study population to achieve sufficient (>50 nmol/L) 25(OH)D concentrations, regardless of the number of risk alleles they carried (Fig. 4). For participants carrying 0–2 or 3 risk alleles, a total daily vitamin D intake between 3 and 7.4 μg seemed to be sufficient for 95 % of the study population to achieve sufficient 25(OH)D concentrations. For participants carrying four risk alleles, a total daily vitamin D intake >7.5 μg seemed to be sufficient for 95 % of the study population to achieve sufficient 25(OH)D concentrations. For participants carrying five risk alleles, a total daily vitamin D intake >10 μg seemed to be sufficient for 95 % of the study population to achieve sufficient 25(OH)D concentrations. For participants carrying 6–8 risk alleles, a total daily vitamin D intake >15 μg was almost enough for 95 % of the study population to achieve sufficient 25(OH)D concentrations.

In addition, we determined the percentage of participants with sufficient 25(OH)D concentrations (Fig. 5). Sufficient 25(OH)D concentrations were achieved for all participants carrying 0–2, 3 or 4 risk alleles and who consumed >15 μg/day of vitamin D. For participants carrying 5 or 6–8 risk alleles, this fell to 86 and 90 %, respectively. Furthermore, sufficient 25(OH)D concentrations were achieved for 87, 90, 83, 84 and 67 % of the participants carrying 0–2, 3, 4, 5 or 6–8 risk alleles and who consumed 10–14.9 μg/day. This fell to 80, 76, 86, 50 and 53 % and 57, 50, 61, 52 and 41 % for participants carrying 0–2, 3, 4, 5 or 6–8 risk alleles and who consumed 7.5–9.9 μg/day or 3.0–7.4 μg/day of vitamin D, respectively.

Fig. 5
figure 5

The prevalence (%) of sufficient 25(OH)D concentrations, defined as >50 nmol/L, for each genetic risk score category stratified by quintile of total vitamin D intake at the end of the study. Total vitamin D intake was estimated as the sum of dietary vitamin D, use of multivitamin and vitamin D supplements and, for the fortification group, intake of vitamin D3-fortified bread and milk. The following quintile stratification was used: quintile 1: 0–2.9 μg/day; quintile 2: 3–7.4 μg/day; quintile 3: 7.5–9.9 μg/day; quintile 4: 10.0–14.9 μg/day; and quintile 5: >15.0 μg/day. Genetic risk score (range 0–8) was calculated as the sum of number of G-alleles of rs10741657, A-alleles of rs10766197, A-alleles of rs4588 and C/A-alleles of rs842999. Individuals carrying 0, 1 or 2 (0–2) risk alleles and individuals carrying 6, 7 or 8 (6–8) risk alleles were combined due to small sample size after quintile stratification by total vitamin D intake

Discussion

In the present study, we show that genetic variation influences 25(OH)D concentrations considerably. Genetically predisposed individuals carrying 6–8 risk alleles of rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC need >15 μg/day or more vitamin D to reach 25(OH)D concentrations >50 nmol/L during winter. Furthermore, there was a statistically significant dose-dependent relationship between 25(OH)D concentration and total vitamin D intake for carriers of 0–5 risk alleles of SNPs rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC. A dose-dependent relationship was also observed for carriers of 6–8 risk alleles, but the increase in 25(OH)D concentrations was not statistically significant.

At baseline, our study showed that there was statistically significant difference in the prevalence of participants presenting with 25(OH)D concentration <50 nmol/L for rs10741657, rs10766197, rs4588 and rs842999. The significant differences in prevalence disappeared during the winter for the control group, but were maintained for rs4588 and rs842999 in the fortification group. For the fortification group, the highest prevalence of 25(OH)D <50 nmol/L was observed for the rs4588-AA genotype. In contrast, in the control group, rs4588-AA carriers had the lowest prevalence of 25(OH)D <50 nmol/L. This indicates that although carriers of the rs4588-AA genotype in the fortification group were more prone to be vitamin D deficient, rs4588-AA carriers in the control group were less prone to be vitamin D deficient. This may indicate that rs4588-AA carriers have a somewhat low but very stable 25(OH)D concentrations. Paradoxically, a recessive effect was observed for rs4588-AA carriers on PTH levels in both the fortification and control group at the end of the study. Participants with the rs4588-AA genotype have the lowest PTH levels and 25(OH)D concentrations compared to rs4588-CC or rs4588-CA carriers. Similar to our findings, Pekkinen et al. (2014) found a dose–response effect of rs4588 on PTH concentrations in 231 Finnish children and adolescents aged 7–19 years, with rs4588-AA carriers having the lowest PTH and 25(OH)D concentrations. Further studies are warranted to investigate the underlying biological mechanism of this observation.

At the end of the study, there was a pronounced positive effect of real-life usage of vitamin D3-fortified bread and milk on 25(OH)D concentrations. For the fortification group, 25(OH)D concentrations were significantly associated with rs10741657 in CYP2R1, and with rs4588 and rs842999 in GC. Furthermore, rs10766197 in CYP2R1 was borderline significantly associated with 25(OH)D concentrations. These winter results resemble the results found at baseline (late summer) and indicate that when vitamin D is received primarily as vitamin D3-food fortification during the winter, the association between 25(OH)D concentrations and genetic variation observed at late summer for rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC is maintained. In contrast, the baseline association between 25(OH)D concentrations and rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC disappeared during the winter for the control group. Our findings are consistent with the findings from two previous studies (Gozdzik et al. 2011; Engelman et al. 2013). Gozdzik et al. (2011) found that rs4588 in GC was associated with 25(OH)D concentrations in Canadians of European descent during the fall (p = 0.009), but not during the winter (p = 0.535). Similarly, Engelman et al. (Engelman et al. 2013) found two SNPs in GC (rs4588 and rs7041) and four SNPs in CYP2R1 (rs105000804, rs11023380, 2060763, 11023374) to be strongly associated with 25(OH)D concentrations in individuals whose blood was drawn in summer but not in individuals whose blood was drawn in winter month.

Engelman et al. (2013) performed a GRS for rs4588 in GC and rs2060793 in CYP2R1. The risk scores were highly significantly associated with 25(OH)D concentrations in individuals with high external source of vitamin D (>10 μg/day) but not in individuals with low external source of vitamin D (<10 μg/day). In addition, Gozdzik et al. (2011) found that vitamin D intake was significantly predictive of 25(OH)D concentrations in individuals carrying the rs4588 (T436 K) or in GC diplotypes during fall and winter. Our results support these findings by Engelman et al. (2013) and Gozdzik et al. (2011). We performed a GRS including the four SNPs rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC. For the fortification group, the GRS was highly significantly associated with 25(OH)D concentrations (p < 0.0001) but not for the control group (p = 0.1428) during winter. In general, children had higher mean 25(OH)D concentrations compared to adults. For the fortification group, an explanation could be that the children consumed more vitamin D3-fortified bread and milk compared to the adults. Approximately 90 % of the total intake of consumed bread and milk was the products provided by the study, with no difference in compliance between children and adults (Madsen et al. 2013). In general, the children were more often multivitamin users compared to the adults (Madsen et al. 2013).

When stratifying total vitamin D intake into quintiles, our data suggest that it is difficult to raise 25(OH)D concentrations to a sufficient level in participants carrying 6–8 risk alleles with vitamin D3-fortified bread and milk. A statistically non-significant increase in 25(OH)D concentrations was found comparing the lowest and highest quintile of vitamin D intake for participants carrying 6–8, but with a much lower rate (+Δ17.6 nmol/L) compared to participants carrying 0–2, 3, 4 or 5 risk alleles (+Δ28.8, 36.5, 24.2 and 33.6 nmol/L), respectively, (Fig. 4). Whether this also applies for vitamin D synthesized in the skin during UVB exposure remains to be further investigated. These increases are similar to the findings by Engelman et al. (2013). They found that among individuals carrying 3–4 risk alleles of GC (rs4588) and CYP2R1 (rs2060793), the lowest increase in 25(OH)D concentrations was observed in individuals carrying 3–4 risk alleles (+Δ16.7 nmol/L) compared to individuals with fewer risk alleles (+Δ27.7 nmol/L).

In our study population, 67 % of the participants carrying 6–8 risk alleles had sufficient 25(OH)D concentrations in contrast to 87, 90, 83 and 84 % for participants carrying 0–2, 3, 4 or 5 risk alleles, respectively, when following IOMs RDA of 15 μg/day for individuals aged 1–70 years. Following the Nordic countries RI of 10 μg/day for individual aged 2–60 years, only 50 and 53 % of the participants carrying 5 or 6–8 risk alleles, respectively, had sufficient 25(OH)D concentrations compared to 80, 76 and 86 % of the participants carrying 0–2, 3 or 4 risk alleles, respectively. This indicates that genetic predisposition may have a large impact on 25(OH)D concentrations. Participants having a high GRS may need a higher amount of vitamin D supplementation than participants carrying a lower GRS in order to reach sufficient 25(OH)D concentrations. We provide evidence that participants with different genetic profiles need different amounts of vitamin D supplementation to achieve sufficient 25(OH)D concentrations. Epidemiological studies have found association between blood levels of vitamin D concentrations and risk of cancer, but the significance of genetically determined low vitamin D concentration is not clear.

In agreement with our findings, Cranney et al. (2007) concluded that vitamin D3-doses of 10–20 μg/day may be insufficient to prevent vitamin D deficiency in at-risk individuals. Cashman et al. (2011) concluded that for a population to achieve 25(OH)D concentrations of 50 nmol/L, an average intake of 9 μg/day vitamin D was needed. Nevertheless, taking inter-individual variation into account 23.5 μg/day of vitamin D3 was needed for 95 % of the population to reach a 25(OH)D concentration of 50 nmol/L. Engelman et al. (2013) found that all of the individuals with no risk alleles of rs4588 and rs2060793 who consumed at least 17 μg/day (670 IU/day) had 25(OH)D >50 nmol/L. This fell to 84, 72 and 62 %, respectively, for individuals carrying 1, 2 or 3–4 risk alleles who also consumed at least 17 μg/day.

Our study has several strengths in that we ensured a large age span (4–60 years), had both genders represented, and both children and adults were included due to the family-based design (Madsen et al. 2013). 25(OH)D concentrations were measured by a specific analytical method (LC–MS/MS). We took into account that non-genetic factors such as vitamin D intake and season are known to influence 25(OH)D concentrations. We estimated total vitamin D intake, and blood samples were drawn during the same seasons for all the participants. A disadvantage is that some of the known predictors of 25(OH)D concentration were quantified by self-reported questionnaire data.

In summary, we found that after consuming vitamin D3-fortified bread and milk during a winter season, the effect of genetic variation in the CYP2R1 and GC genes on 25(OH)D concentrations resembles the results found in late summer. The association with genetic variation observed for CYP2R1 and GC genes in late summer disappeared during the winter season for the control group. We found that carriers of the rs4588-AA genotype had the highest prevalence of 25(OH)D concentration <50 nmol/L at baseline and at the end of the study for the fortification group. In contrast, rs4588-AA carriers in the control group had the lowest prevalence. It seems like rs4588-AA carriers have a low but very stable 25(OH)D concentration, and interestingly, also low PTH level.

In this study, we demonstrated that carriers of a high GRS of CYP2R1 (rs10741657 and rs10766197) and GC (rs4588 and rs842999) are more prone to be vitamin D deficient compared to carriers of a low GRS. Furthermore, carriers of a high GRS may need a higher amount of vitamin D3 supplementation to achieve sufficient 25(OH)D concentrations. Importantly, for public health recommendations, it seems that with increasing vitamin D intake, genetically determined low risk carriers with sufficient 25(OH)D concentrations achieve even higher 25(OH)D concentrations with the used real-life vitamin D3-fortification model.

Abbreviations

DBP:

Vitamin D-binding protein

GC:

Vitamin D-binding protein gene or group-specific component

GRS:

Genetic risk score

GWAS:

Genome-wide association studies

IOM:

Institute of Medicine

LC–MS/MS:

Isotope dilution liquid chromatography tandem mass spectrometry

LD:

Linkage disequilibrium

NNRs:

Nordic nutrition recommendations

PTH:

Parathyroid hormone

RDA:

Recommended dietary allowance

RI:

Recommended intakes

25(OH)D:

25-Hydroxyvitamin D

SNPs:

Single-nucleotide polymorphisms

UVB:

Ultraviolet B radiation

References

  • Abbas S, Linseisen J, Slanger T et al (2008) The Gc2 allele of the vitamin D binding protein is associated with a decreased postmenopausal breast cancer risk, independent of the vitamin D status. Cancer Epidemiol Biomarkers Prev 17:1339–1343. doi:10.1158/1055-9965.EPI-08-0162

    Article  CAS  PubMed  Google Scholar 

  • Ahn J, Yu K, Stolzenberg-Solomon R et al (2010) Genome-wide association study of circulating vitamin D levels. Hum Mol Genet 19:2739–2745. doi:10.1093/hmg/ddq155

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Andersen R, Mølgaard C, Skovgaard LT et al (2005) Teenage girls and elderly women living in northern Europe have low winter vitamin D status. Eur J Clin Nutr 59:533–541. doi:10.1038/sj.ejcn.1602108

    Article  CAS  PubMed  Google Scholar 

  • Arnaud J, Constans J (1993) Affinity differences for vitamin D metabolites associated with the genetic isoforms of the human serum carrier protein (DBP). Hum Genet 92:183–188

    Article  CAS  PubMed  Google Scholar 

  • Bikle D, Gee E, Halloran B et al (1986) Assessment of the free fraction of 25-hydroxyvitamin D in serum and its regulation by albumin and the vitamin D-binding protein. J Clin Endocrinol Metab 63:954–959

    Article  CAS  PubMed  Google Scholar 

  • Bogh MK, Schmedes A, Philipsen P et al (2010) Vitamin D production after UVB exposure depends on baseline vitamin D and total cholesterol but not on skin pigmentation. J Invest Dermatol 130:546–553. doi:10.1038/jid.2009.323

    Article  CAS  PubMed  Google Scholar 

  • Borradale D, Kimlin M (2009) Vitamin D in health and disease: an insight into traditional functions and new roles for the “sunshine vitamin”. Nutr Res Rev 22:118–136. doi:10.1017/S0954422409990102

    Article  CAS  PubMed  Google Scholar 

  • Bu FX, Armas L, Lappe J et al (2010) Comprehensive association analysis of nine candidate genes with serum 25-hydroxy vitamin D levels among healthy Caucasian subjects. Hum Genet 128:549–556. doi:10.1007/s00439-010-0881-9

    Article  CAS  PubMed  Google Scholar 

  • Cashman KD, Fitzgerald AP, Kiely M, Seamans KM (2011) A systematic review and meta-regression analysis of the vitamin D intake-serum 25-hydroxyvitamin D relationship to inform European recommendations. Br J Nutr 106:1638–1648. doi:10.1017/S0007114511005058

    Article  CAS  PubMed  Google Scholar 

  • Cheng JB, Levine M, Bell N et al (2004) Genetic evidence that the human CYP2R1 enzyme is a key vitamin D 25-hydroxylase. Proc Natl Acad Sci USA 101:7711–7715. doi:10.1073/pnas.0402490101

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Cole TJ, Bellizzi MC, Flegal KM, Dietz WH (2000) Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ 320:1240–1243

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Cranney A, Horsley T, O’Donnell S et al (2007) Effectiveness and safety of vitamin D in relation to bone health. Evid Rep Technol Assess (Full Rep) 158:1–235

  • Engelman CD, Fingerlin TE, Langefeld CD et al (2008) Genetic and environmental determinants of 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D levels in Hispanic and African Americans. J Clin Endocrinol Metab 93:3381–3388. doi:10.1210/jc.2007-2702

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Engelman CD, Meyers KJ, Iyengar SK et al (2013) Vitamin D intake and season modify the effects of the GC and CYP2R1 genes on 25-hydroxyvitamin D concentrations. 25:17–26. doi:10.3945/jn.112.169482.17

  • Fu L, Yun F, Oczak M et al (2009) Common genetic variants of the vitamin D binding protein (DBP) predict differences in response of serum 25-hydroxyvitamin D [25(OH)D] to vitamin D supplementation. Clin Biochem 42:1174–1177. doi:10.1016/j.clinbiochem.2009.03.008

    Article  CAS  PubMed  Google Scholar 

  • Gozdzik A, Zhu J, Wong BY et al (2011) Association of vitamin D binding protein (VDBP) polymorphisms and serum 25(OH)D concentrations in a sample of young Canadian adults of different ancestry. J Steroid Biochem Mol Biol 127:405–412. doi:10.1016/j.jsbmb.2011.05.009

    Article  CAS  PubMed  Google Scholar 

  • Grant WB (2011) An estimate of the global reduction in mortality rates through doubling vitamin D levels. Eur J Clin Nutr 65:1016–1026. doi:10.1038/ejcn.2011.68

    Article  CAS  PubMed  Google Scholar 

  • Holick M (2007) Vitamin D deficiency. N Engl J Med 357:266–281. doi:10.1056/NEJMra070553

    Article  CAS  PubMed  Google Scholar 

  • Holick M, Chen T (2008) Vitamin D deficiency : a worldwide problem with health consequences. Am J Clin Nutr 87:1080–1086

    Google Scholar 

  • Kamboh MI, Ferrell RE (1986) Ethnic variation in vitamin D-binding protein (GC): a review of isoelectric focusing studies in human populations. Hum Genet 72:281–293

    Article  CAS  PubMed  Google Scholar 

  • Karohl C, Su S, Kumari M et al (2010) Heritability and seasonal variability of vitamin D concentrations in male twins. Am J Clin Nutr 25:1393–1398. doi:10.3945/ajcn.2010.30176.1

    Article  Google Scholar 

  • Kawakami M, Blum C, Ramakrishnan R et al (1981) Turnover of the plasma binding protein for vitamin D and its metabolites in normal human subjects. J Clin Endocrinol Metab 53:1110–1116

  • Kurylowicz A, Ramos-Lopez E, Bednarczuk T, Badenhoop K (2006) Vitamin D-binding protein (DBP) gene polymorphism is associated with Graves’ disease and the vitamin D status in a Polish population study. Exp Clin Endocrinol Diabetes 114:329–335. doi:10.1055/s-2006-924256

    Article  CAS  PubMed  Google Scholar 

  • Lauridsen AL, Vestergaard P, Hermann AP et al (2005) Plasma concentrations of 25-hydroxy-vitamin D and 1,25-dihydroxy-vitamin D are related to the phenotype of Gc (vitamin D-binding protein): a cross-sectional study on 595 early postmenopausal women. Calcif Tissue Int 77:15–22. doi:10.1007/s00223-004-0227-5

    Article  CAS  PubMed  Google Scholar 

  • Lips P, van Schoor NM (2011) The effect of vitamin D on bone and osteoporosis. Best Pract Res Clin Endocrinol Metab 25:585–591. doi:10.1016/j.beem.2011.05.002

    Article  CAS  PubMed  Google Scholar 

  • Lu L, Sheng H, Li H et al (2012) Associations between common variants in GC and DHCR7/NADSYN1 and vitamin D concentration in Chinese Hans. Hum Genet 131:505–512. doi:10.1007/s00439-011-1099-1

    Article  CAS  PubMed  Google Scholar 

  • Madsen KH, Rasmussen LB, Andersen R et al (2013) Randomized controlled trial of the effects of vitamin D–fortified milk and bread on serum 25-hydroxyvitamin D concentrations in families in Denmark during winter: the VitmaD study 1–3. Am J Clin Nutr 1–9. doi:10.3945/ajcn.113.059469

  • Miller S, Dykes D, Polesky H (1988) A simple salting out procedure for extracting DNA from human nucleated cells. NucleicAcids Res 16:55404

    Google Scholar 

  • Monticielo OA, Teixeira TM, Chies JA et al (2012) Vitamin D and polymorphisms of VDR gene in patients with systemic lupus erythematosus. Clin Rheumatol. doi:10.1007/s10067-012-2021-5

    Google Scholar 

  • National Board of Health (2010) Forebygggelse, diagnistik og behandling af D-vitaminmangel (Prevention, diagnostics and treatment of vitamin D deficiency. In: Natl. Board Heal. http://sundhedsstyrelsen.dk/~/media/FA2FC43A29D146918C9695BEC2716A33.ashx.

  • National Food Institute, Technical University of Denmark. Danish food composition databank. www.foodcomp.dk

  • Nissen J, Rasmussen LB, Ravn-Haren G et al (2014) Common variants in CYP2R1 and GC genes predict vitamin D Concentrations in healthy danish children and adults. PLoS One 9:e89907. doi:10.1371/journal.pone.0089907

    Article  PubMed Central  PubMed  Google Scholar 

  • Nordic Council of Ministers (2004) Nordic Nutrition Recommendations 2004, 4th edn. Copenhagen, Denmark. Norden

  • Nordic Council of Ministers (2014) Nordic Nutrition Recommendations 2012. Part 1. Copenhagen, Denmark. Norden

  • O’Mahony L, Stepien M, Gibney MJ et al (2011) The potential role of vitamin D enhanced foods in improving vitamin D status. Nutrients 3:1023–1041. doi:10.3390/nu3121023

    Article  PubMed Central  PubMed  Google Scholar 

  • Pekkinen M, Saarnio E, Viljakainen HT et al (2014) Vitamin D binding protein genotype is associated with serum 25-hydroxyvitamin d and pth concentrations, as well as bone health in children and adolescents in Finland. PLoS One 9:e87292. doi:10.1371/journal.pone.0087292

    Article  PubMed Central  PubMed  Google Scholar 

  • Ramos-lopez E, Brück P (2007) CYP2R1 (vitamin D 25-hydroxylase) gene is associated with susceptibility to type 1 diabetes and vitamin D levels in Germans. Diabetes Metab Res Rev 1:631–636. doi:10.1002/dmrr

    Article  Google Scholar 

  • Ross A, Manson J, Abrams S et al (2011) The 2011 report on dietary reference intakes for calcium and vitamin D from the Institute of Medicine: what clinicians need to know. J Clin Endocrinol Metab 96:53–58. doi:10.1210/jc.2010-2704

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Shea MK, Benjamin EJ, Dupuis J et al (2009) Genetic and non-genetic correlates of vitamins K and D. Eur J Clin Nutr 63:458–464. doi:10.1038/sj.ejcn.1602959

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Sinotte M, Diorio C, Berube S et al (2009) Genetic polymorphisms of the vitamin D binding protein and plasma concentrations of 25-hydroxyvitamin D in premenopausal women 1–3. Am J Clin Nutr 25:634–640. doi:10.3945/ajcn.2008.26445.INTRODUCTION

    Article  Google Scholar 

  • Speeckaert M, Huang G, Delanghe JR, Taes YEC (2006) Biological and clinical aspects of the vitamin D binding protein (Gc-globulin) and its polymorphism. Clin Chim Acta 372:33–42. doi:10.1016/j.cca.2006.03.011

    Article  CAS  PubMed  Google Scholar 

  • Tetens I, Biltoft-Jensen A, Spagner C et al (2011) Intake of micronutrients among Danish adult users and non-users of dietary supplements. Food Nutr Res 55:1–8. doi:10.3402/fnr.v55i0.7153

    Article  Google Scholar 

  • Thuesen B, Husemoen L, Fenger M et al (2012) Determinants of vitamin D status in a general population of Danish adults. Bone 50:605–610. doi:10.1016/j.bone.2011.12.016

    Article  CAS  PubMed  Google Scholar 

  • Wang TJ, Zhang F, Richards JB et al (2010) Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet 376:180–188. doi:10.1016/S0140-6736(10)60588-0

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Wjst M, Altmüller J, Faus-Kessler T et al (2006) Asthma families show transmission disequilibrium of gene variants in the vitamin D metabolism and signalling pathway. Respir Res 7:60. doi:10.1186/1465-9921-7-60

    Article  PubMed Central  PubMed  Google Scholar 

  • World Health Organization (2000) Obesity: preventing and managing the global epidemic. Report of a WHO consultation. WHO Heal Organ Tech Rep Ser 894:1–253

  • Zhang R, Naughton DP (2010) Vitamin D in health and disease: current perspectives. Nutr J 9:65. doi:10.1186/1475-2891-9-65

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  • Zhang Y, Wang X, Liu Y et al (2012) The GC, CYP2R1 and DHCR7 genes are associated with vitamin D levels in northeastern Han Chinese children. Swiss Med Wkly 142:1–6. doi:10.4414/smw.2012.13636

    Google Scholar 

  • Zhang Z, He JW, Fu WZ et al (2013) An analysis of the association between the vitamin D pathway and serum 25-hydroxyvitamin D levels in a healthy Chinese population. J Bone Miner Res. doi:10.1002/jbmr.1926

    Google Scholar 

Download references

Acknowledgments

The authors would like to acknowledge all the families for their participation. We thank technician Bettina Hansen from Department of Biomedicine, Aarhus University, for performing the genotyping.

Conflict of interest

The authors declare that they have no conflict of interest.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Janna Nissen.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 41 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nissen, J., Vogel, U., Ravn-Haren, G. et al. Real-life use of vitamin D3-fortified bread and milk during a winter season: the effects of CYP2R1 and GC genes on 25-hydroxyvitamin D concentrations in Danish families, the VitmaD study. Genes Nutr 9, 413 (2014). https://0-doi-org.brum.beds.ac.uk/10.1007/s12263-014-0413-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1007/s12263-014-0413-7

Keywords