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.

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.

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).

Figure 1.

Parametric 2-point LOD scores for the genome-wide scan. All participating individuals were genotyped at 414 loci throughout the genome. Maximum 2-point LOD scores were calculated under dominant and recessive models in an affected-only analysis. Results are plotted as a function of order along the chromosomes. Black line, dominant model; gray line, recessive model.

Figure 1.

Parametric 2-point LOD scores for the genome-wide scan. All participating individuals were genotyped at 414 loci throughout the genome. Maximum 2-point LOD scores were calculated under dominant and recessive models in an affected-only analysis. Results are plotted as a function of order along the chromosomes. Black line, dominant model; gray line, recessive model.

Close modal

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.

Table 1.

Fine-scale mapping: linkage analysis on a 1-cM map for particular regions of interest

ChromosomeNumber of markersParametric 2-point LOD scores
Nonparametric
Maximum LOD𝛉Marker (cM*)Flanking markersMaximum NPLPcM*
24 0.57 0.05 D1S2850 (251.41) Positive (one side) 0.19337 0.42 
18 1.48 0.1 D4S1581 (58.8) Positive (both sides) 1.667 0.05 70 
29 1.48 0.0 DG4S1397 (65.2) Positive (both sides)    
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 
ChromosomeNumber of markersParametric 2-point LOD scores
Nonparametric
Maximum LOD𝛉Marker (cM*)Flanking markersMaximum NPLPcM*
24 0.57 0.05 D1S2850 (251.41) Positive (one side) 0.19337 0.42 
18 1.48 0.1 D4S1581 (58.8) Positive (both sides) 1.667 0.05 70 
29 1.48 0.0 DG4S1397 (65.2) Positive (both sides)    
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.

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.

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.

We thank the families for participating in the study.

1
Newton Bishop J, Harland M, Randerson-Moor J, Bishop T. Management of familial melanoma.
Lancet Oncol
2007
;
8
:
46
–54.
2
Hussussian CJ, Struewing JP, Goldstein AM, et al. Germline p16 mutations in familial melanoma.
Nat Genet
1994
;
8
:
15
–21.
3
Zuo L, Weger J, Yang Q, et al. Germline mutations in p16 binding domain of CDK4 in familial melanoma.
Nat Genet
1996
;
12
:
97
–9.
4
Goldstein AM, Chan M, Harland M, et al. Assessment of high-risk melanoma susceptibility genes and their associations with pancreatic cancer, neural system tumors, and uveal melanoma: a Melanoma Genetics Consortium (Genomel) study.
Cancer Res
2006
;
66
:
9818
–28.
5
Balzi D, Bidoli E, Franceschi S, et al. Estimates of cancer incidence and mortality in Italian regions. Aviano, Italia: Centro di Riferimento Oncologico di Aviano 1997;48–9.
6
Jemal A, Siegel R, Ward E, et al. Cancer statistics, 2008.
CA Cancer J Clin
2008
;
58
:
71
–96.
7
Holland EA, Schmid H, Kefford RF, Mann GJ. CDKN2A (p16) and CDK4 mutation analysis in 131 Australian melanoma probands: effect of family history and multiple primary melanomas.
Genes Chromosomes Cancer
1999
;
25
:
339
–48.
8
Jones WO, Harman CR, Ng AK, Shaw JH. Incidence of malignant melanoma in Auckland, New Zealand: highest rates in the world.
World J Surg
1999
;
23
:
732
–5.
9
Calista D, Goldstein AM, Landi MT. Familial melanoma aggregation in north-eastern Italy.
J Invest Dermatol
2000
;
115
:
764
–5.
10
Landi MT, Goldstein AM, Tsang S, et al. Genetic susceptibility in familial melanoma from northeastern Italy.
J Med Genetics
2004
;
41
:
557
–66.
11
Ott J. Computer-simulation methods in human linkage analysis.
Proc Natl Acad Sci U S A
1989
;
86
:
4175
–8.
12
Weeks DE, Ott J, Lathrop GM. SLINK: a general simulation program for linkage analysis.
Am J Hum Genet
1990
;
47
:
A204
.
13
Verdecchia A, Micheli A, Colonna M, et al. A comparative analysis of cancer prevalence in cancer registry areas of France, Italy and Spain.
Ann Oncol
2002
;
13
:
1128
–39.
14
Hodge SE, Anderson CE, Nejswanger K, et al. The search for heterogeneity in insulin-dependent diabetes mellitus (IDDM): linkage studies, two-locus models, and genetic heterogeneity.
Am J Hum Genet
1983
;
35
:
1139
–55.
15
Bale SJ, Dracopoli NC, Tucker MA, et al. Mapping the gene for hereditary cutaneous malignant melanoma-dysplastic nevus to chromosome 1p.
N Engl J Med
1989
;
320
:
1367
–72.
16
Goldstein AM, Dracopoli NC, Ho EC, et al. Further evidence for a locus for cutaneous malignant melanoma–dysplastic nevus (CMM/DN) on chromosome 1p, and evidence for genetic heterogeneity.
Am J Hum Genet
1993
;
52
:
537
–50.
17
Gillanders E, Juo SH, Holland EA, et al. Localization of a novel melanoma susceptibility locus to 1p22.
Am J Hum Genet
2003
;
73
:
301
–13.
18
de Snoo FA, Hayward NK. Cutaneous melanoma susceptibility and progression genes.
Cancer Lett
2005
;
230
:
153
–86.
19
Boi S, Cristofolini M, Micciolo R, Polla E, Dalla Palma P. Epidemiology of skin tumors: data from the cutaneous cancer registry in Trentino, Italy.
J Cutan Med Surg
2003
;
7
:
300
–5.
20
Landi MT, Kanetsky PA, Tsang S, et al. MC1R, ASIP, and DNA repair in sporadic and familial melanoma in a Mediterranean population.
J Natl Cancer Inst
2005
;
97
:
998
–1007.