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.     

A novel CDKN2A variant (p16L117P) in a patient with familial and multiple primary melanomas

Germline mutations in CDKN2A (p16) are commonly found in patients with family history of melanoma or personal history of multiple primary melanomas. The p16 tumor suppressor gene regulates cell cycle progression and senescence through binding of cyclin‐dependent kinases (CDK) and also regulates cellular oxidative stress independently of cell cycle control. We identified a germline missense (c.350T>C, p.Leu117Pro) CDKN2A mutation in a patient who had history of four primary melanomas, numerous nevi, and self‐reported family history of melanoma. This particular CDKN2A mutation has not been previously reported in prior large studies of melanoma kindreds or patients with multiple primary melanomas. Compared with wild‐type p16, the p16^L117P mutant largely retained binding capacity for CDK4 and CDK6 but exhibited impaired capacity for repressing cell cycle progression and inducing senescence, while retaining its ability to reduce mitochondrial reactive oxygen species. Structural modeling predicted that the Leu117Pro mutation disrupts a putative adenosine monophosphate (AMP) binding pocket involving residue 117 in the fourth ankyrin domain. Identification of this new likely pathogenic variant extends our understanding of CDKN2A in melanoma susceptibility and implicates AMP as a potential regulator of p16.     

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.     

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.     

Genomic analysis reveals low tumor mutation burden which may be associated with GNAQ/11 alteration in a series of primary leptomeningeal melanomas

Primary central nervous system melanoma is rare and characterized by a variable prognosis, and no current treatment guidelines exist. We describe the clinical course of a 70‐year‐old female patient diagnosed with primary leptomeningeal melanoma (LMN) whose case represents the diagnostic and management challenges of this tumor. Targeted genomic sequencing of 315 genes from this tumor revealed GNAQ Q209L mutation and low (4 mutations/Megabase) tumor mutation burden (TMB). Wild‐type NRAS, KIT, and BRAF were also observed. A cohort of 4,787 melanomas was subsequently analyzed to identify additional primary central nervous system melanomas, of which 10 additional tumors met pathologic criteria (0.21% of total melanoma cohort). These tumors were genomically assessed according to the same targeted sequencing panel, and 6 of the tumors were also found to harbor a GNAQ mutation. All 10 tumors had low (less than or equal to 2 mutations/Megabase) TMB indicating a potential trend between G‐protein‐coupled receptor (GPCR) alterations and low TMB in LMNs. GPCR alterations were found to significantly correlate with TMB across the cohort of 4,787 melanomas, supporting this potential finding in the limited LMN subset.     

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.     

Functional analysis of RPS27 mutations and expression in melanoma

Next‐generation sequencing has enabled genetic and genomic characterization of melanoma to an unprecedent depth. However, the high mutational background plus the limited depth of coverage of whole‐genome sequencing performed on cutaneous melanoma samples make the identification of novel driver mutations difficult. We sought to explore the somatic mutation portfolio in exonic and gene regulatory regions in human melanoma samples, for which we performed targeted sequencing of tumors and matched germline DNA samples from 89 melanoma patients, identifying known and novel recurrent mutations. Two recurrent mutations found in the RPS27 promoter associated with decreased RPS27 mRNA levels in vitro. Data mining and IHC analyses revealed a bimodal pattern of RPS27 expression in melanoma, with RPS27‐low patients displaying worse prognosis. In vitro characterization of RPS27‐high and RPS27‐low melanoma cell lines, as well as loss‐of‐function experiments, demonstrated that high RPS27 status provides increased proliferative and invasive capacities, while low RPS27 confers survival advantage in low attachment and resistance to therapy. Additionally, we demonstrate that 10 other cancer types harbor bimodal RPS27 expression, and in those, similarly to melanoma, RPS27‐low expression associates with worse clinical outcomes. RPS27 promoter mutation could thus represent a mechanism of gene expression modulation in melanoma patients, which may have prognostic and predictive implications.     

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.     

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.     

Inhibition of melanoma development in the Nras(Q61K)::Ink4a−/− mouse model by the small molecule BI‐69A11

To date, there are no effective therapies for tumors bearing NRAS mutations, which are present in 15–20% of human melanomas. Here we extend our earlier studies where we demonstrated that the small molecule BI‐69A11 inhibits the growth of melanoma cell lines. Gene expression analysis revealed the induction of interferon‐ and cell death‐related genes that were associated with responsiveness of melanoma cell lines to BI‐69A11. Strikingly, the administration of BI‐69A11 inhibited melanoma development in genetically modified mice bearing an inducible form of activated Nras and a deletion of the Ink4a gene (Nras^(Q61K)::Ink4a^?/?). Biweekly administration of BI‐69A11 starting at 10 weeks or as late as 24 weeks after the induction of mutant Nras expression inhibited melanoma development (100 and 36%, respectively). BI‐69A11 treatment did not inhibit the development of histiocytic sarcomas, which constitute about 50% of the tumors in this model. BI‐69A11‐resistant Nras^(Q61K)::Ink4a^?/? tumors exhibited increased CD45 expression, reflective of immune cell infiltration and upregulation of gene networks associated with the cytoskeleton, DNA damage response, and small molecule transport. The ability to attenuate the development of NRAS mutant melanomas supports further development of BI‐69A11 for clinical assessment.     

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.     

