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.

Multi-scale geometric network analysis identifies melanoma immunotherapy response gene modules

Identifying Biomarkers Associated with Immunotherapy Response in Melanoma by Multi-Omics Analysis.

Machine learning-based identification of an immunotherapy-related signature to enhance outcomes and immunotherapy responses in melanoma

Sample-specific network analysis identifies gene coexpression patterns of immunotherapy response in advanced kidney cancer

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

Spatially Resolved Analyses Link Genomic and Immune Diversity and Reveal Unfavorable Neutrophil Activation in Melanoma

Using Machine Learning Algorithms to Predict Immunotherapy Response in Patients with Advanced Melanoma

Immunogenic cell death signatures from on-treatment tumor specimens predict immune checkpoint therapy response in metastatic melanoma

Dissecting immune cell stat regulation network reveals biomarkers to predict ICB therapy responders in melanoma

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

Biological knowledge graph-guided investigation of immune therapy response in cancer with graph neural network

A Gene Signature Predicting Natural Killer Cell Infiltration and Improved Survival in Melanoma Patients

A Comprehensive Proteogenomic and Spatial Analysis of Innate and Acquired Resistance of Metastatic Melanoma to Immune Checkpoint Blockade Therapies

A Cancer Associated Fibroblasts-Related Six-Gene Panel for Anti-PD-1 Therapy in Melanoma Driven by Weighted Correlation Network Analysis and Supervised Machine Learning

Tumor Immune Microenvironment Characterization Identifies Prognosis and Immunotherapy-Related Gene Signatures in Melanoma

Unveiling the Gut Microbiome Network Fingerprint: A Deep Learning Approach to Predict Immunotherapy Response in Melanoma Patients.

Multi-Scale Spatial Analysis of the Tumor Microenvironment Reveals Features of Cabozantinib and Nivolumab Efficacy in Hepatocellular Carcinoma

Integrated analysis reveals the pivotal interactions between immune cells in the melanoma tumor microenvironment

Prediction of Immunotherapy Response in Melanoma through Combined Modeling of Neoantigen Burden and Immune-Related Resistance Mechanisms

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma

Decoding Immune Heterogeneity of Melanoma and identifying immune-prognostic hub genes