Identification of immune-related genes contributing to head and neck squamous cell carcinoma development using weighted gene co-expression network analysis.

Abstract

This study aimed to identify genes related to the degree of immune cell infiltration in head and neck squamous cell carcinoma (HNSCC), explore their new biological functions, and evaluate their diagnostic and prognostic value in HNSCC. Transcriptomic data from The Cancer Genome Atlas (TCGA) HNSCC dataset was used to screen differentially expressed genes between tumors and normal tissues, followed by weighted correlation network analysis (WGCNA) to identify immune-related modules. Differential gene expression, immune cell infiltration, and survival analyses were performed to screen key genes. The expression of these key genes was validated in Oncomine and gene expression omnibus (GEO) datasets and by immunohistochemistry (IHC). 1869 and 1578 genes were significantly upregulated and downregulated in HNSCC. WGCNA showed that the brown module was associated with the most significant number of immune-related genes. PPI network analysis demonstrated that PPL, SCEL, KRT4, KRT24, KRT78, KRT13, SPRR3, TGM3, CRCT1, and CRNN were key components in the brown module. Furthermore, the expression levels of KRT4, KRT78, KRT13, and SPRR3 in HNSCC correlated with infiltration levels of CD8+ T cells and macrophages. Survival analyses revealed that the expression of KRT78, KRT13, and SPRR3 in HNSCC correlated with overall survival (OS). The IHC assay indicated that KRT13 (p=.042), KRT78 (p < .001), and SPRR3 (p=.022) protein expression levels in HNSCC were significantly lower than in normal tissues. Analysis of GSE65858 and GSE41613 datasets showed that a worse OS was associated with low expression of KRT78 (p=.0086, and p=.005) and SPRR3 (p=.017, and p=.02). Our findings suggest that KRT4, KRT78, KRT13, and SPRR3 are related to the occurrence and development of HNSCC. Importantly, KRT78 and SPRR3 might serve as diagnostic and prognostic biomarkers of HNSCC.