Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Int J Cancer. Author manuscript; available in PMC 2013 January 15.
Published in final edited form as:
PMCID: PMC3135781

Vitamin D receptor polymorphisms in patients with cutaneous melanoma


The vitamin D receptor (VDR) gene has been associated with cancer risk, but only a few polymorphisms have been studied in relation to melanoma risk and the results have been inconsistent. We examined 38 VDR gene SNPs in a large international multi-center population-based case-control study of melanoma.

Buccal DNAs were obtained from 1207 people with incident multiple primary melanoma and 2469 with incident single primary melanoma. SNPs with known or suspected impact on VDR activity, htSNPs with ≥10% MAF in Caucasians, and SNPs reported as significant in other association studies were examined. Logistic regression was used to calculate the relative risks conferred by the individual SNP.

Eight of 38 SNPs in the promoter, coding, and 3’ gene regions were individually significantly associated with multiple primary melanoma after adjusting for covariates. The estimated increase in risk for individuals who were homozygous for the minor allele ranged from 25% to 33% for 6 polymorphisms: rs10875712 (OR 1.28; 95%CI, 1.01–1.62), rs4760674 (OR 1.33; 95% CI, 1.06–1.67), rs7139166 (OR 1.26; 95%CI, 1.02–1.56), rs4516035 (OR 1.25; 95%CI, 1.01–1.55), rs11168287 (OR 1.27; 95%CI, 1.03–1.57), rs1544410 (OR 1.30; 95%CI, 1.04–1.63); for 2 polymorphisms, homozygous carriers had a decreased risk: rs7305032 (OR 0.81; 95%CI 0.65–1.02), rs7965281 (OR, 0.78; 95%CI, 0.62–0.99). We recognize the potential false positive findings due to multiple comparisons; however the 8 significant SNPs in this study outnumbered the 2 significant tests expected to occur by chance. The vitamin D receptor may play a role in melanomagenesis.

Keywords: VDR, SNP, melanoma, polymorphism, vitamin D


Vitamin D is hypothesized to lower the risk for or mortality from several cancers, including melanoma. In vitro and in vivo assays show that vitamin D modulates cell proliferation, differentiation and programmed cell death 15, supporting its anti-carcinogenic potential. A number of observational studies in humans have also suggested that vitamin D, or putative surrogates for vitamin D status such as season, geographic latitude and evidence of chronic sun exposure, is associated with better outcomes in melanoma patients, despite the fact that UVB exposure, which activates the precursor of vitamin D present in skin, is a known risk factor for melanomagenesis 69. This adds to the plausibility of the hypothesis that vitamin D affects melanoma risk.

The active form of Vitamin D, 1,25-dihydroxyvitamin D, exerts its biological function by binding to the vitamin D receptor (VDR), a nuclear steroid hormone receptor expressed in many tissues and organs including skin 10. This receptor is encoded by the vitamin D receptor gene (ID 7421, OMIM 601769) which spans approximately 100 Kb on chromosome 12q13.11, with a large non-coding region spanning exons 1F to 1C, and with exons 2 to 9 encoding a 424 aminoacid protein 11. The VDR gene contains numerous variants, some of which are hypothesized to influence the expression and/or function of the VDR protein product, and some of which have been associated with complex traits and diseases such as osteoporosis, diabetes, small stature, and cancer 12, 13.

Despite the large number of SNPs, to date only a few common polymorphisms in the VDR gene have been studied in relation to melanoma risk: FokI, BsmI, TaqI, ApaI, A-1012G, and Cdx2, and the results for some of the SNPs have not been consistent across different studies 1423. For example, the FokI (rs2228570) minor allele T (f, M4) results in production of a VDR protein that is less effective as a transcriptional activator 2426 and has been linked to an increased risk for melanoma in several studies 14,16,21,23, and in pooled analyses 23,27,28, but has been observed to have no effect according to other reports 19,20,22,23. Some of the inconsistencies may be due to differences in the study design, study populations, or in the exposures or adjustments for covariates during the analyses. It is also possible, especially when the size of the VDR gene region is considered, that the causal genetic variant or key group of variants remains to be identified.

Nejentsev et al. (2004)12 identified the minimal set of tag SNPs and then evaluated the informativeness of the most commonly studied VDR SNPs in relation to disease association. The authors found that FokI was not in LD with any other SNP, while ApaI, TaqI, and BsmI predicted only 38% of the SNPs in block B with R2>0.8. Therefore, in studies that genotyped only FokI, TaqI, BsmI or ApaI, information on a large fraction of common SNPs was not captured. The aim of this study was to conduct a comprehensive analysis of common genetic variants including less studied SNPs across the VDR gene and to determine their effects on melanoma risk.


Study population

Subjects were recruited between 1998 and 2003 through an international multi-center population-based case control study of multiple primary melanoma, the Genes, Environment and Melanoma Study (GEM). The GEM population consists of incident cases of multiple primary melanoma identified in eight population- based registries and one hospital center in Australia, Canada, Italy and USA. Patients with incident single primary melanoma were used as controls. This design is, in effect, a case-control study of melanoma conducted in the high risk population of individuals who have already experienced a melanoma diagnosis. Further details of the study design and its rationale have been published 29,30.

The study was approved by the institutional review boards at all centers. All study participants signed informed consent, provided a buccal sample for extraction of germline DNA, and completed a detailed interview, providing information on sun exposure history, sun sensitivity, and phenotypic characteristics.

Selection of VDR SNPs

DNA samples from 1207 individuals with multiple primary melanoma and 2469 with single primary melanoma were available for genotyping. A total of 42 SNPs located on average at 3.6Kb from each other were selected. These included htSNPs with MAF>10% in Caucasians 12,13, SNPs with known or suspected impact on the transcription, stability, and/or activity of the VDR, and SNPs reported as significant in other association studies. Specifically, rs11568820 located on the promoter has been found to affect transcription of VDR by interfering with the binding site of the transcription factor Cdx2 11. SNPs rs7139166 (1A-1521) and rs4516035 (1A-1012), when the haplotype 1521G/-1012A is present, were shown to affect promoter activity 31. The polymorphism rs2228570 (FokI), responsible for an alternative earlier start and a less active protein was included. The commonly studied polymorphisms rs731236 (TaqI) and rs1544410 (BsmI) were also included although their functional impact is less apparent.

