Establishment and validation of a novel risk model based on CD8T cell marker genes to predict prognosis in thyroid cancer by integrated analysis of single-cell and bulk RNA-sequencing.

Abstract

Papillary thyroid cancer (PTC) is a histological type of thyroid cancer, and CD8T is important for the immune response. The single-cell RNA data were acquired from Gene Expression Omnibus. SingleR package was used for cluster identification, and CellChat was exploited to evaluate the interaction among several cell types. Bulk RNA data obtained from the cancer genome atlas were used for determination of prognosis using Kaplan-Meier and Receiver Operating Characteristic curve. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were applied for assessment of function enrichment. The drug sensitivity was calculated in Gene Set Cancer Analysis. The regulatory network was constructed by STRING and Cytoscape. We identified 23 cell clusters and 10 cell types. Cell communication results showed CD8T cell was vital among all immune cell types. Enrichment analysis found the marker genes of CD8T cell was enriched in some signal pathways related to tumor development. Overall, FAM107B and TUBA4A were considered as hub genes and used to construct a risk model. Most immune checkpoint expressions were upregulated in tumor group. Tumor mutation burden results indicated that prognosis of PTC was not related to the mutation of hub genes. Drug sensitivity analysis showed some drugs could be effectively used for the treatment of PTC, and regulatory network identified some targets for the immunotherapy. A 2-gene model of PTC was developed based on the single-cell RNA and bulk RNA data. Besides, we found CD8T was essential for the immune response in PTC.