A Cell Cycle Progression-Derived Gene Signature to Predict Prognosis and Therapeutic Response in Hepatocellular Carcinoma

Hui, Yongfeng; Leng, Junzhi; Jin, Dong; Liu, Di; Wang, Genwang; Wang, Qi; Wang, Yanyang

doi:https://doi.org/10.1155/2021/1986159

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

The Discovery of Nucleotide and Nucleic Acid-Associated Biomarkers and Cancer Personalized Therapies

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 1986159 | https://doi.org/10.1155/2021/1986159

A Cell Cycle Progression-Derived Gene Signature to Predict Prognosis and Therapeutic Response in Hepatocellular Carcinoma

Yongfeng Hui,¹Junzhi Leng,¹Dong Jin,¹Di Liu,¹Genwang Wang,¹Qi Wang,¹and Yanyang Wang²

Academic Editor: Fu Wang

Received23 Aug 2021

Accepted22 Sept 2021

Published21 Oct 2021

Abstract

Objective. Dysregulation of cell cycle progression (CCP) is one of the hallmarks of cancer. Here, our study is aimed at developing a CCP-derived gene signature for predicting high-risk population of hepatocellular carcinoma (HCC). Methods. Our study retrospectively analyzed the transcriptome profiling and clinical information of HCC patients from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) projects. Uni- and multivariate cox regression models were conducted for identifying which hallmarks of cancer were risk factors of HCC. CCP-derived gene signature was developed with LASSO method. The predictive efficacy was verified by ROC curves and subgroup analyses. A nomogram was then generated and validated by ROC, calibration, and decisive curves. Immune cell infiltration was estimated with ssGSEA method. Potential small molecular compounds were predicted via CTRP and CMap analyses. The response to chemotherapeutic agents was evaluated based on the GDSC project. Results. Among hallmarks of cancer, CCP was identified as a dominant risk factor for HCC prognosis. CCP-derived gene signature displayed the favorable predictive efficacy in HCC prognosis independent of other clinicopathological parameters. A nomogram was generated for optimizing risk stratification and quantifying risk evaluation. CCP-derived signature was in relation to immune cell infiltration, HLA, and immune checkpoint expression. Combining CTRP and CMap analyses, fluvastatin was identified as a promising therapeutic agent against HCC. Furthermore, CCP-derived signature might be applied for predicting the response to doxorubicin and gemcitabine. Conclusion. Collectively, CCP-derived gene signature was a promising marker in prediction of survival outcomes and therapeutic responses for HCC patients.

1. Introduction

Hepatocellular carcinoma (HCC), an aggressive malignancy with undesirable prognosis, occupies 85% of liver cancer cases, which usually develops in the context of chronic liver diseases [1]. The complex etiology and highly intratumoral and intertumoral heterogeneity make prognosis prediction more challenging [2]. The current TNM staging system is not sufficient for precise therapy decision-making and prognostic prediction for HCC patients. Emerging evidence demonstrates that genomic signatures can be applied for risk stratification and prognostic prediction in HCC [3–5]. Nevertheless, due to insufficient sample size, interpatient, intertumoral, and intratumoral heterogeneity, and technical bias, most of prognostic models possess low reproducibility, which cannot be utilized for clinical routine practice [6]. Hence, new methods of identifying high-risk subgroups of HCC patients will bring great implications to personalized cancer care.

The hallmarks of cancer include six biological abilities acquired during the multistep progression of cancers, including sustained proliferation, evasion of growth inhibitors, cell death resistance, replicative immortality, angiogenesis, invasion activation, metastases, metabolic reprogramming, and immune escape [7]. Among them, cell cycle progression (CCP) is a key biological event that has a controlled regulation in normal cells, but it nearly generally becomes abnormal or dysregulated in tumor cells [8]. Evidence suggests that CCP could be a useful biomarker in prediction of recurrence and metastasis following surgical resection in several cancer types such as breast cancer and prostate cancer [9–11]. Aberrations of CCP-related genes are frequently found across diverse human neoplastic processes, and their expression profiling possesses wide ranging prognostic significance [12]. In our research, CCP was identified as a dominant risk factor of HCC survival outcomes. Accurate bioinformatic and machine learning methods were adopted for screening reliable candidate variables as well as building a personalized CCP-derived gene signature in prediction of prognosis and the response to immune- and chemotherapies for HCC patients.

2. Materials and Methods

2.1. Data Acquisition

