AUTHOR=Chen Zuoyu , Wang Yizeng , Li Dongyang , Le Yuting , Han Yue , Jia Lanning , Yan Caigu , Tian Zhigang , Song Wenbin , Li Fuxin , Zhao Ke , He Xianghui TITLE=Single-Cell RNA Sequencing Revealed a 3-Gene Panel Predicted the Diagnosis and Prognosis of Thyroid Papillary Carcinoma and Associated With Tumor Immune Microenvironment JOURNAL=Frontiers in Oncology VOLUME=12 YEAR=2022 URL=https://www.frontiersin.org/journals/oncology/articles/10.3389/fonc.2022.862313 DOI=10.3389/fonc.2022.862313 ISSN=2234-943X ABSTRACT=Objective

The 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.

Method

The 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.

Results

The 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.

Conclusion

We 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.