Laboratory methods

The isolation of DNA was previously described 32. We used the Sequenom MassARRAY iPLEX genotyping Platform (Sequenom Inc., San Diego, CA) 33 to test 36 SNPs and pyrosequencing and melting temperature analysis to test 2 SNPs. For iPLEX, experiments were designed with the RealSNP Assay Database and the MassArray software v.3.4. Five to 10 ng of genomic DNA were amplified using specific primers, reagents, and cycling conditions detailed in the Supplementary Material (Table S1). Products were spotted onto 384 SpectroCHIP bioarrays (Sequenom). Cluster plots were evaluated with the TyperAnalyzer application (MassARRAY v.3.4). Assays were considered optimal according to degree of clustering, absence of signal in the blanks, and when sequencing of representative samples within the clusters confirmed the genotypes. Four SNPs were rejected during the assay design or during the wet testing due to non-specific signals (rs2238139, rs739837) or very low signals and call rates (rs2408876, rs3847987). SNPs rs2228570 (267bp) and rs1544410 (209bp) were genotyped using single SNP genotyping methods: melting temperature analysis34 coupled to the LightTyper (Roche Applied Science, Indianapolis, USA) and pyrosequencing 35 with the PSQ™ MA instrument (Biotage AB, Uppsala, Sweden). Specific primers and probes, reaction and cycling conditions are available in Table S1. Quality control (QC) procedures included use of barcodes on samples and plates, inclusion of internal laboratory controls, randomly selected repeats (5% for MassArray and 10% for pyrosequencing and for melting temperature analyses), pre-PCR dedicated materials and working space, and at least two independent readers for review and interpretation of cluster plots, pyrograms, melting curves, and results. A plate containing discordant results for at least one pair of duplicate samples or internal controls was newly assayed. Samples with weak signals and outliers were repeated. QC for data entry included an additional review of 20% of genotype calls. There was 100% agreement in the genotyping calls between two independent laboratory members. Direct sequencing was performed on an independently amplified PCR fragment to confirm the genotype of laboratory control samples and when a new SNP was identified by different sequence pattern in the pyrogram, or different melting profile using standard methods 36.

Statistical analysis

We analyzed a total of 38 VDR SNPs in relation to risk of melanoma using conventional methods for case-control studies. The “variant” allele was defined to be the less frequently occurring allele. For each SNP, logistic regression was used to calculate adjusted odds ratios (OR) with 95% confidence intervals (95% CI). Heterozygous and homozygous carriers of the variant allele were separately compared to homozygous carriers of the most prevalent or “reference” allele within the same SNP-specific logistic regression model. An ordinal variable for the three possible genotypes was also created to evaluate linear trend. A multiplicative age-sex interaction term was included in all models to control for the fact that the age incidence curves for melanoma differ for men and women. In addition to age, sex, and age-sex interaction, study center was included as a covariate in the multivariate analysis to account for potential differences among populations. Analyses further adjusting for known melanoma risk factors, such as mole count, hair and eye color, and other “phenotypic” characteristics, were also conducted, but the results are not materially different and are not presented. Hardy-Weinberg equilibrium (HWE) tests were performed to identify possible gross systematic genotyping problems; since our control group consists of melanoma cases, absence of HWE could be expected in SNPs that confer risk. We evaluated the issue of multiple comparisons using techniques described in Abramovich and Benjamini (2006)37. We also employed quantile-quantile (Q-Q) plots 38 for comparing the distribution of p-values obtained in the logistic regression against those expected if the data occurred randomly and not specifically due to associations between SNPs and melanoma risk. All analyses were performed using SAS statistical software version 9.2 (SAS Institute, Cary, NC).


We conducted a comprehensive, systematic search of the medical literature to identify relevant studies on cutaneous melanoma in which genotyping results were available for each of the SNPs identified as statistically significant in our study. The literature search was limited to studies in humans, included all years up to November 2010, and was performed using the following databases: PubMed, ISI Web of Science, ISI Web of Knowledge, Cochrane, and Embase. Search terms included ‘Vitamin D Receptor’ and separately, ‘VDR’, in combination with ‘melanoma’, ‘allele’, ‘gene’, ‘polymorphism’, and ‘risk’. We searched for publications with frequencies or estimates of relative risk for cutaneous melanoma, with 95% CI. Review articles were considered as sources of additional relevant citations. We extracted information on the characteristics of the study population, study design, year of publication, and source of controls (Table S2). The homozygous occurrence of the most common allele for each polymorphism was considered the reference for our analyses. Hardy-Weinberg tests were performed among controls. For one study in which genotype frequencies were not reported 20, we contacted the authors to obtain this information. Among the articles retrieved, two papers 17,21 presented overlapping data on rs731236 and we therefore included results of the report with the largest number of subjects 21. The variation across different studies attributable to heterogeneity rather than chance was assessed using the I-squared measure 39. We considered that heterogeneity was present when the overall I2 p-value was <0.05. Calculation of summary odds ratios, 95% confidence intervals and study heterogeneity p-values for meta-analysis was done in STATA (version 11.0; STATA Corp, College Station, TX) using the metan function. Publication bias was evaluated using the adjusted rank correlation test within the metabias function of STATA 40.


The SNPs analyzed in this study do not affect the amino acid sequence. Therefore, we evaluated their potential functional relevance on the VDR transcription or activity by searching for overlaps between the SNP genomic location and transcription factor binding sites, splice sites, sequences complementary to miRNA seed regions, and by considering conservation across species. We obtained information from the UCSC’s genome browser tables (accessed October 2009), and each SNP, when information was available, was further annotated with conservation information, and miRNA and transcription factor binding site overlaps. Information on conservation was obtained from the UCSC’s runs of the PHAST package on primates, vertebrates and mammals. The miRNA annotations were obtained from Target Scan, and information on transcription factor binding from the binding site conservation across the human/mouse/rat alignment. Further information can be obtained from: 41. We also used an alternative method for detecting potential overlaps between VDR SNPs and miRNAs. Specifically, we looked for miRNA seed regions by searching for exact matches between the 2–7 nucleotide regions of known mature human miRNAs and the reverse complement of these substrings in the 3’UTR of the VDR gene. Selbach et al. (2008) 42 has shown that this 6-mer miRNA seed is essential for its function. Consequently, a SNP located in the miRNA target site which is complementary to the miRNA seed region is likely to deregulate the expression of the target gene.


