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.

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

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

Exploring Tumor Immune Microenvironment and Its Associations With Molecular Characteristics in Melanoma

The combination of single-cell and RNA sequencing analysis decodes the melanoma tumor microenvironment and identifies novel T cell-associated signature genes

Pivotal factors associated with the immunosuppressive tumor microenvironment and melanoma metastasis

Comprehensive analysis of prognostic immune‐related genes in the tumor microenvironment of cutaneous melanoma

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

Integrative analysis of bulk and single-cell RNA-seq reveals the molecular characterization of the immune microenvironment and oxidative stress signature in melanoma

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

Bioinformatics data combined with single-cell analysis reveals patterns of immunoinflammatory infiltration and cell death in melanoma

The Prognostic Landscape of Adaptive Immune Resistance Signatures and Infiltrating Immune Cells in the Tumor Microenvironment of Uveal Melanoma.

Multiplexed Imaging Mass Cytometry Reveals Distinct Tumor-Immune Microenvironments Linked to Immunotherapy Responses in Melanoma

Advancing immunotherapy for melanoma: the critical role of single-cell analysis in identifying predictive biomarkers

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

Identification of immune-related genes with prognostic significance in the microenvironment of cutaneous melanoma

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

Correlation analysis of gene and immune microenvironment for 87 melanoma cases.

Bioinformatics analysis of the tumor microenvironment in melanoma - Constructing a prognostic model based on CD8+ T cell-related genes: An observational study

Detailed Analysis of Molecular Mechanisms in Primary and Metastatic Melanoma.

Single-Cell Transcriptomic Analysis Reveals the Crosstalk Propensity Between the Tumor Intermediate State and the CD8+ T Exhausted State to be Associated with Clinical Benefits in Melanoma

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