Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 644.
Published online 2017 September 21. doi:  10.1038/s41467-017-00758-3
PMCID: PMC5608693

Combined activation of MAP kinase pathway and β-catenin signaling cause deep penetrating nevi


Deep penetrating nevus (DPN) is characterized by enlarged, pigmented melanocytes that extend through the dermis. DPN can be difficult to distinguish from melanoma but rarely displays aggressive biological behavior. Here, we identify a combination of mutations of the β-catenin and mitogen-activated protein kinase pathways as characteristic of DPN. Mutations of the β-catenin pathway change the phenotype of a common nevus with BRAF mutation into that of DPN, with increased pigmentation, cell volume and nuclear cyclin D1 levels. Our results suggest that constitutive β-catenin pathway activation promotes tumorigenesis by overriding dependencies on the microenvironment that constrain proliferation of common nevi. In melanoma that arose from DPN we find additional oncogenic alterations. We identify DPN as an intermediate stage in the step-wise progression from nevus to melanoma. In summary, we delineate specific genetic alterations and their sequential order, information that can assist in the diagnostic classification and grading of these distinctive neoplasms.


Common nevi are benign melanocytic neoplasms located superficially in the skin. Nevus cells become smaller and less pigmented as their distance from the surface epithelium increases. This characteristic is termed ‘maturation’ and is an important diagnostic feature to distinguish benign nevi from melanoma1, 2. In 1989 Helwig and colleagues described a neoplasm – deep penetrating nevus (DPN) – as an exception to this rule in which the constituent melanocytes maintain their cell size and pigmentation throughout the dermis3. Before their initial description, these microscopic characteristcs resulted in the misdiagnosis of almost one-third of DPN as melanoma. While the establishment of DPN as a diagnostic entity has reduced the number of cases falsely classified as melanoma, tumors with the characteristic phenotypic features can present diagnostic challenges, and occasional tumors classified as DPN metastasize with fatal outcome4.

While the majority of common nevi are clonal proliferations of BRAF V600E mutant melanocytes5, 6, the genetic drivers of DPN are largely unknown. DPN have overlapping phenotypic features with blue nevi in that they are mainly dermal based and contain numerous melanophages7. However, they lack the GNAQ and GNA11 mutations that are found in the majority of blue nevi8. DPN also have overlapping phenotypic features with Spitz nevi in that they have large melanocytes with abundant cytoplasm and have previously been shown to harbor HRAS mutations in 6% of cases8.To investigate oncogenic alterations in DPN, we sequenced DNA from 18 DPN and 10 common nevi, using a platform that includes exons of several hundred cancer genes and introns involved in oncogenic rearrangements. We sequenced additional nevi with overlapping features of DPN and blue nevi. The resulting sequence data was analyzed for single-nucleotide variants, insertion/deletions (indels), structural variants and copy number changes. We identify MAPK activating mutations in combination with β-catenin activating mutations as a characteristic feature of DPN. We demonstrate that CTNNB1 mutations arise after BRAF mutation in DPN arising from common acquired nevi and confer the characteristic phenotypic features of DPN. In DPN-like melanomas, we identify the presence of additional genetic alterations indicating that DPN are genetically intermediate between benign nevi and DPN-like melanomas.


β-catenin and MAPK pathway mutations define DPN

Cases with typical histological features of DPN were genetically distinct from common nevi and nevi whose features overlapped between DPN and blue nevus. 17 out of 18 (94%) DPN had activating mutations of the β-catenin pathway, mostly affecting CTNNB1 (16 cases) with rare APC mutations (Fig. 1a and Supplementary Table 1), whereas no common nevi harbored β-catenin activating mutations (P = 0.007214, binomial test). By contrast, nevi with overlapping features of DPN and blue nevi were devoid of CTNNB1 and APC mutations (P = 0.03, binomial test), instead harboring GNAQ activating mutations in 5 out of 7 (71%) cases (P = 0.02, binomial test). Thus, the latter are better classified as blue nevi based on their genetic profile. Only one DPN harbored GNAQ Q209L without mutation in CTNNB1 or APC. This tumor also had an IDH1 R132C mutation that was not present in the other GNAQ mutant tumors in our series. Perhaps this mutation led to epigenetic changes that altered histopathological features contributing to the apparent morphological misclassification of this lesion that genetically represents a blue nevus.

Fig. 1
Combined MAPK and β-catenin pathway activating mutations define deep penetrating nevi a β-catenin pathway mutations, affecting CTNNB1 or APC, and MAPK pathway mutations co-occur in deep penetrating nevi (DPN). Nevi with overlapping features ...

