Baozhen Sun1
Xingyi Jin- 1Department of Hepatopancreatobiliary, China-Japan Union Hospital of Jilin University, Changchun, China
- 2First Department of Neurosurgery, China-Japan Union Hospital of Jilin University, Changchun, China
Despite studying the various molecular mechanisms of hepatocellular carcinoma (HCC), effective drugs and biomarkers in HCC therapy are still scarce. The present study was designed to investigate dysregulated pathways, novel biomarkers and therapeutic targets for HCC. The gene expression dataset of GSE14520, which included 362 tumor and their paired non-tumor tissues of HCC, was extracted for processing by the Robust multi-array average (RMA) algorithm in the R environment. SAM methods were leveraged to identify differentially expressed genes (DEGs). Functional analysis of DEGs was performed using DAVID. The GeneMania and Cytohubba were used to construct the PPI network. To avoid individual bias, GSEA and survival analysis were employed to verify the results. The results of these analyses indicated that separation of sister chromatids was the most aberrant phase in the progression of HCC, and the most frequently involved genes, EZH2, GINS1, TPX2, CENPF, and BUB1B, require further study to be used as drug targets or biomarkers in diagnosis and treatment of HCC.
Introduction
Hepatocellular carcinoma (HCC) is the third-leading cause of cancer-related deaths worldwide and its incidence continues to rise (Mittal and El-Serag, 2013). It mainly arise from hepatitis B virus (HBV) or hepatitis C virus (HCV) infections, and patients with cirrhosis have more opportunities to get HCC (El-Serag, 2002; Umemura et al., 2009). The limited knowledge on the molecular mechanisms of HCC contribute to poor prognosis and ineffective therapy, which leaves liver transplantation as the best choice of management (Ho et al., 2015; Turato et al., 2017). However, recurrence following transplantation is also associated with an unfavorable prognosis (Spinzi and Paggi, 2008; Sposito et al., 2013). Moreover, surgical intervention is ineffective in patients diagnosed at advanced stages of HCC (Rich et al., 2017). Therefore, new therapies for HCC are direly needed.
Although some molecular events that facilitate the progression to HCC have been investigated, the effective drug targets and potential biomarkers for early treatment and diagnosis of HCC are still unclear (Blum, 2005; Cha and Dematteo, 2005). Therefore, identifying the dysregulated pathways and hub genes involved in this process would allow us to identify patients with HCC as early as possible. Findings from previous studies that have focused on this area are limited due to small sample sizes, resulting in an incomplete understanding of HCC (Jia et al., 2007; Lin et al., 2014; Yin et al., 2017).
In this study, we extracted the gene expression profiles from the GEO database of 362 HBV-related HCC tumors and their paired non-tumor tissues which are mostly accompany with liver cirrhosis. The significance analysis of microarrays (SAMs) algorithm was used to screen the differential expressed genes (DEG), which was performed for pathway enrichment and generation of PPI network. After that, hub genes were identified by GeneMANIA and CytoHubba in Cytoscape. Furthermore, we performed a gene set enrichment analysis (GSEA), which evaluates microarray data at the level of gene sets, to overcome the limitation of individual gene analysis. In the meantime, survival analysis was leveraged, using the TCGA database, to assess the risks of hub gene expression. Finally, five most significantly hub genes were verified by qPCR and IHC in human HCC to confirm the results.
Results
Microarray Data Collection and Processing
The BioConductor package, Simpleaffy, was used for quality control and normalization of the microarray raw files (Supplementary Figure S1).
Identification of DEGs
To identify DEGs from the tumor and non-tumor tissue of HCC patients, we used the SAM method at the delta 14.31, with the FDR < 0.1%. A total of 862 genes were identified as DEGs, including 553 up-regulated and 309 down-regulated genes (Figure 1).
FIGURE 1. Significance analysis of microarrays (SAMs) plot comparing gene expression of HCC. 553 up-regulated genes are colored in red, 309 down-regulated genes are colored in green.
Functional Analysis of DEGs
Gene Ontology (GO) enrichment analysis revealed that GO terms were most significantly enriched in cell division (p = 1.09E-20), sister chromatid cohesion (p = 2.79E-13), mitotic nuclear division (p = 1.62E-11) and DNA replication (p = 1.34E-09; Figure 2). Reactome Pathway enrichment analysis showed that the separation of sister chromatids was the most significantly affected phase in HCC (Figure 3), which was in accordance with the results of the GO enrichment analysis.
FIGURE 2. Bubble plot for visualizing GO Annotation. Term are shown at –log (adj p-value) cut off >10. The z-score predicts existence of a bias in gene regulation.
FIGURE 3. The REACTOME pathway enrichment analysis of DEGs. 24 statistically significant pathways are listed and their colors are shown by p-values.
PPI Network Analysis of DEGs
A PPI network with 873 nodes and 84,272 edges was generated via the GeneMANIA plugin, around the DEGs. A global metric was utilized in the determination of hub proteins, through the Cytohubba plugin. Following this, the relationship between the 20 top-ranked proteins was mapped, based on the MCC as shown in Figure 4. The majority of these were cell division-related genes, such as TOP2A, GINS1, BUB1B, TPX2, and CENPF. The 20 top-ranked proteins were all up-regulated DEGs.
FIGURE 4. PPI network of 20 top-ranked DEGs. The hub genes are screened by cytoHubba and colored by their logFC.
GSEA and Leading-Edge Analysis
In order to further confirm the molecular mechanisms of HCC in the whole transcriptome, GSEA of the gene expression profile data from 362 tumor and paired non-tumor tissue of HCC patients was performed, based on the GO biological process. The results revealed that the most significant biological processes that were enriched were cell division-related processes, including GO_REGULATION_OF_ CELL_DIVISION (FDR = 0.052), GO_REGULATION_ OF_MEIOTIC _CELL_CYCLE (FDR = 0.050), GO_REGULATION_ OF_NUCLEAR_DIVISION (FDR = 0.051), GO_SISTER_ CHROMATID_COHESION (FDR = 0.051), GO_REGULATION_OF_ CHROMOSOME_ORGANIZATION (FDR = 0.043), and GO_SISTER_ CHROMATID_SEGREGATION (FDR = 0.033), which are shown in Figure 5A. Leading-edge analysis was used to find the hub genes which appeared frequently in related gene sets, and the results showed that AURKB appeared in five gene sets, while BUB1B, CDC20, FBXO5, DLGAP5, ESPL1, BIRC5, BUB3, BUB1, CENPE, CENPF, and MAD2L1 appeared in four gene sets (Figures 5B,C).
FIGURE 5. GSEA of gene expression profiling of HCC. (A) GO enrichment plot of six cell division related gene sets. The normalized enrichment score (NES) and the false discovery rates (FDR) are indicated for each gene set. Each bar at the bottom of plot represents a member gene of the respective gene set. (B) Leading edge analysis among gene sets. Color intensity is used to show the overlap between subsets, the darker the color, the greater the overlap between the subsets. Specifically, the intensity of the cell for sets A and B corresponds to an X/Y ratio where X is the number of leading edge genes from set A and Y is the union of leading edge genes in sets A and B. (C) The bar graph shows each gene and the number of subsets in which it appears.
Kaplan–Meier Survival Analysis
SurvExpress was engaged to explore the relationship between the hub genes and the survival of HCC patients in silico. Finally, survival analysis, based on clinical information from the TCGA-liver cancer datasets, revealed that the high expression of EZH2, GINS1, and TPX2 correlated positively with higher risk, CENPF and BUB1B were quite the contrary (Figure 6).
FIGURE 6. Kaplan–Meier survival curve of 422 TCGA liver cancer samples using the SurvExpress database, based on the low or high risk for a poor outcome. Censoring samples are marked with “+.” The Horizontal axis represents time to event and vertical axis represents percentage. High- and low-risk groups are labeled with red and green curves, respectively. (A) High expression of five hub genes is correlated with high risk, poor prognosis and shorter overall survival time. (B–F) Survival curve of five hub genes respectively.
Validation of Hub Genes by qPCR in Human HCC
Five most significantly hub genes, EZH2, GINS1, TPX2, CENPF, and BUB1B, were successfully validated by qPCR in 30 paired human HCC tissues which have no difference with the analysis results of gene expression profiling (Figure 7).
FIGURE 7. Validation of five hub genes by qPCR. HCC sample (red) and their paired non-tumor tissue (green) were validated by qPCR, ∗∗∗∗P < 0.0001.
Immunohistochemistry
IHC was employed to validate the results from bioinformation analysis, which revealed the strong expression of five hub genes in HCC vs. the control group. However, CENPF and BUB1B were not outstandingly in our result, the individual difference may affect the outcome to some extent (Figure 8).
FIGURE 8. Validation of five hub genes by IHC. Representative photomicrographs of immunohistochemical detection and measurement of five hub genes by Image Pro Plus 6.0 between 12 paired tumor and adjacent non-tumor tissues, ∗P < 0.05 and ∗∗∗∗P < 0.0001 (magnification 200×).
Discussion
In the present study, we explored the molecular mechanisms of HCC between tumor and non-tumor tissues, using bioinformatics analysis. Our results indicated that the separation of sister chromatids was the most significantly dysregulated pathway during the transition from cirrhosis to HCC, with the up-regulation of 12 hub genes.
One of the most miraculous events in the human cell cycle is the concurrent separation of 46 pairs of sister chromatids. Since this irreversible separation is highly monitored and regulated, neither damage to the genome, nor errors in chromosome alignment, can be easily rectified after separation (Nasmyth, 1999). The initiation of chromosomal segregation was supervised by the spindle assembly checkpoint (SAC), which ensures the genomic stability during mitosis (Dai, 2009). Defects in sister chromatid segregation could lead to aneuploidy (Panigrahi and Pati, 2009) and promote chromosome instability (CIN) during mitosis (Schvartzman et al., 2010). In addition, CIN may contribute to tumor initiation and/or progression, which has been demonstrated in cell lines (Zhang et al., 2004; Vader and Lens, 2008), mice models (Bernal et al., 2002; Dawlaty et al., 2008), and human tumors (Kronenwett et al., 2004; Carter et al., 2006). In the Sgo1± mouse model, mitotic error-induced CIN was shown to be an important early event in HCC development (Kronenwett et al., 2004). Furthermore, 120 HCC with 195 markers (Nagai et al., 1997), and 48 HCC with 275 chromosomal markers (Boige et al., 1997), indicated that CIN appears widely in HCC. Aneuploidy is a major manifestation of CIN and is seen in over 75% of cancers, and is also considered essential for tumorigenesis, by some biologists (Duesberg and Li, 2003). Therefore, dysregulation of sister chromatids separation might contribute to the initiation and progression of human HCC.
Based on the results of the PPI analysis by GeneMania and Cytohubba, 20 top-ranked proteins from 862 DEGs are thought to participate in the core pathway of HCC, such as TOP2, GINS1, EZH2, TTK2, CDK1, BUB1B, TPX2, CENPF, and MAD2L1. Furthermore, the high expression of GINS1, EZH2, and TPX2, correlates with high-risk in HCC, as confirmed by the survival analysis. BUB1B, TPX2, and CENPF appeared most often in related gene sets, which were demonstrated by GSEA at the level of gene sets. Therefore, EZH2, GINS1, TPX2, CENPF, and BUB1B are thought to be hub genes in HCC and are discussed below.
Epigenetic regulation of gene expression, particularly hypermethylation, plays an important role in tumorigenesis (Jones and Baylin, 2002). EZH2 is the enzymatic subunit of Polycomb repressive complex 2 (PRC2), which methylates H3K27, resulting in silence of the associated tumor suppressor genes (Margueron and Reinberg, 2011; Momparler et al., 2012). EZH2 was up-regulated and expressed in many solid cancers, and YY1 can recruit EZH2 and suppress NFkB function in hepatitis B virus-dependent HCC (Chase and Cross, 2011). Additionally, EZH2 is clinically associated with tumor progression and multiple metastatic features, and epigenetically restrained a subset of miRNA in human HCC (Au et al., 2012). Thus, EZH2 may be regarded as a potential therapeutic target, and a few of compounds have been already investigated as inhibitors of EZH2, in pre-clinical studies (Knutson et al., 2012, 2013; McCabe et al., 2012; Qi et al., 2012; Kim et al., 2013).
GINS1/PSF1 is a subunit of the GINS complex, which is involved in the DNA replication fork and the initiation of chromosome replication (Labib and Gambus, 2007). Research suggests that GINS1 and/or other GINS complex subunits are upregulated in some types of cancers and possess some tumorigenic characteristics (Hokka et al., 2013; Zhang et al., 2015; Zhou et al., 2015). GINS1 is expressed at high levels in HCC tissues, which is associated with more aggressive tumors and worse prognosis. Moreover, in a mouse xenograft model, high levels of GINS1 expression correspond to high proliferative activity, transplantation potential, and metastatic capability (Nagahama et al., 2010). In contrast, knockdown of GINS1 expression led to inhibited tumor growth by disrupting DNA replication and chromosomal segregation, and promoted apoptosis, particularly early apoptosis (Nagahama et al., 2010; Zhou et al., 2015). These findings may make GINS1 a potential theranostic target in the future.
For many years, the function of TPX2 has been studied in mitosis and spindle assembly because of the chromatin-mediated TPX2/Importinα-β/Ran signal and its control of Aurora A kinase (Asteriti et al., 2010). The location of TPX2 is at the long arm of chromosome 20, at position 20q11, which is often amplified in HCC and other tumors (Knuutila et al., 1998; Hodgson et al., 2003; Scotto et al., 2008; Beroukhim et al., 2010). The elevated TPX2 expression results in dysregulation of spindle formation and balanced chromosome segregation, by over activation of Aurora-A, which could lead to unscheduled phosphorylation of downstream targets. TPX2 knockdown inhibits cell proliferation and AKT signaling, and decreases the MMP2 and MMP9 expression in HCC cell lines (Liu et al., 2014). Clinical sample analysis also indicates that TPX2 expression is associated with the tumor–node–metastasis stage, tumor numbers, and tumor differentiation in the HCC tissues (Liu et al., 2015). Moreover, TPX2 inactivation experiments indicated anti-proliferative effects in cancer cells, suggesting the potential value of TPX2 as an anti-cancer target (Fenner et al., 2005; Warner et al., 2009; Li et al., 2010).
CENPF is a large coiled-coil protein whose expression and subcellular localization was cell cycle-dependent, and undergoes rapid degradation at the end of cell division. It is reported that CENPF plays a major role in kinetochore assembly, regulation of chromosome segregation, and control of SAC activity. CENPF may be a potential proliferation marker in the clinical diagnosis of HCC (Ma et al., 2006). CENPF is overexpressed in HCC (Kim et al., 2012) and other tumors (Varis et al., 2006). Additionally, the up-regulated CENPF expression has been shown to contribute to the proliferation of HCC rather than acting as a trigger for malignant cell growth. Consequently, CENPF could be an indicator of tumorigenesis, especially at early stages of HCC (Zhang et al., 2001).
BUB1B is a key component in the SAC protein family, which has been proven to be involved and upregulated in multiple human cancers (Seike et al., 2002; Shichiri et al., 2002; Gupta et al., 2003; Yamamoto et al., 2007; Fu et al., 2016). In mitosis, BUB1B accumulates cyclin B in G2 phase, by binding to CDC20 to inhibit APC/C activity and prolonging the checkpoint signaling by kinase activity at kinetochores (Malureanu et al., 2009). A human study reported that the role of BUB1B was to facilitate accurate chromosome segregation and maintain chromosomal stability, to suppress cancer (Hanks et al., 2004). However, the phosphorylated BUB1B, which is tightly regulated through its own activation and subcellular localization (Bin et al., 1998; Li et al., 1999), was elevated in the SV40 Tag-derived prostate cancer models (Guo et al., 2006). The contradiction of BUB1B, between its role in suppressing cancer and upregulating cancers, is demonstrated in varying reports of cancer-associated missense and nonsense mutations in BUB1B, in several cancers (Cahill et al., 1998, 1999; Myrie et al., 2000; Saeki et al., 2002). However, the true mechanism of BUB1B in cancers remains to be elucidated.
In conclusion, we collectively analyzed the molecular mechanisms of human HCC through interpretation of the functions and PPI network of DEGs, which were confirmed by GSEA and survival analysis. In doing so, we ascertained the molecular genetic differences between tumor and non-tumor tissues of HCC, which suggest that separation of sister chromatids may have the most important influence on initiation and progression of human HCC. Errors in this process contribute to CIN and aneuploidy which were thought to be responsible for tumorigenic progression in human cells. The hub genes we found may be useful as biomarkers for diagnosis and prognosis or in tailoring treatment in human HCC. Finally, additional studies are needed to confirm the findings of these experiments.
Materials and Methods
Microarray Data Collection and Processing
The gene expression dataset of GSE14520 was obtained from NCBI GEO database1 which is based on the Affymetrix Human Genome U133A 2.0 Array. 362 tumor and paired non-tumor tissues of HCC patients were collected for analysis of genome microarrays. The detail of every HCC patient was shown in Supplementary Table S1. The Simpleaffy package was used to read the raw data and perform quality control and normalization by Robust multi-array average (RMA) algorithm in the R environment (Wilson and Miller, 2005). The mean gene expression was considered in multiple probe sets with one name.
Identification of Differentially Expressed Genes (DEGs)
The DEGs in tumor and non-tumor tissues of HCC patients were determined using SAM (Grace and Nacheva, 2012). We used two class unpaired comparison analysis with t-statistics, and permutations of 103. SAM uses permutations of repeated measurements to estimate the percentage of genes identified by chance, and the false discovery rate (FDR).
Functional Analysis of DEGs
Gene Ontology (GO) and Reactome Pathway enrichment analyses were performed using DAVID to explore the biological processes and signaling pathways in which the DEGs were involved (Ashburner et al., 2000; Huang et al., 2009; Fabregat et al., 2018). The enrichment results were visualized in a new R visualization package called GOPlot (Walter et al., 2015).
PPI Network Analysis of DEGs and the Screening of Hub Proteins
A PPI network analysis was performed to evaluate physical relationships between the proteins encoded by the DEGs. The GeneMania (Montojo et al., 2010) and Cytohubba (Chin et al., 2014) were used to construct the PPI network based on co-expression, physical interactions, and genetic interactions.
Gene Set Enrichment Analysis (GSEA) and Leading-Edge Analysis
To further investigate the biological characteristics of HCC, we performed GSEA assay in non-tumor and tumor groups with permutations of 104 in the GO biological process, from the Molecular Signature Database (MSigDB). A leading-edge analysis was performed to elucidate hub genes of HCC according to the results of Reactome Pathway enrichment (Subramanian et al., 2005).
Kaplan–Meier Survival Analysis
Kaplan–Meier analysis was performed with the online multi-cancer biomarker validation tool, SurvExpress, in the TCGA-liver cancer datasets containing 422 samples, using the hub genes as an input (Aguirre-Gamboa et al., 2013).
Sample Collection
Thirty HBV-related HCC and their paired non-tumor tissues were collected from August 2014 to December 2017 at China-Japan Union Hospital of Jilin University. All the specimens were dealing with liquid nitrogen after surgical resection and stored at -80°C. Three independent pathologists made the decision about the diagnosis of HCC and assessed the samples with HE staining.
Quantitative RT-PCR
Total RNA of 30 HCC samples was extracted using Trizol (Invitrogen) as described everywhere, qRT-PCR was performed by One-Step qPCR Kit (Invitrogen) and CFX ConnectTM Real-Time System (BIO-RAD) following manufacturer’s instructions. The data of qPCR were processed by ΔΔCt method with normalizing to GAPDH as the reference gene. The sequence of primers was shown in Supplementary Table S2.
Immunohistochemistry
The 4 μm thick sections cut from formalin-fixed, paraffin-embedded HCC tissue were used for IHC of hub genes as described before (Ho et al., 2012). The primary antibodies against EZH2 (No. 191080, Abcam), BUB1B (No. 183496, Abcam), CENPF (No. 223847, Abcam), TPX2 (11741-1-AP, Proteintech) and GINS1 (No. D161403, Sangon Biotech) were employed for IHC. Image Pro Plus 6.0 (Media Cybernetics, Bethesda, MD, United States) was employed to measure the positive area of hub genes for quantitative analysis.
Statistical Analysis
All data are shown as mean ± SDs except for otherwise indicated. Significance was determined with two-tailed t-test when comparing the variance from HCC to the adjacent non-tumor tissue. GraphPad Prism 6 software (GraphPad Software, La Jolla, CA, United States) was used for analysis. A P < 0.05 is considered significant.
Ethics Statement
This study was carried out in accordance with the recommendations of CIOMS. The protocol was approved by the institutional review boards of the China-Japan Union Hospital of Jilin University. All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Author Contributions
XJ conceived and designed the study. BS collected, analyzed the data, and wrote the manuscript. GL, DJ, SL, and GC revised the manuscript. All authors read and approved the manuscript.
Conflict of Interest Statement
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/fphys.2018.01019/full#supplementary-material
FIGURE S1 |(A) Data before normalization. (B) Data after normalization.
TABLE S1 | Samples detail.
TABLE S2 | qPCR primers used for tissues specific expression analysis.
Footnotes
References
Aguirre-Gamboa, R., Gomez-Rueda, H., Martinez-Ledesma, E., Martinez-Torteya, A., Chacolla-Huaringa, R., Rodriguez-Barrientos, A., et al. (2013). SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS One 8:e74250. doi: 10.1371/journal.pone.0074250
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. the gene ontology consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Asteriti, I. A., Rensen, W. M., Lindon, C., Lavia, P., and Guarguaglini, G. (2010). The Aurora-A/TPX2 complex: a novel oncogenic holoenzyme? Biochim. Biophys. Acta Rev. Cancer 1806, 230–239. doi: 10.1016/j.bbcan.2010.08.001
Au, S. L., Wong, C. C., Lee, J. M., Fan, D. N., Tsang, F. H., Ng, I. O., et al. (2012). Enhancer of zeste homolog 2 epigenetically silences multiple tumor suppressor microRNAs to promote liver cancer metastasis. Hepatology 56, 622–631. doi: 10.1002/hep.25679
Bernal, J. A., Luna, R., Espina, A., Lazaro, I., Ramos-Morales, F., Romero, F., et al. (2002). Human securin interacts with p53 and modulates p53-mediated transcriptional activity and apoptosis. Nat. Genet. 32, 306–311. doi: 10.1038/ng997
Beroukhim, R., Mermel, C. H., Porter, D., Wei, G., Raychaudhuri, S., Donovan, J., et al. (2010). The landscape of somatic copy-number alteration across human cancers. Nature 463, 899–905. doi: 10.1038/nature08822
Bin, O. Y., Lan, Z. D., Meadows, J., Pan, H. Q., Fukasawa, K., Li, W. Q., et al. (1998). Human Bub1: a putative spindle checkpoint kinase closely linked to cell proliferation. Cell Growth Diff. 9, 877–885.
Blum, H. E. (2005). Treatment of hepatocellular carcinoma. Best Pract. Res. Clin. Gastroenterol. 19, 129–145. doi: 10.1016/j.bpg.2004.11.008
Boige, V., Laurentpuig, P., Fouchet, P., Flejou, J. F., Monges, G., Bedossa, P., et al. (1997). Concerted nonsyntenic allelic losses in hyperploid hepatocellular carcinoma as determined by a high-resolution allelotype. Cancer Res. 57, 1986–1990.
Cahill, D. P., Da Costa, L. T., Carson-Walter, E. B., Kinzler, K. W., Vogelstein, B., and Lengauer, C. (1999). Characterization of MAD2B and other mitotic spindle checkpoint genes. Genomics 58, 181–187. doi: 10.1006/geno.1999.5831
Cahill, D. P., Lengauer, C., Yu, J., Riggins, G. J., Willson, J. K. V., Markowitz, S. D., et al. (1998). Mutations of mitotic checkpoint genes in human cancers. Nature 392, 300–303. doi: 10.1038/32688
Carter, S. L., Eklund, A. C., Kohane, I. S., Harris, L. N., and Szallasi, Z. (2006). A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers. Nat. Genet. 38, 1043–1048. doi: 10.1038/ng1861
Cha, C., and Dematteo, R. R. (2005). Molecular mechanisms in hepatocellular carcinoma development. Best Pract. Res. Clin. Gastroenterol. 19, 25–37. doi: 10.1016/j.bpg.2004.11.005
Chase, A., and Cross, N. C. P. (2011). Aberrations of EZH2 in Cancer. Clin. Cancer Res. 17, 2613–2618. doi: 10.1158/1078-0432.Ccr-10-2156
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 8(Suppl. 4):S11. doi: 10.1186/1752-0509-8-S4-S11
Dai, W. (2009). Suppression of genomic instabilities caused by chromosome Mis-segregation: a perspective from studying BubR1 and Sgo1. J. Formos. Med. Assoc. 108, 904–911. doi: 10.1016/s0929-6646(10)60002-2
Dawlaty, M. M., Malureanu, L., Jeganathan, K. B., Kao, E., Sustmann, C., Tahk, S., et al. (2008). Resolution of sister centromeres requires RanBP2-mediated SUMOylation of topoisomerase II alpha. Cell 133, 103–115. doi: 10.1016/j.cell.2008.01.045
Duesberg, P., and Li, R. (2003). Multistep carcinogenesis A chain reaction of aneuploidizations. Cell Cycle 2, 202–210. doi: 10.4161/cc.2.3.382
El-Serag, H. B. (2002). Hepatocellular carcinoma: an epidemiologic view. J. Clin. Gastroenterol. 35, S72–S78. doi: 10.1097/00004836-200211002-00002
Fabregat, A., Jupe, S., Matthews, L., Sidiropoulos, K., Gillespie, M., Garapati, P., et al. (2018). The reactome pathway knowledgebase. Nucleic Acids Res. 46, D649–D655. doi: 10.1093/nar/gkx1132
Fenner, M., Wehrhan, F., Jehle, M., Amann, K., Radespiel-Troger, M., Grabenbauer, G., et al. (2005). Restricted-expressed proliferation-associated protein (repp86) expression in squamous cell carcinoma of the oral cavity. Strahlentherapie Und Onkologie 181, 755–761. doi: 10.1007/s00066-005-1430-7
Fu, X., Chen, G., Cai, Z.-D., Wang, C., Liu, Z.-Z., Lin, Z.-Y., et al. (2016). Overexpression of BUB1B contributes to progression of prostate cancer and predicts poor outcome in patients with prostate cancer. Oncotargets Ther. 9, 2211–2220. doi: 10.2147/ott.S101994
Grace, C., and Nacheva, E. P. (2012). Significance analysis of microarrays (SAM) offers clues to differences between the genomes of adult philadelphia positive all and the lymphoid blast transformation of CML. Cancer Inform. 11, 173–183. doi: 10.4137/CIN.S9258
Guo, C., Wu, G., Chin, J. L., Bauman, G., Moussa, M., Wang, F., et al. (2006). Bub1 up-regulation and hyperphosphorylation promote malignant transformation in SV40 tag-induced transgenic mouse models. Mol. Cancer Res. 4, 957–969. doi: 10.1158/1541-7786.MCR-06-0168
Gupta, A., Inaba, S., Wong, O. K., Fang, G. W., and Liu, J. W. (2003). Breast cancer-specific gene 1 interacts with the mitotic checkpoint kinase BubR1. Oncogene 22, 7593–7599. doi: 10.1038/sj.onc.1206880
Hanks, S., Coleman, K., Reid, S., Plaja, A., Firth, H., Fitzpatrick, D., et al. (2004). Constitutional aneuploidy and cancer predisposition caused by biallelic mutations in BUB1B. Nat. Genet. 36, 1159–1161. doi: 10.1038/ng1449
Ho, D. W. H., Kai, A. K. L., and Ng, I. O. L. (2015). TCGA whole-transcriptome sequencing data reveals significantly dysregulated genes and signaling pathways in hepatocellular carcinoma. Front. Med. 9, 322–330. doi: 10.1007/s11684-015-0408-9
Ho, D. W. Y., Yang, Z. F., Yi, K., Lam, C. T., Ng, M. N. P., Yu, W. C., et al. (2012). Gene expression profiling of liver cancer stem cells by RNA-sequencing. PLoS One 7:e37159. doi: 10.1371/journal.pone.0037159
Hodgson, J. G., Chin, K., Collins, C., and Gray, J. W. (2003). Genome amplification of chromosome 20 in breast cancer. Breast Cancer Res. Treat. 78, 337–345. doi: 10.1023/A:1023085825042
Hokka, D., Maniwa, Y., Tane, S., Nishio, W., Yoshimura, M., Okita, Y., et al. (2013). Psf3 is a prognostic biomarker in lung adenocarcinoma. Lung Cancer 79, 77–82. doi: 10.1016/j.lungcan.2012.10.001
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Jia, H. L., Ye, Q. H., Qin, L. X., Budhu, A., Forgues, M., Chen, Y., et al. (2007). Gene expression profiling reveals potential biomarkers of human hepatocellular carcinoma. Clin. Cancer Res. 13, 1133–1139. doi: 10.1158/1078-0432.CCR-06-1025
Jones, P. A., and Baylin, S. B. (2002). The fundamental role of epigenetic events in cancer. Nat. Rev. Genet. 3, 415–428. doi: 10.1038/nrg816
Kim, H.-E., Kim, D.-G., Lee, K. J., Son, J. G., Song, M.-Y., Park, Y.-M., et al. (2012). Frequent amplification of CENPF, GMNN and CDK13 genes in hepatocellular carcinomas. PLoS One 7:e43223. doi: 10.1371/journal.pone.0043223
Kim, W., Bird, G. H., Neff, T., Guo, G. J., Kerenyi, M. A., Walensky, L. D., et al. (2013). Targeted disruption of the EZH2-EED complex inhibits EZH2-dependent cancer. Nat. Chem. Biol. 9, 643–650. doi: 10.1038/Nchembio.1331
Knutson, S. K., Warholic, N. M., Wigle, T. J., Klaus, C. R., Allain, C. J., Raimondi, A., et al. (2013). Durable tumor regression in genetically altered malignant rhabdoid tumors by inhibition of methyltransferase EZH2. Proc. Natl. Acad. Sci. U.S.A. 110, 7922–7927. doi: 10.1073/pnas.1303800110
Knutson, S. K., Wigle, T. J., Warholic, N. M., Sneeringer, C. J., Allain, C. J., Klaus, C. R., et al. (2012). A selective inhibitor of EZH2 blocks H3K27 methylation and kills mutant lymphoma cells. Nat. Chem. Biol. 8, 890–896. doi: 10.1038/Nchembio.1084
Knuutila, S., Bjorkqvist, A. M., Autio, K., Tarkkanen, M., Wolf, M., Monni, O., et al. (1998). DNA copy number amplifications in human neoplasms - Review of comparative genomic hybridization studies. Am. J. Pathol. 152, 1107–1123.
Kronenwett, U., Huwendiek, S., Ostring, C., Portwood, N., Roblick, U. J., Pawitan, Y., et al. (2004). Improved grading of breast adenocarcinomas based on genomic instability. Cancer Res. 64, 904–909. doi: 10.1158/0008-5472.can-03-2451
Labib, K., and Gambus, A. (2007). A key role for the GINS complex at DNA replication forks. Trends Cell Biol. 17, 271–278. doi: 10.1016/j.tcb.2007.04.002
Li, B., Qi, X. Q., Chen, X., Huang, X., Liu, G. Y., Chen, H. R., et al. (2010). Expression of targeting protein for Xenopus kinesin-like protein 2 is associated with progression of human malignant astrocytoma. Brain Res. 1352, 200–207. doi: 10.1016/j.brainres.2010.06.060
Li, W. Q., Lan, Z. D., Wu, H. Y., Wu, S. C., Meadows, J., Chen, J., et al. (1999). BUBR1 phosphorylation is regulated during mitotic checkpoint activation. Cell Growth Diff. 10, 769–775.
Lin, K. T., Shann, Y. J., Chau, G. Y., Hsu, C. N., and Huang, C. Y. F. (2014). Identification of latent biomarkers in hepatocellular carcinoma by ultra-deep whole-transcriptome sequencing. Oncogene 33, 4786–4794. doi: 10.1038/onc.2013.424
Liu, Q., Tu, K., Zhang, H., Zheng, X., Yao, Y., and Liu, Q. (2015). TPX2 as a novel prognostic biomarker for hepatocellular carcinoma. Hepatol. Res. 45, 906–918. doi: 10.1111/hepr.12428
Liu, Q., Yang, P., Tu, K., Zhang, H., Zheng, X., Yao, Y., et al. (2014). TPX2 knockdown suppressed hepatocellular carcinoma cell invasion via inactivating AKT signaling and inhibiting MMP2 and MMP9 expression. Chin. J. Cancer Res. 26, 410–417. doi: 10.3978/j.issn.1000-9604.2014.08.01
Ma, L., Zhao, X. S., and Zhu, X. L. (2006). Mitosin/CENP-F in mitosis, transcriptional control, and differentiation. J. Biomed. Sci. 13, 205–213. doi: 10.1007/s11373-005-9057-3
Malureanu, L. A., Jeganathan, K. B., Hamada, M., Wasilewski, L., Davenport, J., and Van Deursen, J. M. (2009). BubR1 N terminus acts as a soluble inhibitor of cyclin B Degradation by APC/C-Cdc20 in interphase. Dev. Cell 16, 118–131. doi: 10.1016/j.devcel.2008.11.004
Margueron, R., and Reinberg, D. (2011). The Polycomb complex PRC2 and its mark in life. Nature 469, 343–349. doi: 10.1038/nature09784
McCabe, M. T., Ott, H. M., Ganji, G., Korenchuk, S., Thompson, C., Van Aller, G. S., et al. (2012). EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. Nature 492, 108–112. doi: 10.1038/nature11606
Mittal, S., and El-Serag, H. B. (2013). Epidemiology of hepatocellular carcinoma: consider the population. J. Clin. Gastroenterol. 47(Suppl.), S2–S6. doi: 10.1097/MCG.0b013e3182872f29
Momparler, R. L., Idaghdour, Y., Marquez, V. E., and Momparler, L. F. (2012). Synergistic antileukemic action of a combination of inhibitors of DNA methylation and histone methylation. Leuk. Res. 36, 1049–1054. doi: 10.1016/j.leukres.2012.03.001
Montojo, J., Zuberi, K., Rodriguez, H., Kazi, F., Wright, G., Donaldson, S. L., et al. (2010). GeneMANIA Cytoscape plugin: fast gene function predictions on the desktop. Bioinformatics 26, 2927–2928. doi: 10.1093/bioinformatics/btq562
Myrie, K. A., Percy, M. J., Azim, J. N., Neeley, C. K., and Petty, E. M. (2000). Mutation and expression analysis of human BUB1 and BUB1B in aneuploid breast cancer cell lines. Cancer Lett. 152, 193–199. doi: 10.1016/S0304-3835(00)00340-2
Nagahama, Y., Ueno, M., Miyamoto, S., Morii, E., Minami, T., Mochizuki, N., et al. (2010). PSF1, a DNA replication factor expressed widely in stem and progenitor cells, drives tumorigenic and metastatic properties. Cancer Res. 70, 1215–1224. doi: 10.1158/0008-5472.Can-09-3662
Nagai, H., Pineau, P., Tiollais, P., Buendia, M. A., and Dejean, A. (1997). Comprehensive allelotyping of human hepatocellular carcinoma. Oncogene 14, 2927–2933. doi: 10.1038/sj.onc.1201136
Nasmyth, K. (1999). Separating sister chromatids. Trends Biochem. Sci. 24, 98–104. doi: 10.1016/s0968-0004(99)01358-4
Panigrahi, A. K., and Pati, D. (2009). Road to the crossroads of life and death: linking sister chromatid cohesion and separation to aneuploidy, apoptosis and cancer. Crit. Rev. Oncol. Hematol. 72, 181–193. doi: 10.1016/j.critrevonc.2008.12.002
Qi, W., Chan, H. M., Teng, L., Li, L., Chuai, S. N., Zhang, R. P., et al. (2012). Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation. Proc. Natl. Acad. Sci. U.S.A. 109, 21360–21365. doi: 10.1073/pnas.1210371110
Rich, N. E., Parikh, N. D., and Singal, A. G. (2017). Overdiagnosis: an understudied issue in hepatocellular carcinoma surveillance. Semin. Liver Dis. 37, 296–304. doi: 10.1055/s-0037-1608775
Saeki, A., Tamura, S., Ito, N., Kiso, S., Matsuda, Y., Yabuuchi, I., et al. (2002). Frequent impairment of the spindle assembly checkpoint in hepatocellular carcinoma. Cancer 94, 2047–2054. doi: 10.1002/cncr.10448
Schvartzman, J.-M., Sotillo, R., and Benezra, R. (2010). Mitotic chromosomal instability and cancer: mouse modelling of the human disease. Nat. Rev. Cancer 10, 102–115. doi: 10.1038/nrc2781
Scotto, L., Narayan, G., Nandula, S., Arias-Pulido, H., Subramaniyam, S., Schneider, A., et al. (2008). Identification of copy number gain and overexpressed genes on chromosome arm 20q by an integrative genomic approach in cervical cancer: potential role in progression. Genes Chromosomes Cancer 47, 755–765. doi: 10.1002/gcc.20577
Seike, M., Gemma, A., Hosoya, Y., Hosomi, Y., Okano, T., Kurimoto, F., et al. (2002). The promoter region of the human BUBR1 gene and its expression analysis in lung cancer. Lung Cancer 38, 229–234. doi: 10.1016/S0169-5002(02)00218-0
Shichiri, M., Yoshinaga, K., Hisatomi, H., Sugihara, K., and Hirata, Y. (2002). Genetic and epigenetic inactivation of mitotic checkpoint genes hBUB1 and hBUBR1 and their relationship to survival. Cancer Res. 62, 13–17.
Spinzi, G., and Paggi, S. (2008). Sorafenib in advanced hepatocellular carcinoma. N. Engl. J. Med. 359, 2497–2498. doi: 10.1056/NEJMc081780
Sposito, C., Mariani, L., Germini, A., Flores Reyes, M., Bongini, M., Grossi, G., et al. (2013). Comparative efficacy of sorafenib versus best supportive care in recurrent hepatocellular carcinoma after liver transplantation: a case-control study. J. Hepatol. 59, 59–66. doi: 10.1016/j.jhep.2013.02.026
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjeed, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Turato, C., Balasso, A., Carloni, V., Tiribelli, C., Mastrotto, F., Mazzocca, A., et al. (2017). New molecular targets for functionalized nanosized drug delivery systems in personalized therapy for hepatocellular carcinoma. J. Control. Release 268, 184–197. doi: 10.1016/j.jconrel.2017.10.027
Umemura, T., Ichijo, T., Yoshizawa, K., Tanaka, E., and Kiyosawa, K. (2009). Epidemiology of hepatocellular carcinoma in Japan. J. Gastroenterol. 44(Suppl. 19), 102–107. doi: 10.1007/s00535-008-2251-0
Vader, G., and Lens, S. M. A. (2008). The Aurora kinase family in cell division and cancer. Biochim. Biophys. Acta Rev. Cancer 1786, 60–72. doi: 10.1016/j.bbcan.2008.07.003
Varis, A., Salmela, A.-L., and Kallio, M. J. (2006). Cenp-F (mitosin) is more than a mitotic marker. Chromosoma 115, 288–295. doi: 10.1007/s00412-005-0046-0
Walter, W., Sanchez-Cabo, F., and Ricote, M. (2015). GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics 31, 2912–2914. doi: 10.1093/bioinformatics/btv300
Warner, S. L., Stephens, B. J., Nwokenkwo, S., Hostetter, G., Sugeng, A., Hidalgo, M., et al. (2009). Validation of TPX2 as a potential therapeutic target in pancreatic cancer cells. Clin. Cancer Res. 15, 6519–6528. doi: 10.1158/1078-0432.Ccr-09-0077
Wilson, C., and Miller, C. (2005). Simpleaffy: a BioConductor package for affymetrix quality control and data analysis. Bioinformatics 21, 3683–3685. doi: 10.1093/bioinformatics/bti605
Yamamoto, Y., Matsuyama, H., Chochi, Y., Okuda, M., Kawauchi, S., Inoue, R., et al. (2007). Overexpression of BUBR1 is associated with chromosomal instability in bladder cancer. Cancer Genet. Cytogenet. 174, 42–47. doi: 10.1016/j.cancergencyto.2006.11.012
Yin, L., Chang, C., and Xu, C. (2017). G2/M checkpoint plays a vital role at the early stage of HCC by analysis of key pathways and genes. Oncotarget 8, 76305–76317. doi: 10.18632/oncotarget.19351
Zhang, D. W., Hirota, T., Marumoto, T., Shimizu, M., Kunitoku, N., Sasayama, T., et al. (2004). Cre-loxP-controlled periodic Aurora-A overexpression induces mitotic abnormalities and hyperplasia in mammary glands of mouse models. Oncogene 23, 8720–8730. doi: 10.1038/sj.onc.1208153
Zhang, J., Wu, Q., Wang, Z., Zhang, Y., Zhang, G., Fu, J., et al. (2015). Knockdown of PSF1 expression inhibits cell proliferation in lung cancer cells in vitro. Tumor Biol. 36, 2163–2168. doi: 10.1007/s13277-014-2826-8
Zhang, J. Y., Zhu, W., Imai, H., Kiyosawa, K., Chan, E. K. L., and Tan, E. M. (2001). De-novo humoral immune responses to cancer-associated autoantigens during transition from chronic liver disease to hepatocellular carcinoma. Clin. Exp. Immunol. 125, 3–9. doi: 10.1046/j.1365-2249.2001.01585.x
Keywords: separation of sister chromatids, HCC, biomarker, SAM, GSEA
Citation: Sun B, Lin G, Ji D, Li S, Chi G and Jin X (2018) Dysfunction of Sister Chromatids Separation Promotes Progression of Hepatocellular Carcinoma According to Analysis of Gene Expression Profiling. Front. Physiol. 9:1019. doi: 10.3389/fphys.2018.01019
Received: 23 February 2018; Accepted: 10 July 2018;
Published: 27 July 2018.
Edited by:
Cesare Gargioli, Università degli Studi di Roma Tor Vergata, ItalyReviewed by:
Gabriele Grassi, University of Trieste, ItalyFabiola Moretti, National Council of Italy, Italy
Copyright © 2018 Sun, Lin, Ji, Li, Chi and Jin. 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: Xingyi Jin, jin5756@163.com