Bioinformatics analysis to identify the critical genes, microRNAs and long noncoding RNAs in melanoma
Melanoma, which is usually induced by ultraviolet light exposure and the following DNA damage, is the most dangerous skin cancer. The purpose of the present study was to screen key molecules involved in melanoma.
Microarray data of E-MTAB-1862 were downloaded from the ArrayExpress database, which included 21 primary melanoma samples and 11 benign nevus samples. In addition, the RNASeq version 2 and microRNA (miRNA) sequencing data of cutaneous melanoma were downloaded from The Cancer Genome Atlas database. After identifying the differentially expressed genes (DEGs) using Limma package, enrichment analysis and protein-protein interaction (PPI) network analysis were performed separately for them using DAVID software and Cytoscape software. In addition, survival analysis and regulatory network analysis were further performed by log-rank test and Cytoscape software, respectively. Moreover, real-time reverse transcription polymerase chain reaction (RT-PCR) was performed to further verify the expression patterns of several selected DEGs.
A total of 382 DEGs were identified in primary melanoma samples, including 206 upregulated genes and 176 downregulated genes. Functional enrichment analysis showed that COL17A1 was enriched in epidermis development. In the PPI network, CXCL8 (degree = 29) and STAT1 (degree = 28) had higher degrees and could interact with each other. Survival analysis showed that 21 DEGs, 55 long noncoding RNAs (lncRNAs) and 32 miRNAs were found to be associated with prognosis. Furthermore, several regulatory relationships were found in the lncRNA-gene regulatory network (such as RP11-361L15.4 targeting COL17A1) and the miRNA-gene regulatory network (such as hsa-miR-375 targeting CCL27 and hsa-miR-375 targeting insulin-like growth factor 1 receptor [IGF1R]). Real-time RT-PCR results showed that the overall direction of differential expression was consistent except COL17A1.
CXCL8 interacted with STAT1, CCL27, and IGF1R targeted by hsa-miR-375, and COL17A1 targeted by RP11-361L15.4 might function in the development and progression of melanoma, which should be verified by more detailed experiments.