Development of a prognostic gene signature for hepatocellular carcinoma

Abstract

Accurate prediction of overall survival is important for prognosis and the assignment of appropriate personalized clinical treatment in hepatocellular carcinoma (HCC) patients. The aim of the present study was to establish an optimal gene model for the independent prediction of prognosis associated with common clinical patterns. Gene expression profiles and the corresponding clinical information of the LIHC cohort were obtained from The Cancer Genome Atlas. Differentially expressed genes were found using the R package “limma”. Subsequently, a prognostic gene signature was developed using the LASSO Cox regression model. Kaplan–Meier, log-rank, and receiver operating characteristic (ROC) analyses were performed to verify the predictive accuracy of the prognostic model. Finally, a nomogram and calibration plot were created using the “rms” package. Differentially expressed genes were screened with threshold criteria (FDR < 0.01 and |log FC|>3) and 563 differentially expressed genes were obtained, including 448 downregulated and 115 upregulated genes. Using the LASSO Cox regression model, a prognostic gene signature was developed based on nine genes, IQGAP3, BIRC5, PTTG1, STC2, CDKN3, PBK, EXO1, NEIL3, and HOXD9, the expression levels of which were quantitated using RT-qPCR. According to the risk scores, patients were separated into high-risk and low-risk groups. In conclusion, the prognostic gene signature can be used as a combined biomarker for the independent prediction of overall survival in HCC patients. Moreover, we created a nomogram that can be used to infer prognosis and aid individualized decisions regarding treatment and surveillance.