CCNB2, TOP2A, and ASPM Reflect the Prognosis of Hepatocellular Carcinoma, as Determined by Weighted Gene Coexpression Network Analysis

Zeng, Yuping; He, He; Zhang, Yu; Wang, Xia; Yang, Lidan; An, Zhenmei

doi:https://doi.org/10.1155/2020/4612158

BioMed Research International

On this page

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

Research Article | Open Access

Volume 2020 | Article ID 4612158 | https://doi.org/10.1155/2020/4612158

CCNB2, TOP2A, and ASPM Reflect the Prognosis of Hepatocellular Carcinoma, as Determined by Weighted Gene Coexpression Network Analysis

Yuping Zeng,¹He He,¹Yu Zhang,²Xia Wang,¹Lidan Yang,¹and Zhenmei An²

Academic Editor: Udayan Apte

Received24 Dec 2019

Revised15 May 2020

Accepted02 Jun 2020

Published24 Jun 2020

Abstract

Background. Hepatocellular carcinoma (HCC) is characterized by increased mortality and poor prognosis. We aimed to identify potential prognostic markers by weighted gene coexpression network analysis (WGCNA), to assist clinical outcome prediction and improve treatment decisions for HCC patients. Methods. Prognosis-related gene modules were first established by WGCNA. Venn diagrams obtained intersection genes of module genes and differentially expressed genes. The Kaplan-Meier overall survival curves and disease-free survival curves of intersection genes were further analyzed on the Gene Expression Profiling Interactive Analysis website. Chi-square tests were performed to explore the associations between prognostic gene expressions and clinicopathological features. Results. CCNB2, TOP2A, and ASPM were identified as both prognosis-related genes and differentially expressed genes. TOP2A (HR: 1.7, ) and ASPM (HR: 1.8, ) exhibited a significant difference between the high- and low-expression groups in the overall survival analysis, while CCNB2 (HR: 1.4, ) was not statistically significant. CCNB2 (HR: 1.5, ), TOP2A (HR: 1.7, ), and ASPM (HR: 1.6, ) were all statistically significant in the disease-free survival analysis. All three genes were significantly associated with race and fetoprotein values (). CCNB2 expression was associated with tumor stage (), and ASPM expression was associated with new tumor events (). Conclusion. Overexpression of CCNB2, TOP2A, and ASPM are associated with poor prognosis, and these genes could serve as potential prognostic markers and therapeutic targets for HCC.

1. Introduction

Hepatocellular carcinoma (HCC) is the third most common cause of cancer mortality worldwide [1]. The mechanisms by which hepatitis B virus, hepatitis C virus, alcohol, fatty liver disease, and other environmental factors, such as aflatoxin, cause liver cancer remain unclear. However, advances in genomics provide essential information about tumor initiation and progression [2]. The complex genetic background of HCC makes current clinical staging methods, including the Barcelona Clinic Liver Cancer algorithm or the tumor node metastasis (TNM) staging, insufficient to predict patient prognosis. Indeed, patients with the same HCC stage may have significantly different outcomes. Besides, the high recurrence and metastasis rates after ablation or surgical resection lead to a low survival rate in patients with HCC. These outcomes present an urgent requirement for improved prognostic estimates other than staging [3]. To date, thousands of HCC genomes have been sequenced globally, and most driver gene mutations, structural variants, fusion genes, copy number alterations, and viral integration events have been established [4]. This vast amount of new data provides an opportunity to understand the molecular basis of HCC better [5]. However, to utilize this scientific information, knowledge of the available database resources and bioinformatics tools are indispensable. To improve the treatment decisions and the overall survival of patients with HCC, we summarized the currently available databases supporting HCC research to aid in the identification of potential prognostic markers of clinical outcome prediction.

