Integrated analysis of single-cell, spatial and bulk RNA-sequencing identifies a cell-death signature for predicting the outcomes of head and neck cancer.

Background: Cell death plays an essential role in carcinogenesis, but its function in the recurrence and postoperative prognosis of head and neck cancer (HNC), which ranks as the 7 Methods: Data from five main subtypes of HNC related single-cell RNA sequencing (scRNA-seq) were recruited to establish a single-cell atlas, and the distribution of cell death models (CDMs) across different tissues as well as cell subtypes were analyzed. Bulk RNA-seq from the Cancer Genome Atlas Program (TCGA) dataset was subjected to a machine learning-based integrative procedure for constructing a consensus cell death-related signature risk score (CDRscore) model and validated by external data. The biofunctions including different expression analysis, immune cell infiltration, genomic mutations, enrichment analysis as well as cellchat analysis were compared between the high- and low- risk score groups categorized by this CDRscore model. Finally, samples from laryngeal squamous cell cancer (LSCC) were conducted by spatial transcriptomics (ST) to further validate the results of CDRscore model. Results: T cells from HNC patients manifested the highest levels of cell death while HPV infection attenuates malignant cell death based on single-cell atlas. CDMs are positively correlated with the tumor-cell stemness, immune-related score and T cells are infiltrated. A CDRscore model was established based on the transcription of ten cell death prognostic genes ( Conclusion: The CDRscore model could be utilized as a powerful prognostic indicator for HNC.