Distinct genomic traits of acral and mucosal melanomas revealed by targeted mutational profiling

The incidence of melanoma is rising globally including China. Comparing to Caucasians, the incidence of non‐cutaneous melanomas is significantly higher in Chinese. Herein, we performed genomic profiling of 89 Chinese surgically resected primary melanomas, including acral (n = 54), cutaneous (n = 22), and mucosal (n = 13), by hybrid capture‐based next‐generation sequencing. We show that mucosal melanomas tended to harbor more pathogenic mutations than other types of melanoma, though the biological significance of this finding remains uncertain. Chromosomal arm‐level alterations including 6q, 9p, and 10p/q loss were highly recurrent in all subtypes, but mucosal melanoma was significantly associated with increased genomic instability. Importantly, 7p gain significantly correlated with unfavorable clinical outcomes in non‐cutaneous melanomas, representing an intriguing prognostic biomarker of those subtypes. Furthermore, focal amplification of 4q12 (KIT, KDR, and PDGFRα) and RAD51 deletion were more abundant in mucosal melanoma, while NOTCH2 amplification was enriched in acral melanoma. Additionally, cutaneous melanomas had higher mutation load than acral melanomas, while mucosal melanomas did not differ from other subtypes in mutation burden. Together, our data revealed important features of acral and mucosal melanomas in Chinese including distinctive driver mutation pattern and increased genomic instability. These findings highlight the possibilities of combination therapies in the clinical management of melanoma.     

Molecular pathogenesis of malignant melanoma: a different perspective from the studies of melanocytic nevus and acral melanoma

The Clark model for melanoma progression emphasizes a series of histopathological changes beginning from benign melanocytic nevus to melanoma via dysplastic nevus. Several models of the genetic basis of melanoma development and progression are based on this Clark’s multi-step model, and predict that the acquisition of a BRAF mutation can be a founder event in melanocytic neoplasia. However, our recent investigations have challenged this view, showing the polyclonality of BRAF mutations in melanocytic nevi. Furthermore, it is suggested that many melanomas, including acral and mucosal melanomas, arise de novo, not from melanocytic nevus. While mutations of the BRAF gene are frequent in melanomas on non-chronic sun damaged skin which are prevalent in Caucasians, acral and mucosal melanomas harbor mutations of the KIT gene as well as the amplifications of cyclin D1 or cyclin-dependent kinase 4 gene. Amplifications of the cyclin D1 gene are detected in normal-looking ‘field melanocytes’, which represent a latent progression phase of acral melanoma that precedes the stage of atypical melanocyte proliferation in the epidermis. Based on these observations, we propose an alternative genetic progression model for melanoma.     

Mutations in KIT occur at low frequency in melanomas arising from anatomical sites associated with chronic and intermittent sun exposure

In melanoma, mutations in KIT are most frequent in acral and mucosal subtypes and rarely reported in cutaneous melanomas particularly those associated with intermittent UV exposure. Conversely melanomas arising within chronic sun damaged skin are considered to harbour KIT mutations at higher rates. To characterize the frequency of KIT mutations in a representative melanoma population, 261 patients from two Australian melanoma centres were prospectively screened for mutations in exons 11, 13 and 17 of the KIT gene. A total of 257 patients had cutaneous melanoma arising from non-acral sites and four were acral melanomas. No mucosal or ocular melanomas were analysed. KIT mutations were identified in five tumours (2% of the entire cohort) including two acral melanomas. Two of the three non-acral melanomas with KIT mutations were associated with markers of chronic sun damage as assessed by the degree of skin elastosis. In the remaining cohort, 43% had chronically sun damaged skin. This report confirms that within an Australian population, KIT mutations are infrequent in cutaneous melanomas associated with both intermittent and chronic sun exposed skin.     

Kinase gene fusions in defined subsets of melanoma

Genomic rearrangements resulting in activating kinase fusions have been increasingly described in a number of cancers including malignant melanoma, but their frequency in specific melanoma subtypes has not been reported. We used break‐apart fluorescence in situ hybridization (FISH) to identify genomic rearrangements in tissues from 59 patients with various types of malignant melanoma including acral lentiginous, mucosal, superficial spreading, and nodular. We identified four genomic rearrangements involving the genes BRAF, RET, and ROS1. Of these, three were confirmed by Immunohistochemistry (IHC) or sequencing and one was found to be an ARMC10‐BRAF fusion that has not been previously reported in melanoma. These fusions occurred in different subtypes of melanoma but all in tumors lacking known driver mutations. Our data suggest gene fusions are more common than previously thought and should be further explored particularly in melanomas lacking known driver mutations.     

