Identification of Key Genes and Prognostic Value Analysis in Hepatocellular Carcinoma by Integrated Bioinformatics Analysis

Wang, Meng; Wang, Licheng; Wu, Shusheng; Zhou, Dongsheng; Wang, Xianming

doi:https://doi.org/10.1155/2019/3518378

International Journal of Genomics

On this page

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

Research Article | Open Access

Volume 2019 | Article ID 3518378 | https://doi.org/10.1155/2019/3518378

Identification of Key Genes and Prognostic Value Analysis in Hepatocellular Carcinoma by Integrated Bioinformatics Analysis

Meng Wang,¹Licheng Wang,²Shusheng Wu,³Dongsheng Zhou,^3,4and Xianming Wang^3,4

Academic Editor: Byung-Hoon Jeong

Received02 Jul 2019

Revised07 Aug 2019

Accepted20 Aug 2019

Published22 Nov 2019

Abstract

Emerging evidence indicates that various functional genes with altered expression are involved in the tumor progression of human cancers. This study is aimed at identifying novel key genes that may be used for hepatocellular carcinoma (HCC) diagnosis, prognosis, and targeted therapy. This study included 3 expression profiles (GSE45267, GSE74656, and GSE84402), which were obtained from the Gene Expression Omnibus (GEO). GEO2R was used to analyze the differentially expressed genes (DEGs) between HCC and normal samples. The functional and pathway enrichment analysis was performed by the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network of the identified DEGs was constructed using the Search Tool for the Retrieval of Interacting Gene, and hub genes were identified. ONCOMINE and CCLE databases were used to verify the expression of the hub genes in HCC tissues and cells. Kaplan-Meier plotter was used to assess the effects of the hub genes on the overall survival of HCC patients. A total of 99 DEGs were identified from the 3 expression profiles. These DEGs were enriched with functional processes and pathways related to HCC pathogenesis. From the PPI network, 5 hub genes were identified. The expression of the 5 hub genes was all upregulated in HCC tissues and cells compared with the control tissues and cells. Kaplan-Meier survival curves indicated that high expression of cyclin-dependent kinase (CDK1), cyclin B1 (CCNB1), cyclin B2 (CCNB2), MAD2 mitotic arrest deficient-like 1 (MAD2L1), and topoisomerase IIα (TOP2A) predicted poor overall survival in HCC patients (all log-rank ). These results revealed that the DEGs may serve as candidate key genes during HCC pathogenesis. The 5 hub genes, including CDK1, CCNB1, CCNB2, MAD2L1, and TOP2A, may serve as promising prognostic biomarkers in HCC.

1. Introduction

Hepatocellular carcinoma (HCC) remains a serious health burden and is the second leading cause of cancer mortality worldwide [1]. Statistical data have indicated that the morbidity and mortality of HCC have been increasing in recent years, mainly due to the increased infection of hepatitis C virus [2]. Researchers have identified several established risk factors for the occurrence of HCC, such as liver cirrhosis, viral infection, metabolic disorder, and heavy alcohol consumption [3]. Despite advances in various therapeutic strategies, such as surgery, chemotherapy, radiotherapy, and biologics, the prognosis and outcomes remain poor in patients suffering from HCC [4]. Therefore, efficient diagnosis and prognosis remain great challenges for HCC treatment.

It is generally considered that tumorigenesis is a complex process with a wide spectrum of genetic alterations [5]. These genes typically exhibit aberrant expression patterns and have clinical significance in cancer diagnosis and prognosis [6]. Currently, some molecules have been recognized as diagnostic and prognostic biomarkers in HCC. For example, the high expression of peroxiredoxin 1 (Prdx1) is associated with tumor development and overall survival of HCC patients and serves as a candidate biomarker for the screening and prediction of this malignancy [7]. Sulfite oxidase (SUOX) expression is downregulated during tumorigenesis of HCC and is correlated with HCC diagnosis and prognosis [8]. Upregulated expression of distal-less homeobox gene 4 (DLX4) in HCC samples has been shown to be associated with poor prognosis of HCC patients [9]. Similarly, the altered alpha-fucosidase (AFU) expression has significant prognostic value in HCC patients and acts as a potential target for HCC-targeted therapy [10]. However, the available biomarkers are not suitable for all the HCC cases due to the limitations of sensitivity and specificity. Accordingly, the identification of novel functional genes may contribute to the understanding of tumor pathogenesis and the improvement of diagnosis and prognosis of HCC.

In recent research, differentially expressed genes (DEGs) in tumor samples compared with normal samples can be identified using gene expression profiling arrays [11, 12]. Some key molecules have also been reported in HCC using bioinformatics analysis [13, 14]. However, the number of the identified functional genes is far from sufficient to explain the mechanisms underlying the pathogenesis of HCC. Thus, this study used bioinformatics analyses to further identify key genes in HCC progression from 3 gene expression profiles from the Gene Expression Omnibus (GEO) database and assessed the clinical significance of the DEGs in HCC prognosis. The expression and prognostic value of the identified key genes were further verified using the data from The Cancer Genome Atlas (TCGA) database.

2. Materials and Methods

2.1. Data Collection

