Skip to main content


SNPs related to vitamin D and breast cancer risk: a case-control study



It has been suggested that vitamin D might protect from breast cancer, although studies on levels of vitamin D in association with breast cancer have been inconsistent. Genome-wide association studies (GWASs) have identified several single-nucleotide polymorphisms (SNPs) to be associated with vitamin D. The aim of this study was to investigate such vitamin D-SNP associations in relation to subsequent breast cancer risk. A first step included verification of these SNPs as determinants of vitamin D levels.


The Malmö Diet and Cancer Study included 17,035 women in a prospective cohort. Genotyping was performed and was successful in 4058 nonrelated women from this cohort in which 865 were diagnosed with breast cancer. Levels of vitamin D (25-hydroxyvitamin D) were available for 700 of the breast cancer cases and 643 of unaffected control subjects. SNPs previously associated with vitamin D in GWASs were identified. Logistic regression analyses yielding ORs with 95% CIs were performed to investigate selected SNPs in relation to low levels of vitamin D (below median) as well as to the risk of breast cancer.


The majority of SNPs previously associated with levels of vitamin D showed a statistically significant association with circulating vitamin D levels. Heterozygotes of one SNP (rs12239582) were found to have a statistically significant association with a low risk of breast cancer (OR 0.82, 95% CI 0.68–0.99), and minor homozygotes of the same SNP were found to have a tendency towards a low risk of being in the group with low vitamin D levels (OR 0.72, 95% CI 0.52–1.00). Results from stratified analyses showed diverse associations with breast cancer risk for a few of the tested SNPs, depending on whether vitamin D level was high or low.


SNPs associated with vitamin D may also be associated with the risk of breast cancer. Even if such a risk is small, the allele frequency of the SNP variants is high, and therefore the population attributable risk could be substantial. It is also possible that vitamin D levels may interact with genomic traits with regard to breast cancer risk.


About 5–10% of breast cancers are considered hereditary, from which the known breast cancer genes account for 3–4% [1]. Apart from known breast cancer genes, genome-wide association studies (GWASs) have previously identified more than 170 genetic polymorphisms associated with the risk of breast cancer [2]. It has been suggested that single-nucleotide polymorphisms (SNPs) may add up to 14% of heredity of breast cancer [1].

Ecological and epidemiological studies have suggested a beneficial effect of relatively high vitamin D levels, owing to solar exposure, on breast cancer risk and survival [3,4,5,6]. There are several prospective epidemiological studies on the relationship between serum levels of vitamin D and breast cancer incidence. The results have been conflicting [7,8,9,10], although authors of a meta-analysis concluded that there is an inverse relationship between levels of vitamin D and breast cancer risk [11]. Diverse results may be a result of misclassification of vitamin D status, and a better marker of stable and unconfounded vitamin D status over time may be found in the genotype, as studied using a Mendelian randomization approach [12].

Researchers investigating vitamin D-related SNPs and breast cancer risk have focused mainly on SNPs located in the vitamin D receptor (VDR) gene. Results derived from these studies have also been conflicting [13,14,15,16]. SNPs of the VDR have not been found to be associated with vitamin D in previous GWASs [2]. The aim of the present study was to investigate breast cancer risk in relation to SNPs previously identified in GWASs on vitamin D levels and related phenotypes. The present study was a prospective, nested case-control study with information on SNPs, vitamin D levels and subsequent breast cancer.


Malmö diet and cancer study

Between 1991 and 1996, all residents of Malmö, Sweden, born between 1923 and 1950 were invited to participate in a prospective cohort study. Among invited women, 43% participated, resulting in a female cohort of 17,035 women [17]. At inclusion, written informed consent was obtained from all participants. Baseline examinations included anthropometric measurements by a trained nurse who also drew blood samples. Subjects were included evenly over the calendar year, although there was less recruitment in December and June and none in July. Participants also provided information on reproductive factors and lifestyle via a self-administered questionnaire. Information on previous gynaecological surgery was retrieved from medical records, and menopausal status was defined using these data together with information provided in the questionnaire [18]. The Malmö Diet and Cancer Study (MDCS) (LU 51-90) and the present study (Dnr 153/2004 and Dnr 682/2009) were approved by the regional ethics committee in southern Sweden.

Study population

In a previous case-control study, women diagnosed with breast cancer in the MDCS until December 31, 2006 (n = 764), were matched on age and date of inclusion with control subjects (n = 764), based on 1482 individuals (some were included twice because incidence density matching was applied) [8]. One breast cancer case was classified to be without disease in a later follow-up. Until December 31, 2009, an additional 183 women were diagnosed with breast cancer, leading to inclusion of 946 women with breast cancer in the MDCS cohort. The control group was expanded to include an additional 2658 randomly selected women from the MDCS cohort without a breast cancer diagnosis. These women were also part of the MDCS cardiovascular cohort, a randomized subsample of MDCS [19]. Together with the 704 control subjects from the previous case-control study who had not developed breast cancer, potential control subjects added up to 3362. Material was accessible for genotyping from 901 of the breast cancer cases and 3335 of the control subjects (Fig. 1).