This study gathered the mRNA expression profiling and clinical features of HCC patients from the Cancer Genome Atlas (TCGA; https://cancergenome.nih.gov/) [13] and International Cancer Genome Consortium (ICGC; https://dcc.icgc.org/) [14]. Only patients with complete follow-up time and vital status were included in this study. In total, we finally retrieved 368 HCC patients from TCGA-LIHC cohort (Supplementary table 1) and 260 HCC patients from ICGC (LIRI-JP) cohort for further analysis. Then, FPKM value was converted to TPM value. The somatic mutation data in Mutation Annotation Format (MAF) of HCC were downloaded from TCGA database, which were analyzed through maftools package [15].

2.2. Quantification of Hallmarks of Cancer

The gene sets of hallmarks of cancer raised by Hanahan and Weinberg [7] including glycolysis, CCP, angiogenesis, apoptosis, DNA repair, epithelial-mesenchymal transition (EMT), hypoxia, inflammation, and stemness were acquired from the MSigDB database (http://software.broadinstitute.org/gsea/index.jsp) [16]. Single sample gene set enrichment analysis (ssGSEA) was adopted to quantify the activity of above hallmarks of cancer with gene set variation analysis (GSVA) [17]. Immune infiltration was quantified using Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) computational method [18]. The ssGSEA and immune infiltration scores were scaled with -score. Uni- and multivariate cox regression models were conducted to evaluate the prognostic implications of above hallmarks of cancer in TCGA-LIHC cohort. Supplementary table 2 listed the gene sets of hallmarks of cancer.

2.3. Differential Expression Analysis

CCP-relevant genes were identified through presenting differential expression analysis between high and low CCP score groups via limma package. The criteria were set as and .

2.4. Functional Enrichment Analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of CCP-relevant genes were conducted through adopting clusterProfiler package [19]. GO covered three categories: biological process, cellular component, and molecular function.

2.5. Establishment and Verification of a Prognostic Gene Signature

Univariate Cox regression models were conducted to assess the association between CCP-relevant genes and HCC prognosis in TCGA-LIHC cohort. CCP-relevant genes with were included for least absolute shrinkage and selection operator- (LASSO-) penalized Cox regression analysis utilizing glmnet package [20]. Variable selection and shrinkage were carried out. Penalty parameter () of this model was determined through tenfold crossverification in line with the value that corresponded to the lowest partial likelihood deviance. The risk score (RS) was calculated following the normalized expression of each candidate variable and its matched regression coefficient. The formula was constructed as follows: . HCC patients were stratified into high- and low-risk subgroups in line with the median value of the RS. Overall survival between two subgroups was compared with survminer package. Uni- and multivariate cox regression models were established for assessing the predictive independency of the gene signature through adjustment of clinicopathological parameters (age, gender, grade, and stage). Time-dependent receiver-operating characteristic (ROC) curve analyses were conducted for evaluating the predictive power of the gene signature with survivalROC package. The CCP-derived gene signature was externally verified in the LIRI-JP cohort.

2.6. Establishment of a Prognostic Nomogram

Nomogram was established in TCGA-LIHC cohort with rms package. The 1-, 3-, and 5-year overall survival was assessed through total points, sum points of stage, and RS. The predictive performance of this nomogram was verified by ROC, calibration, and decisive curves.

2.7. Assessment of Tumor Microenvironment (TME)

Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE; https://sourceforge.net/projects/estimateproject/) method [21] was adopted for estimating the fractions of stromal and immune cells in TCGA-LIHC cohort. Through combination of stromal and immune scores, tumor purity was inferred. The ssGSEA algorithm was utilized for quantifying the relative abundance of 28 tumor-infiltrating immune cells in TCGA-LIHC cohort. The marker genes of immune cells were obtained from Charoentong et al. [22, 23]. The enrichment score was determined, which represented the infiltration of immune cells in HCC specimens. The mRNA expression of human leukocyte antigen (HLA) genes and immune checkpoints was analyzed in HCC samples.

2.8. Quantifying the Activity of Known Signaling Pathways

We collected the gene sets of IFN-γ signature, APM signal, base excision repair, cell cycle, DNA replication, Fanconi anemia pathway, homologous recombination, microRNAs in cancer, mismatch repair, nucleotide excision repair, oocyte meiosis, p53 signaling pathway, progesterone-mediated oocyte maturation, proteasome, pyrimidine metabolism, spliceosome, systemic lupus erythematosus, and viral carcinogenesis from published literature [24–26]. Spearman correlation test was presented between RS and above biological pathways.

2.9. Prediction of Potential Therapeutic Agents

This study gathered drug sensitivity profiling of human cancer cell lines from the CTRP (https://portals.broadinstitute.org/ctrp) projects. The area under the curve (AUC) value was utilized for inferring drug sensitivity. Moreover, transcriptome data in Cancer Cell Line Encyclopedia (CCLE) database (https://portals.broadinstitute.org/ccle/) were applied for CTRP analysis [27].

2.10. Screening Small Molecule Compounds

RS-relevant genes with and were screened between high- and low-risk HCC specimens. The up- and downregulated RS-relevant genes were separately uploaded into the Connectivity map (CMap; http://portals.broadinstitute.org/cmap/) project [28]. Potential small molecular compounds with and were discovered. Mechanism of action of these compounds were evaluated with mode-of-action (MoA) analyses.

2.11. Analysis of Chemotherapy Response

The response to chemotherapeutic agents was evaluated based on the Genomics of Drug Sensitivity in Cancer (GDSC, https://www.cancerrxgene.org/) project [29]. Through applying pRRophetic package, the half-maximal inhibitory concentration (IC50) values of doxorubicin, cisplatin, gemcitabine, and sorafenib were calculated in each HCC specimen [30].

2.12. Statistical Analysis

Statistical analysis was conducted with R language (version 4.0.2). Comparisons between two groups were carried out via Student’s -test or Wilcoxon test. Spearman correlation test was utilized for assessment of the correlation between parameters. Two-sided indicated statistical significance.

3. Results

3.1. CCP Acts as a Dominant Risk Factor for HCC Prognosis among Cancer Hallmarks

This study evaluated the prognostic significance of hallmarks of cancer defined by Hanahan and Weinberg. As shown in univariate cox regression models, angiogenesis (HR: 1.55 (0.45-5.29), ), apoptosis (hazard ratio (HR): 0.64 (0.04-10.96), ), EMT (HR: 2.01 (0.63-6.40), ), hypoxia (HR: 13.17 (0.86-200.85), ), inflammation (HR: 1.04 (0.26-4.15), ), and immune infiltration (HR: 0.9997 (0.9996-1), ) did not significantly affect the survival outcomes of HCC patients in TCGA-LIHC cohort. Meanwhile, DNA repair (HR: 1524850.08 (953.23-2439263968.30), ), glycolysis (HR: 3325.14 (67.71-163303.93), ), stemness (HR: 705.26 (35.39-14055.11), ), and CCP (HR: 439.59 (44.74-4319.26), ) were significant risk factors of HCC prognosis. Further multivariate cox regression models showed that glycolysis (HR: 447.32 (4.76-42064.25), ) and CCP (HR: 63.19 (1.33-3001.70), ) not DNA repair (HR: 993.12 (0.08-13051084.68), ) and stemness (HR: 0.30 (0.002-41.35), ) were independent risk factors for HCC. By adjusting clinicopathological parameters, we observed that, in addition to stage (HR: 1.44 (1.17-1.78), ) and age (HR: 1.01 (1.002-1.03), ), CCP (HR: 98.54 (7.10-1366.81), ) was the only significant risk factor of HCC prognosis among hallmarks of cancer according to multivariate cox regression models (Figure 1(a)). The ssGSEA score was used to quantify CCP level in each HCC specimen. Figure 1(b) visualized the associations of CCP score and common clinicopathological parameters. Survival analysis uncovered that high CCP score was in relation to unfavorable clinical outcomes of HCC patients (Figure 1(c)). Above findings suggested that CCP acted as a dominant risk factor for HCC prognosis among cancer hallmarks.

(a)

(b)

(c)

3.2. Identification of CCP-Relevant Genes in HCC

In total, 549 CCP-relevant genes were identified across HCC specimens in TCGA-LIHC cohort (Figures 2(a) and 2(b)). The detailed information was listed in Supplementary table 3. Biological functions of these 549 CCP-relevant genes were explored in depth. In Figure 2(c), biological processes of chromosome segregation, mitotic nuclear division, nuclear division, organelle fission, and sister chromatid segregation were distinctly regulated by the 549 CCP-relevant genes. Furthermore, the genes mainly participated in modulating the cellular components of chromosomal region, chromosome, centromeric region, condensed chromosome, condensed chromosome, centromeric region, and spindle (Figure 2(d)). They had the molecular functions of 3-5 DNA helicase activity, catalytic activity, acting on DNA, DNA helicase activity, DNA replication origin binding, and single-stranded DNA binding (Figure 2(e)). In Figure 2(f), cell cycle, DNA replication, drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, and retinol metabolism were distinctly modulated by the CCP-relevant genes.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2

Identification of CCP-relevant genes in HCC and their biological functions. (a) Volcano plots showing CCP-relevant genes between high and low CCP score groups in HCC specimens from TCGA-LIHC cohort. (b) Heat map visualizing the mRNA expression of CCP-relevant genes in HCC patients with high and low CCP score. (c–e) The networks showing the significant biological process, cellular component, and molecular function enriched by CCP-relevant genes. (f) The gene-pathway network visualizing the KEGG pathways involved in CCP-relevant genes.

3.3. Construction of a CCP-Derived Prognostic Gene Signature in HCC

Through univariate cox regression models, we identified 374 prognostic CCP-relevant genes for HCC in TCGA-LIHC cohort (Supplementary table 4). We adopted glmnet package to establish LASSO Cox regression model. Figure 3(a) depicted the change in trajectory of above prognostic CCP-relevant genes, suggesting that more independent variables possessed regression coefficients approaching zero as gradually increased. Furthermore, tenfold crossverification was presented to establish the prognostic model. The confidence interval corresponding to each was shown in Figure 3(b). Finally, six candidate genes were included in the model, and the formula was as follows: . As depicted in Figure 3(c), there were significantly positive interactions between six candidate genes. In line with the median value of RS, HCC patients were stratified into high- and low-risk groups (Figure 3(d)). More dead patients were found in high-risk group (Figure 3(e)). With the increase in RS, KIF20A, CDCA8, KPNA2, G6PD, NDRG1, and EPS8L3 expression was gradually increased (Figure 3(f)). Survival difference was evaluated between two groups. In Figure 3(g), we observed that low-risk patients possessed the distinct survival advantage. The prognostic significance of RS was further investigated by subgroup analyses. We stratified HCC patients into different subgroups in line with clinicopathological parameters. Our data suggested that patients with high RS had poorer clinical outcomes in comparison to those with low RS in each subgroup of and (Figures 3(h) and 3(i)), female and male (Figures 3(j) and 3(k)), G1-2 and G3-4 (Figures 3(l) and 3(m)), and stages I-II and stages III-IV (Figures 3(n) and 3(o)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

Figure 3

Construction of a CCP-derived prognostic gene signature for HCC in TCGA-LIHC cohort. (a) The change in trajectory of candidate variables. The -axis represented the log (lambda) of variables, and the -axis represented the coefficient of variables. (b) Partial likelihood deviance corresponding to each lambda. (c) The interactions between six candidate genes in the LASSO model. (d) Distribution of CCP-derived RS in each HCC patient. The high- and low-risk groups were separated in line with the median of RS (dotted line). (e) Distribution of survival status in high- and low-risk groups. (f) Heat map visualizing the mRNA expression of six candidate variables in high- and low-risk groups. Yellow, upregulation and blue, downregulation. (g) Prognostic analyses of HCC patients with high and low RS in TCGA-LIHC cohort. (h–o) Comparison of survival outcomes of HCC patients with high and low RS in each subgroup of and , female and male, G1-2 and G3-4, and stages I-II and stages III-IV.

3.4. The CCP-Derived Gene Signature is a Robust Prognostic Predictor of HCC

The predictive efficacy of the CCP-derived gene signature was externally verified in the LIRI-JP cohort. With the median of RS, HCC patients were stratified into two groups (Figure 4(a)). There were more dead cases in high-risk group (Figure 4(b)). Consistently, increased expression of KIF20A, CDCA8, KPNA2, G6PD, NDRG1, and EPS8L3 was found in high-risk patients (Figure 4(c)). Prognostic analyses uncovered the prominent survival advantage of high-risk patients (Figure 4(d)). Univariate cox regression models showed that stage and CCP-derived RS were in relation to HCC prognosis in TCGA-LIHC cohort (Figure 4(e)). As shown in multivariate cox regression models, stage and CCP-derived RS were both independent prognostic indicators of HCC (Figure 4(f)). The AUC values at 1-, 3-, and 5-year survival were 0.776, 0.697, and 0.619 in TCGA-LIHC cohort (Figure 4(g)). Also, in the LIRI-JP cohort, the AUC values at 1-, 3-, and 5-year survival were 0.779, 0.803, and 0.762 (Figure 4(h)). Above findings were indicative of the well-predictive performance of the CCP-derived RS in HCC prognosis.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 4

The CCP-derived gene signature is a robust prognostic predictor of HCC. (a) Distribution of CCP-derived RS across HCC patients in the LIRI-JP cohort. The high- and low-risk groups were separated in line with the median of RS (dotted line). (b) Distribution of survival status in high- and low-risk groups in the LIRI-JP cohort. (c) Heat map visualizing the mRNA expression of six candidate variables in high- and low-risk groups the LIRI-JP cohort. Yellow, upregulation and blue, downregulation. (d) Prognostic analyses of HCC patients with high and low RS in the LIRI-JP cohort. (e, f) Uni- and multivariate cox regression models of CCP-derived RS and clinicopathological parameters in TCGA-LIHC cohort. (g) ROC curves at 1-, 3-, and 5-year survival in TCGA-LIHC cohort. (h) ROC curves at 1-, 3-, and 5-year survival in the LIRI-JP cohort.

3.5. Establishment of a Reliable Nomogram for HCC Prognosis

Two independent prognostic factors (CCP-derived RS and stage) were included for construction of a nomogram in TCGA-LIHC cohort (Figure 5(a)). Each level of RS or stage was assigned one score, and the total score was counted through summing up the scores in each patient. The 1-, 3-, and 5-year survival probabilities were determined by the function conversion relationships of total scores. The AUCs of this nomogram at 1-, 3-, and 5-year survival were 0.718, 0.757, and 0.757, indicative of the excellent predictive performance (Figure 5(b)). The calibration curves showed the high consistency between the nomogram-estimated 1-, 3-, and 5-year clinical outcomes and actual observations (Figure 5(c)). Decisive curves were conducted for evaluation of the guiding implications of the nomogram for clinical application. Our data demonstrated that the nomogram was the best for prediction of clinical outcomes of HCC (Figure 5(d)).

(a)

(b)

(c)

(d)

3.6. CCP-Derived Prognostic Gene Signature is in Relation to TME and Carcinogenic Pathway of HCC

ESTIMATE computational method was utilized for inferring the overall infiltration levels of immune and stromal cells. We observed no significant difference in immune score or tumor purity between high- and low-risk HCC patients (Figure 6(a)). Low-risk patients displayed significantly increased stromal score in comparison to high-risk patients. In Figure 6(b), activated CD4+ T cell, central memory CD4+ T cell, memory B cell, type 17 helper cell, type 2 helper cell, and activated dendritic cell had increased infiltration levels in high-risk than low-risk HCC specimens. Reduced infiltration levels of effector memory CD8+ T cell, CD56 bright natural killer cell, eosinophil, and mast cell were found in high-risk compared to low-risk HCC samples. Furthermore, high RS was in relation to increased mRNA expression of HLA family, including HLA-DOB, HLA-DMA, HLA-DPB1, HLA-DRA, HLA-DOA, HLA-DQA2, HLA-DQA1, HLA-DMB, HLA-DPB2, and HLA-DQB2 (Figure 6(c)). We also compared the mRNA expression of immune checkpoints between high- and low-risk HCC specimens. We observed that CD200, VSIR, CD27, LAG3, CD70, TNFRSF25, CD28, NRP1, CD200R1, TIGIT, VTCN1, TNFSF18, TNFRSF18, ICOS, TNFRSF14, TNFSF9, PDCD1, TNFRSF9, CTLA4, CD86, LGALS9, TNFRSF18, HAVCR2, LAIR1, TNFRSF4, TNFSF15, HHLA2, CD80, TNFSF4, and CD276 expression was gradually increased with the increase of RS (Figure 6(d)). Association between RS and known biological pathways was evaluated across HCC samples. In Figure 6(e), RS displayed positive correlation to carcinogenic pathway activation including IFN-γ signature, APM signal, base excision repair, cell cycle, DNA replication, Fanconi anemia pathway, homologous recombination, microRNAs in cancer, mismatch repair, nucleotide excision repair, oocyte meiosis, p53 signaling pathway, progesterone-mediated oocyte maturation, proteasome, pyrimidine metabolism, spliceosome, systemic lupus erythematosus, and viral carcinogenesis.

(a)

(b)

(c)

(d)

(e)

Figure 6

CCP-derived prognostic gene signature is in relation to TME and carcinogenic pathway of HCC in TCGA-LIHC cohort. (a) Comparison of immune score, stromal score, and tumor purity in high- and low-risk HCC patients with ESTIMATE algorithm. (b) Comparison of the infiltration levels of tumor-infiltrating immune cells in high- and low-risk HCC patients via ssGSEA method. ; ; . (c) Heat map visualizing the mRNA expression of HLA genes in each HCC specimen and their relationships to CCP-derived RS. (d) Heat map showing the mRNA expression of immune checkpoints in each HCC sample as well as their relationships to CCP-derived RS. (e) Association between CCP-derived RS and known biological pathways.

3.7. Prediction of Potential Therapeutic Compounds against HCC Based on CCP-Derived Gene Signature

This study further discovered promising therapeutic agents against HCC based on the CCP-derived gene signature in TCGA-LIHC cohort. In Figure 7(a), 12 CTRP-derived compounds were identified, including paclitaxel (), BI-2536 (), triazolothiadiazine (), ABT-737 (), nakiterpiosin (), mitomycin (), barasertib (), ceranib-2 (), SB-743921 (), clofarabine (), fluvastatin (), and KX2-391 (). AUC values of above compounds were compared between high- and low-risk HCC patients. We observed that high-risk patients had prominently reduced AUC values of each compound. This indicated that high-risk patients were more likely to respond to above compounds. Small molecular compounds were also predicted by CMap project. In total, 33 small molecular compounds with and were identified (Table 1). Shared mechanisms of action of above compounds were shown in Figure 7(b). We observed that digitoxigenin, digoxin, helveticoside, ouabain, and proscillaridin shared the ATPase inhibitor; apigenin and ricinine shared casein kinase inhibitor; 4, 5-dianilinophthalimide and butein shared EGFR inhibitor; scriptaid and vorinostat shared HDAC inhibitor; fluvastatin and simvastatin shared HMGCR inhibitor; and LY-294002 and sirolimus shared mTOR inhibitor. Fluvastatin was predicted both in CTRP and CMap projects, indicative of the potential as a therapeutic agent against HCC.

(a)

(b)

(c)

Figure 7

Prediction of potential therapeutic compounds as well as assessment of chemosensitivity based on CCP-derived gene signature in TCGA-LIHC cohort. (a) Spearman correlation between CTRP-derived compounds and CCP-derived RS as well as comparison of estimated AUC value between high- and low-risk HCC samples. . (b) Shared mechanisms of action of small molecular compounds derived from CMap database. (c) Comparison of the estimated IC50 value of doxorubicin, cisplatin, gemcitabine, and sorafenib in high- and low-risk HCC samples.

3.8. CCP-Derived Gene Signature Predicts Chemotherapeutic Response

The response to chemotherapeutic agents including doxorubicin, cisplatin, gemcitabine, and sorafenib in each HCC sample was evaluated in TCGA-LIHC cohort. We compared the IC50 value of above agents in high- and low-risk HCC patients. As a result, high-risk patients had markedly reduced IC50 values of doxorubicin and gemcitabine than those with low-risk, suggesting that high RS was indicative of increased sensitivity to doxorubicin and gemcitabine (Figure 7(c)). However, no significant difference was investigated between high- and low-risk HCC patients. Thus, CCP-derived prognostic gene signature might be applied for predicting the response to doxorubicin and gemcitabine.

3.9. CCP-Derived Gene Signature Acts as a Reliable Prognostic Indicator Independent of Genetic Mutation

Further analysis revealed that 84.62% occurred somatic mutation across 364 HCC samples (Figure 8(a)). TP53 was the most frequent mutation (30%), followed by CTNNB1 (25%), TTN (24%), and MUC16 (14%). Missense mutation was the main type of mutation. HCC samples were stratified into different subgroups according to whether above genes occurred mutation. We observed that high-RS patients displayed poorer clinical outcomes in comparison to low-RS patients regardless of whether TP53, CTNNB1, TTN, and MUC16 were mutated (Figures 8(b)–8(i)). This indicated that CCP-derived gene signature acted as a reliable prognostic indicator independent of genetic mutation.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

4. Discussion

Nowadays, bioinformatics analysis has become an important tool for cancer research. Herein, our study suggested that CCP acted as a dominant risk factor for HCC prognosis among hallmarks of cancer. The LASSO computational method is a shrinkage estimate, which may be utilized for constructing a penalty function as well as obtaining a relatively refined model, where several coefficients may be shrunk, and several are set to zero. Hence, this model possesses the advantage of subset shrinkage and represents a biased estimate, which may process multiple collinear data. It is capable of estimating parameters and selecting variables, thereby solving the problem of multiple collinearities in regression analyses. Thus, we adopted LASSO method to construct a robust CCP-derived gene signature for HCC prognosis.

This signature contained KIF20A, CDCA8, KPNA2, G6PD, NDRG1, and EPS8L3. KIF20A dysregulation may be independently predictive of unfavorable clinical outcomes and recurrence of HCC patients [31]. KIF20A downregulation may decrease the proliferation and induce the G1 arrest of HCC cells [32]. CDCA8 is in relation to unfavorable stage and survival outcomes of HCC [33]. CDCA8 inhibition weakens HCC growth and stemness through restoring ATF3 and inactivating AKT/β-catenin axis [34]. KPNA2 upregulation accelerates HCC progression via enhancement of migration [35] and proliferation and is indicative of undesirable prognosis [36]. KPNA2 dysregulation is also in relation to early recurrence for patients with small HCC following hepatectomy [37]. G6PD upregulation leads to migration and invasion of HCC cells through enhancing EMT [38]. It can inhibit ferroptosis in HCC through reducing cytochrome P450 oxidoreductase [39]. Hypoxia-mediated NDRG1 upregulation modulates apoptotic levels via driving mitochondrial fission in HCC [40]. NDRG1 acts as a predictor of metastases, relapse, and undesirable prognosis in HCC [41]. EPS8L3 upregulation enhances HCC proliferation through inhibition of the transactivity of FOXO1 [42]. It can also induce HCC proliferation and metastases through modulation of EGFR dimerization and internalization [43]. Also, we observed that the CCP-derived gene signature exhibited positive correlation to carcinogenic pathway activation. Hence, each variable in the CCP-derived gene signature acts as an oncogene during HCC progression.

Our prognostic analysis uncovered that the CCP-derived gene signature can be independently predictive of HCC patients’ clinical outcomes. Patients with high RS indicated unfavorable prognosis. ROC curves at 1-, 3-, and 5-year survival confirmed the well predictive efficacy. The external verification demonstrated the clinical applicability of this CCP-derived gene signature. A recent study has proposed a novel cell cycle-relevant prognostic model for endometrial cancer [44]. This study is the first to construct the CCP-derived gene signature in HCC prognosis. Through integration of stage and CCP-derived signature, we developed a nomogram. Following verification of ROC, calibration, and decisive curves, this nomogram possessed the potential as a clinical tool for predicting HCC prognosis.

Recently, immune checkpoint inhibitors like PD-1, PD-L1, and CTLA-4 antibodies display favorable therapeutic effects against advanced HCC [45]. But current immune-based therapies induce durable response in a minority of HCC subjects. Due to the indispensable role of TME in HCC pathogenesis, it is of significance to uncover the underlying mechanism shaping the unique TME of HCC [46]. Our findings suggested that high-risk HCC possessed the increased infiltration levels of activated CD4+ T cell, central memory CD4+ T cell, memory B cell, type 17 helper cell, type 2 helper cell, and activated dendritic cell as well as had the reduced infiltration levels of effector memory CD8+ T cell, CD56 bright natural killer cell, eosinophil, and mast cell. Furthermore, most of HLA genes and immune checkpoints displayed the significant upregulation in high-risk patients. Hence, the CCP-derived gene signature was in relation to TME reshaping.

Combining CTRP and CMap analyses, we finally identified fluvastatin as a promising therapeutic agent. Evidence suggests that fluvastatin suppresses apoptosis of HCC cells [47]. High RS had increased sensitivity to doxorubicin and gemcitabine for HCC patients. This indicated that CCP-derived gene signature might be utilized for prediction of response to doxorubicin and gemcitabine. Furthermore, we conducted the genetic mutation landscape of HCC. Our data suggested that the CCP-derived gene signature acted as a reliable prognostic predictor independent of genetic mutation. However, several limitations of this study should be pointed out. Although the CCP-derived gene signature has been verified externally, the predictive performance should be observed in a prospective cohort. Moreover, each gene in the CCP-derived gene signature will be validated through in vivo or in vitro experiments.

5. Conclusion

Collectively, our study proposed the CCP-derived gene signature for HCC prognosis. Following verification, this signature acted as a robust prognostic predictor for HCC. Also, a reliable nomogram including the CCP-derived gene signature and stage was developed as a promising clinical tool. Further analysis uncovered that the CCP-derived gene signature was in relation to TME reshaping and response to doxorubicin and gemcitabine. Based on the CCP-derived gene signature, fluvastatin was confirmed as a therapeutic agent against HCC. Despite this, the therapeutic effects of fluvastatin required to be verified in more experiments. Moreover, the predictive efficacy of the CCP-derived gene signature in HCC clinical outcomes will be verified in prospective cohorts.

Abbreviations

HCC:	Hepatocellular carcinoma
TCGA:	The Cancer Genome Atlas
ICGC:	International Cancer Genome Consortium
CCP:	Cell cycle progression
EMT:	Epithelial-mesenchymal transition
ssGSEA:	Single sample gene set enrichment analysis
GSVA:	Gene set variation analysis
GO:	Gene Ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes
LASSO:	Least absolute shrinkage and selection operator
RS:	Risk score
ROC:	Receiver-operating characteristic
TME:	Tumor microenvironment
ESTIMATE:	Estimation of STromal and Immune cells in MAlignant Tumours using Expression data
HLA:	Human leukocyte antigen
AUC:	Area under the curve
CCLE:	Cancer Cell Line Encyclopedia
CMap:	Connectivity map
MoA:	Mode-of-action
GDSC:	Genomics of Drug Sensitivity in Cancer
IC50:	Half-maximal inhibitory concentration.

Data Availability

The data used to support the findings of this study are included within the supplementary information files.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Yongfeng Hui and Junzhi Leng contributed equally to this work.

Acknowledgments

This research was supported by the Ningxia Autonomous Region Key R&D Program Project (2021BEG03067), Natural Science Foundation of Ningxia (2019AAC03192), Key Project of Ningxia Health and Family Planning Commission (2019-NW-015), and General Hospital of Ningxia Medical University (Clinical Medicine Research Center) Open Project.

Supplementary Materials

Supplementary table 1. Clinical information of HCC patients in TCGA-LIHC cohort. Supplementary table 2. The gene sets of hallmarks of cancer. Supplementary table 3. 549 CCP-relevant genes in TCGA-LIHC cohort. Supplementary table 4. Prognostic CCP-relevant genes in TCGA-LIHC cohort. (Supplementary Materials)

References

A. J. Craig, J. von Felden, T. Garcia-Lezana, S. Sarcognato, and A. Villanueva, “Tumour evolution in hepatocellular carcinoma,” Nature Reviews Gastroenterology & Hepatology, vol. 17, no. 3, pp. 139–152, 2020.
View at: Publisher Site | Google Scholar
B. Tang, J. Zhu, J. Li et al., “The ferroptosis and iron-metabolism signature robustly predicts clinical diagnosis, prognosis and immune microenvironment for hepatocellular carcinoma,” Cell Communication and Signaling, vol. 18, no. 1, 2020.
View at: Publisher Site | Google Scholar
W. Hong, L. Liang, Y. Gu et al., “Immune-related lncRNA to construct novel signature and predict the immune landscape of human hepatocellular carcinoma,” Molecular Therapy - Nucleic Acids, vol. 22, pp. 937–947, 2020.
View at: Publisher Site | Google Scholar
J. Y. Liang, D. S. Wang, H. C. Lin et al., “A novel ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma,” International Journal of Biological Sciences, vol. 16, no. 13, pp. 2430–2441, 2020.
View at: Publisher Site | Google Scholar
Y. Liu, X. Zhang, J. Zhang, J. Tan, J. Li, and Z. Song, “Development and validation of a combined ferroptosis and immune prognostic classifier for hepatocellular carcinoma,” Frontiers in Cell and Developmental Biology, vol. 8, p. 596679, 2020.
View at: Publisher Site | Google Scholar
H. Xiao, B. Wang, H. X. Xiong et al., “A novel prognostic index of hepatocellular carcinoma based on immunogenomic landscape analysis,” Journal of Cellular Physiology, vol. 236, no. 4, pp. 2572–2591, 2021.
View at: Publisher Site | Google Scholar
D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.
View at: Publisher Site | Google Scholar
G. I. Evan and K. H. Vousden, “Proliferation, cell cycle and apoptosis in cancer,” Nature, vol. 411, no. 6835, pp. 342–348, 2001.
View at: Publisher Site | Google Scholar
M. R. Cooperberg, J. P. Simko, J. E. Cowan et al., “Validation of a cell-cycle progression gene panel to improve risk stratification in a contemporary prostatectomy cohort,” Journal of Clinical Oncology, vol. 31, no. 11, pp. 1428–1434, 2013.
View at: Publisher Site | Google Scholar
J. D. Mosley and R. A. Keri, “Cell cycle correlated genes dictate the prognostic power of breast cancer gene lists,” BMC Medical Genomics, vol. 1, no. 1, p. 11, 2008.
View at: Publisher Site | Google Scholar
J. Cuzick, G. P. Swanson, G. Fisher et al., “Prognostic value of an RNA expression signature derived from cell cycle proliferation genes in patients with prostate cancer: a retrospective study,” Lancet Oncology, vol. 12, no. 3, pp. 245–255, 2011.
View at: Publisher Site | Google Scholar
J. M. West, D. Ma, S. L. Mott, and J. A. Brown, “Cell cycle progression score has potential prognostic value for stage T1 renal cell carcinomas,” Urologic Oncology, vol. 38, no. 5, pp. 545–552, 2020.
View at: Publisher Site | Google Scholar
Z. Wang, M. A. Jensen, and J. C. Zenklusen, “A practical guide to The Cancer Genome Atlas (TCGA),” Methods in Molecular Biology, vol. 1418, pp. 111–141, 2016.
View at: Publisher Site | Google Scholar
J. Zhang, R. Bajari, D. Andric et al., “The international cancer genome consortium data portal,” Nature Biotechnology, vol. 37, no. 4, pp. 367–369, 2019.
View at: Publisher Site | Google Scholar
A. Mayakonda, D. C. Lin, Y. Assenov, C. Plass, and H. P. Koeffler, “Maftools: efficient and comprehensive analysis of somatic variants in cancer,” Genome Research, vol. 28, no. 11, pp. 1747–1756, 2018.
View at: Publisher Site | Google Scholar
A. Liberzon, C. Birger, H. Thorvaldsdóttir, M. Ghandi, J. P. Mesirov, and P. Tamayo, “The Molecular Signatures Database Hallmark Gene Set Collection,” Cell Systems, vol. 1, no. 6, pp. 417–425, 2015.
View at: Publisher Site | Google Scholar
S. Hänzelmann, R. Castelo, and J. Guinney, “GSVA: gene set variation analysis for microarray and RNA-seq data,” BMC Bioinformatics, vol. 14, no. 1, 2013.
View at: Publisher Site | Google Scholar
A. M. Newman, C. L. Liu, M. R. Green et al., “Robust enumeration of cell subsets from tissue expression profiles,” Nature Methods, vol. 12, no. 5, pp. 453–457, 2015.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
S. Engebretsen and J. Bohlin, “Statistical predictions with glmnet,” Clinical Epigenetics, vol. 11, no. 1, 2019.
View at: Publisher Site | Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, 2013.
View at: Publisher Site | Google Scholar
D. A. Barbie, P. Tamayo, J. S. Boehm et al., “Systematic RNA interference reveals that oncogenic _KRAS_ -driven cancers require TBK1,” Nature, vol. 462, no. 7269, pp. 108–112, 2009.
View at: Publisher Site | Google Scholar
P. Charoentong, F. Finotello, M. Angelova et al., “Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade,” Cell Reports, vol. 18, no. 1, pp. 248–262, 2017.
View at: Publisher Site | Google Scholar
Y. Şenbabaoğlu, R. S. Gejman, A. G. Winer et al., “Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures,” Genome Biol, vol. 17, no. 1, 2016.
View at: Publisher Site | Google Scholar
J. E. Rosenberg, J. Hoffman-Censits, T. Powles et al., “Atezolizumab in patients with locally advanced and metastatic urothelial carcinoma who have progressed following treatment with platinum-based chemotherapy: a single-arm, multicentre, phase 2 trial,” Lancet, vol. 387, no. 10031, pp. 1909–1920, 2016.
View at: Publisher Site | Google Scholar
S. Mariathasan, S. J. Turley, D. Nickles et al., “TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells,” Nature, vol. 554, no. 7693, pp. 544–548, 2018.
View at: Publisher Site | Google Scholar
M. Ghandi, F. W. Huang, J. Jané-Valbuena et al., “Next-generation characterization of the Cancer Cell Line Encyclopedia,” Nature, vol. 569, no. 7757, pp. 503–508, 2019.
View at: Publisher Site | Google Scholar
J. Lamb, E. D. Crawford, and D. Peck, “The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease,” Science, vol. 313, no. 5795, pp. 1929–1935, 2006.
View at: Publisher Site | Google Scholar
W. Yang, J. Soares, P. Greninger et al., “Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells,” Nucleic Acids Research, vol. 41, no. D1, pp. D955–D961, 2012.
View at: Publisher Site | Google Scholar
P. Geeleher, N. Cox, and R. S. Huang, “pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels,” PLoS One, vol. 9, no. 9, article e107468, 2014.
View at: Publisher Site | Google Scholar
M. Lu, X. Huang, Y. Chen et al., “AberrantKIF20Aexpression might independently predict poor overall survival and recurrence-free survival of hepatocellular carcinoma,” IUBMB Life, vol. 70, no. 4, pp. 328–335, 2018.
View at: Publisher Site | Google Scholar
X. Li, W. Huang, W. Huang et al., “Kinesin family members KIF2C/4A/10/11/14/18B/20A/23 predict poor prognosis and promote cell proliferation in hepatocellular carcinoma,” American Journal of Translational Research, vol. 12, pp. 1614–1639, 2020.
View at: Google Scholar
D. Li, Y. Ji, J. Guo, and Q. Guo, “Upregulated expression ofMTFR2as a novel biomarker predicts poor prognosis in hepatocellular carcinoma by bioinformatics analysis,” Future Oncology, vol. 17, no. 24, pp. 3187–3201, 2021.
View at: Publisher Site | Google Scholar
T. Jeon, M. J. Ko, Y. R. Seo et al., “Silencing CDCA8 suppresses hepatocellular carcinoma growth and stemness via restoration of ATF3 tumor suppressor and inactivation of AKT/β-catenin signaling,” Cancers, vol. 13, no. 5, p. 1055, 2021.
View at: Publisher Site | Google Scholar
X. Guo, Z. Wang, J. Zhang et al., “Upregulated KPNA2 promotes hepatocellular carcinoma progression and indicates prognostic significance across human cancer types,” Acta Biochimica et Biophysica Sinica, vol. 51, no. 3, pp. 285–292, 2019.
View at: Publisher Site | Google Scholar
Y. Zan, B. Wang, L. Liang et al., “MicroRNA-139 inhibits hepatocellular carcinoma cell growth through down-regulating karyopherin alpha 2,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, 2019.
View at: Publisher Site | Google Scholar
P. Jiang, Y. Tang, L. He et al., “Aberrant expression of nuclear KPNA2 is correlated with early recurrence and poor prognosis in patients with small hepatocellular carcinoma after hepatectomy,” Medical Oncology, vol. 31, no. 8, p. 131, 2014.
View at: Publisher Site | Google Scholar
M. Lu, L. Lu, Q. Dong et al., “Elevated G6PD expression contributes to migration and invasion of hepatocellular carcinoma cells by inducing epithelial-mesenchymal transition,” Acta Biochimica et Biophysica Sinica, vol. 50, no. 4, pp. 370–380, 2018.
View at: Publisher Site | Google Scholar
F. Cao, A. Luo, and C. Yang, “G6PD inhibits ferroptosis in hepatocellular carcinoma by targeting cytochrome P450 oxidoreductase,” Cell Signal, vol. 87, p. 110098, 2021.
View at: Publisher Site | Google Scholar
D. D. Guo, K. F. Xie, and X. J. Luo, “Hypoxia-induced elevated NDRG1 mediates apoptosis through reprograming mitochondrial fission in HCC,” Gene, vol. 741, p. 144552, 2020.
View at: Publisher Site | Google Scholar
J. Cheng, H. Y. Xie, X. Xu et al., “NDRG1 as a biomarker for metastasis, recurrence and of poor prognosis in hepatocellular carcinoma,” Cancer Letters, vol. 310, no. 1, pp. 35–45, 2011.
View at: Publisher Site | Google Scholar
C. X. ZENG, L. Y. TANG, C. Y. XIE et al., “Overexpression of EPS8L3 promotes cell proliferation by inhibiting the transactivity of FOXO1 in HCC,” Neoplasma, vol. 65, no. 5, pp. 701–707, 2018.
View at: Publisher Site | Google Scholar
Z. Xuan, L. Zhao, Z. Li et al., “EPS8L3 promotes hepatocellular carcinoma proliferation and metastasis by modulating EGFR dimerization and internalization,” American journal of cancer research, vol. 10, no. 1, pp. 60–77, 2020.
View at: Google Scholar
J. Liu, J. Mei, S. Li, Z. Wu, and Y. Zhang, “Establishment of a novel cell cycle-related prognostic signature predicting prognosis in patients with endometrial cancer,” Cancer Cell International, vol. 20, no. 1, 2020.
View at: Publisher Site | Google Scholar
Y. Zongyi and L. Xiaowu, “Immunotherapy for hepatocellular carcinoma,” Cancer Letters, vol. 470, pp. 8–17, 2020.
View at: Publisher Site | Google Scholar
B. Ruf, B. Heinrich, and T. F. Greten, “Immunobiology and immunotherapy of HCC: spotlight on innate and innate-like immune cells,” Cellular & Molecular Immunology, vol. 18, no. 1, pp. 112–127, 2021.
View at: Publisher Site | Google Scholar
L. H. al-Wahaibi, M. S. M. al-Saleem, O. A. A. Ahmed et al., “Optimized conjugation of fluvastatin to HIV-1 TAT displays enhanced pro-apoptotic activity in HepG2 cells,” International Journal of Molecular Sciences, vol. 21, no. 11, p. 4138, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yongfeng Hui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

319

Downloads

580

Citations