Descriptive characteristics of the study participants are provided in Table 1. The distribution of melanoma risk factors such as age, gender, and number of moles on the back differed between GEM cases and GEM controls. Further details on the associations between participants’ characteristics and melanoma risk in this study have been published elsewhere 43.

Table 1
Selected characteristics of study participants

The individual call rates for the SNPs evaluated in this study, their distribution in the population according to HapMap_CEU and other databases, and the genotype frequencies obtained in the GEM cases and controls are shown in Table 2. After excluding four SNPs with poor assay performance, the average call rate was 0.98 for the remaining 38 SNPs. For SNPs rs11568820 and rs7305032 the failure rates were 0.12 and 0.07, respectively. However the proportion of failures did not differ significantly between cases and controls (p=0.33 and p=0.34). The proportion of failures was also similar between cases and controls for SNPs rs2228570 (p=0.85) and rs1544410 (p=0.33). All other SNPs displayed failure rates lower than 0.03. Genotype frequencies were similar to those reported in dbSNP in Caucasians or Europeans (HapMap-CEU or AFD_EUR_PANEL), with few exceptions (rs11168275, rs12370156, and rs7305032). For rs11168275 we and others 44 have found that G, and not A, is the minor allele in Caucasians. All three SNPs are in HWE (Table S1). The following SNPs were in LD with each other (D’ > 0.96): rs7139166 and rs4516035 (D’ = 0.972; R2 = 0.942); rs2544027 and rs2544028 (D’ = 0.968, R2 = 0.708); rs2544027 and rs2544038 (D’ = 0.974, R2 = 0.815); rs2544028 and rs2544038 (D’ = 0.974, R2 = 0.815); rs1544410 and rs731236 (D’=0.957, R2=0.912).

Table 2
Genotype frequency of VDR variants in individuals with multiple primary (GEM cases) and single primary (GEM controls) melanoma and adjusted odds ratios by genotype.

The genotype distribution among cases and controls, and the relative risks based on individual VDR genotypes are presented in Table 2. Nominally significant results (at the 5% significance level) are in bold type. Eight SNPs displayed statistically significant associations (p-value for trend p≤0.05) after adjusting for age, sex, age-sex interaction, and study center. The estimated increase in risk for individuals who were homozygous for the variant allele ranged from 25% to 33% for 6 polymorphisms: rs10875712, rs4760674, rs7139166, rs4516035, rs11168287, rs1544410; for 2 polymorphisms, homozygous variant carriers had a decreased risk: rs7305032, rs7965281. Results of all significant and marginally significant SNPs are displayed in Figure 1, in which the 8 significant SNPs appear in bold typeface. Figure 2 depicts the Q-Q plot, which should track with the 45° line in the absence of association across the locus. The shape of the curve exhibits a strong departure from what we would expect in the absence of association. Five promoter SNPs, two coding SNPs, and one 3’UTR SNP showed significant associations. With respect to previously described haplotype blocks, five of the significant SNPs lie within ‘block C’ and two of these overlap with ‘block 2’ 13; while three SNPs lie within ‘block B’ 12.

Figure 1
Position of the SNPs studied GEM in relation to the VDR structure. SNPs with significant associations in multivariate analyses after adjusting for age, sex, age-sex interaction, center are shown in bold font. The vertical bars represent promoters 1F – ...
Figure 2
Uniform Q-Q Plot for the VDR SNPs risk trend p-values. Quantiles of the obtained p-values were plotted against the theoretical uniform quantiles.

The VDR SNPs tested in this study included promoter, intronic, and 3’UTR variants and therefore do not directly affect the protein sequence. We evaluated their potential functional significance by searching for overlaps between genomic localization of SNPs and putative transcription factor binding sites, and miRNA seed regions. We also evaluated the conservation across species. We found that rs2228570 localizes to a nucleotide conserved among primates, vertebrates, and mammals (conservation scores 392, 418, and 379 respectively); SNPs rs731236 and rs2544028 localize to nucleotides conserved among vertebrates (conservation scores 272, 480) and mammals (conservation scores 275, 443); and SNPs rs4073729 and rs731236 localize to nucleotides conserved among primates (conservation scores of 293 and 407). The genomic locations of two SNPs overlap with transcription factor binding sites: rs4237856 (V$HNF4_01 HNF-4alpha 1 and 2, transcription z-score 1.98; and V$PPARG_01 PPARgamma, transcription z-score 2.58); and rs1544410 (V$IK3_01 Ikaros 3, transcription z-score 2.77; and V$IK1_01 Ikaros 1, transcription z-score 2.14), although the z-scores for HNF-4alpha 1 and 2 and for Ikaros 1 fall below the threshold of 2.33. SNP rs11568820 is known to interfere with the binding site of the transcription factor Cdx2 11. We did not find any overlaps with seed regions of known miRNAs.

At the time of the preparation of this manuscript, thirteen original reports addressed VDR SNPs and the risk of melanoma, and eleven articles containing a total of 13 series reported results for three of our significant or borderline significant SNPs: rs4516035, rs1544410, and rs731236. Of these, two reports overlapped with a subsequent publication from the same investigators and were not considered for the current meta-analyses 17,21. The main study characteristics and the distribution of genotypes among cases and controls are listed in the supplemental material (Table S2; Table S3). The individual unadjusted ORs and summary odds ratio for heterozygotes and for carriers of the homozygote variant compared to the homozygote common allele are presented in Table 3 and in Figure 3.