Fig. 1

Flowchart of cases and control subjects. SNP Single-nucleotide polymorphism

Genotyping was performed using the HumanOmniExpressExome BeadChip and iScan System (Illumina, San Diego, CA, USA) during 2012–2013 at the Broad Institute of MIT and Harvard University (Cambridge, MA, USA) (n = 4058) and at the Clinical Research Centre, Skåne University Hospital, Malmö, Sweden (n = 178). Sample-based quality control (QC) included controls for call rate, > 95% of SNPs per individual, and control for excess heterozygosity, which led to exclusion of 51 individuals. Owing to first- and second-degree relationships between individuals, a further 127 individuals were excluded, and the relative with the highest call rate was kept in the study population. After QC and exclusion of relatives, 865 patients with breast cancer and 3193 control subjects were available for statistical analysis (Fig. 1). PLINK version 1.07 software was used for QC and exclusion of relatives.

Levels of vitamin D

In the previous case-control study of the MDCS cohort mentioned above, researchers investigated, among other things, serum levels of vitamin D in relation to breast cancer risk [8]. For this study, high-pressure liquid chromatography was used to analyse 25-hydroxyvitamin D [25(OH)D] levels, and laboratory analysis was successful in 700 breast cancer cases and 643 control subjects in the present study population (Fig. 1). Amongst subjects with data on SNPs and vitamin D level (1343 in total), 320 women (165 cases, 155 control subjects) were also subjects in the cardiovascular cohort. There is no consensus regarding adequate levels of vitamin D, and the levels vary substantially over the year; thus it is not possible to clinically define a high or low vitamin D level on a single measurement [20]. In this study, the median level of vitamin D of each calendar month of sampling was used as a cut-off between low and high levels of vitamin D.

Selecting SNPs

To find SNPs previously associated with vitamin D, a search in a GWAS catalogue using the search string “vitamin D” identified 20 SNPs [2, 21]. The researchers in these studies had evaluated associations between SNPs and levels of vitamin D (15 SNPs) [22,23,24,25], vitamin D insufficiency (9 SNPs) [26] and levels of vitamin D-binding protein (DBP) (3 SNPs) [27]. Nine of these SNPs were directly genotyped on the BeadChip used in our study. Using the SNAP web-based tool, we identified proxies on the basis of linkage disequilibrium (LD) and physical distance to selected SNPs for an additional eight SNPs [28]. When several proxies were found, the two with the highest LD and closest proximity to the selected SNP were used. Altogether, 20 SNPs from 10 different loci in the genome were tested for associations with vitamin D level and breast cancer risk: 2 from ST6GALNAC3 in chromosome 1, 4 from group-specific component (GC) intron in chromosome 4, 1 from chromosome 6, 2 upstream of NPY in chromosome 7, 1 from proximity to MGMT in chromosome 10, 2 from unidentified genes, 2 from PDE3B intron, 2 from CYP2R1, 3 from NADSYN1 in chromosome 11 and 1 from MTMR4 in chromosome 17. A description of tested SNPs and their original GWAS proxy, as well as previous association with vitamin D, is provided in Additional file 1.

During QC, monomorphic SNPs due to lack of variation in a European population were excluded, as were SNPs deviating from Hardy-Weinberg equilibrium (p < 10−6) or if the variant call rate was < 95% in all samples. PLINK version 1.07 software was used for QC and identification of selected SNPs [29].

Statistical methods

Cases and control subjects were compared in two different sets with regard to established and potential risk factors. The first set contained all cases and control subjects with available SNP data; the other set included only cases and control subjects in whom vitamin D had been analysed in the previous study [8]. Control groups and groups with and without SNPs and vitamin D data were also compared with regard to risk factors in order to investigate the risk of potential selection bias.

Unconditional binary logistic regression analyses were performed to investigate associations between SNPs and low levels of vitamin D (below median), as well as between selected SNPs and risk of breast cancer. ORs with 95% CIs were calculated using the homozygote for the major allele as the reference group. Stratification by low versus high levels of vitamin D was added to see if there were any differences in breast cancer risk between these groups based on genotypes. This analysis was tested for interaction of low versus high vitamin D level. All analyses were also adjusted for year of inclusion, for age at baseline and for established or potential risk factors for breast cancer, including level of education, type of occupation, age at menarche, age at first childbirth, exposure to oral contraceptives, exposure to hormone replacement therapy, height, body mass index (BMI), alcohol consumption and smoking. Missing values for adjustment factors were coded in a separate category and included in multivariable analyses.


