ABSTRACT: BACKGROUND:The immune response within the tumor microenvironment plays a key role in tumorigenesis and determines the clinical outcomes of head and neck squamous cell carcinoma (HNSCC). However, to date, a paucity of robust, reliable immune-related biomarkers has been identified that are capable of estimating prognosis in HNSCC patients. METHODS:High-throughput RNA sequencing was performed in tumors and matched adjacent tissues from five HNSCC patients, and the immune signatures expression of 730 immune-related transcripts selected from the nCounter PanCancer Immune Profiling Panel were assessed. Survival analyzes were performed in a training cohort, consisting of 416 HNSCC cases, retrieved from The Cancer Genome Atlas (TCGA) database. A prognostic signature was built, using elastic net-penalized Cox regression and backward, stepwise Cox regression analyzes. The outcomes were validated by an independent cohort of 115 HNSCC patients, using tissue microarrays and immunohistochemistry staining. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) was also used to estimate the relative fractions of 22 immune-cell types and their correlations coefficients with prognostic biomarkers. RESULTS:Collectively, 248 immune-related genes were differentially expressed in paired tumors and normal tissues using RNA sequencing. After process screening in the training TCGA cohort, four immune-related genes (PVR, TNFRSF12A, IL21R, and SOCS1) were significantly associated with overall survival (OS). Integrating these genes with Path_N stage, a multiplex model was built and suggested better performance in determining 5?years OS (receiver operating characteristic (ROC) analysis, area under the curve (AUC)=0.709) than others. Further protein-based validation was conducted in 115 HNSCC patients. Similarly, high expression of PVR and TNFRSF12A were associated with poor OS (Kaplan-Meier p=0.017 and 0.0032), while high expression of IL21R and SOCS1 indicated favorable OS (Kaplan-Meier p<0.0001 and =0.0018). The integrated model with Path_N stage still demonstrated efficacy in OS evaluation (Kaplan-Meier p<0.0001, ROC AUC=0.893). Besides, the four prognostic genes were significantly correlated with activated CD8+ T cells, CD4+ T cells, follicular helper T cells and regulatory T cells, implying the possible involvement of these genes in the immunoregulation and development of HNSCC. CONCLUSIONS:The well-established model encompassing both immune-related biomarkers and clinicopathological factor might serve as a promising tool for the prognostic prediction of HNSCC.
