Abstract
Recent genome-wide association studies (GWAS) have identified new susceptibility loci for melanoma, but their associations with multiple primary melanoma (MPM) are unclear.
We investigated the associations of 69 SNPs in 39 GWAS-identified loci with odds of MPM relative to single primary melanoma in the international, population-based Genes, Environment, and Melanoma study. Per-minor-allele ORs and 95% confidence intervals (CI) for individuals with MPM “cases” (n = 1,205) relative to single primary melanoma “controls” (n = 2,458) were estimated using multivariable logistic regression, and polygenic risk scores (PRS) were calculated and weighted based on a 2020 GWAS meta-analysis (57 of the 68 independent GWAS SNPs available).
Thirteen SNPs in 11 gene regions (PARP1, CYP1B1/RMDN3, TERT, RAPGEF5, TYRP1, MTAP, CDKN2A/CDKN2B, KLF4, TYR, SOX6, and ASIP) were statistically significantly associated (P < 0.05) with MPM adjusting for age, sex, age-by-sex interaction, and study center. The highest versus lowest PRS quintile was associated with a 2.81-fold higher odds of MPM (95% CI, 2.10–3.78; P = 7.5 × 10−13); this association was attenuated but remained statistically significant after excluding SNPs individually associated with MPM (OR = 1.75, 95% CI, 1.32–2.31).
Inherited genetic variants spanning 11 gene regions were independently associated with MPM. Nonsignificant SNPs were associated with MPM when aggregated into a PRS, indicating that their cumulative effect may influence MPM risk despite lacking individual statistical significance in our study population.
Our findings provide additional evidence that these loci are associated with melanoma risk and estimate the magnitude of their genetic effect on subsequent (multiple) primary melanoma risk.