Fig. 6

High-dimension Weighted Correlation Network Analysis (hdWGCNA) of melanoma cell reveals the crucial roles of C1. A Dendrogram demonstrating nine gene modules in the scale-free network. B Nine gene modules and the associated top 10 hub gene were displayed according to the hdWGCNA pipeline. C Correlation between nine gene modules. Green represents negative correlation and purple represents positive correlation. D Nine module levels in different melanoma cell clusters. E Violin plot demonstrating that Module 1 (M1) activity is significantly higher in Cluster 1 than in the other melanoma cell clusters. F Protein–protein interaction (PPI) network for the top 25 genes in M1. G GO and KEGG enrichmnt analysis of M1. H Expression level of M1 hub genes between melanoma tissues in the TCGA-SKCM cohort (n = 469) and control tissues in in both the TCGA and GTEx databases (n = 813). Wilcoxon test, ***P < 0.05. I Forest plot showing the results of Cox regression analysis of the mean survival (overall survival) of the 10 M1 hub genes in the TCGA-SKCM cohort. J–K All cell types were analyzed by Cellchat and melanoma cells were classified according to stemness level as cluster 1 (high), cluster 0 (low), and other melanoma cells (medium). The results showed both interaction numbers and interaction strengths. L–O Cellchat analysis showed weaker communication strength between melanoma cell cluster 1 and T cells compared to melanoma cell cluster 0