Potential confounders in cases and control subjects

Women with incident breast cancer were more often than control subjects to be non-manual labour workers (Table 1). This was applicable in both sets of study populations, although the difference between groups was smaller in the study population in which vitamin D levels were analysed. Another observation was that women with incident breast cancer in both study populations were more likely to be users of hormone replacement therapy and were also more likely to have been exposed to oral contraceptives (Table 1). In the complete study population, we also saw that women with incident breast cancer were more often included in the MDCS during the summer, had a higher level of education and were peri- or pre-menopausal to a greater extent at inclusion than women in the control group. These differences were not seen in the study population with analysed vitamin D levels, presumably owing to previous age and inclusion date matching of control subjects (Table 1).

Table 1 Baseline characteristics of breast cancer cases and control subjects

When comparing control subjects from the cardiovascular cohort with control subjects used in the previous case-control study, risk factors were similarly distributed, although small differences in a higher level of education and a slightly higher proportion of oral contraceptive users and smokers were seen amongst the previously matched control subjects (Additional file 2). Women whose data were excluded owing to QC of SNPs or relationships with other women in the cohort were less likely to have used oral contraceptives. They also had a lower BMI and had not used tobacco to the same extent as women with included data (Additional file 2). There also were small differences between women with or without available data on vitamin D level. Women with no data on vitamin D used oral contraceptives and hormone replacement therapy less and were also slightly older and less likely to smoke (Additional file 2).

SNPs and risk of low vitamin D level

There were statistically significant associations between a majority of tested SNPs and levels of vitamin D (Table 2). An increased risk of being in the group with lower vitamin D was seen for minor homozygotes of rs705117 (OR 3.38, 95% CI 1.20–9.49), rs7041 (OR 2.12, 95% CI 1.51–2.99) and rs2282679 (OR 2.51, 95% CI 1.56–4.05), all of which are located in the intron of the GC gene on chromosome 4 (Table 2). For rs7041 and rs2282679, the association was seen also for heterozygotes (respectively, OR 1.39, 95% CI 1.09–1.77; OR 1.42, 95% CI 1.13–1.80) (Table 2). Similarly, higher risk of low vitamin D was associated with minor homozygotes of rs12295888 (OR 1.80, 95% CI 1.24–2.63), rs1007392 (OR 1.78, 95% CI 1.26–2.52), rs7944926 (OR 2.39, 95% CI 1.59–3.59) and rs3829251 (OR 3.26, 95% CI 1.72–6.17) on chromosome 11 (Table 2). A statistically significant decreased risk of low vitamin D was seen with heterozygotes (OR 0.78, 95% CI 0.61–1.00) and minor homozygotes (OR 0.58, 95% CI 0.42–0.81) of rs2060793 also located on chromosome 11 but in the CYP2R1 gene (Table 2).

Table 2 Selected single-nucleotide polymorphisms in relation to vitamin D level and breast cancer

When several proxies were used for one GWAS SNP, the results in the analyses were very similar, but only the proxy with an R 2 value closest to 1 is presented in this report.

SNPs and overall breast cancer risk

The tested SNP proxy rs12239582 for the GWAS SNP rs12144344 (located in ST6GALNAC3 in chromosome 1) was found to be statistically significantly associated with a relatively low risk of breast cancer (OR 0.82, 95% CI 0.68–0.99) in the adjusted model when heterozygotes were compared with major homozygotes (Table 2). In the crude analysis, minor homozygotes of SNP proxy rs12791871 for GWAS SNP rs12785878 (located in NADSYN1 in chromosome 11) also had a statistically significant decreased risk of breast cancer (OR 0.73, 95% CI 0.55–0.95) compared with major homozygotes (Table 2). This result was similar in the adjusted analysis, although with only borderline statistical significance (OR 0.77, 95% CI 0.57–1.02) (Table 2).

SNPs and breast cancer risk in different strata of serum vitamin D

A test of interaction of high/low vitamin D levels and breast cancer risk was statistically significant at p < 0.01 for rs198300 (proxy for GWAS rs156299, located upstream of neuropeptide Y [NPY] gene in chromosome 7). Minor homozygotes of rs198300 who were within the group with low vitamin D levels had a statistically significant decrease in breast cancer risk compared with major homozygotes (OR 0.53, 95% CI 0.33–0.87) (Table 3). Looking at minor homozygotes of the same SNP (rs198300) within the group with high vitamin D levels, their risk of breast cancer seemed to be increased (OR 1.38, 95% CI 0.85–2.23) (Table 3). Similar but not statistically significant results were seen for rs2060793 (proxy for GWAS rs10741657 located in gene encoding cytochrome P450 [CYP2R1] in chromosome 11); compared with major homozygotes, heterozygotes in the group with lower levels of vitamin D had a slightly decreased risk of breast cancer (OR 0.95, 95% CI 0.67–1.35), but in the high vitamin D group, a reverse effect with increased breast cancer risk was seen (OR 1.35, 95% CI 0.93–1.97; p = 0.10 for interaction) (Table 3).

