BackgroundAdvances in cancer research have allowed for early diagnosis and improved treatment of cutaneous melanoma (CM). However, its invasiveness and recurrent metastasis, along with rising resistance to newer therapies, have lent urgency to the search for novel biomarkers and the underlying molecular mechanisms of CM.MethodsSingle nucleotide polymorphism- (SNP-) related genes were obtained from the sequencing data of 428 CM samples in The Cancer Genome Atlas. Functional enrichment of these genes was analysed in clusterProfiler. Additionally, a protein-protein interaction (PPI) network was constructed with the Search Tool for the Retrieval of Interacting Gene (STRING) database. Gene Expression Profiling Interactive Analysis (GEPIA) was used to identify the expression and prognostic value of mutated genes. Finally, the Tumour Immune Estimation Resource (TIMER) analysed the relationship between gene expression and immune cell infiltration.ResultsWe constructed a PPI network from the top 60 SNP-related genes. Mutated genes were mainly involved in calcium and oxytocin signalling pathways, as well as circadian entrainment. In addition, three SNP-related genes,BRAF,FLG, andSORL1, were significantly associated with patient prognosis.BRAFandSORL1were positively associated with infiltration abundance of B cells, CD8+ T cells, CD4+ T cells, neutrophils, and dendritic cells, whereasFLGexpression was negatively associated. Furthermore, higher immune cell infiltration was positively correlated with good prognosis.ConclusionsOur study provides vital bioinformatic data and a relevant theoretical basis to further explore the molecular pathogenesis of CM and improve patient prognosis.