Exon 15 BRAF mutations are uncommon in canine oral malignant melanomas

An activating mutation in codon 599 of BRAF has been identified in approximately 60% of human cutaneous nevi and melanomas, but not melanomas of mucosal origin. The purpose of this study was to determine if BRAF mutations occur in canine oral malignant melanomas. The canine BRAF gene was first cloned from normal canine testicular cDNA, and a novel previously unreported splice variant involving exon 5 was identified during this process. To screen canine melanoma samples for BRAF mutation in codon 599, cDNA and genomic DNA were isolated from canine malignant melanoma cell lines and primary tumor samples respectively, all from cases seen at the Veterinary Medical Teaching Hospital at the University of California, Davis. Polymerase chain reaction (PCR) was performed for exon 15 using primers based at the 5 end of exon 15 and the 5 end of intron 15 and the resultant products were directly sequenced. No mutations in codon 599 or exon 15 were identified in any of the 17 samples evaluated. However, all of the melanoma cell lines expressed BRAF and demonstrated high levels of basal ERK phosphorylation suggesting that dysregulation of this pathway is present. Therefore, similar to the case with human mucosal melanomas, canine oral malignant melanomas do not possess codon 599 BRAF mutations commonly identified in human cutaneous melanomas. This finding supports the notion that melanomas arising from non-sun-exposed sites exhibit distinct mechanisms of molecular transformation.     

Clinicopathological characteristics associated with BRAFK601E and BRAFL597 mutations in melanoma

BRAF mutations at codons L597 and K601 occur uncommonly in melanoma. Clinical and pathological associations of these mutations were investigated in a cohort of 1119 patients with known BRAF mutation status. A BRAF mutation was identified in 435 patients; Mutations at L597 and the K601E mutation were seen in 3.4 and 3.2% of these, respectively. K601E melanomas tended to occur in male patients, a median age of 58 yr, were generally found on the trunk (64%) and uncommonly associated with chronically sun‐damaged (CSD) skin. BRAF L597 melanomas occurred in older patients (median 66 yr), but were associated with CSD skin (extremities or head and neck location – 73.3%, P = 0.001). Twenty‐three percent of patients with V600E‐ and 43% of patients with K601E‐mutant melanomas presented with nodal disease at diagnosis compared to just 14% of patients with BRAF wild‐type tumors (P = 0.001 and 0.006, respectively). Overall, these mutations represent a significant minority of BRAF mutations, but have distinct clinicopathological phenotypes and clinical behaviors.     

Tumor genetic heterogeneity analysis of chronic sun‐damaged melanoma