Table 3 Selected single-nucleotide polymorphisms in relation to risk of breast cancer, stratified by vitamin D level

Adverse associations of breast cancer risk dependent on vitamin D levels were seen for heterozygotes of rs10485165 (in chromosome 6). Heterozygotes with low levels of vitamin D were associated with an increased risk of breast cancer (OR 1.34, 95% CI 0.91–1.98), whereas heterozygotes with higher vitamin D levels had a decreased risk of breast cancer (OR 0.83, 95% CI 0.56–1.22), with a p value for interaction of 0.08 (Table 3).


In the present study, a majority of the selected SNPs were associated with vitamin D levels that were in agreement with those reported in previous studies. As for breast cancer risk, a statistically significant association was observed with one of the tested loci (rs12239582 and rs2209458; proxies for GWAS SNP rs12144344). Another locus (rs12791871 and rs7944926; proxies for GWAS SNP rs12785878) showed a tendency towards an association with breast cancer risk. The analyses stratified for level of vitamin D showed diverse associations with breast cancer risk for three of the tested loci: rs156299 (proxies used, rs198300 and rs13245518), rs10741657 (proxies used, rs2060793 and rs1993116) and rs10485165, depending on whether vitamin D level was high or low.

SNPs and risk of low vitamin D

There is no current consensus regarding which level of vitamin D should be considered insufficient; in addition, the level varies over the year owing to sun exposure. Grouping individuals by high or low vitamin D level by calendar month was considered pragmatic to calculate differences in vitamin D level owing to allele variants. SNPs located in the intron of the GC gene in chromosome 4 have previously been associated in several studies with levels of vitamin D and/or vitamin D insufficiency [22, 23, 25,26,27]. Results derived from the present analyses are consistent with previous results because minor alleles of all tested SNPs located in the intron of the GC gene (rs705117, rs7041, rs4588 and rs2282679) were statistically significantly associated with a low level of vitamin D. Equally confirming results were seen for several other SNPs of other genome sites, such as PDE3B intron, CYP2R1, and NADSYN1 in chromosome 11. Borderline statistically significant associations were seen for further SNPs on ST6GALNAC3 in chromosome 1 and on MTMR4 in chromosome 17. SNPs in chromosome 6, chromosome 7 and chromosome 10 did not show any association with low vitamin D, in contrast to previous studies of children in Western Australia [23, 24]. This might depend on the different genetic backgrounds of the two populations, limited statistical power due to small groups with vitamin D levels in our material, or different vitamin D exposures because it is expected that children in Western Australia are likely to have different exposure to sun than middle-aged women in the south of Sweden.

If there is a true association between low levels of vitamin D and an increased risk of breast cancer, a Mendelian randomization study would show that SNPs associated with low vitamin D levels would also be associated with an increased risk of breast cancer [12]. This is somewhat in line with the present results derived from analysing rs12239582 and rs2209458 in ST6GALNAC3, in chromosome 1, where minor homozygotes had a borderline statistically significant association with a decreased risk of low vitamin D level and heterozygotes with a statistically significant decreased risk of breast cancer (Table 2). Contradictory to this, SNPs whose minor alleles showed an association with an increased risk of low vitamin D—rs12791871, rs7944926 and rs3829251—of the NADSYN1 in chromosome 11 were, if anything, associated with a decreased risk of breast cancer. Also, for most SNPs with an association with vitamin D level, no association with breast cancer risk could be found. This might be explained either by no such association being able to be found or by such an association being too small to be observed in our study. The clinical relevance of such small associations is questionable.

SNPs and risk of breast cancer

The association found in our study between rs12144344 (proxies, rs12239582 and rs2209458) in chromosome 1 and breast cancer risk has not been reported previously. This SNP is positioned in ST6GALNAC3, a gene encoding a sialyltransferase which might affect DBP synthesis, concentration and function and that has previously been associated with levels of DBP [27]. No statistically significant association with vitamin D level was seen in our study, although minor homozygotes had a borderline statistically significant association with low levels of vitamin D. Previous associations between rs12144344 and vitamin D levels have not been reported.

Previous studies on vitamin D-associated SNPs and breast cancer risk have been focused mainly on genetic variants of the vitamin D receptor, and a few studies have been done on GC, which is the gene encoding DBP [30]. The present study was focused on SNPs previously associated with vitamin D in GWASs, and because no association was previously found for VDR SNPs, no VDR SNPs were included. Of the SNPs analysed in this study, only rs7041 and rs4588, both located in the GC gene, were previously studied in association with breast cancer risk [31,32,33,34,35]. Our results showing a tendency towards a protective effect of one or two minor alleles of rs7041 and no association of rs4588 are consistent with those of one previous study in which researchers showed a similar statistically significant association [32], as well as one other study in which researchers combined variants of rs7041 and rs4588 and found a protective effect of some combinations for post-menopausal breast cancer [36]. However, other studies have shown no association [30, 33, 35].

