Establishment and validation of a novel CD8+ T cell-associated prognostic signature for predicting clinical outcomes and immunotherapy response in hepatocellular carcinoma via integrating single-cell RNA-seq and bulk RNA-seq

Abstract

CD8+ T lymphocytes are critical in the immune response against neoplasms, yet the prognostic relevance of CD8+ T cell-associated genes in hepatocellular carcinoma (HCC) is not fully understood. We sourced single-cell RNA-sequencing (scRNA-seq) and bulk RNA-seq data for HCC from the GSE98638 dataset and The Cancer Genome Atlas (TCGA) repository. We utilized Weighted Gene Correlation Network Analysis (WGCNA) to identify CD8+ T cell-related genes. A clinical prognostic model for risk stratification was then constructed via Cox-Lasso regression analysis. The Immunophenotypic Score (IPS) was utilized to evaluate the potential of immunotherapeutic interventions in the categorized cohorts. Validation of the expression of CD8+ T cell-associated risk genes was performed using quantitative reverse transcription PCR (qRT-PCR). Integrating scRNA-seq with RNA-seq data, we identified five CD8+ T cell-related signature genes: IKBKE, ATP1B3, MSC, ADA, and BATF. Notably, HCC patients in the high-risk group had markedly decreased overall survival. Elevated infiltration levels of CD8+ T cells, B cells, and macrophages were observed in the high-risk group. Moreover, there was a positive correlation between the risk score and immune checkpoints (ICPs), including PDCD1, CD274, and CTLA4. Patients within the high-risk group subject to PD1 and CTLA4 blockade exhibited higher IPS levels. Additionally, the expression of the five risk genes was upregulated in HCC cell lines and tissues compared to normal cells and tissues. Our findings establish a prognostic signature based on CD8+ T cells, offering a potent predictive model for clinical outcomes and responsiveness to immunotherapy in HCC patients.