Abstract
The study of IRGPs to construct the prognostic signature in head and neck squamous cell carcinoma (HNSCC) has not yet elucidated. The objective of this study was to explore a novel model to predict the prognosis of HNSCC patients. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were set as training and validation cohorts, respectively. The least absolute shrinkage and selection operator (LASSO) and time-dependent ROC were employed to screen the highest frequency immune-related gene pairs (IRGPs) and their best cut-off value. Survival analysis, Cox regression analysis were applied to discover the effects of selected IRGPs signature on survival outcomes. The immune cell proportions were deconvoluted by the CIBERSORT method. After a couple of filtering, we obtained 22 highest frequency IRGPs. The overall survival time of HNSCC patients with a high score of IRGPs was shorter as compared to the ones with a low score in two independent datasets (P < 0.001). Six kinds of immune cells were found to be differentially distributed in the two different risk groups of HNSCC patients (P < 0.001). GO and GSEA analysis showed these differentially expressed genes enriched in multiple molecular functions. The new IRGPs signature probably confers a new insight into the prognosis prediction of HNSCC patients.