Identification of Key Genes in Nasopharyngeal Carcinoma Based on Bioinformatics Analysis

This study aimed to identify key genes associated with the pathogenesis of nasopharyngeal carcinoma (NPC) by bioinformatics analysis.

Datasets (GSE13597 and GSE34573) were screened and downloaded from the comprehensive gene expression database (GEO). GEO2R online tool was adopted to analyze microarray data GSE13597 and GSE34573 related to NPC. Volcano plot was generated using Bioconductor in R software. "Pheatmap" was used to draw heatmaps based on the top 10 regulated genes of GSE13597 and GSE34573. GO and KEGG analyses were conducted via online tool DAVID. We uploaded the DEGs of NPC to STRING software and then used Cytoscape software to draw PPI network of DEGs.

216 DEGs were obtained in GSE13597 between patient and control group (111 up-regulated DEGs and 105 down-regulated DEGs). 1101 DEGs were obtained in GSE34573 (470 up-regulated DEGs and 641 down-regulated DEGs). 63 common differential genes were screened named co-DEGs in the two datasets. These DEGs were mainly associated with defense response to bacterium, cell-matrix adhesion, chemokine-mediated signaling pathway, tissue homeostasis, humoral immune response, cilium movement, cilium organization, cilium assembly, and epithelial cilium movement. KEGG pathway enrichment analysis showed that DEGs were mainly involved in viral protein interaction with cytokine and cytokine receptor, salivary secretion, p53 signaling pathway, IL-17 signaling pathway, cell cycle, PI3K-Akt signaling pathway, and ECM-receptor interaction. We identified seven hub genes, including FN1, MMP-10, MUC1, KIF23, CDK1, MUC5B, and MUC5AC.

Seven hub genes, including FN1, MMP-10, MUC1, KIF23, CDK1, MUC5B, and MUC5AC, might be therapeutic potential biomarkers of NPC.