In the present study, the common homozygote was used as a reference, and relative risks were calculated individually for heterozygotes and less common homozygotes. Others have used other approaches, calculating on the basis of risk alleles and/or grouping several genotypes together. Sometimes it is not clear which allele was used as a reference, which means that comparing results between studies is difficult.

SNPs and breast cancer risk in different strata of serum vitamin D

When the association of analysed SNPs with breast cancer risk was stratified by level of vitamin D, some results differed depending on vitamin D level. This association was found for rs198300 and rs13245518, both proxies for rs156299 located on chromosome 7 upstream of gene encoding neuropeptide Y, a neurotransmitter involved in mediating physiological processes, including food intake and bone homeostasis, previously associated with levels of vitamin D [23]. In non-stratified analyses, no association with either breast cancer risk or vitamin D level was observed. Similar results were found for rs2060793 and rs1993116, proxies for GWAS rs10741657 located in CYP2R1, which encodes cytochrome P450, an enzyme that converts vitamin D to the circulating form 25(OH)D. rs10741657 was previously associated with vitamin D levels and insufficiency [22, 23, 25, 26]. Also, rs10485165, a proxy for rs7763511 (in chromosome 6) previously associated with vitamin D level but with unknown function [24], showed diverse ORs for breast cancer risk dependent on high versus low vitamin D level.

We found no previous reports of associations of this kind. Because this was an exploratory part of the present study, and because vitamin D level was available for only a smaller fraction of cases and control subjects, we suggest that these findings ought to be regarded with caution. The findings must be replicated, but at the same time they indicate that vitamin D level may interact with genetic traits in breast cancer risk.

Methodological issues

All patients diagnosed with cancer within the boundaries of Sweden are reported to the Swedish Cancer Registry, which is regarded as a highly validated source of information [37]. Thanks to the unique civil registration number given to all Swedes at birth, it is possible to link individuals included in the cohort to this registry and retrieve complete and correct information on breast cancer diagnosis. Women who have emigrated and received their diagnosis elsewhere are not included in this registry, which is a limitation of the present study; hence a few cancer cases might be lost, but emigration is not likely to be linked to risk of breast cancer or genetic traits.

An individual’s vitamin D level is influenced by several factors, not only sun exposure. It has also been shown that levels decrease with increasing age and increasing BMI [38,39,40], which was taken into consideration when analyses were adjusted. Physical activity has been shown to influence both level of vitamin D and breast cancer incidence [41, 42], but because the validity of physical activity data in the MDCS has been questioned, this was not included in the statistical model. Moreover, the effect of physical activity on vitamin D levels seems to be limited to outdoor activity, and hence physical activity may be a marker of sun exposure [43]. Vitamin D levels were analysed in only one blood sample per individual, which might be considered insufficient, though it has been shown previously that 25(OH)D levels analysed twice several years apart are individually very highly correlated [44, 45].

Differences noted when study populations were compared with regard to risk (Additional file 2: Table S1) were small. Those differences might depend on the fact that the previously matched control subjects had more in common with breast cancer cases. Considering differences between subjects with or without data regarding vitamin D level, there are of course a larger proportion of breast cancer cases amongst women with available vitamin D levels, which explains some of the differences. Also, because we adjusted analyses for risk factors in which subjects differ, we consider the risk of potential selection affecting our results to be small.

It is not plausible that any confounder included in this study affect an individual’s genetic composition. It might therefore be debatable whether adjustments should be added to the analyses. Because there is still a possibility of a chance association between single SNPs and established and potential risk factors for breast cancer, we chose to include such factors in our statistical model. A limitation of the present study is that there were no data accessible regarding family history of breast cancer or prior benign breast disease; hence those known risk factors were not included when adjustments were made.

In total, 20 SNPs from 10 different loci were tested for associations with breast cancer risk and risk of low vitamin D level in the present study. For each one of these, we used the same assumption: If there is an association with vitamin D, there might also be one with risk of breast cancer. We therefore decided not to increase the CI. With so many tested SNPs, one might suspect associations to be chance findings as a type I error, but when such associations were found, we also saw the same association for proxy SNPs, which strengthens the findings.

Individual SNPs change the risk of getting breast cancer to a very small extent, and large groups are a must in order to find statistically significant associations. Even the present population with 865 cases and 3193 control subjects may be too small to identify all associations. In the population with analysed levels of vitamin D, groups were even smaller, and therefore we expect a higher risk of a type II error in these analyses.


