Identification novel prognostic signatures for Head and Neck Squamous Cell Carcinoma based on ceRNA network construction and immune infiltration analysis.

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.

Zhou H ，He Y ，Li L ，Wu C ，Hu G ... -  《International Journal of Medical Sciences》

Identification of Novel Biomarkers for Predicting Prognosis and Immunotherapy Response in Head and Neck Squamous Cell Carcinoma Based on ceRNA Network and Immune Infiltration Analysis.

Patients with head and neck squamous cell carcinoma (HNSCC) have poor prognosis and show poor responses to immune checkpoint (IC) inhibitor (ICI) therapy. Competing endogenous RNA (ceRNA) networks, tumor-infiltrating immune cells (TIICs), and ICIs may influence tumor prognosis and response rates to ICI therapy. This study is aimed at identifying prognostic and IC-related biomarkers and key TIIC signatures to improve prognosis and ICI therapy response in HNSCC patients. Ninety-five long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and 1746 mRNAs were identified using three independent methods. We constructed a ceRNA network and estimated the proportions of 22 immune cell subtypes. Ten ceRNAs were related to prognosis according to Kaplan-Meier analysis. Two risk signatures based, respectively, on nine ceRNAs (ANLN, CFL2, ITGA5, KDELC1, KIF23, NFIA, PTX3, RELT, and TMC7) and three immune cell types (naïve B cells, neutrophils, and regulatory T cells) via univariate Cox regression, least absolute shrinkage and selection operator, and multivariate Cox regression analyses could accurately and independently predict the prognosis of HNSCC patients. Key mRNAs in the ceRNA network were significantly correlated with naïve B cells and regulatory T cells and with stage, grade, and immune and molecular subtype. Eight IC genes exhibited higher expression in tumor tissues and were correlated with eight key mRNAs in the ceRNA network in HNSCC patients with different HPV statuses according to coexpression and TIMER 2.0 analyses. Most drugs were effective in association with expression of these key signatures (ANLN, CFL2, ITGA5, KIF23, NFIA, PTX3, RELT, and TMC7) based on GSCALite analysis. The prognostic value of key biomarkers and associations between key ceRNAs and IC genes were validated using online databases. Eight key ceRNAs were confirmed to predict response to ICI in other cancers based on TIDE analysis. We constructed two risk signatures to accurately predict prognosis in HNSCC. Key IC-related signatures may be associated with response to ICI therapy. Combinations of ICIs with inhibitors of eight key mRNAs may improve survival outcomes of HNSCC patients.

Guo Y ，Pan WK ，Wang ZW ，Su WH ，Xu K ，Jia H ，Chen J ... -  《-》

Excavating novel diagnostic and prognostic long non-coding RNAs (lncRNAs) for head and neck squamous cell carcinoma: an integrated bioinformatics analysis of competing endogenous RNAs (ceRNAs) and gene co-expression networks.

Yang L ，Lu P ，Yang X ，Li K ，Chen X ，Qu S ... -  《-》

Tumor Mutation Burden, Immune Cell Infiltration, and Construction of Immune-Related Genes Prognostic Model in Head and Neck Cancer.

Background: Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide, and the prognosis of HNSCC remains bleak. Numerous studies revealed that the tumor mutation burden (TMB) could predict the survival outcomes of a variety of tumors. Objectives: This study aimed to investigate the TMB and immune cell infiltration in these patients and construct an immune-related genes (IRGs) prognostic model. Methods: The expression data of 546 HNSCC patients were obtained from The Cancer Genome Atlas (TCGA) database. All patients were divided into high- and low- TMB groups, and the relationship between TMB and clinical relevance was further analyzed. The differentially expressed genes (DEGs) were identified using the R software package, limma. Functional enrichment analyses were conducted to identify the significantly enriched pathways between two groups. CIBERSORT algorithm was adopted to calculate the abundance of 22 leukocyte subtypes. The IRGs prognostic model was constructed via the multivariate Cox regression analysis. Results: Missense mutation and single nucleotide variants (SNV) were the most predominant mutation types in HNSCC. TP53, TTN, and FAT1 were the most frequently mutated genes. Patients with high TMB were observed with worse survival outcomes. The functional analysis of TMB associated DEGs showed that the identified DEGs mainly involved in spliceosome, RNA degradation, proteasome, and RNA polymerase pathways. We observed that macrophages, T cells CD8, and T cells CD4 memory were the most commonly infiltrated subtypes of immune cells in HNSCC. Finally, an IRGs prognostic model was constructed, and the AUC of the ROC curve was 0.635. Conclusions: Our results suggest that high TMB is associated with poor prognosis in HNSCC patients. The constructed model has potential prognostic value for the prognosis of these individuals, and it needs to be further validated in large-scale and prospective studies.

Jiang AM ，Ren MD ，Liu N ，Gao H ，Wang JJ ，Zheng XQ ，Fu X ，Liang X ，Ruan ZP ，Tian T ，Yao Y ... -  《International Journal of Medical Sciences》

Development and validation of a novel ferroptosis-related gene signature for predicting prognosis and immune microenvironment in head and neck squamous cell carcinoma.

Ferroptosis plays an important role across variable cancer types. However, few studies have focused on the prognostic patterns of ferroptosis-related genes in HNSCC. Cohorts with mRNA expression profiles, as well as corresponding clinical data of HNSCC patients from published studies, were collected and consolidated from public databases. We performed random survival forest analysis, Kaplan-Meier (KM) analysis of best combinations, and Cox regression analysis on 231 ferroptosis-related genes to construct a gene signature in the discovery cohort (TCGA), and later validated it in the validation cohort (GEO). The 7-gene signature was constructed to stratify patients into two groups according to their level of risk. Poorer overall survival (OS) was detected in the high risk (HRisk) group than in the low risk (LRisk) group in both the TCGA cohort (P < 0.0001, HR = 1.71, 95%CI:1.41-2.07) and the GEO cohort (P < 0.001, HR = 1.68, 95%CI:1.32-2.13). The risk score was identified as an independent predictive factor of OS in multivariate Cox regression analyses (HR > 1, P < 0.0001) in both cohorts. The signature's predictive capacity was proven by the time-dependent receiver operating characteristic (ROC) curve analysis and nomogram analysis. Functional enrichment analysis revealed that immunosuppressive pathways such as matrix extracellular space, and (transforming growth factor-β)TGF-β were enriched. The HRisk group was strongly associated with upregulation of both cancer-related pathways and stromal scores, while higher proportions of anti-tumor immune cells and immune signatures were enriched in the LRisk group. In conclusion, the signature based on 7 ferroptosis-related genes could be applicable for predicting the prognosis of HNSCC, indicating that ferroptosis may be a potential therapeutic target for HNSCC.

He F ，Chen Z ，Deng W ，Zhan T ，Huang X ，Zheng Y ，Yang H ... -  《-》