Melanoma Genetics: Recent Findings Take Us Beyond Well-Traveled Pathways

Aggregation of melanoma tumour site within Western Australian families

BACKGROUND

Evidence is emerging that melanoma has distinct aetiologic pathways and subtypes, characterized by factors like anatomic site of the tumour. To explore genetic influences on anatomic subtypes, we examined the extent to which melanomas in first-degree relatives shared the same body site of occurrence.

METHODS

Population-level linked data was used to identify the study population of over 1.5 million individuals born in Western Australia between 1945 and 2014, and their first-degree relatives. There were 1009 pairs of invasive tumours from 677 family pairs, each categorised by anatomic site. Greater than expected representation of site-concordant pairs would suggest the presence of genetic factors that predispose individuals to site-specific melanoma.

RESULTS

Comparing observed versus expected totals, we observed a modest increase in site concordance for invasive head/neck and truncal tumours (P=0.02). A corresponding analysis including in situ tumours showed a similar concordance (P=0.05). No further evidence of concordance was observed when stratified by sex.

CONCLUSION

In conclusion, modest evidence of aggregation was observed but with inconsistent patterns between sites. Results suggest that further investigation into the familial aggregation of melanoma by tumour site is warranted, with the inclusion of genetic data in order to disentangle the relative contributions of genetic and environmental factors.

The Driverless Triple-Wild-Type (BRAF, RAS, KIT) Cutaneous Melanoma: Whole Genome Sequencing Discoveries

The genetic makeup of the triple-wild-type melanoma (BRAF, NRAS and NF1) has been known for some time, but those studies grouped together rare histopathological versions with common ones, as well as mucosal and even uveal ones. Here we used whole genome sequencing to genetically characterize the triple-wild-type melanoma (TWM), termed here as BRAF, RAS and KIT wild type (the most frequent oncogenic drivers of skin melanoma), using the most common histological forms and excluding rare ones. All these tumors except one were clearly induced by UV based on the mutational signature. The tumor mutational burden was low in TWM, except in the NF1 mutant forms, and a relatively high frequency of elevated LOH scores suggested frequent homologue recombination deficiency, but this was only confirmed by the mutation signature in one case. Furthermore, all these TWMs were microsatellite-stabile. In this driverless setting, we revealed rare oncogenic drivers known from melanoma or other cancer types and identified rare actionable tyrosine kinase mutations in NTRK1, RET and VEGFR1. Mutations of TWM identified genes involved in antitumor immunity (negative and positive predictors of immunotherapy), Ca++ and BMP signaling. The two regressed melanomas of this cohort shared a 17-gene mutation signature, containing genes involved in antitumor immunity and several cell surface receptors. Even with this comprehensive genomic approach, a few cases remained driverless, suggesting that unrecognized drivers are hiding among passenger mutations.

Modern Concepts in Melanocytic Tumors

The advent of molecular pathology has fueled unprecedented advances in the diagnosis and understanding of melanocytic tumors. These advances, however, have also generated concepts that may be difficult to grasp for clinical practitioners, who are not always conversant with the array of genetic techniques employed in the laboratory. These same practitioners, however, are being increasingly called on to provide treatments that are often based on the latest molecular findings for melanocytic tumors. We review the most recent concepts in the pathway classification of melanocytic tumors, including intermediate lesions known as melanocytomas. We examine the genetic and molecular techniques used to study these tumors, look at where they overlap, and discuss their limitations and some of the most difficult-to-interpret results.

Molecular Pathology of Skin Melanoma: Epidemiology, Differential Diagnostics, Prognosis and Therapy Prediction

Similar to other malignancies, TCGA network efforts identified the detailed genomic picture of skin melanoma, laying down the basis of molecular classification. On the other hand, genome-wide association studies discovered the genetic background of the hereditary melanomas and the susceptibility genes. These genetic studies helped to fine-tune the differential diagnostics of malignant melanocytic lesions, using either FISH tests or the myPath gene expression signature. Although the original genomic studies on skin melanoma were mostly based on primary tumors, data started to accumulate on the genetic diversity of the progressing disease. The prognostication of skin melanoma is still based on staging but can be completed with gene expression analysis (DecisionDx). Meanwhile, this genetic knowledge base of skin melanoma did not turn to the expected wide array of target therapies, except the BRAF inhibitors. The major breakthrough of melanoma therapy was the introduction of immune checkpoint inhibitors, which showed outstanding efficacy in skin melanoma, probably due to their high immunogenicity. Unfortunately, beyond BRAF, KIT mutations and tumor mutation burden, no clinically validated predictive markers exist in melanoma, although several promising biomarkers have been described, such as the expression of immune-related genes or mutations in the IFN-signaling pathway. After the initial success of either target or immunotherapies, sooner or later, relapses occur in the majority of patients, due to various induced genetic alterations, the diagnosis of which could be developed to novel predictive genetic markers.

Follow-up of primary melanoma patients with high risk of recurrence: recommendations based on evidence and consensus

In spite of the good prognosis of patients with early-stage melanoma, there is a substantial proportion of them that develop local or distant relapses. With the introduction of targeted and immune therapies for advanced melanoma, including at the adjuvant setting, early detection of recurrent melanoma and/or second primary lesions is crucial to improve clinical outcomes. However, there is a lack of universal guidelines regarding both frequency of surveillance visits and diagnostic imaging and/or laboratory evaluations. In this article, a multidisciplinary expert panel recommends, after careful review of relevant data in the field, a consensus- and experience-based follow-up strategy for melanoma patients, taking into account prognostic factors and biomarkers and the high-risk periods and patterns of recurrence in each (sub) stage of the disease. Apart from the surveillance intensity, healthcare professionals should focus on patients' education to perform regular self-examinations of the skin and palpation of lymph nodes.

T Cell Repertoire Abnormality in Immunodeficiency Patients with DNA Repair and Methylation Defects

Both DNA damage response and methylation play a crucial role in antigen receptor recombination by creating a diverse repertoire in developing lymphocytes, but how their defects relate to T cell repertoire and phenotypic heterogeneity of immunodeficiency remains obscure. We studied the TCR repertoire in patients with the mutation in different genes (ATM, DNMT3B, ZBTB24, RAG1, DCLRE1C, and JAK3) and uncovered distinct characteristics of repertoire diversity. We propose that early aberrancies in thymus T cell development predispose to the heterogeneous phenotypes of the immunodeficiency spectrum. Shorter CDR3 lengths in ATM-deficient patients, resulting from a decreased number of nucleotide insertions during VDJ recombination in the pre-selected TCR repertoire, as well as the increment of CDR3 tyrosine residues, lead to the enrichment of pathology-associated TCRs, which may contribute to the phenotypes of ATM deficiency. Furthermore, patients with DNMT3B and ZBTB24 mutations who exhibit discrepant phenotypes present longer CDR3 lengths and reduced number of known pathology-associated TCRs.

Cell type identification from single-cell transcriptomes in melanoma

BACKGROUND

Single-cell sequencing approaches allow gene expression to be measured at the single-cell level, providing opportunities and challenges to study the aetiology of complex diseases, including cancer.

