MC1Rgenotype as a predictor of early-onset melanoma, compared with self-reported and physician-measured traditional risk factors: an Australian case-control-family study
BMC Cancer volume 13, Article number: 406 (2013)
Melanocortin-1 receptor (MC1R) gene variants are very common and are associated with melanoma risk, but their contribution to melanoma risk prediction compared with traditional risk factors is unknown. We aimed to 1) evaluate the separate and incremental contribution of MC1R genotype to prediction of early-onset melanoma, and compare this with the contributions of physician-measured and self-reported traditional risk factors, and 2) develop risk prediction models that include MC1R, and externally validate these models using an independent dataset from a genetically similar melanoma population.
Using data from an Australian population-based, case-control-family study, we included 413 case and 263 control participants with sequenced MC1R genotype, clinical skin examination and detailed questionnaire. We used unconditional logistic regression to estimate predicted probabilities of melanoma. Results were externally validated using data from a similar study in England.
When added to a base multivariate model containing only demographic factors, MC1R genotype improved the area under the receiver operating characteristic curve (AUC) by 6% (from 0.67 to 0.73; P < 0.001) and improved the quartile classification by a net 26% of participants. In a more extensive multivariate model, the factors that contributed significantly to the AUC were MC1R genotype, number of nevi and previous non-melanoma skin cancer; the AUC was 0.78 (95% CI 0.75-0.82) for the model with self-reported nevi and 0.83 (95% CI 0.80-0.86) for the model with physician-counted nevi. Factors that did not further contribute were sun and sunbed exposure and pigmentation characteristics. Adding MC1R to a model containing pigmentation characteristics and other self-reported risk factors increased the AUC by 2.1% (P = 0.01) and improved the quartile classification by a net 10% (95% CI 1-18%, P = 0.03).
Although MC1R genotype is strongly associated with skin and hair phenotype, it was a better predictor of early-onset melanoma than was pigmentation characteristics. Physician-measured nevi and previous non-melanoma skin cancer were also strong predictors. There might be modest benefit to measuring MC1R genotype for risk prediction even if information about traditional self-reported or clinically measured pigmentation characteristics and nevi is already available.
Melanoma is one of the most common cancers and a leading cause of cancer death in young adults of European origin [1, 2]. People identified as being at high risk of melanoma could likely benefit from regular skin checks and interventions to improve sun-protection behaviours [3, 4]. Phenotypic characteristics (e.g. hair, eye and skin colour, skin sensitivity to sunlight, number of nevi (moles)), family history, past sun exposure and past history of skin cancer are usually the basis for discriminating individual risk of melanoma [5–7]. However, given the decreasing costs and increasing use of genetic testing, it is becoming more feasible to incorporate genetic risk factors into clinical risk prediction tools. Low penetrant genetic variants for the melanocortin-1 receptor (MC1R) gene [8, 9] are very common in populations of European origin [10, 11] and some of these variants have been associated with a 1.5 to 4-fold increased risk of melanoma [12–15]. MC1R variants are associated with sun-sensitive phenotypes but the association with melanoma appears to be mediated also through non-pigmentary pathways [12, 15]. To date, only three, preliminary, melanoma risk prediction models have included MC1R genotype [16–18]. No study has formally assessed the contribution of MC1R genotype to melanoma risk prediction compared with traditional factors.
The Australian Melanoma Family Study is a multi-centre, population-based, case-control-family study of early-onset melanoma (diagnosis before 40 years of age) that has comprehensive data for MC1R genotype and traditional risk factors including phenotype, UV (ultraviolet) radiation and clinically measured nevus counts . Using data from this study, we aimed to 1) evaluate the separate and incremental contribution of MC1R genotype to prediction of early-onset melanoma, and compare this with the contributions of physician-measured and self-reported traditional risk factors, and 2) develop risk prediction models that include MC1R, and externally validate these models using an independent dataset from a genetically similar melanoma population .
The Australian Melanoma Family Study design, recruitment, data collection and participant characteristics have been previously detailed . Cases and controls were living in Brisbane, Sydney or Melbourne, which comprise about 50% of Australia’s population. Approval for the study was obtained from the ethics committees of The University of Sydney, The University of Melbourne, The University of Queensland, Cancer Council Victoria, Queensland Cancer Register and Cancer Council NSW. All participants provided written, informed consent.
Cases were identified from population-based state cancer registries, diagnosed between 1st July 2000 and 31st December 2002 at ages 18–39 years with incident, histopathologically-confirmed, first-primary invasive cutaneous melanoma. A total of 629 cases were recruited; participation was 54% of those eligible and 76% of those contactable.
Population controls were aged between 18 and 39 years at the time of approach and had no history of invasive or in situ melanoma. They were selected from the electoral roll (registration to vote is compulsory for Australian citizens aged 18 years and over) and were frequency-matched to cases by city, age and sex. A total of 240 population controls were recruited; participation was 23% of those apparently eligible and 42% of those contactable.
Eligible spouse or friend controls were a spouse, partner, or friend nominated by a case as a potential control participant. They were eligible if they were at least 18 years of age and had no history of invasive or in situ melanoma; there were no other age, sex or residency restrictions. A potential control was nominated by 59% of cases. A total of 295 spouse or friend controls were recruited; participation was 80% of those nominated. Population-controls and spouse or friend-controls were combined into one control group as done previously [15, 19].
Data were collected by telephone interview using a structured questionnaire, which included detailed questions on sun exposure, phenotype, residence history, demographic information, ancestry and diagnoses of cancer and non-melanoma skin cancers (basal cell carcinoma and squamous cell carcinoma) [19, 21]. Participants also reported their skin colour and type, eye colour, natural hair colour at age 18 years, usual tanning and sunburn response to prolonged or repeated exposure to sunlight in summer, sunbed use, the number of nevi covering the body (described pictorially as none, few, some, many), freckling in childhood and adulthood, and nevus count on the back. Reported melanoma in relatives was verified where possible .
Clinical skin examinations
All case and control participants were invited to attend clinical skin examinations, which were conducted at dermatology clinics in Brisbane, Sydney, and Melbourne by dermatology trainees trained on the study protocol. A clinical skin examination was completed by 73% of cases, 55% of population controls and 67% of spouse or friend controls.
Measurement of nevi was based on international guidelines . Separate counts were made for melanocytic nevi of 2-5 mm and >5 mm, raised nevi of >2 mm, and clinically atypical nevi of >2 mm, on 30 body sites. The number of solar lentigines on the upper back was recorded by using a 6-level picture scale. Natural hair colour at age 18 and eye colour were recorded using wig hair swatches and eye photographs. Reflected skin colour, a correlate of melanin content , was recorded using a hand-held reflectance spectrophotometer with standard reflectance at 685 nm. The multi-wavelength data quantify colour using the Commission Internationale de l'Éclairage L* a* b* colour space parameters . Inner arm L* values describe base skin colour, b* values describe tanning, and a* values describe erythema [23–25].
MC1Rgenotyping and classification
Blood samples were requested from all participants and were obtained from 597 (95%) cases, 220 (92%) population controls, and 256 (87%) spouse/friend controls. The methods for MC1R genotyping and classification have been described in detail elsewhere . Briefly, we sequenced MC1R and classified variants D84E, R142H, R151C, I155T, R160W, D294H as ‘R’ variants and all other variants excluding synonymous changes and non-coding changes as ‘r’ variants. R variants have been shown to be strongly associated with the presence of ‘red hair colour phenotype’ (red hair, fair skin, freckling, poor sun sensitivity), whereas r variants generally have a weaker association with red-hair colour phenotype . The association of the individual MC1R variants with melanoma risk in this sample has been described previously .
When MC1R genotype was added to the statistical models, it was added together as a group of seven separate variables: one for each of the six ‘R’ variants D84E, R142H, R151C, I155T, R160W, D294H, and one variable for all ‘r’ variants combined. Each of these variables was formatted to indicate the number of variant alleles (i.e. 0,1,2).
In order to compare the contribution of MC1R genotype with all self-reported and clinical traditional risk factors simultaneously, we restricted this analysis to case and control participants who had: complete questionnaire data for the main risk factors examined, a clinical skin examination, MC1R genotype, self-reported exclusive European ancestry, and were aged < 45 years at interview. After exclusions, 676 participants remained for the analysis: 413 cases and 263 controls (115 population-controls and 148 spouse or friend controls). Data were analysed using SAS version 9.2 (SAS Institute, Cary NC) and statistical significance was inferred at two-sided P < 0.05.
In the ‘base’ model, melanoma status was the outcome variable and covariates included demographic factors: age (quadratic), sex, city of recruitment (Brisbane, Sydney, Melbourne), and self-reported European ancestry (British/northern, southern, eastern/mixed/other European) to account for any difference in MC1R allele frequencies across ethnic groups .
We added MC1R genotype and traditional risk factors separately to the base model to evaluate their individual contribution to risk prediction. We also added the risk factors incrementally to the base model in order of their contribution to the area under the receiver operating characteristic (ROC) curve (AUC). To examine the contribution of traditional pigmentation characteristics, we created a pigmentation-related propensity-to-melanoma score (‘pigmentation score’) continuous variable that summarizes the contribution of six correlated, categorical phenotypic variables, including self-reported ability to tan, propensity to sunburn, childhood freckling, skin colour, eye colour, and hair colour . For the more objectively-measured pigmentation score, the last three self-reported variables were replaced with physician-measured skin reflectance, eye colour, and hair colour.
Other self-reported variables that were tested in the models were number of nevi fitted as a categorical variable (none, few, some, many), previous diagnosis of non-melanoma skin cancer (yes, no), and ultraviolet (UV) radiation related exposures: total childhood sun exposure hours (quartiles), childhood blistering sunburns (none, ≤ 8, > 8) and lifetime sunbed use (none, 1–10, >10 sessions). The two childhood measures were chosen over other sun exposure measures such as lifetime, adulthood, weekday and holiday sun exposure, because they were more predictive of melanoma in our study sample. Other physician-measured variables that were tested in the models were separate nevus counts (≥ 2 mm, 2–5 mm, ≥ 5 mm, dysplastic, raised) and solar lentigines. We also included confirmed family history of melanoma in a first-degree relative.
Measures of model performance
As measures of discrimination, i.e. the ability of a model to discriminate those who will develop melanoma from those who will not, we calculated: the AUC, which is equivalent to the concordance (c) statistic; the net reclassification improvement (NRI); discrimination slope; and the integrated discrimination index (IDI) [26–29]. To assess calibration, i.e. the agreement between observed and predicted outcomes, we used the Hosmer-Lemeshow goodness-of-fit test [26, 30]. These measures were based on predicted probabilities of melanoma from the unconditional logistic regression models described above.
The AUC is equal to the probability that, for one case and one control chosen at random from the data set, the predicted probability of melanoma is higher for the case than for the control, and ranges from 0.5 (equivalent to a coin toss) to 1.0 (perfect discrimination). The NRI quantifies overall improvement in model sensitivity and specificity. A net improvement in risk classification implies upward reclassification of case participants and downward reclassification of control participants. The NRI was calculated by first fitting a ‘base model’ which grouped participants into quartiles of their predicted probability of melanoma; these quartile distributions were then compared to the ‘comparison model’. Improvement in sensitivity represents net reclassification of more cases into higher quartiles, improvement in specificity represents net reclassification of controls into lower quartiles, and overall improvement in classification combines the improvements in sensitivity and specificity. In the absence of clinically meaningful cut-points, we used quartiles to define risk categories as done elsewhere . We also calculated the ‘category-free’ NRI, for which the definition of upward or downward movement is simplified to indicate any increase or decrease in probabilities of the outcome . Discrimination slope was calculated as the difference between the mean predicted probability for cases and controls, and the IDI was calculated as the difference between discrimination slopes between the base and comparison models; both of these measures do not require predefined risk categories. As measures of overall model performance, we estimated the Brier score and Nagelkerke’s R2, which are measures of how well future outcomes are likely to be predicted by the model ; a higher Nagelkerke’s R2 and a lower Brier score indicates better predictability of the model.
As a measure of internal validation, we used 100 bootstrap samples to estimate the AUC and Nagelkerke’s R2 for the final models. Odds ratios (OR) for melanoma and their 95% confidence intervals were estimated using unconditional logistic regression models.
We performed external validation of the final regression models using a population-based case–control study of melanoma from a geographically defined area of Yorkshire and the Northern region of the United Kingdom . Case participants had incident pathologically confirmed invasive melanoma diagnosed between September 2000 and December 2005 (67% case participation). Control participants were identified from the cases’ family doctors (55% response) and were frequency-matched to cases by age and sex. A total of 841 case participants and 452 control participants, aged between 18 and 76 years, were included in this analysis. This study was conducted in tandem with the Australian case–control study and used a common protocol for collection of phenotype and sun exposure measures to facilitate comparisons among the datasets. We handled the data variables and analysis in the same way for both datasets.
Characteristics of the study sample
Demographic characteristics and selected risk factors of early-onset cases and controls are shown in Table 1. Fifty-eight percent of cases and 40% of controls had at least one R allele. A previous non-melanoma skin cancer was reported by 8% of cases and 2% of controls.
Separate contribution of MC1Rgenotype and traditional factors
Compared to the base model, the separate addition of MC1R, pigmentation score, nevi, non-melanoma skin cancer and solar lentigines each considerably improved the discriminative ability of the model, whereas inclusion of self-reported sun and sunbed exposure variables (childhood sun exposure, childhood blistering sunburns and lifetime sunbed use) resulted in minimal improvement, and inclusion of family history resulted in no improvement (Table 2).
When added to the base model, MC1R improved the AUC by 6%, sensitivity by 12% (95% CI 5-19%), specificity by 14% (95% CI 6-23%), and improved the quartile classification for a net 26% (95% CI 15-37%) of participants. Further examination showed that the six ‘R’ variants were responsible for most of the improvement to risk prediction, as together they increased the AUC by 4% (P = 0.001) and improved the quartile classification by a net 21% (95% CI 10-31%) of participants whereas the combined ‘r’ variants increased the AUC by less than 1% (P = 0.5) and the net reclassification improvement by 5% (−3-12%) (Table 2).
The contribution of traditional pigmentation characteristics to model improvement was similar for self-reported and the more objectively-measured pigmentation score. Physician-counted number of nevi ≥ 2 mm and 2–5 mm were the nevi variables most predictive of melanoma risk, whereas self-reported number of nevi had a more modest impact. There was no material change to any of our results in this paper when we repeated the models, replacing the single composite ‘pigmentation score’ variable with the six separate variables that comprise the pigmentation score (data not shown). We also tested hair colour as a separate variable and found that it contributed about half as much to the AUC compared to the pigmentation score variable.
Incremental contribution of MC1Rgenotype and traditional factors
In a more extensive multivariate model where each risk factor was added incrementally to the base model in order of their contribution to increasing the AUC, only MC1R, number of nevi and history of non-melanoma skin cancer significantly improved the AUC for both the self-reported and physician-measured models (Table 3). Self-reported pigmentation score weakly increased (by 1%; P = 0.07) the AUC for the self-reported model already containing MC1R, nevi and non-melanoma skin cancer, whereas more objectively-measured pigmentation score did not increase the discrimination of the corresponding physician-measured model. MC1R and number of nevi were the only variables that produced significant quartile reclassification of cases and controls. Measures of sun and sunbed exposure and solar lentigines did not increase the discrimination of the models already containing the other factors; nor did number of dysplastic nevi or raised nevi, once number of nevi ≥ 2 mm (the most predictive nevus variable) was included in the physician-measured model.
Selection and validation of final models, and measures of model performance
Based on improvement to the AUC, the final selected models for both the self-reported and physician-measured models included MC1R, nevi and non-melanoma skin cancer, in addition to demographic factors. Details of the models’ performance and validation are shown in Table 4. The AUC was higher for the physician-measured model (0.83, 95% CI 0.80-0.86) than for the self-reported model (0.78, 95% CI 0.75-0.82), a difference in the AUC of 0.043 (P < 0.001), reflecting better predictive ability of clinically-measured number of nevi than self-reported nevi.
Compared to the base model, the self-reported model improved classification for a net 37% (95% CI 25-48%) of participants based on quartile cut-points and 63% (95% CI 47-78%) using the category-free approach; for the physician-measured model, net reclassification improvement was 53% (95% CI 41-64%) and 85% (95% CI 70-101%), respectively. Overall model performance also improved: Nagelkerke’s R2 increased from 13% in the base model to 32% for the self-reported model and 39% for the physician-measured model, and the Brier score decreased. Internal validation produced similar results for Nagelkerke’s R2 and the AUC. The discrimination slopes for each model (presented as box plots in Additional file 1: Figure S1), show how the physician-measured model achieved the best separation of predicted probabilities between cases and controls.
External validation of the final regression models using data from the English study showed slightly lower discrimination for the self-reported and physician-measured models compared to our Australian study results. However, this appeared to be due to lower discrimination for the baseline model (AUC 0.61 compared to 0.67), as both studies demonstrated similar improvements to the AUC, NRI and Nagelkerke’s R2 for the self-reported and physician-measured models when compared to the respective base model (Table 4). For both studies, the Hosmer-Lemeshow test indicated poor calibration (P <0.05) for the physician-measured model.
Using ROC curves (Figure 1), we estimated the proportion of cases and controls that would be classified as high-risk using different cut-points. Choosing a cut-point value of 90% specificity, equating to 10% of controls being classified as high-risk, the proportion of cases identified as high-risk (i.e. sensitivity) was 23% in the base model, 40% in the self-reported model, and 53% in the physician-measured model. When we chose a different cut-point based on a balance of sensitivity and specificity (by selecting the value in the top, left-hand corner of the ROC plot), sensitivity and specificity were each 62% in the base model, 70% in the self-reported model and 74% in the physician-measured model.
OR estimates for MC1R, nevi and non-melanoma skin cancer are presented in Table 5. Previously, we observed a stronger association between MC1R and melanoma for men than women in this study (P = 0.005) ; in this analysis the multivariate OR for any R allele compared to wild-type consensus alleles was 4.20 (95% CI 1.75-10.10) for men and 2.44 (95% CI 1.31-4.56) for women.
Incremental contribution of MC1Rgenotype to a model including traditional pigmentation variables, nevi and non-melanoma skin cancer
We assessed whether MC1R further improved prediction of melanoma when traditional pigmentation characteristics were already in the model, by replacing MC1R with ‘pigmentation score’ in the final regression models and then testing the addition of MC1R. We found that MC1R increased the AUC by 2.1% (from 0.77 to 0.79, P = 0.01) for the self-reported model and 1.3% (from 0.82 to 0.83, P = 0.04) for the physician-measured model; and improved the quartile classification for a net 10% (95% CI 1-18%, P = 0.03) of participants for the self-reported model and 6% (95% CI −1-13%, P = 0.08) for the physician-measured model (data not shown in tables).
MC1R, nevi and personal history of non-melanoma skin cancer were identified as the strongest predictors of melanoma risk in our study of early-onset melanoma. The contribution of MC1R to prediction of melanoma was similar to that obtained from measuring self-reported nevi, which is considered a strong and discriminative risk factor [5, 33, 34]. When added separately to the base model, MC1R and self-reported nevi increased the AUC by 6% and 5% respectively, and both improved classification for about a quarter of the cases and controls through net movement of cases into higher quartiles and controls into lower quartiles of predicted risk. Total number of nevi measured by physicians (dermatology trainees) was the strongest predictor of risk overall, increasing the AUC by 11% and reclassifying 39% of participants. Although MC1R genotype is strongly associated with skin and hair phenotype [12, 14], it was a better predictor of early-onset melanoma than was pigmentation characteristics.
Our models demonstrated high discrimination: an AUC of 0.78 for the self-reported model and 0.83 for the physician-measured model in the Australian (development) dataset and 0.71 and 0.79 in the English (validation) dataset. The additional predictive value of MC1R, nevi and non-melanoma skin cancer variables when added to the demographic ‘base’ model was similar for both studies, suggesting good generalisability of our results to other genetically similar populations. The differences in the AUC for the base models of the two studies is less important, as it is strongly influenced by how age, sex, ethnic and regional differences have already been accounted for in the study design. It is expected that models will perform better on the development dataset than the validation dataset because of overfitting.
Three other, preliminary, melanoma risk prediction models containing MC1R genotype have been published. Whiteman and Green  published a prototype for a melanoma risk prediction tool but provided no details on predictive performance. In a published conference abstract, Smith et al. showed an AUC of 0.72 (95% CI 0.70-0.75) for a model containing conventional risk factors and 0.75 (95% CI 0.72-0.77) when they added MC1R genotype, outdoor UV and indoor UV exposure, based on data from a case–control study of people aged 25–59 years in Minnesota. In a Greek study with 284 cases and 284 controls, Stefanaki and colleagues  derived a melanoma risk prediction model containing phenotypic traits (except nevi was not available) and 8 single nucleotide polymorphisms (SNPs) from several genes that included the MC1R locus. They found no appreciable change to the AUC after the addition of the 8 SNPs (AUC changed from 0.833 to 0.839), which they suggested may have been partly due to lower risk allele frequencies in their Greek population compared to other European populations . Measurement of single SNPs rather than the causal variants might also underestimate the contribution of genetic variation to melanoma risk.
Other published risk prediction models for melanoma have reported AUCs in the range of 0.54 to 0.86 [34–41]. The highest reported AUC of 0.86 was for a model containing age, hair colour, personal history of melanoma and suspicious melanocytic lesion on dermoscopy, developed using a German cohort . We were not able to include personal history of melanoma in our models because the study eligibility criteria specifically excluded these cases; however previous primary melanoma would be rare in people younger than 40 years and thus would be less relevant to our analysis than in studies with older participants. Previous non-melanoma skin cancer was a very strong risk factor in our study (OR 9), however because it had a low prevalence, the contribution to the AUC and more particularly to NRI was modest. This can be a limitation of these prediction methods, as factors that increase risk substantially can have minimal overall improvement to the model if they are rare. Another risk factor affected in this way is family history of melanoma. This variable did not improve discrimination in our analysis. We used a relatively low threshold for defining positive family history, requiring only one confirmed melanoma in a first-degree relative. This definition is consistent with many other population-based studies  and is associated with approximately two-fold increased melanoma risk [7, 19]. However, a more extensive family history is associated with higher risk estimates  and thus, for some individuals, this factor will strongly influence their personal melanoma risk. This issue demonstrates the different priorities that are placed on the design of risk prediction tools for different settings: one type for stratifying the population into broad risk categories to aid primary prevention strategies, and the other type for more precisely estimating personal risk of melanoma.
Harbauer et al. developed a model containing number of nevi, skin phototype and skin UV damage, adjusted for age and sex, with an AUC of 0.73 (95% CI 0.68-0.77) when measured by self-report and 0.77 (95% CI 0.73-0.83) when measured by a dermatologist. Similarly, our physician-measured model also achieved better discrimination than our self-report model, which was due solely to improved measurement of nevi. Our results indicate that counting all nevi greater than 2 mm is a better predictor of melanoma than counting only large, dysplastic or raised nevi. There was no added benefit to having physicians measure hair colour, eye colour or skin colour, as the discriminative ability of the pigmentation score was similar whether measured by self-report or by a physician.
Including sun and sunbed exposure variables resulted in little improvement to discrimination; when added as a group to the base model, they increased the AUC by 1.5%, but this diminished when the other factors were in the model. Although these sun and sunbed exposures have been shown to be associated with melanoma in our study [21, 43], it has been demonstrated that very strong, independent associations with risk are required to meaningfully increase the AUC [29, 44, 45]. Sun exposure-related factors generally do not have very strong effect estimates, which may be partly due to inherent difficulties measuring past sun exposure  but also reflects the strong influence of genetically determined risk factors.
The Hosmer-Lemeshow test of calibration indicated poor agreement between observed and predicted outcomes for the physician-measured model, for both the development and external validation datasets. We examined the observed and expected values within each decile from the Hosmer-Lemeshow test but there no was consistent pattern describing the differences. Calibration has been rarely reported for other published melanoma risk prediction models. A prospective cohort study with follow-up of individuals for development of melanoma would be the ideal method to evaluate melanoma prediction probabilities.
Several strengths and limitations of the Australian Melanoma Family Study have been discussed previously [15, 19, 21]. We had low participation from cases and population controls, which is a common problem for population-based studies , especially when conducted with highly mobile, young adults . Although poor participation can sometimes lead to selection bias , we did not find strong evidence of this occurring in our study . This analysis was restricted to participants who had a complete set of data for self-reported and clinically-measured risk factors. Having a clinical skin examination was a preferred but not compulsory aspect of participation in the Australian Melanoma Family Study, and thus was only completed by a subset of case and control participants. Although this reduced the sample size for our analysis, it is unlikely to have introduced systematic bias, as a comparison of those with (n = 676) versus without (n = 270) clinical skin examinations showed no statistically significant differences on predictors including sex, self-reported nevi, hair colour, MC1R genotype, previous non-melanoma skin cancer, and childhood blistering sunburns.
Our study focused on prediction of early-onset melanoma. A younger age at diagnosis is more likely to reflect an underlying genetic susceptibility to the disease, thus the contribution of MC1R and traditional factors to risk prediction in our study may differ for those diagnosed at older ages or for other ethnicities. Nevertheless, our external validation results do suggest generalisability of our results to a broader-aged population. The generalisability of our results may also differ by country and the phenotypic and behavioural characteristics of the population. External validation of risk prediction models using an independent dataset has been rarely performed by other studies, and is an important strength of our study.
Physician-measured nevi, MC1R genotype and previous non-melanoma skin cancer were the strongest predictors of early-onset melanoma in this study. Our results suggest that there might be modest benefit to measuring MC1R genotype for risk prediction even if information about traditional self-reported or clinically measured pigmentation characteristics and nevi is already available. Although many nonsynonymous MC1R gene variants exist [10, 48], the six ‘R’ variants were responsible for most of the improvement to risk prediction in this study and thus would be the most important MC1R gene variants to include in melanoma prediction tools if it was not feasible or economically justifiable to measure all variants. We had limited statistical power to evaluate the predictive effect of individual R variants, however both our previous study  and a meta-analysis  attributed the R151C variant with the highest population burden based on prevalence and relative risk of melanoma.
Our study results will help guide the development of melanoma risk prediction tools that incorporate MC1R genotype. Decreasing genotyping costs and increasing use of genetic testing is making it more feasible to incorporate genetic risk factors into clinical risk prediction tools; however, translation into routine clinical practice requires several additional steps [49, 50]. As a screening tool, evidence is needed to show whether or not individuals identified at high risk of melanoma will improve their sun-protection behaviours and perform regular skin checks, that improved outcomes justify the associated costs, and that the benefits of obtaining this information outweigh any disadvantages for patients and their families. For clinical genetic counselling, information on common variants in other melanoma susceptibility genes may need to be incorporated in the risk prediction tool in order to more precisely estimate personal risk of melanoma. Ultimately, the clinical and public health value of such tools will depend not only on its predictive performance, but also on the application setting, feasibility, cost effectiveness, benefits and harms [51, 52].
Receiver operating characteristic curve
Net reclassification improvement
Area under the ROC curve
Integrated discrimination index
Single nucleotide polymorphisms.
Australian Institute of Health and Welfare (AIHW): Cancer in Australia: an overview, 2008. Cancer Series no 46, AIHW cat no CAN 42. 2008, Canberra: AIHW
CancerStats - cancer statistics for the UK. http://info.cancerresearchuk.org/cancerstats,
Aitken JF, Elwood M, Baade PD, Youl P, English D: Clinical whole-body skin examination reduces the incidence of thick melanomas. Int J Cancer. 2010, 126 (2): 450-458. 10.1002/ijc.24747.
Breitbart EW, Waldmann A, Nolte S, Capellaro M, Greinert R, Volkmer B, Katalinic A: Systematic skin cancer screening in Northern Germany. J Am Acad Dermatol. 2012, 66 (2): 201-211. 10.1016/j.jaad.2010.11.016.
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 (1): 28-44. 10.1016/j.ejca.2004.10.015.
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 (1): 45-60. 10.1016/j.ejca.2004.10.016.
Gandini S, Sera F, Cattaruzza MS, Pasquini P, Zanetti R, Masini C, Boyle P, Melchi CF: Meta-analysis of risk factors for cutaneous melanoma: III. Family history, actinic damage and phenotypic factors. Eur J Cancer. 2005, 41 (14): 2040-2059. 10.1016/j.ejca.2005.03.034.
Newton-Bishop J, Gruis N: Melanoma susceptibility genes. Melanoma Res. 2010, 20 (3): 161-162. 10.1097/CMR.0b013e328336b000.
Law MH, Macgregor S, Hayward NK: Melanoma genetics: recent findings take us beyond well-traveled pathways. J Invest Dermatol. 2012, 132 (7): 1763-1774. 10.1038/jid.2012.75.
Savage SA, Gerstenblith MR, Goldstein AM, Mirabello L, Fargnoli MC, Peris K, Landi MT: Nucleotide diversity and population differentiation of the melanocortin 1 receptor gene, MC1R. BMC Genet. 2008, 9: 31-
Gerstenblith MR, Goldstein AM, Fargnoli MC, Peris K, Landi MT: Comprehensive evaluation of allele frequency differences of MC1R variants across populations. Hum Mutat. 2007, 28 (5): 495-505. 10.1002/humu.20476.
Raimondi S, Sera F, Gandini S, Iodice S, Caini S, Maisonneuve P, Fargnoli MC: MC1R variants, melanoma and red hair color phenotype: a meta-analysis. Int J Cancer. 2008, 122 (12): 2753-2760. 10.1002/ijc.23396.
Kanetsky PA, Panossian S, Elder DE, Guerry D, Ming ME, Schuchter L, Rebbeck TR: Does MC1R genotype convey information about melanoma risk beyond risk phenotypes?. Cancer. 2010, 116 (10): 2416-2428.
Williams PF, Olsen CM, Hayward NK, Whiteman DC: Melanocortin-1-receptor and risk of cutaneous melanoma: a meta-analysis and estimates of population burden. Int J Cancer. 2011, 129 (7): 1730-1740. 10.1002/ijc.25804.
Cust AE, Goumas C, Holland EA, Agha-Hamilton C, Aitken JF, Armstrong BK, Giles GG, Kefford RF, Schmid H, Hopper JL, et al: MC1R Genotypes and risk of melanoma before age 40 years: a population-based case–control-family study. Int J Cancer. 2012, 131: E269-E281. 10.1002/ijc.27357.
Whiteman DC, Green AC: A risk prediction tool for melanoma?. Cancer Epidemiol Biomarkers Prev. 2005, 14 (4): 761-763. 10.1158/1055-9965.EPI-14-4-ED.
Smith LA, Qian M, Ng E, Shao Y, Berwick M, Lazovich D, Polsky D: Development of a melanoma risk prediction model incorporating MC1R genotype and indoor tanning exposure. J Clin Oncol. 2012, 30 (suppl; abstr 8574):
Stefanaki I, Panagiotou OA, Kodela E, Gogas H, Kypreou KP, Chatzinasiou F, Nikolaou V, Plaka M, Kalfa I, Antoniou C, et al: Replication and predictive value of SNPs associated with melanoma and pigmentation traits in a Southern European case–control study. PLoS One. 2013, 8 (2): e55712-10.1371/journal.pone.0055712.
Cust AE, Schmid H, Maskiell JA, Jetann J, Ferguson M, Holland EA, Agha-Hamilton C, Jenkins MA, Kelly J, Kefford RF, et al: Population-based, case–control-family design to investigate genetic and environmental influences on melanoma risk: australian melanoma family study. Am J Epidemiol. 2009, 170 (12): 1541-1554. 10.1093/aje/kwp307.
Bishop DT, Demenais F, Iles MM, Harland M, Taylor JC, Corda E, Randerson-Moor J, Aitken JF, Avril MF, Azizi E, et al: Genome-wide association study identifies three loci associated with melanoma risk. Nat Genet. 2009, 41 (8): 920-925. 10.1038/ng.411.
Cust AE, Jenkins MA, Goumas C, Armstrong BK, Schmid H, Aitken JF, Giles GG, Kefford RF, Hopper JL, Mann GJ: Early-life sun exposure and risk of melanoma before age 40 years. Cancer Causes Control. 2011, 22 (6): 885-897. 10.1007/s10552-011-9762-3.
English DR, MacLennan R, Rivers J, Kelly J, Armstrong BK: Epidemiological studies of melanocytic naevi: protocol for identifying and recording naevi. 1990, Lyon (France): International Agency for Research on Cancer (IARC) internal report No. 90/002
Dwyer T, Blizzard L, Ashbolt R: Sunburn associated with increased number of nevi in darker as well as lighter skinned adolescents of northern European descent. Cancer Epidemiol Biomarkers Prev. 1995, 4 (8): 825-830.
Weatherall IL, Coombs BD: Skin color measurements in terms of CIELAB color space values. J Invest Dermatol. 1992, 99 (4): 468-473. 10.1111/1523-1747.ep12616156.
Creech LL, Mayer JA: Ultraviolet radiation exposure in children: a review of measurement strategies. Ann Behav Med. 1997, 19 (4): 399-407. 10.1007/BF02895159.
Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, Pencina MJ, Kattan MW: Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. 2010, 21 (1): 128-138. 10.1097/EDE.0b013e3181c30fb2.
Spitz MR, Amos CI, D'Amelio A, Dong Q, Etzel C: Re: Discriminatory accuracy from single-nucleotide polymorphisms in models to predict breast cancer risk. J Natl Cancer Inst. 2009, 101 (24): 1731-1732. 10.1093/jnci/djp394. author reply 1732
Pepe MS, Gu JW, Morris DE: The potential of genes and other markers to inform about risk. Cancer Epidemiol Biomarkers Prev. 2010, 19 (3): 655-665. 10.1158/1055-9965.EPI-09-0510.
Pencina MJ, D'Agostino RB, D'Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27 (2): 157-172. 10.1002/sim.2929. discussion 207–112
Steyerberg EW: Clinical prediction models: a practical approach to development, validation, and updating. 2009, New York, USA: Springer
Steyerberg EW, Pencina MJ, Lingsma HF, Kattan MW, Vickers AJ, Van Calster B: Assessing the incremental value of diagnostic and prognostic markers: a review and illustration. Eur J Clin Invest. 2012, 42 (2): 216-228. 10.1111/j.1365-2362.2011.02562.x.
Newton-Bishop JA, Chang YM, Iles MM, Taylor JC, Bakker B, Chan M, Leake S, Karpavicius B, Haynes S, Fitzgibbon E, et al: Melanocytic nevi, nevus genes, and melanoma risk in a large case–control study in the United Kingdom. Cancer Epidemiol Biomarkers Prev. 2010, 19 (8): 2043-2054. 10.1158/1055-9965.EPI-10-0233.
English DR, Armstrong BK: Identifying people at high risk of cutaneous malignant melanoma: results from a case–control study in Western Australia. Br Med J (Clin Res Ed). 1988, 296 (6632): 1285-1288. 10.1136/bmj.296.6632.1285.
Fears TR, Guerry D, Pfeiffer RM, Sagebiel RW, Elder DE, Halpern A, Holly EA, Hartge P, Tucker MA: Identifying individuals at high risk of melanoma: a practical predictor of absolute risk. J Clin Oncol. 2006, 24 (22): 3590-3596. 10.1200/JCO.2005.04.1277.
Williams LH, Shors AR, Barlow WE, Solomon C, White E: Identifying persons at highest risk of melanoma using self-assessed risk factors. J Clin Exp Dermatol Res. 2011, 2 (6):
Guther S, Ramrath K, Dyall-Smith D, Landthaler M, Stolz W: Development of a targeted risk-group model for skin cancer screening based on more than 100,000 total skin examinations. J Eur Acad Dermatol Venereol. 2012, 26 (1): 86-94. 10.1111/j.1468-3083.2011.04014.x.
Cho E, Rosner BA, Feskanich D, Colditz GA: Risk factors and individual probabilities of melanoma for whites. J Clin Oncol. 2005, 23 (12): 2669-2675.
Harbauer A, Binder M, Pehamberger H, Wolff K, Kittler H: Validity of an unsupervised self-administered questionnaire for self-assessment of melanoma risk. Melanoma Res. 2003, 13 (5): 537-542. 10.1097/00008390-200310000-00013.
Barbini P, Cevenini G, Rubegni P, Massai MR, Flori ML, Carli P, Andreassi L: Instrumental measurement of skin colour and skin type as risk factors for melanoma: a statistical classification procedure. Melanoma Res. 1998, 8 (5): 439-447. 10.1097/00008390-199810000-00009.
Fortes C, Mastroeni S, Bakos L, Antonelli G, Alessandroni L, Pilla MA, Alotto M, Zappala A, Manoorannparampill T, Bonamigo R, et al: Identifying individuals at high risk of melanoma: a simple tool. Eur J Cancer Prev. 2010, 19 (5): 393-400. 10.1097/CEJ.0b013e32833b492f.
Quereux G, Moyse D, Lequeux Y, Jumbou O, Brocard A, Antonioli D, Dreno B, Nguyen JM: Development of an individual score for melanoma risk. Eur J Cancer Prev. 2011, 20 (3): 217-224. 10.1097/CEJ.0b013e32834474ae.
Hemminki K, Zhang H, Czene K: Familial and attributable risks in cutaneous melanoma: effects of proband and age. J Invest Dermatol. 2003, 120 (2): 217-223. 10.1046/j.1523-1747.2003.12041.x.
Cust AE, Armstrong BK, Goumas C, Jenkins MA, Schmid H, Hopper JL, Kefford RF, Giles GG, Aitken JF, Mann GJ: Sunbed use during adolescence and early adulthood is associated with increased risk of early-onset melanoma. Int J Cancer. 2011, 128 (10): 2425-2435. 10.1002/ijc.25576.
Pepe MS, Janes H, Longton G, Leisenring W, Newcomb P: Limitations of the odds ratio in gauging the performance of a diagnostic, prognostic, or screening marker. Am J Epidemiol. 2004, 159 (9): 882-890. 10.1093/aje/kwh101.
Ware JH: The limitations of risk factors as prognostic tools. N Engl J Med. 2006, 355 (25): 2615-2617. 10.1056/NEJMp068249.
Galea S, Tracy M: Participation rates in epidemiologic studies. Ann Epidemiol. 2007, 17 (9): 643-653. 10.1016/j.annepidem.2007.03.013.
Cardis E, Richardson L, Deltour I, Armstrong B, Feychting M, Johansen C, Kilkenny M, McKinney P, Modan B, Sadetzki S, et al: The INTERPHONE study: design, epidemiological methods, and description of the study population. Eur J Epidemiol. 2007, 22 (9): 647-664. 10.1007/s10654-007-9152-z.
Kanetsky PA, Rebbeck TR, Hummer AJ, Panossian S, Armstrong BK, Kricker A, Marrett LD, Millikan RC, Gruber SB, Culver HA, et al: Population-based study of natural variation in the melanocortin-1 receptor gene and melanoma. Cancer Res. 2006, 66 (18): 9330-9337. 10.1158/0008-5472.CAN-06-1634.
Pearson TA, Manolio TA: How to interpret a genome-wide association study. JAMA. 2008, 299 (11): 1335-1344. 10.1001/jama.299.11.1335.
Collins FS, Green ED, Guttmacher AE, Guyer MS, Institute USNHGR: A vision for the future of genomics research. Nature. 2003, 422 (6934): 835-847. 10.1038/nature01626.
Gail MH: Value of adding single-nucleotide polymorphism genotypes to a breast cancer risk model. J Natl Cancer Inst. 2009, 101 (13): 959-963. 10.1093/jnci/djp130.
Freedman AN, Seminara D, Gail MH, Hartge P, Colditz GA, Ballard-Barbash R, Pfeiffer RM: Cancer risk prediction models: a workshop on development, evaluation, and application. J Natl Cancer Inst. 2005, 97 (10): 715-723. 10.1093/jnci/dji128.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/406/prepub
This work was supported by the National Health and Medical Research Council of Australia (NHMRC) (project grants 566946, 107359, 211172 and program grant number 402761 to GJM and RFK); the Cancer Council New South Wales (project grant 77/00, 06/10), the Cancer Council Victoria and the Cancer Council Queensland (project grant 371); the US National Institutes of Health (via NIH RO1 grant CA-83115-01A2 to the international Melanoma Genetics Consortium - GenoMEL) and a Victorian Cancer Agency Early Career Seed Grant (ECSG07_010). AEC is the recipient of a NHMRC public health postdoctoral fellowship (520018) and a Cancer Institute NSW Early Career Development Fellowship (10/ECF/2-06). BKA’s research is supported by a University of Sydney Medical Foundation Program Grant and JLH is an Australia Fellow of the NHMRC. For the English case–control study (Melanoma Cohort Study), the collection of samples was funded by Cancer Research UK (Project Grant C8216/A6129 and Program awards C588/A4994 and C588/A10589) and by the NIH (R01 CA83115). None of the funding bodies had a role in the design, collection, analysis, or interpretation of data, in the writing of the manuscript, or the decision to submit the manuscript for publication.
We gratefully acknowledge all of the participants, the work and dedication of the research coordinators, interviewers, examiners and data management staff, including Judith Maskiell, Jackie Arbuckle, Steven Columbus, Michaela Lang, Helen Rodais, Caroline Ellis (Centre for MEGA Epidemiology, School of Population Health, University of Melbourne, Melbourne, Australia); Carol El Hayek, Lynne Morgan, Joanne Roland, Emma Tyler, Jodi Barton, Caroline Watts and Lesley Porter (Westmead Institute of Cancer Research, University of Sydney at Westmead Millennium Institute and Melanoma Institute Australia, Sydney, Australia); Jodie Jetann, Megan Ferguson, Michelle Hillcoat, Kellie Holland, Pamela Saunders, Joan Roberts and Sheree Tait (Viertel Centre for Research in Cancer Control, Cancer Council Queensland, Spring Hill, Brisbane, Australia); Anil Kurien, Clare Patterson, Caroline Thoo, Sally de Zwaan, Angelo Sklavos, Shobhan Manoharan, Jenny Cahill and Sarah Brennand (skin examiners). In the UK, recruitment was facilitated by the UK National Cancer Research Network. We gratefully acknowledge the work of all research coordinators, interviewers, data management and laboratory staff involved in the study.
The authors declare that they have no competing interests.
AEC conceived the study idea, participated in its design, coordinated the statistical analysis, and drafted the manuscript. CG conducted the statistical analysis. MAJ and KV participated in the design of the study and helped to draft the manuscript. JRD and JHB conducted the statistical analysis for the external validation. EAH and CA participated in the MC1R sequencing. HS, BKA, RFK, JFA, GGG, DTB, JAN, JLH, and GJM coordinated the acquisition of data. All authors revised the manuscript critically for important intellectual content, and read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Box-and-whisker plots showing the predicted probabilities of melanoma for cases and controls for the A. Base, B. Self-reported and C. Physician-measured final models. The discrimination slope is calculated as the difference between the mean predicted probability of melanoma for cases and controls. The box represents the median and interquartile range, and the bars indicate the range. The base model includes demographic factors age, sex, city of recruitment and European ancestry. Both the self-reported model and the physician-measured model also include MC1R, non-melanoma skin cancer and nevi (self-reported model = none, few, some, many; physician-measured model = number of nevi ≥ 2 mm). (PNG 18 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cust, A.E., Goumas, C., Vuong, K. et al. MC1Rgenotype as a predictor of early-onset melanoma, compared with self-reported and physician-measured traditional risk factors: an Australian case-control-family study. BMC Cancer 13, 406 (2013). https://doi.org/10.1186/1471-2407-13-406