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

Melanoma is the most lethal type of skin cancer. Despite the breakthroughs in the clinical treatment of melanoma using tumor immunotherapy, many patients do not benefit from these immunotherapies because of multiple immunosuppressive mechanisms. Therefore, there is an urgent need to determine the mechanisms of tumor-immune system interactions and their molecular determinants to improve cancer immunotherapy. In this study, combined analysis of microarray data and single-cell RNA sequencing data revealed the key interactions between immune cells in the melanoma microenvironment. First, differentially expressed genes (DEGs) between normal and malignant tissues were obtained using GEO2R. The DEGs were then subjected to downstream analyses, including enrichment analysis and protein–protein interaction analysis, indicating that these genes were associated with the immune response of melanoma. Then, the GEPIA and TIMER databases were used to verify the differential expression and prognostic significance of hub genes, and the relationship between the hub genes and immune infiltration. In addition, we combined single cell analysis from GSE123139 to identify immune cell types, and validated the expression of the hub genes in these immune cells. Finally, cell-to-cell communication analysis of the proteins encoded by the hub genes and their interactions was performed using CellChat. We found that the CCL5-CCR1, SELPLG-SELL, CXCL10-CXCR3, and CXCL9-CXCR3 pathways might play important roles in the communication between the immune cells in tumor microenvironment. This discovery may reveal the communication basis of immune cells in the tumor microenvironment and provide a new idea for melanoma immunotherapy.