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.

Screening and identification of potential biomarkers and therapeutic drugs in melanoma via integrated bioinformatics analysis

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

Identification of potential biomarkers with colorectal cancer based on bioinformatics analysis and machine learning

Bioinformatic Analysis Identifies Potential Key Genes in the Pathogenesis of Melanoma

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

Re-Exploring Biomarkers and Therapeutic Targets in Primary Melanoma Patients: Insights from Network-Based Analysis of Microarray Data

Bioinformatics Analysis Identifies MicroRNAs and Target Genes Associated with Prognosis in Patients with Melanoma

Screening and identification of potential prognostic biomarkers in metastatic skin cutaneous melanoma by bioinformatics analysis

Identification of significant genes with a poor prognosis in skin cutaneous malignant melanoma based on a bioinformatics analysis

Screening Genes Associated with Melanoma Using a Combined Analysis of Mrna and Methylation Microarray

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

Bioinformatical analysis identified significant genes in the PI3K-Akt signaling pathway relating with poor prognosis of skin cutaneous melanoma

Key Candidate Genes of STAT1 and CXCL10 in Melanoma Identified by Integrated Bioinformatical Analysis

Identification of Keratinocyte Differentiation-Involved Genes for Metastatic Melanoma by Gene Expression Profiles

Detailed Analysis of Molecular Mechanisms in Primary and Metastatic Melanoma.

Identification of Core Genes and Pathways in Melanoma Metastasis via Bioinformatics Analysis

Identification of Epigenetically Altered Genes and Potential Gene Targets in Melanoma Using Bioinformatic Methods

Gene Identification and Potential Drug Therapy for Drug-Resistant Melanoma with Bioinformatics and Deep Learning Technology

Network-based Co-Expression Analysis for Exploring the Potential Diagnostic Biomarkers of Metastatic Melanoma

Identification of Hub Genes Associated With Melanoma Development by Comprehensive Bioinformatics Analysis

Prognostic Genes of Melanoma Identified by Weighted Gene Co-Expression Network Analysis and Drug Repositioning Using a Network-Based Method.