Developing an advanced diagnostic model for hepatocellular carcinoma through multi-omics integration leveraging diverse cell-death patterns

Introduction: Hepatocellular carcinoma (HCC), representing more than 80% of primary liver cancer cases, lacks satisfactory etiology and diagnostic methods. This study aimed to elucidate the role of programmed cell death-associated genes (CDRGs) in HCC by constructing a diagnostic model using single-cell RNA sequencing (scRNA-seq) and RNA sequencing (RNA-seq) data. Methods: Six categories of CDRGs, including apoptosis, necroptosis, autophagy, pyroptosis, ferroptosis, and cuproptosis, were collected. RNA-seq data from blood-derived exosomes were sourced from the exoRBase database, RNA-seq data from cancer tissues from the TCGA database, and scRNA-seq data from the GEO database. Subsequently, we intersected the differentially expressed genes (DEGs) of the HCC cohort from exoRBase and TCGA databases with CDRGs, as well as DEGs obtained from single-cell datasets. Candidate biomarker genes were then screened using clinical indicators and a machine learning approach, resulting in the construction of a seven-gene diagnostic model for HCC. Additionally, scRNA-seq and spatial transcriptome sequencing (stRNA-seq) data of HCC from the Mendeley data portal were used to investigate the underlying mechanisms of these seven key genes and their association with immune checkpoint blockade (ICB) therapy. Finally, we validated the expression of key molecules in tissues and blood-derived exosomes through quantitative Polymerase Chain Reaction (qPCR) and immunohistochemistry experiments. Results: Collectively, we obtained a total of 50 samples and 104,288 single cells. Following the meticulous screening, we established a seven-gene diagnostic model for HCC, demonstrating high diagnostic efficacy in both the exoRBase HCC cohort (training set: AUC = 1; testing set: AUC = 0.847) and TCGA HCC cohort (training set: AUC = 1; testing set: AUC = 0.976). Subsequent analysis revealed that HCC cluster 3 exhibited a higher stemness index and could serve as the starting point for the differentiation trajectory of HCC cells, also displaying more abundant interactions with other cell types in the microenvironment. Notably, key genes TRIB3 and NQO1 displayed elevated expression levels in HCC cells. Experimental validation further confirmed their elevated expression in both tumor tissues and blood-derived exosomes of cancer patients. Additionally, stRNA analysis not only substantiated these findings but also suggested that patients with high TRIB3 and NQO1 expression might respond more favorably to ICB therapy. Conclusions: The seven-gene diagnostic model demonstrated remarkable accuracy in HCC screening, with TRIB3 emerging as a promising diagnostic tool and therapeutic target for HCC.