So far, two familial melanoma genes have been identified, accounting for a minority of genetic risk in families. Mutations in CDKN2A account for approximately 40% of familial cases1, and predisposing mutations in CDK4 have been reported in a very small number of melanoma kindreds2. To identify other familial melanoma genes, here we conducted whole-genome sequencing of probands from several melanoma families, identifying one individual carrying a novel germline variant (coding DNA sequence c.G1075A; protein sequence p.E318K; rs149617956) in the melanoma-lineage-specific oncogene microphthalmia-associated transcription factor (MITF). Although the variant co-segregated with melanoma in some but not all cases in the family, linkage analysis of 31 families subsequently identified to carry the variant generated a log odds ratio (lod) score of 2.7 under a dominant model, indicating E318K as a possible intermediate risk variant. Consistent with this, the E318K variant was significantly associated with melanoma in a large Australian case–control sample. Likewise, it was similarly associated in an independent case–control sample from the United Kingdom. In the Australian sample, the variant allele was significantly over-represented in cases with a family history of melanoma, multiple primary melanomas, or both. The variant allele was also associated with increased naevus count and non-blue eye colour. Functional analysis of E318K showed that MITF encoded by the variant allele had impaired sumoylation and differentially regulated several MITF targets. These data indicate that MITF is a melanoma-predisposition gene and highlight the utility of whole-genome sequencing to identify novel rare variants associated with disease susceptibility.

Selumetinib (AZD6244, ARRY-142886) is a selective, non–ATP-competitive inhibitor of mitogen-activated protein/extracellular signal–regulated kinase kinase (MEK)-1/2. The range of antitumor activity seen preclinically and in patients highlights the importance of identifying determinants of response to this drug. In large tumor cell panels of diverse lineage, we show that MEK inhibitor response does not have an absolute correlation with mutational or phospho-protein markers of BRAF/MEK, RAS, or phosphoinositide 3-kinase (PI3K) activity. We aimed to enhance predictivity by measuring pathway output through coregulated gene networks displaying differential mRNA expression exclusive to resistant cell subsets and correlated to mutational or dynamic pathway activity. We discovered an 18-gene signature enabling measurement of MEK functional output independent of tumor genotype. Where the MEK pathway is activated but the cells remain resistant to selumetinib, we identified a 13-gene signature that implicates the existence of compensatory signaling from RAS effectors other than PI3K. The ability of these signatures to stratify samples according to functional activation of MEK and/or selumetinib sensitivity was shown in multiple independent melanoma, colon, breast, and lung tumor cell lines and in xenograft models. Furthermore, we were able to measure these signatures in fixed archival melanoma tumor samples using a single RT-qPCR–based test and found intergene correlations and associations with genetic markers of pathway activity to be preserved. These signatures offer useful tools for the study of MEK biology and clinical application of MEK inhibitors, and the novel approaches taken may benefit other targeted therapies.

To identify microRNAs potentially involved in melanomagenesis, we compared microRNA expression profiles between melanoma cell lines and cultured melanocytes. The most differentially expressed microRNA between the normal and tumor cell lines was miR-211. We focused on this pigment-cell-enriched miRNA as it is derived from the microphthalmia-associated transcription factor (MITF)-regulated gene, TRPM1 (melastatin). We find that miR-211 expression is greatly decreased in melanoma cells and melanoblasts compared to melanocytes. Bioinformatic analysis identified a large number of potential targets of miR-211, including POU3F2 (BRN2). Inhibition of miR-211 in normal melanocytes resulted in increased BRN2 protein, indicating that endogenous miR-211 represses BRN2 in differentiated cells. Over-expression of miR-211 in melanoma cell lines changed the invasive potential of the cells in vitro through directly targeting BRN2 translation. We propose a model for the apparent non-overlapping expression levels of BRN2 and MITF in melanoma, mediated by miR-211 expression.

