- Primary research
- Open Access
Three novel circRNAs upregulated in tissue and plasma from hepatocellular carcinoma patients and their regulatory network
Cancer Cell International volume 21, Article number: 72 (2021)
- The Correction to this article has been published in Cancer Cell International 2021 21:86
The regulatory roles of circular RNAs (circRNAs) in tumorigenesis have attracted increasing attention. However, novel circRNAs with the potential to be used as serum/plasma biomarkers and their regulatory mechanism in the pathogenesis of hepatocellular carcinoma (HCC) remain explored.
CircRNA expression profiles of tumor tissues and plasma samples from HCC patients were compiled and jointly analyzed. CircRNA–miRNA–mRNA interactions were predicted by bioinformatics tools. The expression of interacting miRNAs and mRNA was verified in independent datasets. Survival analysis and pathway enrichment analysis were conducted on hub genes.
We identified three significantly up-regulated circRNAs (hsa_circ_0009910, hsa_circ_0049783, and hsa_circ_0089172) both in HCC tissues and plasma samples. Two of them were validated to be indeed circular and could be excreted from hepatoma cells. We further revealed four miRNAs (hsa-miR-455-5p, hsa-miR-615-3p, hsa-miR-18a-3p, hsa-miR-4524a-3p) that targeting circRNAs and expressed in human HCC samples, and 95 mRNAs targeted by miRNAs and significantly up-regulated in two HCC cohorts. A protein-protein interaction network revealed 19 hub genes, 12 of them (MCM6, CCNB1, CDC20, NDC80, ZWINT, ASPM, CENPU, MCM3, MCM5, ECT2, CDC7, and DLGAP5) were associated with reduced survival in two HCC cohorts. KEGG, Reactome, and Wikipathway enrichment analysis indicated that the hub genes mainly functioned in DNA replication and cell cycle.
Our study uncovers three novel deregulated circRNAs in tumor and plasma from HCC patients and provides an insight into the pathogenesis from the circRNA–miRNA–mRNA regulatory network.
Hepatocellular carcinoma (HCC) accounts for approximately 90 % of primary liver cancer, the sixth most commonly diagnosed tumor and the fourth leading cause of cancer-related death worldwide . Due to the early asymptomatic, most HCC patients are diagnosed at advanced stages with metastasis, thus losing the opportunity for radical resection [2, 3]. With the improvement of diagnosis and treatment technologies, the survival of HCC patients has increased. However, the 5-year survival rate is only 18 % as high frequencies of tumor metastasis and recurrence . To improve patients’ survival and quality of life, it is urgent to explore novel biomarkers for early diagnosis, prognosis, and treatment of liver cancer [2, 4, 5].
Circular RNA (circRNA) is a type of stable transcripts that do not have a free 5′ end cap and 3′ end poly (A) tail and forms a circular structure by covalent bonds [6, 7]. CircRNA was initially considered to be a by-product of splicing. Nevertheless, with the advancement of high-throughput sequencing technology and bioinformatics analysis, the cellular functions and potential biomedical applications of circRNAs become a new research hotspot . CircRNAs are abundant in the cytoplasm of eukaryotic cells and cell-type specific [9,10,11]. Recent studies have shown that circRNAs are essential to modulate gene expression in the nucleus, act as decoys for miRNAs and proteins, and serve as templates for translation [3, 12,13,14,15]. The regulatory roles of circRNAs in tumorigenesis have also attracted more attention [5, 16, 17]. Increasing studies suggest that circRNAs play essential roles in the occurrence and development of HCC, and are expected to have the potential to become diagnostic biomarkers and therapeutic targest [18, 19]. However, novel circRNAs with the potential to be used as serum/plasma biomarkers and their regulatory mechanism in liver cancer pathogenesis remain to be explored.
In this study, we jointly analyzed three HCC tissues-related microarray datasets from the Gene Expression Omnibus (GEO). Such a relatively large sample size allows for higher statistical power in the identification of differentially expressed circRNAs. CircRNAs that were also increased in HCC plasma were obtained by overlapping with an HCC plasma samples-related microarray dataset. Subsequently, a circRNA-miRNA-mRNA network was constructed to provide an insight into the potential pathogenesis of the identified circRNAs in HCC. The flow chart recapitulating the present work is shown in Fig. 1.
CircRNA expression processing of HCC microarray datasets
The four publicly available HCC-related circRNA microarray datasets used in this study were downloaded from the GEO database, including three gene expression profiles on platform GPL19978 and one from the platform GPL26925. Data normalization and combination were performed by pooling the three datasets of the GPL19978 platform using the Bioconductor Limma package. The principal component analysis was performed to ensure a minimal batch effect in the meta-cohort. Then, differential analysis of the pooled expression data was also conducted by Limma with the criteria of |log2(foldchange)| ≥ 0.585 and P-value < 0.05. Subsequently, the processed plasma circRNA data of GSE135806 based on the thresholds of |log2(foldchange)| ≥ 1 and P-value < 0.05 was downloaded from the GEO, and further intersected with the differentially expressed circRNAs in HCC tissues with an online tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
Cells and cell culture
Human HCC cell lines HepG2, Huh7, and Hep3B were purchased from the Cell Bank of Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences. Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Gibco) supplemented with 10 % fetal bovine serum (Gibco) and 1 % penicillin/streptomycin at 37 °C in a humidified 5 % CO2 incubator.
The culture medium of hepatoma cells was collected, centrifuged at 1,000g for 10 min and 14,000g for 2 min to remove cells and debris . Total RNA from cells or culture medium was extracted using Total RNA Kit I (OMEGA) or TRIzol LS Reagent (Invitrogen) according to the manufacturer’s instructions. cDNA was synthesized from total RNA treated with or without RNase R exonuclease (Geneseed) for real-time PCR on a CFX96 system (Bio-Rad). PCR products were further visualized by agarose gel electrophoresis. β-actin was used as an endogenous control. Primer sequences were shown in Table 1.
Prediction of miRNA binding and miRNA target genes
The Cancer-specific circRNAs database (CSCD)  and circbank (http://www.circbank.cn/index.html) were used to predict miRNA response elements (MRE or binding sites). Overlapped miRNAs in the two databases were considered as potential target miRNAs of circRNAs. The expression of the candidate miRNAs was further verified base on the miRNA-seq dataset GSE76903 . The miRNA–mRNA interactions were predicted with miRWalk , with only the results from binding to 3′ UTR filtered.
Collection of differently expressed genes in HCC
RNA-Seq data of the TCGA-LIHC project consisting of 374 liver cancers and 50 normal tissues were downloaded from the GDC Data Portal. The DESeq2 package was used to screen differentially expressed genes with thresholds of |log 2 (fold change)| ≥ 1 and P-value < 0.01. With the same filter criteria, differently expressed genes from expression profiling by an array (GSE14520, GPL3921 Affymetrix HT Human Genome U133A Array) containing 225 HBV-related HCC and 220 non-tumor tissues were also determined.
Establishment of circRNA–miRNA–mRNA network and identification of hub genes
A circRNA–miRNA–mRNA regulatory network was established by integrating the circRNA-miRNA and miRNA-mRNA interactions and visualized using the Cytoscape 3.7.2 software . The candidate mRNAs were imported to the STRING database (https://string-db.org), and protein-protein interaction (PPI) network was created. Then, MCODE, a clustering algorithm that detects densely connected regions in large protein interaction networks , in Cytoscape was used to identify hub genes in the PPI network. The highest-scoring nodes recognized as hub genes are listed in Table 2.
Survival analysis of hub genes
Survival data of patients from TCGA-LIHC and patients from GSE14520 were obtained from UCSC Xena (https://xenabrowser.net/) and GEO. Patients were divided into two groups based on the median expression values for TCGA-LIHC and based on the minimum P-value cut off values for GSE14520. The Kaplan-Meier’s survival curves for the hub genes were visualized using GraphPad Prism.
Functional enrichment analysis of hub genes
KEGG, Reactome, and Wiki pathway enrichment analysis of the 12 hub genes were carried out using Webgestalt online tool (http://www.webgestalt.org/). A P-value of less than 0.05 was set as the cut-off criterion.
Statistical analysis was performed using GraphPad Prism (version: 8.0.2) and R software. The survival curves were calculated using the Kaplan-Meier method, and a log-rank test assessed the differences. Statistical significance was indicated by P values less than 0.05. *P < 0.05, ** P < 0.01, *** P < 0.001, **** P < 0.0001.
Three circRNAs were significantly up‐regulated both in HCC tissues and plasma samples
To identify differentially expressed circRNAs in HCC tissues and plasma samples, four publicly available microarray datasets providing circRNA expression profile (GSE78520, GSE94508, GSE97332, and GSE135806) were collected in this study. GSE78520, GSE94508, and GSE97332 were HCC tissues-related and were all from platform GPL19978 Agilent-069978 Arraystar Human CircRNA microarray V1, whereas GSE135806 was HCC plasma samples-related and was from the platform of GPL26925 Agilent-084217 CapitalBio Technology Human CircRNA Array v2. A summary of the four datasets was presented in Table 3. First, we performed a quality check, cohort-specific normalization, and combined the three HCC tissues-related datasets by Limma, resulting in a meta-cohort of 15 pairs of HCC and matched non-tumor tissues (GSE78520&GSE94508&GSE97332). The principal component analysis plot indicated that the batch effect due to cohort and processing centers was effectively minimized while retaining the disease state’s differences in the meta-cohort (Fig. 2a, b). A total of 687 differentially expressed circRNAs (421 up-regulated and 266 down-regulated) were identified in the pooled 15 pairs of samples (Fig. 2c, d). Previously reported differentially expressed circRNAs, including hsa_circ_0004913, hsa_circRNA_0008514, hsa_circ_0067934, hsa_circRNA_0006461, hsa_circ_0017639, and hsa_circ_0074854 [26, 27] were also identified in our analysis (Fig. 2c), proving the robustness of the above HCC meta-cohort compilation. On the other hand, a total of 371 circRNAs were differentially expressed (326 up-regulated and 45 down-regulated) in HCC plasma samples (GSE135806 ). Venn diagrams and volcano plots showed three circRNAs (hsa_circ_0009910, hsa_circ_0049783, and hsa_circ_0089172) were consistently up-regulated, and none were down-regulated in the integrated datasets (Fig. 2d, e). Scatter plots confirmed the increased expression of these three circRNAs in HCC tissues and plasma samples compared with their respective controls (Fig. 2f). The essential characteristics of the three identified circRNAs are displayed in Table 4. Intriguingly, regulation directions of their parent genes were different from the significantly up-regulated circRNAs, indicating that the three increased circRNAs in HCC were not by-products of splicing but were suggestive of functionality .
We used divergent primers specifically amplifying circRNA in cDNA but not genomic DNA, and convergent primers amplifying both circRNA and linear RNA and in cDNA and parental gene in genomic DNA to confirm that hsa_circ_0009910 and hsa_circ_0089172 were indeed circular (Fig. 2g). Following RNase R treatment, expression levels of hsa_circ_0009910 and hsa_circ_0089172 were almost unaffected. In contrast, linear β-actin was substantially digested, demonstrating that hsa_circ_0009910 and hsa_circ_0089172 were resistant to RNase R treatment and were stable transcripts (Fig. 2h). Agarose gel electrophoresis following circRNAs amplification showed that these two circRNAs existed in the culture medium of tumor cells but not in the control group, indicating hsa_circ_0009910 and hsa_circ_0089172 excretion from hepatoma cells (Fig. 2i).
Identification of four circRNA–miRNA interactions
To depict whether the three up-regulated circRNAs could regulate gene expression by functioning as miRNA sponges in HCC, we collected their basic structural patterns from the CSCD database, exhibiting numerous miRNA response elements (MREs; Fig. 3a). We further collected their potential miRNA bindings by intersecting the MREs predicted by two web tools, CSCD and circbank (Fig. 3b). A total of five miRNAs (hsa-miR-455-5p, hsa-miR-615-3p, hsa-miR-3926, hsa-miR-5197-3p, and hsa-miR-6836-3p) were identified for hsa_circ_0009910, two miRNAs (hsa-miR-18a-3p and hsa-miR-8071) for hsa_circ_0049783, and three miRNAs (hsa-miR-4524a-3p, hsa-miR-3154, hsa-miR-3190-5p) for hsa_circ_0089172. Among these miRNAs, 4 (hsa-miR-455-5p, hsa-miR-615-3p, hsa-miR-18a-3p, hsa-miR-4524a-3p) are expressed in primary HCC or matched adjacent normal tissues (average normalized count ≥ 1) according to a recent study providing genome-wide miRNA expression profile (Fig. 3c; Table 5). Thus, a total of four circRNA–miRNA interactions, including three circRNAs and four miRNAs were identified.
Construction of circRNA–miRNA–mRNA network
A total of 12,908 target genes of the four miRNAs mentioned above were obtained from miRWalk, a web tool for predicting miRNA binding sites. Additionally, 6648 genes were significantly up-regulated in tumors compared with normal tissues from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) with mixed Asian/European origin, while 337 increased genes from GSE14520 dataset from a Chinese cohort (Fig. 4a). To further narrow the scope, we intersected the miRNA targets predicted by miRWalk with the significantly up-regulated genes determined from TCGA-LIHC and GSE14520, resulting in a total of 95 candidate genes for further analysis (Fig. 4b). By integrating the circRNA-miRNA and miRNA-mRNA interactions, we established and visualized a circRNA–miRNA–mRNA regulatory network in HCC (Fig. 4c).
Identification of 19 hub genes
Figuring out the biological function of potential circRNA-related genes could predict the corresponding function of circRNAs. Interactions among the 95 candidate genes were predicted with the STRING tool, leading to a network containing 67 nodes and 358 edges after removing unconnected nodes (Fig. 5a). Subsequently, we employed the MCODE clustering algorithm to find hub genes from the interaction network. One subnetwork with 19 significantly up-regulated nodes and 166 edges was identified with the k-core = 2 (Fig. 5b, c, and Table 2), suggesting the 19 hub genes’ critical roles in HCC.
Twelve hub genes with prognostic values
To investigate the prognostic values of the hub genes, survival analysis was performed on patients from TCGA-LIHC and GSE14520. Kaplan-Meier’s survival curves showed that among the 19 hub genes, high expression of 12 genes, including MCM6, CCNB1, CDC20, NDC80, ZWINT, ASPM, CENPU, MCM3, MCM5, ECT2, CDC7, and DLGAP5 was significantly associated with poorer overall survival (Log-rank P < 0.05) in patients from TCGA-LIHC and GSE14520 (Fig. 6; Table 2). Correlations between expression levels of these 12 hub genes and reduced disease-free interval, disease-specific survival, and progression-free interval were also observed in patients from TCGA-LIHC (hazard ratio > 1; Table 6). A circRNA–miRNA–hub gene network consisting of four regulatory modules (hsa_circ_0009910–hsa-miR-455-5p–DLGAP5/MCM5 axis, hsa_circ_0009910–hsa-miR-615-3p–MCM5/MCM6/MCM3/CDC20/CCNB1/CDC7 axis, hsa_circ_0049783–hsa-miR-18a-3p–ZWINT axis, and hsa_circ_0089172–hsa-miR-4524a-3p–CDC7/CCNB1/CENPU/ASPM/ECT2/NDC80 axis) was then built to delineate the links among the identified circRNAs, miRNAs and hub genes (Fig. 6a).
Functional enrichment analysis of hub genes
Enrichment analysis was performed to illustrate the functional annotations of the 12 survival-related hub genes. KEGG pathway enrichment analysis indicated that hub genes were mainly enriched in DNA replication and cell cycle. For the Reactome pathways, the hub genes were enriched in activation of the pre-replicative complex, G1/S transition, resolution of sister chromatid cohesion, and DNA replication. For the Wikipathway, enrichment of the hub genes in regulation of sister chromatid separation at the metaphase-anaphase transition, DNA replication, and retinoblastoma gene in cancer was identified (Fig. 7b-d).
In recent years, the application of large-scale gene expression profiling and the development of bioinformatics has allowed us to systematically uncover the molecular regulation mechanisms associated with tumorigenesis and development. However, large sample sizes are required to identify circRNA expression differences with greater statistical power . In this study, we obtained three up-regulated circRNAs, namely hsa_circ_0009910, hsa_circ_0049783, and hsa_circ_0089172, through jointly analyzing of three tissue microarray datasets and a blood microarray dataset, suggesting these deregulated circRNAs in HCC tissues could be secreted into plasma . Recent studies found that hsa_circ_0009910 was overexpressed in osteosarcoma and ovarian cancer cells, acting as a sponge of miR-449a and promoting the expression of miR-449a target IL6R in osteosarcoma while suppressing miR-145 in ovarian cancer cells [31, 32]. The expression of hsa_circ_0089172 was up-regulated in Hashimoto’s thyroiditis patients and might function via sponging miR-125a-3p . However, the role of hsa_circ_0049783 has not been reported.
To explore the interaction of the candidate circRNAs with miRNA and mRNA, we constructed a circRNA–miRNA–mRNA regulatory network based on the mechanism that circRNAs harboring miRNA binding sites might act as miRNA sponges to regulate downstream gene expression . A total of four miRNAs, namely hsa-miR-455-5p, hsa-miR-615-3p, hsa-miR-18a-3p, and hsa-miR-4524a-3p, were screened out to interact with the candidate circRNAs. Cumulative evidence indicates that hsa-miR-455-5p functions as a tumor suppressor in cancer progression of stomach and prostate [35, 36]. Hsa-miR-615-3p has also been reported as a tumor suppressor in cancers of lung and esophagus . Downstream hub genes, including MCM6, CCNB1, CDC20, NDC80, ZWINT, ASPM, CENPU, MCM3, MCM5, ECT2, CDC7, and DLGAP5, were up-regulated and associated with reduced survival in two HCC cohorts. Our enrichment analysis showed that these hub genes were mainly enriched in DNA replication and cell cycle, indicating their critical roles in the pathogenesis of HCC. Consistently, Jia et al. suggested that MCM6 could be an optimal diagnostic biomarker and a potential therapeutic target for HCC in a Southern Chinese Zhuang population . Ju et al. proved that the expression level of NDC80 was significantly elevated in HCC tissues and contributed to tumor progression . Ying et al. reported that ZWINT was overexpressed in HCC samples and associated with poor overall survival and tumor recurrence . Lin et al. suggested that ASPM overexpression was a marker for enhanced vascular invasive, metastatic potential, and poor prognosis of HCC . Chen et al. identified that the ECT2 upregulation was associated with early HCC recurrence and poor survival . Liao et al. reported that DLGAP5 was significantly up-regulated in HCC specimens and promoted cell proliferation .
In summary, we screened out three novel circRNAs up-regulated both in tissue and plasma from HCC patients and constructed a regulatory mechanism of circRNA–miRNA–mRNA interaction by integrating multiple expression profiles and bioinformatics tools. Our study provides a new perspective on the effectiveness of the candidate circRNAs as biomarkers and clinical treatment targets for HCC.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Villanueva A. Hepatocellular Carcinoma. N Engl J Med. 2019;380(15):1450–62.
Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, Gores G. Hepatocellular carcinoma. Nature Rev Dis Primers. 2016;2:16018.
Wang M, Yu F, Li P. Circular RNAs: Characteristics, Function and Clinical Significance in Hepatocellular Carcinoma. Cancers 2018, 10(8).
Llovet JM, Montal R, Sia D, Finn RS. Molecular therapies and precision medicine for hepatocellular carcinoma. Nature Rev Clin Oncol. 2018;15(10):599–616.
Qiu LP, Wu YH, Yu XF, Tang Q, Chen L, Chen KP. The Emerging Role of Circular RNAs in Hepatocellular Carcinoma. J Cancer. 2018;9(9):1548–59.
Vicens Q, Westhof E. Biogenesis of Circular RNAs. Cell. 2014;159(1):13–4.
Chen LL. The biogenesis and emerging roles of circular RNAs. Nature Rev Molecular Cell Biol. 2016;17(4):205–11.
Ji P, Wu W, Chen S, Zheng Y, Zhou L, Zhang J, Cheng H, Yan J, Zhang S, Yang P, et al. Expanded expression landscape and prioritization of circular RNAs in mammals. Cell Rep. 2019;26(12):3444-60 e3445.
Bahn JH, Zhang Q, Li F, Chan TM, Lin X, Kim Y, Wong DT, Xiao X. The landscape of microRNA, Piwi-interacting RNA, and circular RNA in human saliva. Clinical Chem. 2015;61(1):221–30.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, Zhong G, Yu B, Hu W, Dai L, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Quan G, Li J. Circular RNAs: biogenesis, expression and their potential roles in reproduction. Journal Ovarian Res. 2018;11(1):9.
Zhong Z, Huang M, Lv M, He Y, Duan C, Zhang L, Chen J. Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway. Cancer Lett. 2017;403:305–17.
Chen I, Chen CY, Chuang TJ. Biogenesis, identification, and function of exonic circular RNAs. Wiley Interdis Rev RNA. 2015;6(5):563–79.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, Luo Y, Lyu D, Li Y, Shi G, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nature Commun. 2016;7:11215.
Li X, Yang L, Chen LL. The Biogenesis, Functions, and Challenges of Circular RNAs. Molecular cell . 2018;71(3):428-442.
Zhang Z, Yang T, Xiao J. Circular RNAs: Promising Biomarkers for Human Diseases. EBioMedicine. 2018;34:267–74.
Shang Q, Yang Z, Jia R, Ge S. The novel roles of circRNAs in human cancer. Mol Cancer. 2019;18(1):6.
Bai N, Peng E, Qiu X, Lyu N, Zhang Z, Tao Y, Li X, Wang Z. circFBLIM1 act as a ceRNA to promote hepatocellular cancer progression by sponging miR-346. J Exp Clin Cancer Res. 2018;37(1):172.
Yao Z, Luo J, Hu K, Lin J, Huang H, Wang Q, Zhang P, Xiong Z, He C, Huang Z, et al. ZKSCAN1 gene and its related circular RNA (circZKSCAN1) both inhibit hepatocellular carcinoma cell growth, migration, and invasion but through different signaling pathways. Molecular oncol. 2017;11(4):422–37.
Wan C, Sun Y, Tian Y, Lu L, Dai X, Meng J, Huang J, He Q, Wu B, Zhang Z, et al. Irradiated tumor cell-derived microparticles mediate tumor eradication via cell killing and immune reprogramming. Sci Adv. 2020;6(13):eaay9789.
Xia S, Feng J, Chen K, Ma Y, Gong J, Cai F, Jin Y, Gao Y, Xia L, Chang H, et al. CSCD: a database for cancer-specific circular RNAs. Nucleic Acids Res. 2018;46(D1):D925-9.
Yang Y, Chen L, Gu J, Zhang H, Yuan J, Lian Q, Lv G, Wang S, Wu Y, Yang YT, et al. Recurrently deregulated lncRNAs in hepatocellular carcinoma. Nature Commun. 2017;8:14421.
Sticht C, De La Torre C, Parveen A, Gretz N. miRWalk: An online resource for prediction of microRNA binding sites. PloS ONE. 2018;13(10):e0206239.
Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Current Protocols Bioinform. 2014;47:8.13.11-24.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003;4:2.
Li Q, Pan X, Zhu D, Deng Z, Jiang R, Wang X. Circular RNA MAT2B promotes glycolysis and malignancy of hepatocellular carcinoma through the miR-338-3p/PKM2 Axis under hypoxic stress. Hepatology. 2019;70(4):1298–316.
Xiong DD, Dang YW, Lin P, Wen DY, He RQ, Luo DZ, Feng ZB, Chen G. A circRNA-miRNA-mRNA network identification for exploring underlying pathogenesis and therapy strategy of hepatocellular carcinoma. J Translational Med. 2018;16(1):220.
Yu J, Ding WB, Wang MC, Guo XG, Xu J, Xu QG, Yang Y, Sun SH, Liu JF, Qin LX, et al. Plasma circular RNA panel to diagnose hepatitis B virus-related hepatocellular carcinoma: A large-scale, multicenter study. Int J Cancer. 2020;146(6):1754–63.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, Sun SH, Yang F, Zhou WP. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68(6):1214–27.
Tan TZ, Rouanne M, Tan KT, Huang RY, Thiery JP. Molecular Subtypes of Urothelial Bladder Cancer: Results from a Meta-cohort Analysis of 2411 Tumors. European Urol. 2019;75(3):423–32.
Deng N, Li L, Gao J, Zhou J, Wang Y, Wang C, Liu Y. Hsa_circ_0009910 promotes carcinogenesis by promoting the expression of miR-449a target IL6R in osteosarcoma. Biochem Biophys Res Commun. 2018;495(1):189–96.
Li Y, Lin S, An N. Hsa_circ_0009910: oncogenic circular RNA targets microRNA-145 in ovarian cancer cells. Cell Cycle. 2020;19(15):1857–68.
Xiong S, Peng H, Ding X, Wang X, Wang L, Wu C, Wang S, Xu H, Liu Y. Circular RNA expression profiling and the potential role of hsa_circ_0089172 in Hashimoto’s thyroiditis via sponging miR125a-3p. Molecular Therapy Nucleic Acids. 2019;17:38–48.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Xing Q, Xie H, Zhu B, Sun Z, Huang Y. MiR-455-5p suppresses the progression of prostate cancer by targeting CCR5. BioMed Research Int. 2019;2019:6394784.
Liu J, Zhang J, Li Y, Wang L, Sui B, Dai D. MiR-455-5p acts as a novel tumor suppressor in gastric cancer by down-regulating RAB18. Gene. 2016;592(2):308–15.
Godinez-Rubi M, Ortuno-Sahagun D. miR-615 Fine-Tunes Growth and Development and Has a Role in Cancer and in Neural Repair. Cells 2020, 9(7).
Jia W, Xie L, Wang X, Zhang Q, Wei B, Li H, Qin S, Chen S, Liu J, Tan Y, et al. The impact of MCM6 on hepatocellular carcinoma in a Southern Chinese Zhuang population. Biomedicine Pharmacotherapy 2020;127:110171.
Ju LL, Chen L, Li JH, Wang YF, Lu RJ, Bian ZL, Shao JG. Effect of NDC80 in human hepatocellular carcinoma. World J Gastroenterol. 2017;23(20):3675–83.
Ying H, Xu Z, Chen M, Zhou S, Liang X, Cai X. Overexpression of Zwint predicts poor prognosis and promotes the proliferation of hepatocellular carcinoma by regulating cell-cycle-related proteins. Onco Targets Therapy. 2018;11:689–702.
Lin SY, Pan HW, Liu SH, Jeng YM, Hu FC, Peng SY, Lai PL, Hsu HC. ASPM is a novel marker for vascular invasion, early recurrence, and poor prognosis of hepatocellular carcinoma. Clinical Cancer Res. 2008;14(15):4814–20.
Chen J, Xia H, Zhang X, Karthik S, Pratap SV, Ooi LL, Hong W, Hui KM. ECT2 regulates the Rho/ERK signalling axis to promote early recurrence in human hepatocellular carcinoma. J Hepatol. 2015;62(6):1287–95.
Liao W, Liu W, Yuan Q, Liu X, Ou Y, He S, Yuan S, Qin L, Chen Q, Nong K, et al. Silencing of DLGAP5 by siRNA significantly inhibits the proliferation and invasion of hepatocellular carcinoma cells. PloS ONE. 2013;8(12):e80789.
This work was supported by the National Natural Science Foundation of China (No. 82060513, 82060297, 82060579, and 81902850), the Youth Science and Technology Innovation Leading Talents Project of Corps (2020CB015), and the Youth Innovation Talents Project of Shihezi University (CXBJ201907).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, L., Zhou, L., Hou, J. et al. Three novel circRNAs upregulated in tissue and plasma from hepatocellular carcinoma patients and their regulatory network. Cancer Cell Int 21, 72 (2021). https://doi.org/10.1186/s12935-021-01762-w
- Liver cancer