Most β-catenin pathway activating mutations were missense mutations in exon 3 of CTNNB1 that disrupt phosphorylation of β-catenin and its ubiquitin-mediated degradation, leading to increased β-catenin levels (Fig. 1b)911. In one DPN, two inactivating mutations of APC were present, representing an alternative mechanism of β-catenin activation.

16 out of 18 DPN had mutations of the mitogen-activated protein kinase (MAPK) pathway in addition to the activating mutations of the β-catenin pathway, affecting BRAF, MAP2K1 and HRAS in a mutually exclusive pattern. 9 out of 18 DPN (50%) had BRAF mutations: BRAF V600E in six cases, BRAF A598_T599insI, BRAF K601E and AKAP9-BRAF fusion in one case each. One DPN harbored an HRAS Q61R mutation. MAP2K1 alterations were present in 6 of 18 DPN (33%) consisting of small indels near codons 102 to 107, within helix C of MEK1, which is displaced by the activation loop in the active conformation12, 13. In one case, a small deletion affected the negative regulatory region (Fig. 1c)13. Similar MAP2K1 mutations occur in Langerhans cell histiocytosis and hairy cell leukemia, also in a mutually exclusive pattern with BRAF V600E 14, 15.

Three of the five MAP2K1 indels have been previously observed in melanoma or other cancers, and two have been experimentally demonstrated to constitutively activate MAP2K114, 16, 17. We expressed 4 of the observed MAP2K1 mutations in 293FT cells and found that they activate MAP kinase signaling, which could be inhibited by the MEK inhibitor trametinib (Fig. 1d). Notably, the MAP2K1 mutations present in DPN are distinct from those at positions P124 and E203, which often occur in combination with BRAF or NRAS mutations in melanoma18, 19.

By contrast, common nevi harbored BRAF V600E in 7/10, NRAS Q61R in 3/10, with no additional identified mutations. Thus DPN are genetically defined by the combination of MAPK activating and β-catenin activating mutations and are distinct from common and blue nevi.

To further investigate the role of CTNNB1 mutation in DPN, we used immunohistochemistry to evaluate β-catenin expression level and localization in DPN and common nevus controls. In common nevi, nuclear β-catenin was expressed in superficial melanocytes situated near the epidermis, whereas the smaller melanocytes in the deeper portions of the nevus had considerably lower expression levels, consistent with prior reports (Fig. 2)20. Only deep dermal melanocytes situated near hair follicles also had increased levels of nuclear β-catenin, suggesting that proximity to epithelial structures - epidermis or hair follicle epithelium - may increase β-catenin signaling (Supplementary Fig. 1). In DPN, β-catenin staining was more intense and by contrast did not diminish with distance from epithelia (Fig. 2, Supplementary Fig. 2). We also measured the expression of AXIN2, a marker of β-catenin transcription21, by RNA in situ hybridization. In DPN, melanocytes demonstrated constant levels of AXIN2 expression throughout the neoplasm, whereas in common nevi AXIN2 expression decreased with distance from epithelium (Supplementary Table 2 and Supplementary Fig. 3).

Fig. 2
Cell size and pigmentation decrease with β-catenin and cyclin D1 levels in common nevi in contrast with deep penetrating nevi. a BRAFV600E mutant common nevus. Scale bars: 300 microns, inset 12.5 microns. b BRAFV600E/CTNNB1T41A mutant DPN. Scale ...

Cyclin D1 is a direct transcriptional target of β-catenin22 and is frequently amplified in melanoma23. Melanocytes in DPN demonstrated strong and uniform nuclear expression of cyclin D1 by immunohistochemistry. In contrast, immunoreactivity for cyclin D1 in common nevi was limited to melanocytes near epithelia, mimicking the expression pattern of CTNNB1 (Fig. 2, Supplementary Figs. 1, 2).

CTNNB1 mutations and confer phenotypic characteristics of DPN

Some DPN arise juxtaposed with a common nevus2426. We hypothesized genetic progression of a BRAF V600E melanocyte within a common nevus by acquisition of a β-catenin activating mutation results in an adjacent DPN. In 66 out of 68 (97%) bi-phenotypic cases the DPN displayed increased nuclear β-catenin and cyclin D1 compared to the adjacent common nevus. β-catenin and cyclin D1 levels in the DPN remained uniform, whereas they decreased in the adjacent common nevus with distance from the epithelium (Fig. 3). We genotyped BRAF, MAP2K1, and CTNNB1 in the DPN of 11 biphasic lesions. Activating CTNNB1 mutations were identified in 10 out of 11 DPN within common nevi, usually in combination with BRAF V600E (Supplementary Table 3). In 2 cases we genotyped the adjacent common nevus component and found identical BRAF V600E mutations as in the adjacent DPN but no CTNNB1 mutations. These results indicate that common nevi can evolve to DPN by acquiring β-catenin activating mutations.