Chronic sun‐damaged (CSD) melanoma represents 10%–20% of cutaneous melanomas and is characterized by infrequent BRAF V600E mutations and high mutational load. However, the order of genetic events or the extent of intra‐tumor heterogeneity (ITH) in CSD^high melanoma is still unknown. Ultra‐deep targeted sequencing of 40 cancer‐associated genes was performed in 72 in situ or invasive CMM, including 23 CSD^high cases. In addition, we performed whole exome and RNA sequencing on multiple regions of primary tumor and multiple in‐transit metastases from one CSD^high melanoma patient. We found no significant difference in mutation frequency in melanoma‐related genes or in mutational load between in situ and invasive CSD^high lesions, while this difference was observed in CSD^low lesions. In addition, increased frequency of BRAF V600K, NF1, and TP53 mutations (p < .01, Fisher's exact test) was found in CSD^high melanomas. Sequencing of multiple specimens from one CSD^high patient revealed strikingly limited ITH with >95% shared mutations. Our results provide evidence that CSD^high and CSD^low melanomas are distinct molecular entities that progress via different genetic routes.     

Melanomas of unknown primary have a mutation profile consistent with cutaneous sun‐exposed melanoma

Melanoma of unknown primary (MUP) is an uncommon phenomenon whereby patients present with metastatic disease without an evident primary site. To determine their likely site of origin, we combined exome sequencing from 33 MUPs to assess the total rate of somatic mutations and degree of UV mutagenesis. An independent cohort of 91 archival MUPs was also screened for 46 hot spot mutations highly prevalent in melanoma including BRAF, NRAS, KIT, GNAQ, and GNA11. Results showed that the majority of MUPs exhibited high somatic mutation rates, high ratios of C>T/G>A transitions, and a high rate of BRAF (45 of 101, 45%) and NRAS (32 of 101, 32%) mutations, collectively indicating a mutation profile consistent with cutaneous sun‐exposed melanomas. These data suggest that a significant proportion of MUPs arise from regressed or unrecognized primary cutaneous melanomas or arise de novo in lymph nodes from nevus cells that have migrated from the skin.     

High frequency of PTEN mutations in nevi and melanomas from xeroderma pigmentosum patients

We examined nevi and melanomas in 10 xeroderma pigmentosum (XP) patients with defective DNA repair. The lesions had a lentiginous appearance with markedly increased numbers of melanocytes. Using laser capture microdissection, we performed DNA sequencing of 18 benign and atypical nevi and 75 melanomas (melanoma in situ and invasive melanomas). The nevi had a similar high frequency of PTEN mutations as melanomas [61% (11/18) versus 53% (39/73)]. Both had a very high proportion of UV‐type mutations (occurring at adjacent pyrimidines) [91% (10/11) versus 92% (36/39)]. In contrast to melanomas in the general population, the frequency of BRAF mutations (11%, 7/61), NRAS mutations (21%, 13/62), and KIT mutations (21%, 6/28) in XP melanomas was lower than for PTEN. Phospho‐S6 immunostaining indicated activation of the mTOR pathway in the atypical nevi and melanomas. Thus, the clinical and histological appearances and the molecular pathology of these UV‐related XP nevi and melanomas were different from nevi and melanomas in the general population.     

BRAF internal deletions and resistance to BRAF/MEK inhibitor therapy

BRAF and MEK inhibitors have improved clinical outcomes in advanced, BRAF^V600‐mutated melanomas. Acquired resistance occurs in most patients, with numerous and diverse drivers. We obtained pretreatment and progression biopsies from a patient who progressed on dabrafenib and trametinib. In addition to a preserved BRAF^V600E mutation, an internal deletion (rearrangement) of BRAF was observed in the progression sample. This deletion involved exons 2–8, which includes the Ras‐binding domain, and is analogous to previously documented BRAF fusions and splice variants known to reactivate RAS‐RAF‐MEK‐ERK signaling. In a large cohort of melanomas, 10 additional internal deletions were identified (0.4% of all melanomas; nine of which had concurrent BRAF mutations), as well as sporadically in other tumor types. Thus, we describe a novel mechanism of resistance to BRAF and MEK inhibition.     

Genomic analysis and clinical management of adolescent cutaneous melanoma

Melanoma in young children is rare; however, its incidence in adolescents and young adults is rising. We describe the clinical course of a 15‐year‐old female diagnosed with AJCC stage IB non‐ulcerated primary melanoma, who died from metastatic disease 4 years after diagnosis despite three lines of modern systemic therapy. We also present the complete genomic profile of her tumour and compare this to a further series of 13 adolescent melanomas and 275 adult cutaneous melanomas. A somatic BRAF^V600E mutation and a high mutational load equivalent to that found in adult melanoma and composed primarily of C>T mutations were observed. A germline genomic analysis alongside a series of 23 children and adolescents with melanoma revealed no mutations in known germline melanoma‐predisposing genes. Adolescent melanomas appear to have genomes that are as complex as those arising in adulthood and their clinical course can, as with adults, be unpredictable.     

NRAS and BRAF Mutations in Melanoma-Associated Nevi and Uninvolved Nevi

According to the prevailing multistep model of melanoma development, oncogenic BRAF or NRAS mutations are crucial initial events in melanoma development. It is not known whether melanocytic nevi that are found in association with a melanoma are more likely to carry BRAF or NRAS mutations than uninvolved nevi. By laser microdissection we were able to selectively dissect and genotype cells either from the nevus or from the melanoma part of 46 melanomas that developed in association with a nevus. In 25 cases we also genotyped a control nevus of the same patients. Available tissue was also immunostained using the BRAF^V600E-mutation specific antibody VE1. The BRAF^V600E mutation was found in 63.0% of melanomas, 65.2% of associated nevi and 50.0% of control nevi. No significant differences in the distribution of BRAF or NRAS mutations could be found between melanoma and associated nevi or between melanoma associated nevi and control nevi. In concordant cases immunohistochemistry showed a higher expression (intensity of immunohistochemistry) of the mutated BRAF^V600E-protein in melanomas compared to their associated nevi. In this series the presence of a BRAF- or NRAS mutation in a nevus was not associated with the risk of malignant transformation. Our findings do not support the current traditional model of stepwise tumor progression.     

Genetic and morphologic features for melanoma classification

Melanoma is comprised of biologically distinct subtypes. The defining clinical, histomorphologic, and molecular features are not fully established. This study sought to validate the association between genetic and histomorphologic features previously described and to determine their reproducibility and association with important clinical variables. Detailed clinical and histomorphologic features of 365 primary cutaneous melanomas were assessed by 11 pathologists and correlated with mutation status of BRAF and NRAS. There was substantial agreement in the quantitative assessment of histomorphologic features showing similar or better interobserver reproducibility than the established World Health Organization classification scheme. We confirmed that melanomas with BRAF mutations showed characteristic morphologic features (P < 0.0001) and metastasized more frequently to regional lymph nodes (P = 0.046). Importantly, melanomas without mutations were a heterogeneous group, with a subset having very similar clinical and morphological features as those with BRAF mutation raising the possibility that they are biologically related. Our study confirms an association between histomorphologic features, mutation status, and pattern of metastasis, providing criteria for a refined melanoma classification aimed at defining biologically homogeneous disease subgroups.     

Clinical and pathological associations of the activating RAC1 P29S mutation in primary cutaneous melanoma

Activating mutations in the GTPase RAC1 are a recurrent event in cutaneous melanoma. We investigated the clinical and pathological associations of RAC1^P29S in a cohort of 814 primary cutaneous melanomas with known BRAF and NRAS mutation status. The RAC1^P29S mutation had a prevalence of 3.3% and was associated with increased thickness (OR=1.6 P = 0.001), increased mitotic rate (OR=1.3 P = 0.03), ulceration (OR=2.4 P = 0.04), nodular subtype (OR=3.4 P = 0.004), and nodal disease at diagnosis (OR=3.3 P = 0.006). BRAF mutant tumors were also associated with nodal metastases (OR=1.9 P = 0.004), despite being thinner at diagnosis than BRAF WT (median 1.2 mm versus 1.6 mm, P < 0.001). Immunohistochemical analysis of 51 melanomas revealed that 47% were immunoreactive for RAC1. Melanomas were more likely to show RAC1 immunoreactivity if they were BRAF mutant (OR=5.2 P = 0.01). RAC1 may therefore be important in regulating the early migration of BRAF mutant tumors. RAC1 mutations are infrequent in primary melanomas but may accelerate disease progression.     

The metabolic microenvironment of melanomas: Prognostic value of MCT1 and MCT4

BRAF mutations are known drivers of melanoma development and, recently, were also described as players in the Warburg effect, while this reprogramming of energy metabolism has been identified as a possible strategy for treating melanoma patients. Therefore, the aim of this work was to evaluate the expression and prognostic value of a panel of glycolytic metabolism-related proteins in a series of melanomas. The immunohistochemical expression of MCT1, MCT4, GLUT1, and CAIX was evaluated in 356 patients presenting melanoma and 20 patients presenting benign nevi. Samples included 20 benign nevi, 282 primary melanomas, 117 lymph node and 54 distant metastases samples. BRAF mutation was observed in 29/92 (31.5%) melanoma patients and 17/20 (85%) benign nevi samples. NRAS mutation was observed in 4/36 (11.1%) melanoma patients and 1/19 (5.3%) benign nevi samples. MCT4 and GLUT1 expression was significantly increased in metastatic samples, and MCT1, MCT4 and GLUT1 were significantly associated with poor prognostic variables. Importantly, MCT1 and MCT4 were associated with shorter overall survival. In conclusion, the present study brings new insights on metabolic aspects of melanoma, paving the way for the development of new-targeted therapies.     

Naturally occurring melanomas in dogs as models for non‐UV pathways of human melanomas

Spontaneously occurring melanomas are frequent in dogs. They appear at the same localizations as in humans, i.e. skin, mucosal sites, nail matrix and eyes. They display variable behaviors: tumors at oral localizations are more frequent and aggressive than at other anatomical sites. Interestingly, dog melanomas are associated with strong breed predispositions and overrepresentation of black‐coated dogs. Epidemiological analysis of 2350 affected dogs showed that poodles are at high risk of developing oral melanoma, while schnauzers or Beauce shepherds mostly developped cutaneous melanoma. Clinical and histopathological analyses were performed on a cohort of 153 cases with a 4‐yr follow‐up. Histopathological characterization showed that most canine tumors are intradermal and homologous to human rare morphological melanomas types – ‘nevocytoid type’ and ‘animal type’‐. Tumor cDNA sequencing data, obtained from 95 dogs for six genes, relevant to human melanoma classification, detected somatic mutations in oral melanoma, in NRAS and PTEN genes, at human hotspot sites, but not in BRAF. Altogether, these findings support the relevance of the dog model for comparative oncology of melanomas, especially for the elucidation of non‐UV induced pathways.     

Identification of PLX4032‐resistance mechanisms and implications for novel RAF inhibitors

BRAF inhibitors improve melanoma patient survival, but resistance invariably develops. Here we report the discovery of a novel BRAF mutation that confers resistance to PLX4032 employing whole‐exome sequencing of drug‐resistant BRAF^V600K melanoma cells. We further describe a new screening approach, a genome‐wide piggyBac mutagenesis screen that revealed clinically relevant aberrations (N‐terminal BRAF truncations and CRAF overexpression). The novel BRAF mutation, a Leu505 to His substitution (BRAF^L505H), is the first resistance‐conferring second‐site mutation identified in BRAF mutant cells. The mutation replaces a small nonpolar amino acid at the BRAF‐PLX4032 interface with a larger polar residue. Moreover, we show that BRAF^L505H, found in human prostate cancer, is itself a MAPK‐activating, PLX4032‐resistant oncogenic mutation. Lastly, we demonstrate that the PLX4032‐resistant melanoma cells are sensitive to novel, next‐generation BRAF inhibitors, especially the ‘paradox‐blocker’ PLX8394, supporting its use in clinical trials for treatment of melanoma patients with BRAF‐mutations.     

Novel Somatic Mutations to PI3K Pathway Genes in Metastatic Melanoma

Background

BRAF^V600 inhibitors have offered a new gateway for better treatment of metastatic melanoma. However, the overall efficacy of BRAF^V600 inhibitors has been lower than expected in clinical trials, and many patients have shown resistance to the drug’s effect. We hypothesized that somatic mutations in the Phosphoinositide 3-Kinase (PI3K) pathway, which promotes proliferation and survival, may coincide with BRAF^V600 mutations and contribute to chemotherapeutic resistance.

Methods

We performed a somatic mutation profiling study using the 454 FLX pyrosequencing platform in order to identify candidate cancer genes within the MAPK and PI3K pathways of melanoma patients. Somatic mutations of theses candidate cancer genes were then confirmed using Sanger sequencing.

Results

As expected, BRAF^V600 mutations were seen in 51% of the melanomas, whereas NRAS mutations were seen in 19% of the melanomas. However, PI3K pathway mutations, though more heterogeneous, were present in 41% of the melanoma, with PTEN being the highest mutated PI3K gene in melanomas (22%). Interestingly, several novel PI3K pathway mutations were discovered in MTOR, IRS4, PIK3R1, PIK3R4, PIK3R5, and NFKB1. PI3K pathway mutations co-occurred with BRAF^V600 mutations in 17% of the tumors and co-occurred with 9% of NRAS mutant tumors, implying cooperativity between these pathways in terms of melanoma progression.

Conclusions

These novel PI3K pathway somatic mutations could provide alternative survival and proliferative pathways for metastatic melanoma cells. They therefore may be potential chemotherapeutic targets for melanoma patients who exhibit resistance to BRAF^V600 inhibitors.     

Clinical and therapeutic implications of BRAF mutation heterogeneity in metastatic melanoma

Heterogeneity of BRAF mutation in melanoma has been a controversial subject. Quantitative data on BRAF allele frequency (AF) are sparse, and the potential relationship with response to BRAF inhibitors (BRAFi) in patients with metastatic melanoma is unknown. We quantitatively measured BRAF AF in a cohort of treatment naïve metastatic melanoma samples by pyrosequencing and correlated with survival data in patients treated with BRAFi as part of their clinical care. Fifty‐two samples from 50 patients were analysed. BRAF V600E mutations were detected in 71.1% of samples followed by V600K (25%) and V600R (3.9%). There was a wide range of AF from 3.9% to 80.3% (median 41.3%). In 33 patients treated with BRAFi, there was no difference in overall or progression‐free survival when the patients were categorized into high or low AF groups. There was no correlation between AF and degree of response, and no difference in survival based on genotype.