- Open Access
A circRNA–miRNA–mRNA network identification for exploring underlying pathogenesis and therapy strategy of hepatocellular carcinoma
Journal of Translational Medicinevolume 16, Article number: 220 (2018)
Circular RNAs (circRNAs) have received increasing attention in human tumor research. However, there are still a large number of unknown circRNAs that need to be deciphered. The aim of this study is to unearth novel circRNAs as well as their action mechanisms in hepatocellular carcinoma (HCC).
A combinative strategy of big data mining, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and computational biology was employed to dig HCC-related circRNAs and to explore their potential action mechanisms. A connectivity map (CMap) analysis was conducted to identify potential therapeutic agents for HCC.
Six differently expressed circRNAs were obtained from three Gene Expression Omnibus microarray datasets (GSE78520, GSE94508 and GSE97332) using the RobustRankAggreg method. Following the RT-qPCR corroboration, three circRNAs (hsa_circRNA_102166, hsa_circRNA_100291 and hsa_circRNA_104515) were selected for further analysis. miRNA response elements of the three circRNAs were predicted. Five circRNA–miRNA interactions including two circRNAs (hsa_circRNA_104515 and hsa_circRNA_100291) and five miRNAs (hsa-miR-1303, hsa-miR-142-5p, hsa-miR-877-5p, hsa-miR-583 and hsa-miR-1276) were identified. Then, 1424 target genes of the above five miRNAs and 3278 differently expressed genes (DEGs) on HCC were collected. By intersecting the miRNA target genes and the DEGs, we acquired 172 overlapped genes. A protein–protein interaction network based on the 172 genes was established, with seven hubgenes (JUN, MYCN, AR, ESR1, FOXO1, IGF1 and CD34) determined from the network. The Gene Oncology, Kyoto Encyclopedia of Genes and Genomes and Reactome enrichment analyses revealed that the seven hubgenes were linked with some cancer-related biological functions and pathways. Additionally, three bioactive chemicals (decitabine, BW-B70C and gefitinib) based on the seven hubgenes were identified as therapeutic options for HCC by the CMap analysis.
Our study provides a novel insight into the pathogenesis and therapy of HCC from the circRNA–miRNA–mRNA network view.
Circular RNA (circRNA), with a complete closed loop structure, is firstly identified in 1976 . However, these transcripts without poly-A tail are ignored for a long time due to the limitation of traditional RNA detection methods. In recent years, with the development of high-throughput sequencing technology, a myriad of circRNAs has been found in the eukaryotic transcriptome . With the features of cell-type specific  and highly conserved across species , circRNAs are thought to be new star RNAs which play important roles in various diseases, including human cancers [5, 6].
Competing endogenous RNAs (ceRNAs) are transcripts that act as miRNA sponges, modulating each other at post-transcriptional level via competely binding to shared miRNAs . Recently, circRNAs have become new hotspots in ceRNA family since they have been demonstrated to harbor abundant conserved miRNA response elements (MREs) . Increasing study has revealed that some circRNAs are involved in tumor initiation and progression by the ceRNA mechanism. For example, a classic cirRNA, ciRs7, is deemed to be a miRNA sponge, absorbing miR-7 and liberating the latter inhibitory effect on its target gene in many human cancers [9, 10]. CircRNAs as ceRNAs mediating pathological processes has also been reported in HCC [11, 12]. However, many unknown circRNAs still remain to be explored.
In this study, we employed a combinative strategy of gene chip and computational biology to investigate novel circRNAs and their potential action mechanisms in HCC. The flow chart recapitulating the present work is shown in Fig. 1 as follows: First, we collected HCC-related microarray datasets providing expression profile of circRNAs from the Gene Expression Omnibus (GEO), obtaining differently expressed circRNAs (DECs) with RobustRankAggreg method and corroborating their expression using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). To depict whether the DECs function as ceRNAs in HCC, we collected their sponge miRNAs and miRNA target genes, constructing a circRNA–miRNA–mRNA network. A protein–protein interaction (PPI) network was subsequently established and the hubgenes were identified. Then, Gene Oncology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome enrichment analyses on the hubgenes were performed to elucidate the potential pathogenesis of HCC. Furthermore, we conducted a connectivity map (CMap) analysis to acquire bioactive compounds for the treatment of HCC, which provide a new insight into the latent therapeutic capacity of circRNAs in HCC.
Screening of DECs in HCC from GEO
Microarray datasets providing circRNA expression profile in HCC were achieved from the GEO database . All raw expression data were normalized and log2-transformed. First, we used Limma, a Bioconductor package for differential analysis of microarray data, to determine DECs in each dataset with the criteria of |log2(foldchange)| > 1 and P-value < 0.05 . Then, we integrated and ranked all of the DECs with R package RobustRankAggreg .
Validation of DECs with RT-qPCR method
Sixteen paired fresh frozen HCC tissues and corresponding adjacent non-tumor tissues were obtained from patients diagnosed with HCC at the First Affiliated Hospital of Guangxi Medical University (Nanning, Guangxi, People’s Republic of China). No patient received any radiotherapy or chemotherapy before surgery. The Ethics Committee of the First Affiliated Hospital of Guangxi Medical University approved this study.
Total RNA was isolated with TRIzol® Reagent (Life technologies, Thermo Fisher Scientific, USA) following the manufacturer’s instruction. Then, 1 µg total RNA was reversed into 20 µl complementary DNA (cDNA) with Geneseed® II First Strand cDNA Synthesis Kit (Geneseed, Guangzhou, China). RT-qPCR was conducted using Geneseed® qPCR SYBR® Green Master Mix (Geneseed) on ABI7500 system (Applied Biosystems, CA, USA) in line with the manufacturer’s procedure. Beta-actin (Geneseed) was set as the endogenous reference. All of the primer sequences used in this study were synthesized by Geneseed and are displayed in Table 1. CircRNA expression was determined using the 2−ΔCT method. Significance between groups was analyzed by paired T-test with SPSS 22.0 (IBM, New York, USA). A P-value < 0.05 denotes a statistical significance.
Prediction of MREs
The miRNA binding sites, also known as MREs, of those selected DECs, were predicted with two web tools, Cancer-Specific CircRNA (CSCD)  and Circular RNA Interactome (CircInteractome) . We identified overlapped miRNAs of the two algorithms as potential target miRNAs of the DECs.
Verification of miRNA expression based on data from GEO and TCGA
Microarray datasets providing miRNA expression profile in HCC were obtained from the public databases GEO and The Cancer Genome Atlas (TCGA) . The retrieval terms were as follows: (hepatocellular OR hepatic OR liver OR HCC) and (“cancer” OR “tumor” OR “tumour” OR “carcinoma” OR “neoplasm” OR “malignan*”) and (miRNA OR microRNA OR miR OR “non-coding RNA” OR ncRNA OR “noncoding RNA” OR “non coding RNA”).
We screened available datasets based on the listed inclusion criteria: (1) all of the patients were diagnosed with HCC; (2) the studies must contain circRNA expression data both in cancerous and normal liver tissues; and (3) the sample sizes in tumor and non-tumor group were at least three.
Two revivers extracted the basic information of each included record: miRNA type, first author and published year, region, data source, platform, number of case, and expression level of miRNA. Any divergences were settled via discussion with a third investigator.
The combined standard mean difference (SMD) and 95% confidence interval (95% CI) were computed by STATA 12.0 (StataCorp, College Station, TX, USA). A SMD > 0 represents high expression of miRNA in HCC samples than in normal controls. The corresponding 95% CI do not cross 1 and a P-value < 0.05 suggest a statistical significance.
Prediction of miRNA target genes
The miRNA–mRNA interactions were predicted with miRWalk 2.0 , which involves 12 predicted algorithms (Targetscan, RNAhybrid, RNA22, PITA, Pictar2, miRWalk, Microt4, miRNAMap, miRDB, mirbridge, miRanda and miRMap). Target genes forecasted by at least eight algorithms were selected for further analysis.
Collection of differently expressed genes (DEGs) of HCC form TCGA
RNA-sequencing (RNA-seq) data containing 374 HCC samples and 50 normal controls was downloaded from the TCGA. DEGs were determined by the edgeR package  in Bioconductor with the filter criteria of |log2(foldchange)| > 1 and adjust P-value < 0.05.
Construction of circRNA–miRNA–mRNA network
The overlapping genes between the predicted miRNA target genes and the DEGs were obtained for circRNA–miRNA–mRNA network construction. The Cytoscape 3.6.1 software  was used to visualize the regulatory network.
Establishment of PPI network and identification of hub-genes
A PPI network was established by the STRING (v10.5)  and visualized by the Cytoscape 3.6.1. Then, the “Molecular Complex Detection” (MCODE), a clustering algorithm identifying locally densely connected regions in a large PPI network based on a node-weighting arithmetic , was employed to recognize highly interacted hubgene clustering.
GO, KEGG and Reactome enrichment analyses
GO annotation and KEGG pathway analyses were conducted by clusterProfiler, an R package for functional classification and enrichment of gene clusters using hypergeometric distribution . Reactome pathway analysis was performed by Reactome FI, a plugin of Cytoscape for network and pathway analysis .
Hubgenes consisting of two lists of up- and down-regulated tags were uploaded into the CMap web tool, matching against over 7000 gene expression profiles following treatment of 1309 bioactive compounds in human cell lines . The match between the signatures of interest and chemicals from CMap was assessed by a connectivity score from − 1 to 1: a positive score denotes a stimulative effect of compound on the query signatures; while a negative score implicates a repressed effect of a compound on the query signatures.
Identification of six DECs in HCC based on RobustRankAggreg method
Three microarray datasets (GSE78520, GSE94508 and GSE97332) were included in our study. All of the three gene chips were from the platform of Agilent-069978 Arraystar Human CircRNA microarray V1. The basic information of the three datasets is concluded in Table 2. A total of 259 DECs with 211 up-regulated circRNAs and 48 down-regulated circRNAs were found in gene chip GSE78520 (Fig. 2a); 299 DECs with 46 up-regulated circRNAs and 253 down-regulated circRNAs were determined in gene chip GSE94508 (Fig. 2b); 882 DECs with 429 up-regulated circRNAs and 453 down-regulated circRNAs were identified in gene chip GSE97332 (Fig. 2c). Subsequently, we integrated the DECs of the three datasets and ranked them with a robust method. A total of six circRNAs, including two up-regulated circRNAs (hsa_circRNA_103510 and hsa_circRNA_100542) and four down-regulated circRNAs (hsa_circRNA_102166, hsa_circRNA_104515, hsa_circRNA_105031 and hsa_circRNA_100291), were found to be in the top rankings with an adjust P-value < 0.05 (Fig. 3). The expression of the six circRNAs in each dataset is shown in Fig. 4. The essential characteristics of the six circRNAs are displayed in Table 3. The basic structural patterns of the six circRNAs are exhibited in Fig. 5.
Corroboration of the six circRNAs with RT-qPCR
RT-qPCR was conducted in 16 pairs HCC samples and adjacent non-cancerous tissues to corroborate the expression of the six circRNAs. Among the six circRNAs, hsa_circRNA_103510 and hsa_circRNA_100542 cannot be detected in these in-house tissues due to their relative low expression levels. The low expression of hsa_circRNA_102166 in HCC tissues was corroborated by RT-qPCR (P = 0.047, Fig. 6a). Hsa_circRNA_100291 and hsa_circRNA_104515 showed a down-regulated tendency in HCC, though the P-values were greater than 0.05 (hsa_circRNA_100291: P = 0.069, Fig. 6b; hsa_circRNA_104515: P = 0.059, Fig. 6c). For hsa_circRNA_105031, its expression in HCC tissues was similar to that in non-cancerous tissues (P = 0.588, Fig. 6d). The head-to-tail splicing in the RT-qPCR product of hsa_circRNA_102166, hsa_circRNA_100291, hsa_circRNA_104515 and hsa_circRNA_105031 was confirmed by Sanger sequencing (Fig. 7).
Identification of five circRNA–miRNA interactions
Hsa_circRNA_102166, hsa_circRNA_100291 and hsa_circRNA_104515 were selected for further analysis. Increasing evidence has demonstrated that some circRNAs play critical roles in tumors by functioning as “decoys” to sponge miRNAs. To depict whether the three circRNAs perform the similar role in HCC, we collected their potential target miRNAs from two online databases, CSCD and CircInteractome. A total of five circRNA–miRNA interactions including two circRNAs (hsa_circRNA_104515 and hsa_circRNA_100291) and five miRNAs (hsa-miR-1303, hsa-miR-142-5p, hsa-miR-877-5p, hsa-miR-583 and hsa-miR-1276) were identified (Fig. 8). DIANA-miRPath  was exploited to explore the signaling pathways in which the five miRNAs may be involved. As shown in Fig. 9, all of the five miRNAs were closely linked with some cancer-related pathways.
Expression of the five miRNAs based on the data from GEO and TCGA
In all, 11 microarray and RNA-seq datasets were included. Among them, 7 datasets with 339 HCC samples and 113 normal controls were for miR-1303; 11 datasets with 881 HCC samples and 351 normal controls were for miR-142-5p; 7 datasets with 671 HCC samples and 149 normal controls were for miR-877-5p; 5 datasets with 307 HCC samples and 97 normal controls were for miR-583; and 5 datasets with 534 HCC samples and 110 normal controls were for miR-1276. The essential properties of the 11 records are concluded in Table 4. According to the pooled results, miR-142-5p was down-regulated (SMD = − 0.53, 95% CI − 0.88 to − 0.19, P = 0.003; Fig. 10a) and miR-877 was up-regulated (SMD = 0.80, 95% CI 0.25–1.34, P = 0.004; Fig. 10b) in HCC. MiR-1303 and miR-1276 showed high expression tendencies in HCC tissues compared to in normal liver tissues, while the P-values were greater than 0.05 (miR-1303: SMD = 0.49, 95% CI − 0.19 to 1.16, P = 0.157; Fig. 10c; miR-1276: SMD = 0.21, 95% CI − 0.30 to 0.72, P = 0.422; Fig. 10d). For miR-583, its expression in HCC tissues was similar to that in normal liver tissues (SMD = 0.03, 95% CI − 0.25 to 0.31, P = 0.851, Fig. 10e).
Construction of circRNA–miRNA–mRNA network
Total 1424 target genes of the aforementioned five miRNAs were obtained from the miRWalk. Additionally, 3278 DEGs in HCC were gained from the TCGA (Fig. 11a). By intersecting the predicted target genes and DEGs, we identified 172 target genes that exert momentous roles in HCC (Fig. 11b).
We integrated the circRNA–miRNA interactions and miRNA–mRNA interactions to construct a circRNA–miRNA–mRNA network (Fig. 12), which provided a preliminary insight into the links between the DECs (hsa_circRNA_104515 and hsa_circRNA_100291), the five miRNAs (hsa-miR-1303, hsa-miR-142-5p, hsa-miR-877-5p, hsa-miR-583 and hsa-miR-1276) and the 172 mRNAs.
Identification of seven hubgenes with MCODE algorithm from PPI network
Removing unconnected nodes, we established a PPI network consisting of 91 nodes and 131 edges to view the interactions among the 172 mRNAs (Fig. 13a). Considering the importance of hubgene in a network, we employed a MCODE approach to screen hubgenes from the PPI network. With the k-core = 2, one subnetwork with 7 nodes and 18 edges was identified (Fig. 13b), which unveiled the critical roles of the seven genes (JUN, MYCN, AR, ESR1, FOXO1, IGF1 and CD34) in HCC. The expression levels of the seven genes in HCC are exhibited in Fig. 14. A circRNA-miRNA-hubgene network was then built to delineate the links among the DECs, miRNAs and hubgenes (Fig. 15). Eight circRNA–miRNA–mRNA regulatory modules, including hsa_circRNA_100291/hsa-miR-1276/FOXO1 regulatory axis, hsa_circRNA_100291/hsa-miR-583/ESR1 regulatory axis, hsa_circRNA_100291/hsa-miR-583/JUN regulatory axis, hsa_circRNA_100291/hsa-miR-583/AR regulatory axis, hsa_circRNA_104515/hsa-miR-877-5p/AR regulatory axis, hsa_circRNA_104515/hsa-miR-142-5p/MYCN regulatory axis, hsa_circRNA_104515/has-miR-142-5p/IGF1 regulatory axis and has_circRNA_104515/has-miR-142-5p/CD34 axis, were found from the network.
GO annotation, KEGG pathway and Reactome pathway analyses of the seven hubgenes
GO analysis was carried out to illustrate the functional annotations of the seven hubgenes. The top five highly enriched GO terms of biological process (BP), cellular component (CC) and molecular function (MF) are shown in Fig. 16. The most enriched GO terms in BP was “epithelial cell proliferation” (P < 0.0001), that in CC was “chromatin” (P < 0.0001), and that in MF was “beta-catenin binding” (P < 0.0001). KEGG pathway analysis was conducted to ascertain the signaling cascade that the seven genes participate. With an adjust P-value < 0.05, 12 significantly enriched pathways were obtained (Fig. 17a). Among the 12 pathways, “AMPK signaling pathway”, “FoxO signaling pathway” and “Estrogen signaling pathway” are linked with the progression of HCC [27,28,29]. Additionally, some other pathways such as “Prostate cancer”, “Breast cancer” and “Transcriptional misregulation in cancer” were also tumor-related pathways. Reactome pathway analysis was further performed to delineate the metabolic pathways that the seven hubgenes related to. With a P-value < 0.05, a total of 29 pathways were identified (data not shown). The top 10 significantly enriched Reactome pathways are displayed in Fig. 17b.
Identification of three bioactive compounds for the treatment of HCC based on CMap analysis
The seven hubgenes consisting of two up-regulated genes (CD34, MYCN) and five down-regulated genes (AR, JUN, ESR1, FOXO1, IGF1) were loaded into the CMap web tool as up-regulated tags and down-regulated tags, respectively. Following the signature query, three compounds (decitabine, BW-B70C, gefitinib) with the highest negative enrichment score were determined as the potential therapeutic agents for HCC (Table 5). The chemical structures of the three compounds are presented in Fig. 18.
Due to the lack of 5′ caps and 3′polyadenylated tails, circRNAs are ignored in classical polyadenylated transcriptome studies for a long time. In the past few decades, with the development of high-throughput sequencing, biochemical and computational biology methods, a large number of circRNAs are unmasked the veil in various tissues and cells . Increasing study has unveiled the important roles of circRNAs in a myriad of human diseases, including malignant tumors [31,32,33]. CircRNAs usually serve as diagnostic and prognostic biomarkers because of their relative tolerance to exonucleases, which benefits from the covalently closed loop structures [34, 35]. In addition, owing to the high cell- and tissue- specificity of circRNAs, their roles in different neoplasms are not complete accord. In HCC, an increasing number of circRNAs, such as ciRS-7 , has_circ_0067934 , CDR1as  and circMTO1 , have been reported to exert momentous roles in regulating pathophysiological process and guiding clinical diagnosis and treatment. However, there are still a lot of circRNAs that need to be unearthed.
In this study, we firstly collected three gene chips (GSE78520, GSE94508 and GSE97332) from the GEO database and screened six DECs (hsa_circRNA_103510,hsa_circRNA_100542, hsa_circRNA_102166, hsa_circRNA_104515, hsa_circRNA_105031 and hsa_circRNA_100291) with the RobustRankAggreg approach, an algorithm for integration of different gene expression spectra, ranking genes by their fold changes and aggregating these rankings to achieve the final robust and rigorous results. Following the RT-qPCR validation of the six DECs, three circRNAs (hsa_circRNA_102166, hsa_circRNA_100291 and hsa_circRNA_104515) were selected for further analysis.
As highly conserved endogenous RNAs, many circRNAs harbors abundant miRNA binding sites, indicating that they can sponge corresponding miRNAs and thus function as ceRNAs to regulate gene expression [8, 38, 39]. To ascertain whether the aforementioned three circRNAs function as ceNAs in HCC, we predicted their MREs via two online tools, CSCD and CircInteractome. The former web tool predicted MREs within 50 bp upstream and downstream of circRNA junction point . The latter web tool predicted MREs based on the TargetScane algorithm, which forecasts MREs by searching for 7mer or 8mer complementarity to the 3′ end and the seed region of each miRNA [17, 40]. We chose miRNA predicted by both algorithms as the putative target miRNA for the three circRNAs. Finally, five circRNA–miRNA interactions consisting of two circRNAs (hsa_circRNA_104515 and hsa_circRNA_100291) and five miRNAs (hsa-miR-1303, hsa-miR-142-5p, hsa-miR-877-5p, hsa-miR-583 and hsa-miR-1276) were determined. The expression of the five miRNAs was then verified based on the microarray and RNA-seq data from the GEO and TCGA. The results showed that miR-142-5p was down-regulated in HCC, which was consistent with previous studies [41, 42]. No relevant study reported the expression of miR-877, miR-1303, miR-1276 and miR-583 in HCC. According to our results, miR-877 was up-regulated in HCC. MiR-1303 and miR-1276 showed high expression tendencies in HCC, while the differences were not statistical significant. For miR-583, its expression in HCC tissues was similar to that in normal liver tissues. Due to the inter-study heterogeneity that cannot be ignored, the reliable of the pooled results is reduced. Thereby, further rigorous studies are necessary to validate these findings.
Following the collection of the 172 overlapped genes between the target genes of the five miRNAs and the DEGs in HCC, we constructed a circRNA–miRNA–mRNA regulatory network. We found that hsa_circRNA_104515 and hsa_circRNA_100291 may act as ceRNAs to capture miR-1276, miR-583, miR-877-5p or miR-142-5p, and subsequently regulate the 172 genes expression. Our results provide an evidence of the ceRNA regulatory mechanism of hsa_circRNA_104515 and hsa_circRNA_100291 in HCC. To further elucidate the action mechanism of the ceRNA network, we constructed a PPI network, screening seven hubgenes (JUN, MYCN, AR, ESR1, FOXO1, IGF1 and CD34) from the PPI network. The functional annotations and pathway analyses showed that the seven hubgenes genes were involved in many critical tumor-related biological functions and metabolic pathways, such as “epithelial cell proliferation”, “AMPK signaling pathway”, “FoxO signaling pathway”, “Estrogen signaling pathway”, “MAPK6/MAPK4 signaling” and “MAPK family signaling cascades”. The important roles of the seven genes in HCC have also been demonstrated previously [43,44,45,46,47,48,49]. However, their links with circRNAs have not yet been explored. Here, we identified seven circRNA–miRNA–mRNA axes (hsa_circRNA_100291/miR-1276/FOXO1, hsa_circRNA_100291/miR-583/ESR1, hsa_circRNA_100291/miR-583/JUN, hsa_circRNA_100291/miR-583/AR, hsa_circRNA_104515/miR-877-5p/AR, hsa_circRNA_104515/miR-142-5p/MYC, hsa_circRNA_104515/miR-142-5p/MYCN, hsa_circRNA_104515/miR-142-5p/IGF1 and hsa_circRNA_104515/miR-142-5p/CD34), indicating competitive regulatory relationships of hsa_circRNA_100291 and hsa_circRNA_104515 with the seven genes in HCC. However, given that the results are on the basis of computational biology, further in-depth studies are indispensable to verify the possible roles of the seven axes in HCC.
To date drug control is an important treatment for patients with HCC . Digging effective and sensitive drugs against HCC helps improve patients’ outcomes. We therefore implemented the CMap analysis of the seven hubgenes to explore usable drugs for the treatment of HCC. Based on the genome-wide expression profiling of transcripts technology, CMap provides a comprehensive and accurate data resource for exploration of novel drug or relocation of existing drug . Drugs available in CMap are all licensed for human use by the Food and Drug Administration, thus it is an ideal and reliable approach to obtain therapeutic agents for human diseases . Three chemicals (decitabine, BW-B70C and gefitinib) were identified as the treatment options for HCC. As a cytidine antimetabolite analogue, decitabine represses DNA methylation, arresting cells into G1/S phase and inhibiting tumor cell proliferation. Its antitumor activity in solid tumors including HCC has been elucidated in previous studies . Xing et al.  have demonstrated that decitabine could facilitate the expression of miR-122 via inhibition of methylation in HCC cells. More importantly, Skårn et al.  have demonstrated that miR-142 is epigenetically restrained by DNA methylation. Treating decitabine in mesenchymal cells promotes mature miR-142-5p/3p expression and thus depresses cell proliferation. We thus hypothesize that decitabine may exert its anti-HCC effect by augmenting miR-142-5p via demethylation and subsequently regulating the downstream target genes of miR-142-5p. Further well-design study is necessary to validate the conclusion. BW-B70C is an inhibitor of arachidonic acid 5-lipoxygenase. Previous study has reported its anti-neoplastic activity in leukemic cells by suppressing the NOTCH1-P13K-AKT-eNOS axis . However, its anti-HCC effect and action mechanism have not been elucidated as of yet. In this study, we found its potential as therapeutic agent for HCC. More studies are needed to verify this finding. Gefitinib is a selective inhibitor of tyrosine kinase receptor used in clinic for the treatment of locally advanced or metastatic non-small cell lung cancer (NSCLC) . Its antineoplastic effect on the other solid tumors including HCC has also been reported [58, 59]. However, the responsiveness of different HCC patients to gefitinib varies greatly and most patients even develop gefitinib resistance . Multiple non-coding RNAs including miRNAs and lncRNAs have been reported to contribute to gefitinib resistance and sensitivity in NSCLC [61,62,63]. We speculate that gefitinib share the similar resistance mechanisms in HCC and lung cancer. Our study provides a theoretical basis for studying gefitinib resistance mechanism and enhancing gefitinib sensitivity in patients with HCC from the perspective of circRNA–miRNA–mRNA network.
In conclusion, by employing a comprehensive strategy of big data mining, RT-qPCR and computational biology, we constructed a circRNA–miRNA–mRNA network and found that hsa_circRNA_100291 and hsa_circRNA_104515 may function as ceRNAs to exert critical roles in HCC. In addition, three bioactive chemicals (decitabine, BW-B70C and gefitinib) based on the CMap analysis was determined as therapeutic agents for HCC. Our study provides a novel insight into the pathogenesis and therapy for HCC from the circRNA–miRNA–mRNA view.
reverse transcription-quantitative polymerase chain reaction
competing endogenous RNA
miRNA response element
Gene Expression Omnibus
differently expressed circRNA
Kyoto Encyclopedia of Genes and Genomes
Circular RNA Interactome
The Cancer Genome Atlas
standard mean difference
- 95% CI:
95% confidence interval
Molecular Complex Detection
differently expressed gene
Food and Drug Administration
non-small cell lung cancer
Sanger HL, Klotz G, Riesner D, Gross HJ, Kleinschmidt AK. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc Natl Acad Sci U S A. 1976;73:3852–6.
Wang PL, Bao Y, Yee MC, Barrett SP, Hogan GJ, Olsen MN, et al. Circular RNA is expressed across the eukaryotic tree of life. PLoS ONE. 2014;9:e90859.
Salzman J, Chen RE, Olsen MN, Wang PL, Brown PO. Cell-type specific features of circular RNA expression. PLoS Genet. 2013;9:e1003777.
Vidal AF, Sandoval GT, Magalhaes L, Santos SE, Ribeiro-dos-Santos A. Circular RNAs as a new field in gene regulation and their implications in translational research. Epigenomics. 2016;8:551–62.
Wang Y, Mo Y, Gong Z, Yang X, Yang M, Zhang S, et al. Circular RNAs in human cancer. Mol Cancer. 2017;16:25.
Kristensen LS, Hansen TB, Veno MT, Kjems J. Circular RNAs in cancer: opportunities and challenges in the field. Oncogene. 2018;37:555–65.
Qi X, Zhang DH, Wu N, Xiao JH, Wang X, Ma W. ceRNA in cancer: possible functions and clinical implications. J Med Genet. 2015;52:710–8.
Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer. 2018;17:79.
Peng L, Yuan XQ, Li GC. The emerging landscape of circular RNA ciRS-7 in cancer (Review). Oncol Rep. 2015;33:2669–74.
Geng HH, Li R, Su YM, Xiao J, Pan M, Cai XX, et al. The circular RNA Cdr1as promotes myocardial infarction by mediating the regulation of miR-7a on its target genes expression. PLoS ONE. 2016;11:e0151753.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66:1151–64.
Zhu Q, Lu G, Luo Z, Gui F, Wu J, Zhang D, et al. CircRNA circ_0067934 promotes tumor growth and metastasis in hepatocellular carcinoma through regulation of miR-1324/FZD5/Wnt/beta-catenin axis. Biochem Biophys Res Commun. 2018;497:626–32.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–5.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28:573–80.
Xia S, Feng J, Chen K, Ma Y, Gong J, Cai F, et al. CSCD: a database for cancer-specific circular RNAs. Nucleic Acids Res. 2018;46:D925–9.
Dudekula DB, Panda AC, Grammatikakis I, De S, Abdelmohsen K, Gorospe M. CircInteractome: a web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016;13:34–42.
Deng M, Bragelmann J, Schultze JL, Perner S. Web-TCGA: an online platform for integrated analysis of molecular cancer data sets. BMC Bioinform. 2016;17:72.
Dweep H, Gretz N, Sticht C. miRWalk database for miRNA-target interactions. Methods Mol Biol. 2014;1182:289–305.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinform. 2014;47:8.13.1–8.1324.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003;4:2.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Wu G, Dawson E, Duong A, Haw R, Stein L. ReactomeFIViz: a Cytoscape app for pathway and network-based data analysis. F1000Res. 2014;3:146.
Musa A, Ghoraie LS, Zhang SD, Glazko G, Yli-Harja O, Dehmer M, et al. A review of connectivity map and computational approaches in pharmacogenomics. Brief Bioinform. 2017;18:903.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43:W460–6.
Huang GM, Jiang QH, Cai C, Qu M, Shen W. SCD1 negatively regulates autophagy-induced cell death in human hepatocellular carcinoma through inactivation of the AMPK signaling pathway. Cancer Lett. 2015;358:180–90.
Hou YQ, Yao Y, Bao YL, Song ZB, Yang C, Gao XL, et al. Juglanthraquinone C induces intracellular ROS increase and apoptosis by activating the Akt/Foxo signal pathway in HCC cells. Oxid Med Cell Longev. 2016;2016:4941623.
You H, Meng K, Wang ZY. The ER-alpha36/EGFR signaling loop promotes growth of hepatocellular carcinoma cells. Steroids. 2018;134:78–87.
Chen LL, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12:381–8.
Yang Y, Gao X, Zhang M, Yan S, Sun C, Xiao F, Huang N, Yang X, Zhao K, Zhou H, Huang S, Xie B, Zhang N. Novel role of FBXW7 circular RNA in repressing glioma tumorigenesis. J Natl Cancer Inst. 2018. https://doi.org/10.1093/jnci/djx166.
Hu W, Bi ZY, Chen ZL, Liu C, Li LL, Zhang F, Zhou Q, Zhu W, Song YY, Zhan BT, Zhang Q, Bi YY, Sun CC, Li DJ. Emerging landscape of circular RNAs in lung cancer. Cancer Lett. 2018;427:18–27.
Li P, Yang X, Yuan W, Yang C, Zhang X, Han J, Wang J, Deng X, Yang H, Li P, Tao J, Lu Q, Gu M. CircRNA-Cdr1as exerts anti-oncogenic functions in bladder cancer by sponging MicroRNA-135a. Cell Physiol Biochem. 2018;46:1606–16.
Li P, Chen S, Chen H, Mo X, Li T, Shao Y, Xiao B, Guo J. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clin Chim Acta. 2015;444:132–6.
Jiang XM, Li ZL, Li JL, Xu Y, Leng KM, Cui YF, Sun DJ. A novel prognostic biomarker for cholangiocarcinoma: circRNA Cdr1as. Eur Rev Med Pharmacol Sci. 2018;22:365–71.
Xu L, Zhang M, Zheng X, Yi P, Lan C, Xu M. The circular RNA ciRS-7 (Cdr1as) acts as a risk factor of hepatic microvascular invasion in hepatocellular carcinoma. J Cancer Res Clin Oncol. 2017;143:17–27.
Yang X, Xiong Q, Wu Y, Li S, Ge F. Quantitative proteomics reveals the regulatory networks of circular RNA CDR1as in hepatocellular carcinoma cells. J Proteome Res. 2017;16:3891–902.
Zhang Y, Liang W, Zhang P, Chen J, Qian H, Zhang X, et al. Circular RNAs: emerging cancer biomarkers and targets. J Exp Clin Cancer Res. 2017;36:152.
Cortes-Lopez M, Miura P. Emerging functions of circular RNAs. Yale J Biol Med. 2016;89:527–37.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120:15–20.
Lou K, Chen N, Li Z, Zhang B, Wang X, Chen Y, et al. MicroRNA-142-5p overexpression inhibits cell growth and induces apoptosis by regulating FOXO in hepatocellular carcinoma cells. Oncol Res. 2017;25:65–73.
Tsang FH, Au SL, Wei L, Fan DN, Lee JM, Wong CC, et al. MicroRNA-142-3p and microRNA-142-5p are downregulated in hepatocellular carcinoma and exhibit synergistic effects on cell motility. Front Med. 2015;9:331–43.
Fu L, Balasubramanian M, Shan J, Dudenhausen EE, Kilberg MS. Auto-activation of c-JUN gene by amino acid deprivation of hepatocellular carcinoma cells reveals a novel c-JUN-mediated signaling pathway. J Biol Chem. 2011;286:36724–38.
Hsu TY, Fourel G, Etiemble J, Tiollais P, Buendia MA. Integration of hepatitis virus DNA near c-myc in woodchuck hepatocellular carcinoma. Gastroenterol Jpn. 1990;25(Suppl 2):43–8.
Ma WL, Jeng LB, Lai HC, Liao PY, Chang C. Androgen receptor enhances cell adhesion and decreases cell migration via modulating beta1-integrin-AKT signaling in hepatocellular carcinoma cells. Cancer Lett. 2014;351:64–71.
Hishida M, Nomoto S, Inokawa Y, Hayashi M, Kanda M, Okamura Y, et al. Estrogen receptor 1 gene as a tumor suppressor gene in hepatocellular carcinoma detected by triple-combination array analysis. Int J Oncol. 2013;43:88–94.
Yang XW, Shen GZ, Cao LQ, Jiang XF, Peng HP, Shen G, et al. MicroRNA-1269 promotes proliferation in human hepatocellular carcinoma via downregulation of FOXO1. BMC Cancer. 2014;14:909.
Lei T, Ling X. IGF-1 promotes the growth and metastasis of hepatocellular carcinoma via the inhibition of proteasome-mediated cathepsin B degradation. World J Gastroenterol. 2015;21:10137–49.
Paschoal JP, Bernardo V, Canedo NH, Ribeiro OD, Caroli-Bottino A, Pannain VL. Microvascular density of regenerative nodule to small hepatocellular carcinoma by automated analysis using CD105 and CD34 immunoexpression. BMC Cancer. 2014;14:72.
zniewski P, Tournigand C. Unravelling the pharmacologic opportunities and future directions for targeted therapies in gastro-intestinal cancers part 2: neuroendocrine tumours, hepatocellular carcinoma, and gastro-intestinal stromal tumours. Pharmacol Ther. 2018;181:49–75.
Qu XA, Rajpal DK. Applications of connectivity map in drug discovery and development. Drug Discov Today. 2012;17:1289–98.
Lamb J. The connectivity map: a new tool for biomedical research. Nat Rev Cancer. 2007;7:54–60.
Tao SF, Zhang CS, Guo XL, Xu Y, Zhang SS, Song JR, et al. Anti-tumor effect of 5-aza-2′-deoxycytidine by inhibiting telomerase activity in hepatocellular carcinoma cells. World J Gastroenterol. 2012;18:2334–43.
Xing TJ, Xu HT, Yu WQ, Jiang DF. Methylation regulation of liver-specific microRNA-122 expression and its effects on the proliferation and apoptosis of hepatocellular carcinoma cells. Genet Mol Res. 2013;12:3588–97.
Skarn M, Baroy T, Stratford EW, Myklebost O. Epigenetic regulation and functional characterization of microRNA-142 in mesenchymal cells. PLoS ONE. 2013;8:e79231.
Villegas SN, Gombos R, Garcia-Lopez L, Gutierrez-Perez I, Garcia-Castillo J, Vallejo DM, et al. PI3K/Akt cooperates with oncogenic notch by inducing nitric oxide-dependent inflammation. Cell Rep. 2018;22:2541–9.
Rahman AF, Korashy HM, Kassem MG. Gefitinib. Profiles Drug Subst Excip Relat Methodol. 2014;39:239–64.
Gu HR, Park SC, Choi SJ, Lee JC, Kim YC, Han CJ, et al. Combined treatment with silibinin and either sorafenib or gefitinib enhances their growth-inhibiting effects in hepatocellular carcinoma cells. Clin Mol Hepatol. 2015;21:49–59.
Shao J, Xu Z, Peng X, Chen M, Zhu Y, Xu L, et al. Gefitinib synergizes with irinotecan to suppress hepatocellular carcinoma via antagonizing Rad51-mediated DNA-repair. PLoS ONE. 2016;11:e0146968.
Rossi L, Zoratto F, Papa A, Iodice F, Minozzi M, Frati L, et al. Current approach in the treatment of hepatocellular carcinoma. World J Gastrointest Oncol. 2010;2:348–59.
Yu F, Liu JB, Wu ZJ, Xie WT, Zhong XJ, Hou LK, et al. Tumor suppressive microRNA-124a inhibits stemness and enhances gefitinib sensitivity of non-small cell lung cancer cells by targeting ubiquitin-specific protease 14. Cancer Lett. 2018;427:74–84.
Cao X, Lai S, Hu F, Li G, Wang G, Luo X, et al. miR-19a contributes to gefitinib resistance and epithelial mesenchymal transition in non-small cell lung cancer cells by targeting c-Met. Sci Rep. 2017;7:2939.
Wang Z, Pan L, Yu H, Wang Y. The long non-coding RNA SNHG5 regulates gefitinib resistance in lung adenocarcinoma cells by targeting miR-377/CASP1 axis. Biosci Rep. 2018. https://doi.org/10.1042/BSR20180400.
DDX and YWD collected circRNA-related datasets; screened differently expressed circRNAs; conducted RT-qPCR; collected and analyzed miRNA-related data from TCGA and GEO; and were the major contributors in writing the manuscript. PL and DYW collected miRNA target genes as well as differently expressed genes in HCC; construed PPI network; and conducted GO, KEGG and Reactome analyses. RQH conducted CMap analysis and checked all of the data used in the manuscript. ZBF and GC guided the design of all of the experiments and the writing of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The Ethical Committee of the First Affiliated Hospital of Guangxi Medical University approved the present study.
This study was supported by the National Natural Science Foundation of China (NSFC81560489 and NSFC81560386), the Natural Science Foundation of Guangxi, China (2017GXNSFAA198017), the Guangxi Medical University Training Program for Distinguished Young Scholars (2017) and the Medical Excellence Award Funded by the Creative Research Development Grant from the First Affiliated Hospital of Guangxi Medical University (Gang Chen and Yi-wu Dang).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.