METHODS

Based on single-cell gene and lncRNA expression levels, we proposed a computational framework for cell type identification that fully considers cell dropout characteristics. First, we defined the dropout features of the cells and identified the dropout clusters. Second, we constructed a differential co-expression network and identified differential modules. Finally, we identified cell types based on the differential modules.

RESULTS

The method was applied to single-cell melanoma data, and eight cell types were identified. Enrichment analysis of the candidate cell marker genes for the two key cell types showed that both key cell types were closely related to the physiological activities of the major histocompatibility complex (MHC); one key cell type was associated with mitosis-related activities, and the other with pathways related to ten diseases.

CONCLUSIONS

Through identification and analysis of key melanoma-related cell types, we explored the molecular mechanism of melanoma, providing insight into melanoma research. Moreover, the candidate cell markers for the two key cell types are potential therapeutic targets for melanoma.

Association of Melanoma-Risk Variants with Primary Melanoma Tumor Prognostic Characteristics and Melanoma-Specific Survival in the GEM Study

Genome-wide association studies (GWAS) and candidate pathway studies have identified low-penetrant genetic variants associated with cutaneous melanoma. We investigated the association of melanoma-risk variants with primary melanoma tumor prognostic characteristics and melanoma-specific survival. The Genes, Environment, and Melanoma Study enrolled 3285 European origin participants with incident invasive primary melanoma. For each of 47 melanoma-risk single nucleotide polymorphisms (SNPs), we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and odds ratios for presence of ulceration, mitoses, and tumor-infiltrating lymphocytes (TILs). We also used Cox proportional hazards regression modeling to estimate the per allele hazard ratios for melanoma-specific survival. Passing the false discovery threshold (p = 0.0026) were associations of IRF4 rs12203592 and CCND1 rs1485993 with log of Breslow thickness, and association of TERT rs2242652 with presence of mitoses. IRF4 rs12203592 also had nominal associations (p < 0.05) with presence of mitoses and melanoma-specific survival, as well as a borderline association (p = 0.07) with ulceration. CCND1 rs1485993 also had a borderline association with presence of mitoses (p = 0.06). MX2 rs45430 had nominal associations with log of Breslow thickness, presence of mitoses, and melanoma-specific survival. Our study indicates that further research investigating the associations of these genetic variants with underlying biologic pathways related to tumor progression is warranted.

Association between a 46-SNP Polygenic Risk Score and melanoma risk in Dutch patients with familial melanoma

BACKGROUND

Familial clustering of melanoma suggests a shared genetic predisposition among family members, but only 10%-40% of familial cases carry a pathogenic variant in a known high-risk melanoma susceptibility gene. We investigated whether a melanoma-specific Polygenic Risk Score (PRS) is associated with melanoma risk in patients with genetically unexplained familial melanoma.

METHODS

Dutch familial melanoma cases (n=418) were genotyped for 46 SNPs previously identified as independently associated with melanoma risk. The 46-SNP PRS was calculated and standardised to 3423 healthy controls (sPRS) and the association between PRS and melanoma risk was modelled using logistic regression. Within the case series, possible differences were further explored by investigating the PRS in relation to (1) the number of primary melanomas in a patient and (2) the extent of familial clustering of melanoma.

RESULTS

The PRS was significantly associated with melanoma risk, with a per-SD OR of 2.12 (95% CI 1.90 to 2.35, p<0.001), corresponding to a 5.70-fold increased risk (95% CI 3.93 to 8.28) when comparing the top 90th to the middle 40-60th PRS percentiles. The mean PRS was significantly higher in cases with multiple primary melanomas than in cases with a single melanoma (sPRS 1.17 vs 0.71, p=0.001). Conversely, cases from high-density melanoma families had a lower (but non-significant) mean PRS than cases from low-density families (sPRS 0.60 vs 0.94, p=0.204).

CONCLUSION

Our work underlines the significance of a PRS in determining melanoma susceptibility and encourages further exploration of the diagnostic value of a PRS in genetically unexplained melanoma families.

FRAMe: Familial Risk Assessment of Melanoma-a risk prediction tool to guide CDKN2A germline mutation testing in Australian familial melanoma

Germline mutations in CDKN2A greatly increase risk of developing cutaneous melanoma. We have constructed a risk prediction model, Familial Risk Assessment of Melanoma (FRAMe), for estimating the likelihood of carrying a heritable CDKN2A mutation among Australian families, where the prevalence of these mutations is low. Using logistic regression, we analysed characteristics of 299 Australian families recruited through the Sydney site of GenoMEL (international melanoma genetics consortium) with at least three cases of cutaneous melanoma (in situ and invasive) among first-degree blood relatives, for predictors of the presence of a pathogenic CDKN2A mutation. The final multivariable prediction model was externally validated in an independent cohort of 61 melanoma kindreds recruited through GenoMEL Queensland. Family variables independently associated with the presence of a CDKN2A mutation in a multivariable model were number of individuals diagnosed with melanoma under 40 years of age, number of individuals diagnosed with more than one primary melanoma, and number of individuals blood related to a melanoma case in the first degree diagnosed with any cancer excluding melanoma and non-melanoma skin cancer. The number of individuals diagnosed with pancreatic cancer was not independently associated with mutation status. The risk prediction model had an area under the receiver operating characteristic curve (AUC) of 0.851 (95% CI 0.793, 0.909) in the training dataset, and 0.745 (95%CI 0.612, 0.877) in the validation dataset. This model is the first to be developed and validated using only Australian data, which is important given the higher rate of melanoma in the population. This model will help to effectively identify families suitable for genetic counselling and testing in areas of high ambient ultraviolet radiation. A user-friendly electronic nomogram is available at www.melanomarisk.org.au .

The prognostic role of BRAF and WNT pathways activation in kinase inhibitors-naïve clinical stage III cutaneous melanoma

The results of local-regional advanced melanoma (stage III) management are still not satisfactory. Particularly, there is no personalized treatment in stage III melanoma patients due to the lack of useful classical pathological markers for prognostication of indolent or aggressive course of the disease. The aim of this study was to explore melanoma genomic landscape by means of the mutational profiling of 50 genes influencing carcinogenesis pathways in the randomly selected 93 kinase inhibitor-naïve (KI-naïve) stage III patients. The genomic alterations were found in 27 out of 50 tested genes and at least one pathogenic variant was detected in 77 out of 93 cases (82.7%). Survival rate was negatively affected by the presence of the somatic mutations in AKT1, ATM, CDH1 and SMARCB1, while the BRAF+ status in KI-naïve stage III population correlated with the longer median overall survival. Genomic alterations in WNT pathway correlated with extranodal adipocyte tissue involvement (P = 0.027) and higher number of metastatic lymph nodes (P = 0.045). In terms of survival, the Cox model confirmed the worse prognosis in patients with mutation in the WNT pathway [hazard ratio (HR) = 2.9, P = 0.017], and better prognosis in cases with mutations in BRAF pathway (HR = 0.5, P = 0.004). WNT/β-catenin pathway alteration was associated with more advanced/aggressive disease. From this perspective, the concept of blocking the activity of the WNT pathway in selected cases appears promising and complementary to the BRAF inhibition therapeutic option for the future.

