Abstract
It has been found that progression from leukoplakia to head and neck squamous cell carcinoma (HNSCC) is a long-term process that may involve changes in the multicellular ecosystem. We acquired scRNA-seq samples information from gene expression omnibus and UCSC Xena database. The BEAM function was used to construct the pseudotime trajectory and analyze the differentially expressed genes in different branches. We used the ssGSEA method to explore the correlation between each cell subgroup and survival time, and obtained the cell subgroup related to prognosis. During the progression from leukoplakia to HNSCC, we found several prognostic cell subgroups, such as AURKB + epithelial cells, SFRP1 + fibroblasts, SLC7A8 + macrophages, FCER1A + CD1C + dendritic cells, and TRGC2 + NK/T cells. All cell subgroups had two different fates, one tending to cell proliferation, migration, and enhancement of angiogenesis capacity, and the other tending to inflammatory immune response, leukocyte chemotaxis, and T cell activation. Tumor-promoting genes such as CD163 and CD209 were highly expressed in the myeloid cells, and depletion marker genes such as TIGIT, LAG3 were highly expressed in NK/T cells. Our study may provide a reference for the molecular mechanism of HNSCC and theoretical basis for the development of new therapeutic strategies.