A coexpression analysis is an efficient method to describe free-scale gene coexpression networks. The weighted gene coexpression network analysis (WGCNA), an algorithm based on large-scale datasets and modules of highly correlated genes, was used to explore associations between gene sets and clinical features and to identify candidate biomarkers [6, 7]. This approach has been successfully applied in multiple tumors, such as clear cell renal cell carcinoma [8], glioblastoma [6], pancreatic carcinoma [9], adrenocortical carcinoma [10], breast cancer [11], and HCC [12, 13]. We mined prognostic markers by constructing a coexpression network and performing differential gene analysis and survival analysis to verify their prognostic values.

2. Materials and Methods

2.1. Data Sources

HCC-related gene expression profiles were downloaded from Gene Expression Omnibus (GEO) datasets (https://www.ncbi.nlm.nih.gov/geo/).The search details used were as follows: (“2010”(UDAT): “3000”(UDAT)) and (((“carcinoma, hepatocellular”(MeSH terms) or hepatocellular carcinoma (all fields)) and (“liver neoplasms”(MeSH terms) or liver cancer (all fields)) and (“mortality”(subheading) or “survival”(MeSH terms) or survival (all fields))) and “Homo sapiens”(porgn) and “Expression profiling by array”(filter)) and “Expression profiling by array”(filter). Four datasets contained HCC and paracancerous nontumor tissues, and samples over 30 were included for our analysis. GSE54236 [14] was utilized for WGCNA analysis because there were clinical features such as survival time in the dataset. GSE60502 [15], GSE64041 [16], and GSE45114 [17] were used for differential expression gene (DEG) analysis to obtain reliable DEGs. Standardized series matrix file and annotation platform were extracted through the R package “GEOquery” [18]. The data uniformity between samples was judged by boxplot and log2 conversion was performed to standardize datasets if necessary.

2.2. WGCNA Analysis

To identify hub genes associated with phenotypes, we analyzed GSE54236 by the R package “WGCNA” as follows [19, 20]. (1) Data procession: the median absolute deviation (MAD) of each gene expression between samples was calculated. We selected genes with MAD greater than the 70% quantile interval of the MAD of all genes and greater than 0.01 for the subsequent analysis. Missing values were checked by the “goodSamplesGenes” function, and outlier samples were identified by the “hclust” function to eliminate possible sample interference. (2) Network construction: soft power was determined with a threshold of 0.85 by “pickSoftThreshold” function to make the constructed network in line with a scale-free network. Then, we constructed the network by “blockwiseModules” function and showed each module by a hierarchical clustering tree. (3) Modules related to phenotypes: To explore phenotype-related modules, we analyzed the associations between modules and phenotypes. We calculated the correlations of modules and genes (module membership), and phenotypes and genes (gene significance) by the Pearson correlation analysis. Genes highly associated with phenotypes in modules, with and were selected for further analysis. (4) Hub genes identified: the topological overlap matrix was calculated from module genes of interest, and hub nodes with an edge-adjacency threshold of 0.2 were exported to visualize in Cytoscape (version 3.7.2) by “exportNetworkToCytoscape” function.

2.3. Pathway Enrichment Analysis

We performed Gene Ontology (GO), including biological processes, cellular components, and molecular functions, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of module genes of interest by the “clusterProfiler” package [21]. Next, the top 10 most significant enrichment pathways were graphically displayed. Moreover, protein-protein interaction (PPI) analysis was performed on the Metascape website (http://metascape.org/gp/index.html#/main/step1), which determined the PPI network, Molecular Complex Detection (MCODE) components, and the top three enrichment pathways of each component.

2.4. DEG Analysis

We analyzed DEGs of GSE60502, GSE64041, and GSE45114 between HCC and paracancerous nontumor tissues by the “limma” [22] package. and adjusted values (Benjaminiand-Hochberg adjustment) < 0.05 were set as the screening cutoff for DEGs. The intersection of WGCNA phenotype-related module genes and DEGs was obtained and visualized in Venn diagrams (http://bioinfogp.cnb.csic.es/tools/venny/index.html), representing both prognosis-related and differentially expressed genes.

2.5. Survival Analysis

We performed the Kaplan-Meier overall survival curves and disease-free survival curves of intersection genes on the Gene Expression Profiling Interactive Analysis (GEPIA) website [23], a web-based tool for analyzing Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) project data from tumors and normal samples.

2.6. Correlation Analysis of Clinicopathological Features

To further explore the associations between prognosis gene expressions and clinicopathological features, RNA-seq-counts of the TCGA liver cancer () and phenotypes () were downloaded from the UCSC Xena (https://xenabrowser.net/datapages/). Genes were divided into high- and low-expression groups, according to the median of gene expressions. We compared the differences between the two groups in clinicopathological features by Chi-square tests. All analyses were performed by R (version 3.6.0). Bilateral values less than 0.05 were considered statistically significant.

3. Results

3.1. Weighted Coexpression Network Construction and Key Module Identification

To mine coexpression module genes for HCC, the top 30% MAD of all genes in GSE54236 were extracted, which provided 5,612 genes for further analysis. No outlier samples were identified in the sample clustering tree, so no samples were removed from the analysis (Figure 1(a)). Next, we constructed the gene-gene similarity network. We set the soft power as 7 to ensure a scale-free network (Figure 1(b)). By the “blockwiseModules” function, we divided the network into 11 modules with similar gene expressions (Figure 1(c)). Then, we analyzed the relationships between modules and phenotypes. Results showed that the red module had a significant negative correlation with doubling time and survival time (Figure 1(d)). Figures 1(e) and 1(f) further proved the associations between the red module genes and doubling time and survival time with and . Finally, we identified 153 hub genes in the red module as key genes by Cytoscape.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1

Weighted gene coexpression network analysis of GSE54236. (a) Sample clustering tree to detect outliers. (b) Identification of soft power with a threshold of 0.85. (c) Hierarchical clustering tree of each module. Different colors represent different modules. (d) Correlations of modules and doubling time and survival time. The correlation coefficient varied from -1 (blue) to 1 (red) and was annotated. (e) Correlations between gene-module and gene-survival time. (f) Correlations between gene-module and gene-doubling time.

3.2. Gene Function and Pathway Enrichment Analyses

To further understand the biological functions and signaling pathways in which red module genes participate, we performed enrichment analysis for 153 hub genes by GO and KEGG. Biological processes were mainly enriched in chromosome segregation, nuclear division, and organelle fission (Figure 2(a)). Associated cell components were mainly present in the chromosomal region, chromosome, centromeric region, and condensed chromosome (Figure 2(b)). Molecular functions occurred in microtubule motor activity, tubulin binding, and microtubule activity (Figure 2(c)). KEGG pathways were considerably enriched in the cell cycle and oocyte meiosis (Figure 2(d)). PPI and MCODE analyses revealed that red module genes were divided into five components (Figures 2(e) and 2(f)). The top three enrichment pathways for each component were shown in Supplementary Table 1. These genes mainly participated in the cell cycle and divisions.

(a)

(b)

(c)

(d)

(e)

(f)

3.3. Identification of Intersection Genes

We explored the intersection genes both related to prognosis and differentially expressed in tumors and paracancerous nontumor tissues by the Venn diagram. There were 510, 75, and 415 DEGs for GSE60502, GSE64041, and GSE45114, respectively. We found three intersection genes (CCNB2, TOP2A, and ASPM) which were both red module genes and DEGs in the Venn diagram (Figure 3(a)).

(a)

(b)

(c)

(d)

3.4. Identification and Validation of Hub Genes

To further validate hub genes, we analyzed expression levels (Figure 3(b)), overall survival curves (Figure 3(c)), and disease-free survival curves (Figure 3(d)) of CCNB2, TOP2A, and ASPM on the GEPIA website. All three genes were differentially expressed in HCC and normal tissues with . Genes were divided into high- and low-expression groups based on the median gene expression levels. The differences between the two groups were compared by the logrank test. TOP2A (HR: 1.7, ) and ASPM (HR: 1.8, ) had a statistical significance between the two groups in the overall survival analysis, while CCNB2 (HR: 1.4, ) was not statistically significant. CCNB2 (HR: 1.5, ), TOP2A (HR: 1.7, ), and ASPM (HR: 1.6, ) were statistically significant in the disease-free survival analysis.

3.5. Associations between Gene Expressions and Clinicopathological Features

To explore the associations between CCNB2, TOP2A, ASPM, and clinicopathological features, we divided three genes into the low- and high-expression groups according to the median of gene expressions in the TCGA. Through Chi-square tests, we found that all three genes were significantly associated with race and fetoprotein (AFP) values (). CCNB2 was significantly associated with tumor stage (), while there was no significant association between TOP2A or ASPM and tumor stage. Moreover, a significant association between ASPM and new tumor events was observed (), as shown in Table 1. We did not found significant associations between CCNB2, TOP2A, ASPM, and BMI, adjacent hepatic tissue inflammation, fibrosis Ishak score, residual tumor, vascular tumor, and Child-Pugh classification, as shown in Supplementary Table 2.

4. Discussion

Clinical prognostic information is limited for patients of HCC. Several staging systems that attempt to predict patient prognosis and guide treatment modality have been proposed for HCC. However, the clinical features of these systems have drawbacks concerning the efficacy of treatment guidance and prognostic accuracy. Staging systems based on clinical features lack universal applicability, having been designed using data mostly acquired from a Western population or by being too rigid in suggesting therapy for some problematic prognostic factors [24]. Therefore, it is necessary to identify prognostic markers for patients with HCC. Patients with HCC present with features of high heterogeneity, metastasis, recurrence, and mortality and provide potential targets for prognosis prediction and treatment monitoring [25]. On the other hand, there have been very few reports detailing data mining by combining coexpressed and differentially expressed genes to explore HCC prognostic markers.

In this study, we explored genes associated with HCC prognosis using WGCNA, a method more frequently used to investigate phenotype-related coexpression module genes than to analyze DEGs. We identified 153 hub genes in the red module that were significantly negatively associated with doubling time and survival time, symbolizing a poor prognosis. Besides, the strong correlations of modules and genes and phenotypes and genes further proved the robust prognostic efficacy of red module genes. Second, gene function and pathway enrichment analyses indicated that red module genes were mainly engaged in the cell division and proliferation of HCC. DEGs were mined by combining GSE60502, GSE64041, and GSE45114, making results more credible. Further combining overall and disease-free survival analysis using the GEPIA website, we identified that TOP2A, ASPM, and CCNB2 were promising prognostic markers and potential therapeutic targets for HCC. It is worth noting that the overall survival analysis was based on patient all-cause mortality, and disease-free survival analysis was based on disease progression or death. Therefore, overall survival analysis endpoints were optimal observation endpoints without subjective judgment bias. In our study, both overall and disease-free survival analyses were evaluated by public databases.

We further investigated associations between TOP2A, ASPM, CCNB2, and clinicopathological features by Chi-square tests. We found significant associations between three prognostic genes and AFP levels; therefore, we speculated that these three genes were potential markers of the development and progression of HCC. However, whether these genes might become favorable diagnostic markers for HCC and specific mechanisms required further experimental verification. Importantly, high CCNB2 expression was strongly associated with tumor stage (), and high ASPM expression was associated with new tumor events (). These findings indicated that the overexpression of CCNB2 and ASPM accelerated the possibility of HCC deterioration and metastasis, suggesting a worse prognosis for HCC patients; hence, more aggressively integrated treatment and assessment were needed when these genes were found to be overexpressed.

TOP2A, DNA topoisomerase II alpha encoding a DNA topoisomerase, regulates the topologic states of DNA and controls tumor cell response. Wong et al. [26] showed that TOP2A overexpression in HCC tumors, relative to adjacent nontumors, correlated with histological grading, microvascular invasion, tumor deterioration, shorter survival, and chemoresistance, which was consistent with our clinicopathological features correlation analysis [27]. TOP2A targeted therapeutic drugs, such as doxorubicin, were suggested when TOP2A was overexpressed in HCC. Overexpression of microRNA-23a potentiated the response of HCC to the TOP2A targeted drugs etoposide and doxorubicin in vivo and in vitro, while suppressing another topoisomerase, TOP1, but not altering TOP2A expression levels [28]. Another study showed that TOP2A was targeted for proteasomal degradation by histone deacetylase inhibitors by activating casein kinase 2α and GSK3β double phosphorylation, highlighting a novel potential mechanism of HCC treatment [29]. The significance of TOP2A in tumor aggressiveness and chemoresistance indicate that the mechanisms through which it functions require further exploration.

ASPM encodes an abnormal spindle-like, microcephaly-associated protein that participates in mitosis and is reported to be a recurrence, invasion, metastasis, and prognostic marker for HCC. However, the specific molecular mechanism and pathways through which ASPM function have not been reported [30, 31]. CCNB2 encodes cyclin B2, the cell cycle mediated by transforming growth factor β, which, according to the GEPIA website, did not have a statistically significant influence on overall survival. However, we further analyzed the influence on CCNB2 on the UALCAN website [32] (http://ualcan.path.uab.edu/index.html) and found a statistical significance with . The variation in these results could be attributed to the different sample sizes and sequencing databases. A review of the literature revealed that CCNB2 expression is increased in HCC tissues compared to adjacent nontumor tissues [33], and CCNB2 is a target molecule following knockdown of XPOT [34], TPX2 [35], KPNA2 [36], and TMEM9 [37]. The prognostic value and target therapy of the hub genes for HCC patients remain to be confirmed by further clinical studies.

There were some limitations in our study. All of these sequencing data were from tissue samples and may not be suitable for unresectable patients whose liver tissues may be unavailable. Given that bioinformatics predictions are obtained in silico, further experimental validation is required. Additionally, we analyzed expression profiling by the array in the GEO datasets and concluded that TOP2A, ASPM, and CCNB2 were promising prognostic markers and potential therapeutic targets. However, other prognostic markers, such as methylation, microRNA, lncRNA, and circRNA, should be explored, and the data should be integrated with the findings of the present study to propose a new treatment algorithm that incorporates these markers in determining how aggressive is the treatment approach.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

This article does not contain any studies with human participants performed by any of the authors.

Conflicts of Interest

All authors declare that they have no conflict of interest.

Authors’ Contributions

Yuping Zeng and He He contributed equally to this work.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (81902142), Science and Technology Department of Sichuan Province (No. 2018SZ0271), and Health Department of Sichuan Provincial (No. 17PJ514).

Supplementary Materials

Supplementary Table 1: the top three enrichment pathways of five MCODE components. Supplementary Table 2: association between CCNB2, TOP2A, and ASPM expressions and other clinical features. (Supplementary Materials)

References

A. Forner, M. Reig, and J. Bruix, “Hepatocellular carcinoma,” The Lancet, vol. 391, no. 10127, pp. 1301–1314, 2018.
View at: Publisher Site | Google Scholar
R. Dhanasekaran, J.-C. Nault, L. R. Roberts, and J. Zucman-Rossi, “Genomic Medicine and Implications for Hepatocellular Carcinoma Prevention and Therapy,” Gastroenterology, vol. 156, no. 2, pp. 492–509, 2019.
View at: Publisher Site | Google Scholar
K. M. Kim, D. H. Sinn, S. H. Jung et al., “The recommended treatment algorithms of the BCLC and HKLC staging systems: does following these always improve survival rates for HCC patients?” Liver International, vol. 36, no. 10, pp. 1490–1497, 2016.
View at: Publisher Site | Google Scholar
H. Nakagawa, M. Fujita, and A. Fujimoto, “Genome sequencing analysis of liver cancer for precision medicine,” Seminars in Cancer Biology, vol. 55, pp. 120–127, 2019.
View at: Publisher Site | Google Scholar
A. Teufel, “Bioinformatics and database resources in hepatology,” Journal of Hepatology, vol. 62, no. 3, pp. 712–719, 2015.
View at: Publisher Site | Google Scholar
P. Xu, J. Yang, J. Liu et al., “Identification of glioblastoma gene prognosis modules based on weighted gene co-expression network analysis,” BMC Medical Genomics, vol. 11, no. 1, 2018.
View at: Publisher Site | Google Scholar
M. Giulietti, G. Occhipinti, A. Righetti et al., “Emerging Biomarkers in Bladder Cancer Identified by Network Analysis of Transcriptomic Data,” Frontiers in Oncology, vol. 8, 2018.
View at: Publisher Site | Google Scholar
Y. Wang, L. Chen, G. Wang et al., “Fifteen hub genes associated with progression and prognosis of clear cell renal cell carcinoma identified by co-expression analysis,” Journal of Cellular Physiology, vol. 234, no. 7, pp. 10225–10237, 2018.
View at: Publisher Site | Google Scholar
Z. Zhou, Y. Cheng, Y. Jiang et al., “Ten hub genes associated with progression and prognosis of pancreatic carcinoma identified by co-expression analysis,” International Journal of Biological Sciences, vol. 14, no. 2, pp. 124–136, 2018.
View at: Publisher Site | Google Scholar
L. Yuan, G. Qian, L. Chen et al., “Co-expression network analysis of biomarkers for adrenocortical carcinoma,” Frontiers in Genetics, vol. 9, 2018.
View at: Publisher Site | Google Scholar
J. Tang, D. Kong, Q. Cui et al., “Prognostic genes of breast cancer identified by gene co-expression network analysis,” Frontiers in Oncology, vol. 8, 2018.
View at: Publisher Site | Google Scholar
X. Xu, Y. Zhou, R. Miao et al., “Transcriptional modules related to hepatocellular carcinoma survival: coexpression network analysis,” Frontiers of Medicine, vol. 10, no. 2, pp. 183–190, 2016.
View at: Publisher Site | Google Scholar
W. Xu, Q. Rao, Y. An, M. Li, and Z. Zhang, “Identification of biomarkers for Barcelona Clinic Liver Cancer staging and overall survival of patients with hepatocellular carcinoma,” PLoS ONE, vol. 13, no. 8, article e0202763, 2018.
View at: Publisher Site | Google Scholar
E. Villa, R. Critelli, B. Lei et al., “Neoangiogenesis-related genes are hallmarks of fast-growing hepatocellular carcinomas and worst survival. Results from a prospective study,” Gut, vol. 65, no. 5, pp. 861–869, 2016.
View at: Publisher Site | Google Scholar
Y. H. Wang, T. Y. Cheng, T. Y. Chen, K. M. Chang, V. P. Chuang, and K. J. Kao, “Plasmalemmal vesicle associated protein (PLVAP) as a therapeutic target for treatment of hepatocellular carcinoma,” BMC Cancer, vol. 14, no. 1, 2014.
View at: Publisher Site | Google Scholar
Z. Makowska, T. Boldanova, D. Adametz et al., “Gene expression analysis of biopsy samples reveals critical limitations of transcriptome-based molecular classifications of hepatocellular carcinoma,” The Journal of Pathology: Clinical Research, vol. 2, no. 2, pp. 80–92, 2016.
View at: Publisher Site | Google Scholar
L. Wei, B. Lian, Y. Zhang et al., “Application of microRNA and mRNA expression profiling on prognostic biomarker discovery for hepatocellular carcinoma,” BMC Genomics, vol. 15, Suppl 1, 2014.
View at: Publisher Site | Google Scholar
S. Davis and P. S. Meltzer, “GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor,” Bioinformatics, vol. 23, no. 14, pp. 1846-1847, 2007.
View at: Publisher Site | Google Scholar
P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, 2008.
View at: Publisher Site | Google Scholar
P. Langfelder and S. Horvath, “FastRFunctions for Robust Correlations and Hierarchical Clustering,” Journal of statistical software, vol. 46, no. 11, 2012.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, article e47, 2015.
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
S. Tellapuri, P. D. Sutphin, M. S. Beg, A. G. Singal, and S. P. Kalva, “Staging systems of hepatocellular carcinoma: a review,” Indian Journal of Gastroenterology, vol. 37, no. 6, pp. 481–491, 2018.
View at: Publisher Site | Google Scholar
G. Bidkhori, R. Benfeitas, M. Klevstig et al., “Metabolic network-based stratification of hepatocellular carcinoma reveals three distinct tumor subtypes,” Proceedings of the National Academy of Sciences, vol. 115, no. 50, pp. E11874–E11883, 2018.
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
R. Panvichian, A. Tantiwetrueangdet, N. Angkathunyakul, and S. Leelaudomlipi, “TOP2A Amplification and Overexpression in Hepatocellular Carcinoma Tissues,” BioMed Research International, vol. 2015, Article ID 381602, 8 pages, 2015.
View at: Publisher Site | Google Scholar
N. Wang, M. Zhu, S. W. Tsao, K. Man, Z. Zhang, and Y. Feng, “MiR-23a-mediated inhibition of topoisomerase 1 expression potentiates cell response to etoposide in human hepatocellular carcinoma,” Molecular Cancer, vol. 12, no. 1, 2013.
View at: Publisher Site | Google Scholar
M. C. Chen, C. H. Chen, H. C. Chuang, S. K. Kulp, C. M. Teng, and C. S. Chen, “Novel mechanism by which histone deacetylase inhibitors facilitate topoisomerase IIα degradation in hepatocellular carcinoma cells,” Hepatology, vol. 53, no. 1, pp. 148–159, 2011.
View at: Publisher Site | Google Scholar
S. Y. Lin, H. W. Pan, S. H. Liu et al., “ASPM Is a Novel Marker for Vascular Invasion, Early Recurrence, and Poor Prognosis of Hepatocellular Carcinoma,” Clinical Cancer Research, vol. 14, no. 15, pp. 4814–4820, 2008.
View at: Publisher Site | Google Scholar
F. Wang, Y. Chang, J. Li et al., “Strong correlation between ASPM gene expression and HCV cirrhosis progression identified by co-expression analysis,” Digestive and Liver Disease, vol. 49, no. 1, pp. 70–76, 2017.
View at: Publisher Site | Google Scholar
D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.
View at: Publisher Site | Google Scholar
S. Liu, X. Yao, D. Zhang et al., “Analysis of transcription factor-related regulatory networks based on bioinformatics analysis and validation in hepatocellular carcinoma,” BioMed Research International, vol. 2018, Article ID 1431396, 16 pages, 2018.
View at: Publisher Site | Google Scholar
J. Lin, Y. Hou, S. Huang et al., “Exportin-T promotes tumor proliferation and invasion in hepatocellular carcinoma,” Molecular Carcinogenesis, vol. 58, no. 2, pp. 293–304, 2019.
View at: Publisher Site | Google Scholar
C.-W. Hsu, Y.-C. Chen, H.-H. Su et al., “Targeting TPX2 suppresses the tumorigenesis of hepatocellular carcinoma cells resulting in arrested mitotic phase progression and increased genomic instability,” Journal of Cancer, vol. 8, no. 8, pp. 1378–1394, 2017.
View at: Publisher Site | Google Scholar
H. Gao, A. Kerr, H. Jiao et al., “Long non-coding RNAs associated with metabolic traits in human white adipose tissue,” EBioMedicine, vol. 30, pp. 248–260, 2018.
View at: Publisher Site | Google Scholar
Y. Zhang, Y. Ran, Y. Xiong et al., “Effects of TMEM9 gene on cell progression in hepatocellular carcinoma by RNA interference,” Oncology Reports, vol. 36, no. 1, pp. 299–305, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Yuping Zeng 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

489

Downloads

1121

Citations