Abstract
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
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
FIGURE 3
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
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
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
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
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
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.
Statements
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
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 S1Samples detail.
TABLE S1qPCR primers used for tissues specific expression analysis.
Footnotes
References
1
Aguirre-GamboaR.Gomez-RuedaH.Martinez-LedesmaE.Martinez-TorteyaA.Chacolla-HuaringaR.Rodriguez-BarrientosA.et al (2013). SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis.PLoS One8:e74250. 10.1371/journal.pone.0074250
2
AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene ontology: tool for the unification of biology. the gene ontology consortium.Nat. Genet.2525–29. 10.1038/75556
3
AsteritiI. A.RensenW. M.LindonC.LaviaP.GuarguagliniG. (2010). The Aurora-A/TPX2 complex: a novel oncogenic holoenzyme?Biochim. Biophys. Acta Rev. Cancer1806230–239. 10.1016/j.bbcan.2010.08.001
4
AuS. L.WongC. C.LeeJ. M.FanD. N.TsangF. H.NgI. O.et al (2012). Enhancer of zeste homolog 2 epigenetically silences multiple tumor suppressor microRNAs to promote liver cancer metastasis.Hepatology56622–631. 10.1002/hep.25679
5
BernalJ. A.LunaR.EspinaA.LazaroI.Ramos-MoralesF.RomeroF.et al (2002). Human securin interacts with p53 and modulates p53-mediated transcriptional activity and apoptosis.Nat. Genet.32306–311. 10.1038/ng997
6
BeroukhimR.MermelC. H.PorterD.WeiG.RaychaudhuriS.DonovanJ.et al (2010). The landscape of somatic copy-number alteration across human cancers.Nature463899–905. 10.1038/nature08822
7
BinO. Y.LanZ. D.MeadowsJ.PanH. Q.FukasawaK.LiW. Q.et al (1998). Human Bub1: a putative spindle checkpoint kinase closely linked to cell proliferation.Cell Growth Diff.9877–885.
8
BlumH. E. (2005). Treatment of hepatocellular carcinoma.Best Pract. Res. Clin. Gastroenterol.19129–145. 10.1016/j.bpg.2004.11.008
9
BoigeV.LaurentpuigP.FouchetP.FlejouJ. F.MongesG.BedossaP.et al (1997). Concerted nonsyntenic allelic losses in hyperploid hepatocellular carcinoma as determined by a high-resolution allelotype.Cancer Res.571986–1990.
10
CahillD. P.Da CostaL. T.Carson-WalterE. B.KinzlerK. W.VogelsteinB.LengauerC. (1999). Characterization of MAD2B and other mitotic spindle checkpoint genes.Genomics58181–187. 10.1006/geno.1999.5831
11
CahillD. P.LengauerC.YuJ.RigginsG. J.WillsonJ. K. V.MarkowitzS. D.et al (1998). Mutations of mitotic checkpoint genes in human cancers.Nature392300–303. 10.1038/32688
12
CarterS. L.EklundA. C.KohaneI. S.HarrisL. N.SzallasiZ. (2006). A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers.Nat. Genet.381043–1048.10.1038/ng1861
13
ChaC.DematteoR. R. (2005). Molecular mechanisms in hepatocellular carcinoma development.Best Pract. Res. Clin. Gastroenterol.1925–37. 10.1016/j.bpg.2004.11.005
14
ChaseA.CrossN. C. P. (2011). Aberrations of EZH2 in Cancer.Clin. Cancer Res.172613–2618. 10.1158/1078-0432.Ccr-10-2156
15
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):S11. 10.1186/1752-0509-8-S4-S11
16
DaiW. (2009). Suppression of genomic instabilities caused by chromosome Mis-segregation: a perspective from studying BubR1 and Sgo1.J. Formos. Med. Assoc.108904–911. 10.1016/s0929-6646(10)60002-2
17
DawlatyM. M.MalureanuL.JeganathanK. B.KaoE.SustmannC.TahkS.et al (2008). Resolution of sister centromeres requires RanBP2-mediated SUMOylation of topoisomerase II alpha.Cell133103–115. 10.1016/j.cell.2008.01.045
18
DuesbergP.LiR. (2003). Multistep carcinogenesis A chain reaction of aneuploidizations.Cell Cycle2202–210. 10.4161/cc.2.3.382
19
El-SeragH. B. (2002). Hepatocellular carcinoma: an epidemiologic view.J. Clin. Gastroenterol.35S72–S78. 10.1097/00004836-200211002-00002
20
FabregatA.JupeS.MatthewsL.SidiropoulosK.GillespieM.GarapatiP.et al (2018). The reactome pathway knowledgebase.Nucleic Acids Res.46D649–D655. 10.1093/nar/gkx1132
21
FennerM.WehrhanF.JehleM.AmannK.Radespiel-TrogerM.GrabenbauerG.et al (2005). Restricted-expressed proliferation-associated protein (repp86) expression in squamous cell carcinoma of the oral cavity.Strahlentherapie Und Onkologie181755–761. 10.1007/s00066-005-1430-7
22
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.Oncotargets Ther.92211–2220. 10.2147/ott.S101994
23
GraceC.NachevaE. 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.11173–183. 10.4137/CIN.S9258
24
GuoC.WuG.ChinJ. L.BaumanG.MoussaM.WangF.et al (2006). Bub1 up-regulation and hyperphosphorylation promote malignant transformation in SV40 tag-induced transgenic mouse models.Mol. Cancer Res.4957–969. 10.1158/1541-7786.MCR-06-0168
25
GuptaA.InabaS.WongO. K.FangG. W.LiuJ. W. (2003). Breast cancer-specific gene 1 interacts with the mitotic checkpoint kinase BubR1.Oncogene227593–7599. 10.1038/sj.onc.1206880
26
HanksS.ColemanK.ReidS.PlajaA.FirthH.FitzpatrickD.et al (2004). Constitutional aneuploidy and cancer predisposition caused by biallelic mutations in BUB1B.Nat. Genet.361159–1161. 10.1038/ng1449
27
HoD. W. H.KaiA. K. L.NgI. O. L. (2015). TCGA whole-transcriptome sequencing data reveals significantly dysregulated genes and signaling pathways in hepatocellular carcinoma.Front. Med.9322–330. 10.1007/s11684-015-0408-9
28
HoD. W. Y.YangZ. F.YiK.LamC. T.NgM. N. P.YuW. C.et al (2012). Gene expression profiling of liver cancer stem cells by RNA-sequencing.PLoS One7:e37159. 10.1371/journal.pone.0037159
29
HodgsonJ. G.ChinK.CollinsC.GrayJ. W. (2003). Genome amplification of chromosome 20 in breast cancer.Breast Cancer Res. Treat.78337–345. 10.1023/A:1023085825042
30
HokkaD.ManiwaY.TaneS.NishioW.YoshimuraM.OkitaY.et al (2013). Psf3 is a prognostic biomarker in lung adenocarcinoma.Lung Cancer7977–82. 10.1016/j.lungcan.2012.10.001
31
HuangD. W.ShermanB. T.LempickiR. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat. Protoc.444–57. 10.1038/nprot.2008.211
32
JiaH. L.YeQ. H.QinL. X.BudhuA.ForguesM.ChenY.et al (2007). Gene expression profiling reveals potential biomarkers of human hepatocellular carcinoma.Clin. Cancer Res.131133–1139. 10.1158/1078-0432.CCR-06-1025
33
JonesP. A.BaylinS. B. (2002). The fundamental role of epigenetic events in cancer.Nat. Rev. Genet.3415–428. 10.1038/nrg816
34
KimH.-E.KimD.-G.LeeK. J.SonJ. G.SongM.-Y.ParkY.-M.et al (2012). Frequent amplification of CENPF, GMNN and CDK13 genes in hepatocellular carcinomas.PLoS One7:e43223. 10.1371/journal.pone.0043223
35
KimW.BirdG. H.NeffT.GuoG. J.KerenyiM. A.WalenskyL. D.et al (2013). Targeted disruption of the EZH2-EED complex inhibits EZH2-dependent cancer.Nat. Chem. Biol.9643–650. 10.1038/Nchembio.1331
36
KnutsonS. K.WarholicN. M.WigleT. J.KlausC. R.AllainC. J.RaimondiA.et al (2013). Durable tumor regression in genetically altered malignant rhabdoid tumors by inhibition of methyltransferase EZH2.Proc. Natl. Acad. Sci. U.S.A.1107922–7927. 10.1073/pnas.1303800110
37
KnutsonS. K.WigleT. J.WarholicN. M.SneeringerC. J.AllainC. J.KlausC. R.et al (2012). A selective inhibitor of EZH2 blocks H3K27 methylation and kills mutant lymphoma cells.Nat. Chem. Biol.8890–896. 10.1038/Nchembio.1084
38
KnuutilaS.BjorkqvistA. M.AutioK.TarkkanenM.WolfM.MonniO.et al (1998). DNA copy number amplifications in human neoplasms - Review of comparative genomic hybridization studies.Am. J. Pathol.1521107–1123.
39
KronenwettU.HuwendiekS.OstringC.PortwoodN.RoblickU. J.PawitanY.et al (2004). Improved grading of breast adenocarcinomas based on genomic instability.Cancer Res.64904–909. 10.1158/0008-5472.can-03-2451
40
LabibK.GambusA. (2007). A key role for the GINS complex at DNA replication forks.Trends Cell Biol.17271–278. 10.1016/j.tcb.2007.04.002
41
LiB.QiX. Q.ChenX.HuangX.LiuG. Y.ChenH. R.et al (2010). Expression of targeting protein for Xenopus kinesin-like protein 2 is associated with progression of human malignant astrocytoma.Brain Res.1352200–207. 10.1016/j.brainres.2010.06.060
42
LiW. Q.LanZ. D.WuH. Y.WuS. C.MeadowsJ.ChenJ.et al (1999). BUBR1 phosphorylation is regulated during mitotic checkpoint activation.Cell Growth Diff.10769–775.
43
LinK. T.ShannY. J.ChauG. Y.HsuC. N.HuangC. Y. F. (2014). Identification of latent biomarkers in hepatocellular carcinoma by ultra-deep whole-transcriptome sequencing.Oncogene334786–4794. 10.1038/onc.2013.424
44
LiuQ.TuK.ZhangH.ZhengX.YaoY.LiuQ. (2015). TPX2 as a novel prognostic biomarker for hepatocellular carcinoma.Hepatol. Res.45906–918. 10.1111/hepr.12428
45
LiuQ.YangP.TuK.ZhangH.ZhengX.YaoY.et al (2014). TPX2 knockdown suppressed hepatocellular carcinoma cell invasion via inactivating AKT signaling and inhibiting MMP2 and MMP9 expression.Chin. J. Cancer Res.26410–417. 10.3978/j.issn.1000-9604.2014.08.01
46
MaL.ZhaoX. S.ZhuX. L. (2006). Mitosin/CENP-F in mitosis, transcriptional control, and differentiation.J. Biomed. Sci.13205–213. 10.1007/s11373-005-9057-3
47
MalureanuL. A.JeganathanK. B.HamadaM.WasilewskiL.DavenportJ.Van DeursenJ. M. (2009). BubR1 N terminus acts as a soluble inhibitor of cyclin B Degradation by APC/C-Cdc20 in interphase.Dev. Cell16118–131. 10.1016/j.devcel.2008.11.004
48
MargueronR.ReinbergD. (2011). The Polycomb complex PRC2 and its mark in life.Nature469343–349. 10.1038/nature09784
49
McCabeM. T.OttH. M.GanjiG.KorenchukS.ThompsonC.Van AllerG. S.et al (2012). EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations.Nature492108–112. 10.1038/nature11606
50
MittalS.El-SeragH. B. (2013). Epidemiology of hepatocellular carcinoma: consider the population.J. Clin. Gastroenterol.47(Suppl.), S2–S6. 10.1097/MCG.0b013e3182872f29
51
MomparlerR. L.IdaghdourY.MarquezV. E.MomparlerL. F. (2012). Synergistic antileukemic action of a combination of inhibitors of DNA methylation and histone methylation.Leuk. Res.361049–1054. 10.1016/j.leukres.2012.03.001
52
MontojoJ.ZuberiK.RodriguezH.KaziF.WrightG.DonaldsonS. L.et al (2010). GeneMANIA Cytoscape plugin: fast gene function predictions on the desktop.Bioinformatics262927–2928. 10.1093/bioinformatics/btq562
53
MyrieK. A.PercyM. J.AzimJ. N.NeeleyC. K.PettyE. M. (2000). Mutation and expression analysis of human BUB1 and BUB1B in aneuploid breast cancer cell lines.Cancer Lett.152193–199. 10.1016/S0304-3835(00)00340-2
54
NagahamaY.UenoM.MiyamotoS.MoriiE.MinamiT.MochizukiN.et al (2010). PSF1, a DNA replication factor expressed widely in stem and progenitor cells, drives tumorigenic and metastatic properties.Cancer Res.701215–1224. 10.1158/0008-5472.Can-09-3662
55
NagaiH.PineauP.TiollaisP.BuendiaM. A.DejeanA. (1997). Comprehensive allelotyping of human hepatocellular carcinoma.Oncogene142927–2933. 10.1038/sj.onc.1201136
56
NasmythK. (1999). Separating sister chromatids.Trends Biochem. Sci.2498–104. 10.1016/s0968-0004(99)01358-4
57
PanigrahiA. K.PatiD. (2009). Road to the crossroads of life and death: linking sister chromatid cohesion and separation to aneuploidy, apoptosis and cancer.Crit. Rev. Oncol. Hematol.72181–193. 10.1016/j.critrevonc.2008.12.002
58
QiW.ChanH. M.TengL.LiL.ChuaiS. N.ZhangR. P.et al (2012). Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation.Proc. Natl. Acad. Sci. U.S.A.10921360–21365. 10.1073/pnas.1210371110
59
RichN. E.ParikhN. D.SingalA. G. (2017). Overdiagnosis: an understudied issue in hepatocellular carcinoma surveillance.Semin. Liver Dis.37296–304. 10.1055/s-0037-1608775
60
SaekiA.TamuraS.ItoN.KisoS.MatsudaY.YabuuchiI.et al (2002). Frequent impairment of the spindle assembly checkpoint in hepatocellular carcinoma.Cancer942047–2054. 10.1002/cncr.10448
61
SchvartzmanJ.-M.SotilloR.BenezraR. (2010). Mitotic chromosomal instability and cancer: mouse modelling of the human disease.Nat. Rev. Cancer10102–115. 10.1038/nrc2781
62
ScottoL.NarayanG.NandulaS.Arias-PulidoH.SubramaniyamS.SchneiderA.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 Cancer47755–765. 10.1002/gcc.20577
63
SeikeM.GemmaA.HosoyaY.HosomiY.OkanoT.KurimotoF.et al (2002). The promoter region of the human BUBR1 gene and its expression analysis in lung cancer.Lung Cancer38229–234. 10.1016/S0169-5002(02)00218-0
64
ShichiriM.YoshinagaK.HisatomiH.SugiharaK.HirataY. (2002). Genetic and epigenetic inactivation of mitotic checkpoint genes hBUB1 and hBUBR1 and their relationship to survival.Cancer Res.6213–17.
65
SpinziG.PaggiS. (2008). Sorafenib in advanced hepatocellular carcinoma.N. Engl. J. Med.3592497–2498. 10.1056/NEJMc081780
66
SpositoC.MarianiL.GerminiA.Flores ReyesM.BonginiM.GrossiG.et al (2013). Comparative efficacy of sorafenib versus best supportive care in recurrent hepatocellular carcinoma after liver transplantation: a case-control study.J. Hepatol.5959–66. 10.1016/j.jhep.2013.02.026
67
SubramanianA.TamayoP.MoothaV. K.MukherjeedS.EbertB. L.GilletteM. 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.10215545–15550. 10.1073/pnas.0506580102
68
TuratoC.BalassoA.CarloniV.TiribelliC.MastrottoF.MazzoccaA.et al (2017). New molecular targets for functionalized nanosized drug delivery systems in personalized therapy for hepatocellular carcinoma.J. Control. Release268184–197. 10.1016/j.jconrel.2017.10.027
69
UmemuraT.IchijoT.YoshizawaK.TanakaE.KiyosawaK. (2009). Epidemiology of hepatocellular carcinoma in Japan.J. Gastroenterol.44(Suppl. 19), 102–107. 10.1007/s00535-008-2251-0
70
VaderG.LensS. M. A. (2008). The Aurora kinase family in cell division and cancer.Biochim. Biophys. Acta Rev. Cancer178660–72. 10.1016/j.bbcan.2008.07.003
71
VarisA.SalmelaA.-L.KallioM. J. (2006). Cenp-F (mitosin) is more than a mitotic marker.Chromosoma115288–295. 10.1007/s00412-005-0046-0
72
WalterW.Sanchez-CaboF.RicoteM. (2015). GOplot: an R package for visually combining expression data with functional analysis.Bioinformatics312912–2914. 10.1093/bioinformatics/btv300
73
WarnerS. L.StephensB. J.NwokenkwoS.HostetterG.SugengA.HidalgoM.et al (2009). Validation of TPX2 as a potential therapeutic target in pancreatic cancer cells.Clin. Cancer Res.156519–6528. 10.1158/1078-0432.Ccr-09-0077
74
WilsonC.MillerC. (2005). Simpleaffy: a BioConductor package for affymetrix quality control and data analysis.Bioinformatics213683–3685. 10.1093/bioinformatics/bti605
75
YamamotoY.MatsuyamaH.ChochiY.OkudaM.KawauchiS.InoueR.et al (2007). Overexpression of BUBR1 is associated with chromosomal instability in bladder cancer.Cancer Genet. Cytogenet.17442–47. 10.1016/j.cancergencyto.2006.11.012
76
YinL.ChangC.XuC. (2017). G2/M checkpoint plays a vital role at the early stage of HCC by analysis of key pathways and genes.Oncotarget876305–76317. 10.18632/oncotarget.19351
77
ZhangD. W.HirotaT.MarumotoT.ShimizuM.KunitokuN.SasayamaT.et al (2004). Cre-loxP-controlled periodic Aurora-A overexpression induces mitotic abnormalities and hyperplasia in mammary glands of mouse models.Oncogene238720–8730. 10.1038/sj.onc.1208153
78
ZhangJ.WuQ.WangZ.ZhangY.ZhangG.FuJ.et al (2015). Knockdown of PSF1 expression inhibits cell proliferation in lung cancer cells in vitro.Tumor Biol.362163–2168. 10.1007/s13277-014-2826-8
79
ZhangJ. Y.ZhuW.ImaiH.KiyosawaK.ChanE. K. L.TanE. M. (2001). De-novo humoral immune responses to cancer-associated autoantigens during transition from chronic liver disease to hepatocellular carcinoma.Clin. Exp. Immunol.1253–9. 10.1046/j.1365-2249.2001.01585.x
80
ZhouL.SunX.-J.LiuC.WuQ.-F.TaiM.-H.WeiJ.-C.et al (2015). Overexpression of PSF1 is correlated with poor prognosis in hepatocellular carcinoma patients.Int. J. Biol. Mark.30E56–E64. 10.5301/jbm.5000105
Summary
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
Volume
9 - 2018
Edited by
Cesare Gargioli, Università degli Studi di Roma Tor Vergata, Italy
Reviewed by
Gabriele Grassi, University of Trieste, Italy; Fabiola Moretti, National Council of Italy, Italy
Updates
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
This article was submitted to Integrative Physiology, a section of the journal Frontiers in Physiology
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.