Abstract:Background Melanoma is a malignant tumor of melanocytes, and the incidence has increased faster than any other cancer over the past half century. Most primary melanoma can be cured by local excision, but metastatic melanoma has a poor prognosis. Cutaneous melanoma (CM) is prone to metastasis, so the research on the mechanism of melanoma occurrence and metastasis will be beneficial to diagnose early, improve treatment, and prolong life survival. In this study, we compared the gene expression of normal skin (N), primary cutaneous melanoma (PM) and metastatic cutaneous melanoma (MM) in the Gene Expression Omnibus (GEO) database. Then we identified the key genes and molecular pathways that may be involved in the development and metastasis of cutaneous melanoma, thus to discover potential markers or therapeutic targets. Methods Three gene expression profiles ( GSE7553 , GSE15605 and GSE46517 ) were downloaded from the GEO database, which contained 225 tissue samples. R software identified the differentially expressed genes (DEGs) between pairs of N, PM and MM samples in the three sets of data. Subsequently, we analyzed the gene ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of the DEGs, and constructed a protein-protein interaction (PPI) network. MCODE was used to seek the most important modules in PPI network, and then the GO function and KEGG pathway of them were analyzed. Finally, the hub genes were calculated by the cytoHubba in Cytoscape software. The Cancer Genome Atlas (TCGA) data were analyzed using UALCAN and GEPIA to validate the hub genes and analyze the prognosis of patients. Results A total of 134, 317 and 147 DEGs were identified between N, PM and MM in pair. GO functions and KEGG pathways analysis results showed that the upregulated DEGs mainly concentrated in cell division, spindle microtubule, protein kinase activity and the pathway of transcriptional misregulation in cancer. The downregulated DEGs occurred in epidermis development, extracellular exosome, structural molecule activity, metabolic pathways and p53 signaling pathway. The PPI network obtained the most important module, whose GO function and KEGG pathway were enriched in oxidoreductase activity, cell division, cell exosomes, protein binding, structural molecule activity, and metabolic pathways. 14, 18 and 18 DEGs were identified respectively as the hub genes between N, PM and MM, and TCGA data confirmed the expression differences of hub genes. In addition, the overall survival curve of hub genes showed that the differences in these genes may lead to a significant decrease in overall survival of melanoma patients. Conclusions In this study, several hub genes were found from normal skin, primary melanoma and metastatic melanoma samples. These hub genes may play an important role in the production, invasion, recurrence or death of CM, and may provide new ideas and potential targets for its diagnosis or treatment.

Comparative quantitative proteomic analysis of melanoma subtypes, nevus-associated melanoma, and corresponding nevi

Single-Cell Characterization of the Immune Microenvironment of Melanoma Brain and Leptomeningeal Metastases.

Whole-genome landscapes of major melanoma subtypes

Network-based Molecular Subtyping of Acral Melanoma

New panel of biomarkers to discriminate between amelanotic and melanotic metastatic melanoma

Identification, Validation, And Functional Annotations Of Genome-Wide Profile Variation Between Melanocytic Nevus And Malignant Melanoma

Comprehensive Genomic Profiling Reveals Distinct Patterns of Driver Mutations and Chromosomal Alterations in Acral and Mucosal Melanomas.

A Study on the Expression of miRNA Gene Subtypes in Tumoral vs. Benign Nevi

Identification of two molecularly and prognostically distinct subtypes in acral melanoma using network prediction method

Proteogenomic insights into the biology and treatment of pan-melanoma

FIGURE 3 from Comparative Epigenetic Profiling Reveals Distinct Features of Mucosal Melanomas Associated with Immune Cell Infiltration and Their Clinical Implications

Bioinformatics Analysis to Identify the Critical Genes, Micrornas and Long Noncoding RNAs in Melanoma

Analysis of differential gene expression in human melanocytic tumour lesions by custom made oligonucleotide arrays

Biology and genetics of acquired and congenital melanocytic naevi

Figure S1 from Comparative Epigenetic Profiling Reveals Distinct Features of Mucosal Melanomas Associated with Immune Cell Infiltration and Their Clinical Implications

Deconvoluted methylation profiles discriminate between closely related melanocytic nevi

Comprehensive analysis and identification of key genes and signaling pathways in the occurrence and metastasis of cutaneous melanoma

Proteomic analysis of pleomorphic dermal sarcoma reveals a fibroblastic cell of origin and distinct immune evasion mechanisms

Single-cell RNA sequencing unveils tumor heterogeneity and immune microenvironment between subungual and plantar melanoma

Whole-Genome-Sequenzierung von Schleimhautmelanomen zeigt diverse Treiber und therapeutische Targets auf

Figure S3 from Comparative Epigenetic Profiling Reveals Distinct Features of Mucosal Melanomas Associated with Immune Cell Infiltration and Their Clinical Implications