The 2018 World Health Organization Classification of Cutaneous, Mucosal, and Uveal Melanoma: Detailed Analysis of 9 Distinct Subtypes Defined by Their Evolutionary Pathway

CONTEXT.—

There have been major advances in the understanding of melanoma since the last revision of the World Health Organization (WHO) classification in 2006.

OBJECTIVE.—

To discuss development of the 9 distinct types of melanoma and distinguishing them by their epidemiology, clinical and histologic morphology, and genomic characteristics. Each melanoma subtype is placed at the end of an evolutionary pathway that is rooted in its respective precursor, wherever appropriate and feasible, based on currently known data. Each precursor has a variable risk of progression culminating in its fully evolved, invasive melanoma.

DATA SOURCES.—

This review is based on the "Melanocytic Tumours" section of the 4th edition of the WHO Classification of Skin Tumours, published in 2018.

CONCLUSIONS.—

Melanomas were divided into those etiologically related to sun exposure and those that are not, as determined by their mutational signatures, anatomic site, and epidemiology. Melanomas on the sun-exposed skin were further divided by the histopathologic degree of cumulative solar damage (CSD) of the surrounding skin, into low and high CSD, on the basis of degree of associated solar elastosis. Low-CSD melanomas include superficial spreading melanomas and high-CSD melanomas incorporate lentigo maligna and desmoplastic melanomas. The "nonsolar" category includes acral melanomas, some melanomas in congenital nevi, melanomas in blue nevi, Spitz melanomas, mucosal melanomas, and uveal melanomas. The general term melanocytoma is proposed to encompass "intermediate" tumors that have an increased (though still low) probability of disease progression to melanoma.

Association Between Parkinson's Disease and Melanoma: Putting the Pieces Together

Patients with Parkinson’s disease (PD) generally have reduced risk of developing many types of cancers, except melanoma—a malignant tumor of melanin-producing cells in the skin. For decades, a large number of epidemiological studies have reported that the occurrence of melanoma is higher than expected among subjects with PD, and the occurrence of PD is reciprocally higher than expected among patients with melanoma. More recent epidemiological studies further indicated a bidirectional association, not only in the patients themselves but also in their relatives. This association between PD and melanoma offers a unique opportunity to understand PD. Here, we summarize epidemiological, clinical, and biological evidence in regard to shared risk factors and possible underlying mechanisms for these two seemingly distinct conditions.

Circular RNA ITCH downregulates GLUT1 and suppresses glucose uptake in melanoma to inhibit cancer cell proliferation

BACKGROUND

Circular RNA ITCH has been proved as a tumor suppressor in several types of cancer. Genetic alterations are known to play important roles in melanoma development. However, the function of circular RNAs and the role of ITCH in melanoma are still hardly known.

METHODS

We investigated the expression levels of ITCH and the correlation between ITCH expression and GLUT1 expression, which showed that ITCH was downregulated in melanoma tissues and inversely correlated with GLUT1.

RESULTS

ITCH overexpression resulted in the downregulation of GLUT1 and suppressed glucose uptake in melanoma. GLUT1 overexpression promoted glucose uptake but failed to affect ITCH. ITCH overexpression resulted in increased, while GLUT1 overexpression resulted in decreased rate of melanoma cell proliferation. In addition, GLUT1 overexpression reduced the effects of ITCH overexpression on cell proliferation.

CONCLUSIONS

We concluded that ITCH may downregulate GLUT1 and suppresses glucose uptake in melanoma to inhibit cancer cell proliferation.

Immunomodulatory germline variation associated with the development of multiple primary melanoma (MPM)

Multiple primary melanoma (MPM) has been associated with a higher 10-year mortality risk compared to patients with single primary melanoma (SPM). Given that 3–8% of patients with SPM develop additional primary melanomas, new markers predictive of MPM risk are needed. Based on the evidence that the immune system may regulate melanoma progression, we explored whether germline genetic variants controlling the expression of 41 immunomodulatory genes modulate the risk of MPM compared to patients with SPM or healthy controls. By genotyping these 41 variants in 977 melanoma patients, we found that rs2071304, linked to the expression of SPI1, was strongly associated with MPM risk reduction (OR = 0.60; 95% CI = 0.45–0.81; p = 0.0007) when compared to patients with SPM. Furthermore, we showed that rs6695772, a variant affecting expression of BATF3, is also associated with MPM-specific survival (HR = 3.42; 95% CI = 1.57–7.42; p = 0.0019). These findings provide evidence that the genetic variation in immunomodulatory pathways may contribute to the development of secondary primary melanomas and also associates with MPM survival. The study suggests that inherited host immunity may play an important role in MPM development.

Association between risk factors and detection of cutaneous melanoma in the setting of a population-based skin cancer screening

Early detection is considered to improve the prognosis of cutaneous melanoma. The value of population-based screening for melanoma, however, is still controversial. The aim of this study was to evaluate the predictive power of established risk factors in the setting of a population-based screening and to provide empirical evidence for potential risk stratifications. We reanalyzed data (including age, sex, risk factors, and screening results) of 354 635 participants in the Skin Cancer Research to provide Evidence for Effectiveness of Screening in Northern Germany (SCREEN)project conducted in the German state of Schleswig-Holstein (2003-2004). In multivariable analysis, atypical nevi [odds ratio (OR): 17.4; 95% confidence interval (CI): 14.4-20.1], personal history of melanoma (OR: 5.3; 95% CI: 3.6-7.6), and multiple (≥40) common nevi (OR: 1.3; 95% CI: 1.1-1.6) were associated with an increased risk of melanoma detection. Family history and congenital nevi were not significantly associated with melanoma detection in the SCREEN. The effects of several risk-adapted screening strategies were evaluated. Hypothesizing a screening of individuals aged more than or equal to 35 years, irrespective of risk factors (age approach), the number needed to screen is 559 (95% CI: 514-612), whereas a screening of adults (aged ≥20) with at least one risk factor (risk approach) leads to a number needed to screen of 178 (95% CI: 163-196). Converted into one screen-detected melanoma, the number of missed melanomas is 0.15 (95% CI: 0.12-0.18) with the age approach and 0.22 (95% CI: 0.19-0.26) with the risk approach. The results indicate that focusing on individuals at high risk for melanoma may improve the cost-effectiveness and the benefit-to-harm balance of melanoma screening programs.

New susceptibility loci for cutaneous melanoma risk and progression revealed using a porcine model

