Immune Infiltration-Related ceRNA Network Revealing Potential Biomarkers for Prognosis of Head and Neck Squamous Cell Carcinoma

Abstract

Background. Head and neck squamous cell carcinoma (HNSCC) is a frequently lethal malignancy, and the mortality is considerably high. The tumor microenvironment (TME) has been identified as a critical participation in cancer development, treatment, and prognosis. However, competing endogenous RNA (ceRNA) networks grouping with immune/stromal scores of HNSCC patients need to be further illustrated. Therefore, our study aimed to provide clues for searching promising prognostic markers of TME in HNSCC. Materials and Methods. ESTIMATE algorithm was used to calculate immune scores and stromal scores of the enrolled HNSCC patients. Differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified by comparing the expression difference between high and low immune/stromal scores. Then, a ceRNA network and protein-protein interaction (PPI) network were constructed for selecting hub regulators. In addition, survival analysis was performed to access the association between immune scores, stromal scores, and differentially expressed RNAs in the ceRNA network and the overall survival (OS) of HNSCC patients. Then, the GSE65858 datasets from Gene Expression Omnibus (GEO) database was used for verification. At last, the difference between the clinical characteristics and immune cell infiltration in different expression groups of IL10RA, PRF1, and IL2RA was analyzed. Results. Survival analysis showed a better OS in the high immune score group, and then we constructed a ceRNA network composed of 97 DEGs, 79 DELs and 22 DEMs. Within the ceRNA network, FOXP3, IL10RA, STAT5A, PRF1, IL2RA, miR-148a-3p, miR-3065-3p, and lncRNAs, including CXCR2P1, HNRNPA1P21, CTA-384D8.36, and IGHV1OR15-2, were closely correlated with the OS of HNSCC patients. Especially, using the data from GSE65858, we successfully verified that IL10RA, PRF1, and IL2RA were not only significantly upregulated in patients high immune scores, but also their high expressions were associated with longer survival time. In addition, stratified analysis showed that PRF1 and IL2RA might be involved in the mechanism of tumor progress. Conclusion. In conclusion, we constructed a ceRNA network related to the TME of HNSCC, which provides candidates for therapeutic intervention and prognosis evaluation.