Many previous findings of SNPs associated with vitamin D levels were reassuringly replicated in this study. One SNP (rs12239582) previously associated with levels of DBP was associated with a slightly decreased risk of breast cancer, as was a tendency towards a decreased risk of low vitamin D level, in the present study. The allele frequency of the SNP variants is high, and therefore even a small increase in risk per individual may be a substantial population attributable risk. Further results of the present study show vitamin D level to be an effect modifier of the risk of breast cancer associated with certain SNPs. This suggests that an individual’s composition of SNPs may affect the extent to which levels of vitamin D are associated with breast cancer risk.



Body mass index


Vitamin D-binding protein


Group-specific component


Genome-wide association study


Hormone replacement therapy


Linkage disequilibrium


Malmö Diet and Cancer Study


Quality control


Single-nucleotide polymorphism


Vitamin D receptor


  1. 1.

    Eccles SA, Aboagye EO, Ali S, Anderson AS, Armes J, Berditchevski F, Blaydes JP, Brennan K, Brown NJ, Bryant HE, et al. Critical research gaps and translational priorities for the successful prevention and treatment of breast cancer. Breast Cancer Res. 2013;15(5):R92.

  2. 2.

    Burdett T, Hall PN, Hastings E, Hindorff LA, Junkins HA, Klemm AK, MacArthur J, Manolio TA, Morales J, Parkinson H, Welter D. The NHGRI-EBI Catalog of published genome-wide association studies. Accessed 11 Apr 2016 and April 23, 2017, version 1.0.

  3. 3.

    Grant WB. An ecological study of cancer incidence and mortality rates in France with respect to latitude, an index for vitamin D production. Dermatoendocrinol. 2010;2(2):62–7.

  4. 4.

    Porojnicu A, Robsahm TE, Berg JP, Moan J. Season of diagnosis is a predictor of cancer survival. Sun-induced vitamin D may be involved: a possible role of sun-induced vitamin D. J Steroid Biochem Mol Biol. 2007;103(3-5):675–8.

  5. 5.

    Robsahm TE, Tretli S, Dahlback A, Moan J. Vitamin D3 from sunlight may improve the prognosis of breast-, colon- and prostate cancer (Norway). Cancer Causes Control. 2004;15(2):149–58.

  6. 6.

    Rohan T. Epidemiological studies of vitamin D and breast cancer. Nutr Rev. 2007;65 Suppl 2:S80–3.

  7. 7.

    Colston KW. Vitamin D, and breast cancer risk. Best Pract Res Clin Endocrinol Metab. 2008;22(4):587–99.

  8. 8.

    Almquist M, Bondeson AG, Bondeson L, Malm J, Manjer J. Serum levels of vitamin D, PTH and calcium and breast cancer risk-a prospective nested case-control study. Int J Cancer. 2010;127(9):2159–68.

  9. 9.

    Gorham ED, Mohr SB, Garland FC, Garland CF. Vitamin D for cancer prevention and survival. Clin Rev Bone Miner Metab. 2009;7(2):159–75.

  10. 10.

    Shao T, Klein P, Grossbard ML. Vitamin D and breast cancer. Oncologist. 2012;17(1):36–45.

  11. 11.

    Wang D, Velez de-la-Paz OI, Zhai JX, Liu DW. Serum 25-hydroxyvitamin D and breast cancer risk: a meta-analysis of prospective studies. Tumour Biol. 2013;34(6):3509–17.

  12. 12.

    Ong JS, Cuellar-Partida G, Lu Y, Fasching PA, Hein A, Burghaus S, Beckmann MW, Lambrechts D, Van Nieuwenhuysen E, Vergote I, et al. Association of vitamin D levels and risk of ovarian cancer: a Mendelian randomization study. Int J Epidemiol. 2016;45(5):1619–30.

  13. 13.

    Khan MI, Bielecka ZF, Najm MZ, Bartnik E, Czarnecki JS, Czarnecka AM, Szczylik C. Vitamin D receptor gene polymorphisms in breast and renal cancer: current state and future approaches. Int J Oncol. 2014;44(2):349–63.

  14. 14.

    Shan JL, Dai N, Yang XQ, Qian CY, Yang ZZ, Jin F, Li M, Wang D. FokI polymorphism in vitamin D receptor gene and risk of breast cancer among Caucasian women. Tumour Biol. 2014;35(4):3503–8.

  15. 15.

    Xu Y, He B, Pan Y, Deng Q, Sun H, Li R, Gao T, Song G, Wang S. Systematic review and meta-analysis on vitamin D receptor polymorphisms and cancer risk. Tumour Biol. 2014;35(5):4153–69.

  16. 16.

    Zhang K, Song L. Association between vitamin D receptor gene polymorphisms and breast cancer risk: a meta-analysis of 39 studies. PLoS One. 2014;9(4):e96125.

  17. 17.

    Manjer J, Carlsson S, Elmstahl S, Gullberg B, Janzon L, Lindstrom M, Mattisson I, Berglund G. The Malmo diet and cancer study: representativity, cancer incidence and mortality in participants and non-participants. Eur J Cancer Prev. 2001;10(6):489–99.

  18. 18.

    Manjer J, Elmståhl S, Janzon L, Berglund G. Invitation to a population-based cohort study: differences between subjects recruited using various strategies. Scand J Public Health. 2002;30(2):103–12.

  19. 19.

    Fernandez C, Sandin M, Sampaio JL, Almgren P, Narkiewicz K, Hoffmann M, Hedner T, Wahlstrand B, Simons K, Shevchenko A, et al. Plasma lipid composition and risk of developing cardiovascular disease. PLoS One. 2013;8(8):e71846.

  20. 20.

    Lucas RM, Gorman S, Black L, Neale RE. Clinical, research, and public health implications of poor measurement of vitamin D status. J AOAC Int. 2017;100(5):1225–9.

  21. 21.

    Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, Klemm A, Flicek P, Manolio T, Hindorff L, et al. The NHGRI GWAS catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42(Database issue):D1001–6.

  22. 22.

    Ahn J, Yu K, Stolzenberg-Solomon R, Simon KC, McCullough ML, Gallicchio L, Jacobs EJ, Ascherio A, Helzlsouer K, Jacobs KB, et al. Genome-wide association study of circulating vitamin D levels. Hum Mol Genet. 2010;19(13):2739–45.

  23. 23.

    Anderson D, Holt BJ, Pennell CE, Holt PG, Hart PH, Blackwell JM. Genome-wide association study of vitamin D levels in children: replication in the Western Australian pregnancy cohort (raine) study. Genes Immun. 2014;15(8):578–83.

  24. 24.

    Benjamin EJ, Dupuis J, Larson MG, Lunetta KL, Booth SL, Govindaraju DR, Kathiresan S, Keaney Jr JF, Keyes MJ, Lin JP, et al. Genome-wide association with select biomarker traits in the Framingham heart study. BMC Med Genet. 2007;8 Suppl 1:S11.

  25. 25.

    Jorde R, Schirmer H, Wilsgaard T, Joakimsen RM, Mathiesen EB, Njølstad I, Løchen ML, Figenschau Y, Berg JP, Svartberg J, et al. Polymorphisms related to the serum 25-hydroxyvitamin D level and risk of myocardial infarction, diabetes, cancer and mortality: the Tromsø study. PLoS One. 2012;7(5):e37295.

  26. 26.

    Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, Kiel DP, Streeten EA, Ohlsson C, Koller DL, et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet. 2010;376(9736):180–8.

  27. 27.

    Moy KA, Mondul AM, Zhang H, Weinstein SJ, Wheeler W, Chung CC, Mannisto S, Yu K, Chanock SJ, Albanes D. Genome-wide association study of circulating vitamin D-binding protein. Am J Clin Nutr. 2014;99(6):1424–31.

  28. 28.

    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.

  29. 29.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.

  30. 30.

    Tagliabue E, Raimondi S, Gandini S. Meta-analysis of vitamin D-binding protein and cancer risk. Cancer Epidemiol Biomarkers Prev. 2015;24(11):1758–65.

  31. 31.

    Abbas S, Linseisen J, Slanger T, Kropp S, Mutschelknauss EJ, Flesch-Janys D, Chang-Claude J. 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. 2008;17(6):1339–43.

  32. 32.

    Anderson LN, Cotterchio M, Cole DE, Knight JA. Vitamin D-related genetic variants, interactions with vitamin D exposure, and breast cancer risk among Caucasian women in Ontario. Cancer Epidemiol Biomarkers Prev. 2011;20(8):1708–17.

  33. 33.

    Jorde R, Schirmer H, Wilsgaard T, Bøgeberg Mathiesen E, Njølstad I, Lochen ML, Joakimsen RM, Grimnes G. The DBP phenotype Gc-1f/Gc-1f is associated with reduced risk of cancer: the Tromsø Study. PLoS One. 2015;10(5):e0126359.

  34. 34.

    Maneechay W, Boonpipattanapong T, Kanngurn S, Puttawibul P, Geater SL, Sangkhathat S. Single nucleotide polymorphisms in the Gc gene for vitamin D binding protein in common cancers in Thailand. Asian Pac J Cancer Prev. 2015;16(8):3339–44.

  35. 35.

    Shi J, Grundy A, Richardson H, Burstyn I, Schuetz JM, Lohrisch CA, SenGupta SK, Lai AS, Brooks-Wilson A, Spinelli JJ, et al. Genetic variation in vitamin D-related genes and risk of breast cancer among women of European and East Asian descent. Tumour Biol. 2016;37(5):6379–87.

  36. 36.

    Abbas S, Nieters A, Linseisen J, Slanger T, Kropp S, Mutschelknauss EJ, Flesch-Janys D, Chang-Claude J. Vitamin D receptor gene polymorphisms and haplotypes and postmenopausal breast cancer risk. Breast Cancer Res. 2008;10(2):R31.

  37. 37.

    Barlow L, Westergren K, Holmberg L, Talbäck M. The completeness of the Swedish cancer register: a sample survey for year 1998. Acta Oncol. 2009;48(1):27–33.

  38. 38.

    Brock K, Huang WY, Fraser DR, Ke L, Tseng M, Stolzenberg-Solomon R, Peters U, Ahn J, Purdue M, Mason RS, et al. Low vitamin D status is associated with physical inactivity, obesity and low vitamin D intake in a large US sample of healthy middle-aged men and women. J Steroid Biochem Mol Biol. 2010;121(1-2):462–6.

  39. 39.

    Lips P. Vitamin D, physiology. Prog Biophys Mol Biol. 2006;92(1):4–8.

  40. 40.

    Need AG, O’Loughlin PD, Morris HA, Horowitz M, Nordin BE. The effects of age and other variables on serum parathyroid hormone in postmenopausal women attending an osteoporosis center. J Clin Endocrinol Metab. 2004;89(4):1646–9.

  41. 41.

    Romieu II, Amadou A, Chajes V. The role of diet, physical activity, body fatness, and breastfeeding in breast cancer in young women: epidemiological evidence. Rev Invest Clin. 2017;69(4):193–203.

  42. 42.

    Scragg R, Camargo Jr CA. Frequency of leisure-time physical activity and serum 25-hydroxyvitamin D levels in the US population: results from the third national health and nutrition examination survey. Am J Epidemiol. 2008;168(6):577–91.

  43. 43.

    Pilch W, Tota L, Sadowska-Krepa E, Piotrowska A, Kepinska M, Palka T, Maszczyk A. The effect of a 12-week health training program on selected anthropometric and biochemical variables in middle-aged women. Biomed Res Int. 2017;2017:9569513.

  44. 44.

    Meng JE, Hovey KM, Wactawski-Wende J, Andrews CA, Lamonte MJ, Horst RL, Genco RJ, Millen AE. Intraindividual variation in plasma 25-hydroxyvitamin D measures 5 years apart among postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2012;21(6):916–24.

  45. 45.

    Platz EA, Leitzmann MF, Hollis BW, Willett WC, Giovannucci E. Plasma 1,25-dihydroxy- and 25-hydroxyvitamin D and subsequent risk of prostate cancer. Cancer Causes Control. 2004;15(3):255–65.

