Dissecting order amidst chaos of programmed cell deaths: construction of a diagnostic model for KIRC using transcriptomic information in blood-derived exosomes and single-cell multi-omics data in tumor microenvironment

Abstract

Background: Kidney renal clear cell carcinoma (KIRC) is the most frequently diagnosed subtype of renal cell carcinoma (RCC); however, the pathogenesis and diagnostic approaches for KIRC remain elusive. Using single-cell transcriptomic information of KIRC, we constructed a diagnostic model depicting the landscape of programmed cell death (PCD)-associated genes, namely cell death-related genes (CDRGs).

Methods: In this study, six CDRG categories, including apoptosis, necroptosis, autophagy, pyroptosis, ferroptosis, and cuproptosis, were collected. RNA sequencing (RNA-seq) data of blood-derived exosomes from the exoRBase database, RNA-seq data of tissues from The Cancer Genome Atlas (TCGA) combined with control samples from the GTEx databases, and single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database were downloaded. Next, we intersected the differentially expressed genes (DEGs) of the KIRC cohort from exoRBase and the TCGA databases with CDRGs and DEGs obtained from single-cell datasets, further screening out the candidate biomarker genes using clinical indicators and machine learning methods and thus constructing a diagnostic model for KIRC. Finally, we investigated the underlying mechanisms of key genes and their roles in the tumor microenvironment using scRNA-seq, single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq), and the spatial transcriptomics sequencing (stRNA-seq) data of KIRC provided by the GEO database.

Result: We obtained 1,428 samples and 216,155 single cells. After the rational screening, we constructed a 13-gene diagnostic model for KIRC, which had high diagnostic efficacy in the exoRBase KIRC cohort (training set: AUC = 1; testing set: AUC = 0.965) and TCGA KIRC cohort (training set: AUC = 1; testing set: AUC = 0.982), with an additional validation cohort from GEO databases presenting an AUC value of 0.914. The results of a subsequent analysis revealed a specific tumor epithelial cell of TRIB3^high subset. Moreover, the results of a mechanical analysis showed the relatively elevated chromatin accessibility of TRIB3 in tumor epithelial cells in the scATAC data, while stRNA-seq verified that TRIB3 was predominantly expressed in cancer tissues.

Conclusions: The 13-gene diagnostic model yielded high accuracy in KIRC screening, and TRIB3^high tumor epithelial cells could be a promising therapeutic target for KIRC.

Keywords: biomarkers; exosomes; kidney renal clear cell carcinoma; prognosis; programmed cell death; single-cell RNA sequencing; spatial transcriptome.