Abstract
There has been an increase in the mortality rate of thyroid cancer (THCA), which is the most common endocrine malignancy. We identified six distinct cell types in the THAC microenvironment by analyzing single-cell RNA sequencing (Sc-RNAseq) data from 23 THCA tumor samples, indicating high intratumoral heterogeneity. Through re-dimensional clustering of immune subset cells, myeloid cells, cancer-associated fibroblasts, and thyroid cell subsets, we deeply reveal differences in the tumor microenvironment of thyroid cancer. Through an in-depth analysis of thyroid cell subsets, we identified the process of thyroid cell deterioration (normal, intermediate, malignant cells). Through cell-to-cell communication analysis, we found a strong link between thyroid cells and fibroblasts and B cells in the MIF signaling pathway. In addition, we found a strong correlation between thyroid cells and B cells, TampNK cells, and bone marrow cells. Finally, we developed a prognostic model based on differentially expressed genes in thyroid cells from single-cell analysis. Both in the training set and the testing set, it can effectively predict the survival of thyroid patients. In addition, we identified significant differences in the composition of immune cell subsets between high-risk and low-risk patients, which may be responsible for their different prognosis. Through in vitro experiments, we identify that knockdown of NPC2 can significantly promote thyroid cancer cell apoptosis, and NPC2 may be a potential therapeutic target for thyroid cancer. In this study, we developed a well-performing prognostic model based on Sc-RNAseq data, revealing the cellular microenvironment and tumor heterogeneity of thyroid cancer. This will help to provide more accurate personalized treatment for patients in clinical diagnosis.