Abstract
Recent studies, including genome-wide association studies, have identified several putative low-penetrance susceptibility loci for melanoma. We sought to determine their generalizability to genetic predisposition for multiple primary melanoma in the international population-based Genes, Environment, and Melanoma (GEM) Study. GEM is a case–control study of 1,206 incident cases of multiple primary melanoma and 2,469 incident first primary melanoma participants as the control group. We investigated the odds of developing multiple primary melanoma for 47 SNPs from 21 distinct genetic regions previously reported to be associated with melanoma. ORs and 95% confidence intervals were determined using logistic regression models adjusted for baseline features (age, sex, age by sex interaction, and study center). We investigated univariable models and built multivariable models to assess independent effects of SNPs. Eleven SNPs in 6 gene neighborhoods (TERT/CLPTM1L, TYRP1, MTAP, TYR, NCOA6, and MX2) and a PARP1 haplotype were associated with multiple primary melanoma. In a multivariable model that included only the most statistically significant findings from univariable modeling and adjusted for pigmentary phenotype, back nevi, and baseline features, we found TERT/CLPTM1L rs401681 (P = 0.004), TYRP1 rs2733832 (P = 0.006), MTAP rs1335510 (P = 0.0005), TYR rs10830253 (P = 0.003), and MX2 rs45430 (P = 0.008) to be significantly associated with multiple primary melanoma, while NCOA6 rs4911442 approached significance (P = 0.06). The GEM Study provides additional evidence for the relevance of these genetic regions to melanoma risk and estimates the magnitude of the observed genetic effect on development of subsequent primary melanoma. Cancer Epidemiol Biomarkers Prev; 24(6); 992–7. ©2015 AACR.
This article is featured in Highlights of This Issue, p. 887
Introduction
Clinically, melanoma is one of a small group of cancers where patients are at increased risk of potentially life-threatening subsequent primaries (1, 2), but the underlying genetic predispositions to multiple primaries are relatively unexplored. Recent genome-wide association studies (GWAS) and candidate pathways studies have identified several low-penetrant genetic variants associated with cutaneous melanoma (3, 4). The majority of these variants are in gene regions associated with fair pigmentation, such as TYRP1, TYR, HERC2/OCA2, SLC45A2, and ASIP; nevi, such as PLA2G6, MTAP, and NID1; or both, such as IRF4 (4–12). More recent GWAS have identified melanoma risk–associated variants in genes, including ATM, MX2, PARP1, ARNT, and CASP8, which may not be associated with phenotypic risk (8, 13). However, the risk associated with these low-penetrant genetic variants in relation to multiple primary melanomas (MPM) has rarely been evaluated.
We studied these variants in the Genes, Environment and Melanoma (GEM) Study, a large, international population-based case–control study, in which the “cases” are patients with MPM and the “controls” are patients with single primary melanoma (SPM; refs. 14, 15). Participants' germline DNAs were genotyped for 47 polymorphisms from 21 distinct genomic regions. We compared the odds of carrying the genotypes and haplotypes in MPM relative to SPM patients in univariable and multivariable analyses and assessed effect modification.
Materials and Methods
Study population
The GEM Study is a population-based case–control study that enrolled 1,206 cases diagnosed with MPM (a second or higher-order invasive or in situ primary melanoma) between 1998 and 2003 and 2,469 controls diagnosed with invasive SPM in 2000. In situ melanomas were eligible as MPM in order to take into account surveillance when the patient had a previous invasive melanoma. Patients were recruited from eight population-based cancer registries in the United States (New Jersey, North Carolina, California), Australia (New South Wales, Tasmania), Canada (Ontario, British Columbia), and Italy (Turin), and one hospital center in Michigan. GEM recruitment procedures and data collection have been described (14, 15). The Institutional Review Boards of all participating institutions approved the protocol; informed consent was obtained from each participant.
Demographic and melanoma risk factors, including hair and eye color, ability to tan, and number of back nevi, were collected from telephone interview and self-administered questionnaire. Using a glossy-colored guide to aid in differentiating between nevi and other skin lesions, subjects had the nevi on their backs counted by a family member or friend; and back nevi counts were categorized as 0 to 10 or >10 for this article. Back nevus counts were significantly correlated with whole-body nevus diagrams in GEM (data not shown). A phenotypic index variable was derived from additively combining: hair color (black or dark brown = 0; light brown or blond = 1; red = 2), eye color (black or brown = 0; green, gray, or hazel = 1; blue = 2), and ability to tan in response to sun exposure (deeply or moderately = 0; occasionally or none = 1). Those with index scores of 0 or 1 were defined as very low/low, 2 as medium low, 3 as medium high, and 4 or 5 as high/very high risk.
Genotyping
DNA was collected from buccal swab kits. SNPs were genotyped using the MassArray iPLEX platform (Sequenom Inc.) with quality control measures as previously reported (16). Two SNPs of interest were not compatible with the Sequenom design, and proxy SNPs (r2 > 0.95) were chosen (1000 Genomes, CEU population; Proxy SNP; Broad Institute).
Statistical analyses
Logistic regression models were used to estimate the OR and 95% confidence interval (95% CI) for each SNP assuming an additive model of inheritance of the variant allele. All models were adjusted for baseline features: age, sex, an age by sex interaction, and study center. For each locus with multiple associated SNPs, we applied stepwise logistic regression to determine the SNP with the strongest association from among the significantly associated SNPs, keeping baseline variables fixed.
For the genes with at least two SNPs genotyped, we determined their haplotype blocks using the Haploview software algorithm (17) based on the pair-wise linkage disequilibrium information of the GEM population in combination with the HapMap CEU population (18). Within each haplotype block, we inferred the haplotypes in terms of probabilities for each individual from the SNP genotype input data using the PHASE algorithm, a Bayesian method in which the prior was chosen to approximate the coalescent process (19). For each haplotype block, the haplotype associations with MPM were assessed through haplotype trend regressions (HTR; ref. 20). The HTR method effectively took into account the haplotype phase uncertainty and reduced bias by incorporating the inferred individual haplotype probabilities as the predictor variables in the regression models. Haplotypes with low estimated frequencies (<0.01) were grouped together, reducing the number of haplotype categories and increasing the efficiency and power of haplotype analysis. Each haplotype or grouped rare haplotypes were then compared with the most common haplotypes in our study population.
Genotype and haplotype associations with phenotypic index were estimated using multinomial models and with back nevi using logistic regression models. In subsequent analyses, we limited the participants to those with no missing data for phenotypic index, back nevus counts, genotypes, and haplotypes of interest. Baseline models for genotypes and haplotypes were adjusted for baseline features; we then also adjusted the models for phenotypic index and back nevi, and finally also included all genotypes and haplotypes of interest in a multivariable model. Further, in exploratory stratified analyses, we assessed effect measure modification by phenotypic index and back nevus counts. The likelihood ratio test was used to test interactions, comparing models with main effects to models with main effects and interaction terms. All statistical tests were two-sided with P < 0.05 considered statistically significant. All data were analyzed using R (http://www.r-project.org/) or SAS 9.3 programs.
Results
The SPM and MPM patients' age, sex, race, number of back nevi, phenotypic index, and tumors' Breslow thicknesses are in Supplementary Table S1. Twelve non-Caucasian patients were excluded from analyses. Forty-seven SNPs within 21 genetic loci previously reported to be in association with melanoma were genotyped. SNP locations, minor allele frequencies, numbers of cases and controls genotyped, and literature references are in Supplementary Table S2. Proxy SNPs rs6735656 and rs12278954 were used, respectively, for CASP8 rs10931936 and ATM rs1801516 identified by Barrett and colleagues (8).
Eleven SNPs in 6 gene neighborhoods (TERT/CLPTM1L, TYRP1, MTAP, TYR, NCOA6, and MX2) were significantly (P < 0.05) associated with MPM compared with SPM using additive models adjusted for baseline features (Table 1). The MTAP and TYR loci each had more than one significantly associated SNP. MTAP rs1335510 and TYR rs10830253 were brought forward for subsequent analyses because they each were the only SNP that remained in the stepwise logistic regression model for their locus.
Association of genotypes with MPM (N = 1,206) compared with SPM (N = 2,469) in the GEM Study
Of the haplotypes examined (Table 2 and Supplementary Table S3), a PARP1 haplotype (rs3219090, rs2695238) was significantly (P = 0.03) associated with MPM, as were haplotypes in MTAP and TYR, when adjusting for baseline features. However, the statistical significance of the MTAP and TYR haplotype associations with MPM was weaker than the respective single MTAP rs1335510 and TYR rs10830253 SNP associations (Supplementary Table S3). Thus, only the PARP1 haplotype was selected for further analysis.
Association of a PARP1 haplotype with MPM in the GEM Study
Several of the MPM-associated genotypes were associated with phenotypic index or back nevus counts (Supplementary Table S4), indicating that these SNPs are correlated with these phenotypes and may increase risk of MPM via these phenotypes, which are known melanoma risk factors.
When limiting the dataset to patients with no missing data for genotypes, haplotype, or traits of interest (Table 3), the TERT/CLPTM1L, TYRP1, MTAP, TYR, NCOA6, and MX2 genotypes, but not the PARP1 haplotype (P = 0.14), remained significantly associated with MPM after adjusting for baseline features. After additionally adjusting for phenotypic index and back nevi, the ORs did not appreciably change, although the association with NCOA6 rs4911442 became insignificant (P = 0.07). In a multivariable model further adjusting for genotypes and the PARP1 haplotype, the TERT/CLPTM1L, TYRP1, MTAP, TYR, and MX2 genotypes remained significant, but the NCOA6 genotype (P = 0.06) and PARP1 haplotype (P = 0.22) did not; none of the ORs appreciably changed.
Genotype and haplotype associations with MPM (N = 1,074) compared with SPM (N = 2,137) for participants with no missing data in the GEM Study
In an exploratory stratified analysis adjusted for baseline features, no evidence was found of effect modification by phenotypic index or back nevus counts on the association of genotypes with MPM (Supplementary Table S5). However, there was evidence of effect modification by back nevi, but not phenotypic index, on the association between the PARP1 haplotype and MPM (P for interaction = 0.01). The PARP1 haplotype AG was negatively associated (P = 0.02) with MPM when 0 to 10 back nevi were present, whereas both the AG (P = 0.03) and AC (P = 0.01) haplotypes were negatively associated with MPM when >10 back nevi were present.
Discussion
In the international GEM Study, we found that SNPs in TERT/CLPTM1L, TYRP1, MTAP, TYR, NCOA6, and MX2 and a PARP1 haplotype were associated with the occurrence of MPM. TERT/CLPTM1L rs401681, TYRP1 rs1408799, MTAP rs1335510, TYR rs10830253, and MX2 rs45430 were associated with MPM independently of each other and of phenotypic index and back nevi. NCOA6 rs4911442 and the PARP1 haplotype were not significant in the multivariable model, possibly as a result of diminished statistical power as there was little change in the ORs. There was no evidence for effect modification of SNP associations with MPM by patient phenotype; however, back nevi did modify the association of the PARP1 haplotype with MPM.
The single SNP associations reported in GEM are in the same direction as those reported in the literature (Supplementary Table S2). In a recent large meta-analysis, Chatzinasiou and colleagues (21) found variants in 8 of the 21 low-penetrant loci examined in our study to have strong epidemiologic credibility (Venice criteria overall grade, A) as associated with melanoma, although this study did not specify whether the associations were specific to MPM, SPM, or both. Variants in each of these loci reached significance in our study except for PIGU rs910873 (P = 0.07), MYH7B rs1885120 (P = 0.11), and SLC45A2 rs16891982 (P = 0.14), the ORs of which were in the same direction as previously reported. We are not aware of another group examining the PARP1 rs3219090, rs2695238 haplotype. Pena-Chilet and colleagues (22) found a similar trend toward protection from melanoma with a PARP1 haplotype containing the minor alleles of rs1136410 and rs3219090.
We know of only one other study addressing associations with MPM of common genetic variants in the loci discussed here. Helsing and colleagues (23) found no association of ASIP rs1015362 and rs4911414, TYR rs1126809, and TYRP1 rs1408799 with MPM in patients identified through the Norwegian Cancer Registry compared with melanoma-free blood donors. Although this design differs from GEM, we also found no association of ASIP rs4911414 and TYRP1 rs1408799 with MPM. We did not genotype ASIP rs1015362 and TYR rs1126809.
Strengths of our study include its large size, population-based case ascertainment, homogeneous questionnaire administration with respect to cases and controls, phenotypic trait adjustments, multivariable model, and examination of risk stratified by phenotypic traits. There may be a tendency for effect estimates of individual risk factors to be attenuated in our high baseline risk population (15). As evidence, several genotypes previously reported as associated with melanoma risk (PARP1 rs3219090, CASP8 rs13016963, TYRP1 rs2733832, and PIGU rs910873) had borderline associations (0.10 < P ≥ 0.05) in the same direction as the literature in GEM. However, the GEM Study estimates are more relevant for survivors and their risk of subsequent melanoma than studies investigating lower-risk populations. A limitation is that GEM may not have had sufficient power to detect associations of SNPs with lower minor allele frequencies (e.g., SLC45A2 rs16891982, MAF = 0.017).
Subsequent melanomas are a major problem for melanoma patients, but few studies have explored their genetic predisposition. Our results provide evidence that several putative low-penetrance susceptibility loci for melanoma are generalizable to risk of subsequent melanoma. Also, validation of genetic associations in the large international population-based GEM study adds further credibility that these loci are melanoma risk–associated. Knowledge of genetic risk factors for subsequent melanoma could inform screening algorithms, future risk estimation modeling, and future prevention studies for melanoma survivors.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: D.C. Gibbs, I. Orlow, A. Kricker, B.K. Armstrong, H. Anton-Culver, S.B. Gruber, L.D. Marrett, R.P. Gallagher, S. Rosso, N.E. Thomas
Development of methodology: D.C. Gibbs, I. Orlow, H. Anton-Culver, L.D. Marrett, R.P. Gallagher, N.E. Thomas
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): I. Orlow, A. Kricker, B.K. Armstrong, H. Anton-Culver, S.B. Gruber, L.D. Marrett, R.P. Gallagher, S. Rosso, T. Dwyer, A. Sharma, E. La Pilla, L. From, K.J. Busam, D.W. Ollila, M. Berwick, N.E. Thomas
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D.C. Gibbs, I. Orlow, P.A. Kanetsky, L. Luo, R. Zanetti, A.E. Cust, C.B. Begg, M. Berwick, N.E. Thomas
Writing, review, and/or revision of the manuscript: D.C. Gibbs, I. Orlow, P.A. Kanetsky, L. Luo, A. Kricker, B.K. Armstrong, H. Anton-Culver, S.B. Gruber, R.P. Gallagher, R. Zanetti, T. Dwyer, A. Sharma, E. La Pilla, L. From, K.J. Busam, A.E. Cust, D.W. Ollila, C.B. Begg, M. Berwick, N.E. Thomas
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): D.C. Gibbs, R.P. Gallagher, M. Berwick, N.E. Thomas
Study supervision: D.C. Gibbs, B.K. Armstrong, H. Anton-Culver, R.P. Gallagher, S. Rosso, M. Berwick, N.E. Thomas
Grant Support
This study was supported by NCI grants R01CA112243 (to N.E. Thomas), R01CA112524 (to M. Berwick), R01CA112243-05S1 (to N.E. Thomas), R01CA112524-05S2 (to M. Berwick), CA098438 (to C.B. Begg), U01CA83180 (to M. Berwick), R33CA10704339 (to N.E. Thomas), P30CA016086 (to S.H. Earp), P30CA014089 (to S.B. Gruber), and P30CA008748 (to C. Thompson); National Institute of Environmental Health Sciences (P30ES010126, to J. Swenberg); University of Sydney Medical Foundation Program grant (to B.K. Armstrong); and Michael Smith Foundation for Health Research Infrastructure Award (to R.P. Gallagher).
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.
Acknowledgments
The study was conducted by The GEM Study Group: Coordinating Center, Memorial Sloan-Kettering Cancer Center, New York, NY: Marianne Berwick, M.P.H., Ph.D. [Principal Investigator (PI), currently at the University of New Mexico], Colin B. Begg, Ph.D. (co-PI), Irene Orlow, Ph.D. (co-Investigator), Klaus J. Busam, M.D. (Dermatopathologist), Anne S. Reiner, M.P.H. (Biostatistician), Pampa Roy, Ph.D. (Laboratory Technician), Ajay Sharma, M.S. (Laboratory Technician), Emily La Pilla (Laboratory Technician). University of New Mexico, Albuquerque: Marianne Berwick, M.P.H., Ph.D. (PI), Li Luo, Ph.D. (Biostatistician), Kirsten White, MSc (Laboratory Manager), Susan Paine, M.P.H. (Data Manager). Study centers included the following: The University of Sydney and The Cancer Council New South Wales, Sydney, Australia: Bruce K. Armstrong, M.B.B.S., D.Phil. (PI), Anne Kricker, Ph.D. (co-PI), Anne E. Cust, Ph.D. (co-Investigator); Menzies Research Institute Tasmania, University of Tasmania, Hobart, Australia: Alison Venn, Ph.D. (current PI), Terence Dwyer, M.D. (PI, currently at University of Oxford, United Kingdom), Paul Tucker, M.D. (Dermatopathologist); British Columbia Cancer Research Centre, Vancouver, Canada: Richard P. Gallagher, M.A. (PI), Donna Kan (Coordinator); Cancer Care Ontario, Toronto, Canada: Loraine D. Marrett, Ph.D. (PI), Elizabeth Theis, M.Sc. (co-Investigator), Lynn From, M.D. (Dermatopathologist); CPO, Center for Cancer Prevention, Torino, Italy: Roberto Zanetti, M.D (PI), Stefano Rosso, M.D., M.Sc. (co-PI); University of California, Irvine, CA: Hoda Anton-Culver, Ph.D. (PI), Argyrios Ziogas, Ph.D. (Statistician); University of Michigan, Ann Arbor, MI: Stephen B. Gruber, M.D., M.P.H., Ph.D. (PI, currently at University of Southern California, Los Angeles, CA), Timothy Johnson, M.D. (Director of Melanoma Program), Duveen Sturgeon, M.S.N. (co-Investigator, joint at USC-University of Michigan); University of North Carolina, Chapel Hill, NC: Nancy E. Thomas, M.D., Ph.D. (PI), Robert C. Millikan, Ph.D. (previous PI, deceased), David W. Ollila, M.D. (co-Investigator), Kathleen Conway, Ph.D. (co-Investigator), Pamela A. Groben, M.D. (Dermatopathologist), Sharon N. Edmiston, B.A. (Research Analyst), Honglin Hao (Laboratory Specialist), Eloise Parrish, MSPH (Laboratory Specialist), David C. Gibbs, B.S. (Research Assistant), Jill S. Frank, M.S. (Research Assistant), Jennifer I. Bramson (Research Assistant); University of Pennsylvania, Philadelphia, PA: Timothy R. Rebbeck, Ph.D. (PI), Peter A. Kanetsky, M.P.H., Ph.D. (co-Investigator, currently at H. Lee Moffitt Cancer Center & Research Institute, Tampa, Florida); UV data consultants: Julia Lee Taylor, Ph.D. and Sasha Madronich, Ph.D., National Centre for Atmospheric Research, Boulder, CO.
Footnotes
Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).
- Received December 30, 2014.
- Revision received March 9, 2015.
- Accepted March 13, 2015.
- ©2015 American Association for Cancer Research.