Identification and validation of a novel signature based on cell-cell communication in head and neck squamous cell carcinoma by integrated analysis of single-cell transcriptome and bulk RNA-sequencing

BACKGROUND: The heterogeneous crosstalk between tumor cells and other cells in their microenvironment means a notable difference in clinical outcomes of head and neck squamous cell carcinoma (HNSCC). CD8+ T cells and macrophages are effector factors of the immune system, which have direct killing and phagocytosis effects on tumor cells. How the evolution of their role in the tumor microenvironment influences patients clinically remains a mystery. This study aims to investigate the complex communication networks in the HNSCC tumor immune microenvironment, elucidate the interactions between immune cells and tumors, and establish prognostic risk model. METHODS: 20 HNSCC samples single-cell rna sequencing (scRNA-seq) data and bulk rna-seq data were derived from public databases. The “cellchat” R package was used to identify cell-to-cell communication networks and prognostic related genes, and then cell-cell communication (ccc) molecular subtypes were constructed by unsupervised clustering. Kaplan-Meier(K-M) survival analysis, clinical characteristics analysis, immune microenvironment analysis, immune cell infiltration analysis and CD8+T cell differentiation correlation analysis were performed. Finally, the ccc gene signature including APP, ALCAM, IL6, IL10 and CD6 was constructed based on univariate Cox analysis and multivariate Cox regression. Kaplan-Meier analysis and time-dependent receiver operating characteristic (ROC) analysis were used to evaluate the model in the train group and the validation group, respectively. RESULTS: With CD8+T cells from naive to exhaustion state, significantly decreased expression of protective factor (CD6 gene) is associated with poorer prognosis in patients with HNSCC. The role of macrophages in the tumor microenvironment has been identified as tumor-associated macrophage (TAM), which can promote tumor proliferation and help tumor cells provide more nutrients and channels to facilitate tumor cell invasion and metastasis. In addition, based on the strength of all ccc in the tumor microenvironment, we identified five prognostic ccc gene signatures (cccgs), which were identified as independent prognostic factors by univariate and multivariate analysis. The predictive power of cccgs was well demonstrated in different clinical groups in train and test cohorts. CONCLUSION: Our study highlights the propensity for crosstalk between tumors and other cells and developed a novel signature on the basis of a strong association gene for cell communication that has a powerful ability to predict prognosis and immunotherapy response in patients with HNSCC. This may provide some guidance for developing diagnostic biomarkers for risk stratification and therapeutic targets for new therapeutic strategies.