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.

CRISPR-Cas9 screening identified novel subtypes of cutaneous melanoma based on essential cancer genes

The Identification of the Metabolism Subtypes of Skin Cutaneous Melanoma Associated with the Tumor Microenvironment and the Immunotherapy

Analysis of immune subtypes based on immunogenomic profiling identifies prognostic signature for cutaneous melanoma

New Prognostic Biomarkers and Drug Targets for Skin Cutaneous Melanoma Via Comprehensive Bioinformatic Analysis and Validation

The molecular subtypes and clinical prognosis characteristic of tertiary lymphoid structures-related gene of cutaneous melanoma

Identification and Validation of Three Risk Models in Skin Melanoma Based on Bioinformatics Analysis

Analysis of prognostic immune and stromal- related genes in melanoma microenvironment

Integrated transcriptomic and immunological profiling reveals new diagnostic and prognostic models for cutaneous melanoma

A novel risk-scoring model developed based on tertiary lymphoid structure related genes predict the prognosis and immunotherapeutic efficacy in cutaneous melanoma.

Analyzing Prognostic Hub Genes in the Microenvironment of Cutaneous Melanoma by Computer Integrated Bioinformatics

Identification and validation of TME-related signatures to predict prognosis and response to anti-tumor therapies in skin cutaneous melanoma

Molecular Classification and Subtype-Specific Characterization of Skin Cutaneous Melanoma by Aggregating Multiple Genomic Platform Data

Identification of potential therapeutic targets for skin cutaneous melanoma on the basic of transcriptomics

Identification of Genes Related to Immune Infiltration in the Tumor Microenvironment of Cutaneous Melanoma

Identification of Immune-Related Biomarkers Associated with Tumorigenesis and Prognosis in Skin Cutaneous Melanoma.

Identification of anoikis-related subtypes and development of risk stratification system in skin cutaneous melanoma

Melanoma Molecular Subtypes and Development of Prognostic and Immunotherapy-Related Genetic Characteristics by Ferroptosis Gene Analysis

Prognostic biomarker CCR6 and its correlation with immune infiltration in cutaneous melanoma

Prognostic Signature and Immune Efficacy of M1a‐, M5c‐ and M6a‐related Regulators in Cutaneous Melanoma

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

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