Identification of novel hub genes associated with lymph node metastasis of head and neck squamous cell carcinoma by completive bioinformatics analysis

BACKGROUND: Head and neck squamous cell carcinoma (HNSCC) is one of the most serious diseases affecting populations worldwide and lymph node metastasis is a key pathological feature of HNSCC which predicts poor survival. However, the molecular mechanisms associated with the development of lymph node metastasis in HNSCC have not been fully elucidated. METHODS: Differentially expressed genes (DEGs) were identified in two HNSCC datasets (GES6631 and GES58911). Functional annotation analysis was constructed via Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Meanwhile, the protein-protein interaction (PPI) network and module analysis using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape were carried out to identify the hub genes. The expression differences, overall survival (OS), and disease-free survival (DFS) of hub genes were analyzed by Gene Expression Profiling Interactive Analysis 2 (GEPIA2) and verified by immunohistochemistry (IHC) from Human Protein Atlas (HPA). Moreover, receiver operating characteristic (ROC) curve analysis was conducted to judge whether those hub genes had good diagnostic and prognostic ability, and the web tool Tumor Immune Estimation Resource (TIMER) was used to analyze the correlation of hub genes’ expression and immune infiltration. RESULTS: In total, 913 DEGs including 476 upregulated and 437 downregulated genes were identified. The genes Aurora kinase A (AURKA), CyclinB1 (CCNB1), Cyclin-dependent kinase regulatory subunit 1B (CKS1B), Serpin Family H Member 1(SERPINH1), and Transforming growth factor-beta-induced protein (TGFBI) were screened out as hub genes and were associated with lymph node metastasis, showing notably abnormal expression in HNSCC samples, and the high expression of all the hub genes in HNSCC patients was related to worse overall survival. CONCLUSIONS: The genes AURKA, CCNB1, CKS1B, SERPINH1, and TGFBI may be involved in the lymph node metastasis of HNSCC and reveal the potential to serve as molecular biomarkers in the diagnosis of HNSCC. This study may help to elucidate the molecular mechanisms of the development of lymph node metastasis and facilitate the selection of targets for the treatment and diagnosis of HNSCC.