Abstract Background Hepatocellular carcinoma (HCC) is one of the most common human malignant tumors. The prognosis of HCC patients is still unsatisfying. In this study, we performed the integrated bioinformatics analysis to identify potential biomarkers and biological pathways in HCC. Methods Gene expression profiles were obtained from the Gene Expression Omnibus database ([32]GSE55048, [33]GSE55758, and [34]GSE56545) for the screening of the common differentially expressed genes (DEGs) between HCC tissues and matched non-tumor tissues. DEGs were subjected to Gene Ontology, KEGG pathway, and Reactome pathway analysis. The hub genes were identified by using protein–protein interaction (PPI) network analysis. The hub genes in HCC were further subjected to overall survival analysis of HCC patients. The hub genes were further validated by in vitro functional assays. Results A total of 544 common differentially expressed genes were screened from three datasets. Gene Ontology, KEGG and Reactome analysis results showed that DEGs are significantly associated with the biological process of cell cycle, cell division, and DNA replication. PPI network analysis identified 20 hub genes from the DEGs. These hub genes except CENPE were all significantly up-regulated in the HCC tissues when compared to non-tumor tissues. The Kaplan–Meier survival analysis results showed that the high expression of the 20 hub genes was associated with shorter survival of the HCC patients. Further validation studies showed that knockdown of KIF14 and KIF23 both suppressed the proliferative potential, increased the caspase-3/-7 activity, up-regulated Bax expression, and promoted the invasive and migratory abilities in the HCC cells. In addition, knockdown of KIF14 and KIF23 enhanced chemosensitivity to cisplatin and sorafenib in the HCC cells. Finally, the high expression of KIF14 and KIF23 was associated with shorter progression-free survival, recurrence-free survival, and disease-specific survival of patients with HCC. Conclusion In conclusion, the present study performed the integrated bioinformatics analysis and showed that KIF14 and KIF23 silence attenuated cell proliferation, invasion, and migration, and promoted chemosensitivity of HCC cells. KIF14 and KIF23 may serve as potential biomarkers for predicting the worse prognosis of patients with HCC. Keywords: hepatocellular carcinoma, differentially expressed genes, hub genes, KIF14, KIF23, survival analysis Introduction Hepatocellular carcinoma (HCC) is one of the most common human malignancies and is the second leading cause of cancer-related deaths,[35]^1 and nearly 50% of the new cases were diagnosed in China.[36]^2 The high incidence of HCC in China has been suggested to be associated with Hepatitis B virus infection in the Chinese populations.[37]^2 The prognosis of the HCC patients is relatively poor, due to the insufficient tools for the early diagnosis of HCC, and a substantial proportion of HCC patients was diagnosed at a late stage.[38]^3 So far, surgical intervention and radio-/chemotherapy were the main strategy for the treatment of HCC. However, the high recurrence of HCC significantly affects the overall survival of HCC patients after the different treatments.[39]^3 Therefore, it is urgent for us to further explore the molecular mechanisms underlying the pathophysiology of HCC, which will be of clinical significance to development of novel diagnostic and therapeutic targets for the management of HCC. Recently, the high throughput techniques have been intensively applied in the field of cancer research, and microarray analysis and RNA-seq techniques have been able to generate a large amount of datasets.[40]^4 With the aid of integrated bioinformatics analysis, a large number of novel genes have been identified in various types of cancers.[41]^5 However, the public available datasets including microarray data and RNA-seq data are still not fully explored by the researchers, and the re-analysis of these datasets will enable us to identify novel genes related to the disease of interest. For examples, Ma et al[42]^6 employed the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) for screening differentially expressed genes (DEGs) and identified 11 hub genes associated with HCC progression; Yan et al[43]^7 performed the integrated bioinformatics analysis and suggested that MSTO2P has significant diagnostic and prognostic value in HCC; and Xu et al[44]^8 analyzed the three datasets from the GEO database and found that ten hub genes were associated with the survival of the HCC patients. Due to the large amounts of the datasets in the public database, more hub genes related to the HCC development and progression are yet to be identified. In this study, three GEO datasets including [45]GSE55048, [46]GSE55758, and [47]GSE56545 were used to analyze the differentially expressed genes between non-tumor and HCC groups. The common DEGs from these three datasets were subjected to Gene Ontology (GO), KEGG pathway, Reactome pathway, and protein–protein interaction (PPI) network analysis. A total of 20 hub genes were identified based on the PPI network analysis. Furthermore, the hub genes were further analyzed in the TGCA database to determine the association between the expression of these hub genes and overall survival of the HCC patients. Finally, the biological functions of the potential hub genes were validated by the in vitro functional assays. The integrated bioinformatics analysis of the three GEO datasets may provide us with novel insights into understanding the pathophysiology of HCC. Materials and Methods Collection of Microarray Data The microarray data was downloaded from the GEO database ([48]http://www.ncbi.nlm.nih.gov/geo/) under the accession number [49]GSE55048, [50]GSE55758, and [51]GSE56545. For the [52]GSE55048 dataset, the [53]GPL9115 Illumina Genome Analyzer II platform was used, and a total of four human HCC tissues and four matched non-tumor tissues were included in this dataset; for the [54]GSE55758 dataset, the [55]GPL9052 Illumina Genome Analyzer platform was used, and a total of eight human HCC tissues and eight matched non-tumor tissues were included in this dataset; for the [56]GSE56545 dataset, a total of 12 human HCC tissues and 12 matched non-tumor tissues were included in this dataset. The details of three included datasets are summarized in [57]Table 1. Table 1. Details of the Analyzed GEO Datasets GEO Datasets Platform Samples (Tumor versus Normal) HBV-Positive Tumor Information Country [58]GSE55048 [59]GPL9115 Illumina Genome Analyzer II platform 4 versus 4 4 out of 4 patients Not applicable China [60]GSE55758 [61]GPL9052 Illumina Genome Analyzer 8 versus 8 7 out of 8 patients All patients had a single tumor and most of the primary tumors (5 of 8) had moderately differentiated histology; 6 of 8 had stage II tumor China [62]GSE56545 [63]GPL15433 Illumina HiSeq 1000 12 versus 12 12 out of 12 patients All the patients were male; 3 patients had stage I tumor; 2 patients had stage I–II tumor; 4 patients had stage II tumor; 3 patients had tumor II–III tumor USA [64]Open in a new tab Data Reprocessing and DEG Screening Geo RNA-seq experiments Interactive Navigator (GREIN; [65]http://www.ilincs.org/apps/grein/) was used to retrieve the expression datasets between non-tumor group and HCC group for the [66]GSE55048, [67]GSE55758, and [68]GSE56545 datasets. The DEGs were analyzed using the iGEAK tool based on edgeR. Principle component analysis (PCA) analysis, clustering heatmap plot, volcano plot, and Venn diagram were processed by iGEAK tool. The common DEGs were then exported into excel for further analysis. Functional Analysis of DEGs The ClusterProfiler tool was used to perform Gene Ontology (GO), GO enrichment, KEGG pathway, and Reactome pathway enrichment analysis for the DEGs. GO enrichment analysis including biological processes, cellular components, and molecular function, identified which GO terms were over or underrepresented within a given set of genes.[69]^9 The KEGG knowledge database, an integrated database resource, is generally used to identify functional and metabolic pathways.[70]^9 PPI Network Analysis The Search Tool for the Retrieval of Interacting Genes (STRING)[71]^10 is a biological database for predicting pairs of PPIs. The interactions between DEGs were evaluated using STRING and genes with a combined score >0.7 were defined as key DEGs. Subsequently, Cytoscape (version 3.6.1; [72]https://cytoscape.org/)[73]^11 was used to develop the PPI network of the key DEGs that were identified. Molecular complex detection (MCODE),[74]^12 a Cytoscape plugin, was used with the default parameters to identify the most important modules of the PPI network. Analysis of Hub Gene Expression Level In this study, the Gene Expression Profiling Interactive Analysis (GEPIA) database was sued to visualize the expression of 20 hub genes by using the boxplot. The protein levels of hub genes were analyzed using the Human Protein Atlas database ([75]https://www.proteinatlas.org/). P<0.05 was considered statistically significant. Survival Analysis Survival analysis of the DEGs was performed by using the Gene Expression Profiling Interactive Analysis (GEPIA) database[76]^13 and the Kaplan–Meier (KM) plotter database.[77]^14 The Kaplan–Meier curves were plotted in the GEPIA database including 364 HCC patients. P<0.05 was considered statistically significant. Cell Lines, Cell Culture, Transfection with Small Interfering RNAs, and Drug Treatment LO2, HepG2, SMMC7221, and HCC-LM3 cells were obtained from the Cell Bank of Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences (Shanghai, China). The cells were cultured in the DMEM medium (Thermo Fisher Scientific, Waltham, USA) supplemented with fetal bovine serum (FBS; Thermo Fisher Scientific) and were kept in a humidified incubator at 37°C with 5% CO[2]. The KIF14 and KIF23 siRNAs (si-KIF14 and siKIF-23) and the scrambled control siRNA (si-NC) were obtained from RiboBio (Guangzhou, China). SMMC7721 and HCC-LM3 cells were transfected with siRNAs using Lipofectamine 2000 reagent (Invitrogen, Carlsbad, CA, USA). At 24 hours after transfection, cells were processed for further experimental assays. The cisplatin and sorafenib were purchased from the Sigma (St. Louis, MO, USA). For the cisplatin and sorafenib treatments, the cells were exposed to different concentrations of cisplatin or sorafenib for 48 hours before further experimental assays. RNA Isolation and RT-qPCR The extraction of RNA was performed using the MiniBEST Universal RNA extraction kit (TaKaRa, Dalian, China). A total of 2 μg RNA was reversely transcribed into cDNA using the PrimeScript 1st strand cDNA Synthesis Kit (Takara). The real-time PCR was performed on Applied Biosystems (Applied Biosystems, Foster City, CA, China) using One Step PrimeScript™ RT-PCR Kit (Takara) in accordance with the manufacture’s protocol. GAPDH was used as the reference control for KIF14, KIF23, and Bax expression. The mRNA expression levels were calculated using the 2^−ΔΔCt method. Cell Counting Kit-8 (CCK-8) Assay Cell proliferative ability was measured by the CCK-8 kit (Beyotime, Beijing, China). Briefly, SMMC7721 and HCC-LM3 cells with different treatments were seeded in a 96-well plate and were cultured for different time durations. After that, cells were incubated with 10 µL of CCK-8 reagent at room temperature for 4 hours. The cell proliferation was detected by measuring the absorbance value at 450 nm using a microplate reader. BrdU Incorporation Assay The cell proliferative potential of SMMC7721 and HCC-LM3 cells with different treatments was determined by the BrdU incorporation assay kit (Beyotime) according to the manufacturer’s protocol. Caspase-3/-7 Activity Assay The caspase-3/-7 activity of the HCC cells were determined using the Caspase-Glo assay kit (Promega, Madison, WI, USA). In brief, the HCC cells after different treatments were taken out from the incubator and equilibrated to room temperature for 30 minutes. After that, the cells were incubated with 100 μL of Caspase-Glo reagent for 2 hours at room temperature. The caspase-3/-7 activity of HCC cells was determined by measuring the luminescence using a plate-reading luminometer (Thermo Fisher Scientific). Transwell Invasion Assay The invasive potential of the HCC cells was assessed by the transwell invasion assay. HCC cells after different treatments were seeded on the upper chamber with transwell inserts (Matrigel-coated membranes, pore size: 8 µm; Millipore, Burlington, VT, USA). The upper chamber was filled with culture medium without FBS, and the lower chamber was filled with culture medium with FBS. After 24 hours, the invaded cells were fixed with 4% paraformaldehyde for 30 minutes followed by staining with 0.1% crystal violet solution for 30 minutes, and the invaded cells were counted under a light microscope by randomly choosing five fields. Wound Healing Assay The migratory ability of the HCC cells was determined by the wound healing assay. Briefly, HCC cells after different treatments were seeded onto the 6-well plates. After the cells reached 90% confluence, the wound was created by scratching the cellular monolayer with a sterile pipette tip, and the cell debris was washed away by phosphate buffered saline (PBS). The wound healing was observed for 48 hours. The wound width was determined at 0 and 48 hours, respectively. Statistical Analysis Statistical analysis was performed using GraphPad Prism software V7.0 (GraphPad Software, La Jolla, CA, USA). One-way ANOVA or unpaired Student’s t-test was used for determination of significance of differences among/between different treatment groups. A P-value less than 0.05 was considered statistically significant. Results Analysis of DEGs in [78]GSE55048, [79]GSE55758, and [80]GSE56545 Datasets Firstly, PCA analysis was performed on the expression datasets, and the PCA results showed that the non-tumor group and tumor group were divided into two distinctive groups for these three datasets ([81]Figure 1A). Secondly, the Pearson correlation coefficients between samples were determined by RNA-seq correlation analysis, and the Pearson correlation coefficient of each sample was greater than 0.7, indicating good repeatability ([82]Figure 1B). Furthermore, the DEGs from the three datasets were analyzed using edgeR, and 2,439 DEGs (1,050 genes were up-regulated and 1,389 genes were down-regulated), 2,002 DEGs (1,228 genes were up-regulated and 774 genes were down-regulated), and 3,383 DEGs (1,997 genes were up-regulated and 1,386 genes were down-regulated) were identified in the [83]GSE55048, [84]GSE55758, and [85]GSE56545 datasets, respectively. The DEGs were illustrated in the heatmap clustering and volcano plot ([86]Figure 1C and [87]D). Based on the Veen diagram results, a total of 544 common DEGs were identified in these three datasets ([88]Figure 1E). The gene list of the 544 common DEGs was supplied in [89]Supplemental Table S1. Figure 1. [90]Figure 1 [91]Open in a new tab Analysis of DEGs in [92]GSE55048, [93]GSE55758, and [94]GSE56545 datasets. (A) PCA analysis of [95]GSE55048, [96]GSE55758, and [97]GSE56545 expression datasets. (B) Correlation analysis of [98]GSE55048, [99]GSE55758, and [100]GSE56545 expression datasets. (C) Heatmap clustering of the DEGs between non-tumor and tumor groups in [101]GSE55048, [102]GSE5575, and [103]GSE56545 datasets. (D) Volcano plot of the DEGs between non-tumor and tumor groups in [104]GSE55048, [105]GSE5575, and [106]GSE56545 datasets. (E) Common DEGs shared by the three datasets. GO Analysis, KEGG Pathway, and Reactome Pathway Analysis The GO terms annotation of the common DEGs were shown in [107]Figure 2A, the top terms of the biological process include “cell metabolic process”, “nitrogen compound metabolic process”, “regulation of metabolic process”, etc. ([108]Figure 2A); the top terms of the cellular component include “extracellular region part”, “membrane part”, “plasma membrane”, etc ([109]Figure 2A); the top terms of molecular function include “hydrolase activity”, “DNA-binding transcription factor activity”, “transferase activity”, and so on ([110]Figure 2A). The top GO enrichment results are shown in [111]Figure 2B. In the biological process group, the majority of DEGs were enriched in “organelle fission”, “nuclear division”, “chromosome segregation”, “mitotic nuclear division”, “region of cell cycle phase transition” ([112]Figure 2B); in the cellular component group, the majority of DEGs were enriched in “chromosomal region”, “spindle”, “chromosome, centromeric region”, condensed chromosome” ([113]Figure 2B); in the molecular function group, the majority of the DEGs were enriched in “catalytic activity, acting on DNA”, “tubulin binding”, “ATPasae activity”, “coenzyme binding”, and “helicase activity” ([114]Figure 2B). [115]Figure 2C illustrated the enrichment map of the enriched GO terms in the different groups. Figure 2. [116]Figure 2 [117]Open in a new tab GO annotations and enrichment analysis of the 544 common DEGs identified from [118]GSE55048, [119]GSE55758, and [120]GSE56545 datasets. (A) GO classification of the common DEGs into the biological process, cellular component and molecular function. (B) Top enriched GO terms for common DEGs classified into biological process, cellular component, and molecular function groups. (C) Enrichment map of enriched GO terms for common DEGs. The size of the circle represents the number of genes enriched in the pathway; the color of the circle represents the P-value. The KEGG pathway analysis results are shown in [121]Figure 3A. The top enriched pathways of the DEGs were “Cell cycle”, “oocyte meiosis”, “cellular senescence”, and “DNA replication”. [122]Figure 3B illustrated the enrichment map of the enriched KEGG pathways. Furthermore, the Reactome pathway analysis results are shown in [123]Figure 3C. The top enriched pathways of the DEGs were “cell cycle checkpoints”, “M phase”, “mitotic prometaphase”, mitotic G1-G1/S phases”, “mitotic anaphase” ([124]Figure 3C), and the enrichment map of the enriched Reactome pathways is shown in [125]Figure 3D. Figure 3. [126]Figure 3 [127]Open in a new tab KEGG pathway and Reactome pathway analysis of the 544 common DEGs identified from [128]GSE55048, [129]GSE55758, and [130]GSE56545 datasets. (A) Top enriched KEGG pathways for the 544 DEGs. (B) Enrichment map of enriched KEGG pathways for common DEGs. (C) Top enriched Reactome pathways for the 544 DEGs. (D) Enrichment map of enriched Reactome pathways for the 544 DEGs. The size of the circle represents the number of genes enriched in the pathway; the color of the circle represents the P-value. PPI Network and Identification of Hub Genes The PPI network of the DEGs was constructed using the STRING base, as shown in [131]Figure 4. A total of 530 nodes and 3,086 edges were shown in the PPI network. Subsequently, when “score greater than 5” was defined as the cut-off criterion in MCODE, five clusters were identified from the PPI network, and the most significant clusters consisted of 55 nodes and 167 edges ([132]Figure 5). Furthermore, the top 20 genes (CENPE, DTL, FANCI, MKI67, ZWINT, CENPF, CDC45, ASPM, TPX2, NEK2, KIF14, UBE2C, AURKB, CCNB1, KIF23, RACGAP1, BUB1, HJURP, NCAPG, CDK1) with highest MCODE score were identified as hub genes for subsequent analysis. Figure 4. [133]Figure 4 [134]Open in a new tab PPI network of the common DEGs identified from [135]GSE55048, [136]GSE55758, and [137]GSE56545 datasets. Figure 5. [138]Figure 5 [139]Open in a new tab The most significant cluster identified based on MCODE from the PPI network. Expression Profiles of Hub Genes in the HCC By querying the TGCA database, we evaluated the expression of these hub genes in non-tumor tissues and HCC tissues. As shown in [140]Figure 6, CENPE gene was not significantly up-regulated in the HCC tissues when compared the non-tumor tissues; while all the other hub genes were significantly up-regulated in the HCC tissues when compared to the non-tumor tissues ([141]Figure 6). Figure 6. [142]Figure 6 [143]Open in a new tab Expression of top 20 hub genes in non-tumor tissues and HCC tissues. Red column represents tumor tissues; grey column represents normal tissues; *P<0.05. Abbreviations: num (T), number of tumor tissues; num (N), number of normal tissues. Association Between the Expression of Hub Genes and the Overall Survival of Patients with HCC By querying the TGCA database from the GEPIA tool, we assessed the association between the expression of hub genes and the overall survival of HCC patients. As shown in [144]Figure 7, all the hub genes showed the significant correlation with overall survival of the HCC patients. Furthermore, the analysis using the KM plotter database showed the consistent results ([145]Figure 8). The results indicated that the high expression of these hub genes were significantly associated with the poorer overall survival of the HCC patients. Figure 7. [146]Figure 7 [147]Open in a new tab Association between the expression of hub genes and the overall survival of patients with HCC was analyzed using GEPIA tool. Figure 8. [148]Figure 8 [149]Open in a new tab Association between the expression of hub genes and the overall survival of patients with HCC was analyzed using KM plotter database. KIF14 and KIF23 Knockdown Suppressed the Proliferation of HCC Cells The qRT-PCR analysis showed that the mRNA expression levels of KIF14 and KIF23 were significantly up-regulated in the HCC cell lines (SMMC7721 and HCC-LM3) when compared to LO2 cells ([150]Figure 9A and [151]B). The knockdown of KIF14 and KIF23 in SMMC7721 and HCC-LM3 cells were achieved by transfecting the HCC cells with respective siRNAs targeting KIF14 and KIF23 ([152]Figure 9C and [153]D). The BrdU assay results showed that knockdown of KIF14 and KIF23 both significantly reduced the BrdU-positive SMMC7721 and HCC-LM3 cells ([154]Figure 9E and [155]F). The CCK-8 results showed that KIF14 and KIF23 silence both significantly suppressed the proliferative potential of SMMC7721 and HCC-LM3 cells ([156]Figure 9G–[157]J). In order to assess the effects of KIF14 and KIF23 on the cell apoptosis, we measured the caspase-3/-7 activity and the expression of Bax in the HCC cells. As shown in [158]Figure 9K and [159]L, knockdown of KIF14 and KIF23 both markedly increased the caspase-3/-7 activities of SMMC7721 and HCC-LM3 cells ([160]Figure 9K and [161]L). Consistently, the mRNA expression levels of Bax were elevated in these HCC cells with KIF14 or KIF23 knockdown ([162]Figure 9M and [163]N). Furthermore, the transwell invasion and wound healing assays showed that knockdown of KIF14 and KIF23 both remarkably suppressed the invasive and migratory abilities of SMMC7721 and HCC-LM3 cells ([164]Figure 9O–[165]R). Figure 9. [166]Figure 9 [167]Open in a new tab KIF14 and KIF23 knockdown suppressed the proliferation of HCC cells. (A, B) The mRNA expression of KIF14 (A) and KIF23 (B) in LO2, SMMC7721, and HCC-LM3 was determined by qRT-PCR. (C) The mRNA expression of KIF14 in SMMC7721 and HCC-LM3 cells transfected with si-NC or si-KIF14 was determined by qRT-PCR. (D) The mRNA expression of KIF23 in SMMC7721 and HCC-LM3 cells transfected with si-NC or si-KIF23 was determined by qRT-PCR. (E) The cell proliferation of SMMC7721 and HCC-LM3 cells transfected with si-NC or si-KIF14 was determined by BrdU assay. (F) The cell proliferation of SMMC7721 and HCC-LM3 cells transfected with si-NC or si-KIF23 was determined by BrdU assay. (G, H) The cell viability of the SMMC7721 (G) and HCC-LM3 cells (H) transfected with si-NC or si-KIF14 was determined by CCK-8 assay. (I, J) The cell viability of the SMMC7721 (I) and HCC-LM3 cells (J) transfected with si-NC or si-KIF23 was determined by CCK-8 assay. (K, L) The caspase-3/-7 activities of the SMMC7721 (K) and HCC-LM3 cells (L) were determined by the caspase-3/-7 activity assay kit. (M, N) The mRNA expression levels of Bax in SMMC7721 (M) and HCC-LM3 cells (N) were determined by qRT-PCR. (O, P) The invasive potential of the SMMC7721 (O) and HCC-LM3 cells (P) were measured by transwell invasion assay. (Q, R) The migratory abilities of the SMMC7721 (Q) and HCC-LM3 cells (R) were assessed by transwell invasion assay. N=3; significant difference relative to si-NC groups was shown as *P<0.05, **P<0.01, and ***P<0.001. KIF14 and KIF23 Knockdown Promoted the Chemosensitivity to Cisplatin and Sorafenib The effects of cisplatin and sorafenib on the cell viability of HCC cells were determined by CCK-8 assay. As show in [168]Figure 10A and B, KIF14 and KIF23 knockdown both significantly promoted the sensitivity of HCC cells to cisplatin ([169]Figure 10A and B, [170]Table 2). Consistently, KIF14 and KIF23 silence both also enhanced the sensitivity of HCC cells to sorafenib ([171]Figure 10C and D, [172]Table 2). Figure 10. [173]Figure 10 [174]Open in a new tab KIF14 and KIF23 knockdown promoted the chemosensitivity to cisplatin and sorafenib. (A) The cell viability of the si-NC or si-KIF14-transfected SMMC7721 cells after being treated with different concentrations of cisplatin was determined by CCK-8 assay. (B) The cell viability of the si-NC or si-KIF23-transfected SMMC7721 cells after being treated with different concentrations of cisplatin was determined by CCK-8 assay. (C) The cell viability of the si-NC or si-KIF14-transfected SMMC7721 cells after being treated with different concentrations of sorafenib was determined by CCK-8 assay. (D) The cell viability of the si-NC or si-KIF23-transfected SMMC7721 cells after being treated with different concentrations of sorafenib was determined by CCK-8 assay. N=3. Table 2. IC50 Values of Cisplatin and Sorafenib in SMMC7721 Cells si-NC si-KIF14 si-NC si-KIF23 IC50 of cisplatin (µM) 26.03±2.734 17.75±2.117* 24.20±5.575 13.29±1.560* IC50 of sorafenib (µM) 6.272±1.670 3.240±0.9836* 6.337±0.2505 2.528±1.286* [175]Open in a new tab Notes: Data were presented as mean±standard deviation. N=3; *P<0.05 compared to the si-NC group. The High Expression of KIF14 and KIF23 Predicts Worse Prognosis of Patients with HCC In order to further confirm the prognostic role of KIF14 and KIF23 in the HCC, we analyzed the association between KIF14/KIF23 expression and PFS/RFS/DSS of patients with HCC. As shown in [176]Figure 11A–C, high expression of KIF14 was associated with shorter progression-free survival (PFS; [177]Figure 11A), recurrence-free survival (RFS; [178]Figure 11B), and disease-specific survival (DSS; [179]Figure 11C) of patients with HCC. Consistently, high expression of KIF23 was also associated with shorter PFS ([180]Figure 11D), RFS ([181]Figure 11E), and DSS ([182]Figure 11F) of patients with HCC. In addition, the protein expression levels of KIF14 and KIF23 were significantly higher in the HCC tissues than that in the normal liver tissues, as revealed by the immunohistochemistry from the Human Protein Atlas ([183]Figure 11G and H). Figure 11. [184]Figure 11 [185]Open in a new tab The high expression of KIF14 and KIF23 predicts worse prognosis of patients with HCC. (A–C) The association between the expression of KIF14 and the PFS (A), RFS (B), and DSS (C) of patients with HCC was analyzed using KM plotter database. (D–F) The association between the expression of KIF23 and the PFS (D), RFS (E), and DSS (F) of patients with HCC was analyzed using KM plotter database. (G) The protein expression of KIF14 in normal and HCC tissues was illustrated by immunohistochemistry from Human Protein Atlas. (H) The protein expression of KIF23 in normal and HCC tissues was illustrated by immunohistochemistry from Human Protein Atlas. Discussion The prognosis of HCC patients is still not satisfying due to lack of a sufficient diagnostic tool for early stage HCC and effective therapies for late stage HCC. Because of the complicated biological processes in the development and progression of HCC, the understanding of HCC pathophysiology is still limited. Using the high throughput techniques, researchers have been more efficiently identifying novel targets in the HCC pathophysiology. In this study, we employed three GEO datasets including [186]GSE55048, [187]GSE55758, and [188]GSE5645, and identified the 544 common DEGs in these three datasets. The GO enrichment, KEGG pathway, and Reactome pathway analysis revealed key pathways of these DEGs in the biological functions. Further PPI network analysis identified 20 hub genes from these DEGs. The analysis of the TGCA database showed that the high expression of these hub genes were significantly associated with the poorer overall survival of the HCC patients. The functional assays showed that KIF14 and KIF23 silence attenuated cell proliferation and promoted chemosensitivity of HCC cells. KIF14 and KIF23 may serve as potential biomarkers for predicting the worse prognosis of patients with HCC. The newly identified hub genes in the present study may play an important role in the pathophysiology of HCC. In the present study, by combining these large datasets, we for the first time identified 544 common DEGs, and the GO enrichment analysis showed that the DEGs were enriched in the biological process of nuclear division, mitotic nuclear division, and regulation of cell phase transition, where these biological processes are key for the progression of cancer cells. In the KEGG pathway and Reactome pathway, analysis revealed that the DEGs were enriched in “Cell cycle” and “DNA replication” pathways, suggesting that these DEGs may play a regulatory role in the cell cycle and DNA replication in HCC cells. The PPI network analysis has been commonly used to identify the hub genes from the DEGs. Jin et al[189]^7 used the PPI network and identified sixhub genes associated with HCC progression. Zhou et al[190]^15 identify 15 hub genes related with HCC patients survival by using network analysis; by using PPI network analysis, Teng et al[191]^16 found that ten top hub genes were associated with HCC progression. In the present study, we used the MCODE method from cystoscope and identified 20 hub genes. CENPE functioned as an oncogene and promoted the progression of various types of cancers including lung cancer, esophageal adenocarcinoma, and prostate cancer.[192]^17–19 DTL was found to promote cancer progression by PDCD4 ubiquitin-dependent degradation[193]^20 and DTL depletion inhibited liver cancer cell growth, increased senescence, and reduced tumorigenesis.[194]^21 FANCI was functioned as a DNA-repair protein and has been found to be involved in the breast cancer progression.[195]^22 MKI67 functioned as an oncogene and could promote the HCC cell progression.[196]^23 ZWINT acted as an oncogene in HCC, and overexpression of ZWINT predicted poor prognosis and promoted HCC cell proliferation.[197]^24 For CENPF, frequent amplification of this gene was detected in HCC[198]^25 and CENPF functioned as an oncogene to promote HCC cell proliferation and tube formation.[199]^26 CDC45 was identified as an oncogene in other types of cancers including lung cancer.[200]^27 ASPM has been suggested as a novel marker for vascular invasion, early recurrence, and poor prognosis of hepatocellular carcinoma.[201]^28 TPX2, NEK2, UBE2C, AURKB, RACGAP1, BUB1, HJURP, and NCAPG functioned as oncogenes in HCC.[202]^29–37 CCNB1 is a key regulator in the cell cycle process and can be targeted by oncogenes, which in turn promote HCC cell progression.[203]^38^,[204]^39 Among these hub genes, the kinesin-like protein family has emerged as important regulators in the tumor development and progression.[205]^40–45 In this regard, we further examined the biological actions of KIF14 and KIF23 in HCC cells. KIF14 and KIF23 belong to the kinesin-like protein family, which has been shown to play an important role in cancer biology.[206]^46 The KIF14 protein has a citron kinase binding region, an N-terminal extension for the binding of PCR1 and a C-terminal motor domain, and it functions to modulate midbody formation, cell cytokinesis, mitotic spindle formation and chromosome segregation.[207]^47 KIF14 has been found to be up-regulated in various types of tumors.[208]^47–50 Knockdown of KIF14 has been found to suppress the progression of HCC.[209]^31 KIF23 involves the regulation of cytokinesis and plays a critical role in the process of cytoplasm separation in mitosis.[210]^43 KIF23 has been identified as an oncogene in gastric cancer and glioma,[211]^43^,[212]^51 and recent studies showed that KIF23 variant 1 expression may represent a novel prognostic factor in HCC patients, while the detailed molecular functions of KIF23 remain unknown.[213]^52 In the present study, we found that KIF14 and KIF23 were up-regulated in the HCC cell lines; while KIF14 and KIF23 knockdown suppressed cell proliferation, invasion, and migration, but promoted chemosensitivity to cisplatin and sorafenib in HCC cells. Moreover, mechanistic studies showed that KIF14 and KIF23 knockdown increased caspase-3/-7 activity and Bax mRNA expression level, implying that the effects of KIF14 and KIF23 on the HCC cell proliferation and chemosensitivity may be associated with the enhanced apoptosis. In addition, the high expression of KIF14 and KIF23 was associated with shorter PFS, RFS, and DSS of patients with HCC, suggesting that KIF14 and KIF23 may server as two potential biomarkers for worse prognosis of patients with HCC. The experimental results regarding the role of KIF14 and KIF23 in the pathophysiology of HCC in the present are still at the preliminary stages, which may require further consideration. The results of the study are limited to the in vitro assays, and future studies should examine KIF14- and KIF23-mediated tumor growth and metastasis in vivo. As there are emerging GEO datasets generated by different studies, future studies may consider to analyze more GEO datasets to confirm the prognostic of KIF14 and KIF23 in HCC. In the overall survival analysis using the TCGA cohort, we were unable to access the clinical characteristics (tumor stage, etiology, cirrhosis, etc) of the HCC patients, and, in the future studies, the correlation between the expression of KIF14/23 an the clinical parameters of HCC patients should be determined by using a more detailed TCGA database. More importantly, the up-stream and down-stream mediators of KIF14 and KIF23 have not been determined in our study, which may be considered in our future plans. Conclusion In summary, the present study performed the integrated bioinformatics analysis and identified 20 hub genes that were closely correlated with the overall survival of the HCC patients. In addition, KIF14 and KIF23 silence attenuated cell proliferation, invasion, and migration, and promoted chemosensitivity of HCC cells. KIF14 and KIF23 may serve as potential biomarkers for predicting the worse prognosis of patients with HCC. However, more experimental studies are warranted in the future studies. Funding Statement There is no funding to report. Author Contributions All authors made substantial contributions to the conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; agreed to submit to the current journal; gave final approval of the version to be published; and agree to be accountable for all aspects of the work. Disclosure The authors declare that they have no conflicts of interest. References