A prognostic model based on Scissor plus cancer associated fibroblasts identified from bulk and single cell RNA sequencing data in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is one of the most lethal diseases in the world, which often recur after multimodality treatment approaches, leading to a poor prognosis. Fibroblasts, a heterogeneous component of the tumor microenvironment, can modulate numerous aspects of tumor biology and have been increasingly acknowledged in dictating the clinical outcome of patients with HNSCC. However, the subpopulation of fibroblasts that are related to the prognosis of HNSCC has not yet been fully explored. To do so, we combined a single-cell RNA sequencing (scRNA-seq) dataset and bulk RNA-sequencing dataset with clinical information, identifying the fibroblast population that are related to poor prognosis of HNSCC. We found these specific population of fibroblasts are less differentiated. In addition, to identify the prognostic signatures of HNSCC, bioinformatics analysis included least absolute shrinkage and selection operator (LASSO) analyses and univariate cox and were performed. We selected 12 prognosis-related genes for constructing a risk model using The Cancer Genome Atlas (TCGA). The AUC values and calibration plots of this model indicated good prognostic prediction efficacy. This model also was validated in two Gene Expression Omnibus (GEO) datasets. In conclusion, we constructed an optimal model that was derived from single cell RNA-seq and bulk RNA-seq to predict the survival probability of HNSCC patients. Among this model, AKR1C3 higher expression in cancer associated fibroblasts (CAFs) of HNSCC has been confirmed by preliminary experiments.