Abstract
Nicotinic acetylcholine receptors (CHRNs) expression and their critical role in various types of cancer have been reported. However, it is still unclear which CHRNs and their associated genes play essential roles in metastasis in melanoma patients. Here, we performed bioinformatics analyses on publicly available bulk RNA sequencing (RNA-seq) data of patients with melanoma to identify the CHRNs highly expressed in metastatic melanoma. We found that CHRNA1 was highly expressed in metastatic melanoma samples compared to primary melanoma samples and was strongly associated with CHRNB1 and CHRNG. These muscle-type CHRNs (CHRNA1, CHRNB1, and CHRNG) were correlated with the ZEB1 and Rho/ROCK pathway-related genes in metastatic melanoma samples. Pairwise correlations and enrichment analyses revealed that CHRNA1 was significantly associated with myogenesis/muscle contraction and cell cycle genes. Kaplan-Meier curves illustrated the involvement of CHRNA1, four of its correlated genes (DES, FLNC, CDK1, and CDC20), and the myogenesis gene signature in the prognosis of melanoma patients. Following the bulk RNA-seq analysis, single-cell RNA-seq (scRNA-seq) analysis showed that the CHRNA1-expressing melanoma cells are primarily metastatic and had high expression levels of CHRNB1, CHRNG, and myogenesis/cell cycle-related genes. Our bioinformatics analyses of the bulk RNA-seq and scRNA-seq data of patients with melanoma revealed that CHRNA1 and its correlated myogenesis/cell-related cycle genes are critical prognosis-related markers of metastatic melanoma.