Identification of new biomarkers and pathways associated with treatment and prognosis in melanoma patients
Abstract Background Skin cutaneous melanoma (SKCM) is known as the most malignancy and treatment-resistant in human tumor, causing about 72% of deaths in skin carcinoma. However, the potential mechanism and new effective targets remain to be further elucidated. Available datasets such as Gene Expression Omnibus (GEO) can be utilized to search for novel therapeutic targets and prognostic biomarkers. Methods Three data sets were downloaded from GEO database . The differentially expressed genes (DEGs) were identified via Venn software. Protein‐protein interaction network of DEGs was developed and the module hub genes analysis was constructed by Cytoscape. Subsequently, multiple online tools and Kaplan-Meier survival curves were analyzed to detect underlying signaling pathways, gene expression, drug-gene interaction and prognostic value of hub genes. In addition, we explored the correlation between hub genes and immune cell infiltration. At last, the related miRNA, lncRNA networks were constructed by R software. Results A total of 308 DEGs and 12 hub genes were identified. Function and pathway enrichment results demonstrated a correlation between DEGs and the tumor microenvironment, immune response and melanoma tumorigenesis. Subsequently, we focused on assessing potential value of 12 hub genes. Seven hub genes ( CCL4, CCL5, NMU, GAL, CXCL9, CXCL10, CXCL13 ) were identified with significant overall survival for prognosis. What’s more, five of these seven hub genes were found to be related to clinical stages (P values<0.05). In addition, the most important pathways of hub genes include interleukin-10 signaling, peptide ligand-binding receptors, which play important roles in tumor microenvironment for immune activation or immunosuppressive by regulating the infiltration of immune cells. Our results revealed a strong positive correlation between gene expression (CCL4, CCL5, CXCL9, CXCL10 and CXCL13) and immune cell infiltration (B-cell, CD8+ T cells, CD4+ T cells, macrophages, Neutrophils, Dendritic cells). Interestingly, 8 of 12 hub genes (CXCL10, CCL4, CCL5, IL6, CXCL2, PTGER3, GAL, NPY1R) were also found in the predicted drug-gene interaction. The related miRNA, lncRNA for diagnosis and prognosis were found in networks. Conclusion In conclusion, CCL4, CCL5, NMU, GAL, CXCL9, CXCL10, CXCL13 were of high prognostic value and may be potential targets for the diagnosis and therapy of patients with melanoma.