| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Epidemiology and Prevention |
1 Department of Biostatistics and Epidemiology and Center for Clinical Epidemiology and Biostatistics, University of Pennsylvania, Philadelphia, Pennsylvania; 2 Department of Epidemiology and Biostatistics, Memorial Sloan-Kettering Cancer Center, New York, New York; 3 University of Sydney, Sydney, New South Wales, Australia; 4 Cancer Care Ontario; 5 Women's College Hospital, Toronto, Ontario, Canada; 6 University of North Carolina, Chapel Hill, North Carolina; 7 University of Michigan, Ann Arbor, Michigan; 8 University of California, Irvine, California; 9 Centro per la Prevenzione Oncologia Torino, Piemonte, Italy; 10 British Columbia Cancer Agency, Vancouver, British Columbia, Canada; 11 Murdoch Childrens Research Institute, Royal Children's Hospital, Melbourne, Victoria, Australia; 12 New Jersey Department of Health and Senior Services, Trenton, New Jersey; and 13 Department of Epidemiology, University of New Mexico, Albuquerque, New Mexico
Requests for reprints: Peter A. Kanetsky, Center for Clinical Epidemiology and Biostatistics and Epidemiology, University of Pennsylvania, 903 Blockley Hall, Philadelphia, PA 19104-6021. Phone: 215-573-3282; Fax: 215-573-2260; E-mail: pkanetsk{at}cceb.med.upenn.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Further evidence of the role of MC1R variants in melanoma etiology comes from studies of melanoma-prone families that show a younger average age of onset of melanoma in families inheriting mutations in both MC1R and cyclin-dependent kinase inhibitor 2A (CDKN2A) than in families with mutations in CDKN2A alone (1013).
In this study, we examine the frequencies of MC1R variants and their associations with melanoma in a large population-based investigation. The study involved population-based recruitment of patients with an incident second primary melanoma or a higher-order primary melanoma (i.e., third, fourth, etc.). The "control" group represents the population at risk for these subsequent primary melanomas (i.e., patients who are diagnosed with a first primary melanoma). This design enriches the sampling for rare variants that are associated with the disease (14). The validity of the design for evaluating risk factors for melanoma has been examined in detail in a recent article (15). By sequencing the entire MC1R coding region in this large population-based sample, the study provides the most comprehensive data on MC1R in melanoma to date.
| Materials and Methods |
|---|
|
|
|---|
MC1R genotyping. All MC1R genotypes were determined at a single center as previously described in Kanetsky et al. (17). Genomic DNA was isolated from buccal cells using either a modified version of the Richards' protocol (18) or by Puregene DNA Purification kits for buccal cells (Gentra Systems, Minneapolis, MN). DNA was isolated from blood lymphocytes using Qiagen QIAamp DNA Blood kit (Qiagen, Inc., Valencia, CA). Using a MJ Research (Waltham, MA) PTC-100 thermocycler, standard PCR technique was used to amplify the entire 951-nucleotide MC1R coding region. All amplified products were directly sequenced on an ABI Prism 377 or 3100 (Applied Biosystems, Foster City, CA) using BigDye Terminators (Applied Biosystems) according to the manufacturer's specifications. Sequencing primers consisted of either a set of two or four oligonucleotides: 5'-GCCATGAGCACCAGCATAG-3' and 5'-GACCACACAAATATCACCACCT or 5'-TCGTCTTCAGCACGCTCTTC-3', 5'-TTTAAGGCCAAAGCCCTGGT-3', 5'-AACCTGCACTCACCCATGTA-3', and 5'-CTGCAGGTGATCACGTCAAT, respectively. MC1R chromatograms were read with the aid of Sequencher software version 4.0.5 (Gene Codes Corp., Ann Arbor, MI) and/or SeqScape software version 1.0 to 2.1.1 (Applied Biosystems). A minority of novel variants observed only once were verified independently. All MC1R genotypes were double entered into a customized Microsoft Office Access 2003 database before delivery to the coordinating center for statistical analysis.
Data coding. Information obtained from the self-administered questionnaire instrument (19) was used to create a phenotypic index that was based on participant responses to hair color (black or dark brown hair = 1; light brown or blond hair = 2; red hair = 3), eye color (black or brown eyes = 0; hazel, green, gray, or blue eyes = 1), and inability to tan in response to sun exposure (no = 0; yes = 1). An index score of 1 or 2 signifies overall darker cutaneous phenotype and indicates a low phenotypic risk, an index score of 3 indicates medium phenotypic risk, and an index score of 4 or 5 signifies overall fairer cutaneous phenotype and indicates high phenotypic risk. Using a glossy colored guide to aid in differentiating between nevi and other skin lesions, subjects were asked to have the nevi on their backs counted by a family member or friend; the number of self-reported nevi on the back was categorized as <5, 5 to 10, 11 to 25, and >26.
We calculated allele frequencies for the MC1R consensus sequence (Genbank accession no. AF326275) and observed MC1R variants and used the nomenclature and definitions suggested by Sturm et al. (9) to group MC1R variants as higher-risk R variants (D84E, R151C, R160W, and D294H) or lower-risk r variants (all other variants excluding synonymous changes), acknowledging that these risk categories are inexact because the precise functional status of many MC1R variants is unknown. We also categorized MC1R R and r variants using two alternative coding schemes: the first additionally included the R142H and the g.86_87insA variant (along with the other four observed insertion/deletion variants) as R variants based on prior work, showing the importance of these variants in pigmentation or MC1R function (5, 20, 21), and the second used previously published results from a Sorting Intolerant from Tolerant analysis, which identified MC1R protein positions that are predicted intolerant to amino acid substitutions, thus indicating those variants expected to be higher-risk based on putative functional importance (17).
Statistical analysis. Contingency table analysis was used to compare univariate associations [odds ratios (OR)] of individual or categories of MC1R variants and case/control status. Unconditional logistic regression was used to obtain adjusted estimates. These models were adjusted for age of diagnosis of most recent melanoma (continuous variable), sex, study center (eight indicator variables), and a term for the interaction of age and sex. Adjusted ORs (aOR) and corresponding 95% confidence intervals (95% CI) are presented.
| Results |
|---|
|
|
|---|
|
|
There were three instances of an ambiguous translation from the nucleotide change to corresponding amino acid designation because three individuals carried two variants that fell within the same codon. Although it is not possible with genotyping data alone to determine the phase of these nucleotide changes (i.e., whether they occur in cis or trans orientation), analysis of inferred MC1R haplotypes in our study population indicates that, with few exceptions, MC1R variants are inherited individually on an otherwise genetic background of the consensus sequence (results not presented). Therefore, we assumed that these nucleotide changes were in trans orientation and assigned amino acid changes accordingly. The first individual carried the R151C (g.451C>T) and R151R (g.453C>G) variants, the second carried the R160W (g.478C>T) and R160Q (g.479G>A) variants, and the third carried the T272A (g.814A>G) and T272K (g.815C>A) variants.
Overall, 86% of individuals carried at least one MC1R variant; 84% of all individuals carried nonsynonymous and insertion/deletion variants. Carriage of synonymous variants alone is an uncommon event. One thousand four hundred ten individuals (43.8%) carried one, 1,288 (40.0%) carried two, and 6 (0.19%) carried three nonsynonymous or insertion/deletion variants.
MC1R variants and melanoma. Table 3 shows the crude and aOR for associations between MC1R variants and melanoma. Despite our large sample size, there were few homozygous carriers of individual variants or compound heterozygous carriers with rare variants. Hence, for variant-specific associations, we present the OR for carriage of at least one variant versus homozygous carriage of the MC1R consensus. There was no association between the V60L, V92M, I155T, and R163Q variants and case status as evident by aORs near 1.0 and nonsignificant 95% CIs. In contrast, cases were more likely to carry the R variants and the R142H variant compared with controls as indicated by modest aORs in the range of 1.2 to 1.4 and statistically or borderline significant 95% CIs. In aggregate, case status was not associated with carriage of nonsynonymous variants that occur at a MAF <0.01 (aOR, 0.94; 95% CI, 0.60-1.5) nor with carriage of insertion/deletion variants (aOR, 0.96; 95% CI, 0.51-1.8).
|
|
MC1R and melanoma by levels of pigmentation and nevi. We did stratified analysis of MC1R-melanoma associations based on categories of phenotypic risk (low, medium, and high) as determined from the phenotypic index. With few exceptions, most stratum-specific associations did not reach statistical significance (Fig. 1A ). Comparing the three phenotypic index groups, the pattern of aORs across MC1R genotype categories seemed similar.
|
| Discussion |
|---|
|
|
|---|
This study represents by far the largest investigation of MC1R variants and melanoma to date. Because all participants in this study have melanoma, the observed MC1R allele frequencies are higher than previous reports among populations of healthy individuals yet are generally consistent with allele frequencies reported among other melanoma case groups (5, 7). We detected 46 previously unreported variants in MC1R; all were rare. Although we had a very large sample size, we had small numbers of rare variants. Thus, their contribution to melanoma risk was difficult to evaluate and is unlikely to be large.
Although our associations between MC1R variants and melanoma are consistent with previous studies, there are notable differences in the strengths of association observed in our study compared with those observed by other investigators. In most previously published work, associations between MC1R variants and melanoma are stronger than those found in our study. Palmer et al. (4) reported a 2.2- and 4.1-fold odds of melanoma among Australians who carry one or two copies of the R151C, R160W, and D294H variants. ORs reported by Kennedy et al. (5) among Dutch study participants were 3.1 and 4.9 for carriage of one or two copies, respectively, of any MC1R variant. A recent report among a French clinical population noted strong association between melanoma and carriage of one (OR, 4.3) or two (OR, 6.8) MC1R variants (7), whereas that from an Italian sample found that melanoma cases were over twice as likely to carry at least one MC1R and nearly four times more likely to carry two variants (8).
There are several methodologic differences between our study and other published investigations that may explain, in part, the difference in the observed strengths of association. Three studies used clinic- and/or hospital-based recruitment of melanoma cases for which individuals with early-stage melanoma may have been underrepresented (5, 7, 8). Indeed, a recent finding by Landi et al. (8) showed that the R151C, R160W, and D294H variants were more common among individuals with thicker melanomas (>1.06 mm). Thus, exclusion of thinner (early stage) lesions could lead to overestimation of MC1R associations. Further, in two studies, the case group consisted, in part, of individuals selected for increased likelihood of an underlying genetic susceptibility to melanoma (4, 7). This selection was severe in one study for which selected case subjects accounted for 73% of the case group (7). Specialized selection of case subjects may have lead to overestimates of true associations between MC1R and melanoma status.
We must, however, consider the possibility that our study design contributed to the possible underestimation of the true MC1R-melanoma effect. Our study design is adapted from traditional case-control methodology, in which both cases and controls arise from a theoretical joint source population. Here, cases (individuals with multiple melanomas) arise from a population of at-risk persons with incident single primary melanoma (controls). The rationale for implementing this design was motivated by the desire to create an efficient population-based study for investigating genetic variants in CDKN2A and CDK4 that are highly penetrant in the arena of familial melanoma (25) but rare at the population level (26). The degree to which this design may affect associations related to common exposures, including MC1R variants, is unknown. It is plausible that, among a population already at increased risk for melanoma (e.g., controls with single primary melanoma), the relative effect of exposures may be less than among a population at baseline risk (e.g., healthy controls with no history of melanoma), although the absolute effect on risk may be the same or greater (5).
Our stratified analysis by phenotypic index revealed that the overall pattern of association between melanoma and MC1R variants is not different across levels of cutaneous phenotype. Because of differences in the definition of "high-risk" MC1R variants, it was not possible to directly compare our results with those of Palmer et al. (4), Landi et al. (8), or Dwyer et al. (27), who found stronger associations between melanoma and MC1R variants among individuals with darker pigmentation.
For a subset of individuals, DNA did not amplify. Individuals whose DNA amplified (n = 3,218) compared with those whose DNA did not amplify (n = 366) were slightly older at their most recent melanoma diagnosis (60 versus 58 years; P = 0.05), more likely to be male (64% versus 55%; P < 0.001), and more likely from Australia (59% versus 31%) than from North America (37% versus 67%) or Italy (4% versus 2%; P < 0.001, for comparison of the three geographic regions). There was not a statistically significant difference between the groups for report of moles on the back (P = 0.06, for comparison of the four mole categories) or distribution of the phenotypic index (P = 0.18, for comparison of the three categories of phenotypic index). Importantly, we found no difference in genotyping success by case status (P = 0.12). Therefore, it is unlikely that the addition of these 366 individuals to our analyses would have had a significant effect on our point estimates.
Because our control group consisted of individuals with incident first primary melanoma rather than those with prevalent first primary melanoma, we also assessed the effect of MC1R variants on survival by examining time from ultimate to penultimate diagnosis of melanoma among multiple primary cases within individual variants and by carriage number (0 versus
1 variants). There were no statistically significant differences in the mean number of years between melanoma diagnoses, although difference for the R142H and D294H variants approached borderline statistical significance (P = 0.07 and 0.10, respectively).
This study of >3,200 individuals ascertained predominately from population-based sources and for whom complete genotype information at the MC1R locus was available is the largest molecular epidemiologic investigation of melanoma susceptibility to date. Although the point estimates of melanoma-MC1R associations were smaller than previously observed and our attenuated estimates may have affected our ability to definitively confirm differences in the association between MC1R variants and melanoma among individuals with different phenotypic profiles, it is noteworthy that our finding of a positive association between melanoma and MC1R variants and a trend toward increasingly stronger estimates with increasing MC1R risk category based on carriage number and variant (R or r) type is consistent with previous investigations.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
The Genes Environment and Melanoma Study Group includes the following. Coordinating center: Memorial Sloan-Kettering Cancer Center, Marianne Berwick (principal investigator, currently at the University of New Mexico, Albuquerque, NM), Colin B. Begg (coprincipal investigator), Irene Orlow (coinvestigator), Urvi Mujumdar (project coordinator), Amanda J. Hummer (biostatistician), Klaus Busam (dermatopathologist), Pampa Roy (laboratory technician), Rebecca Canchola (laboratory technician), Brian Clas (laboratory technician), Javiar Cotignola (laboratory technician), and Yvette Monroe (interviewer). Study centers: The University of Sydney and The Cancer Council New South Wales (Sydney, New South Wales, Australia), Bruce K. Armstrong (principal investigator), Anne Kricker (coprincipal investigator), and Melisa Litchfield (study coordinator); Menzies Centre for Population Health Research, University of Tasmania (Hobart, Tasmania, Australia), Terence Dwyer (principal investigator), Paul Tucker (dermatopathologist), and Nicola Stephens (study coordinator); British Columbia Cancer Agency (Vancouver, British Columbia, Canada), Richard P. Gallagher (principal investigator) and Teresa Switzer (coordinator); Cancer Care Ontario (Toronto, Ontario, Canada), Loraine D. Marrett (principal investigator), Elizabeth Theis (coinvestigator), Lynn From (dermatopathologist), Noori Chowdhury (coordinator), Louise Vanasse (coordinator), Mark Purdue (research officer), and David Northrup (manager for computer-assisted telephone interviewing); Centro per la Prevenzione Oncologia Torino (Piemonte, Italy), Roberto Zanetti (principal investigator), Stefano Rosso (data manager), and Carlotta Sacerdote (coordinator); University of California (Irvine, CA), Hoda Anton Culver (principal investigator), Nancy Leighton (coordinator), and Maureen Gildea (data manager); University of Michigan (Ann Arbor, MI), Stephen B. Gruber (principal investigator), Joe Bonner (data manager), and Joanne Jeter (coordinator); New Jersey Department of Health and Senior Services (Trenton, NJ), Judith Klotz (principal investigator), Homer Wilcox (coprincipal investigator), and Helen Weiss (coordinator); University of North Carolina (Chapel Hill, NC), Robert C. Millikan (principal investigator), Nancy Thomas (coinvestigator), Dianne Mattingly (coordinator), Jon Player (laboratory technician), and Chiu-Kit Tse (data analyst); University of Pennsylvania (Philadelphia, PA), Timothy R. Rebbeck (principal investigator), Peter A. Kanetsky (coinvestigator), Amy Walker (laboratory technician), and Saarene Panossian (laboratory technician). Consultants: Harvey Mohrenweiser (University of California at Irvine) and Richard Setlow (Brookhaven National Laboratory, Upton, NY).
Received 5/ 4/06. Revised 7/24/06. Accepted 7/24/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. K. Oleksyk, M. W. Smith, and S. J. O'Brien Genome-wide scans for footprints of natural selection Phil Trans R Soc B, January 12, 2010; 365(1537): 185 - 205. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |