Abstract

Background: Head and neck squamous cell carcinoma (HNSCC) was the seventh most common cancer worldwide in 2018. Lymphatic metastasis (LM) is closely related to HNSCC prognosis and recurrence. However, the underlying mechanism of LM remains unclear. Therefore, this study aimed to identify the key genes in the LM of HNSCC. Methods: We used The Cancer Genome Atlas (TCGA) to identify differentially expressed genes (DEGs) between LM and non-LM cases. A random forest model, the Search Tool for the Retrieval of Interacting Genes, Cytoscape, and cytoHubba were used to identify hub genes among DEGs, including KRT20 (Cytokeratins 20). We analyzed the survival of KRT20 in TCGA, and we overexpressed KRT20 in HNSCC cell lines to investigate its effects on migration and invasion. We also correlated the expression of KRT20 in HNSCC tissue microarrays with survival and clinicopathological features. Results: We identified 243 DEGs—143 upregulated genes and 100 downregulated genes. Further analysis revealed that KRT20 is a potential key gene associated with LM and overall survival rates among patients with HNSCC. Overexpression of KRT20 increased the migration and invasion ability of HNSCC cell lines Tu686 and FD-LSC-1. Tissue microarray studies demonstrated an overexpression of KRT20 among N1+ patients (including N1-N3 patients). Survival analysis results and the clinicopathological features of HNSCC tissue microarrays were consistent with our analysis of TCGA. Thus, a high KRT20 expression level might suggest an adverse HNSCC prognosis. Our gene set enrichment analysis showed that KRT20 participates in many metabolic pathways, including those related to tumorigenesis and cancer development. Conclusions: We propose that KRT20 may be a key gene in HNSCC with LM.