Abstract
Background: Head and neck squamous cell carcinoma (HNSCC) is a common malignancy with high mortality and morbidity worldwide, but the underlying biological mechanisms of molecules and tumor infiltrating-immune cells (TIICs) are still unknown. Methods and Results: We obtained mRNAs, lncRNAs, and miRNAs expression profiles of 546 HNSCC from The Cancer Genome Atlas (TCGA) database to develop a ceRNA network. CIBERSORT was employed to estimate the fraction of 22 types of TIICs in HNSCC. Univariate and multivariate Cox regression and lasso regression analyses were used to develop prognostic signatures. Then, two novel risk signatures were constructed respectively based on six ceRNAs (ANLN, KIT, PRKAA2, NFIA, PTX3 and has-miR-148a-3p) and three immune cells (naïve B cells, regulatory T cells and Neutrophils). Kaplan-Meier (K-M) analysis and Cox regression analysis further proved that these two signatures were significant prognostic factors independent of multiple clinicopathological characteristics. Two nomograms were built based on ceRNAs-riskScore and TIICs-riskScore that could be used to predict the prognosis of HNSCC. Co-expression analysis showed significant correlations between miR-148a-3p and naive B cells, naive B cells and plasmas cells. Conclusion: Through construction of the ceRNA network and estimation of TIICs, we established two risk signatures and their nomograms with excellent utility, which indicated the potential molecular and cellular mechanisms, and predicted the prognosis of HNSCC.