Research Paper Volume 11, Issue 2 pp 423—447
Expression and gene regulation network of RBM8A in hepatocellular carcinoma based on data mining
- 1 Department of Medical Oncology, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, Guangxi 530021, People’s Republic of China
- 2 Maternal and Child Health Hospital and Obstetrics and Gynecology Hospital of Guangxi Zhuang Autonomous Region, Guangxi 530021, People’s Republic of China
- 3 Department of Medical Oncology, Affiliated Langdong Hospital of Guangxi Medical University, Nanning, Guangxi 530021, People’s Republic of China
- 4 The Fifth Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi 530021, People’s Republic of China
- 5 Department of Biology, Pennsylvania State University, University Park, PA 16802, USA
Received: October 12, 2018 Accepted: December 25, 2018 Published: January 22, 2019https://doi.org/10.18632/aging.101749
How to Cite
Copyright: Lin et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
RNA binding motif protein 8A (RBM8A) is an RNA binding protein in a core component of the exon junction complex. Abnormal RBM8A expression is associated with carcinogenesis. We used sequencing data from the Cancer Genome Atlas database and Gene Expression Omnibus, analyzed RBM8A expression and gene regulation networks in hepatocellular carcinoma (HCC). Expression was analyzed using OncomineTM and Gene Expression Profiling Interactive Analysis tools, while RBM8A alterations and related functional networks were identified using cBioPortal. LinkedOmics was used to identify differential gene expression with RBM8A and to analyze Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. Gene enrichment analysis examined target networks of kinases, miRNAs and transcription factors. We found that RBM8A is overexpressed and the RBM8A gene often amplified in HCC. Expression of this gene is linked to functional networks involving the ribosome and RNA metabolic signaling pathways. Functional network analysis suggested that RBM8A regulates the spliceosome, ribosome, DNA replication and cell cycle signaling via pathways involving several cancer-related kinases, miRNAs and E2F Transcription Factor 1. Our results demonstrate that data mining efficiently reveals information about RBM8A expression and potential regulatory networks in HCC, laying a foundation for further study of the role of RBM8A in carcinogenesis.
Hepatocellular carcinoma (HCC) is one of the top ten malignant tumors and the third leading cause of cancer-related death in the world . China accounts for 55% of new HCC cases and HCC-related deaths every year . Due to the high recurrence and metastasis, the 5-year survival rate of patients with advanced HCC does not exceed 5% . The development of various targeted drugs have prolonged the survival of patients and made a revolutionary breakthrough in the treatment of advanced HCC. However, existing targeted drugs show unsatisfactory efficacy. Even with sorafenib or regorafenib therapy, the overall life expectancy of HCC patients is less than 1 year [4, 5]. The pathogenesis of HCC is extremely complex, involving processes such as cell cycle regulation and signal transduction, which reflects the function and interaction of multiple genes at multiple steps . It may be possible to identify new drug targets for HCC by screening gene networks for changes related to tumor formation and progression.
RNA binding motif protein 8A (RBM8A), also known as Y14, is a conserved protein widely expressed in cells, and is an RNA binding motif protein . RBM8A serves as a core factor of the RNA surveillance machinery for the exon junction complex (EJC). It is also known to be the core protein of nonsense-mediated mRNA decay (NMD), which monitors abnormal mRNA in eukaryotes [8, 9]. RBM8A is related to RNA transcription, translation, cell cycle regulation and apoptosis regulation [10, 11], and it is also involved in several crucial cell signaling pathways [12–15], in which it plays an important role in tumorigenesis and development. Abnormal expression of RBM8A was first detected in cervical cancer , and its overexpression was subsequently detected in various malignant tumors, including non-small cell lung cancer and myeloma [16–18].
In our studies with patient samples from Guangxi, one of the regions with the highest incidence of HCC in China, we found that RBM8A was overexpressed in HCC tumor tissues compared to normal liver tissues, and this overexpression was associated with the surface antigen of the hepatitis B virus (HBsAg) and Edmondson pathological grading. Moreover, Kaplan-Meier survival analysis showed that high RBM8A expression was associated with poor overall survival and progression-free survival of HCC patients. Gain- and loss-of-function experiments further demonstrated that RBM8A promotes invasion and metastasis via the signaling pathway involved in the endothelial-to-mesenchymal (EMT) transition, while loss of RBM8A induces apoptosis .
These results suggest that RBM8A is a novel proto-oncogene. Thus, we studied RBM8A expression and mutations in data from patients with HCC in The Cancer Genome Atlas (TCGA) and various public databases. Using multi-dimensional analysis methods, we analyzed genomic alterations and functional networks related to RBM8A in HCC. Thus, our results could potentially reveal new targets and strategies for HCC diagnosis and treatment.
RBM8A expression in HCC
We initially evaluated RBM8A transcription levels in multiple HCC studies from TCGA and the Gene Expression Omnibus (GEO). Data in the Oncomine 4.5 database revealed that mRNA expression and DNA copy number variation (CNV) of RBM8A were significantly higher in HCC tissues than in normal tissues (p < 0.01). Although the fold differences were within 2, RBM8A ranked within the top 25% based on mRNA expression and within the top 5% based on DNA CNVs (Figure 1). Further sub-group analysis of multiple clinic pathological features of 371 liver hepatocellular carcinoma (LIHC) samples in the TCGA consistently showed high transcription of RBM8A. The transcription level of RBM8A was significantly higher in HCC patients than healthy people in subgroup analyses based on gender, age, ethnicity, disease stages and tumor grade (Figure 2). Thus, RBM8A expression may serve as a potential diagnostic indicator in HCC.
Figure 1. RBM8A transcription in hepatocellular carcinoma (Oncomine). Levels of RBM8A mRNA and DNA copy number were significantly higher in hepatocellular carcinoma than in normal tissue. Shown are fold change, associated p values, and overexpression gene rank, based on Oncomine 4.5 analysis. (A–D) Box plot showing RBM8A mRNA levels in, respectively, the Roessler Liver 2, Chen Liver, Wurmbach Liver and Mas Liver datasets. (E–F) Box plot showing RBM8A copy number in The Cancer Genome Atlas (TCGA) Liver and Guichard Liver datasets, respectively.
Figure 2. RBM8A transcription in subgroups of patients with hepatocellular carcinoma, stratified based on gender, age and other criteria (UALCAN). (A) Boxplot showing relative expression of RBM8A in normal and LIHC samples. (B) Boxplot showing relative expression of RBM8A in normal individuals of either gender or male or female LIHC patients. (C) Boxplot showing relative expression of RBM8A in normal individuals of any age or in LIHC patients aged 21–40, 41–60, 61–80, or 81–100 yr. (D) Boxplot showing relative expression of RBM8A in normal individuals of any ethnicity or in LIHC patients of Caucasian, African-American or Asian ethnicity. (E) Boxplot showing relative expression of RBM8A in normal individuals or in LIHC patients in stages 1, 2, 3 or 4. (F) Boxplot showing relative expression of RBM8A in normal individuals or LIHC patients with grade 1, 2, 3 or 4 tumors. Data are mean ± SE. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Genomic alterations of RBM8A in HCC
Frequency and type of RBM8A alterations in HCC
We then used the cBioPortal to determine the types and frequency of RBM8A alterations in HCC based on sequencing data from LIHC patients in the TCGA database. RBM8A was altered in 80 of 370 (22%) LIHC patients (Figure 3A). These alterations were mRNA upregulation in 60 cases (16%), amplification in 38 cases (10%), mutation in 1 case (0.3%), and multiple alterations in 19 cases (5%) (Table 1). Thus, amplification is the most common type of RBM8A CNV in HCC.
Figure 3. Visual summary of RBM8A alterations and biological interaction network in hepatocellular carcinoma (cBioPortal). (A) OncoPrint of RBM8A alterations in LIHC. The OncoPrint provides an overview of genomic alterations in RBM8A affecting individual samples (columns) in LIHC from the TCGA. The different types of genetic alterations are highlighted in different colors. (B) Network view of the RBM8A neighborhood in LIHC. RBM8A are seed genes (indicated with thick border), and all other genes are automatically identified as altered in LIHC. Darker red indicates increased frequency of alteration in LIHC. The interaction types are derived from the Biological Pathway Exchange (BioPAX): the blue connection indicates that the first protein controls a reaction that changes the state of the second protein; the red connection indicates that the proteins are members of the same complex.
Table 1. The type and frequency of RBM8A neighbor gene alterations in hepatocellular carcinoma (cBioPortal)
|Gene Symbol||Amplification||Homozygous Deletion||Up-regulation||Down-regulation||Mutation||Total Alteration|
Biological interaction network of RBM8A alterations in HCC
We next wanted to determine the biological interaction network of RBM8A in HCC. To do this, we used the tab Network in cBioPortal to show RBM8A-neighboring genes that were altered at frequencies >20% (Figure 3B and Table 1). The neighbor genes of RBM8A with the most frequent alterations were POLR2K (40.7%), RPL30 (38.8%) and CPSF1 (37.2%). The 50 most frequently altered neighbor genes were determined using Network. Analysis of significantly enriched gene ontology (GO) terms indicated that these genes encode proteins localized mainly to the cytosol, ribosome, and ribosome subunits. These proteins are primarily involved in viral gene expression and RNA catabolism, while they also serve as structural constituents of ribosomes and mRNA binding (Figure 4A–4C). Similarly, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed enrichment in ribosome signaling, RNA transport, mRNA surveillance and spliceosome signaling pathways (Figure 4D, 4E). Thus, the biological interaction network of RBM8A alterations is involved in the ribosome and several RNA metabolic processes.
Figure 4. Enrichment analysis of the genes altered in the RBM8A neighborhood in hepatocellular carcinoma. The bubble diagrams display the enrichment results of the top 50 genes altered in the RBM8A neighborhood in LIHC. (A) Cellular components. (B) Biological processes. (C) Molecular functions. (D) KEGG pathway analysis. (E) KEGG pathway annotations of the ribosome pathway. The red star denotes altered genes.
Enrichment analysis of RBM8A functional networks in HCC
GO and KEGG pathway analyses of co-expression genes correlated with RBM8A in HCC
The Function module of LinkedOmics was used to analyze mRNA sequencing data from 371 LIHC patients in the TCGA. As shown in the volcano plot (Figure 5A), 2,596 genes (dark red dots) showed significant positive correlations with RBM8A, whereas 3,050 genes (dark green dots) showed significant negative correlations (false discovery rate [FDR] < 0.01). The 50 significant gene sets positively and negatively correlated with RBM8A as shown in the heat map (Figure 5B, 5C). This result suggests a widespread impact of RBM8A on the transcriptome. The statistical scatter plots for individual genes are shown in Supplementary Figure 1A–1C. RBM8A expression showed a strong positive association with expression of POLR3C (positive rank #1, Pearson correlation = 0.68, p = 1.48e–52), VPS72 (Pearson correlation = 0.61, p = 6.64e–39), and MRPL9 (Pearson correlation = 0.60, p = 3.64e–37), which reflect changes in RNA polymerase III transcription initiation, transcriptional regulation/DNA repair/apoptosis and mitochondrial ribosome composition. Significant GO term analysis by gene set enrichment analysis (GSEA) showed that genes differentially expressed in correlation with RBM8A were located mainly in the condensed chromosome, replication fork and spliceosome complex, where they participate primarily in cell cycle checkpoint control, DNA replication and mRNA processing. They act as structural constituents in ribosomes and monooxygenase (Figure 6A–6C and Supplementary Tables 1–3). KEGG pathway analysis showed enrichment in the spliceosome, ribosome, DNA replication and cell cycle pathways (Figure 6D, 6E and Supplementary Table 4).
Figure 5. Genes differentially expressed in correlation with RBM8A in hepatocellular carcinoma (LinkedOmics). (A) A Pearson test was used to analyze correlations between RBM8A and genes differentially expressed in LIHC. (B–C) Heat maps showing genes positively and negatively correlated with RBM8A in LIHC (TOP 50). Red indicates positively correlated genes and green indicates negatively correlated genes.
Figure 6. Significantly enriched GO annotations and KEGG pathways of RBM8A in hepatocellular carcinoma. The significantly enriched GO annotations and KEGG pathways of RBM8A co-expression genes in LIHC were analyzed using GSEA. (A) Cellular components. (B) Biological processes. (C) Molecular functions. (D) KEGG pathway analysis. The blue column represents the LeadingEdgeNum, and the orange represents the false discovery rate (FDR). The FDR from GSEA in the figure is 0. (E) KEGG pathway annotations of the cell cycle pathway. Red marked nodes are associated with the LeadingEdgeGene.
RBM8A networks of kinase, miRNA or transcription factor targets in HCC
To further explore the targets of RBM8A in HCC, we analyzed the kinase, miRNA and transcription factor target networks of positively correlated gene sets generated by GSEA. The top 5 most significant target networks were the kinase-target networks related primarily to the kinases ataxia telangiectasia-mutated and Rad3-related (ATR), Aurora kinase B (AURKB), cyclin-dependent kinase 1 (CDK1), cyclin-dependent kinase 2 (CDK2) and checkpoint kinase 1 (CHEK1) (Table 2 and Supplementary Tables 5–7). The miRNA-target network was associated with (CTTTGCA) MIR-527, (CCCACAT) MIR-299-3P and (GTGCAAT) MIR-25, MIR-32, MIR-92, MIR-363, and MIR-367. The transcription factor-target network was related mainly to the E2F Transcription Factor (E2F) family, including E2F1DP1RB_01, E2F1DP1_01 and E2F1DP2_01. The protein-protein interaction network constructed by GeneMANIA revealed correlation among genes for the kinases ATR, miRNA-527 and TF E2F1DP1RB_01. The gene set enriched for kinase ATR and transcription factor E2F1DP1RB is responsible mainly for regulating DNA replication, DNA repair and the cell cycle checkpoint (Figure 7 and Supplementary Figure 3). In addition, the gene set enriched for transcription factor E2F1DP1RB is also involved in the G1/S transition in mitosis, telomere maintenance and the MCM complex. The gene set enriched for miRNA-527 is involved mainly in regulation of the tricarboxylic acid cycle, organ development and angiogenesis (Supplementary Figure 2).
Table 2. The Kinase, miRNA and transcription factor-target networks of RBM8A in in hepatocellular carcinoma (LinkedOmics)
|miRNA Target||CTTTGCA, MIR-527||71||0.025|
|GTGCAAT, MIR-25, MIR-32, MIR-92, MIR-363, MIR-367||78||0.029|
|Transcription Factor Target||V$E2F1DP1RB_01||76||0|
Figure 7. Protein-protein interaction network of ATR kinase-target networks (GeneMANIA). Protein-protein interaction (PPI) network and functional analysis indicating the gene set that was enriched in the target network of ATR kinases. Different colors of the network edge indicate the bioinformatics methods applied: co-expression, website prediction, pathway, physical interactions and co-localization. The different colors for the network nodes indicate the biological functions of the set of enrichment genes.
Differential expression and dysfunction of the exon junction complex (EJC) core protein has been reported in various cancers[7, 20–22]. RBM8A, a core component of this complex, is involved in multiple steps of transcription [8, 23]. In previous work, we found that RBM8A was overexpressed in 105 HCC tissue samples from patients in Guangxi, and that high RBM8A expression predicted poor prognosis . In that work, we also showed in vitro that RBM8A promotes tumor cell migration and invasion in HCC by activating the EMT signaling pathway. To gain more detailed insights into the potential functions of RBM8A in HCC and its regulatory network, we performed bioinformatics analysis of public sequencing data to guide future research in HCC.
Early detection of HCC is a difficult problem that perplexes clinicians. Alpha-fetoprotein (AFP) has long been used as an indicator for early screening of HCC, but only approximately 70% of HCC patients are AFP-positive . Therefore, new HCC markers are needed to improve early diagnosis. Analysis of transcriptional sequencing data from more than 1000 clinical samples in the GEO and TCGA databases comprising six geographic regions and ethnic HCC studies [25–29] confirmed that RBM8A mRNA levels and CNVs are significantly higher in HCC than in normal liver tissue. The fold change is similar across the various HCC studies and, while not large, ranks RBM8A among the top 3–4% of all genes upregulated in HCC, based on CNVs. We suggest that RBM8A overexpression occurs in many cases of HCC and deserves further clinical validation as a potential diagnostic and prognostic marker.
CNVs can have major genomic implications, disrupting genes and altering genetic content, leading to phenotypic differences . Our study found that the copy number of RBM8A was increased in HCC, and that the major type of RBM8A alteration was amplification, which was associated with shorter survival. We speculate that altered RBM8A expression and RBM8A dysfunction in HCC may result from alterations in chromosomal structure. Since RBM8A plays several important physiological functions, its alteration may cause changes in various downstream signaling pathways. Indeed, neighboring gene networks close to RBM8A generally show different degrees of amplification in HCC. Related functional networks are involved in ribosome signaling, RNA transport, mRNA surveillance and spliceosome signaling. Thus, the network of RBM8A alterations is involved in the core node of post-transcriptional regulation, which is closely related to RNA splicing and protein translation. This is consistent with the physiological function of RBM8A [23, 31].
Enrichment analysis of target gene sets using GSEA can help reveal important networks of target kinases, miRNAs and transcription factors. Our results suggest that the functional network of RBM8A participates primarily in the spliceosome, ribosome, DNA replication and cell cycle. Like the mutation network, the functional network of RBM8A transcription is involved in genomic stability, gene expression and the cell cycle. These findings are consistent with the fact that RBM8A is critical for efficient and faithful splicing of a specific subset of short introns in mitotic cell cycle-related genes . It is critical to understand how alteration in a protein important for ensuring normal transcription can lead to major dysfunction and even to cancer such as HCC.
Genomic instability and mutagenesis are fundamental characteristics of cancer cells, and kinases and their associated signaling pathways help stabilize and repair genomic DNA [33, 34]. We found that RBM8A in HCC is associated with a network of kinases including ATR, AURKB, and CDK1. These kinases regulate genomic stability, mitosis and the cell cycle [35, 36]. In fact, ATR is one of the core kinase regulators of genomic stability; it initiates cellular responses to genomic instability and repair, directly phosphorylating more than 1000 important substrates, including tumor suppressor gene p53 protein and cell cycle regulatory proteins . ATR kinase inhibitors can kill tumor cells and synergize with the cell-killing effects of chemoradiotherapy . In HCC, RBM8A may regulate DNA replication, repair, and cell cycle progression via ATR kinase.
In 2011, Hanahan et al. described the 10 hallmark features of tumors, with “continuous proliferation” at the top . The most important reason is that abnormal expression of cell cycle regulatory factors in tumor cells leads to cell cycle disorder, which results in aberrant proliferation, decreased differentiation, decreased apoptosis, and rapid multiplication and development. E2F1 is one of the key links in the cell cycle regulation network . Abnormal E2F1 expression is actively involved in the occurrence and development of HCC. Studies have shown that increased expression of E2F1 is associated with poor prognosis in HCC patients , and studies with a transgenic mouse model of E2F1-induced HCC showed that the network of transcription factors targeted by RBM8A is related to E2F1 . Moreover, we showed previously that RBM8A activates the EMT signaling pathway, which E2F1 alters in a variety of solid tumors, including HCC [43, 44]. Therefore, our analyses suggest that E2F1 is an important target of RBM8A, and that RBM8A acts through this factor to regulate the cell cycle and proliferation capacity of HCC. Further studies should test this hypothesis.
Our study identified several miRNAs that were associated with RBM8A. These short (20–24 nt) non-coding RNAs, normally involved in post-transcriptional regulation of gene expression, can contribute to human carcinogenesis . The particular miRNAs in our study have been linked to tumor proliferation, apoptosis, cell cycle, invasion, metastasis, drug resistance and angiogenesis [46–49]. In fact, miR-363, miR-25 and miR-299 can be used as diagnostic and prognostic markers of HCC [50–52]. While miR-363 regulates E2F Transcription Factor 3 to inhibit HCC migration and invasion , miR-367 and miR-32 participate in EMT progression [54, 55]. Dysregulation of these miRNAs is consistent with the phenotype of RBM8A overexpression in HCC from our previous work .
This study provides multi-level evidence for the importance of RBM8A in hepatocarcinogenesis and its potential as a marker in HCC. Our results suggest that RBM8A overexpression in HCC has far-reaching effects in genomic stability and at multiple steps of gene expression (DNA replication, RNA splicing and protein translation) and of the cell cycle. RBM8A is specifically related to several tumor-associated kinases (such as ATR), miRNAs (such as miRNA-363), and transcription factors (such as E2F1). This study uses online tools based on the most popular bioinformatics theories to perform target gene analyses on tumor data from public databases. Compared with traditional chip screening, this method has the advantages of large sample size, low cost, and simplicity. This enables large-scale HCC genomics research and subsequent functional studies.
At the same time, the TCGA database has limitations. One is that the TCGA LIHC samples contain three ethnic groups. The genetic background and etiology of HCC can differ significantly across ethnic groups. Another limitation is that the LIHC samples contain relatively few patients in stage 4, yet the clinical reality is that most HCC patients are first diagnosed when their disease is advanced and prognosis is extremely poor. Therefore, our results should be verified in clinical samples showing sufficient coverage of different ethnic groups and HCC stages. The third limitation is that transcriptome sequencing can detect only static mutations; it cannot directly provide information on protein activity or expression level. These questions should be addressed in follow-up studies using molecular biology techniques.
Materials and Methods
The mRNA expression and DNA copy number of RBM8A in HCC were analyzed within the Oncomine 4.5 database. Oncomine (www.oncomine.org), currently the world’s largest oncogene chip database and integrated data mining platform, contains 715 gene expression data sets and data from 86,733 cancer tissues and normal tissues . This analysis drew on a series of HCC studies, including Roessler Liver 2, Chen Liver, Wurmbach Liver, Mas Liver, TCGA Liver and Guichard Liver studies [25–29]. RBM8A expression was assessed in HCC tissue relative to its expression in normal tissue, and differences associated with p < 0.01 were considered significant.
UALCAN uses TCGA level 3 RNA-seq and clinical data from 31 cancer types, which is an interactive web-portal to perform in-depth analyses of TCGA gene expression data . One of the portal’s user-friendly features is that it allows analysis of relative expression of a query gene(s) across tumor and normal samples, as well as in various tumor sub-groups based on individual cancer stages, tumor grade or other clinicopathological features. UALCAN is publicly available at http://ualcan.path.uab.edu.
The cBio Cancer Genomics Portal (http://cbioportal.org) is an open-access resource for interactive exploration of multidimensional cancer genomics data sets, currently containing 225 cancer studies . We used c-BioPortal to analyze RBM8A alterations in the TCGA LIHC sample. The search parameters included mutation, CNVs, and mRNA expression. The tab OncoPrint displays an overview of genetic alterations per sample in RBM8A. The tab Network visualizes the biological interaction network of RBM8A derived from public pathway databases, with color-coding and filter options based on the frequency of genomic alterations in each gene. Neighboring genes with alteration frequencies greater than 20% were included. We then performed GO and KEGG pathway enrichment analyses of the most frequently altered neighbor genes using the “clusterProfiler” package in R . The GO annotation had three parts: cellular component (CC), biological process (BP), and molecular function (MF).
The LinkedOmics database (http://www.linkedomics.org/login.php) is a Web-based platform for analyzing 32 TCGA cancer-associated multi-dimensional datasets . The LinkFinder module of LinkedOmics was used to study genes differentially expressed in correlation with RBM8A in the TCGA LIHC cohort (n=371). Results were analyzed statistically using Pearson’s correlation coefficient. The LinkFinder also created statistical plots for individual genes. All results were graphically presented in volcano plots, heat maps or scatter plots. The LinkInterpreter module of LinkedOmics performs pathway and network analyses of differentially expressed genes. The comprehensive functional category database in the Web-based Gene SeT AnaLysis Toolkit (WebGestalt)  was applied. Data from the LinkFinder results were signed and ranked, and GSEA was used to perform analyses of GO (CC, BP and MF), KEGG pathways, kinase-target enrichment, miRNA-target enrichment and transcription factor-target enrichment. The latter two network analyses were based on the Molecular Signatures Database (MSigDB) . The rank criterion was an FDR < 0.05, and 500 simulations were performed.
GeneMANIA (http://www.genemania.org) is a flexible, user-friendly web interface for constructing protein-protein interaction (PPI) network, generating hypotheses about gene function, analyzing gene lists and prioritizing genes for functional assays . The website can set the source of the edge of the network, and it features several bioinformatics methods: physical interaction, gene co-expression, gene co-location, gene enrichment analysis and website prediction. We used GeneMANIA to visualize the gene networks and predict function of genes that GSEA identified as being enriched in HCC: kinase ATR, mi-RNA 527 and transcription factor E2F1.
Y.L, R.L. and Y.Q. designed and performed the research, analyzed and interpreted data, and drafted the manuscript. Y.L., J. Z., Y. L., C.Y., and G.Q. participated in data analysis and figure preparation. Z.L. and Y. L. reviewed the manuscript. D.Z and Y.M. conceived the study and participated in research design and data interpretation. All authors read and approved the final manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This study was supported by the National Natural Science Foundation of China (81660498), Guangxi Natural Science Foundation (2018GXNSFBA281091, 2016GXNSFBA380090 and 2015GXNSFCA139014), the Basic Ability Enhancement Program for Young and Middle-aged Teachers of Guangxi (2017KY0120), the Guangxi Medical and Health Appropriate Technology Development and Application Project (S2018 062 and S2017 101), the Project of Nanning Scientific Research and Technology Development Plan (20183040-2), the China Scholarship Council (201608455001 and 201708455059), the American Heart Association (13SDG14670072), and the Guangxi High-Level Medical Talents Training Project.
- 1. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012.Int J Cancer.2015;136:E359-86. https://doi.org/10.1002/ijc.29210 [PubMed]
- 2. Chen W. Cancer statistics: updated cancer burden in China.Chin J Cancer Res.2015;27:1 https://doi.org/10.3978/j.issn.1000-9604.2015.02.07 [PubMed]
- 3. Louafi S, Boige V, Ducreux M, Bonyhay L, Mansourbakht T, de Baere T, Asnacios A, Hannoun L, Poynard T, Taieb J. Gemcitabine plus oxaliplatin (GEMOX) in patients with advanced hepatocellular carcinoma (HCC): results of a phase II study.Cancer.2007;109:1384-90. https://doi.org/10.1002/cncr.22532 [PubMed]
- 4. Bruix J, Qin S, Merle P, Granito A, Huang YH, Bodoky G, Pracht M, Yokosuka O, Rosmorduc O, Breder V, Gerolami R, Masi G, Ross PJ, et al. Regorafenib for patients with hepatocellular carcinoma who progressed on sorafenib treatment (RESORCE): a randomised, double-blind, placebo-controlled, phase 3 trial.Lancet.2017;389:56-66. https://doi.org/10.1016/S0140-6736(16)32453-9 [PubMed]
- 5. Rimassa L and Santoro A. Sorafenib therapy in advanced hepatocellular carcinoma: the SHARP trial.Expert Rev Anticancer Ther.2009;9:739-45. https://doi.org/10.1586/era.09.41 [PubMed]
- 6. Levrero M and Zucman-Rossi J. Mechanisms of HBV-induced hepatocellular carcinoma.J Hepatol.2016;64:S84-S101. https://doi.org/10.1016/j.jhep.2016.02.021 [PubMed]
- 7. Salicioni AM, Xi M, Vanderveer LA, Balsara B, Testa JR, Dunbrack RL Jr, Godwin AK. Identification and structural analysis of human RBM8A and RBM8B: two highly conserved RNA-binding motif proteins that interact with OVCA1, a candidate tumor suppressor.Genomics.2000;69:54-62. https://doi.org/10.1006/geno.2000.6315 [PubMed]
- 8. Chuang TW, Lee KM, Tarn WY. Function and pathological implications of exon junction complex factor Y14.Biomolecules.2015;5:343-55. https://doi.org/10.3390/biom5020343 [PubMed]
- 9. Fatscher T, Boehm V, Gehring NH. Mechanism, factors, and physiological role of nonsense-mediated mRNA decay.Cell Mol Life Sci.2015;72:4523-44. https://doi.org/10.1007/s00018-015-2017-9 [PubMed]
- 10. Ishigaki Y, Nakamura Y, Tatsuno T, Hashimoto M, Shimasaki T, Iwabuchi K, Tomosugi N. Depletion of RNA-binding protein RBM8A (Y14) causes cell cycle deficiency and apoptosis in human cells.Exp Biol Med (Maywood).2013;238:889-97. https://doi.org/10.1177/1535370213494646 [PubMed]
- 11. Michelle L, Cloutier A, Toutant J, Shkreta L, Thibault P, Durand M, Garneau D, Gendron D, Lapointe E, Couture S, Le Hir H, Klinck R, Elela SA, et al. Proteins associated with the exon junction complex also control the alternative splicing of apoptotic regulators.Mol Cell Biol.2012;32:954-67. https://doi.org/10.1128/MCB.06130-11 [PubMed]
- 12. Lu CC, Lee CC, Tseng CT, Tarn WY. Y14 governs p53 expression and modulates DNA damage sensitivity.Sci Rep.2017;7:45558 https://doi.org/10.1038/srep45558 [PubMed]
- 13. Muromoto R, Taira N, Ikeda O, Shiga K, Kamitani S, Togi S, Kawakami S, Sekine Y, Nanbo A, Oritani K, Matsuda T. The exon-junction complex proteins, Y14 and MAGOH regulate STAT3 activation.Biochem Biophys Res Commun.2009;382:63-8. https://doi.org/10.1016/j.bbrc.2009.02.127 [PubMed]
- 14. Togi S, Shiga K, Muromoto R, Kato M, Souma Y, Sekine Y, Kon S, Oritani K, Matsuda T. Y14 positively regulates TNF-alpha-induced NF-kappaB transcriptional activity via interacting RIP1 and TRADD beyond an exon junction complex protein.J Immunol.2013;191:1436-44. https://doi.org/10.4049/jimmunol.1300501 [PubMed]
- 15. Wang C, Cigliano A, Delogu S, Armbruster J, Dombrowski F, Evert M, Chen X, Calvisi DF. Functional crosstalk between AKT/mTOR and Ras/MAPK pathways in hepatocarcinogenesis: implications for the treatment of human liver cancer.Cell Cycle.2013;12:1999-2010. https://doi.org/10.4161/cc.25099 [PubMed]
- 16. Carrasco DR, Tonon G, Huang Y, Zhang Y, Sinha R, Feng B, Stewart JP, Zhan F, Khatry D, Protopopova M, Protopopov A, Sukhdeo K, Hanamura I, et al. High-resolution genomic profiles define distinct clinico-pathogenetic subgroups of multiple myeloma patients.Cancer Cell.2006;9:313-25. https://doi.org/10.1016/j.ccr.2006.03.019 [PubMed]
- 17. Kim TJ, Choi JJ, Kim WY, Choi CH, Lee JW, Bae DS, Son DS, Kim J, Park BK, Ahn G, Cho EY, Kim BG. Gene expression profiling for the prediction of lymph node metastasis in patients with cervical cancer.Cancer Sci.2008;99:31-8. https://doi.org/10.1111/j.1349-7006.2007.00652.x [PubMed]
- 18. Petroziello J, Yamane A, Westendorf L, Thompson M, McDonagh C, Cerveny C, Law CL, Wahl A, Carter P. Suppression subtractive hybridization and expression profiling identifies a unique set of genes overexpressed in non-small-cell lung cancer.Oncogene.2004;23:7734-45. https://doi.org/10.1038/sj.onc.1207921 [PubMed]
- 19. Liang R, Lin Y, Ye JZ, Yan XX, Liu ZH, Li YQ, Luo XL, Ye HH. High expression of RBM8A predicts poor patient prognosis and promotes tumor progression in hepatocellular carcinoma.Oncol Rep.2017;37:2167-76. https://doi.org/10.3892/or.2017.5457 [PubMed]
- 20. Dvinge H, Kim E, Abdel-Wahab O, Bradley RK. RNA splicing factors as oncoproteins and tumour suppressors.Nat Rev Cancer.2016;16:413-30. https://doi.org/10.1038/nrc.2016.51 [PubMed]
- 21. Han D, Gao X, Wang M, Qiao Y, Xu Y, Yang J, Dong N, He J, Sun Q, Lv G, Xu C, Tao J, Ma N. Long noncoding RNA H19 indicates a poor prognosis of colorectal cancer and promotes tumor growth by recruiting and binding to eIF4A3.Oncotarget.2016;7:22159-73. https://doi.org/10.18632/oncotarget.8063 [PubMed]
- 22. Xu L, Dai W, Li J, He L, Wang F, Xia Y, Chen K, Li S, Liu T, Lu J, Zhou Y, Wang Y, Guo C. Methylation-regulated miR-124-1 suppresses tumorigenesis in hepatocellular carcinoma by targeting CASC3.Oncotarget.2016;7:26027-41. https://doi.org/10.18632/oncotarget.8266 [PubMed]
- 23. Chuang TW, Lee KM, Lou YC, Lu CC, Tarn WY. A Point Mutation in the Exon Junction Complex Factor Y14 Disrupts Its Function in mRNA Cap Binding and Translation Enhancement.J Biol Chem.2016;291:8565-74. https://doi.org/10.1074/jbc.M115.704544 [PubMed]
- 24. Luo P, Yin P, Hua R, Tan Y, Li Z, Qiu G, Yin Z, Xie X, Wang X, Chen W, Zhou L, Wang X, Li Y, et al. A Large-scale, multicenter serum metabolite biomarker identification study for the early detection of hepatocellular carcinoma.Hepatology.2018;67:662-75. https://doi.org/10.1002/hep.29561 [PubMed]
- 25. Chen X, Cheung ST, So S, Fan ST, Barry C, Higgins J, Lai KM, Ji J, Dudoit S, Ng IO, Van De Rijn M, Botstein D, Brown PO. Gene expression patterns in human liver cancers.Mol Biol Cell.2002;13:1929-39. https://doi.org/10.1091/mbc.02-02-0023 [PubMed]
- 26. Guichard C, Amaddeo G, Imbeaud S, Ladeiro Y, Pelletier L, Maad IB, Calderaro J, Bioulac-Sage P, Letexier M, Degos F, Clement B, Balabaud C, Chevet E, et al. Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma.Nat Genet.2012;44:694-8. https://doi.org/10.1038/ng.2256 [PubMed]
- 27. Mas VR, Maluf DG, Archer KJ, Yanek K, Kong X, Kulik L, Freise CE, Olthoff KM, Ghobrial RM, McIver P, Fisher R. Genes involved in viral carcinogenesis and tumor initiation in hepatitis C virus-induced hepatocellular carcinoma.Mol Med.2009;15:85-94. https://doi.org/10.2119/molmed.2008.00110 [PubMed]
- 28. Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, Thorgeirsson SS, Sun Z, Tang ZY, Qin LX, Wang XW. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients.Cancer Res.2010;70:10202-12. https://doi.org/10.1158/0008-5472.CAN-10-2607 [PubMed]
- 29. Wurmbach E, Chen YB, Khitrov G, Zhang W, Roayaie S, Schwartz M, Fiel I, Thung S, Mazzaferro V, Bruix J, Bottinger E, Friedman S, Waxman S, et al. Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma.Hepatology.2007;45:938-47. https://doi.org/10.1002/hep.21622 [PubMed]
- 30. Urrutia E, Chen H, Zhou Z, Zhang NR, Jiang Y. Integrative pipeline for profiling DNA copy number and inferring tumor phylogeny.Bioinformatics.2018;34:2126-8. https://doi.org/10.1093/bioinformatics/bty057 [PubMed]
- 31. Le Hir H, Sauliere J, Wang Z. The exon junction complex as a node of post-transcriptional networks.Nat Rev Mol Cell Biol.2016;17:41-54. https://doi.org/10.1038/nrm.2015.7 [PubMed]
- 32. Fukumura K, Wakabayashi S, Kataoka N, Sakamoto H, Suzuki Y, Nakai K, Mayeda A, Inoue K. The Exon Junction Complex Controls the Efficient and Faithful Splicing of a Subset of Transcripts Involved in Mitotic Cell-Cycle Progression.Int J Mol Sci.2016;17 https://doi.org/10.3390/ijms17081153 [PubMed]
- 33. Karimian A, Ahmadi Y, Yousefi B. Multiple functions of p21 in cell cycle, apoptosis and transcriptional regulation after DNA damage.DNA Repair (Amst).2016;42:63-71. https://doi.org/10.1016/j.dnarep.2016.04.008 [PubMed]
- 34. Yogosawa S and Yoshida K. Tumor suppressive role for kinases phosphorylating p53 in the DNA damage-induced apoptosis.Cancer Sci.2018; https://doi.org/10.1111/cas.13792 [PubMed]
- 35. Fang H, Niu K, Mo D, Zhu Y, Tan Q, Wei D, Li Y, Chen Z, Yang S, Balajee AS, Zhao Y. RecQL4-Aurora B kinase axis is essential for cellular proliferation, cell cycle progression, and mitotic integrity.Oncogenesis.2018;7:68 https://doi.org/10.1038/s41389-018-0080-4 [PubMed]
- 36. Herrera MC, Chymkowitch P, Robertson JM, Eriksson J, Boe SO, Alseth I, Enserink JM. Cdk1 gates cell cycle-dependent tRNA synthesis by regulating RNA polymerase III activity.Nucleic Acids Res.2018; https://doi.org/10.1093/nar/gky846 [PubMed]
- 37. Cortez D, Guntuku S, Qin J, Elledge SJ. ATR and ATRIP: partners in checkpoint signaling.Science.2001;294:1713-6. https://doi.org/10.1126/science.1065521 [PubMed]
- 38. Hafsi H, Dillon MT, Barker HE, Kyula JN, Schick U, Paget JT, Smith HG, Pedersen M, McLaughlin M, Harrington KJ. Combined ATR and DNA-PK Inhibition Radiosensitizes Tumor Cells Independently of Their p53 Status.Front Oncol.2018;8:245 https://doi.org/10.3389/fonc.2018.00245 [PubMed]
- 39. Hanahan D and Weinberg RA. Hallmarks of cancer: the next generation.Cell.2011;144:646-74. https://doi.org/10.1016/j.cell.2011.02.013 [PubMed]
- 40. Polager S and Ginsberg D. E2F - at the crossroads of life and death.Trends Cell Biol.2008;18:528-35. https://doi.org/10.1016/j.tcb.2008.08.003 [PubMed]
- 41. Chen YL, Uen YH, Li CF, Horng KC, Chen LR, Wu WR, Tseng HY, Huang HY, Wu LC, Shiue YL. The E2F transcription factor 1 transactives stathmin 1 in hepatocellular carcinoma.Ann Surg Oncol.2013;20:4041-54. https://doi.org/10.1245/s10434-012-2519-8 [PubMed]
- 42. Conner EA, Lemmer ER, Omori M, Wirth PJ, Factor VM, Thorgeirsson SS. Dual functions of E2F-1 in a transgenic mouse model of liver carcinogenesis.Oncogene.2000;19:5054-62. https://doi.org/10.1038/sj.onc.1203885 [PubMed]
- 43. Qi P, Lin WR, Zhang M, Huang D, Ni SJ, Zhu XL, Bai QM, Sheng WQ, Du X, Zhou XY. E2F1 induces LSINCT5 transcriptional activity and promotes gastric cancer progression by affecting the epithelial-mesenchymal transition.Cancer Manag Res.2018;10:2563-71. https://doi.org/10.2147/CMAR.S171652 [PubMed]
- 44. Chen J, Gong C, Mao H, Li Z, Fang Z, Chen Q, Lin M, Jiang X, Hu Y, Wang W, Zhang X, Chen X, Li H. E2F1/SP3/STAT6 axis is required for IL-4-induced epithelial-mesenchymal transition of colorectal cancer cells.Int J Oncol.2018;53:567-78. https://doi.org/10.3892/ijo.2018.4429 [PubMed]
- 45. Hayes J, Peruzzi PP, Lawler S. MicroRNAs in cancer: biomarkers, functions and therapy.Trends Mol Med.2014;20:460-9. https://doi.org/10.1016/j.molmed.2014.06.005 [PubMed]
- 46. Zhang Z, Hong Y, Xiang D, Zhu P, Wu E, Li W, Mosenson J, Wu WS. MicroRNA-302/367 cluster governs hESC self-renewal by dually regulating cell cycle and apoptosis pathways.Stem Cell Reports.2015;4:645-57. https://doi.org/10.1016/j.stemcr.2015.02.009 [PubMed]
- 47. Zhao L, Han T, Li Y, Sun J, Zhang S, Liu Y, Shan B, Zheng D, Shi J. The lncRNA SNHG5/miR-32 axis regulates gastric cancer cell proliferation and migration by targeting KLF4.FASEB J.2017;31:893-903. https://doi.org/10.1096/fj.201600994R [PubMed]
- 48. Xu J, Lin H, Li G, Sun Y, Chen J, Shi L, Cai X, Chang C. The miR-367-3p Increases Sorafenib Chemotherapy Efficacy to Suppress Hepatocellular Carcinoma Metastasis through Altering the Androgen Receptor Signals.EBioMedicine.2016;12:55-67. https://doi.org/10.1016/j.ebiom.2016.07.013 [PubMed]
- 49. Cai H, Liu X, Zheng J, Xue Y, Ma J, Li Z, Xi Z, Li Z, Bao M, Liu Y. Long non-coding RNA taurine upregulated 1 enhances tumor-induced angiogenesis through inhibiting microRNA-299 in human glioblastoma.Oncogene.2017;36:318-31. https://doi.org/10.1038/onc.2016.212 [PubMed]
- 50. Zhang J, Fan J, Zhou C, Qi Y. miR-363-5p as potential prognostic marker for hepatocellular carcinoma indicated by weighted co-expression network analysis of miRNAs and mRNA.BMC Gastroenterol.2017;17:81 https://doi.org/10.1186/s12876-017-0637-2 [PubMed]
- 51. Wang F, Dai M, Chen H, Li Y, Zhang J, Zou Z, Yang H. Prognostic value of hsa-mir-299 and hsa-mir-7706 in hepatocellular carcinoma.Oncol Lett.2018;16:815-20. https://doi.org/10.3892/ol.2018.8710 [PubMed]
- 52. Li LM, Hu ZB, Zhou ZX, Chen X, Liu FY, Zhang JF, Shen HB, Zhang CY, Zen K. Serum microRNA profiles serve as novel biomarkers for HBV infection and diagnosis of HBV-positive hepatocarcinoma.Cancer Res.2010;70:9798-807. https://doi.org/10.1158/0008-5472.CAN-10-1001 [PubMed]
- 53. Ye J, Zhang W, Liu S, Liu Y, Liu K. miR-363 inhibits the growth, migration and invasion of hepatocellular carcinoma cells by regulating E2F3.Oncol Rep.2017;38:3677-84. https://doi.org/10.3892/or.2017.6018 [PubMed]
- 54. Wang Y, Liu Z, Yao B, Li Q, Wang L, Wang C, Dou C, Xu M, Liu Q, Tu K. Long non-coding RNA CASC2 suppresses epithelial-mesenchymal transition of hepatocellular carcinoma cells through CASC2/miR-367/FBXW7 axis.Mol Cancer.2017;16:123 https://doi.org/10.1186/s12943-017-0702-z [PubMed]
- 55. Yang D, Ma M, Zhou W, Yang B, Xiao C. Inhibition of miR-32 activity promoted EMT induced by PM2.5 exposure through the modulation of the Smad1-mediated signaling pathways in lung cancer cells.Chemosphere.2017;184:289-98. https://doi.org/10.1016/j.chemosphere.2017.05.152 [PubMed]
- 56. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, Barrette TR, Anstet MJ, Kincead-Beal C, Kulkarni P, Varambally S, Ghosh D, Chinnaiyan AM. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles.Neoplasia.2007;9:166-80. https://doi.org/10.1593/neo.07112 [PubMed]
- 57. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses.Neoplasia.2017;19:649-58. https://doi.org/10.1016/j.neo.2017.05.002 [PubMed]
- 58. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, Cerami E, Sander C, Schultz N. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal.Sci Signal.2013;6:pl1 https://doi.org/10.1126/scisignal.2004088 [PubMed]
- 59. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters.OMICS.2012;16:284-7. https://doi.org/10.1089/omi.2011.0118 [PubMed]
- 60. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types.Nucleic Acids Res.2018;46:D956-D63. https://doi.org/10.1093/nar/gkx1090 [PubMed]
- 61. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017:a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit.Nucleic Acids Res.2017;45:W130-W7. https://doi.org/10.1093/nar/gkx356 [PubMed]
- 62. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0.Bioinformatics.2011;27:1739-40. https://doi.org/10.1093/bioinformatics/btr260 [PubMed]
- 63. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function.Nucleic Acids Res.2010;38:W214-20. https://doi.org/10.1093/nar/gkq537 [PubMed]