Exploration and validation of metastasis-associated genes for skin cutaneous melanoma

Skin cutaneous melanoma is a malignant and highly metastatic skin tumor, and its morbidity and mortality are still rising worldwide. However, the molecular mechanisms that promote melanoma metastasis are unclear. Two datasets (GSE15605 and GSE46517) were retrieved to identify the differentially expressed genes (DEGs), including 23 normal skin tissues (N), 77 primary melanoma tissues (T) and 85 metastatic melanoma tissues (M). Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed to explore the functions of the DEGs. We constructed protein–protein interaction network using the STRING database and Cytoscape software. Using the cytoHubba plugin of Cytoscape, we identified the most significant hub genes by five analytical methods (Degree, Bottleneck, MCC, MNC, and EPC). Hub gene expression was validated using the UALCAN website. Clinical relevance was investigated using The Cancer Genome Atlas resources. Finally, we explored the association between metastasis-associated genes and immune infiltrates through the Tumor Immune Estimation Resource (TIMER) database and performed drug–gene interaction analysis using the Drug-Gene Interaction database. A total of 294 specific genes were related to melanoma metastasis and were mainly involved in the positive regulation of locomotion, mitotic cell cycle process, and epithelial cell differentiation. Four hub genes (CDK1, FOXM1, KIF11, and RFC4) were identified from the cytoHubba plugin of Cytoscape. CDK1 was significantly upregulated in metastatic melanoma compared with primary melanoma, and high CDK1 expression was positively correlated with worse overall survival. Immune infiltration analysis revealed that CDK1 expression negatively correlated with macrophage infiltration (Rho = − 0.164, P = 2.02e−03) and positively correlated with neutrophil cells (Rho = 0.269, P = 2.72e−07) in SKCM metastasis. In addition, we identified that CDK1 had a close interaction with 10 antitumor drugs. CDK1 was identified as a hub gene involved in the progression of melanoma metastasis and may be regarded as a therapeutic target for melanoma patients to improve prognosis and prevent metastasis in the future.