We sequenced 8 melanoma exomes to identify novel somatic mutations in metastatic melanoma. Focusing on the MAP3K family, we found that 24% of melanoma cell lines have mutations in the protein-coding regions of either MAP3K5 or MAP3K9. Structural modelling predicts that mutations in the kinase domain may affect the activity and regulation of MAP3K5/9 protein kinases. The position of the mutations and loss of heterozygosity of MAP3K5 and MAP3K9 in 85% and 67% of melanoma samples, respectively, together suggest that the mutations are likely inactivating. In vitro kinase assay shows reduction in kinase activity in MAP3K5 I780F and MAP3K9 W333X mutants. Overexpression of MAP3K5 or MAP3K9 mutant in HEK293T cells reduces phosphorylation of downstream MAP kinases. Attenuation of MAP3K9 function in melanoma cells using siRNA leads to increased cell viability after temozolomide treatment, suggesting that decreased MAP3K pathway activity can lead to chemoresistance in melanoma.

Our findings add to the growing support for telomere dysregulation as a key process associated with melanoma susceptibility.

To identify 'melanoma-specific' microRNAs (miRNAs) we used an unbiased microRNA profiling approach to comprehensively study cutaneous melanoma in relation to other solid malignancies, which revealed 233 differentially expressed (≥ 2 fold, p < 0.05) miRNAs. Among the top 20 most significantly different miRNAs was hsa-miR-514a-3p. miR-514a is a member of a cluster of miRNAs (miR-506-514) involved in initiating melanocyte transformation and promotion of melanoma growth.We found miR-514a was expressed in 38/55 (69%) melanoma cell lines but in only 1/34 (3%) other solid cancers. To identify miR-514a regulated targets we conducted a miR-514a-mRNA 'pull-down' experiment, which revealed hundreds of genes, including: CTNNB1, CDK2, MC1R, and NF1, previously associated with melanoma. NF1 was selected for functional validation because of its recent implication in acquired resistance to BRAFV600E-targeted therapy. Luciferase-reporter assays confirmed NF1 as a direct target of miR-514a and over-expression of miR-514a in melanoma cell lines inhibited NF1 expression, which correlated with increased survival of BRAFV600E cells treated with PLX4032. These data provide another mechanism for the dysregulation of the MAPK pathway which may contribute to the profound resistance associated with current RAF-targeted therapies.

Melanocytic behavior, survival, and proliferation are regulated through a complex system of cell-cell adhesion molecules. Pathologic changes leading to development of malignant melanoma, upset the delicate homeostatic balance between melanocytes and keratinocytes and can lead to altered expression of cell-cell adhesion and cell-cell communication molecules. Malignant transformation of melanocytes frequently coincides with loss of E-cadherin expression. We now show loss of another member of the superfamily of classical cadherins, H-cadherin (CDH13), which may be involved in the development of malignant melanoma. The provided data show that H-cadherin expression is lost in nearly 80% of the analyzed melanoma cell lines. Knockdown of H-cadherin using siRNA increases invasive capacity in melanocytes. Functional assays show that the re-expression of H-cadherin decreases migration and invasion capacity, as well as anchorage-independent growth in comparison to control melanoma cells. Furthermore, melanoma cells, which re-express H-cadherin via stable transfection show a reduction in rate of tumor growth in a nu/nu mouse tumor model in comparison to the parental control transfected cell lines. Our study presents for the first time the down-regulation of H-cadherin in malignant melanomas and its possible functional relevance in maintenance healthy skin architecture.

Here, we showed that the secretome of senescent melanoma cells drives basal melanoma cells towards a mesenchymal phenotype, with characteristic of stems illustrated by increased level of the prototype genes FN1, SNAIL, OCT4 and NANOG. This molecular reprogramming leads to an increase in the low-MITF and slow-growing cell population endowed with melanoma-initiating cell features. The secretome of senescent melanoma cells induces a panel of 52 genes, involved in cell movement and cell/cell interaction, among which AXL and ALDH1A3 have been implicated in melanoma development. We found that the secretome of senescent melanoma cells activates the STAT3 pathway and STAT3 inhibition prevents secretome effects, including the acquisition of tumorigenic properties. Collectively, the findings provide insights into how the secretome of melanoma cells entering senescence upon chemotherapy treatments increases the tumorigenicity of naïve melanoma cells by inducing, through STAT3 activation, a melanoma-initiating cell phenotype that could favor chemotherapy resistance and relapse.