Despite major advances, it is estimated that a large part of melanoma predisposing genes remains to be discovered. Animal models of spontaneous diseases are valuable tools and experimental crosses can be used to identify and fine-map new susceptibility loci associated with melanoma. We performed a Genome-Wide Association Study (GWAS) of melanoma occurrence and progression (clinical ulceration and presence of metastasis) in a porcine model of spontaneous melanoma, the MeLiM pig. Five loci on chromosomes 2, 5, 7, 8 and 16 showed genome-wide significant associations (p < 5 × 10^–6) with either one of these phenotypes. Suggestive associations (p < 5 × 10^–5) were also found at 16 additional loci. Moreover, comparison of the porcine results to those reported by human melanoma GWAS indicated shared association signals notably at CDKAL1 and TERT loci but also nearby CCND1, FTO, PLA2G6 and TMEM38B-RAD23B loci. Extensive search of the literature revealed a potential key role of genes at the identified porcine loci in tumor invasion (DST, PLEKHA5, CBY1, LIMK2 and ETV5) and immune response modulation (ETV5, HERC3 and DICER1) of the progression phenotypes. These biological processes are consistent with the clinico-pathological features of MeLiM tumors and can open new routes for future melanoma research in humans.

The melanoma genomics managing your risk study: A protocol for a randomized controlled trial evaluating the impact of personal genomic risk information on skin cancer prevention behaviors

BACKGROUND

Reducing ultraviolet radiation (UV) exposure and improving early detection may reduce melanoma incidence, mortality and health system costs. This study aims to evaluate the efficacy and cost-effectiveness of providing information on personal genomic risk of melanoma in reducing UV exposure at 12 months, according to low and high traditional risk.

METHODS

In this randomized controlled trial, participants (target sample = 892) will be recruited from the general population, and randomized (1:1 ratio, intervention versus control). Intervention arm participants provide a saliva sample, receive personalized melanoma genomic risk information, a genetic counselor phone call, and an educational booklet on melanoma prevention. Control arm participants receive only the educational booklet. Eligible participants are aged 18-69 years, have European ancestry and no personal history of melanoma. All participants will complete a questionnaire and wear a UV dosimeter to objectively measure their sun exposure at baseline, 1- and 12-month time-points, except 1-month UV dosimetry will be limited to ~250 participants. The primary outcome is total daily Standard Erythemal Doses at 12 months. Secondary outcomes include objectively measured UV exposure for specific time periods (e.g. midday hours), self-reported sun protection and skin-examination behaviors, psycho-social outcomes, and ethical considerations surrounding offering genomic testing at a population level. A within-trial and modelled economic evaluation will be undertaken from an Australian health system perspective to assess the intervention costs and outcomes.

DISCUSSION

This trial will inform the clinical and personal utility of introducing genomic testing into the health system for melanoma prevention and early detection at a population-level.

TRIAL REGISTRATION

Australian New Zealand Clinical Trials Registry ACTRN12617000691347.

Polyunsaturated fatty acids and risk of melanoma: A Mendelian randomisation analysis

Melanoma is the deadliest form of skin cancer, mainly affecting populations of European ancestry. Some observational studies suggest that particular diets reduce melanoma risk, putatively through an increase in polyunsaturated fatty acid (PUFA) consumption. However, interpretation of these observational findings is difficult due to residual confounding or reverse causality. To date, a randomized controlled trial has not been carried out to examine the relationship between PUFAs and melanoma. Hence, we performed a Mendelian randomisation (MR) study to evaluate the link between PUFAs and melanoma. To perform MR, we used summary results from the largest risk genome-wide association study (GWAS) meta-analysis of melanoma, consisting of 12,874 cases and 23,203 controls. As instrumental variables we selected SNPs associated with PUFA levels from a GWAS meta-analysis of PUFA levels, from the CHARGE consortium. We used the inverse variance weighted method to estimate a causal odds ratio. To aid interpretation, we established a benchmark "large" predicted change in PUFAs in which, for example, an increase in docosahexaenoic acid (DPA) of 0.17 units (equal to 1 standard deviation) moves a person from the 17^th percentile to the median. Raising PUFA levels by a large amount (increasing DPA by 0.17 units) only negligibly changed melanoma risk: odds ratio [OR] = 1.03 (95% confidence interval [CI] = 0.96-1.10). Other PUFAs yielded similar results as DPA. Our MR analysis suggests that the effect of PUFA levels on melanoma risk is either zero or very small.

Integrated case-control and somatic-germline interaction analyses of melanoma susceptibility genes

While a number of genes have been implicated in melanoma susceptibility, the role of protein-coding variation in melanoma development and progression remains underexplored. To better characterize the role of germline coding variation in melanoma, we conducted a whole-exome case-control and somatic-germline interaction study involving 322 skin cutaneous melanoma cases from The Cancer Genome Atlas and 3607 controls of European ancestry. We controlled for cross-platform technological stratification using XPAT and conducted gene-based association tests using VAAST 2. Four established melanoma susceptibility genes achieved nominal statistical significance, MC1R (p = .0014), MITF (p = .0165) BRCA2 (p = .0206), and MTAP (p = .0393). We also observed a suggestive association for FANCA (p = .002), a gene previously implicated in melanoma survival. The association signal for BRCA2 was driven primarily by likely gene disrupting (LGD) variants, with an Odds Ratio (OR) of 5.62 (95% Confidence Interval (CI) 1.03-30.1). In contrast, the association signals for MC1R and MITF were driven primarily by predicted pathogenic missense variants, with estimated ORs of 1.4 to 3.0 for MC1R and 4.1 for MITF. MTAP exhibited an excess of both LGD and predicted damaging missense variants among cases, with ORs of 5.62 and 3.72, respectively, although neither category was significant. For individuals with known or predicted damaging variants, age of disease onset was significantly lower for two of the four genes, MC1R (p = .005) and MTAP (p = .035). In an analysis of germline carrier status and overlapping copy number alterations, we observed no evidence to support a two-hit model of carcinogenesis in any of the four genes. Although MC1R carriers were represented proportionally among the four molecular tumor subtypes, these individuals accounted for 69% of ultraviolet (UV) radiation mutational signatures among triple-wild type tumors (p = .040), highlighting the increased sensitivity to UV exposure among individuals with loss-of-function variants in MC1R.

Rare germline variants in known melanoma susceptibility genes in familial melanoma

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in <40% of melanoma-prone families, suggesting the existence of additional high-risk genes or perhaps a polygenic mechanism involving multiple genetic modifiers. The goal of this study was to systematically characterize rare germline variants in 42 established melanoma genes among 144 CMM patients in 76 American CMM families without known mutations using data from whole-exome sequencing. We identified 68 rare (<0.1% in public and in-house control datasets) nonsynonymous variants in 25 genes. We technically validated all loss-of-function, inframe insertion/deletion, and missense variants predicted as deleterious, and followed them up in 1, 559 population-based CMM cases and 1, 633 controls. Several of these variants showed disease co-segregation within families. Of particular interest, a stopgain variant in TYR was present in five of six CMM cases/obligate gene carriers in one family and a single population-based CMM case. A start gain variant in the 5'UTR region of PLA2G6 and a missense variant in ATM were each seen in all three affected people in a single family, respectively. Results from rare variant burden tests showed that familial and population-based CMM patients tended to have higher frequencies of rare germline variants in albinism genes such as TYR, TYRP1, and OCA2 (P < 0.05). Our results suggest that rare nonsynonymous variants in low- or intermediate-risk CMM genes may influence familial CMM predisposition, warranting further investigation of both common and rare variants in genes affecting functionally important pathways (such as melanogenesis) in melanoma risk assessment.

