Identification of Genes Potentially Associated with Melanoma Tumorigenesis Through Co-Expression Network Analysis

Background: Melanoma is one of the most malignant and aggressive skin tumors, and its incidence is increasing worldwide. However, few studies have investigated the process of tumorigenesis from normal skin to melanoma. Methods: Several bioinformatics analyses, including GEO databases, Oncomine database, TCGA database, STRING, MCODE and cytoHubba plug-ins, GEPIA, TIMER and TRRUST and DGIdb, were performed to disclose the hub genes and immunology implicated in primary melanoma tumorigenesis. Finally, immunohistochemistry (IHC) and quantitative real-time PCR (qPCR) were used to validate the results of bioinformatics analysis in vitro. Results: A total of 295 overlapping DEGs (ODEGs) (157 upregulated and 138 downregulated) and 9 hub genes were identified between primary melanoma and normal skin tissues. Functional analysis of these 9 hub genes indicated that the genes were primarily enriched in cell chemotaxis, the chemokine-mediated signaling pathway, the extracellular region, the extracellular space, chemokine activity and CXCR3 chemokine receptor binding. KEGG pathway enrichment showed that these genes were primarily involved in the chemokine signaling pathway, cytokine-cytokine receptor interaction, the toll-like receptor signaling pathway, the cytosolic DNA sensing pathway and the TNF signaling pathway. Upregulated CCL5, CCL4, CXCL9 and CXCL10 demonstrated good overall survival (OS), and most of them have a higher expression in stage 0 and 1 of melanoma. Moreover, immune infiltration analysis showed that the above hub genes showed a strong positive correlation between their expression and infiltration of the six immune cell subsets. Transcription factor regulation network suggested that RELA and NFKB1 are the transcription factors of CCL4, CCL5, CXCL10 and CXCL2, while IRF7, IRF3 and IRF1 are the transcription factors of CCL5 and CXCL10. Drug-gene interaction analysis identified 46 drug-gene interactions. In vitro data demonstrated that the level of CCL4, CCL5, CXCL9 and CXCL10 is higher in melanoma than that in normal skin tissues, either at tissue or cell lines level. Conclusion: In summary, we identified 4 key chemokine members related to tumorigenesis and progression in primary melanoma, and these results may help to elucidate melanoma tumorigenesis and facilitate its treatment.