Concordance of somatic mutational profile in multiple primary melanomas

This study aimed to determine the frequency and concordance of BRAF and NRAS mutation in tumours arising in patients with multiple primary melanoma (MPM ). Patients with MPM managed at one of three tertiary referral centres in Melbourne, Australia, from 2010 to 2015 were included. Incident and subsequent melanomas underwent mutation testing. Cohen's kappa (κ) coefficient assessed agreement between incident and subsequent primary melanomas for both BRAF and NRAS mutation status (mutant versus wild‐type). Mutation testing of at least two primary tumours from 64 patients was conducted. There was poor agreement for both BRAF and NRAS mutation status between incident and subsequent melanomas (κ = 0.10, 95% CI ?0.10 to 0.42; κ = 0.06, 95% CI ?0.10 to 0.57, respectively). In view of the low concordance in BRAF mutation status between incident and subsequent melanomas, mutational analysis of metastatic tissue, rather than of a primary melanoma, in patients with MPM should be used to guide targeted therapy.     

Selection for a dominant oncogene and large male size as a risk factor for melanoma in the Xiphophorus animal model

Adult height is a risk factor in numerous human cancers that involve aberrant receptor tyrosine kinase (RTK) signalling. However, its importance is debated due to conflicting epidemiological studies and the lack of useful in vivo models. In Xiphophorus fishes (Platyfishes/Swordtails), a functional RTK, Xiphophorus melanoma receptor kinase (Xmrk), serves as the dominant oncogene and has been maintained for several million years despite being deleterious and in an extremely unstable genomic region. Here we show that the Xmrk genotype is positively correlated with standard length in male and female wild caught Xiphophorus cortezi sampled throughout their phylogeographic distribution. Histopathology confirms the occurrence of malignant melanomas in both sexes; however, melanoma incidence was extremely male biased. Furthermore, males collected with malignant melanomas in the field were significantly larger than both Xmrk males collected without melanomas and wildtype (Xmrk deficient) males. These results not only provide a novel selective mechanism for the persistence of the germline Xmrk oncogene but also create an innovative avenue of melanoma research within the Xiphophorus fishes. Wildlife cancer in natural systems is a growing concern, therefore, future research investigating life history characteristics associated with certain phenotypes and genotypes that predispose an individual to cancer will be fundamental to increasing our understanding of the evolutionary biology of cancer in nature as well as in humans.     

Recurrent BRAF kinase fusions in melanocytic tumors offer an opportunity for targeted therapy

BRAF is the most prevalent oncogene and an important therapeutic target in melanoma. In some cancers, BRAF is activated by rearrangements that fuse its kinase domain to 5′ partner genes. We examined 848 comparative genomic hybridization profiles of melanocytic tumors and found copy number transitions within BRAF in 10 tumors, of which six could be further characterized by sequencing. In all, the BRAF kinase domain was fused in‐frame to six N‐terminal partners. No other mutations were identified in melanoma oncogenes. One of the seven melanoma cell lines without known oncogenic mutations harbored a similar BRAF fusion, which constitutively activated the MAP kinase pathway. Sorafenib, but not vemurafenib, could block MAP kinase pathway activation and proliferation of the cell line at clinically relevant concentrations, whereas BRAF^V600E mutant melanoma cell lines were significantly more sensitive to vemurafenib. The patient from whom the cell line was derived showed a durable clinical response to sorafenib.     

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.     

Sporadic naturally occurring melanoma in dogs as a preclinical model for human melanoma

Melanoma represents a significant malignancy in humans and dogs. Different from genetically engineered models, sporadic canine melanocytic neoplasms share several characteristics with human disease that could make dogs a more relevant preclinical model. Canine melanomas rarely arise in sun‐exposed sites. Most occur in the oral cavity, with a subset having intra‐epithelial malignant melanocytes mimicking the in situ component of human mucosal melanoma. The spectrum of canine melanocytic neoplasia includes benign lesions with some analogy to nevi, as well as invasive primary melanoma, and widespread metastasis. Growing evidence of distinct subtypes in humans, differing in somatic and predisposing germ‐line genetic alterations, cell of origin, epidemiology, relationship to ultraviolet radiation and progression from benign to malignant tumors, may also exist in dogs. Canine and human mucosal melanomas appear to harbor BRAF, NRAS, and c‐kit mutations uncommonly, compared with human cutaneous melanomas, although both species share AKT and MAPK signaling activation. We conclude that there is significant overlap in the clinical and histopathological features of canine and human mucosal melanomas. This represents opportunity to explore canine oral cavity melanoma as a preclinical model.