Genetic variation in IRF4 expression modulates growth characteristics, tyrosinase expression and interferon-gamma response in melanocytic cells

A SNP within intron4 of the interferon regulatory factor4 (IRF4) gene, rs12203592*C/T, has been independently associated with pigmentation and age-specific effects on naevus count in European-derived populations. We have characterized the cis-regulatory activity of this intronic region and using human foreskin-derived melanoblast strains, we have explored the correlation between IRF4 rs12203592 homozygous C/C and T/T genotypes with TYR enzyme activity, supporting its association with pigmentation traits. Further, higher IRF4 protein levels directed by the rs12203592*C allele were associated with increased basal proliferation but decreased cell viability following UVR, an etiological factor in melanoma development. Since UVR, and accompanying IFNγ-mediated inflammatory response, is associated with melanomagenesis, we evaluated its effects in the context of IRF4 status. Manipulation of IRF4 levels followed by IFNγ treatment revealed a subset of chemokines and immuno-evasive molecules that are sensitive to IRF4 expression level and genotype including CTLA4 and PD-L1.

The MeLiM Minipig: An Original Spontaneous Model to Explore Cutaneous Melanoma Genetic Basis

Melanoma is the deadliest skin cancer and is a major public health concern with a growing incidence worldwide. As for other complex diseases, animal models are needed in order to better understand the mechanisms leading to pathology, identify potential biomarkers to be used in the clinics, and eventually molecular targets for therapeutic solutions. Cutaneous melanoma, arising from skin melanocytes, is mainly caused by environmental factors such as UV radiation; however a significant genetic component participates in the etiology of the disease. The pig is a recognized model for spontaneous development of melanoma with features similar to the human ones, followed by a complete regression and a vitiligo-like depigmentation. Three different pig models (MeLiM, Sinclair, and MMS-Troll) have been maintained through the last decades, and different genetic studies have evidenced a complex inheritance of the disease. As in humans, pigmentation seems to play a prominent role, notably through MC1R and MITF signaling. Conversely, cell cycle genes as CDKN2A and CDK4 have been excluded as predisposing for melanoma in MeLiM. So far, only sparse studies have focused on somatic changes occurring during oncogenesis, and have revealed major cytological changes and a potential dysfunction of the telomere maintenance system. Finally, the spontaneous tumor progression and regression occurring in these models could shed light on the interplay between endogenous retroviruses, melanomagenesis, and adaptive immune response.

Integrative genetic analysis suggests that skin color modifies the genetic architecture of melanoma

Melanoma is the deadliest form of skin cancer and presents a significant health care burden in many countries. In addition to ultraviolet radiation in sunlight, the main causal factor for melanoma, genetic factors also play an important role in melanoma susceptibility. Although genome-wide association studies have identified many single nucleotide polymorphisms associated with melanoma, little is known about the proportion of disease risk attributable to these loci and their distribution throughout the genome. Here, we investigated the genetic architecture of melanoma in 1,888 cases and 990 controls of European non-Hispanic ancestry. We estimated the overall narrow-sense heritability of melanoma to be 0.18 (P < 0.03), indicating that genetics contributes significantly to the risk of sporadically-occurring melanoma. We then demonstrated that only a small proportion of this risk is attributable to known risk variants, suggesting that much remains unknown of the role of genetics in melanoma. To investigate further the genetic architecture of melanoma, we partitioned the heritability by chromosome, minor allele frequency, and functional annotations. We showed that common genetic variation contributes significantly to melanoma risk, with a risk model defined by a handful of genomic regions rather than many risk loci distributed throughout the genome. We also demonstrated that variants affecting gene expression in skin account for a significant proportion of the heritability, and are enriched among melanoma risk loci. Finally, by incorporating skin color into our analyses, we observed both a shift in significance for melanoma-associated loci and an enrichment of expression quantitative trait loci among melanoma susceptibility variants. These findings suggest that skin color may be an important modifier of melanoma risk. We speculate that incorporating skin color and other non-genetic factors into genetic studies may allow for an improved understanding of melanoma susceptibility and guide future investigations to identify melanoma risk genes.

Role Played by Signalling Pathways in Overcoming BRAF Inhibitor Resistance in Melanoma

The discovery of the BRAF^V600E mutation led to the development of vemurafenib (PLX4032), a selective BRAF inhibitor specific to the kinase, for the treatment of metastatic melanomas. However, initial success of the drug was dampened by the development of acquired resistance. Melanoma was shown to relapse in patients following treatment with vemurafenib which eventually led to patients' deaths. It has been proposed that mechanisms of resistance can be due to (1) reactivation of the mitogen-activated protein kinase (MAPK) signalling pathway via secondary mutations, amplification or activation of target kinase(s), (2) the bypass of oncogenic pathway via activation of alternative signalling pathways, (3) other uncharacterized mechanisms. Studies showed that receptor tyrosine kinases (RTK) such as PDGFRβ, IGF1R, EGFR and c-Met were overexpressed in melanoma cells. Along with increased secretion of growth factors such as HGF and TGF-α, this will trigger intracellular signalling cascades. This review discusses the role MAPK and Phosphatidylinositol-3-kinase-protein kinase B-mammalian target of rapamycin (PI3K-AKT-mTOR) pathways play in the mechanism of resistance of melanomas.

DNA methylation/hydroxymethylation in melanoma

Melanoma is a malignant tumor of melanocytes and is considered to be the most aggressive cancer among all skin diseases. The pathogenesis of melanoma has not been well documented, which may restrict the research and development of biomarkers and therapies. To date, several genetic and epigenetic factors have been identified as contributing to the development and progression of melanoma. Besides the findings on genetic susceptibilities, the recent progress in epigenetic studies has revealed that loss of the DNA hydroxymethylation mark, 5-hydroxymethylcytosine (5-hmC), along with high levels of DNA methylation at promoter regions of several tumor suppressor genes in melanoma, may serve as biomarkers for melanoma. Moreover, 5-Aza-2′-deoxycytidine, an epigenetic modifier causing DNA demethylation, and ten-eleven translocation family dioxygenase (TET), which catalyzes the generation of 5-hmC, demonstrate therapeutic potential in melanoma treatment. In this review, we will summarize the latest progress in research on DNA methylation/hydroxymethylation in melanoma, and we will discuss and provide insight for epigenetic biomarkers and therapies for melanoma. Particularly, we will discuss the role of DNA hydroxymethylation in melanoma infiltrating immune cells, which may also serve as a potential target for melanoma treatment.

The Pathophysiological Impact of HLA Class Ia and HLA-G Expression and Regulatory T Cells in Malignant Melanoma: A Review

