Single-Cell RNA Sequencing Revealed a 3-Gene Panel Predicted the Diagnosis and Prognosis of Thyroid Papillary Carcinoma and Associated With Tumor Immune Microenvironment

ObjectiveThe objective of this research was to screen prognostic related genes of thyroid papillary carcinoma (PTC) by single-cell RNA sequencing (scRNA-seq), to construct the diagnostic and prognostic models based on The Cancer Genome Atlas Thyroid Cancer (TCGA-THCA) data, and to evaluate the association between tumor immune microenvironment and the prognostic model. MethodThe differentially expressed genes (DEGs) and tumor evolution were analyzed by scRNA-seq based on public databases. The potential regulatory networks of DEGs related to prognosis were analyzed by multi-omics data in the THCA. Logistic regression and Cox proportional hazards regression were utilized to construct the diagnosis and prognostic model of PTC. The performance of the diagnostic model was verified by bulk RNA sequencing (RNA-seq) of our cohort. The tumor immune microenvironment associated with the prognostic model was evaluated using multi-omics data. In addition, qRT-PCR was performed on tumor tissues and adjacent normal tissues of 20 patients to verify the expression levels of DEGs. ResultsThe DEGs screened by scRNA-seq can distinguish between tumor and healthy samples. DEGs play different roles in the evolution from normal epithelial cells to malignant cells. Three DEGs ((FN1, CLU, and ANXA1)) related to prognosis were filtered, which may be regulated by DNA methylation, RNA methylation (m6A) and upstream transcription factors. The area under curve (AUC) of the diagnostic model based on 3-gene in the validation of our RNA-seq was 1. In the prognostic model based on 3-gene, the overall survival (OS) of high-risk patients was shorter. Combined with the clinical information of patients, a nomogram was constructed by using tumor size (pT) and risk score to quantify the prognostic risk. The age and tumor size of high-risk patients in the prognostic model were greater. In addition, the increase of tumor mutation burden (TMB) and diversity of T cell receptor (TCR), and the decrease of CD8(+) T cells in high-risk group suggest the existence of immunosuppressive microenvironment. ConclusionWe applied the scRNA-seq pipeline to focus on epithelial cells in PTC, simulated the process of tumor evolution, and revealed a prognostic prediction model based on 3 genes, which is related to tumor immune microenvironment.