|Home | About | Journals | Submit | Contact Us | Français|
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.
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 1–5, 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 6–9. 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 14–23. 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 24–26 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.
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.
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.
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.
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: http://genome.ucsc.edu/cgi-bin/hgTables 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.
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).
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.
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.
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.
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.
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,21–23. 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,27–28, 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.
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.
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.