Figure 3
Meta-analysis of VDR SNPs reported in other investigations of melanoma in relation to risk. (A) SNP rs4516035 (variant allele G), (B) SNP rs1544410 (variant allele A or ‘B’); and (C) SNP rs731236 (variant allele C or ‘t’). ...
Table 3
Individual odds ratios and summary odds ratios for heterozygote and homozygote variants versus common alleles for SNPs rs4516035 and rs731236

For studies that addressed rs4516035 (1A-1012), the adjustment for covariates varied substantially. Specifically, four studies 18,19,22,23 matched cases and controls by age and gender; three studies 14,15,21 included age and sex in multivariate analyses; while one study 20 adjusted for phenotypic features. The variant allele G produces a null association with melanoma. There was evidence of heterogeneity for the observed risks of the homozygous variant (p=0.04) but not for heterozygote (p=0.31), and no evidence of publication bias (p=0.88 and p=0.65, respectively). For rs1544410 (BsmI) the allele A (‘B’) showed a statistically negative association with melanoma (Figure 3-B, Table 3). There was no evidence of heterogeneity for the observed risks of the heterozygous variant (p=0.78) or for the homozygote variant (p=0.18), and no evidence of publication bias (p=0.45 and p=0.45, respectively). The rs731236 (TaqI) allele C (‘t’) showed a statistically non-significant negative association with melanoma (Figure 3-C, Table 3). There was evidence of heterogeneity for the observed risks of the heterozygous variant (p=0.03) but not for the heterozygote variant (p=0.10), and no evidence of publication bias (p=0.19 and p=0.85, respectively). We did not observe departure from HWE in any of the series (Tables S3).


In this large population based case-control study comprising 3676 people with incident single and multiple primary melanoma we evaluated 38 VDR gene single nucleotide polymorphisms for their impact in the development of skin melanoma and found that in addition to the BsmI (rs1544410), seven SNPs (rs10875712, rs4760674, rs7139166, rs4516035, rs11168287, rs7305032, rs7965281) were significantly associated with risk. Because we conducted 38 individual statistical tests we recognize the increased risk of false positive findings due to these “multiple comparisons”; we would expect about 2 of the tests to be significant simply due to chance. The fact that we observed 8 significant SNPs, reflected by the skewness in the Q-Q plot (Figure 2), suggests that the evidence broadly favors the hypothesis that the gene harbors at least some causal variants.

Promoter SNPs

Most of the VDR SNPs identified in the last several years are of unknown functional significance; however in vitro studies suggest that at least some SNPs affect the VDR mRNA and/or protein levels 13,24,25,26. Such is the case for two of the SNPs possibly associated with multiple primary melanoma in our study: rs7139166 and rs4516035, located in the promoter region in positions 1A-1521 and 1A-1012. The genotype frequencies in our study are similar to those reported in dbSNP in Caucasians and to those reported by others for rs4516035 15,18,22,23,45. These two neighboring SNPs are in strong LD (D’ = 0.968), therefore it is not surprising that the genotype frequencies and the ORs conferred by each of the minor alleles are very similar in our study (Table 2). A previous report suggested a link between these SNPs and stature 31 although results from a larger study did not support this finding 46. We found that the G allele of rs4516035 appears to confer risk. Recently, Randerson-Moor et al. (2009) 23 in a population-based study of 1028 melanoma cases and 402 controls found no association for rs4516035, while in a second case-control study of 299 melanoma cases and 560 unaffected females, the authors reported an increased risk associated with the same allele (G), though the results fall short of statistical significance. We further compared our results to those obtained in our meta-analysis of a total of six published reports on this SNP and melanoma and found no evidence overall that this SNP is associated with melanoma risk, although the results are inconsistent, with significant heterogeneity in the observed relative risks across studies.

Our findings are consistent with the in vitro data for rs4516035 13,31. Specifically, d’Alesio et al. (2005)31 observed that one base-change in any of the two variant sites located in positions 1A-1521 or 1A-1012 led to a dramatic change in protein-DNA complex formation, and that the promoter activity of the VDR was nearly doubled in cells carrying the haplotype rs7139166-G/rs4516035-A when compared to those with the rs7139166-C/rs4516035-G haplotype.

It has been recently hypothesized that the A allele in position 1A-1012 (rs4516035) indirectly promotes a GATA-3 driven T-cell switch of naïve T cells to T-helper 2 cells and that this allele may play different roles in susceptibility and in metastasis possibly as a function of the transcription factors secreted by different cellular backgrounds 13,45,47. Our results do not support an adverse role of the allele ‘A’ in position 1A-1012 in relation to melanoma risk, and in support of our findings, the qualitative analyses of the meta-analyses conducted on published melanoma case-control series suggests that ‘G’ may represent the risk allele.

We found four significant VDR promoter SNPs (rs10875712, rs4760674, rs7139166, rs11168287) that have not been previously reported in relation to melanoma risk or progression.

Coding region and 3’UTR SNPs

