Identifying Apoptosis-Related Transcriptomic Aberrations and Revealing Clinical Relevance as Diagnostic and Prognostic Biomarker in Hepatocellular Carcinoma

Abstract

In view of the unsatisfactory treatment outcome of liver cancer under current treatment, where the mortality rate is high and the survival rate is poor, in this study we aimed to use RNA sequencing data to explore potential molecular markers that can be more effective in predicting diagnosis and prognosis of hepatocellular carcinoma. RNA sequencing data and corresponding clinical information were obtained from multiple databases. After matching with the apoptotic genes from the Deathbase database, 14 differentially expressed human apoptosis genes were obtained. Using univariate and multivariate Cox regression analyses, two apoptosis genes (BAK1 and CSE1L) were determined to be closely associated with overall survival (OS) in HCC patients. And subsequently experiments also validated that knockdown of BAK1 and CSE1L significantly inhibited cell proliferation and promoted apoptosis in the HCC. Then the two genes were used to construct a prognostic signature and diagnostic models. The high-risk group showed lower OS time compared to low-risk group in the TCGA cohort (P < 0.001, HR = 2.11), GSE14520 cohort (P = 0.003, HR = 1.85), and ICGC cohort (P < 0.001, HR = 4). And the advanced HCC patients showed higher risk score and worse prognosis compared to early-stage HCC patients. Moreover, the prognostic signature was validated to be an independent prognostic factor. The diagnostic models accurately predicted HCC from normal tissues and dysplastic nodules in the training and validation cohort. These results indicated that the two apoptosis-related signature effectively predicted diagnosis and prognosis of HCC and may serve as a potential biomarker and therapeutic target for HCC.