Abstract
Background:
Hepatocellular carcinoma (HCC) is a type of primary liver tumor with poor prognosis and high mortality, and its molecular mechanism remains incompletely understood. This study aimed to use bioinformatics technology to identify differentially expressed genes (DEGs) in HCC pathogenesis, hoping to identify novel biomarkers or potential therapeutic targets for HCC research.
Methods:
The bioinformatics analysis of our research mostly involved the following two datasets: Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). First, we screened DEGs based on the R packages (limma and edgeR). Using the DAVID database, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were carried out. Next, the protein-protein interaction (PPI) network of the DEGs was built in the STRING database. Then, hub genes were screened through the cytoHubba plug-in, followed by verification using the GEPIA and Oncomine databases. We demonstrated differences in levels of the protein in hub genes using the Human Protein Atlas (HPA) database. Finally, the hub genes prognostic values were analyzed by the GEPIA database. Additionally, using the Comparative Toxicogenomics Database (CTD), we constructed the drug-gene interaction network.
Results:
We ended up with 763 DEGs, including 247 upregulated and 516 downregulated DEGs, that were mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process, and metabolism-related pathways. Through the constructed PPI network, it can be concluded that the P53 signaling pathway and the cell cycle are the most obvious in module analysis. From the PPI, we filtered out eight hub genes, and these genes were significantly upregulated in HCC samples, findings consistent with the expression validation results. Additionally, survival analysis showed that high level gene expression of CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2 were connected with the poor overall survival of HCC patients. Toxicogenomics analysis showed that only topotecan, oxaliplatin, and azathioprine could reduce the gene expression levels of all seven hub genes.
Conclusion:
The present study screened out the key genes and pathways that were related to HCC pathogenesis, which could provide new insight for the future molecularly targeted therapy and prognosis evaluation of HCC.
Introduction
Accounting for 75-85% of all primary liver cancer, hepatocellular carcinoma (HCC) is the main histological classification of liver cancer, which is the fourth most frequent cause of cancer-related death globally (; ). The liver is the second most common cancer-prone organ, after the lungs, as was shown by the recent cancer study in China (). On the whole, the estimated morbidity of HCC per 100,000 world standard population is 40.0 in males and 15.3 in females (). Major risk factors for HCC include genetic predisposition, epigenetic variation, chronic hepatitis B infection, hepatitis C virus infection, smoking, obesity, aflatoxin exposure, and diabetes (; ). Transplantation is the most useful way to treat HCC; however, after the transplantation process, the tumor recurrence and metastasis rates are high (; ). More than 70% of patients at advanced stage are not suitable for transplantation, whether due to the tumor burden or liver dysfunction (). Therefore, it is urgent to recognize new biomarkers that can act as molecular targets for therapy, and predictors of the prognosis of HCC. With the development of times and technological progress, microarray and high-throughput sequencing technologies have matured and become more reliable, and public databases are improving, such as the Gene Expression Omnibus (GEO)1 and the Cancer Genome Atlas (TCGA)2. The advancement of microarray () and high throughput sequencing technologies () has provided a highly efficient tools to explore key genetic or epigenetic changes in disease to identify biological markers that can be applied to disease diagnosis, therapy, and prognosis (; ; ; ). Additionally, the application of integrated bioinformatics methods in cancer research can solve the problem of different results due to errors caused by different technical platforms or small sample size, thus finding much valuable biological information (; ; ; ).
In this research, by analyzing and distinguishing genes in human HCC samples and normal hepatocyte samples using TCGA and GEO datasets, we firstly screened out differentially expressed genes (DEGs). Then, GO and KEGG pathway enrichment analyses were applied in the further exploration of the main biological functions, which regulated by the DEGs. After that, the final step is to utilize a protein–protein interaction (PPI) network, survival analyses and drug-gene interaction network analyses to ascertain crucial genes and pathways which affecting the pathogenic mechanism and prognosis of HCC patients.
Materials and Methods
Gene Expression Datasets
The microarray gene expression dataset of GSE121248, which comprises 70 hepatocellular carcinoma samples and 37 normal liver samples, was obtained from the GEO website and exploited as discovery dataset to identify DEGs. The included dataset met the following criteria: (1) dataset included human HCC samples and normal liver samples. (2) they contained at least ten samples. (3) dataset was obtained from the Affymetrix Human Genome U133 Plus 2.0 Array [HG-U133_Plus_2] microarray platform. The raw RNA sequencing data, which comprises 374 HCC samples and 50 normal liver tissue samples, was selected from the TCGA liver hepatocellular carcinoma (TCGA-LIHC) dataset and used as a validation dataset.
Identification of DEGs
We used the R language to analyze the original CEL files of the GSE121248 dataset.
The preprocessing procedures: using the “affy” R package to RMA background correction, Log2conversion, Quantile normalization, and Median polish algorithm summarization (; ). Using the bioconductor annotation package to convert microarray data probes into gene symbol. If multiple probes were mapped to a gene symbol, take the average value as the final expression value of the gene (). Next, | log2fold change (FC)| > 1 and adjusted p value <0.05 were used to select the DEGs between tumor and normal tissues using the LIMMA package (; ).
DEGs Validation Using the TCGA Dataset
The DEGs from the GSE121248 dataset were validated using the TCGA-LIHC dataset.
The edgeR package of R software was applied to normalize and analyze the TCGA-LIHC dataset (). | log2fold change (FC)| > 1 and p-value <0.05 were considered significant differences. The overlapping DEGs between GSE121248 and TCGA-LIHC datasets were clustered using the pheatmap and were retained for further study. The overlapping DEGs were analyzed using VennDiagram and ggplot2 packages in R software to draw Venn diagrams and volcano plots, to visualize the identified DEGs ().
Functional Enrichment Analysis of Overlapping DEGs
We used the Database for Annotation, Visualization and Integrated Discovery (DAVID version 6.8)3 to elucidate potential GO function [including biological processes (BP), molecular functions (MF), cellular components (CC)] and signaling pathways (KEGG) related to the overlapping DEGs (; ). We used threshold p-value 0.05.
Protein–Protein Interaction Network Construction and Module Analysis
The Search Tool for the Retrieval of Interacting Genes (STRING version 11)4 database was one of the largest online databases of known protein-protein interactions covering the largest number of species (). The parameter of interactions was set with a confidence score >0.7. The confidence score refers to the strength of data support in terms of the thickness of the line. Confidence score >0.7 means high confidence. Overlapping DEGs were entered into Cytoscape software (version 3.7.2)5 to construct and analyze PPI network (). Moreover, the Cytoscape plug-in MCODE was used to screen crucial clustering modules in the entire network ().
Identification of Hub Genes
The Cytoscape plug-in CytoHubba was used to calculate the protein node degree (; ). The top three methods [(Maximal Clique Centrality (MCC), Maximum Neighborhood Component (MNC), and Density of Maximum Neighborhood Component (DMNC)] were selected to provide the analyzed results. Each method displayed their top ten genes. A Venn diagram was generated to visualize common hub genes based on these three methods.
Expression Analysis of Hub Genes in Multiple Databases
The hub genes mRNA expression levels were finally validated in two databases, Gene Expression Profiling Interactive Analysis (GEPIA)6 () and Oncomine. Oncomine (Version4.5)7 is an online database that has the comprehensive cancer mutation spectrum, gene expression data and related clinical information, which can be used to discover new biomarkers or new therapeutic targets (). In addition to detecting the mRNA expression levels of the hub genes, we also investigated the protein levels in HCC tissues and normal liver tissues using the human protein atlas database (HPA v19)8 ().
Survival Analysis
Gene Expression Profiling Interactive Analysis is a newly developed online database for cancer and normal gene expression profiling. In the current study, the overall survival of each hub gene was analyzed using LIHC dataset in the GEPIA database. The patients were divided into two groups (the high- and low-expression group) according to the median expression level of each hub gene. This division method could evaluate the difference in overall survival probability between these two groups. We were drawn the overall survival curves of each hub gene using the GEPIA database, with a p-value <0.05.
Drug-Gene Interaction Network Analysis
The Comparative Toxicogenomics Database (CTD)9, an online database providing information on the interactions between gene products and chemotherapeutic drugs, and their relationships to diseases) was used to construct the chemotherapeutic drug-gene interaction network (). The networks were visualized by Cytoscape software 3.7.210.
Results
Identification of DEGs
The gene expression dataset of GSE121248, which contains 70 LIHC samples and 37 normal liver samples, was analyzed in the limma package using | logFC| > 1 and corrected p-value <0.05 of R software. In total, 1,518 DEGs (557 high expression genes and 961 low expression genes) were identified between HCC tissue samples and normal liver tissue samples. The volcano map and heatmap of all DEGs are shown in Figures 1A,C. Additionally, compared with normal liver tissues in the TCGA-LIHC dataset, 2,898 DEGs were obtained in LIHC tissues, comprising 1,299 upregulated genes and 1,599 downregulated genes (Figure 1B). Furthermore, 763 overlapping DEGs (247 high expression genes and 516 low expression genes) were identified between the GSE121248 and TCGA-LIHC datasets using a Venn diagram (Figure 1D). Figure 1E shows clustering analysis results of the 763 overlapping DEGs based on the TCGA-LIHC dataset.
FIGURE 1
Enrichment Analysis of Overlapping DEGs
We conducted GO and KEGG pathway enrichment analysis to further elucidate potential biological functions associated with the 763 overlapping DEGs of HCC. The GO analysis results of the DEGs were classified into molecular functions, biological processes and cellular components. For molecular functions, the overlapping DEGs were mainly associated with oxidoreductase activity, monooxygenase activity, heme binding and oxygen binding (Figure 2A). In the BP category, the epoxygenase P450 pathway, oxidation-reduction process, response to drug and cell division were enriched (Figure 2B). In the CC category, they were enriched in extracellular regions, such as extracellular exosomes and the extracellular space (Figure 2C). The pathway enrichment analysis results showed that overlapping DEGs mainly participated in multiple metabolism pathways, such as fatty acid degradation, glycine, serine and threonine metabolism, and tryptophan metabolism (Figure 2D).
FIGURE 2
PPI Network Establishment and Module Analysis
To further reveal the potential relationships between proteins encoded by DEGs, a PPI network was constructed using the STRING database. Network analysis of overlapping DEGs revealed 526 nodes and 4,173 edges in the PPI network. Additionally, we conducted module analysis using the MCODE plug-in to detect crucial clustering modules. In total, 29 clusters were obtained in MCODE, and the top three modules with the highest scores were selected as hub modules. Module 1 contained 63 nodes and 1,752 edges with the highest score of 56.516 and was mainly enriched in cell cycle, oocyte meiosis, P53 signaling pathway and progesterone-mediated oocyte maturation (Figure 3A). Module 2 contained 17 nodes and 80 edges with a score of 10 and mainly participated in PPAR signaling pathway and glycerolipid metabolism (Figure 3B). Module 3 comprised 28 nodes and 100 edges with a score of 7.407 and was mainly implicated in Chemical carcinogenesis, Peroxisome, Metabolic pathways and Drug metabolism cytochrome P450 (Figure 3C).
FIGURE 3
Hub Genes Selection From the PPI Network
The Cytoscape plug-in cytoHubba including the top three algorithms (MCC, MNC, and DMNC) was applied to select hub genes, and the top 10 genes were selected by each of the three methods. The Venn diagram identified eight overlapping hub genes based on these three methods (Figure 3D): cell division cycle protein 20 homolog (CDC20), cyclin-dependent kinase1 (CDK1), mitotic spindle assembly checkpoint protein MAD2A (MAD2L1), threonine-protein kinase BUB1 (BUB1), threonine-protein kinase BUB1 beta (BUB1B), mitotic-specific cyclin-B1 (CCNB1), mitotic-specific cyclin-B2 (CCNB2) and cyclin-A2 (CCNA2). These eight hub genes were used for further analysis.
Validation of Hub Genes in Multiple Databases
Oncomine and GEPIA were applied to validate the differentially expression levels of 8 hub genes between HCC tissues and normal liver tissues in HCC. These eight hub genes were all remarkably overexpressed in HCC samples (Figure 4). Moreover, a summary of hub genes in multiple tumors indicated that hub genes were significantly overexpressed in HCC (Figure 5). Furthermore, we also investigated the protein expression levels in HCC tissue samples and normal liver tissue samples using the human protein atlas database. Because the HPA dataset could not provide immunohistochemical information on BUB1 and BUB1B, we showed the results of the remaining six staining pairs in Figure 6. The protein expression levels of hub genes were agreed with the mRNA expression results, and most genes were overexpressed in HCC tissue (Figure 7). These findings indicate that the overexpression of these hub genes may play a critical role in HCC mechanism.
FIGURE 4
FIGURE 5
FIGURE 6
FIGURE 7
Survival Analysis
We further used the GEPIA database to analyze the prognostic value of these 8 hub genes in HCC patients. The survival analysis of patients in the GEPIA database was based on the TCGA-LIHC data set. We used threshold p-value 0.05 and calculated the hazards ratio based on Cox PH Model (). The relatively higher expression of CDC20 (HR = 2.3; P = 3.8e-06), CDK1 (HR = 2; P = 0.00017), MAD2L1 (HR = 1.7; P = 0.0047), BUB1 (HR = 1.8; P = 0.001), BUBIB (HR = 1.7; P = 0.0028), CCNB1 (HR = 2; P = 0.00015), and CCNA2 (HR = 1.7; P = 0.0037) were associated with a poor prognosis in HCC patients, while only CCNB2 (HR = 1.4; P = 0.052) showed no statistical significance in the overall survival of patients (Figure 6).
Drug-Gene Interaction Network Analysis
To investigate the potential information on the interactions between hub genes and cancer chemotherapeutics drugs, we used the CTD database to construct chemotherapeutics drug-gene interaction network. Various drugs could influence the mRNA expression level of seven hub genes, namely, CDC20, CDK1, MAD2L1, CCNA2, CCNB1, BUB1, and BUB1B (Figure 8). However, only topotecan, oxaliplatin and azathioprine could reduce expression levels of all seven hub genes.
FIGURE 8
Discussion
Hepatocellular carcinoma is a type of primary liver tumor with poor prognosis and high mortality, and the progress in its diagnosis and treatment has always attracted widespread attention from researchers around the world. Because the high recurrence and metastasis rate of HCC remains a challenge, identifying new molecules as biological markers is urgently needed. Integrated bioinformatics analysis, which focuses on screening of DEGs, discovering hub node of network-based and doing survival analysis, which has been diffusely used to recognize latent biological markers related to cancer diagnosis, therapy, and prognosis estimation. In recent years, increasing researches have demonstrated that abnormal gene expression is a factor in the tumorigenesis and development, so it is feasible to screen differential genes as biomarkers to assist diagnosis and treatment. In 2017, by developing an integrated approach including GO and KEGG analysis, PPI network creation, hub gene identification, and overall survival analysis, picked out 16 hub genes for HCC from three GEO datasets, five of which may be playing a part in the occurrence, development, invasion, metastasis or recurrence of HCC. In 2018, used bioinformatics methods to select 10 genes from the GEO dataset GSE64041 for the identification of hub genes and pathways of HCC. recognized 13 crucial genes correlated with progression and prognosis of HCC from the TCGA-LIHC dataset by weighted gene coexpression network analysis. Compared with previous similar studies, our study not only integrated a large sample size of mRNA expression data from the GEO database but also analyzed RNA sequencing result and clinical data from the TCGA-LIHC database to screen out potential hub genes in HCC. And in the second place, this study validated the DEGs through multiple databases. Finally, we explored the relationship between seven hub genes and existing drugs for cancer therapy, which may provide some guidance for the molecular targeting therapy of HCC in the future.
In our research, DEGs in HCC based on the GEO expression profile of GSE121248 (70 HCC samples and 37 normal samples) and TCGA-LIHC RNA sequencing data (374 HCC samples and 50 normal samples) were identified by bioinformatics analysis. In total, 763 significantly robust DEGs, including 247 upregulated DEGs and 516 downregulated DEGs, were identified. The enrichment analysis results of GO indicated that the DEGs were mostly relevant to “oxidoreductase activity, acting on paired donors,” “monooxygenase activity,” “arachidonic acid epoxygenase activity,” “oxidation-reduction process,” “epoxygenase P450 pathway,” “response to drug,” “extracellular exosome,” “extracellular region,” and “cytosol.” The analysis of KEGG pathway showed that the DEGs were mainly concentrated in the following: “fatty acid degradation pathway,” “metabolic pathways,” “chemical carcinogenesis pathway,” “cell cycle pathway,” and “biosynthesis of antibiotics pathway.” Previous studies have reported that the arachidonic acid-derived metabolites and cytochrome P450 epoxygenase CYP2J2 possibly play vital roles in regulating malignant tumor, stimulating tumor cell growth, and inhibiting tumor cell apoptosis (; ; ). Additionally, metabolic pathways are important for cancer cell survival because the metabolic demands of cancer cells are often expressed as increased, and HCC shows a significant alteration in lipid metabolism (). Moreover, dysregulation of the cell cycle processes and mitotic cell cycle plays a vital role in the tumorigenesis and progression (; ). These theories are consistent with our results in GO and KEGG enrichment analysis.
Through building PPI network and analyzing it, we identified crucial hub genes in the PPI network, including CDC20, CDK1, MAD2L1, BUB1, BUB1B CCNB1, CCNB2, and CCNA2. Using Oncomine and GEPIA validation, the mRNA expression of these eight hub genes in HCC samples was higher than normal liver samples, the finding that was in accord with the microarray results. Subsequently, HPA database data displayed that the protein and mRNA expression of hub genes were consistent, and most genes were overexpressed in HCC tissue. To inquire prognostic biological markers of HCC, we applied the GEPIA to analyze the influence of hub genes expression level on survival of HCC patients and found that, except CCNB2, the high level gene expression of CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2 were related to HCC patients poor overall survival. Therefore, these seven genes may be functional in HCC occurrence and development.
It was reported that high expression of CDC20 (cell division cycle protein 20) is associated with poor survival in astrocytoma (), cutaneous squamous cell carcinoma () and pancreatic ductal adenocarcinoma (). CDC20 promotes the progression of prostate cancer by stabilizing hypo-catenin in tumor-like dry cells (). However, the expression of cell division cycle protein 20 in HCC still lacks accurate experimental data. As a part of the Ser/Thr protein kinase family, CDK1 (cyclin-dependent kinase 1) is a key molecule that controls the eukaryotic cell cycle. By phosphorylating Bora, Cyclin A/cdk1 could facilitate the phosphorylation, activation and mitotic entry of Aurora A-dependent Plk1 (). It is reported that CDK1 overexpression has been found in colorectal cancer, pancreatic ductal adenocarcinoma and thyroid cancer (; ; ). It was also reported that CDK1 amplification rate in HCC tissues was usually up to 46% (18/39), which was meaningfully related to poor overall survival (p = 0.008) (). These results were in accord with our study findings.
As a pro-oncogene upregulated in gastric cancer, MAD2L1 (mitotic arrest deficient 2-like protein 1) can be downregulated expression by miR-30a-3p, resulting in inhibition of the proliferation of gastric cancer cells (). Besides, by restraining MAD2L1, miR-200c-5p can inhibit HCC cells proliferation, migration and invasion (), suggesting that MAD2L1 can be used in HCC patients prognostic evaluation and targeted therapy. As a cyclin controlling the G1/S and G2/M phases in the cell cycle, CCNA2 (cyclin-A2) is more expressed in CRC samples than in normal samples. The reduction of CCNA2 gene expression would disrupt cell cycle progression and induce apoptosis, thus significantly inhibiting the growth of CRC cells (). By maintaining the expression of CCNA2 protein and the production of arginine, arginine metabolic enzyme argininosuccinate lyase (ASL) can promote the production of nitric oxide synthase, thus promoting the formation of HCC ().
As a mitotic checkpoint serine/threonine kinase, BUB1 is related to tumorigenesis in many cancers. shRNA silencing inhibits the expression of BUB1 gene in glioblastoma tumor cells, thereby reducing the proliferation and tumorigenicity of tumor cells in vivo and in vitro (). Increased BUB1 expression signally facilitates cell proliferation, while decreased BUB1 expression restrains liver cancer cells proliferation (). The proliferation, migration, and invasion of PCa cell lines can be enhanced via BUB1B overexpression (). Worse OS and DFS of HCC patients can be predicted by the high expression of BUB1B (). CCNB1, an important protein regulating the G2/M (mitotic) cell cycle, is activated by Chk1, exerting its oncogenic role in colorectal cancer cells growth in vivo and in vitro (). Abnormal FOXM1 expression can transcriptionally activate CCNB1 expression, thereby promoting the proliferation of HCC cells ().
After exploring the potential information about the interactions between the seven hub genes and existing chemotherapeutic drugs, we found that various drugs could influence the expression levels of these hub genes. However, only topotecan, oxaliplatin and azathioprine could simultaneously reduce the expression level all seven hub genes. And it should be noted that further experiments are needed to support whether HCC patients with hub gene overexpression can benefit from hub gene inhibition or whether these key genes may be targets of drug treatment of tumor need ulteriorly biological experiments support.
In the current study, we have discussed that the development of HCC is associated with the overexpression of seven hub genes, which lead to poor overall survival, indicating that they may be considered as potential prognostic biomarkers for HCC. However, our study has several limitations: (1) some important clinical information (for example, different age, tumor size, TNM stage and grade) were not considered; (2) biological experiments must be carried out in the future to verify the results of our research; (3) the molecular mechanism of hub gene upregulation remains unclear. Therefore, the verification of hub genes will be the focus of our next work.
Conclusion
Adopting a series of bioinformatics analysis methods, the current study identified 763 DEGs and seven hub genes (CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2) that may be involved in hepatocellular carcinoma tumorigenesis and progression. Additionally, multiple database analysis and survival analysis demonstrated that these seven hub genes may regard as a latent prognostic biomarker and the overexpression of these seven hub genes might lead to reduced overall survival in HCC patients. These results provide a theoretical basis for the molecularly targeted therapy and prognosis evaluation of HCC.
Statements
Data availability statement
The data in this study were obtained from the GEO (GSE121248) and TCGA public databases, and the acquisition and application methods complied with the corresponding database guidelines and policies.
Author contributions
ZL conceived the study, analyzed the raw data, and wrote the manuscript. YL prepared the diagrams in the manuscript. BC and QZ provided the useful suggestions in methodology. YC designed and instructed the research. All the authors have read and approved the final manuscript.
Funding
This study was funded by the Shantou Medical and health plan (Grant No. 180404094011013).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.571231/full#supplementary-material
Supplementary Table 1List of differentially expressed genes in GSE121248.
Supplementary Table 2List of differentially expressed genes between HCC tissue and normal liver tissue in the TCGA dataset.
Supplementary Table 3List of common differentially expressed genes both in GSE121248 and the TCGA-LIHC dataset.
Supplementary Table 4Gene ontology enrichment analysis of the common differentially expressed genes both in GSE121248 and the TCGA-LIHC dataset.
Supplementary Table 5Pathway enrichment analysis of the common DEGs both in GSE121248 and the TCGA-LIHC dataset.
Supplementary Table 6STRING significance scores of common differentially expressed genes.
Footnotes
1.^https://www.ncbi.nlm.nih.gov/geo/
2.^http://tcga-data.nci.nih.gov
6.^http://gepia.cancer-pku.cn/
References
1
AuK. P.ChokK. S. H. (2018). Multidisciplinary approach for post-liver transplant recurrence of hepatocellular carcinoma: a proposed management algorithm.World J. Gastroenterol.245081–5094. 10.3748/wjg.v24.i45.5081
2
AufhauserD. D.Jr.SadotE.MurkenD. R.EddingerK.HoteitM.AbtP. L.et al (2018). Incidence of occult intrahepatic metastasis in hepatocellular carcinoma treated with transplantation corresponds to early recurrence rates after partial hepatectomy.Ann. Surg.267922–928. 10.1097/sla.0000000000002135
3
BaderG. D.HogueC. W. (2003). An automated method for finding molecular complexes in large protein interaction networks.BMC Bioinformatics4:2. 10.1186/1471-2105-4-2
4
BaeckerA.LiuX.La VecchiaC.ZhangZ. F. (2018). Worldwide incidence of hepatocellular carcinoma cases attributable to major risk factors.Eur. J. Cancer Prev.27205–212. 10.1097/cej.0000000000000428
5
BolstadB. M.IrizarryR. A.AstrandM.SpeedT. P. (2003). A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.Bioinformatics19185–193. 10.1093/bioinformatics/19.2.185
6
CaoL.ChenY.ZhangM.XuD. Q.LiuY.LiuT.et al (2018). Identification of hub genes and potential molecular mechanisms in gastric cancer by integrated bioinformatics analysis.PeerJ6:e5180. 10.7717/peerj.5180
7
ChaiN.XieH. H.YinJ. P.SaK. D.GuoY.WangM.et al (2018). FOXM1 promotes proliferation in human hepatocellular carcinoma cells by transcriptional activation of CCNB1.Biochem. Biophys. Res. Commun.500924–929. 10.1016/j.bbrc.2018.04.201
8
ChenH.BoutrosP. C. (2011). VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R.BMC Bioinformatics12:35. 10.1186/1471-2105-12-35
9
ChinC. H.ChenS. H.WuH. H.HoC. W.KoM. T.LinC. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome.BMC Syst. Biol.8 Suppl. 4(Suppl. 4):S11. 10.1186/1752-0509-8-s4-s11
10
DavisA. P.GrondinC. J.JohnsonR. J.SciakyD.McMorranR.WiegersJ.et al (2019). The comparative toxicogenomics database: update 2019.Nucleic Acids Res.47D948–D954. 10.1093/nar/gky868
11
DengJ. L.XuY. H.WangG. (2019). Identification of potential crucial genes and key pathways in breast cancer using bioinformatic analysis.Front. Genet.10:695. 10.3389/fgene.2019.00695
12
DennisG.Jr.ShermanB. T.HosackD. A.YangJ.GaoW.LaneH. C.et al (2003). DAVID: database for annotation, visualization, and integrated discovery.Genome Biol.4:3.
13
DingY.YuS.BaoZ.LiuY.LiangT. (2017). CDC20 with malignant progression and poor prognosis of astrocytoma revealed by analysis on gene expression.J. Neurooncol.13387–95. 10.1007/s11060-017-2434-8
14
DongS.HuangF.ZhangH.ChenQ. (2019). Overexpression of BUB1B, CCNA2, CDC20, and CDK1 in tumor tissues predicts poor survival in pancreatic ductal adenocarcinoma.Biosci. Rep.391–10. 10.1042/bsr20182306
15
FangY.YuH.LiangX.XuJ.CaiX. (2014). Chk1-induced CCNB1 overexpression promotes cell proliferation and tumor growth in human colorectal cancer.Cancer Biol. Ther.151268–1279. 10.4161/cbt.29691
16
FuJ.WangH. (2018). Precision diagnosis and treatment of liver cancer in China.Cancer Lett.412283–288. 10.1016/j.canlet.2017.10.008
17
FuX.ChenG.CaiZ. D.WangC.LiuZ. Z.LinZ. Y.et al (2016). Overexpression of BUB1B contributes to progression of prostate cancer and predicts poor outcome in patients with prostate cancer.Onco Targets Ther.92211–2220. 10.2147/ott.S101994
18
GanY.LiY.LiT.ShuG.YinG. (2018). CCNA2 acts as a novel biomarker in regulating the growth and apoptosis of colorectal cancer.Cancer Manag. Res.105113–5124. 10.2147/cmar.S176833
19
GautierL.CopeL.BolstadB. M.IrizarryR. A. (2004). affy–analysis of Affymetrix GeneChip data at the probe level.Bioinformatics20307–315. 10.1093/bioinformatics/btg405
20
GuY.LiJ.GuoD.ChenB.LiuP.XiaoY.et al (2020). Identification of 13 key genes correlated with progression and prognosis in hepatocellular carcinoma by weighted gene co-expression network analysis.Front. Genet.11:153. 10.3389/fgene.2020.00153
21
HarrisP. S.HansenR. M.GrayM. E.MassoudO. I.McGuireB. M.ShoreibahM. G. (2019). Hepatocellular carcinoma surveillance: an evidence-based approach.World J. Gastroenterol.251550–1559. 10.3748/wjg.v25.i13.1550
22
HungY. H.HuangH. L.ChenW. C.YenM. C.ChoC. Y.WengT. Y.et al (2017). Argininosuccinate lyase interacts with cyclin A2 in cytoplasm and modulates growth of liver tumor cells.Oncol. Rep.37969–978. 10.3892/or.2016.5334
23
KanehisaM.FurumichiM.TanabeM.SatoY.MorishimaK. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs.Nucleic Acids Res.45D353–D361. 10.1093/nar/gkw1092
24
LiC.QinF.HongH.TangH.JiangX.YangS.et al (2019). Identification of Flap endonuclease 1 as a potential core gene in hepatocellular carcinoma by integrated bioinformatics analysis.PeerJ7:e7619. 10.7717/peerj.7619
25
LiL.LeiQ.ZhangS.KongL.QinB. (2017). Screening and identification of key biomarkers in hepatocellular carcinoma: evidence from bioinformatic analysis.Oncol. Rep.382607–2618. 10.3892/or.2017.5946
26
LiY.BaiW.ZhangJ. (2017). MiR-200c-5p suppresses proliferation and metastasis of human hepatocellular carcinoma (HCC) via suppressing MAD2L1.Biomed. Pharmacother.921038–1044. 10.1016/j.biopha.2017.05.092
27
LiuL.ChenC.GongW.LiY.EdinM. L.ZeldinD. C.et al (2011). Epoxyeicosatrienoic acids attenuate reactive oxygen species level, mitochondrial dysfunction, caspase activation, and apoptosis in carcinoma cells treated with arsenic trioxide.J. Pharmacol. Exp. Ther.339451–463. 10.1124/jpet.111.180505
28
LiuX.WuJ.ZhangD.BingZ.TianJ.NiM.et al (2018). Identification of potential key genes associated with the pathogenesis and prognosis of gastric cancer based on integrated bioinformatics analysis.Front. Genet.9:265. 10.3389/fgene.2018.00265
29
NagyÁLánczkyA.MenyhártO.GyőrffyB. (2018). Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets.Sci. Rep.8:9227. 10.1038/s41598-018-27521-y
30
PiaoJ.ZhuL.SunJ.LiN.DongB.YangY.et al (2019). High expression of CDK1 and BUB1 predicts poor prognosis of pancreatic ductal adenocarcinoma.Gene70115–22. 10.1016/j.gene.2019.02.081
31
PopeE. D.IIIKimbroughE. O.VemireddyL. P.SurapaneniP. K.CoplandJ. A.IIIModyK. (2019). Aberrant lipid metabolism as a therapeutic target in liver cancer.Expert Opin. Ther. Targets23473–483. 10.1080/14728222.2019.1615883
32
PuszykW. M.TrinhT. L.ChappleS. J.LiuC. (2013). Linking metabolism and epigenetic regulation in development of hepatocellular carcinoma.Lab. Invest.93983–990. 10.1038/labinvest.2013.94
33
RhodesD. R.YuJ.ShankerK.DeshpandeN.VaramballyR.GhoshD.et al (2004). ONCOMINE: a cancer microarray database and integrated data-mining platform.Neoplasia61–6. 10.1016/s1476-5586(04)80047-2
34
RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies.Nucleic Acids Res.43:e47. 10.1093/nar/gkv007
35
RobinsonM. D.McCarthyD. J.SmythG. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.Bioinformatics26139–140. 10.1093/bioinformatics/btp616
36
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504. 10.1101/gr.1239303
37
SzklarczykD.MorrisJ. H.CookH.KuhnM.WyderS.SimonovicM.et al (2017). The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible.Nucleic Acids Res.45D362–D368. 10.1093/nar/gkw937
38
TangZ.LiC.KangB.GaoG.LiC.ZhangZ. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses.Nucleic Acids Res.45W98–W102. 10.1093/nar/gkx247
39
ThulP. J.LindskogC. (2018). The human protein atlas: a spatial map of the human proteome.Protein Sci.27233–244. 10.1002/pro.3307
40
VigneronS.SundermannL.LabbéJ.-C.PintardL.RadulescuO.CastroA.et al (2018). Cyclin A-cdk1-dependent phosphorylation of bora is the triggering factor promoting mitotic entry.Dev. Cell45637–650.e7. 10.1016/j.devcel.2018.05.005
41
WangJ.LiH.XiaC.YangX.DaiB.TaoK.et al (2019). Downregulation of CENPK suppresses hepatocellular carcinoma malignant progression through regulating YAP1.Onco Targets Ther.12869–882. 10.2147/ott.S190061
42
WangX. B.CuiN. H.LiuX. N.MaJ. F.ZhuQ. H.GuoS. R.et al (2018). Identification of promoter hypermethylation as a biomarker for intra-epithelial lesion and cervical cancer: a meta-analysis of published studies, TCGA, and GEO datasets.Front. Genet.9:258. 10.3389/fgene.2018.00258
43
WangY.WangF.HeJ.DuJ.ZhangH.ShiH.et al (2019). miR-30a-3p targets MAD2L1 and regulates proliferation of gastric cancer cells.Onco Targets Ther.1211313–11324. 10.2147/ott.S222854
44
WeinsteinJ. N.CollissonE. A.MillsG. B.ShawK. R.OzenbergerB. A.EllrottK.et al (2013). The cancer genome Atlas pan-cancer analysis project.Nat. Genet.451113–1120. 10.1038/ng.2764
45
WilliamsG. H.StoeberK. (2012). The cell cycle and cancer.J. Pathol.226352–364. 10.1002/path.3022
46
WlodarchakN.XingY. (2016). PP2A as a master regulator of the cell cycle.Crit. Rev. Biochem. Mol. Biol.51162–184. 10.3109/10409238.2016.1143913
47
WuC. X.WangX. Q.ChokS. H.ManK.TsangS. H. Y.ChanA. C. Y.et al (2018). Blocking CDK1/PDK1/β-Catenin signaling by CDK1 inhibitor RO3306 increased the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma.Theranostics83737–3750. 10.7150/thno.25487
48
XuJ. S.WangX.HeJ.LiaoK. L.WangX. Z. (2020). Combining bioinformatics techniques to explore the molecular mechanisms involved in pancreatic cancer metastasis and prognosis.J. Cell. Mol. Med.2414128–14138. 10.1111/jcmm.16023
49
XuX.ZhangX. A.WangD. W. (2011). The roles of CYP450 epoxygenases and metabolites, epoxyeicosatrienoic acids, in cardiovascular and malignant diseases.Adv. Drug Deliv. Rev.63597–609. 10.1016/j.addr.2011.03.006
50
YanH.ZhengG.QuJ.LiuY.HuangX.ZhangE.et al (2019). Identification of key candidate genes and pathways in multiple myeloma by integrated bioinformatics analysis.J. Cell. Physiol.23423785–23797. 10.1002/jcp.28947
51
YangJ. D.HainautP.GoresG. J.AmadouA.PlymothA.RobertsL. R. (2019). A global view of hepatocellular carcinoma: trends, risk, prevention and management.Nat. Rev. Gastroenterol. Hepatol.16589–604. 10.1038/s41575-019-0186-y
52
YangK.ZhangS.ZhangD.TaoQ.ZhangT.LiuG.et al (2019). Identification of SERPINE1, PLAU and ACTA1 as biomarkers of head and neck squamous cell carcinoma based on integrated bioinformatics analysis.Int. J. Clin. Oncol.241030–1041. 10.1007/s10147-019-01435-9
53
YangX.ZhuS.LiL.ZhangL.XianS.WangY.et al (2018). Identification of differentially expressed genes and signaling pathways in ovarian cancer by integrated bioinformatics analysis.Onco Targets Ther.111457–1474. 10.2147/ott.S152238
54
YarlaN. S.BishayeeA.SethiG.ReddannaP.KalleA. M.DhananjayaB. L.et al (2016). Targeting arachidonic acid pathway by natural products for cancer prevention and therapy.Semin. Cancer Biol.448–81. 10.1016/j.semcancer.2016.02.001
55
YuH.ZhangS.IbrahimA. N.DengZ.WangM. (2019). Serine/threonine kinase BUB1 promotes proliferation and radio-resistance in glioblastoma.Pathol. Res. Pract.215:152508. 10.1016/j.prp.2019.152508
56
ZhangL.HuangY.LingJ.ZhuoW.YuZ.ShaoM.et al (2018). Screening and function analysis of hub genes and pathways in hepatocellular carcinoma via bioinformatics approaches.Cancer Biomark.22511–521. 10.3233/cbm-171160
57
ZhangP.KawakamiH.LiuW.ZengX.StrebhardtK.TaoK.et al (2018). Targeting CDK1 and MEK/ERK overcomes apoptotic resistance in BRAF-mutant human colorectal cancer.Mol. Cancer Res.16378–389. 10.1158/1541-7786.Mcr-17-0404
58
ZhangQ.HuangH.LiuA.LiJ.LiuC.SunB.et al (2019). Cell division cycle 20 (CDC20) drives prostate cancer progression via stabilization of beta-catenin in cancer stem-like cells.EBioMedicine42397–407. 10.1016/j.ebiom.2019.03.032
59
ZhangX.FengH.LiZ.LiD.LiuS.HuangH.et al (2018). Application of weighted gene co-expression network analysis to identify key modules and hub genes in oral squamous cell carcinoma tumorigenesis.Onco Targets Ther.116001–6021. 10.2147/ott.S171791
60
ZhengH. P.HuangZ. G.HeR. Q.LuH. P.DangY. W.LinP.et al (2019). Integrated assessment of CDK1 upregulation in thyroid cancer.Am. J. Transl. Res.117233–7254.
61
ZhuL. J.PanY.ChenX. Y.HouP. F. (2020). BUB1 promotes proliferation of liver cancer cells by activating SMAD2 phosphorylation.Oncol. Lett.193506–3512. 10.3892/ol.2020.11445
62
ZhuR. X.SetoW. K.LaiC. L.YuenM. F. (2016). Epidemiology of hepatocellular carcinoma in the Asia-pacific region.Gut Liver10332–339. 10.5009/gnl15257
63
ZhuangL.YangZ.MengZ. (2018). Upregulation of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in tumor tissues predicted worse overall survival and disease-free survival in hepatocellular carcinoma patients.Biomed. Res. Int.2018:7897346. 10.1155/2018/7897346
Summary
Keywords
hepatocellular carcinoma, bioinformatics, differentially expressed genes, survival, biomarker, GEO, TCGA
Citation
Li Z, Lin Y, Cheng B, Zhang Q and Cai Y (2021) Identification and Analysis of Potential Key Genes Associated With Hepatocellular Carcinoma Based on Integrated Bioinformatics Methods. Front. Genet. 12:571231. doi: 10.3389/fgene.2021.571231
Received
10 June 2020
Accepted
18 February 2021
Published
09 March 2021
Volume
12 - 2021
Edited by
Ancha Baranova, George Mason University, United States
Reviewed by
Yuriy L. Orlov, I.M. Sechenov First Moscow State Medical University, Russia; Evgenii Chekalin, Michigan State University, United States
Updates
Copyright
© 2021 Li, Lin, Cheng, Zhang and Cai.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yingmu Cai, st_ymcai@163.com
This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.