This article has Open Peer Review reports available.
Serum 25-hydroxyvitamin D3 levels and vitamin D receptor variants in melanoma patients from the Mediterranean area of Barcelona
© Ogbah et al.; licensee BioMed Central Ltd. 2013
Received: 23 July 2012
Accepted: 12 February 2013
Published: 16 February 2013
Serum 25-hydroxyvitamin D3 (Vitamin D) insufficiency and single-nucleotide polymorphisms (SNPs) on its receptor, Vitamin D receptor (VDR), have been reported to be involved in melanoma susceptibility in populations mostly from northern countries.
To investigate 25-hydroxyvitamin D3 levels and VDR SNPs in melanoma patients from sunny area of Barcelona, two studies were carried out. The first study evaluated the levels of Vitamin D at time of melanoma diagnosis and the second one analyzed the association between VDR genetic variants and risk of having a high nevus number, the strongest phenotypic risk factor for melanoma.
The levels of 25-hydroxyvitamin D3 in 81 melanoma patients at diagnosis were measured. In a second group of melanoma patients, including 150 with low and 113 with high nevus number, 11 VDR SNPs were analyzed for their association with nevus number.
In the first study, 68% of patients had insufficient levels of 25-hydroxyvitamin D3 (<25 ng/ml). Autumn-winter months and fair phototype were associated with 25-hydroxyvitamin D3 insufficiency; after multivariate analysis, season of sampling remained the only independent predictor of 25-hydroxyvitamin D3 levels. In the second study, VDR variant rs2189480 (P = 0.006) was associated with risk of high nevus number whereas rs2239179 (P = 0.044) and rs7975128 (P = 0.0005) were protective against high nevus number. After Bonferroni adjustment only rs7975128 remained significant. In stratified analysis, SNP rs7975128 was found protective against multiple melanomas (P = 0.021).
This study showed that even in Barcelona, a sunny Mediterranean area, 25-hydroxyvitamin D3 levels were sub-optimal in the majority of melanoma patients at diagnosis. The involvement of VDR in nevi and, in turn, in melanoma susceptibility has also been suggested. Larger studies are needed to confirm our findings.
Serum 25-hydroxyvitamin D3 (Vitamin D) is a hormone classically known for its regulation of bone and calcium metabolism. In recent years interest in 25-hydroxyvitamin D3 has grown as evidence has shown it to play a role in cancer susceptibility and outcome, postulated to be related to its anti-proliferative pro-apoptotic and pro-differentiation properties [1, 2]. Suboptimal levels of 25-hydroxyvitamin D3 have been shown to be associated with increased risk of several cancers including melanoma [3, 4]. Patients diagnosed with melanoma had suboptimal levels of 25-hydroxyvitamin D3 in the North of England within a few months of diagnosis, especially in winter months  and a reduced level of 25-hydroxyvitamin D3 was associated with thicker melanomas, poorer outcome  and aggressive melanoma . There are no reports of 25-hydroxyvitamin D3 levels in melanoma patients living in sunny Mediterranean countries and we therefore herein report a small study in melanoma patients from Barcelona.
Sunburn and intermittent sun exposure are the most important environmental risk factors for melanoma  but sun exposure is also the main source for 25-hydroxyvitamin D3 photosynthesis. 25-hydroxyvitamin D3 derives from cholesterol metabolites in skin after exposure to ultraviolet B radiation (UVB) . The skin exposure to UVB induces the photolytic conversion of 7-dehydrocholesterol (7-DHC) to previtamin D and from there to vitamin D3. The metabolic activation of vitamin D3 requires two hydroxylation steps . The first hydroxylation occurs in the liver where vitamin D3 is converted in 25-hydroxyvitamin D3, the major circulating metabolite of vitamin D3. The second hydroxylation transformed in the kidney the 25-hydroxyvitamin D3 in 1,25-dihydroxyvitamin D3. While circulating 1,25-dihydroxyvitamin D3 levels mainly reflect the synthesis of 1,25-dihydroxyvitamin D3 in the kidney, 25-hydroxyvitamin D3 serum levels are a good parameter of whole body vitamin D status and for this reason is the metabolite most measured for the classification of vitamin D status . In vitro, 25-hydroxyvitamin D3 has been shown to inhibit the growth of normal melanocytes and melanocytes derived from melanoma [11, 12]. The 25-hydroxyvitamin D3 mediates its biological effects by binding to its specific nuclear receptor, Vitamin D Receptor (VDR).
The VDR gene, located on chromosome 12q13.11, has more than 200 single-nucleotide polymorphisms (SNPs)  but only a few SNPs, considered to be functional , have been studied for their potential associations with melanoma susceptibility [16–21]. Meta-analysis of all reports published by 2009 supported the association between VDR functional SNPs BsmI (rs1544410) and FokI (rs2228570) and melanoma risk [22, 23]. More recently a large SNP study showed that common inherited variants at VDR influence melanoma risk , providing further evidence supporting a role for 25-hydroxyvitamin D3 in melanoma susceptibility.
Common acquired nevi are melanocytes in proliferation and high nevus number represents the strongest phenotypic risk factor for melanoma development [25, 26]. Genome-wide association studies have identified common susceptibility genes for melanoma and increased nevus number [27, 28].
If VDR SNPs are associated with melanoma risk therefore, elucidating the relationship between VDR variants and nevi is an important step towards understanding the involvement of VDR and, indirectly, of 25-hydroxyvitamin D3, in melanoma carcinogenesis.
The association between 25-hydroxyvitamin D3 levels and/or inherited variation in VDR and melanoma risk has been studied in few populations, mostly from northern countries such as England, USA, Poland and Germany [4, 6, 21–23]. Only two reports on VDR variants have been published from southern European countries, specifically Italy and Spain [18, 29].
We presented the results of two studies performed in Barcelona melanoma patients. The first study evaluated the levels of 25-hydroxyvitamin D3 at time of melanoma diagnosis and the second one analyzed the association between VDR genetic variants and risk of having a high nevus number.
Two cohort studies both including Caucasian melanoma patients were carried out. The two studies were mutually exclusive (each patient is present in only one study) and included groups of patients belonging to two different studies. In the first, a retrospective cohort study, 25-hydroxyvitamin D3 levels were measured in serum samples stored from melanoma patients at diagnosis for the investigation of biomarkers. In the second, the association between VDR variants and nevus number was evaluated in a second cohort of melanoma patients recruited to a genetic epidemiological study of nevi. Patients were selected according to the total nevus number at the time of melanoma diagnosis and classified into two groups: 150 cases with a low nevus number (<50 nevi) and 113 cases with a high nevus number (>100 nevi). We selected two groups of patients where there would be no possible overlap between them, as <50 nevi were closer to the normal nevus number in our population  and >100 nevi being a well- documented risk factor for melanoma .
Both cohort studies included melanoma patients from the urban area of Barcelona recruited from the Melanoma Unit of Hospital Clinic of Barcelona. This is a referral centre for Melanoma in the Catalonia region.
Both studies were approved by the Ethics Committee of the Hospital Clinic of Barcelona and written informed consent was obtained from all participants in the study.
The first set of melanoma patients (25-hydroxyvitamin D3 study)
The patients had all been diagnosed with melanoma between 2004–2008. Patients were those recruited to the epidemiological study from which serum samples had been cryopreserved in the Biobank of Hospital Clínic and where the blood was collected within a maximum of 3 months after melanoma diagnosis. All patients had localized tumors (stage I and stage II) except one (stage IIIa) which had lymph node micrometastasis. Data on gender, age at melanoma diagnosis, nevus number, Fitzpatrick skin type, eye color, hair color, examined actinic damage, number of melanoma, Breslow thickness, reported solar exposure before 10 years of age, between 10–18 and over 18 years were collected through a questionnaire or extracted by medical records. Participants were asked about weight, height and habits of photoprotection via an interview.
The second set of melanoma patients (VDRpolymorphisms study)
A total of 263 melanoma patients were recruited between 2004–2007. The patients were selected according to the total number of nevi at the time of melanoma diagnosis and classified into two groups: 150 cases with a low nevus number (cases with <50) and 113 cases with a high nevus number (cases with >100). Total body nevus number has been performed by dermatologists and trained nurses. Data on gender, age at melanoma diagnosis, nevus number, Fitzpatrick skin type, eye color, hair color, examined actinic damage, number of melanoma, reported solar exposure as before 10 years of age, between 10–18 and over 18 years were collected through a questionnaire or extracted by medical records.
25-hydroxyvitamin D3 measurement
The levels of 25-hydroxyvitamin D3 (ng/ml) were analyzed by a chemiluminescence immunoassay method using the Liaison ® kit (DiaSorin, Minnesota, USA) as previously described . This method has a detection limit range from 4 ng/ml to 150 ng/ml. Insufficient level of 25-hydroxyvitamin D3 was defined as less than 25 ng/ml (60 nmol/L).
In total, 11 VDR SNPs rs7136534, rs11574027, 11168287, rs2238136, rs3782905, rs2189480, rs2239179, rs11574077, rs11168267 and rs7975128 were analyzed. All SNPs were selected as having minor allele frequencies (MAF) ≥ 0.05 in HapMap (http://www.hapmap.org) European individuals and at least one independent validation criterion, as established in dbSNP (http://www.ncbi.nlm.nih.gov/snp). The SNPs were selected using data from the HapMap project and the Tagger tool  with linkage disequilibrium r2 > 0.5.
Genomic DNA was isolated from peripheral blood samples using the salting-out method. Genotyping was performed with GenomeLab SNPstream genotyping platform (Beckman & Coulter Inc. Fullerton, CA) and its accompanying SNPstream software, according to the manufacturer’s protocols. Primers were designed using Autoprimer (http://www.autoprimer.com). The design of the primers failed in the presence of adjacent SNPs or of repeat sequence which can interfere with primers annealing. If a selected SNP failed assay primer design, if possible, an alternative tagging SNP in complete LD was chosen. Two SNPs, FokI (rs2228570) and BsmI (rs1544410), failed primer design. By HapMap, BsmI was replaced with rs7975128 as in strong LD (r2 = 1) in the Caucasian population while no other variants in LD were found for FokI. Samples and SNPs with call rates below 90% were excluded from the analysis. Failed genotypes were not repeated. As quality controls, 10% of samples were added as duplicates so that concordance between genotype calls could be assessed. The genotyping assays were performed at the LIMM (Leeds) genotyping facilities.
Quantitative variables were described using median and interquartile range (percentiles 25th and 75th). Inferential analyses for these variables were performed by U Mann–Whitney or Wilcoxon test for unpaired data. Qualitative variables were tabulated by absolute frequencies and percentages and evaluated using χ 2 test or Fisher’s exact test where appropriate.
The analysis was performed using SPSS v15.0 (SPSS Inc., Chicago IL). The level of statistical significance was set at P < 0.05. Logistic regression was used to estimate Odds Ratio (OR) and 95% Confidence Interval (CI) for association between 25-hydroxyvitamin D3 levels (<25 ng/ml) and gender, age diagnosis of melanoma (<55/>55), BMI (normal weight/overweight/obese), Fitzpatrick skin type (dark/fair), eye color (dark/light), hair color (dark/light), actinic damage (no/yes); habits of photoprotection (usually/never), solar exposure (none to mild, that is under 50 hour per year (h/y) or moderate to severe that is more than 50 (h/y)) before 10 years of age, between 10–18 and over 18 years, number of melanoma (single/multiple) and Breslow (<0.75 mm/0.75-1 mm/>1 mm). For the purpose of the analysis, Breslow thickness was also analyzed as a continuous variable. The lineal trend of association between 25-hydroxyvitamin D3 levels and Breslow thickness was evaluated by Spearman correlation. In this cohort, cases with an intermediate nevus number between 51 and 99 (51–99) were present; nevus number variable included three groups (≤50 nevi), (51–99) and (≥100 nevi).
To test the effect of the season on blood sampling as an independent factor, a multivariate approach was performed adjusting by phototype and gender.
Statistical analyses were performed using SNPassoc package in R (version 2.10.1) . Departure from Hardy–Weinberg equilibrium (HWE) for each SNP was assessed by χ 2 test. Being a study including only cases, SNPs found outside HWE were analyzed because HWE could be informative for SNPs that confer risk rather than genotyping error. To measure the association between VDR variant alleles and risk for increased number of nevi, cases with less than 50 nevi were compared with cases with more than 100 nevi, where cases <50 were considered as the referent group.
The association between genotype and nevus number was estimated using logistic regression and analyzed under co-dominant, dominant, recessive and log-additive models. Significance was defined as P < 0.05. The best inheritance model was selected using the Akaike criteria. The reported ORs and 95% CI were adjusted for age and gender. Significant associations between SNPs and nevus number were also tested by the Bonferroni test. To control confounder variables, regression analysis was repeated adjusting for whether there was a single or multiple primaries. Logistic regression was also used to estimate the association between polymorphisms and clinical characteristics such as eye color (dark/light), hair color (dark/light), Fitzpatrick skin phototype (dark/fair), actinic damage (no/yes) and number of melanomas (single/multiple). Associations between VDR SNPs and solar exposure (none-mild < 50 h/y/severe-moderate > 50 h/y) before 10, between 10–18 and over 18 years was evaluated by interaction analysis.
Descriptive analysis of characteristics and medians of 25-hydroxyvitamin D3 levels at time of diagnosis in melanoma patients
Vitamin D levels ng/ml
Age at diagnosis of melanoma (years)
Fitzpatrick skin type
20.2 (17.8- 24.3)
Habits of solar exposure and photoprotection
Solar exposure before 10 years old
Solar exposure 10–18 years old
Solar exposure over 18 years old
Photoprotection before 10 years old
Photoprotection 10–18 years old
Photoprotection over 18 years old
Number of Melanoma
> = 1
It is known that BMI influences 25-hydroxyvitamin D3 levels. The 25-hydroxyvitamin D3 levels were shown to be lower (median 15.7 ng/ml) in obese patients with higher BMI (BMI ≥ 30) than in overweight patients with BMI between 25–29.9 (median 19.6 ng/ml) and in normal weight patients with BMI ≤ 24.9 (median 22.9 ng/ml) (P = 0.06).
Patients with phototype I/II had lower levels of 25-hydroxyvitamin D3 (median 18 ng/ml) than patients with phototype III/IV (median 22.9 ng/ml) (P = 0.005) in univariate analysis. Those patients who reported moderate/severe solar exposure before 10 years of age (≥ 50 h/y) had higher levels of 25-hydroxyvitamin D3 (median 22.5 ng/ml) than those with none/mild solar exposure (< 50 h/y) (median 16.8 ng/ml) (P = 0.019).
A higher percentage of patients with fair skin had none/mild solar exposure (<50 h/y) than dark-skinned patients in age groups <10 years (35% vs 12%; P = 0.06) (unpublished observations). The use of sunscreen was reported in 70% of patients. Skin type was also related with sunscreen: 62% of patients with light skin phototype usually used sunscreen compared with 29% in patients with darker skin phototype (P = 0.023) (unpublished observations).
As regards nevus number, patients with >100 showed higher 25-hydroxyvitamin D3 levels (26.3 ng/ml) than patients with 51–99 (20.2 ng/ml) or <50 nevi (21.2 ng/ml) but no significant differences were found among the three groups (P = 0.266).
Most of our patients (74%) had melanomas with Breslow thickness ≤ 2 mm (range from 0.2 to 7.5 mm).
Breslow was analyzed both as continuous and categorical variable. No statistically significant differences were found between Breslow thickness and 25-hydroxyvitamin D3 levels (r2 = 0.009, P = 0.42). The same result was found when examining the association between Vitamin D levels and Breslow as a categorical variable (P = 0.544).
No associations in levels of 25-hydroxyvitamin D3 were detected for actinic damage (P = 0.161), eye color (P = 0.777), hair color (P = 0.912), photoprotection habits before 10 (P = 0.569), between 10–18 (P = 0.672), over 18 years old (P = 0.824) and solar exposure between 10–18 (P = 0.08), over 18 years old (P = 0.267) and number of melanomas (P = 0.992) (Table 1).
Association between clinical characteristics and suboptimal 25-hydroxyvitamin D3 levels at time of diagnosis in melanoma patients
VitD level sufficient >25 ng/mL n = 26(32%)
VitD level insufficient ≤25 ng/mL n = 55(68%)
Age at diagnosis of melanoma**
0.99 (0.96- 1.01)
1.9 (0.74- 4.88)
BMI basal (kg/m2)
0.56 (0.16- 1.92)
0.61 (0.16- 2.49)
Season blood sampling
3.6 (1.32- 9.81)
Fitzpatrick skin type
2.59 (0.99- 6.76)
0.89 (0.32- 2.49)
0.67 (0.26- 1.76)
0.89 (0.27- 2.89)
0.97 (0.34- 2.73)
Habits of solar exposure and photoprotection
Solar exposure before 10 years old
1.55 (0.49- 4.9)
Solar exposure 10–18 years old
3.83 (1.01- 14.5)
Solar exposure > 18 years old
6.11 (0.68- 54.5)
Photoprotection before 10 years old
1.20 (0.43- 3.37).33)
Photoprotection 10–18 years old
1.62 (0.46- 5.73)
Photoprotection over 18 years old
1.79 (0.66- 4.80)
Number of Melanoma
0.68 (0.13- 3.63)
0.89 (0.22- 3.54)
> = 1
1.4 (0.43- 4.48)
0.8 (0.52- 1.3)
0.8 (0.52- 1.2)
1.1 (0.78- 1.54)
The relationship between 25-hydroxyvitamin D3 levels, as a dichotomous variable, and clinical characteristics was showed in Table 2. The group with insufficient 25-hydroxyvitamin D3 level (≤25 ng/ml) differed from those with appropriate 25-hydroxyvitamin D3 (>25 ng/ml) in season of blood sampling (P = 0.012), solar exposure between 10–18 years (P = 0.049) and skin type (P = 0.05) (Table 2). Since patient recollection of solar exposure may be subject to response and recall bias, solar exposure variable was not included in multivariate analysis.
Multivariate models for prediction of presence of low levels of 25-Hydroxyvitamin D3
Season blood sampling (October- March)
Skin phototype (Fair)
2.21 (0.81- 6.04)
Season blood sampling (October- March)
3.3 (1.17- 9.27)
Skin phototype (Fair)
2.33 (0.84- 6.46)
1.88 (0.68- 5.21)
Clinical characteristics of patients with low (<50) and high (>100) nevus number
Cases <50(%) n = 150
Cases >100(%) n = 113
Age at melanoma diagnosis*
Fitzpatrick skin type
Habits of solar exposure
Solar exposure before 10 years old
Solar exposure 10–18 years old
Solar exposure over 18 years old
Number of melanoma
A total of 10 out of 11 SNPs were successfully genotyped in 263 melanoma patients. SNP rs739837 (BglI) failed genotype analysis. All the 10 SNPs followed HWE in the reference group (cases with <50) (Table 4). Genotyping quality was confirmed with 100% concordance between duplicates.
Association between VDR polymorphisms and nevus number
P additive a
P additive b
1.15 (0.65- 2.02)
0.78 (0.21- 2.94)
1.03 (0.65- 1.62)
1.18 (0.73- 1.92)
1.19 (0.64- 2.21)
1.30 (0.67- 2.52)
0.99 (0.55- 1.76)
0.58 (0.27- 1.24)
0.79 (0.55- 1.14)
0.83 (0.56- 1.23)
1.61 (0.93- 2.79)
0.47 (0.12- 1.82)
1.11 (0.72- 1.72)
1.05 (0.66- 1.67)
1.07 (0.62- 1.83)
0.22 (0.07- 0.71)
0.70 (0.46- 1.06)
0.69 (0.44- 1.08)
0.21 (0.07- 0.66)
1.92 (1.09- 3.39)
2.97 (1.12- 7.87)
1.80 (1.17- 2.76)
1.66 (1.05- 2.60)
0.70 (0.39- 1.26)
0.44 (0.20- 0.99)
0.67 (0.46- 0.99)
0.67 (0.44- 1.02)
0.54 (0.21- 1.38)
0.52 (0.21- 1.27)
0.45 (0.17- 1.18)
0.68 (0.32- 1.48)
0.84 (0.41- 1.72)
0.84 (0.38- 1.84)
0.47 (0.27- 0.82)
0.29 (0.12- 0.69)
0.51 (0.35- 0.76)
0.58 (0.38- 0.87)
Sub-group analysis revealed that rs7975128 was associated with decreased risk of multiple primary melanoma (OR = 0.64 95%CI (0.43-0.94); P additive = 0.021) and with light eye color (OR = 1.48 95%CI (1.00-2.19); P additive = 0.049).
Our work is the first study to describe 25-hydroxyvitamin D3 insufficiency in melanoma patients from Spain. Furthermore, the association between VDR SNPs and nevus number, an established risk factor for melanoma, was also evaluated.
The optimal level of 25-hydroxyvitamin D3 remains contentious4. Several meta-analyses support the evidence that low 25-hydroxyvitamin D3 serum level, which is highly prevalent worldwide, is a risk factor for many chronic diseases including cancer . The minimum desirable serum level of 25-hydroxyvitamin D3 has been suggested to be between 20 ng/ml and 30 ng/ml . In 2011 The IOM committee reported guidelines stating that a level of 25-hydroxyvitamin D3 >20 ng/ml is needed for good bone and general health for practically all individuals . However, the IOM report did not consider clinical and demographic variables (ie: race/ethnicity, adiposity, body composition, sun exposure, etc.) by which 25-hydroxyvitamin D3 levels are notoriously affected. Recently, a meta-analysis involving 59 231 healthy subjects from eleven prospective cohort studies suggested that the optimal 25-hydroxyvitamin D3 concentrations of between 30 ng/ml and 40 ng/ml (75 nmol/L and 100 nmol/L) are needed to reduce mortality . The majority of authors site the cut off as <30 ng/ml as this level is reported to correlate with leveling off of the parathyroid hormone (PTH) level. The level chosen as cut off varied in studies in European populations . In our study a conservative level of ≤25 ng/ml as 25-hydroxyvitamin D3 insufficiency was chosen as reported in a previous study in melanoma patients . Using this definition, most of our patients (68%) at the time of melanoma diagnosis had suboptimal levels of 25-hydroxyvitamin D3. Newton-Bishop et al. determined the levels of 25-hydroxyvitamin D3 in 1132 patients with melanoma from the North of England (Leeds), where insufficient levels of 25-hydroxyvitamin D3 was found in 64% of patients, a proportion comparable to our study. In the same study, 25-hydroxyvitamin D3 levels were associated with Breslow thickness. In contrast to the UK study, we did not found any association between 25-hydroxyvitamin D3 and Breslow thickness as the majority of our melanoma patients has Breslow thickness <2 mm.
The prevalence of insufficiency in Spain, a Mediterranean country with many hours of sunshine and where tanning is considered indicative of a healthy and attractive status, was therefore unexpected. 25-hydroxyvitamin D3 levels in general Spanish population have been reported to be insufficiency by studies tended to be small and in different sub-groups of the population [40–42]. The largest study analyzed 25-hydroxyvitamin D3 level in 1262 individuals from Spain and reported 25-hydroxyvitamin D3 insufficiency in one-third (33%) of the population .
Consistent with previous finding , there was seasonal variation in serum levels of 25-hydroxyvitamin D3 where 25-hydroxyvitamin D3 levels were lower in autumn-winter months.
The role of solar exposure in melanoma is complex as, on the one hand, it is a strong risk factor for melanoma, and, on the other, lack of solar exposure leads to 25-hydroxyvitamin D3 deficiency. It has been recently reported that regular weekend solar exposure is protective against melanoma  and is associated with increased 25-hydroxyvitamin D3 levels . We have observed a correlation between less solar exposure during adolescence and insufficient 25-hydroxyvitamin D3 levels. However this result has to be interpreted with caution as self- reported past solar exposure is subject to recall bias.
In our study, 25-hydroxyvitamin D3 levels at the time of diagnosis were lower in fair-skinned patients than darker ones. Although this relationship was no longer statistically significant when corrected for season, this may reflect low statistical power. However, three other studies have in fact also reported that within white-skinned populations, fair skinned people have lower levels of 25-hydroxyvitamin D3 [45–47].
It was postulated that lower levels of 25-hydroxyvitamin D3 in individuals with fair skin reflect their habits of sun protection and in fact, this group of patients did report fewer hours of daily sun exposure and greater use of sunscreen than patients with darker skin type (data not shown). Our results support the hypothesis that increased photoprotection practices may contribute to increase 25-hydroxyvitamin D3 deficiency although the results are still contradictory [48–50].
The effects of 25-hydroxyvitamin D3 are mediated by VDR and SNPs in this gene are expected to influence the anti-proliferative activity of 25-hydroxyvitamin D3 in melanocytes. The potential role of functional VDR variants such as FokI and BsmI in melanoma risk has been assessed by meta-analysis studies [22, 23]. In our study, significant associations have been detected between nevus number, the most potent phenotypic risk factor of melanoma and 3 out of 11 VDR SNPs. The SNP rs2189480 was associated with risk of having a high number of nevi whereas SNPs rs7975128 and rs2239179 with decreased risk of having a high number of nevi. All 3 SNPs are localized in the VDR coding region between SNPs FokI and BsmI and have been previously described in other pathologies [50–52]. After Bonferroni correction, only SNP rs7975128 remained significant.
When we analyzed the effect of VDR variant and melanoma susceptibility, we found that SNP rs7975218 was also associated with decreased risk of multiple primary melanomas. Recently, in a large case-only study, 8 VDR SNPs, 6 not previously studied, were associated with melanoma susceptibility . As Orlow et al. report, the present work confirmed the efficiency of the tagSNP method and showed that information about the genetic variability of VDR is still incomplete.
According to HapMap, rs7975218 is in complete LD (r2 = 1) with VDR functional variants BsmI and TaqI (rs731236) which formed a block in tight LD, including also ApaI (rs7975232), located at the 3’UTR region of VDR. Consistent with previous findings in melanoma where, both TaqI and BsmI have been shown to be protective against melanoma risk [21–23], SNP rs7975128 showed a protective tendency for high nevus number and number of melanoma.
By functional study, Carling et al. have demonstrated that baT (haplotype of BsmI/ApaI/TaqI) is associated with low VDR mRNA expression. The 3’UTR variants are thought to affect VDR mRNA stability  but the functionality of these variants is still unclear .
We also reported an association between SNP rs7975128 and light eye color, which, with pigmentary features such as light hair and fair skin, is considered a melanoma risk factor related to sun exposure. Although the significance of this association is not clear, the presence of this SNP among patients with light eye color could suggest the involvement of the VDR gene in melanoma predisposition among a subgroup of patients at risk of low 25-hydroxyvitamin D3. Relationship between VDR genetic variants and phenotype features among melanoma patients has been previously described .
The only report in a Spanish population found significant associations between VDR SNPs and clinical characteristics such as fair skin, absence of childhood sunburns and tumor in head-neck and trunk but not with melanoma risk, suggesting that VDR may modulate melanoma susceptibility .
Studies including only melanoma patients have been shown to be valid in identifying rare genetic variants associated with melanoma predisposition as rare variants have a higher prevalence among cases [24, 54].
The limitations of our study should be considered. The small sample size of the two groups of patients may affect our results; however the proportion of patients with insufficient 25-hydroxyvitamin D3 is comparable to a previous study in a larger melanoma population . Information regarding dietary and 25-hydroxyvitamin D3 supplement sources, conditions that may affect 25-hydroxyvitamin D3 levels, were not taken into account in the analysis as not available.
Information on both serum 25-hydroxyvitamin D3 levels and VDR variants are required to understand the role of 25-hydroxyvitamin D3 in melanoma risk. High prevalence of sub-optimal levels of 25-hydroxyvitamin D3 was found among melanoma patients at diagnosis. The relationship between VDR SNP rs7975128 and, respectively, low nevus number and less multiple primary melanoma, suggests the involvement of VDR, and in turn of 25-hydroxyvitamin D3, in melanoma susceptibility.
Our results provide further evidence that 25-hydroxyvitamin D3 could be involved in melanoma etiology in a population from Barcelona, a Mediterranean area characterized by a sunny climate. However this finding needs to be validated in a larger sample.
The Barcelona team’s research has been supported by Fondo de Investigaciones Sanitarias (grants 05/0302, 06/0265 and 09/1393) and by the European Commission under the 6th Framework Programme, Contract nr: LSHC-CT-2006-018702 (GenoMEL) and by the National Cancer Institute (NCI) of the US National Institute of Health (NIH) (CA83115). Zighereda Ogbah was supported by financial support provided by Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR).We also thank the principal investigators of the GenoMEL founding grants David Elder, Julia Newton Bishop and Nelleke Gruis for their support and the leader of the exchange Programme Giovanna Bianchi-Scarrá and also Linda Battistuzzi. Thanks to research nurses from the Melanoma Unit of Hospital Clínic of Barcelona, Daniel Gabriel and Pablo Iglesias for collecting phenotype data of patients; to Michael Churchman from LIMM (Leeds) for genotyping support; and to Helena Kruyer for helping with the English editing and correction of the manuscript.
This work was developed in the Melanoma Units of Hospital Clinic I Provincial de Barcelona, IDIBAPS; the Bioinformatics and Biostatistics Unit, and Translational Research Laboratory, Catalan Institute of Oncology, Bellvitge Biomedical Research Institute (IDIBELL), L’Hospitalet, Barcelona, Spain; and in Epidemiology and Biostatistics, LIMM, University of Leeds.
- Garland CF, Garland FC, Gorham ED, Lipkin M, Newmark H, Mohr SB, Holick MF: The role of vitamin D in cancer prevention. Am J Public Health. 2006, 96: 252-261. 10.2105/AJPH.2004.045260.View ArticlePubMedPubMed CentralGoogle Scholar
- Borradale D, Kimlin M: Vitamin D in health and disease: an insight into traditional functions and new roles for the ’sunshine vitamin’. Nutr Res Rev. 2009, 22: 118-136.View ArticlePubMedGoogle Scholar
- Holick MF: Vitamin D deficiency in 2010: health benefits of vitamin D and sunlight: a D-bate. Nat Rev Endocrinol. 2011, 7: 73-75. 10.1038/nrendo.2010.234.View ArticlePubMedGoogle Scholar
- Field S, Newton-Bishop JA: Melanoma and vitamin D. Mol Oncol. 2011, 5: 197-214. 10.1016/j.molonc.2011.01.007.View ArticlePubMedGoogle Scholar
- 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, Jackson S, Gerry E, Nolan C, Bertram C, Marsden J, Elder DE, Barrett JH, Bishop DT: Serum 25-hydroxyvitamin D3 levels are associated with breslow thickness at presentation and survival from melanoma. J Clin Oncol. 2009, 27: 5439-5444. 10.1200/JCO.2009.22.1135.View ArticlePubMedPubMed CentralGoogle Scholar
- Nurnberg B, Graber S, Gartner B, Geisel J, Pfohler C, Schadendorf D, Tilgen W, Reichrath J: Reduced serum 25-hydroxyvitamin D levels in stage IV melanoma patients. Anticancer Res. 2009, 29: 3669-3674.PubMedGoogle Scholar
- Gandini S, Sera F, Cattaruzza MS, Pasquini P, Picconi O, Boyle P, Melchi CF: Meta-analysis of risk factors for cutaneous melanoma: II Sun exposure. Eur J Cancer. 2005, 41: 45-60. 10.1016/j.ejca.2004.10.016.View ArticlePubMedGoogle Scholar
- Holick MF: Sunlight and vitamin D for bone health and prevention of autoimmune diseases, cancers, and cardiovascular disease. Am J Clin Nutr. 2004, 80 (Suppl 6): 1678S-1688S.PubMedGoogle Scholar
- Holick MF: Vitamin D, deficiency. N Engl J Med. 2007, 357: 266-281. 10.1056/NEJMra070553.View ArticlePubMedGoogle Scholar
- Dusso AS, Brown AJ, Slatopolsky E: Vitamin D. Am J Physiol Renal Physiol. 2005, 289: F8-F28. 10.1152/ajprenal.00336.2004.View ArticlePubMedGoogle Scholar
- Seifert M, Rech M, Meineke V, Tilgen W, Reichrath J: Differential biological effects of 1,25-dihydroxyVitamin D3 on melanoma cell lines in vitro. J Steroid Biochem Mol Biol. 2004, 89–90: 375-379.View ArticlePubMedGoogle Scholar
- 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. 10.1111/1523-1747.ep12477126.View ArticlePubMedGoogle Scholar
- Whitfield GK, Hsieh JC, Jurutka PW, Selznick SH, Haussler CA, MacDonald PN, Haussler MR: Genomic actions of 1,25-dihydroxyvitamin D3. Journal Nutr. 1995, 125 (Suppl 6): 1690S-1694S.Google Scholar
- Nejentsev S, Godfrey L, Snook H, Rance H, Nutland S, Walker NM, Lam AC, Guja C, Ionescu-Tirgoviste C, Undlien DE, Rønningen KS, Tuomilehto-Wolf E, Tuomilehto J, Newport MJ, Clayton DG, Todd JA: 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. 10.1093/hmg/ddh169.View ArticlePubMedGoogle Scholar
- 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. 10.1016/j.gene.2004.05.014.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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, Zuppi C, Capoluongo E: Association between cutaneous melanoma, breslow thickness and vitamin D receptor BsmI polymorphism. Br J Dermatol. 2007, 156: 277-282. 10.1111/j.1365-2133.2006.07620.x.View ArticlePubMedGoogle Scholar
- Gapska P, Scott RJ, Serrano-Fernandez P, Mirecka A, Rassoud I, Górski B, Cybulski C, Huzarski T, Byrski T, Nagay L, Maleszka R, Sulikowski M, Lubinski J, Debniak T: Vitamin D receptor variants and the malignant melanoma risk: a population-based study. Cancer Epidemiol. 2009, 33: 103-107. 10.1016/j.canep.2009.06.006.View ArticlePubMedGoogle Scholar
- 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. 10.1002/ijc.23357.View ArticlePubMedPubMed CentralGoogle Scholar
- 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, Bale L, Bertram C, Field H, Barth JH, Silva Idos S, Swerdlow A, Kanetsky PA, Barrett JH, Bishop DT, Bishop JA: 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. 10.1016/j.ejca.2009.06.011.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1002/cncr.23867.View ArticlePubMedGoogle Scholar
- 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. 10.1016/j.ejca.2008.10.003.View ArticlePubMedGoogle Scholar
- Orlow I, Roy P, Reiner AS, Yoo S, Patel H, Paine S, Armstrong BK, Kricker A, Marrett LD, Millikan RC, Thomas NE, Gruber SB, Anton-Culver H, Rosso S, Gallagher RP, Dwyer T, Kanetsky PA, Busam K, From L, Begg CB, Berwick M, GEM Study Group: Vitamin D receptor polymorphisms in patients with cutaneous melanoma. Int J Cancer. 2012, 130: 405-418. 10.1002/ijc.26023.View ArticlePubMedGoogle Scholar
- Chang YM, Newton-Bishop JA, Bishop DT, Armstrong BK, Bataille V, Bergman W, Berwick M, Bracci PM, Elwood JM, Ernstoff MS, Green AC, Gruis NA, Holly EA, Ingvar C, Kanetsky PA, Karagas MR, Le Marchand L, Mackie RM, Olsson H, Østerlind A, Rebbeck TR, Reich K, Sasieni P, Siskind V, Swerdlow AJ, Titus-Ernstoff L, Zens MS, Ziegler A, Barrett JH: A pooled analysis of melanocytic nevus phenotype and the risk of cutaneous melanoma at different latitudes. Int J Cancer. 2009, 124: 420-428. 10.1002/ijc.23869.View ArticlePubMedPubMed CentralGoogle Scholar
- Gandini S, Sera F, Cattaruzza MS, Pasquini P, Abeni D, Boyle P, Melchi CF: Meta-analysis of risk factors for cutaneous melanoma: I. Common and atypical naevi. Eur J Cancer. 2005, 41: 28-44. 10.1016/j.ejca.2004.10.015.View ArticlePubMedGoogle Scholar
- Bishop DT, Demenais F, Iles MM, Harland M, Taylor JC, Corda E, Randerson-Moor J, Aitken JF, Avril MF, Azizi E, Bakker B, Bianchi-Scarrà G, Bressac-de Paillerets B, Calista D, Cannon-Albright LA, Chin-A-Woeng T, Debniak T, Galore-Haskel G, Ghiorzo P, Gut I, Hansson J, Hocevar M, Höiom V, Hopper JL, Ingvar C, Kanetsky PA, Kefford RF, Landi MT, Lang J, Lubiński J: Genome-wide association study identifies three loci associated with melanoma risk. Nature Genet. 2009, 41: 920-925. 10.1038/ng.411.View ArticlePubMedPubMed CentralGoogle Scholar
- Falchi M, Bataille V, Hayward NK, Duffy DL, Bishop JA, Pastinen T, Cervino A, Zhao ZZ, Deloukas P, Soranzo N, Elder DE, Barrett JH, Martin NG, Bishop DT, Montgomery GW, Spector TD: Genome-wide association study identifies variants at 9p21 and 22q13 associated with development of cutaneous nevi. Nature Genet. 2009, 41: 915-919. 10.1038/ng.410.View ArticlePubMedGoogle Scholar
- 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, Lazaro P, Benitez J, Ribas G: Genetic analysis of the vitamin D receptor Gene in two epithelial cancers: melanoma and breast cancer case–control studies. BMC Cancer. 2008, 8: 385-10.1186/1471-2407-8-385.View ArticlePubMedPubMed CentralGoogle Scholar
- Aguilera P, Puig S, Guilabert A, Julia M, Romero D, Vicente A, Gonzalez-Ensenat MA, Malvehy J: Prevalence study of nevi in children from Barcelona. Dermoscopy, constitutional and environmental factors. Dermatology. 2009, 218: 203-214. 10.1159/000183179.View ArticlePubMedGoogle Scholar
- Ersfeld DL, Rao DS, Body JJ, Sackrison JL, Miller AB, Parikh N, Eskridge TL, Polinske A, Olson GT, MacFarlane GD: Analytical and clinical validation of the 25 OH vitamin D assay for the LIAISON automated analyzer. Clin Biochem. 2004, 37: 867-874. 10.1016/j.clinbiochem.2004.06.006.View ArticlePubMedGoogle Scholar
- de Bakker PI, Yelensky R, Pe’er I, Gabriel SB, Daly MJ, Altshuler D: Efficiency and power in genetic association studies. Nature Genet. 2005, 37: 1217-1223. 10.1038/ng1669.View ArticlePubMedGoogle Scholar
- Gonzalez JR, Armengol L, Sole X, Guino E, Mercader JM, Estivill X, Moreno V: SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007, 23: 644-645.PubMedGoogle Scholar
- Lips P: Vitamin D physiology. Prog Biophys Mol Biol. 2006, 92: 4-8. 10.1016/j.pbiomolbio.2006.02.016.View ArticlePubMedGoogle Scholar
- Guerrieri-Gonzaga A, Gandini S: Vitamin D and overall mortality. Pigment Cell Melanoma Res. 2013, 26: 16-28. 10.1111/pcmr.12036.View ArticlePubMedGoogle Scholar
- Norman AW, Bouillon R, Whiting SJ, Vieth R, Lips P: 13th Workshop consensus for vitamin D nutritional guidelines. J Steroid Biochem Mol Biol. 2007, 103: 204-5. 10.1016/j.jsbmb.2006.12.071.View ArticlePubMedPubMed CentralGoogle Scholar
- Ross AC, Manson JE, Abrams SA, Aloia JF, Brannon PM, Clinton SK, Durazo-Arvizu RA, Gallagher JC, Gallo RL, Jones G, Kovacs CS, Mayne ST, Rosen CJ, Shapses SA: The 2011 report on dietary reference intakes for calcium and vitamin D from the institute of medicine: what clinicians need to know. J Clin Endocrinol Metab. 2011, 96: 53-8. 10.1210/jc.2010-2704.View ArticlePubMedGoogle Scholar
- Zittermann A, Iodice S, Pilz S, Grant WB, Bagnardi V, Gandini S: Vitamin D deficiency and mortality risk in the general population: a meta-analysis of prospective cohort studies. Am J Clin Nutr. 2012, 95: 91-100. 10.3945/ajcn.111.014779.View ArticlePubMedGoogle Scholar
- Ovesen L, Andersen R, Jakobsen J: Geographical differences in vitamin D status, with particular reference to European countries. Proc Nutr Soc. 2003, 62: 813-821. 10.1079/PNS2003297.View ArticlePubMedGoogle Scholar
- Vaqueiro M, Bare ML, Anton E, Andreu E, Gimeno C, Group DA: Evaluation assessment of the cut-off point of vitamin D in the population older than 64 years old. Med Clin (Barc). 2006, 127: 648-650. 10.1157/13094820.View ArticleGoogle Scholar
- Mata-Granados JM, Luque de Castro MD, Quesada Gomez JM: Inappropriate serum levels of retinol, alpha-tocopherol, 25 hydroxyvitamin D3 and 24,25 dihydroxyvitamin D3 levels in healthy Spanish adults: simultaneous assessment by HPLC. Clin Biochem. 2008, 41: 676-680. 10.1016/j.clinbiochem.2008.02.003.View ArticlePubMedGoogle Scholar
- Calatayud M, Jodar E, Sanchez R, Guadalix S, Hawkins F: Prevalence of deficient and insufficient vitamin D levels in a young healthy population. Endocrinol Nutr. 2009, 56: 164-169. 10.1016/S1575-0922(09)70980-5.View ArticlePubMedGoogle Scholar
- González-Molero I, Morcillo S, Valdés S, Pérez-Valero V, Botas P, Delgado E, Hernández D, Olveira G, Rojo G, Gutierrez-Repiso C, Rubio-Martín E, Menéndez E, Soriguer F: Vitamin D deficiency in Spain: a population-based cohort study. Eur J Clin Nutr. 2011, 65: 321-328. 10.1038/ejcn.2010.265.View ArticlePubMedGoogle Scholar
- Newton-Bishop JA, Chang YM, Elliott F, Chan M, Leake S, Karpavicius B, Haynes S, Fitzgibbon E, Kukalizch K, Randerson-Moor J, Elder DE, Bishop DT, Barrett JH: Relationship between sun exposure and melanoma risk for tumours in different body sites in a large case–control study in a temperate climate. Eur J Cancer. 2011, 47: 732-741. 10.1016/j.ejca.2010.10.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Davies JR, Chang YM, Snowden H, Chan M, Leake S, Karpavicius B, Haynes S, Kukalizch K, Randerson-Moor J, Elliott F, Barth J, Kanetsky PA, Harland M, Bishop DT, Barrett JH, Newton-Bishop JA: The determinants of serum vitamin D levels in participants in a melanoma case–control study living in a temperate climate. Cancer Causes Control. 2011, 22: 1471-1482. 10.1007/s10552-011-9827-3.View ArticlePubMedPubMed CentralGoogle Scholar
- Malvy DJ, Guinot C, Preziosi P, Galan P, Chapuy MC, Maamer M, Arnaud S, Meunier PJ, Hercberg S, Tschachler E: Relationship between vitamin D status and skin phototype in general adult population. Photochem Photobiol. 2000, 71: 466-469. 10.1562/0031-8655(2000)071<0466:RBVDSA>2.0.CO;2.View ArticlePubMedGoogle Scholar
- Glass D, Lens M, Swaminathan R, Spector TD, Bataille V: Pigmentation and vitamin D metabolism in Caucasians: low vitamin D serum levels in fair skin types in the UK. PLoS One. 2009, 4: e6477-10.1371/journal.pone.0006477.View ArticlePubMedPubMed CentralGoogle Scholar
- Springbett P, Buglass S, Young AR: Photoprotection and vitamin D status. J Photochem Photobiol B. 2010, 101: 160-168. 10.1016/j.jphotobiol.2010.03.006.View ArticlePubMedGoogle Scholar
- Berwick M: The good, the bad, and the ugly of sunscreens. Clin Pharmacol Ther. 2011, 89: 31-33. 10.1038/clpt.2010.227.View ArticlePubMedGoogle Scholar
- Kiel DP, Demissie S, Dupuis J, Lunetta KL, Murabito JM, Karasik D: Genome-wide association with bone mass and geometry in the Framingham heart study. BMC Med Genet. 2007, 8 (Suppl 1): S14-10.1186/1471-2350-8-S1-S14.View ArticlePubMedPubMed CentralGoogle Scholar
- Karami S, Brennan P, Rosenberg PS, Navratilova M, Mates D, Zaridze D, Janout V, Kollarova H, Bencko V, Matveev V, Szeszenia-Dabrowska N, Holcatova I, Yeager M, Chanock S, Menashe I, Rothman N, Chow WH, Boffetta P, Moore LE: Analysis of SNPs and haplotypes in vitamin D pathway genes and renal cancer risk. PLoS One. 2009, 4: e7013-10.1371/journal.pone.0007013.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonilla C, Hooker S, Mason T, Bock CH, Kittles RA: Prostate cancer susceptibility loci identified on chromosome 12 in African Americans. PLoS One. 2011, 6: e16044-10.1371/journal.pone.0016044.View ArticlePubMedPubMed CentralGoogle Scholar
- Carling T, Rastad J, Akerstrom G, Westin G: Vitamin D receptor (VDR) and parathyroid hormone messenger ribonucleic acid levels correspond to polymorphic VDR alleles in human parathyroid tumors. J Clin Endocrinol Metab. 1998, 83: 2255-2259. 10.1210/jc.83.7.2255.PubMedGoogle Scholar
- Millikan RC, Hummer A, Begg C, Player J, de Cotret AR, Winkel S, Mohrenweiser H, Thomas N, Armstrong B, Kricker A, Marrett LD, Gruber SB, Culver HA, Zanetti R, Gallagher RP, Dwyer T, Rebbeck TR, Busam K, From L, Mujumdar U, Berwick M: Polymorphisms in nucleotide excision repair genes and risk of multiple primary melanoma: the genes environment and melanoma study. Carcinogenesis. 2006, 27: 610-618.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/14/26/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.