Identification of CXCL13 as an Immune-Related Biomarker Associated with Tumorigenesis and Prognosis in Cutaneous Melanoma Patients

BACKGROUND: Melanoma is one of the most lethal tumors and its treatment is still challenging. It is urgent to detect novel therapy targets in melanoma. MATERIAL/METHODS: The GEO dataset was used to obtain a list of DEGS (differentially-expressed genes). Integrative bioinformatics analyses, including HPRD database, TCGA data, and TIMER, were performed to determine the role of CXCL13 in SKCM (skin cutaneous melanoma) progression and the immune environment. Furthermore, Pearson correlation coefficient analysis was used to measure correlations between CXCL13 and its co-expressed genes. Survival analysis, GO, and KEGG enrichment analysis were performed to investigate the role of CXCL13 in SKCM. RESULTS: A total of 41 DEGs were identified in 3 GEO datasets, and 4 out of 41 DEGs are hub genes. Among the 4 hub genes, CXCL13 is involved in the most KEGG terms. CXCL13 is co-expressed with well-known immune checkpoint blockade targets, and it was associated with better overall survival. In addition, CXCL13 levels in infiltrating immune cells (neutrophil and myeloid dendritic cells) affect prognosis and survival in SKCM. Functional enrichment analysis clarified that CXCL13-co-expressed top 30 genes were associated with immune signaling pathways. Network analysis identified CXCL13 as a hub gene that interacts with CXCR5 to participate in immune-related biological process. CONCLUSIONS: This study found that CXCL13 is associated with SKCM tumorigenesis and prognosis and immune infiltrations. Our result suggests that CXCL13 has great potential in development of novel immunotherapy targets in melanoma.