Three adjacent RFLP SNPs for BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) respectively, at the 3’ end of the gene, have been the most frequently studied in cancer and in melanoma. The functional characterization of these SNPs has produced inconsistent results across studies, and if these SNPs are non-functional, the effect observed in some investigations may be due to truly relevant SNPs in strong LD located elsewhere in the gene. Here we report that in this gene region, in addition to the minor allele of the BsmI polymorphism (rs1544410) which increases risk, the minor allele of rs7965281, a SNP which has never been studied in relation to melanoma prior to this work, seems to confer some reduction in risk. We also report that two additional SNPs in this location may modify risk: rs7305032, which was not previously reported in relation to melanoma, and rs731236 (TaqI), although the latter, which is in high LD with BsmI and maps to a nucleotide conserved across species, exhibited only borderline significance (p-value 0.06, Table 2). Some investigations targeted to determine the expression differences for polymorphisms in the 3’ end of the VDR gene have shown higher levels of mRNA expression for the BsmI-ApaI-TaqI haplotype ‘BAt’ (rs1544410-A/rs7975232-A/rs731236-C) than for the haplotype ‘baT’(rs1544410-G/rs7975232-C/rs731236-T), however, the opposite results have been found not only for the VDR mRNA expression but also for the mRNA stability and transactivation in human fibroblasts, and in leukemia and prostate cancer cell lines 11. Our findings for rs1544410 suggest that the A (‘B’) might be the adverse allele in relation to melanoma risk which differs from what other studies have found with one exception 16,19,2123. Specifically, three studies were null, one study reported reduction of risk in heterozygous bB individuals but not in homozygotes BB, and one small hospital-based study conducted in 101 cases and 101 controls found a protective effect for the B allele; it is of note that the reported MAF among controls was higher than the one reported in dbSNP and in other studies including ours. A more recent population-based study of 763 cases and 763 controls conducted in a Polish population found that haplotypes containing the B allele increased risk for melanoma 22. Results from our meta-analysis conducted on five previously published studies, like in previous meta-analyses 23,2728, suggest a protective role for the B allele; however, considering the results from our large study after adjusting for covariates, we argue that the results from the meta-analyses should be considered with caution and that BsmI should be further characterized in large population-ascertained cohorts and in functional assays preferentially using melanoma and/or skin models. Our findings for rs731236 suggest that the C (‘t’) allele might be the adverse allele in relation to melanoma risk. Recent reports have found similar genotype frequencies to those in the present work, however investigators found either no significant associations with melanoma for rs731236 14,20,22,23, or, in one study, a significant association for the contralateral allele T (‘T’) 21. Our meta-analysis conducted on five publications does not provide sufficient evidence for an effect of this SNP and melanoma, however it is important to note that our meta-analyses were unadjusted for covariates, and that the results may differ in a multivariate logistic regression.

In relation to bioinformatics predictions, there are two opposing theories concerning conservation of sequences across species in regards to functional relevance. The first supports the idea of a plausible biological function, and as an extension of this idea, one would expect to uncover melanoma risk SNPs in this region. A second opposing hypothesis proposes that variants in a conserved region are subjected to negative selection, and those that persist will more likely be non functional. Interestingly, two of the investigated SNPs (rs731236 and rs2544028) are located on the boundary of the conserved region and although the associations do not reach statistical significance, the variants increase the risk by 21–27% (p-values 0.06 and 0.08, respectively). SNP rs11568820 has been extensively characterized in in vitro studies and overlaps with the binding site of Cdx2. One would expect to find a significant role for this polymorphism in melanoma risk; however it is possible that its effect is cell-specific. In fact, Cdx2 is differentially expressed in normal stomach and intestine and seems progressively decreased in gastric intestinal metaplasia, dysplasia and cancer 17,48,49. SNP rs4237856 overlaps with the binding site of PPAR-gamma. The PPARs are known to heterodimerize with the retinoid X receptor and with the vitamin D bound to its receptor. PPARs may be of importance in cell growth regulation in melanocytes; however expression of PPAR-gamma is much weaker in melanocytes, at least based on observations in MeWo cells 50; therefore, one could speculate that even if this SNP can facilitate or hinder the binding of PPAR-gamma, the contribution of this binding in melanoma risk may not be as significant as with other PPARs.

A major limitation so far of association studies using VDR SNPs in relation to complex-diseases such as cancer is the limited number of variants studied. The VDR gene is very polymorphic; however information on the existing gene variants was initially available from the rather insensitive search of polymorphic patterns within a limited region using restriction enzymes. Additional polymorphisms have been identified during the search for mutation in tumors 24. In more recent years, the studies of Nejentsev et al. (2004)12 and Fang et al. (2005)13 have contributed to the field by systematically scanning and characterizing variations across the entire VDR gene, and by describing patterns of linkage disequilibrium between alleles in different ethnic groups. The main strengths of our study are the comprehensive evaluation of VDR and the large number of study subjects, as well as the availability of data on known melanoma risk confounders. The overall effect on melanoma risk found here was small. Furthermore the results overall do not reach significance if evaluated either by the Bonferroni correction or by restricting the false discovery rate to 5% 37. However eight of 38 (~21%) investigated SNPs revealed statistically significant associations at the 5% significance level (versus 1.9 expected), and 16 of 38 are significant at the 10% level (versus 3.8 expected), proportions considerably higher than would be expected if there is in fact no association between VDR genotype and melanoma risk. Thus we believe that the results, collectively, are strongly suggestive that VDR is a risk locus for melanoma.

Our case-control study design was novel since the controls have melanoma and the cases are required to have experienced melanoma twice. The premise is that risk factors that affect melanoma must affect the occurrence of a second melanoma in a population of individuals with melanoma, and this is borne out empirically in earlier results from this and other studies, although there is evidence that the magnitude of the observed relative risks may be attenuated29. The rationale for this design is that it delivers much greater statistical power for the evaluation of rare, strong risk factors than a conventional case-control design with healthy controls 30. In such a design the double primary cases are typically older than the single primary controls, and in the case of melanoma the gender distribution is altered because of the strong interaction between age and gender on melanoma incidence, and so we adjusted for these factors in our analyses. A further complication is that many second primaries are detected as in situ lesions due to the careful clinical surveillance that melanoma patients experience, and we elected to include these individuals as cases on the assumption that they would ultimately have been diagnosed with an invasive lesion.

In summary, we found that eight VDR SNPs, located in the promoter, coding and 3’ gene regions, were associated with melanoma and confer a modest but statistically significant increased (rs10875712, rs4760674, rs7139166, rs4516035, rs11168287, rs1544410) or decreased (rs7305032, rs7965281) risk of becoming a multiple primary melanoma. To our knowledge, this constitutes the largest study of VDR polymorphisms in melanoma to date. Our results provide some evidence in support of the hypothesis that the vitamin D pathway may play an important role in melanoma genesis.

Supplementary Material

Supp Table S1

Supp Table S2

Supp Table S3


We thank Dr. Laetitia Borsu for assistance with the Sequenom assays; Javier Cotignola for help provided during the initial stages of the iPLEX assay design and genotyping; Brian Clas and Concetta Perretta for their help during the individual SNP genotyping; Aaron Gabow and the MSKCC Bioinformatics Core for providing assistance with the computational characterization of SNPs; Joseph Kanik for assistance with graphics; and all the individuals who agreed to participate in the study.

