Introduction
Melanoma-prone families are often defined as having two or more 1st-degree relatives with melanoma, with the number of cases per family varying according to the incidence rates in the respective populations and the presence of founder mutations (1). The germline mutations of two genes, the CDKN2A tumor suppressor gene on chromosome 9p21 (2) and the CDK4 oncogene on chromosome 12q14 (3), account for melanoma susceptibility in ∼20% to 25% of families with multiple cases of cutaneous malignant melanoma worldwide (4), and mutation frequency increases with the number of melanoma cases in the family (1). The estimated cumulative risk of developing melanoma over a lifetime in the Italian population is ∼0.5% (5), unlike in the United States, Australia, and New Zealand, where the risks are ∼2.0 (6), ∼3.3 (7), and ∼5.7% (8), respectively. Thus, a familial aggregation in a Mediterranean country like Italy is suggestive of an inherited predisposition (9). Moreover, melanoma in a population with a wide range of pigmentary phenotypes, small sized nevi, and intense sun exposure may reveal susceptibility pathways specific to this population. Therefore, we did linkage analysis in CDKN2A-negative and CDK4-negative melanoma-prone families from Italy to take advantage of the increased statistical power from the genetically enriched cluster of individuals to identify novel genetic loci linked to melanoma susceptibility.
Materials and Methods
Our study population included 51 CDKN2A-negative and CDK4-negative melanoma-prone families with 2 to 5 melanoma patients (ages 18-87 y) and consisted of 280 relatives with 85 affected individuals (10). Ten additional families from central and northern Italy, including 2 to 4 affected individuals per family (ages 26-82 y), were added for the final fine-scale mapping.
We used the Applied Biosystems Inc. linkage mapping set with 414 microsatellite markers at an average density of 10 cM. PCR amplification of genomic DNA and marker interrogation was conducted as previously reported (10). Mendelian inconsistencies were identified with the use of PedCheck (version 1.1), and inconsistent genotypes were removed. Fine mapping was done at deCODE Genetics with an average marker density of one marker every 1 cM.
Power calculation for linkage analysis based on these families was done with the use of SLINK (11, 12). Parametric 2-point and multipoint linkage analyses, assuming an autosomal dominant inheritance model, were done with the use of the MLINK program of the FASTLINK software package (version 4.1P) and GENEHUNTER (2.1 release 6), respectively. The disease allele frequency was calculated as 0.0004 based on the prevalent sex averaged rate of melanoma in the Italian population of 60 per 100,000 (13). An affected-only model was used to reduce the effect of incomplete penetrance of the trait on the analysis by coding all unaffected individuals as “unknown.” Nonparametric multipoint linkage analysis was done with GENEHUNTER.
The predivided sample test (14) was used to test for heterogeneity between families based on age at diagnosis. Families were classified as “early” (age of onset below the study population mean of 50.8 y), “late” (age of onset above the mean), or “mixed” (affected individuals in family with both early and late ages of onset).
Parametric 2-point linkage analysis was done under a dominant model, with a redefined affection status (i.e., the presence of dysplastic nevi regardless of melanoma status); all other subjects were coded as unknown. This analysis included 125 subjects with dysplastic nevi.
Results
Based on our initial 51 families, we had 85% power to detect a logarithm of the odds (LOD) score > 3 under a dominant model with 50% penetrance and the assumption of no heterogeneity, with the use of a simulated marker with six equally frequent alleles.
The plots of the 2-point parametric maximum LOD scores are shown in Fig. 1. The highest LOD score observed was 1.83 at marker D11S908 on chromosome 11. No significant P values were observed with the nonparametric analyses (data not shown).
To minimize genetic heterogeneity, we did a predivided sample test based on the mean age at diagnosis for affected individuals, hypothesizing that families with a lower age of onset would more likely have a genetic basis for their disease. There were 8, 27, and 16 families in the early, mixed, and late onset classes, respectively. The predivided sample test revealed no significant genetic heterogeneity among the groups (data not shown). The highest LOD score observed was 3.27 on chromosome 8 at marker D8S277 (8.34 cM) in the late onset classification of families. However, the flanking markers ∼7 and 13 cM away were strongly negative. Moreover, there was no haplotype segregating with the disease, and the melanoma's related alleles had similar frequencies in affected and unaffected individuals.
Linkage analyses with dysplastic nevi as the affection status showed no significant evidence for linkage (the maximum LOD score was 0.98 on chromosome 10).
We selected 5 regions on 4 chromosomes with LOD scores above 1.5 with flanking positive LOD scores for fine-scale mapping. The corresponding LOD scores are shown in Table 1. In this follow-up analysis, no significant evidence of linkage was observed in the expanded population.
Chromosome . | Number of markers . | Parametric 2-point LOD scores . | . | . | . | Nonparametric . | . | . | |||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | Maximum LOD . | 𝛉 . | Marker (cM*) . | Flanking markers . | Maximum NPL . | P . | cM* . | |||||
1 | 24 | 0.57 | 0.05 | D1S2850 (251.41) | Positive (one side) | 0.19337 | 0.42 | 0 | |||||
4 | 18 | 1.48 | 0.1 | D4S1581 (58.8) | Positive (both sides) | 1.667 | 0.05 | 70 | |||||
4 | 29 | 1.48 | 0.0 | DG4S1397 (65.2) | Positive (both sides) | ||||||||
7 | 10 | 1.33 | 0.1 | D7S661 (155) | Positive (one side) | 0.96 | 0.16 | 155 | |||||
11 | 39 | 1.06 | 0.1 | D11S908 (118.6) | Positive (one side) | 1.28 | 0.10 | 105.63 |
Chromosome . | Number of markers . | Parametric 2-point LOD scores . | . | . | . | Nonparametric . | . | . | |||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | Maximum LOD . | 𝛉 . | Marker (cM*) . | Flanking markers . | Maximum NPL . | P . | cM* . | |||||
1 | 24 | 0.57 | 0.05 | D1S2850 (251.41) | Positive (one side) | 0.19337 | 0.42 | 0 | |||||
4 | 18 | 1.48 | 0.1 | D4S1581 (58.8) | Positive (both sides) | 1.667 | 0.05 | 70 | |||||
4 | 29 | 1.48 | 0.0 | DG4S1397 (65.2) | Positive (both sides) | ||||||||
7 | 10 | 1.33 | 0.1 | D7S661 (155) | Positive (one side) | 0.96 | 0.16 | 155 | |||||
11 | 39 | 1.06 | 0.1 | D11S908 (118.6) | Positive (one side) | 1.28 | 0.10 | 105.63 |
cM location provided by deCODE Genetics.
Discussion
A number of molecular studies on melanoma families with multiple affected individuals of Northern European descent have shown linkage to a couple of genetic loci on chromosome 1 (15-17). However, the gene(s) responsible for these associations are still unknown (18). We found no evidence for linkage to melanoma susceptibility genes in either chromosome 1 or other loci in 61 Italian melanoma-prone families with 2 to 5 cutaneous malignant melanoma that were mutation negative for both CDKN2A and CDK4.
Although possible, it is unlikely that familial aggregation is due to chance alone in this population because the incidence rates of melanoma in Italy are relatively low, with 8.5 cases per 100,000 for men and 10.2 cases per 100,000 for women (19). Aggregation might also be the result of intense sun exposure in susceptible individuals. As shown for the melanocortin-1 receptor (MC1R) gene (20), genetic variants contributing to this aggregation may be common, have low penetrance, and contribute to small increments of melanoma risks, making it difficult to identify them in linkage analyses. A larger number of families and population-based association studies may be required to identify susceptibility loci for melanoma in Mediterranean populations.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Grant support: NIH Intramural Research Program, National Cancer Institute, Division of Cancer Epidemiology and Genetics (M.T. Landi, K.F. Kerstann, P.T. Bradford, A.M. Goldstein, and D. Calista); the Italian Ministry of Health Oncology Ordinary Project 2006 (G. Bianchi-Scarra); and Fondazione Italiana per la Ricerca sul Cancro fellowship (M.C. Scaini).
Note: K.F. Kerstann and P.T. Bradford contributed equally to this work.
Acknowledgments
We thank the families for participating in the study.