Fig. 3
Secondary CTNNB1 mutations in common acquired nevi with BRAF mutations confer the phenotypic characteristics of deep penetrating nevus. a Clinical presentation of ‘combined DPN’ with common nevus and DPN areas. left A dark papule grew ...

We transduced immortalized mouse melanocytes (melan-a) with BRAF V600E either alone or in combination with CTNNB1 S33 F. Melanocytes expressing BRAF V600E and CTNNB1 S33F were larger, more pigmented, and expressed higher levels of cyclin D1 than melanocytes expressing BRAF V600E alone (Fig. 4), recapitulating the characteristic features of DPN. These experimental findings indicate that CTNNB1 mutations are sufficient to induce the phenotypic changes characteristic of DPN in BRAF mutant melanocytes.

Fig. 4
Mutant CTNNB1 confers increased cell volume, pigmentation and cyclin D1 levels in vitro. a Stably transduced melanocytes (melan-a) with BRAFV600E and CTNNB1S33F demonstrate increased β-catenin activity as measured by TOP flash luciferase expression ...

DPN-like melanomas harbor additional oncogenic mutations

Rare DPN metastasize with lethal outcomes4, 27. It is not known whether DPN that metastasize are genetically different from those that do not. We genotyped two metastases that originated from previously diagnosed DPN and identified that in addition to the characteristic combination of MAPK activating (NRAS Q61R or MAP2K1 I103D/E203K) and CTNNB1 mutation (Fig. 5, Supplementary Fig. 4, Supplementary Table 4) both harbored multiple DNA copy number alterations, a hallmark of melanoma28. One also harbored TERT promoter (C228T) and TP53 T86M mutations.

Fig. 5
Additional genetic alterations in DPN-like melanoma. Distant metastases of DPN harbored a combination of MAPK and b-catenin pathway activating mutations in addition to oncogenic point mutations and copy number alterations, left. Half of primary melanomas ...

Some melanomas have features of DPN such as melanocytes with pigmented cytoplasm and numerous stromal melanophages. To determine whether DPN-like melanomas share genetic characteristics with DPN, we analyzed 6 cases. One melanoma harbored a GNAQ R183Q activating mutation and truncating BAP1 mutation with loss of heterozygosity and thus was re-classified as a blue nevus-like melanoma29. The five remaining melanomas all harbored MAPK pathway activating mutations in BRAF or NRAS and three had activating mutations in the β-catenin pathway (Fig. 5). In all cases, copy number alterations and additional oncogenic alterations were identified, including TERT promoter mutation or biallelic loss of CDKN2A.


In summary, our study shows that DPN result from the combined mutational activation of the MAP kinase and β-catenin pathways. Some DPN evolve from pre-existing nevi by acquiring a β-catenin pathway mutation secondary to the initiating MAP kinase pathway mutation. Prior work demonstrates that β-catenin signaling in melanocytes results in reduced T-cell response30, which could contribute to the deeply invasive phenotype of DPN. The differences in β-catenin pathway activation between common nevi and DPN reveal another possible mechanism by which β-catenin pathway activation contributes to neoplastic progression. In DPN, melanocytes remain large and pigmented throughout the entire neoplasm, whereas in common nevi without mutational activation of the β-catenin pathway these cytomorphological features are constrained to areas in proximity to epithelia, and melanocytes become smaller and less pigmented with increasing distance from epithelia. Epithelia secrete WNTs31, which may explain the decrease in β-catenin signaling with distance from epithelia. A dependence on WNT signaling may prevent common nevi from invading the deeper dermis. Constitutive β-catenin pathway activation in DPN may override this limitation, allowing these tumors to ‘deeply penetrate’.

Our data further indicates that β-catenin and MAPK pathway mutations alone are insufficient to fully transform melanocytes, and that additional mutations such as immortalizing mutations of TERT and loss of CDKN2A are required for DPN to progress to melanoma, (Fig. 6) consistent with prior work in mouse models32, 33. Our results identify DPN as an intermediate melanocytic neoplasm, with a progression stage positioned between benign nevus and DPN-like melanoma.

Fig. 6
Model of step-wise progression in DPN-like melanoma. BRAF mutation leads to a common nevus. Subsequent CTNNB1 mutation results in the phenotypic switch to DPN. Additional genetic alterations result full transformation to DPN-like melanoma


Clinical Case Selection and Histopathologic Classification