In this study, we firstly downloaded 3 gene expression profiles from GEO database (http://www.ncbi.nlm.nih.gov/geo), including GSE45267, GSE74656, and GSE84402. The inclusion criteria for the expression profiles were as follows: (1) the samples detected are tissues, (2) all tissues are diagnosed with HCC tissues and normal tissues, (3) gene expression profiling of mRNA, (4) samples collected from the same racial population, (4) probes can be converted into the corresponding gene symbols, and (5) complete information for our analyses. The array data of GSE45267 included 49 HCC tumor tissues and 38 normal tissues. GSE74656 contained 10 samples, including 5 HCC tumors and 5 adjacent normal tissues. GSE84402 was comprised of 14 tumor tissues and 14 adjacent noncancerous tissues [15].

2.2. Data Processing

The DEGs between the HCC samples and normal samples were analyzed by GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r), which is a built-in online tool of GEO [16]. Adjusted value and |log fold change| (|log FC|) were used to evaluate the significance of DEGs, and adjusted and were set as the cutoff criteria.

2.3. Functional and Pathway Enrichment Analysis

The Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.ncifcrf.gov/) is an essential program for the comprehensive gene function analysis, which aids the researchers to understand the biological significance of abundant genes [17]. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway enrichment analysis were performed for the obtained DEGs. A result with a was considered statistically significant.

2.4. PPI Network Construction and Module Selection

