A signature of immune-related gene pairs (IRGPs) for risk stratification and prognosis of oral cancer patients

Background With low response to present immunotherapy, it is imperative to identify new immune-related biomarkers for more effective immunotherapies for oral cancer. Methods RNA profiles for 390 oral cancer patients and 32 normal samples were downloaded from The Cancer Genome Atlas (TCGA) database and differentially expressed genes (DEGs) were analyzed. Immune genesets from ImmPort repository were overlapped with DEGs. After implementing univariate Cox analysis and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis, key immune-related gene pairs (IRGPs) among the overlapped DEGs for predicting the survival risk were obtained. Then, the cutoff of risk score was calculated by the receiver operating characteristic (ROC) curve to stratify oral cancer patients into high and low-risk groups. Multivariate Cox analysis was used to analyze independent prognostic indicators for oral cancer. Besides, infiltration of immune cells, functional annotation, and mutation analysis of IRGPs were conducted. Biological functions correlated with IRGPs were enriched by Gene Set Enrichment Analysis (GSEA) method. Results We identified 698 differentially expressed genes (DEGs) in response to oral cancer. 17 IRGPs among the DEGs were identified and integrated into a risk score model. Patients in the high-risk group have a significantly worse prognosis than those in the low-risk group in both training (P