The study was conducted by the GEM Study Group: Marianne Berwick (PI, University of New Mexico), Memorial Sloan-Kettering Cancer Center, New York, NY, USA: Colin Begg (Co-PI), Irene Orlow (Co-Investigator), Urvi Mujumdar (Project Coordinator), Klaus Busam (Dermatopathologist), Pampa Roy (Laboratory Technician). Study Centers: The University of Sydney and The Cancer Council New South Wales, Sydney (Australia): Bruce Armstrong (PI), Anne Kricker (co-PI), Melisa Litchfield (Study Coordinator). Menzies Research Institute, University of Tasmania, Hobart (Australia): Terence Dwyer (PI, currently at the Murdoch Childrens Research Institute, Melbourne, Victoria), Paul Tucker (Dermatopathologist), Alison Venn (co-Investigator), Nicola Stephens (Study Coordinator). British Columbia Cancer Agency, Vancouver (Canada): Richard Gallagher (PI), Teresa Switzer (Coordinator). Cancer Care Ontario, Toronto (Canada): Loraine Marrett (PI), Elizabeth Theis (Co-Investigator), Lynn From (Dermatopathologist), Noori Chowdhury (Coordinator), Louise Vanasse (Coordinator). Centro per la Prevenzione Oncologia Torino, Piemonte (Italy): Stefano Rosso (PI), Roberto Zanetti (co-PI), Carlotta Sacerdote (Coordinator). University of California, Irvine (USA): Hoda Anton-Culver (PI), Nancy Leighton (Coordinator). University of Michigan, Ann Arbor (USA): Stephen Gruber (PI), Joanne Jeter (Coordinator). New Jersey Department of Health and Senior Services, Trenton (USA): Judith Klotz (PI), Homer Wilcox (Co-PI), Helen Weiss (Coordinator). University of North Carolina, Chapel Hill (USA): Robert Millikan (PI), Nancy Thomas (Co-Investigator), Dianne Mattingly (Coordinator), Jon Player (Laboratory Technician). University of Pennsylvania, Philadelphia, PA (USA): Timothy Rebbeck (PI), Peter Kanetsky (Co-Investigator), Amy Walker (Laboratory Manager), Saarene Panossian (Laboratory Technician). Consultants: Julia Lee Taylor and Sasha Madronich, National Centre for Atmospheric Research, Boulder, Colorado (USA).

This work is supported by the National Cancer Institute, award CA112524. The MSKCC Sequenom facility is supported by the Anbinder Fund.

Abbreviations used

single nucleotide polymorphism
vitamin D receptor
Genes, environment and melanoma study
reference sequence
minor allele frequency
linkage disequilibrium
haplotype tagging SNP


Novelty and impact: We analyzed 38 common vitamin D receptor (VDR) gene polymorphisms in relation to melanoma risk in a large population-based case control study comprised of 3676 individuals with incident primary melanoma and found modest but statistically significant associations between 8 SNPs and the risk for developing subsequent new primary melanomas, including 6 SNPs investigated for the first time in relation to melanoma. These results support the hypothesis that the vitamin D pathway plays an important role in the genesis of melanomas. To our knowledge, this is the largest and most comprehensive study of VDR polymorphisms in melanoma to date.