We searched the archives at the University of California Dermatopathology Unit and the personal consultation files of AdlF for cases annotated as having features of DPN. The study was approved by the Committee on Human Research of the University of California, San Francisco, with a waiver of patient consent. Consent to publish patient photos was obtained. The histopathology of all cases were reviewed by at least 2 expert dermatopathologists (AdlF, BCB, IY, ED) and categorized as either DPN, nevus with overlapping features of DPN and blue nevus, ‘combined nevi’ with common nevus and DPN cytomorphology, or melanoma with features of DPN (markedly pigmented with defective maturation and occasionally deep invasion). Samples were excluded if the amount of tumor did not fall into one of the above predefined histopathologic subtypes or the tumor volume was deemed insufficient for further analysis. The number of samples in our cohort for sequencing was limited by the availability of DPN with enough residual material in our archives.

DNA extraction and targeted next-generation sequencing

Areas of tumor were microdissected from 20 µm sections of formalin-fixed, paraffin-embedded (FFPE) tumor. After deparaffinization by washing with Safeclear and ethanol DNA was extracted by phenol chloroform extraction. Multiplex library preparation was performed using the Ovation Ultralow Library System (NuGEN, San Carlos, CA, p/n 0331-32), Nextflex (Bioo Scientific, Austin, TX, p/n No. 5140-53) or Kapa Hyper Prep Kit (Kapa Biosystems, Wilmington, MA, p/n 07962363001) according to the manufacturer’s specifications with up to 200 ng of sample DNA. Hybridization-capture of pooled libraries was performed using custom-designed bait libraries (Nimblegen SeqCap EZ Choice, p/n 06588786001) including the exons of BRAF, NRAS, HRAS, KIT, GNAQ and select introns of BRAF. Three versions of the custom-designed bait libraries were used, and the genes targeted and their overlap between versions are detailed in Supplementary Data 16. The target intervals cover mostly exonic but also some intronic and untranslated regions of 293 (version 1), 538 (version 2), and 480 (version 3) target genes. The target genes were curated to comprise common cancer genes with particular relevance to melanoma.

Captured libraries were sequenced as paired-end 100 bp reads on a HiSeq 2000 or HiSeq 2500 instrument (Illumina). Sequence reads were mapped to the reference human genome (hg19) using the Burrows-Wheeler aligner (BWA)34. Recalibration of reads was performed using the Genome Analysis Toolkit (GATK)35. Mutation calling was performed with FreeBayes36. Coverage and sequencing statistics were determined using Picard CalculateHsMetrics and Picard CollectInsertSizeMetrics37 (Supplementary Table 5). Variant annotation was performed with Annovar38 and variants with frequency in the 1000 genomes or exome sequencing project datasets of > 0.001 were excluded from further analysis. For fusion detection, read pairs with one or more reads unaligned, insert sizes greater than 1000 bp, or with soft clipping of at least one read were extracted and re-aligned using BWA-SW39 and used as input to CREST40. Structural variants predicted by CREST were reviewed by visual inspection in the Integrative Genomics Viewer41. We predicted the resulting fusion transcript by joining the exon directly upstream from the genomic breakpoint with the exon directly downstream. Predicted protein sequences were then determined from the predicted transcripts. Copy number analysis was performed using CNVkit42.

Statistical analysis

We used the binomial test with the null hypothesis that β-catenin activating mutations are equally distributed between DPN and common nevi, and DPN and nevi with features of DPN and blue nevus. We also tested the null hypothesis that GNAQ mutations distinguish nevi with features of DPN and blue nevus from DPN.

Plasmid construction