Download references


The authors thank the data manager and registered nurse, Anna Hwasser, who kindly prepared original files and assisted in data management.


The present study was supported by The Einar and Inga Nilsson Foundation, The Gunnar Nilsson Cancer Foundation, The Ernhold Lundström Foundation, The Swedish Cancer Society, The Swedish Research Council (VR), The Malmö University Hospital Cancer Research Fund, The Skåne University Hospital Funds and Donations, The Breast Cancer network at Lund University (BCLU), and The Region Skåne (ALF).

Availability of data and materials

The authors have full access to the data, and the dataset used for the present study is kept by the corresponding author.

Author information

LH selected SNPs to analyse, extracted information of such SNPs from a PLINK file and did all statistical analyses using SPSS software. LH was the main author of the manuscript and constructed the tables and figures. STB and JM supervised the project. PA performed quality control of SNPs and exclusion of relatives and supervised PLINK software use. SB and JB revised and added material to the manuscript. AF supervised selection of SNPs and provided information on handling SNP data. AF revised the manuscript. OM performed SNP analyses. All authors read and approved the final manuscript.

Correspondence to Linnea Huss.

Ethics declarations

Authors’ information

LH is a general surgeon in the Department of Surgery at Skane University Hospital and is also a doctoral student.

Ethics approval and consent to participate

The MDCS and the present study were approved by the ethics committee in Lund, Sweden (LU 51-90, Dnr 153/2004 and Dnr 682/2009). At the time of inclusion in MDCS, written informed consent was obtained from all participants. Follow-up studies have been advertised in local media with an option to withdraw participation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Appendix. Description of GWAS SNPs associated with vitamin D, SNP proxies and SNPs analysed. (DOCX 27 kb)

Additional file 2: Table S1.

Baseline characteristics of control subjects and subjects with and without data on SNPs/vitamin D. (DOCX 24 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark


  • Vitamin D
  • Breast cancer
  • Polymorphism
  • SNP