Malignant melanoma, a very common type of cancer, is a rapidly growing cancer of the skin with an increase in incidence among the Caucasian population. The disease is seen through all age groups and is very common in the younger age groups. Several studies have examined the risk factors and pathophysiological mechanisms of malignant melanoma, which have enlightened our understanding of the development of the disease, but we have still to fully understand the complex immunological interactions. The examination of the interaction between the human leucocyte antigen (HLA) system and prognostic outcome has shown interesting results, and a correlation between the down- or upregulation of these antigens and prognosis has been seen through many different types of cancer. In malignant melanoma, HLA class Ia has been seen to influence the effects of pharmaceutical drug treatment as well as the overall prognosis, and the HLA class Ib and regulatory T cells have been correlated with tumor progression. Although there is still no standardized immunological treatment worldwide, the interaction between the human leucocyte antigen (HLA) system and tumor progression seems to be a promising focus in the way of optimizing the treatment of malignant melanoma.

Current Data on Risk Factor Estimates Does Not Explain the Difference in Rates of Melanoma between Hispanics and Non-Hispanic Whites

United States Hispanics have seven times lower melanoma incidence rates than non-Hispanic whites (NHW). It is unclear whether this difference can be explained solely by phenotypic risk factors, like darker skin, or whether modifiable risk factors, like sun exposure, also play a role. The purpose of this paper is to summarize what is currently known about melanoma risk factors among Hispanics and NHWs, and whether or not those differences could explain the difference in melanoma incidence. Through literature review, relative risks and prevalence of melanoma risk factors in Hispanics and NHWs were identified and used to calculate the expected rate in Hispanics and rate ratio compared to NHWs. We found that melanoma risk factors either have similar frequency in Hispanics and NHWs (e.g., many large nevi) or are less frequent in Hispanics but do not explain a high proportion of disease variation (e.g., red hair). Considering current knowledge of risk factor prevalence, we found that melanoma incidence rates in the two groups should actually be similar. Sun exposure behavior among Hispanics may contribute to the explanation for the 7-fold difference in melanoma rates. Currently, limited data exist on sun exposure behavior among Hispanics, but possibilities for improving primary prevention by further studying these practices are substantial.

Shorter telomere length in Europeans than in Africans due to polygenetic adaptation

Leukocyte telomere length (LTL), which reflects telomere length in other somatic tissues, is a complex genetic trait. Eleven SNPs have been shown in genome-wide association studies to be associated with LTL at a genome-wide level of significance within cohorts of European ancestry. It has been observed that LTL is longer in African Americans than in Europeans. The underlying reason for this difference is unknown. Here we show that LTL is significantly longer in sub-Saharan Africans than in both Europeans and African Americans. Based on the 11 LTL-associated alleles and genetic data in phase 3 of the 1000 Genomes Project, we show that the shifts in allele frequency within Europe and between Europe and Africa do not fit the pattern expected by neutral genetic drift. Our findings suggest that differences in LTL within Europeans and between Europeans and Africans is influenced by polygenic adaptation and that differences in LTL between Europeans and Africans might explain, in part, ethnic differences in risks for human diseases that have been linked to LTL.

Association of Interferon Regulatory Factor-4 Polymorphism rs12203592 With Divergent Melanoma Pathways

BACKGROUND

Solar elastosis and neval remnants are histologic markers characteristic of divergent melanoma pathways linked to differences in age at onset, host phenotype, and sun exposure. However, the association between these pathway markers and newly identified low-penetrance melanoma susceptibility loci remains unknown.

METHODS

In the Genes, Environment and Melanoma (GEM) Study, 2103 Caucasian participants had first primary melanomas that underwent centralized pathology review. For 47 single-nucleotide polymorphisms (SNPs) previously identified as low-penetrant melanoma risk variants, we used multinomial logistic regression to compare melanoma with solar elastosis and melanoma with neval remnants simultaneously to melanoma with neither of these markers, excluding melanomas with both markers. All statistical tests were two-sided.

RESULTS

IRF4 rs12203592 was the only SNP to pass the false discovery threshold in baseline models adjusted for age, sex, and study center. rs12203592*T was associated positively with melanoma with solar elastosis (odds ratio [OR] = 1.47, 95% confidence interval [CI] = 1.18 to 1.82) and inversely with melanoma with neval remnants (OR = 0.65, 95% CI = 0.48 to 0.87) compared with melanoma with neither marker (P global = 3.78 x 10(-08)). Adjusting for phenotypic characteristics and total sun exposure hours did not materially affect rs12203592's associations. Distinct early- and late-onset age distributions were observed in patients with IRF4 rs12203592 [CC] and [TT] genotypes, respectively.

CONCLUSIONS

Our findings suggest a role of IRF4 rs12203592 in pathway-specific risk for melanoma development. We hypothesize that IRF4 rs12203592 could underlie in part the bimodal age distribution reported for melanoma and linked to the divergent pathways.

Melanoma Epidemiology and Prevention

The epidemiology of melanoma is complex, and individual risk depends on sun exposure, host factors, and genetic factors, and in their interactions as well. Sun exposure can be classified as intermittent, chronic, or cumulative (overall) exposure, and each appears to have a different effect on type of melanoma. Other environmental factors, such as chemical exposures-either through occupation, atmosphere, or food-may increase risk for melanoma, and this area warrants further study. Host factors that are well known to be important are the numbers and types of nevi and the skin phenotype. Genetic factors are classified as high-penetrant genes, moderate-risk genes, or low-risk genetic polymorphisms. Subtypes of tumors, such as BRAF-mutated tumors, have different risk factors as well as different therapies. Prevention of melanoma has been attempted using various strategies in specific subpopulations, but to date optimal interventions to reduce incidence have not emerged.

High incidence of primary melanomas in an MC1R RHC homozygote/CDKN2A mutant genotype patient

Melanoma incidence in Australia remains the highest in the world; hence understanding its causation is paramount for future therapeutic developments. Multiple primary melanomas are also common occurrences among the Australian population with identified risk factors such as personal and family history of melanoma, fair skin type, dysplastic naevus syndrome and history of significant ultraviolet exposure. The roles of both environmental and genetic factors have been elucidated in melanoma development, but the synergy of interactions between the two remains complex given the heterogeneous nature of the disease. We present a rare case of a 57-year-old female with 20 cutaneous melanomas and review the role of genetic and environmental factors in development of her multiple primary melanomas.

Peripheral T-cell and NK cell lymphoproliferative disorders: cell of origin, clinical and pathological implications

T-cell lymphoproliferative disorders are a heterogeneous group of neoplasms with distinct clinical-biological properties. The normal cellular counterpart of these processes has been postulated based on functional and immunophenotypic analyses. However, T lymphocytes have been proven to be remarkably capable of modulating their properties, adapting their function in relationship with multiple stimuli and to the microenvironment. This impressive plasticity is determined by the equilibrium among a pool of transcription factors and by DNA chromatin regulators. It is now proven that the acquisition of specific genomic defects leads to the enforcement/activation of distinct pathways, which ultimately alter the preferential activation of defined regulators, forcing the neoplastic cells to acquire features and phenotypes distant from their original fate. Thus, dissecting the landscape of the genetic defects and their functional consequences in T-cell neoplasms is critical not only to pinpoint the origin of these tumors but also to define innovative mechanisms to re-adjust an unbalanced state to which the tumor cells have become addicted and make them vulnerable to therapies and targetable by the immune system. In our review, we briefly describe the pathological and clinical aspects of the T-cell lymphoma subtypes as well as NK-cell lymphomas and then focus on the current understanding of their pathogenesis and the implications on diagnosis and treatment.