Since the interactions between proteins represent the pivotal events during cellular biological processes, we constructed a protein-protein interaction (PPI) network of the identified DEGs using the Search Tool for the Retrieval of Interacting Gene (STRING, http://string.embl.de/) database [18]. The PPI network was visualized using Cytoscape (version 3.7.0) [19], and a was used as the cutoff criterion. Subsequently, the modules of the PPI network were screened by the Molecular Complex Detection (MCODE) with the following parameters: degree , node score , , and maximum [20].

2.5. Expression Analysis

The mRNA expression levels of the hub genes between HCC tissues and normal controls were analyzed using the ONCOMINE (http://www.oncomine.org) database [21], and the data were collected from three literatures [22–24]. In addition, the expression results were further confirmed in HCC cells by the CCLE (http://portals.broadinstitute.org/ccle/home) database [25]. The differences between two groups were analyzed using Student’s -test. and were set as the cutoff criteria.

2.6. Survival Analysis of DEGs

The prognostic value of the identified hub genes in HCC was further assessed by the Kaplan-Meier plotter (KM plotter, http://www.kmplot.com/analysis/) [26]. The analysis included 364 patients, and their KM survival curves were conducted. In addition, the KM survival curves for patients with different tumor stages were separately plotted. However, only 5 patients were diagnosed with tumor stage 4, and the curve could not be performed due to the limited sample size. The actual number in the other stages can be lower due to missing expression values and/or incomplete survival data. The gene expression was grouped using a cutoff value that is located between the lower and upper quartiles and computed by the Kaplan-Meier plotter with a best performing threshold. The log-rank value for the different survival distribution between the low and high expression group was assessed, and the hazard radio (HR) with 95% confidence interval (95% CI) was calculated and plotted on the webpage.

2.7. Verification of Expression and Prognostic Value of DEGs Using TCGA Data

To confirm the clinical significance of the 5 hub genes in the prognosis of HCC, data from TCGA database were further assessed using the Gene Expression Profiling Interactive Analysis (GEPIA), which is a web-based tool to deliver fast and customizable functionalities based on TCGA data [27]. The expression patterns in HCC tissues and the Kaplan-Meier survival curves were all performed using the GEPIA.

3. Results

3.1. Identification of DEGs in HCC

According to the GEO2R analysis, a total of 352, 249, and 455 DEGs were, respectively, identified in GSE45267, GSE74656, and GSE84402. Among these DEGs, 99 genes with significant aberrant expression were extracted from all the three datasets (Figure 1), including 38 upregulated genes and 61 downregulated genes (Table 1).

3.2. GO Analysis and Pathway Enrichment Analysis of DEGs in HCC

The potential biological function of the identified 99 DEGs was assessed using GO analysis. As shown in Table 2, these genes were mainly enriched in biological processes related to cell division and mitotic nuclear division. Moreover, the potential signaling pathways which these DEGs involved were examined using KEGG analysis. From the results in Table 2, we found that the DEGs were mostly enriched in cell cycle and mineral absorption processes.

3.3. PPI Network Construction and Significant Module Analysis

The 99 identified DEGs were all filtered into the PPI network complex, which contained 99 nodes and 298 edges (Figure 2(a)). Furthermore, the most significant module was extracted from the PI network, containing 32 nodes and 78 edges (Figure 2(b)). In this module, 5 nodes with a were identified as hub genes, including cyclin-dependent kinase (CDK1), cyclin B1 (CCNB1), cyclin B2 (CCNB2), MAD2 mitotic arrest deficient-like 1 (MAD2L1), and topoisomerase IIα (TOP2A) (Table 3).

(a)

(b)

3.4. CDK1 Expression Validation and Prognostic Value in HCC

To further confirm the expression patterns of the 5 hub genes in HCC, we obtained 4 datasets from the ONCOMINE database to analyze the differential expression between HCC tissues and normal tissues. As shown in Figures 3(a)–3(d), the expression of CDK1 was significantly upregulated in HCC tissues compared with the normal controls in each dataset (all ), and this difference was also statistically significant combined with the 4 datasets (, Figure 3(e)). Additionally, the mRNA expression of CDK1 in HCC cells was also analyzed using the CCLE database. The results shown in Figure 3(f) revealed that the CDK1 expression was also elevated in HCC cells. Furthermore, the Kaplan-Meier survival curves were constructed based on CDK1 expression in HCC patients. As shown in Figure 3(g), we considered that the high CDK1 expression was associated with poor overall survival compared with the low CDK1 expression in HCC patients (log-rank ). In addition, survival curves for HCC patients with different tumor stages were also plotted, which showed that the high CDK1 predicts poor overall survival in patients with tumor stage 2 (log-rank , Figure 3(i)) and tumor stage 3 (log-rank , Figure 3(j)). However, no significantly different survival times were observed between patients with high CDK1 expression levels and patients with low CDK1 expression levels at tumor stage 1 (log-rank , Figure 3(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 3

Expression and prognostic value of CDK1 in HCC. (a–e) The expression data collected from 4 datasets from ONCOMINE indicated that CDK1 expression was upregulated in HCC tissues compared with the normal control tissues and cells (all ). (f) Expression of CDK1 was increased in HCC cells based on the data from CCLE. (g) The Kaplan-Meier survival curves revealed that the high CDK1 expression predicted worse overall survival compared with the low CDK1 expression in HCC patients (log-rank ). (h) No significantly different survival times were observed between patients with high CDK1 expression and patients with low CDK1 expression at tumor stage 1 (log-rank ). (i) High CDK1 expression predicted worse overall survival compared with the low CDK1 expression in HCC at tumor stage 2 (log-rank ). (j) High CDK1 expression predicted worse overall survival compared with the low CDK1 expression in HCC patients at tumor stage 3 (log-rank ).

3.5. CCNB1 Expression Validation and Prognostic Value in HCC

According to the expression investigation, we observed that the expression of CCNB1 was upregulated in both tissues and cells compared with the control tissues and cells (all , Figures 4(a)–4(e)). Furthermore, the prognostic value of CCNB1 was examined using the KM plotter. As shown in Figure 4(f), we considered that patients with the high CCNB1 expression level had poor overall survival compared with those with low CCNB1 expression level (log-rank ). Moreover, the survival analysis for patients with different tumor stage revealed that the high CCNB1 expression level was associated with shorter survival time compared with the low CCNB1 expression level in HCC patients with tumor stage 1 (log-rank , Figure 4(g)), tumor stage 2 (log-rank , Figure 4(h)), and tumor stage 3 (log-rank , Figure 4(i)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 4

Expression and prognostic value of CCNB1 in HCC. (a–d) The expression data collected from 3 datasets from ONCOMINE indicated that CCNB1 expression was upregulated in HCC tissues compared with the normal controls (all ). (e) Expression of CCNB1 was increased in HCC cells based on the data from CCLE. (f) The Kaplan-Meier survival curves revealed that the high CCNB1 expression predicted worse overall survival compared with the low CCNB1 expression in HCC patients (log-rank ). (g) High CCNB1 expression predicted worse overall survival compared with the low CCNB1 expression in HCC patients at tumor stage 1 (log-rank ). (h) High CCNB1 expression predicted worse overall survival compared with the low CCNB1 expression in HCC patients at tumor stage 2 (log-rank ). (i) High CCNB1 expression predicted worse overall survival compared with the low CCNB1 expression in HCC patients at tumor stage 3 (log-rank ).

3.6. CCNB2 Expression Validation and Prognostic Value in HCC

The expression patterns of CCNB2 in both HCC tissues and cells were greater than those in the normal control tissues and cells (all , Figures 5(a)–5(e)). To investigate the clinical significance of CCNB2 in HCC prognosis, the KM plotter was used to plot the survival curves for HCC patients. Patients with low CCNB2 expression had longer survival times compared with those with high CCNB2 expression (log-rank , Figure 5(f)). Additionally, the effect of CCNB2 on overall survival of patients with different tumor stages were also assessed. The curves indicated that the high CCNB2 expression was associated with shorter survival times compared with the low CCNB2 expression in HCC patients with tumor stage 2 (log-rank , Figure 5(h)) and tumor stage 3 (log-rank , Figure 5(i)). However, no significantly different survival times were observed between patients with high CCNB2 expression levels and patients with low CCNB2 expression levels at tumor stage 1 (log-rank , Figure 5(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 5

Expression and prognostic value of CCNB2 in HCC. (a–d) The expression data collected from 3 datasets from ONCOMINE indicated that CCNB2 expression was upregulated in HCC tissues compared with the normal controls (all ). (e) Expression of CCNB2 was increased in HCC cells based on the data from CCLE. (f) The Kaplan-Meier survival curves revealed that the high CCNB2 expression predicted worse overall survival compared with the low CCNB2 expression in HCC patients (log-rank ). (g) No significantly different survival times were observed between patients with high CDK1 expression and patients with low CDK1 expression at tumor stage 1 (log-rank ). (h) High CCNB2 expression predicted worse overall survival compared with the low CCNB2 expression in HCC patients at tumor stage 2 (log-rank ). (i) High CCNB2 expression predicted worse overall survival compared with the low CCNB2 expression in HCC patients at tumor stage 3 (log-rank ).

3.7. MAD2L1 Expression Validation and Prognostic Value in HCC

The expression of MAD2L1 was analyzed using the ONCOMINE database and the CCLE database and was proved to be upregulated in HCC tissues and cells compared with the control tissues and cells (all , Figures 6(a)–6(f)). Furthermore, the KM plotter was used to plot survival curves based on MAD2L1 expression in HCC patients. As shown in Figure 6(g), HCC patients with high MAD2L1 expression had poor overall survival compared with those with low MAD2L1 expression (log-rank ). To explore the effect of MAD2L1 expression on HCC tumors with different tumor stages, survival analysis was performed for patients with tumor stages 1-3. The results indicated that the high MAD2L1 expression predicted shorter survival times compared with the low MAD2L1 expression in HCC patients with tumor stage 1 (log-rank , Figure 6(h)), tumor stage 2 (log-rank , Figure 6(i)), and tumor stage 3 (log-rank , Figure 6(j)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 6

Expression and prognostic value of MAD2L1 in HCC. (a–e) The expression data collected from 4 datasets from ONCOMINE indicated that MAD2L1 expression was upregulated in HCC tissues compared with the normal controls (all ). (f) Expression of MAD2L1 was increased in HCC cells based on the data from CCLE. (g) The Kaplan-Meier survival curves revealed that the high MAD2L1 expression predicted worse overall survival compared with the low MAD2L1 expression in HCC patients (log-rank ). (h) High MAD2L1 expression predicted worse overall survival compared with the low MAD2L1 expression in patients with HCC at tumor stage 1 (log-rank ). (i) High MAD2L1 expression predicted worse overall survival compared with the low MAD2L1 expression in HCC patients at tumor stage 2 (log-rank ). (j) High MAD2L1 expression predicted worse overall survival compared with the low MAD2L1 expression in patients with HCC at tumor stage 3 (log-rank ).

3.8. TOP2A Expression Validation and Prognostic Value in HCC

The expression of TOP2A was investigated using the ONCOMINE database and the CCLE database, and the expression level was observed to be elevated in both HCC tissues and cells compared with the normal control tissues and cells (all , Figures 7(a)–7(f)). The survival analysis indicated that HCC patients with high TOP2A expression had poor overall survival compared with those with low TOP2A expression (log-rank , Figure 7(g)). Moreover, we also found that the high TOP2A expression predicted poor overall survival compared with the low TOP2A expression in HCC patients with tumor stage 2 (log-rank , Figure 7(i)) and tumor stage 3 (log-rank , Figure 7(j)). In HCC patients with tumor stage 1, we observed that the trend of the high TOP2A expression was related to a shorter survival time, but the difference was not statistically significant (log-rank , Figure 7(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 7

Expression and prognostic value of TOP2A in HCC. (a–e) The expression data collected from 4 datasets from ONCOMINE indicated that TOP2A expression was upregulated in HCC tissues compared with the normal controls (all ). (f) Expression of TOP2A was increased in HCC cells based on the data from CCLE. (g) The Kaplan-Meier survival curves revealed that the high TOP2A expression predicted worse overall survival compared with the low TOP2A expression in HCC patients (log-rank ). (h) No significantly different survival times were found between patients with high TOP2A and patients with low TOP2A at tumor stage 1 (log-rank ). (i) High TOP2A expression predicted worse overall survival compared with the low TOP2A expression in HCC patients at tumor stage 2 (log-rank ). (j) High TOP2A expression predicted worse overall survival compared with the low TOP2A expression in patients with HCC at tumor stage 3 (log-rank ).

3.9. Expression and Prognostic Value Verification Using TCGA Data

By using the GEPIA, the expression patterns and prognostic value of the 5 hub genes were verified based on the data from TCGA database. Consistent with the expression results analyzed by ONCOMINE, the expression levels of CDK1, CCNB1, CCNB2, MAD2L1, and TOP2A assessed by TCGA data were all upregulated in tumor tissues compared with the normal controls (all , Figure 8(a)). Furthermore, the survival curves shown in Figure 8(b) indicated that high CDK1, CCNB1, MAD2L1, and TOP2A expression predicted poor overall survival (all log-rank ). Although high expression of CCNB2 was also associated with shorter survival time, the difference of survival distribution between high and low expression groups was not statistically significant (log-rank ).

(a)

(b)

4. Discussion

Accurate diagnosis and prognosis remain the great challenges for the improvement of HCC outcomes. To meet the clinical requirements of HCC treatment, various therapeutic methods have been developed in recent decades [28]. Moreover, targeted therapy, which is mainly dependent on genes that have pivotal roles during tumor pathogenesis, has attracted increasing attention [29]. These key genes are involved in tumor progression and typically have considerable clinical significance in the diagnosis and prognosis of various human cancers, including HCC. Ba and colleagues have reported that the serum expression of Golgi protein-73 (GP73) was higher in HCC patients compared with healthy individuals and serves as a novel biomarker for HCC diagnosis [30]. Similarly, upregulated expression of lysine specific demethylase 1 (LSD1) has been proven to be associated with poor prognosis of HCC [31]. To identify novel key genes that might be involved in HCC pathogenesis, we performed a systematic analysis of 3 expression profiles from GEO database using bioinformatics analysis. The DEGs in the expression profiles and the prognostic value of the key genes were assessed in the present study.

A total of 68 HCC samples and 57 normal control samples were included in the 3 expression profiles, and 99 DEGs were screened for further analyses. According to the functional and pathway enrichment analysis, the identified DEGs were shown to be enriched in biological processes that related to cell division and mitotic nuclear division and in signaling pathways that associated with cell cycle and mineral absorption. It is generally considered that cell division, mitotic nuclear division, and cell cycle are important cell processes in both normal and tumor cells [32]. Tumor-related key genes are typically involved in tumor progression by the regulation of these cell processes [33, 34].

Two interesting results were presented in our study. Firstly, the 99 DEGs included 6 members of the cytochrome P450 proteins (CYPs), including CYP1A2, CYP26A1, CYP2B6, CYP2C9, CYP2C19, and CYP2C18. CYPs represent a large group of enzymes with critical roles in the molecular metabolism [35]. They act critical roles in the development of various human cancers, including HCC, and mediate the metabolism of most of the procarcinogens [36]. The members of CYPs in our study were found to be downregulated in HCC tissues, indicating that the CYPs were involved in the progression of HCC and might inhibit the drug sensitivity. Secondly, the DEGs were found to be enriched in the mineral absorption pathway in this study. A previous research has demonstrated that mineral supplementation could improve the status of essential trace elements in biological samples collected from patients with liver cirrhosis and cancer [37]. Therefore, we speculated that the mineral absorption pathway might have effects on the maintenance of HCC tumor microenvironment, which needs to be analyzed and confirmed in future studies.

Furthermore, the PPI network of the DEGs was constructed, and 5 hub genes were extracted from a significant module, including CDK1, CCNB1, CCNB2, MAD2L1, and TOP2A. A study by Xing et al. [38] also focused on the DEGs in HCC tissues compared with the normal controls, and a same expression profile GSE45267 was analyzed and CCNB2 and TOP2A were identified as two of the hub genes, which was consistent with our corresponding data. Furthermore, the expression patterns of the 5 hub genes were found all upregulated in HCC tissues and cells compared with the normal controls, and their prognostic value was evaluated by plotting the Kaplan-Meier survival curves. The analysis results indicated that the expression levels of CDK1, CCNB1, CCNB2, MAD2L1, and TOP2A were associated with the overall survival of HCC patients. Additionally, the relationships between the hub genes and overall survival of HCC cases at different tumor stages were also observed, suggesting that these genes might serve as promising prognostic biomarkers in HCC.

CDK1 belongs to a serine/threonine kinase family and serves as a critical cell cycle-regulating protein. It has been widely investigated in human malignancies and has been found to be involved in tumor progression. In epithelial ovarian cancer, upregulated expression of CDK1 has been observed in cancer cells and promotes cancer growth and has a significant effect on the overall survival of patients [39]. In addition, a study scheduled by Luo et al. revealed that CDK1 had comprehensive effects on gene interaction networks in the tumor progression of cervical cancer and thus indicated the potential role of CDK1 as a therapeutic target [40]. In HCC, the aberrant CDK1 expression could regulate the apoptin-induced apoptosis with a pivotal role in tumor progression [41]. We also found the increased CDK1 expression in HCC samples and proved its prognostic value for cancer patients. The molecular mechanisms underlying the role of CDK1 in human cancers await more research.

CCNB1 and CCNB2 are two important cyclins that are closely correlated with the cell cycle and cell growth. Overexpression of CCNB1 and CCNB2 has been observed in some human cancer samples, and CCNB1 and CCNB2 possess clinical significance in the diagnosis and prognosis of various cancers, such as lung cancer [42] and pancreatic cancer [43]. Our study also showed the upregulated expression of CCNB1 and CCNB2 in both HCC tissues and cells and reported their prognostic value for the patients. However, the clinical significance verification using TCGA data showed that the difference between the survival distributes of low and high CCNB2 expression groups was not statistically significant, which might be due to the limited sample size and the incomplete survival information. Thus, although CCNB1 and CCNB2 have been previously reported to act as therapeutic target genes in HCC [44, 45], the clinical significance of CCNB1 and CCNB2 needs to be investigated in cancer-related research.

MAD2L1 plays an important role in spindle checkpoints during mitosis. Dysregulation of MAD2L1 induces the instability of chromosomes and chromosomal aneuploidy, which are common events in cancer [46]. It has been determined as a useful prognostic biomarker in some cancers, such as breast cancer [47] and lung adenocarcinoma [48]. An increased expression level of MAD2L1 was observed in HCC samples in the present study, which was consistent with the results from a study by Li et al. [49], which also found the overexpression of MAD2L1 in HCC. Collectively, the role of MAD2L1 in cancer pathogenesis and the related molecular mechanisms need to be assessed with in-depth studies.

TOP2A is an enzyme that is closely correlated with DNA replication, recombination, transcription, and chromatin remodeling [50]. The functional and clinical roles of TOP2A have been demonstrated in human cancers, including prostate cancer [51], breast cancer [52], and nasopharyngeal carcinoma [53]. This study showed elevated TOP2A expression levels in HCC tissues and cells and demonstrated its potential as a prognostic biomarker in this malignancy. The upregulation of TOP2A and its prognostic value have been reported in HCC in previous studies, which also revealed its correlation with tumor onset and chemoresistance [54]. Further studies should be carried out to explore the mechanisms underlying the role of TOP2A during cancer pathogenesis.

5. Conclusion

In conclusion, this study identified 99 DEGs from 3 expression profiles by integrated bioinformatics analysis. These DEGs may contain key genes involved in HCC pathogenesis. In addition, CDK1, CCNB1, CCNB2, MAD2L1, and TOP2A were the top five hub genes and serve as candidate prognostic biomarkers in HCC. The results of this study further enrich the number of key genes that may be involved in the pathogenesis of HCC and give in silico evidence for the key genes in the prognosis of HCC. However, our study fails to evaluate the clinical significance and biological function of the key genes in tumor samples by in vitro and in vivo analyses. Thus, further studies are needed to confirm the prognostic value and functional roles of these key genes in HCC.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Meng Wang, Licheng Wang, and Shusheng Wu contributed equally to this work.

Acknowledgments

This study was funded by the National Natural Science Foundation (81702752) and the Shandong Province Natural Science Foundation (ZR2017LH008).

References

L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, J. Lortet-Tieulent, and A. Jemal, “Global cancer statistics, 2012,” CA: A Cancer Journal for Clinicians, vol. 65, no. 2, pp. 87–108, 2015.
View at: Publisher Site | Google Scholar
B. Njei, Y. Rotman, I. Ditah, and J. K. Lim, “Emerging trends in hepatocellular carcinoma incidence and mortality,” Hepatology, vol. 61, no. 1, pp. 191–199, 2015.
View at: Publisher Site | Google Scholar
D. L. White, F. Kanwal, and H. B. El-Serag, “Association between nonalcoholic fatty liver disease and risk for hepatocellular cancer, based on systematic review,” Clinical Gastroenterology and Hepatology, vol. 10, no. 12, pp. 1342–1359.e2, 2012.
View at: Publisher Site | Google Scholar
A. A. Ashhab, H. Rodin, J. Powell, and J. D. Debes, “Hepatocellular carcinoma diagnosis and surveillance: socioeconomic factors don’t seem to matter, unless you are an immigrant,” Journal of Hepatology, vol. 67, no. 3, pp. 648-649, 2017.
View at: Publisher Site | Google Scholar
M. Y. Huang, H. C. Chen, I. P. Yang et al., “Tumorigenesis and tumor progression related gene expression profiles in colorectal cancer,” Cancer Biomarkers, vol. 13, no. 4, pp. 269–279, 2013.
View at: Publisher Site | Google Scholar
F. J. Xie, H. Y. Lu, Q. Q. Zheng et al., “The clinical pathological characteristics and prognosis of FGFR1 gene amplification in non-small-cell lung cancer: a meta-analysis,” OncoTargets and Therapy, vol. 9, pp. 171–181, 2016.
View at: Publisher Site | Google Scholar
Q. K. Sun, J. Y. Zhu, W. Wang et al., “Diagnostic and prognostic significance of peroxiredoxin 1 expression in human hepatocellular carcinoma,” Medical Oncology, vol. 31, no. 1, p. 786, 2014.
View at: Publisher Site | Google Scholar
G. Z. Jin, W. L. Yu, H. Dong et al., “SUOX is a promising diagnostic and prognostic biomarker for hepatocellular carcinoma,” Journal of Hepatology, vol. 59, no. 3, pp. 510–517, 2013.
View at: Publisher Site | Google Scholar
Y. Gao, Z. Li, X. Guo, Y. Liu, and K. Zhang, “DLX4 as a prognostic marker for hepatocellular carcinoma,” Neoplasma, vol. 61, no. 3, pp. 318–323, 2014.
View at: Google Scholar
K. Wang, W. Guo, N. Li et al., “Alpha-1-fucosidase as a prognostic indicator for hepatocellular carcinoma following hepatectomy: a large-scale, long-term study,” British Journal of Cancer, vol. 110, no. 7, pp. 1811–1819, 2014.
View at: Publisher Site | Google Scholar
W. Sui, Z. Shi, W. Xue et al., “Circular RNA and gene expression profiles in gastric cancer based on microarray chip technology,” Oncology Reports, vol. 37, no. 3, pp. 1804–1814, 2017.
View at: Publisher Site | Google Scholar
J. C. Tantai, X. F. Pan, and H. Zhao, “Network analysis of differentially expressed genes reveals key genes in small cell lung cancer,” European Review for Medical and Pharmacological Sciences, vol. 19, no. 8, pp. 1364–1372, 2015.
View at: Google Scholar
C. Zhang, L. Peng, Y. Zhang et al., “The identification of key genes and pathways in hepatocellular carcinoma by bioinformatics analysis of high-throughput data,” Medical Oncology, vol. 34, no. 6, p. 101, 2017.
View at: Publisher Site | Google Scholar
T. Mou, D. Zhu, X. Wei et al., “Identification and interaction analysis of key genes and microRNAs in hepatocellular carcinoma by bioinformatics analysis,” World Journal of Surgical Oncology, vol. 15, no. 1, p. 63, 2017.
View at: Publisher Site | Google Scholar
H. Wang, X. Huo, X. R. Yang et al., “STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4,” Molecular Cancer, vol. 16, no. 1, p. 136, 2017.
View at: Publisher Site | Google Scholar
T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, pp. D991–D995, 2013.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, pp. D447–D452, 2015.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
G. D. Bader and C. W. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, p. 2, 2003.
View at: Publisher Site | Google Scholar
D. R. Rhodes, J. Yu, K. Shanker et al., “ONCOMINE: A Cancer Microarray Database and Integrated Data-Mining Platform,” Neoplasia, vol. 6, no. 1, pp. 1–6, 2004.
View at: Publisher Site | Google Scholar
X. Chen, S. T. Cheung, S. So et al., “Gene expression patterns in human liver cancers,” Molecular Biology of the Cell, vol. 13, no. 6, pp. 1929–1939, 2002.
View at: Publisher Site | Google Scholar
S. Roessler, H. L. Jia, A. Budhu et al., “A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients,” Cancer Research, vol. 70, no. 24, pp. 10202–10212, 2010.
View at: Publisher Site | Google Scholar
E. Wurmbach, Y. B. Chen, G. Khitrov et al., “Genome‐wide molecular profiles of HCV‐induced dysplasia and hepatocellular carcinoma,” Hepatology, vol. 45, no. 4, pp. 938–947, 2007.
View at: Publisher Site | Google Scholar
R. Petryszak, M. Keays, Y. A. Tang et al., “Expression Atlas update—an integrated database of gene and protein expression in humans, animals and plants,” Nucleic Acids Research, vol. 44, no. D1, pp. D746–D752, 2016.
View at: Publisher Site | Google Scholar
B. Győrffy, P. Surowiak, J. Budczies, and A. Lánczky, “Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer,” PLoS One, vol. 8, no. 12, article e82241, 2013.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017.
View at: Publisher Site | Google Scholar
D. Graf, D. Vallböhmer, W. T. Knoefel et al., “Multimodal treatment of hepatocellular carcinoma,” European Journal of Internal Medicine, vol. 25, no. 5, pp. 430–437, 2014.
View at: Publisher Site | Google Scholar
A. L. Fonseca, V. L. da Silva, M. M. da Fonsêca et al., “Bioinformatics analysis of the human surfaceome reveals new targets for a variety of tumor types,” International Journal of Genomics, vol. 2016, Article ID 8346198, 7 pages, 2016.
View at: Publisher Site | Google Scholar
M. C. Ba, H. Long, Y. Q. Tang, and S. Z. Cui, “GP73 expression and its significance in the diagnosis of hepatocellular carcinoma: a review,” International Journal of Clinical and Experimental Pathology, vol. 5, no. 9, pp. 874–881, 2012.
View at: Google Scholar
Z. K. Zhao, H. F. Yu, D. R. Wang et al., “Overexpression of lysine specific demethylase 1 predicts worse prognosis in primary hepatocellular carcinoma patients,” World Journal of Gastroenterology, vol. 18, no. 45, pp. 6651–6656, 2012.
View at: Publisher Site | Google Scholar
S. Y. Hyun, E. M. Rosen, and Y. J. Jang, “Novel DNA damage checkpoint in mitosis: mitotic DNA damage induces re- replication without cell division in various cancer cells,” Biochemical and Biophysical Research Communications, vol. 423, no. 3, pp. 593–599, 2012.
View at: Publisher Site | Google Scholar
H. O. Kaseb, D. W. Lewis, W. S. Saunders, and S. M. Gollin, “Cell division patterns and chromosomal segregation defects in oral cancer stem cells,” Genes, Chromosomes & Cancer, vol. 55, no. 9, pp. 694–709, 2016.
View at: Publisher Site | Google Scholar
X. Wang, Q. Yu, Y. Zhang, Z. Ling, and P. Yu, “Tectonic 1 accelerates gastric cancer cell proliferation and cell cycle progression in vitro,” Molecular Medicine Reports, vol. 12, no. 4, pp. 5897–5902, 2015.
View at: Publisher Site | Google Scholar
E. A. Korobkova, “Effect of natural polyphenols on CYP metabolism: implications for diseases,” Chemical Research in Toxicology, vol. 28, no. 7, pp. 1359–1390, 2015.
View at: Publisher Site | Google Scholar
J. Zhou, Q. Wen, S. F. Li et al., “Significant change of cytochrome P450s activities in patients with hepatocellular carcinoma,” Oncotarget, vol. 7, no. 31, pp. 50612–50623, 2016.
View at: Publisher Site | Google Scholar
T. G. Kazi, N. F. Kolachi, H. I. Afridi et al., “Effects of mineral supplementation on liver cirrhotic/cancer male patients,” Biological Trace Element Research, vol. 150, no. 1-3, pp. 81–90, 2012.
View at: Publisher Site | Google Scholar
T. Xing, T. Yan, and Q. Zhou, “Identification of key candidate genes and pathways in hepatocellular carcinoma by integrated bioinformatical analysis,” Experimental and Therapeutic Medicine, vol. 15, no. 6, pp. 4932–4942, 2018.
View at: Publisher Site | Google Scholar
W. Yang, H. Cho, H. Y. Shin et al., “Accumulation of cytoplasmic Cdk1 is associated with cancer growth and survival rate in epithelial ovarian cancer,” Oncotarget, vol. 7, no. 31, pp. 49481–49497, 2016.
View at: Publisher Site | Google Scholar
Y. Luo, Y. Wu, Y. Peng, X. Liu, J. Bie, and S. Li, “Systematic analysis to identify a key role of CDK1 in mediating gene interaction networks in cervical cancer development,” Irish Journal of Medical Science, vol. 185, no. 1, pp. 231–239, 2016.
View at: Publisher Site | Google Scholar
J. Zhao, S.-X. Han, J.-L. Ma et al., “The role of CDK1 in apoptin-induced apoptosis in hepatocellular carcinoma cells,” Oncology Reports, vol. 30, no. 1, pp. 253–259, 2013.
View at: Publisher Site | Google Scholar
X. Qian, X. Song, Y. He et al., “CCNB2 overexpression is a poor prognostic biomarker in Chinese NSCLC patients,” Biomedicine & Pharmacotherapy, vol. 74, pp. 222–227, 2015.
View at: Publisher Site | Google Scholar
L. Zhou, J. Li, Y. P. Zhao et al., “The prognostic value of Cyclin B1 in pancreatic cancer,” Medical Oncology, vol. 31, no. 9, p. 107, 2014.
View at: Publisher Site | Google Scholar
B. Wang, Xunsun, J. Y. Liu et al., “The effect of cell cycle and expression of cyclin B1 and cyclin C protein in hepatocellular carcinoma cell line HepG2 and SMMC-7721 after of silencing β-catenin gene,” Hepatogastroenterology, vol. 59, no. 114, pp. 515–518, 2012.
View at: Publisher Site | Google Scholar
C. L. Gao, G. W. Wang, G. Q. Yang, H. Yang, and L. Zhuang, “Karyopherin subunit-α 2 expression accelerates cell cycle progression by upregulating CCNB2 and CDK1 in hepatocellular carcinoma,” Oncology Letters, vol. 15, no. 3, pp. 2815–2820, 2018.
View at: Publisher Site | Google Scholar
Y. Li and R. Benezra, “Identification of a human mitotic checkpoint gene: hsMAD2,” Science, vol. 274, no. 5285, pp. 246–248, 1996.
View at: Publisher Site | Google Scholar
Z. Wang, D. Katsaros, Y. Shen et al., “Biological and clinical significance of MAD2L1 and BUB1, genes frequently appearing in expression signatures for breast cancer prognosis,” PLoS One, vol. 10, no. 8, article e0136246, 2015.
View at: Publisher Site | Google Scholar
Y. X. Shi, T. Zhu, T. Zou et al., “Prognostic and predictive values of CDK1 and MAD2L1 in lung adenocarcinoma,” Oncotarget, vol. 7, no. 51, pp. 85235–85243, 2016.
View at: Publisher Site | Google Scholar
Y. Li, W. Bai, and J. Zhang, “MiR-200c-5p suppresses proliferation and metastasis of human hepatocellular carcinoma (HCC) via suppressing MAD2L1,” Biomedicine & Pharmacotherapy, vol. 92, pp. 1038–1044, 2017.
View at: Publisher Site | Google Scholar
E. L. Baldwin and N. Osheroff, “Etoposide, topoisomerase II and cancer,” Current Medicinal Chemistry Anti-Cancer Agents, vol. 5, no. 4, pp. 363–372, 2005.
View at: Publisher Site | Google Scholar
M. F. de Resende, S. Vieira, L. T. D. Chinen et al., “Prognostication of prostate cancer based on TOP2A protein and gene assessment: TOP2A in prostate cancer,” Journal of Translational Medicine, vol. 11, no. 1, p. 36, 2013.
View at: Publisher Site | Google Scholar
E. Moretti, C. Desmedt, C. Biagioni et al., “TOP2A protein by quantitative immunofluorescence as a predictor of response to epirubicin in the neoadjuvant treatment of breast cancer,” Future Oncology, vol. 9, no. 10, pp. 1477–1487, 2013.
View at: Publisher Site | Google Scholar
J. Lan, H. Y. Huang, S. W. Lee et al., “TOP2A overexpression as a poor prognostic factor in patients with nasopharyngeal carcinoma,” Tumour Biology, vol. 35, no. 1, pp. 179–187, 2014.
View at: Publisher Site | Google Scholar
N. Wong, W. Yeo, W. L. Wong et al., “TOP2A overexpression in hepatocellular carcinoma correlates with early age onset, shorter patients survival and chemoresistance,” International Journal of Cancer, vol. 124, no. 3, pp. 644–652, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Meng Wang 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

913

Downloads

1198

Citations