Total RNA from 293FT cells (purchased from Life Technologies) was reverse transcribed with random primers, MAP2K1 cDNA was amplified with primer pairs 5′-caccatgcccaagaagaagccgac-2′ and 5′-ttagacgccagcagcatggg, and cloned into pENTR™/D-TOPO vector (ThermoFisher, p/n K240020). BRAF V600E and CTNNB1 cDNAs were cloned into pENTR™/D-TOPO using primer pairs 5′-caccatggcggcgctgagc and 5′-tcagtggacaggaaacgcacc-3′, or 5′-caccatggctactcaagctgatttgatgg-3′ and 5′-ttacaggtcagtatcaaaccaggc-3′, respectively, and subsequently cloned into pLenti6.3/V5-DEST (i.e., BRAF, MAP2K1, and GFP constructs) or pLenti CMV/TO Puro DEST (670-1), a gift from Eric Campeau (Addgene plasmid #17293), (CTNNB1 construct)43. CTNNB1 S33F and MAP2K1 mutations were introduced by site-directed mutagenesis using the following primer pairs: CTNNB1 S33F: 5′-agcaacagtcttacctggactttggaatccattctg-3′ and 5′-cagaatggattccaaagtccaggtaagactgttgct; MAP2K1 Q58_E62del, 5-ctttcttacccagaagctgaaggatgacgact-3′ and 5′-agtcgtcatccttcagcttctgggtaagaaag; MAP2K1 E102_I103del, 5′-aaagctaattcatctgaaacccgcaatccgga-3′ and 5′-tccggattgcgggtttcagatgaattagcttt-3′; MAP2K1 I103_K104del, 5′-aagctaattcatctgg aacccgcaatccggaa-3′ and 5′-ttccggattgcgggttccagatgaattagctt-3′; and MAP2K1 P105_I107delinsS, 5′-tcatctggagatcaaatcccggaaccagatcata-3′ and 5′-tatgatctggttccgggatttgatctccagatga-3′.

Drug studies

293FT cells were transfected with MAP2K1 wild-type and mutant expression constructs using Jetprime (Polyplus-transfection, Illkirch) and 2 µg plasmid DNA as recommended by the manufacturer. Twenty four hours after transfection, the cells were treated for 4 h with 100 nM of trametinib (Selleckchem).


Immunohistochemical analysis was performed on archival FFPE tumor specimens using the following antibodies: CTNNB1 clone B-catenin 1 (Dako catalog# M3539 or IR702), CCND1 clone EP-12 (Dako, catalog #IR083) or clone SP4 (Thermo Scientific, catalog #RM-9104-R7).

RNA in situ hybridization

Formalin-fixed, paraffin-embedded tissue sections were pretreated with 13 min of boiling followed by 12 min of protease digestion (RNAscope Protease Plus). Hybridization for hs-Axin2 (catalog #400249) was performed in parallel with DapB (catalog # 310043) and PPIB (catalog # 313901) as negative and positive controls respectively. Scoring for each case was performed in melanocytes near the epidermis and at a distance from epithelium (epidermis and adnexal epithelium). A semi-quantitative scoring method was used for each area (Supplementary Table 6)44.

Sanger sequencing

Sanger sequencing was performed for CTNNB1 exon 3, BRAF exon 15, MAP2K1 exons 2 and 3. After selection of the relevant paraffin blocks or hematoxylin-eosin stained sections, genomic DNA was extracted from macrodissected FFPE material (using the QIAamp DNA FFPE extraction kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s guidelines. Sanger sequencing analysis of BRAF exon 15, CTNNB1 exon 3, and MAP2K1 exons 2 and 3 were performed.






Generation of stably transduced cell lines

Melan-a cells were generously provided by Dr Dorothy C. Bennett (St George’s Hospital, University of London, London, UK)45 and maintained in glutamine-containing RPMI-1640 supplemented with 10% heat-inactivated fetal bovine serum, 200 nM of 12-O-tetradecanoylphorbol-13-acetate (TPA), penicillin (100 units/ml) and streptomycin (50 mg/ml). 293FT cells were maintained in DME-H21 medium containing 10% heat-inactivated fetal bovine serum, MEM Non-Essential Amino Acids (0.1 mM), sodium pyruvate (1 mM), penicillin (100 units/ml) and streptomycin (50 mg/ml). Lentiviruses were produced by co-transfection of pLenti6.3/V5-DEST or pLenti CMV/TO Puro DEST (670-1) plasmid constructs with packaging vectors pCMV-VSV-G and pCMV delta R8.2 in 293FT cells. The supernatant containing viral particles were filtered on 0.45-μm polyvinylidene difluoride (PVDF) membranes and used to transduce melan-a cells in the presence of 10 µg/ml polybrene (Santa Cruz Biotechnology). Cells were selected using either 5 µg/ml blasticidin or 1 µg/ml puromycin for 1-3 weeks. Mycoplasma contamination was tested using MycoSensor PCR Assay Kit (Agilent Technologies).

Western blotting

Cells were lysed in RIPA buffer supplemented with Halt protease and phosphatase inhibitor cocktail (Thermo Scientific) and 5 mM EDTA, pH 8. Protein samples (20 µg) were electrophoresed in a NuPAGE 4–12% gradient Bis-Tris gel, transferred to a PVDF membrane. Antibodies were incubated in 5% milk-Tris-buffered saline, 0.1% tween-20 (TBST), washed in TBST, and bands were visualized by an enhanced chemiluminescent reagent and autoradiography. The antibodies used were: phospho-ERK (Cell Signaling, #9101;;, CTNNB1 (BD Biosciences, #610153;, BRAF (Santa Cruz Biotechnology, sc166;, MEK (Santa Cruz Biotechnology, sc436;, MEK (Santa Cruz Biotechnology, sc436;, phospho-MEK (Cell Signaling, #9121;;, CCND1 (Santa Cruz Biotechnology, sc718;, and HSP60 (Santa Cruz Biotechnology, sc1722). All antibodies were used at 1:1000 dilution. Complete blots in Supplementary Fig. 5.

TOPflash assay

Stably transduced Melan-a cell lines cultured in 12-well plates were transfected with 1 ug of β-catenin reporter construct (M50 Super 8x TOPFlash, a gift from Randall Moon, Addgene plasmid, #12456)46 and renilla control reporter (pRL-CMV Vector, Promega, E2231). At 24 h after transfection, luciferase activity was assayed in 10 μl of lysate using the Dual-Luciferase reporter assay system (Promega, E1910) and a Glomax-Multi luminometer (Promega).

Melanin quantitation

Cells were trypsinized and centrifuged and the resulting cell pellet was dissolved in 1N NaOH-10% DMSO, and incubated at 74 °C for 2 h. The absorbance was measured at 420 nm. Melanin content was determined using a standard curve of melanin from Sepia officinalis (Sigma, #M2649).

HoloMonitor imager measurements

Stably transduced cell lines were plated in 6-well plates (Griener EK-27160) with 20,000 cells/well in serum free conditions. Digital holograms of the cells were generated 24–48 h after plating when they were at approximately 30% cell density using a HoloMonitor M4 Digital Holography Cytometer (Phase Holographic Imaging PHI AB, Lund, Sweden). Cell volume was calculated using Hstudio M4 software (Phase Holographic Imaging PHI AB, Lund, Sweden)47.

Data availability

Sequence data that support the findings of this study have been deposited in the Sequence Read Archive (SRA), with the accession code BioProject ID: PRJNA384770 ( All other remaining data are available within the Article and Supplementary Files, or available from the authors upon request.

Electronic supplementary material


We thank Richard Jordan for critical review of the manuscript and Tyrrell Nelson, Alexander Gagnon and Arthur DeLance for excellent technical assistance. We thank Dr Pauline Guyot for providing clinical and dermoscopic photographs. We thank Dr Andy Chien for generously providing reagents.

Author contributions

Author contributions

I.Y., B.C.B. and A.d.l. F. designed the study, critically reviewed the data and wrote the manuscript. I.Y., A.H.S., V.H., D.P. and A.d.l.F. performed DNA sequence analysis. U.E.L., M.K.T., A.J., X.C. and R.L.J. performed in vitro studies. I.Y., B.C.B., A.d.l.F., E.D., L.C., P.E.L. and T.H.M. identified and reviewed the histopathological features of the study cohort. I.Y., A.d.l.F. and E.D. performed the data collection and analysis of immunohistochemical and in situ hybridization results.


Competing interests

The authors declare no competing financial interests.


Ursula E. Lang, Emeline Durieux and Meng Kian Tee contributed equally to this work.

Boris C. Bastian and Arnaud de la Fouchardière jointly supervised this work.

Electronic supplementary material

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-00758-3.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Barnhill, R. L., Piepkorn, M. & Busam, K. J. Pathology of Melanocytic Nevi and Malignant Melanoma (Springer, 2004).
2. Massi, G. & LeBoit, P. E. Histological Diagnosis of Nevi and Melanoma (Springer, 2014).
3. Seab JA, Graham JH, Helwig EB. Deep penetrating nevus. Am. J. Surg. Pathol. 1989;13:39–44. doi: 10.1097/00000478-198901000-00005. [PubMed] [Cross Ref]
4. Cerroni L, et al. Melanocytic tumors of uncertain malignant potential. Am. J. Surg. Pathol. 2010;34:314–326. doi: 10.1097/PAS.0b013e3181cf7fa0. [PubMed] [Cross Ref]
5. Pollock PM, et al. High frequency of BRAF mutations in nevi. Nat. Genet. 2003;33:19–20. doi: 10.1038/ng1054. [PubMed] [Cross Ref]
6. Yeh I, von Deimling A, Bastian BC. Clonal BRAF mutations in melanocytic nevi and initiating role of BRAF in melanocytic neoplasia. J. Natl. Cancer Inst. 2013;105:917–919. doi: 10.1093/jnci/djt119. [PMC free article] [PubMed] [Cross Ref]
7. Murali R, McCarthy SW, Scolyer RA. Blue nevi and related lesions: a review highlighting atypical and newly described variants, distinguishing features and diagnostic pitfalls. Adv. Anat. Pathol. 2009;16:365–382. doi: 10.1097/PAP.0b013e3181bb6b53. [PubMed] [Cross Ref]
8. Bender RP, McGinniss MJ, Esmay P, Velazquez EF, Reimann JD. Identification of HRAS mutations and absence of GNAQ or GNA11 mutations in deep penetrating nevi. Mod. Pathol. Off. J. US Can. Acad. Pathol. Inc. 2013;26:1320–1328. [PubMed]
9. Liu C, et al. Control of β-catenin phosphorylation/degradation by a dual-kinase mechanism. Cell. 2002;108:837–847. doi: 10.1016/S0092-8674(02)00685-2. [PubMed] [Cross Ref]
10. Rubinfeld B, et al. Stabilization of β-catenin by genetic defects in melanoma cell lines. Science. 1997;275:1790–1792. doi: 10.1126/science.275.5307.1790. [PubMed] [Cross Ref]
11. Provost E, et al. Functional correlates of mutation of the Asp32 and Gly34 residues of beta-catenin. Oncogene. 2004;24:2667–2676. doi: 10.1038/sj.onc.1208346. [PubMed] [Cross Ref]
12. Ohren JF, et al. Structures of human MAP kinase kinase 1 (MEK1) and MEK2 describe novel noncompetitive kinase inhibition. Nat. Struct. Mol. Biol. 2004;11:1192–1197. doi: 10.1038/nsmb859. [PubMed] [Cross Ref]
13. Fischmann TO, et al. Crystal structures of MEK1 binary and ternary complexes with nucleotides and inhibitors. Biochemistry. (Mosc). 2009;48:2661–2674. doi: 10.1021/bi801898e. [PubMed] [Cross Ref]
14. Chakraborty R, et al. Mutually exclusive recurrent somatic mutations in MAP2K1 and BRAF support a central role for ERK activation in LCH pathogenesis. Blood. 2014;124:3007–3015. doi: 10.1182/blood-2014-05-577825. [PubMed] [Cross Ref]
15. Waterfall JJ, et al. High prevalence of MAP2K1 mutations in variant and IGHV4-34-expressing hairy-cell leukemias. Nat. Genet. 2014;46:8–10. doi: 10.1038/ng.2828. [PMC free article] [PubMed] [Cross Ref]
16. Hodis E, et al. A landscape of driver mutations in melanoma. Cell. 2012;150:251–263. doi: 10.1016/j.cell.2012.06.024. [PMC free article] [PubMed] [Cross Ref]
17. Imielinski M, et al. Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing. Cell. 2012;150:1107–1120. doi: 10.1016/j.cell.2012.08.029. [PMC free article] [PubMed] [Cross Ref]
18. Nikolaev SI, et al. Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma. Nat. Genet. 2012;44:133–139. doi: 10.1038/ng.1026. [PubMed] [Cross Ref]
19. Akbani R, et al. Genomic classification of cutaneous melanoma. Cell. 2015;161:1681–1696. doi: 10.1016/j.cell.2015.05.044. [PMC free article] [PubMed] [Cross Ref]
20. Pawlikowski JS, et al. Wnt signaling potentiates nevogenesis. Proc. Natl Acad. Sci. USA. 2013;110:16009–16014. doi: 10.1073/pnas.1303491110. [PubMed] [Cross Ref]
21. Jho E, et al. Wnt/β-Catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway. Mol. Cell. Biol. 2002;22:1172–1183. doi: 10.1128/MCB.22.4.1172-1183.2002. [PMC free article] [PubMed] [Cross Ref]
22. Shtutman M, et al. The cyclin D1 gene is a target of the beta-catenin/LEF-1 pathway. Proc. Natl Acad. Sci. USA. 1999;96:5522–5527. doi: 10.1073/pnas.96.10.5522. [PubMed] [Cross Ref]
23. Sauter ER, et al. Cyclin D1 is a candidate oncogene in cutaneous melanoma. Cancer Res. 2002;62:3200–3206. [PubMed]
24. Ball NJ, Golitz LE. Melanocytic nevi with focal atypical epithelioid cell components: a review of seventy-three cases. J. Am. Acad. Dermatol. 1994;30:724–729. doi: 10.1016/S0190-9622(08)81502-2. [PubMed] [Cross Ref]
25. Scolyer RA, Zhuang L, Palmer AA, Thompson JF, McCarthy SW. Combined naevus: a benign lesion frequently misdiagnosed both clinically and pathologically as melanoma. Pathology. 2004;36:419–427. doi: 10.1080/00313020412331283879. [PubMed] [Cross Ref]
26. Baran JL, Duncan LM. Combined melanocytic nevi: histologic variants and melanoma mimics. Am. J. Surg. Pathol. 2011;35:1540–1548. doi: 10.1097/PAS.0b013e31822e9f5e. [PubMed] [Cross Ref]
27. Luzar B, Calonje E. Deep penetrating nevus: a review. Arch. Pathol. Lab. Med. 2011;135:321–326. [PubMed]
28. Bastian BC, LeBoit PE, Hamm H, Bröcker EB, Pinkel D. Chromosomal gains and losses in primary cutaneous melanomas detected by comparative genomic hybridization. Cancer Res. 1998;58:2170–2175. [PubMed]
29. Costa S, et al. Melanomas associated with blue nevi or mimicking cellular blue nevi: clinical, pathologic, and molecular study of 11 cases displaying a high frequency of GNA11 mutations, BAP1 expression loss, and a predilection for the scalp. Am. J. Surg. Pathol. 2016;40:368–377. doi: 10.1097/PAS.0000000000000568. [PubMed] [Cross Ref]
30. Spranger S, Bao R, Gajewski TF. Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity. Nature. 2015;523:231–235. doi: 10.1038/nature14404. [PubMed] [Cross Ref]
31. Clevers H, Loh KM, Nusse R. An integral program for tissue renewal and regeneration: Wnt signaling and stem cell control. Science. 2014;346:1248012. doi: 10.1126/science.1248012. [PubMed] [Cross Ref]
32. Delmas V, et al. Beta-catenin induces immortalization of melanocytes by suppressing p16INK4a expression and cooperates with N-Ras in melanoma development. Genes Dev. 2007;21:2923–2935. doi: 10.1101/gad.450107. [PubMed] [Cross Ref]
33. Damsky WE, et al. β-catenin signaling controls metastasis in Braf-activated Pten-deficient melanomas. Cancer. Cell. 2011;20:741–754. doi: 10.1016/j.ccr.2011.10.030. [PMC free article] [PubMed] [Cross Ref]
34. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–1760. doi: 10.1093/bioinformatics/btp324. [PMC free article] [PubMed] [Cross Ref]
35. McKenna A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–1303. doi: 10.1101/gr.107524.110. [PubMed] [Cross Ref]
36. Garrison, E. & Gabor, M. Haplotype-based variant detection from short-read sequencing. arXiv 1207.3907 [q-GN],
37. Broad Institute. Picard.
38. Yang H, Wang K. Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR. Nat. Protoc. 2015;10:1556–1566. doi: 10.1038/nprot.2015.105. [PMC free article] [PubMed] [Cross Ref]
39. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26:589–595. doi: 10.1093/bioinformatics/btp698. [PMC free article] [PubMed] [Cross Ref]
40. Wang J, et al. CREST maps somatic structural variation in cancer genomes with base-pair resolution. Nat. Methods. 2011;8:652–654. doi: 10.1038/nmeth.1628. [PMC free article] [PubMed] [Cross Ref]
41. Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–192. doi: 10.1093/bib/bbs017. [PMC free article] [PubMed] [Cross Ref]
42. Talevich E, Shain AH, Botton T, Bastian BC. CNVkit: genome-wide copy number detection and visualization from targeted DNA sequencing. PLoS Comput. Biol. 2016;12:e1004873. doi: 10.1371/journal.pcbi.1004873. [PMC free article] [PubMed] [Cross Ref]
43. Campeau E, et al. A versatile viral system for expression and depletion of proteins in mammalian cells. PLoS ONE. 2009;4:e6529. doi: 10.1371/journal.pone.0006529. [PMC free article] [PubMed] [Cross Ref]
44. Yan D, et al. Elevated expression of axin2 and hnkd mRNA provides evidence that Wnt/β-catenin signaling is activated in human colon tumors. Proc. Natl Acad. Sci. USA. 2001;98:14973–14978. doi: 10.1073/pnas.261574498. [PubMed] [Cross Ref]
45. Bennett DC, Cooper PJ, Hart IR. A line of non-tumorigenic mouse melanocytes, syngeneic with the B16 melanoma and requiring a tumour promoter for growth. Int. J. Cancer J. Int. Cancer. 1987;39:414–418. doi: 10.1002/ijc.2910390324. [PubMed] [Cross Ref]
46. Veeman MT, Slusarski DC, Kaykas A, Louie SH, Moon RT. Zebrafish prickle, a modulator of noncanonical Wnt/Fz signaling, regulates gastrulation movements. Curr. Biol. CB. 2003;13:680–685. doi: 10.1016/S0960-9822(03)00240-9. [PubMed] [Cross Ref]
47. Hejna, M., Jorapur, A., Song, J. S. & Judson, R. L. High accuracy label-free classification of kinetic cell states from holographic cytometry. bioRxiv 127449 doi:10.1101/127449 (2017).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group