Interactions between ultraviolet light and MC1R and OCA2 variants are determinants of childhood nevus and freckle phenotypes

BACKGROUND

Melanocytic nevi (moles) and freckles are well known biomarkers of melanoma risk, and they are influenced by similar UV light exposures and genetic susceptibilities to those that increase melanoma risk. Nevertheless, the selective interactions between UV exposures and nevus and freckling genes remain largely undescribed.

METHODS

We conducted a longitudinal study from ages 6 through 10 years in 477 Colorado children who had annual information collected for sun exposure, sun protection behaviors, and full body skin exams. MC1R and HERC2/OCA2 rs12913832 were genotyped and linear mixed models were used to identify main and interaction effects.

RESULTS

All measures of sun exposure (chronic, sunburns, and waterside vacations) contributed to total nevus counts, and cumulative chronic exposure acted as the major driver of nevus development. Waterside vacations strongly increased total nevus counts in children with rs12913832 blue eye color alleles and facial freckling scores in those with MC1R red hair color variants. Sunburns increased the numbers of larger nevi (≥2 mm) in subjects with certain MC1R and rs12913832 genotypes.

CONCLUSIONS

Complex interactions between different UV exposure profiles and genotype combinations determine nevus numbers and size, and the degree of facial freckling.

IMPACT

Our findings emphasize the importance of implementing sun-protective behavior in childhood regardless of genetic make-up, although children with particular genetic variants may benefit from specifically targeted preventive measures to counteract their inherent risk of melanoma. Moreover, we demonstrate, for the first time, that longitudinal studies are a highly powered tool to uncover new gene-environment interactions that increase cancer risk.

Epigenetic marks responsible for cadmium-induced melanoma cell overgrowth

Cadmium (Cd) is a human carcinogen that likely acts via epigenetic mechanisms. However, the precise role of Cd in melanoma remains to be defined. The goals of this study are to: (i) examine the effect of Cd on the proliferation rate of cutaneous and uveal melanoma cells; (ii) identify the genes affected by Cd exposure; (iii) understand whether epigenetic changes are involved in the response to Cd. The cell growth capacity increased at 48 h after Cd treatment at doses ranging from 0.5 to 10 μM. The research on the key genes regulating proliferation has shown that aberrant methylation is responsible for silencing of p16(INK4A) and caspase 8 in uveal and cutaneous melanoma cells, respectively. The methylation and expression patterns of p14(ARF), death receptors 4/5, and E-cadherin remained unmodified after Cd treatment in all the cell lines analyzed. Ectopic expression of p16(INK4A) abolished the overgrowth of uveal melanoma cells in response to Cd and the overexpression of caspase 8 drastically increased the apoptotic rate of Cd-treated cutaneous melanoma cells. In conclusion, our data suggest that hypermethylation of p16(INK4A) and caspase 8 represents the most common event linked to Cd-induced stimulation of cell growth and inhibition of cell death pathway in melanoma.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Genetic epidemiology of malignant melanoma susceptibility

Germline CDKN2A mutations were the first to be associated with familial melanoma. MC1R polymorphisms are associated, in conformity with epidemiological observations, with fair skin phenotype and a moderately increased risk for melanoma. The wider implementation of genome-wide association studies along with improved whole exome sequencing techniques made possible the identification of novel high-penetrant mutations (TERT, MITF, POT1, BAP1) beyond the established pathways of pigmentation and nevus count suggesting an additional role for pathways involved in cell cycle control and DNA repair. A multitude of common polymorphisms in the general population have been associated through candidate gene studies with a low risk for melanoma, supporting the hypothesis of a complex disease.

Epidemiological and genetic factors underlying melanoma development in Italy

Among human cancers, melanoma remains one of the malignancies with an ever-growing incidence in white populations. Recent advances in biological and immunological therapeutic approaches as well as increased efforts for secondary prevention are contributing to improve the survival rates. It is likely that a significant fall in mortality rates for melanoma will be achieved by further increase of the early detection through a more accurate selection of the higher-risk individuals (i.e., carriers of predisposing genetic alterations). A similar scenario occurs in Italy. In the present review, we have considered data on incidence, survival and mortality rates of melanoma in Italian population, including evaluation of the main risk factors and genetic mutations underlying disease susceptibility.

DNA Methylation Levels of Melanoma Risk Genes Are Associated with Clinical Characteristics of Melanoma Patients

In melanoma development, oncogenic process is mediated by genetic and epigenetic mutations, and few studies have so far explored the role of DNA methylation either as predisposition factor or biomarker. We tested patient samples for germline CDKN2A methylation status and found no evidence of inactivation by promoter hypermethylation. We have also investigated the association of clinical characteristics of samples with the DNA methylation pattern of twelve genes relevant for melanomagenesis. Five genes (BAP1, MGMT, MITF, PALB2, and POT1) presented statistical association between blood DNA methylation levels and either CDKN2A-mutation status, number of lesions, or Breslow thickness. In tumors, five genes (KIT, MGMT, MITF, TERT, and TNF) exhibited methylation levels significantly different between tumor groups including acral compared to nonacral melanomas and matched primary lesions and metastases. Our data pinpoint that the methylation level of eight melanoma-associated genes could potentially represent markers for this disease both in peripheral blood and in tumor samples.

Inherited genetic variants associated with occurrence of multiple primary melanoma

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.

Updated field synopsis and systematic meta-analyses of genetic association studies in cutaneous melanoma: the MelGene database

We updated a field synopsis of genetic associations of cutaneous melanoma (CM) by systematically retrieving and combining data from all studies in the field published as of August 31, 2013. Data were available from 197 studies, which included 83,343 CM cases and 187,809 controls and reported on 1,126 polymorphisms in 289 different genes. Random-effects meta-analyses of 81 eligible polymorphisms evaluated in >4 data sets confirmed 20 single-nucleotide polymorphisms across 10 loci (TYR, AFG3L1P, CDK10, MYH7B, SLC45A2, MTAP, ATM, CLPTM1L, FTO, and CASP8) that have previously been published with genome-wide significant evidence for association (P<5 × 10(-8)) with CM risk, with certain variants possibly functioning as proxies of already tagged genes. Four other loci (MITF, CCND1, MX2, and PLA2G6) were also significantly associated with 5 × 10(-8) Collapse

Key Words Collapse

MESH Headings

• Chromosome Mapping
• Databases, Genetic
• Gene Expression Regulation, Neoplastic
• Genetic Predisposition to Disease
• Genetic Variation
• Genome-Wide Association Study
• Germ-Line Mutation
• Humans
• Linkage Disequilibrium
• Melanoma/genetics
• Odds Ratio
• Polymorphism, Single Nucleotide
• Skin Neoplasms/genetics
• Melanoma, Cutaneous Malignant

Collapse

Grants Collapse

Affiliation(s)

