Identification of MicroRNAs Associated with Invasive and Aggressive Phenotype in Cutaneous Melanoma by Next-generation Sequencing

A comprehensive repertoire of human microRNAs (miRNAs) that could be involved in early melanoma invasion into the dermis remains unknown. To this end, we sequenced small RNAs (18-30 nucleotides) isolated from an annotated series of invasive melanomas (average invasive depth, 2.0 mm), common melanocytic nevi, and matched normal skin (n=28). Our previously established bioinformatics pipeline identified 765 distinct mature known miRNAs and defined a set of top 40 list that clearly segregated melanomas into thin (0.75 mm) and thick (2.7 mm) groups. Among the top, miR-21-5p, let-7b-5p, let-7a-5p, miR-424-5p, miR-423-5p, miR-21-3p, miR-199b-5p, miR-182-5p, and miR-205-5p were differentially expressed between thin and thick melanomas. In a validation cohort (n=167), measured expression of miR-21-5p and miR-424-5p, not previously reported in melanoma, were significantly increased in invasive compared with in situ melanomas (P<0.0001). Increased miR-21-5p levels were significantly associated with invasive depth (P=0.038), tumor mitotic index (P=0.038), lymphovascular invasion (P=0.0036), and AJCC stage (P=0.038). In contrast, let-7b levels were significantly decreased in invasive and in situ melanomas compared with common and dysplastic nevi (P<0.0001). Decreased let-7b levels were significantly associated with invasive depth (P=0.011), Clark's level (P=0.013), ulceration (P=0.0043), and AJCC stage (P=0.011). These results define a distinct set of miRNAs associated with invasive and aggressive melanoma phenotype.