1. Colston K, Colston MJ, Feldman D. 1,25-dihydroxyvitamin D3 and malignant melanoma: the presence of receptors and inhibition of cell growth in culture. Endocrinology. 1981;108:1083–1086. [PubMed]
2. Danielsson C, Fehsel K, Polly P, Carlberg C. Differential apoptotic response of human melanoma cells to 1 alpha,25-dihydroxyvitamin D3 and its analogues. Cell Death Differ. 1998;5:946–952. [PubMed]
3. Evans SR, Houghton AM, Schumaker L, Brenner RV, Buras RR, Davoodi F, Nauta RJ, Shabahang M. Vitamin D receptor and growth inhibition by 1,25-dihydroxyvitamin D3 in human malignant melanoma cell lines. J Surg Res. 1996;61:127–133. [PubMed]
4. Ranson M, Posen S, Mason RS. Human melanocytes as a target tissue for hormones: in vitro studies with 1 alpha-25, dihydroxyvitamin D3, alpha-melanocyte stimulating hormone, and beta-estradiol. J Invest Dermatol. 1988;91:593–598. [PubMed]
5. Yudoh K, Matsuno H, Kimura T. 1alpha,25-dihydroxyvitamin D3 inhibits in vitro invasiveness through the extracellular matrix and in vivo pulmonary metastasis of B16 mouse melanoma. J Lab Clin Med. 1999;133:120–128. [PubMed]
6. Berwick M, Armstrong BK, Ben-Porat L, Fine J, Kricker A, Eberle C, Barnhill R. Sun exposure and mortality from melanoma. J Natl Cancer Inst. 2005;97:195–199. [PubMed]
7. Moan J, Porojnicu AC, Dahlback A, Setlow RB. Addressing the health benefits and risks, involving vitamin D or skin cancer, of increased sun exposure. Proc Natl Acad Sci U S A. 2008;105:668–673. [PubMed]
8. Rosso S, Sera F, Segnan N, Zanetti R. Sun exposure prior to diagnosis is associated with improved survival in melanoma patients: results from a long-term follow-up study of Italian patients. Eur J Cancer. 2008;44:1275–1281. [PubMed]
9. Newton-Bishop JA, Beswick S, Randerson-Moor J, Chang YM, Affleck P, Elliott F, Chan M, Leake S, Karpavicius B, Haynes S, Kukalizch K, Whitaker L, et al. Serum 25-hydroxyvitamin D3 levels are associated with breslow thickness at presentation and survival from melanoma. J Clin Oncol. 2009;27:5439–5444. [PMC free article] [PubMed]
10. Stumpf WE, Sar M, Reid FA, Tanaka Y, DeLuca HF. Target cells for 1,25-dihydroxyvitamin D3 in intestinal tract, stomach, kidney, skin, pituitary, and parathyroid. Science. 1979;206:1188–1190. [PubMed]
11. Uitterlinden AG, Fang Y, Van Meurs JB, Pols HA, Van Leeuwen JP. Genetics and biology of vitamin D receptor polymorphisms. Gene. 2004;338:143–156. [PubMed]
12. Nejentsev S, Godfrey L, Snook H, Rance H, Nutland S, Walker NM, Lam AC, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Tuomilehto-Wolf E, et al. Comparative high-resolution analysis of linkage disequilibrium and tag single nucleotide polymorphisms between populations in the vitamin D receptor gene. Hum Mol Genet. 2004;13:1633–1639. [PubMed]
13. Fang Y, van Meurs JB, d'Alesio A, Jhamai M, Zhao H, Rivadeneira F, Hofman A, van Leeuwen JP, Jehan F, Pols HA, Uitterlinden AG. Promoter and 3'-untranslated-region haplotypes in the vitamin d receptor gene predispose to osteoporotic fracture: the Rotterdam study. Am J Hum Genet. 2005;77:807–823. [PubMed]
14. Hutchinson PE, Osborne JE, Lear JT, Smith AG, Bowers PW, Morris PN, Jones PW, York C, Strange RC, Fryer AA. Vitamin D receptor polymorphisms are associated with altered prognosis in patients with malignant melanoma. Clin Cancer Res. 2000;6:498–504. [PubMed]
15. Halsall JA, Osborne JE, Potter L, Pringle JH, Hutchinson PE. A novel polymorphism in the 1A promoter region of the vitamin D receptor is associated with altered susceptibilty and prognosis in malignant melanoma. Br J Cancer. 2004;91:765–770. [PMC free article] [PubMed]
16. Han J, Colditz GA, Hunter DJ. Polymorphisms in the MTHFR and VDR genes and skin cancer risk. Carcinogenesis. 2007;28:390–397. [PubMed]
17. Li C, Liu Z, Zhang Z, Strom SS, Gershenwald JE, Prieto VG, Lee JE, Ross MI, Mansfield PF, Cormier JN, Duvic M, Grimm EA, et al. Genetic variants of the vitamin D receptor gene alter risk of cutaneous melanoma. J Invest Dermatol. 2007;127:276–280. [PubMed]
18. Povey JE, Darakhshan F, Robertson K, Bisset Y, Mekky M, Rees J, Doherty V, Kavanagh G, Anderson N, Campbell H, MacKie RM, Melton DW. DNA repair gene polymorphisms and genetic predisposition to cutaneous melanoma. Carcinogenesis. 2007;28:1087–1093. [PubMed]
19. Santonocito C, Capizzi R, Concolino P, Lavieri MM, Paradisi A, Gentileschi S, Torti E, Rutella S, Rocchetti S, Di Carlo A, Di Stasio E, Ameglio F, et al. Association between cutaneous melanoma, Breslow thickness and vitamin D receptor BsmI polymorphism. Br J Dermatol. 2007;156:277–282. [PubMed]
20. Barroso E, Fernandez LP, Milne RL, Pita G, Sendagorta E, Floristan U, Feito M, Aviles JA, Martin-Gonzalez M, Arias JI, Zamora P, Blanco M, et al. Genetic analysis of the vitamin D receptor gene in two epithelial cancers: melanoma and breast cancer case-control studies. BMC Cancer. 2008;8:385. [PMC free article] [PubMed]
21. Li C, Liu Z, Wang LE, Gershenwald JE, Lee JE, Prieto VG, Duvic M, Grimm EA, Wei Q. Haplotype and genotypes of the VDR gene and cutaneous melanoma risk in non-Hispanic whites in Texas: a case-control study. Int J Cancer. 2008;122:2077–2084. [PMC free article] [PubMed]
22. Gapska P, Scott RJ, Serrano-Fernandez P, Mirecka A, Rassoud I, Gorski B, Cybulski C, Huzarski T, Byrski T, Nagay L, Maleszka R, Sulikowski M, et al. Vitamin D receptor variants and the malignant melanoma risk: a population-based study. Cancer Epidemiol. 2009;33:103–107. [PubMed]
23. Randerson-Moor JA, Taylor JC, Elliott F, Chang YM, Beswick S, Kukalizch K, Affleck P, Leake S, Haynes S, Karpavicius B, Marsden J, Gerry E, et al. Vitamin D receptor gene polymorphisms, serum 25-hydroxyvitamin D levels, and melanoma: UK case-control comparisons and a meta-analysis of published VDR data. Eur J Cancer. 2009;45:3271–3281. [PMC free article] [PubMed]
24. Arai H, Miyamoto KI, Yoshida M, Yamamoto H, Taketani Y, Morita K, Kubota M, Yoshida S, Ikeda M, Watabe F, Kanemasa Y, Takeda E. The polymorphism in the caudal-related homeodomain protein Cdx-2 binding element in the human vitamin D receptor gene. J Bone Miner Res. 2001;16:1256–1264. [PubMed]
25. Jurutka PW, Remus LS, Whitfield GK, Thompson PD, Hsieh JC, Zitzer H, Tavakkoli P, Galligan MA, Dang HT, Haussler CA, Haussler MR. The polymorphic N terminus in human vitamin D receptor isoforms influences transcriptional activity by modulating interaction with transcription factor IIB. Mol Endocrinol. 2000;14:401–420. [PubMed]
26. Colin EM, Weel AE, Uitterlinden AG, Buurman CJ, Birkenhager JC, Pols HA, Van Leeuwen JP. Consequences of vitamin D receptor gene polymorphisms for growth inhibition of cultured human peripheral blood mononuclear cells by 1, 25-dihydroxyvitamin D3. Clin Endocrinol (Oxf) 2000;52:211–216. [PubMed]
27. Mocellin S, Nitti D. Vitamin D receptor polymorphisms and the risk of cutaneous melanoma: a systematic review and meta-analysis. Cancer. 2008;113:2398–2407. [PubMed]
28. Gandini S, Raimondi S, Gnagnarella P, Dore JF, Maisonneuve P, Testori A. Vitamin D and skin cancer: a meta-analysis. Eur J Cancer. 2009;45:634–641. [PubMed]
29. Begg CB, Hummer AJ, Mujumdar U, Armstrong BK, Kricker A, Marrett LD, Millikan RC, Gruber SB, Culver HA, Zanetti R, Gallagher RP, Dwyer T, et al. A design for cancer case-control studies using only incident cases: experience with the GEM study of melanoma. Int J Epidemiol. 2006;35:756–764. [PubMed]
30. Begg CB, Berwick M. A note on the estimation of relative risks of rare genetic susceptibility markers. Cancer Epidemiol Biomarkers Prev. 1997;6:99–103. [PubMed]
31. d'Alesio A, Garabedian M, Sabatier JP, Guaydier-Souquieres G, Marcelli C, Lemacon A, Walrant-Debray O, Jehan F. Two single-nucleotide polymorphisms in the human vitamin D receptor promoter change protein-DNA complex formation and are associated with height and vitamin D status in adolescent girls. Hum Mol Genet. 2005;14:3539–3548. [PubMed]
32. Begg CB, Orlow I, Hummer AJ, Armstrong BK, Kricker A, Marrett LD, Millikan RC, Gruber SB, Anton-Culver H, Zanetti R, Gallagher RP, Dwyer T, et al. Lifetime risk of melanoma in CDKN2A mutation carriers in a population-based sample. J Natl Cancer Inst. 2005;97:1507–1515. [PubMed]
33. Nakai K, Habano W, Fujita T, Schnackenberg J, Kawazoe K, Suwabe A, Itoh C. Highly multiplexed genotyping of coronary artery disease-associated SNPs using MALDI-TOF mass spectrometry. Hum Mutat. 2002;20:133–138. [PubMed]
34. Bennett CD, Campbell MN, Cook CJ, Eyre DJ, Nay LM, Nielsen DR, Rasmussen RP, Bernard PS. The LightTyper: high-throughput genotyping using fluorescent melting curve analysis. Biotechniques. 34(6):1288–1292. 1294-2003. [PubMed]
35. Ronaghi M, Uhlen M, Nyren P. A sequencing method based on real-time pyrophosphate. Science. 1998;281(5375):363–365. [PubMed]
36. Cotignola J, Reva B, Mitra N, Ishill N, Chuai S, Patel A, Shah S, Vanderbeek G, Coit D, Busam K, Halpern A, Houghton A, Sander C, Berwick M, Orlow I. Matrix Metalloproteinase-9 (MMP-9) polymorphisms in patients with cutaneous malignant melanoma. BMC Med Genet. 2007 Mar 8;8:10. [PMC free article] [PubMed]
37. Abramovich F, Benjamini Y. False Discovery Rate. In: Balakrishnan N, editor. Encyclopedia of Statistical Sciences. Vol. 4. New Jersey: Wiley & Sons, Inc.; 2006. pp. 2240–2243.
38. Wilk MB, Gnanadesikan R. Probability plotting methods for the analysis of data. Biometrika. 1968 Mar;55(1):1–17. [PubMed]
39. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–560. [PMC free article] [PubMed]
40. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50:1088–1101. [PubMed]
41. Rhead B, Karolchik D, Kuhn RM, Hinrichs AS, Zweig AS, Fujita PA, Diekhans M, Smith KE, Rosenbloom KR, Raney BJ, Pohl A, Pheasant M, et al. The UCSC Genome Browser database: update 2010. Nucleic Acids Res. 2010;38:D613–D619. [PMC free article] [PubMed]
42. Selbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63. [PubMed]
43. Berwick M, Orlow I, Hummer AJ, Armstrong BK, Kricker A, Marrett LD, Millikan RC, Gruber SB, Anton-Culver H, Zanetti R, Gallagher RP, Dwyer T, et al. The prevalence of CDKN2A germ-line mutations and relative risk for cutaneous malignant melanoma: an international population-based study. Cancer Epidemiol Biomarkers Prev. 2006;15:1520–1525. [PubMed]
44. Martin RJ, McKnight AJ, Patterson CC, Sadlier DM, Maxwell AP. A rare haplotype of the vitamin D receptor gene is protective against diabetic nephropathy. Nephrol Dial Transplant. 2010;25:497–503. [PubMed]
45. Halsall JA, Osborne JE, Epstein MP, Pringle JH, Hutchinson PE. The unfavorable effect of the A allele of the vitamin D receptor promoter polymorphism A-1012G has different mechanisms related to susceptibility and outcome of malignant melanoma. Dermatoendocrinol. 2009;1:54–57. [PMC free article] [PubMed]
46. Macgregor S, Hottenga JJ, Lind PA, Suchiman HE, Willemsen G, Slagboom PE, Montgomery GW, Martin NG, Visscher PM, Boomsma DI. Vitamin D receptor gene polymorphisms have negligible effect on human height. Twin Res Hum Genet. 2008;11:488–494. [PubMed]
47. Rengarajan J, Szabo SJ, Glimcher LH. Transcriptional regulation of Th1/Th2 polarization. Immunol Today. 2000;21:479–483. [PubMed]
48. Yamamoto H, Miyamoto K, Li B, Taketani Y, Kitano M, Inoue Y, Morita K, Pike JW, Takeda E. The caudal-related homeodomain protein Cdx-2 regulates vitamin D receptor gene expression in the small intestine. J Bone Miner Res. 1999;14:240–247. [PubMed]
49. Rozek LS, Lipkin SM, Fearon ER, Hanash S, Giordano TJ, Greenson JK, Kuick R, Misek DE, Taylor JM, Douglas JA, Rennert G, Gruber SB. CDX2 polymorphisms, RNA expression, and risk of colorectal cancer. Cancer Res. 2005;65:5488–5492. [PubMed]
50. Sertznig P, Dunlop T, Seifert M, Tilgen W, Reichrath J. Cross-talk between vitamin D receptor (VDR)- and peroxisome proliferator-activated receptor (PPAR)-signaling in melanoma cells. Anticancer Res. 2009;29:3647–3658. [PubMed]