Kyriaki Antonopoulou Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece
Irene Stefanaki Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece
Christina M Lill Neuropsychiatric Genetics Group, Department of Vertebrate Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany; Department of Neurology, Focus Program Translational Neuroscience, University Medical Center of the Johannes Gutenberg University Mainz, Mainz, Germany
Foteini Chatzinasiou Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece
Katerina P Kypreou Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece
Fani Karagianni Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece
Emmanouil Athanasiadis Center of Systems Biology, Biomedical Research Foundation, Academy of Athens, Athens, Greece
George M Spyrou Center of Systems Biology, Biomedical Research Foundation, Academy of Athens, Athens, Greece
John P A Ioannidis Department of Medicine, Stanford Prevention Research Center, Stanford University School of Medicine, Stanford, California, USA
Lars Bertram Neuropsychiatric Genetics Group, Department of Vertebrate Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany; Department of Medicine, School of Public Health, Imperial College London, London, UK
Evangelos Evangelou Department of Hygiene and Epidemiology, Clinical and Molecular Epidemiology Unit, School of Medicine, University of Ioannina, Ioannina, Greece; Department of Epidemiology and Biostatistics, Imperial College London, St Mary's Campus, London, UK
Alexander J Stratigos Department of Dermatology, University of Athens School of Medicine, Andreas Sygros Hospital, Athens, Greece.

Collapse

DCAF4, a novel gene associated with leucocyte telomere length

BACKGROUND

Leucocyte telomere length (LTL), which is fashioned by multiple genes, has been linked to a host of human diseases, including sporadic melanoma. A number of genes associated with LTL have already been identified through genome-wide association studies. The main aim of this study was to establish whether DCAF4 (DDB1 and CUL4-associated factor 4) is associated with LTL. In addition, using ingenuity pathway analysis (IPA), we examined whether LTL-associated genes in the general population might partially explain the inherently longer LTL in patients with sporadic melanoma, the risk for which is increased with ultraviolet radiation (UVR).

RESULTS

Genome-wide association (GWA) meta-analysis and de novo genotyping of 20 022 individuals revealed a novel association (p=6.4×10(-10)) between LTL and rs2535913, which lies within DCAF4. Notably, eQTL analysis showed that rs2535913 is associated with decline in DCAF4 expressions in both lymphoblastoid cells and sun-exposed skin (p=4.1×10(-3) and 2×10(-3), respectively). Moreover, IPA revealed that LTL-associated genes, derived from GWA meta-analysis (N=9190), are over-represented among genes engaged in melanoma pathways. Meeting increasingly stringent p value thresholds (p<0.05, <0.01, <0.005, <0.001) in the LTL-GWA meta-analysis, these genes were jointly over-represented for melanoma at p values ranging from 1.97×10(-169) to 3.42×10(-24).

CONCLUSIONS

We uncovered a new locus associated with LTL in the general population. We also provided preliminary findings that suggest a link of LTL through genetic mechanisms with UVR and melanoma in the general population.

Height and bone mineral density are associated with naevus count supporting the importance of growth in melanoma susceptibility

Naevus count is the strongest risk factor for melanoma. Body Mass Index (BMI) has been linked to melanoma risk. In this study, we investigate the link between naevus count and height, weight and bone mineral density (BMD) in the TwinsUK cohort (N = 2119). In addition we adjusted for leucocyte telomere length (LTL) as LTL is linked to both BMD and naevus count. Naevus count was positively associated with height (p = 0.001) but not with weight (p = 0.187) despite adjusting for age and twin relatedness. This suggests that the previously reported melanoma association with BMI may be explained by height alone. Further adjustment for LTL did not affect the significance of the association between height and naevus count so LTL does not fully explain these results. BMD was associated with naevus count at the spine (coeff 18.9, p = 0.01), hip (coeff = 18.9, p = 0.03) and forearm (coeff = 32.7, p = 0.06) despite adjusting for age, twin relatedness, weight, height and LTL. This large study in healthy individuals shows that growth via height, probably in early life, and bone mass are risk factors for melanoma via increased naevus count. The link between these two phenotypes may possibly be explained by telomere biology, differentiation genes from the neural crests but also other yet unknown factors which may influence both bones and melanocytes biology.

Is Mc1r an important regulator of non-pigmentary responses to UV radiation?

MC1R is recognized for its role in the regulation of melanin pigmentation. In addition, many investigators believe that it also plays a crucial role in immunomodulation (immunosuppression) and in melanogenesis-independent protective responses against ultraviolet radiation (UVR). Surprisingly, Wolnicka-Glubisz et al. have shown that loss of function in the MC1R has no effect on inflammatory responses and immunosuppression induced by UVR in C57BL/6 mice as well as on the degree of UVA-induced DNA damage in the epidermis and dermis. These findings, by challenging the existing dogmas on the precise role of MC1R in non-pigmentary responses to the UVR, mandate further research to either validate the presented data or to define to which degree these phenomena are restricted to the C57BL/6 mouse model or are applicable to other species including humans. The alternative target for immunomodulation is represented by MC3R. However, cutaneous expression of MC3R remains to be demonstrated.

POT1 loss-of-function variants predispose to familial melanoma

Deleterious germline variants in CDKN2A account for around 40% of familial melanoma cases, and rare variants in CDK4, BRCA2, BAP1 and the promoter of TERT have also been linked to the disease. Here we set out to identify new high-penetrance susceptibility genes by sequencing 184 melanoma cases from 105 pedigrees recruited in the UK, The Netherlands and Australia that were negative for variants in known predisposition genes. We identified families where melanoma cosegregates with loss-of-function variants in the protection of telomeres 1 gene (POT1), with a proportion of family members presenting with an early age of onset and multiple primary tumors. We show that these variants either affect POT1 mRNA splicing or alter key residues in the highly conserved oligonucleotide/oligosaccharide-binding (OB) domains of POT1, disrupting protein-telomere binding and leading to increased telomere length. These findings suggest that POT1 variants predispose to melanoma formation via a direct effect on telomeres.

Epistatic and combinatorial effects of pigmentary gene mutations in the domestic pigeon

Understanding the molecular basis of phenotypic diversity is a critical challenge in biology, yet we know little about the mechanistic effects of different mutations and epistatic relationships among loci that contribute to complex traits. Pigmentation genetics offers a powerful model for identifying mutations underlying diversity and for determining how additional complexity emerges from interactions among loci. Centuries of artificial selection in domestic rock pigeons (Columba livia) have cultivated tremendous variation in plumage pigmentation through the combined effects of dozens of loci. The dominance and epistatic hierarchies of key loci governing this diversity are known through classical genetic studies, but their molecular identities and the mechanisms of their genetic interactions remain unknown. Here we identify protein-coding and cis-regulatory mutations in Tyrp1, Sox10, and Slc45a2 that underlie classical color phenotypes of pigeons and present a mechanistic explanation of their dominance and epistatic relationships. We also find unanticipated allelic heterogeneity at Tyrp1 and Sox10, indicating that color variants evolved repeatedly though mutations in the same genes. These results demonstrate how a spectrum of coding and regulatory mutations in a small number of genes can interact to generate substantial phenotypic diversity in a classic Darwinian model of evolution.