Skip to main content
  • Primary research
  • Open access
  • Published:

Construction of a circRNA-miRNA-mRNA network based on competitive endogenous RNA reveals the function of circRNAs in osteosarcoma

Abstract

Background

Osteosarcoma (OS) is a common primary malignant bone tumour. Growing evidence suggests that circular RNAs (circRNAs) are closely related to the development of tumours. However, the function of circRNAs in OS remains unknown. Here, we aimed to determine the regulatory mechanisms of circRNAs in OS.

Methods

The expression profiles of OS circRNA (GSE96964), microRNA (GSE65071) and mRNA (GSE33382) were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed circRNAs, miRNAs and mRNAs in OS. A ceRNA network was constructed based on circRNA-miRNA pairs and miRNA-mRNA pairs. MRNAs with significant prognostic differences were identified by the TARGET database in the network. Functional and pathway enrichment analyses were performed, and interactions between proteins were predicted using Cytoscape. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to elucidate the possible functions of these differentially expressed circRNAs.

Results

A total of 15 downregulated circRNAs, 136 upregulated miRNAs and 52 downregulated mRNAs were identified in OS. Finally, a circRNA-miRNA-mRNA network was constructed in OS based on 14 circRNAs, 24 miRNAs, and 52 mRNAs. GO and KEGG pathway analyses suggested that the mRNAs in the network may be involved in the pathogenesis and progression of OS. Four mRNAs identified by the TARGET database were significantly associated with OS survival prognosis. A circRNA-miRNA-mRNA subnetwork was constructed based on these four mRNAs.

Conclusion

Our results provide a deeper understanding of the regulatory mechanisms by which circRNAs compete for endogenous RNAs in OS.

Background

Osteosarcoma (OS) is one of the most common bone malignancies [1]. It mainly occurs during adolescence, and its incidence decreases with age [2, 3]. At present, the treatment of OS mainly involves surgery and adjuvant chemotherapy, but the prognosis of patients is still not optimistic [4, 5]. Studies have found that the occurrence and development of OS is closely related to gene mutations and gene expression disorders [6,7,8]. Currently, gene-targeted therapy has been widely studied in various cancers and may be a great potential treatment strategy [9,10,11]. Therefore, studying the molecular mechanism of OS regarding its occurrence and development processes and finding effective molecular targets are crucial for the diagnosis and treatment of OS.

Circular RNAs (circRNAs) are a new class of non-coding RNAs with covalently closed loop structures [12, 13]. Due to the lack of a 5′ cap structure and a 3′ poly (A) tail, circRNAs are resistant to exonuclease and are more stable than linear RNAs [14]. A growing body of research indicates that circRNAs play an important regulatory role in the development of cancer in vivo. Zhou et al. found that by sponging miR-203, circRNA-0008717 inhibits the proliferation and invasion of OS cells [15]. Song et al. found that hsa_circ_0001564 acts as a miR-29c-3p sponge to cause tumourigenesis and may serve as a potential biomarker for OS [16]. In addition, a related study found that circRNA CDR1as/miR-7 signalling promotes tumour growth in OS, with potential therapeutic and diagnostic value [17]. These studies indicate that circRNAs play an important role in the development and progression of OS.

In this study, we collected the expression profiles of circRNAs, miRNAs and mRNAs in OS from the Gene Expression Omnibus (GEO) database and the TARGET database. Differentially expressed circRNAs, miRNAs and mRNAs were identified. CircRNA-miRNA pairs and miRNA-mRNA pairs were screened, and a circRNA-miRNA-mRNA network was constructed. The functional pathways of the mRNAs in this network were analysed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment. To better understand the pathogenesis of OS, TARGET data were used to screen mRNAs related to OS prognosis, and a circRNA-miRNA-mRNA subnetwork related to OS prognosis was constructed. Our research provides new insights into the mechanisms of the occurrence and development of OS. The flowchart for this procedure is shown in Fig. 1.

Fig. 1
figure 1

Flow chart of the approach utilized in the present study

Methods

Data acquisition

GEO (http://www.ncbi.nlm.nih.gov/geo) is a public functional genomic database. Three datasets (GSE96964, GSE65071, and GSE33382; analysed on Agilent-069978 Arraystar Human CircRNA microarray V1, Exiqon human V3 microRNA PCR panel I + II, and Illumina human-6 v2.0 expression beadchip, respectively) were downloaded from GEO. The probes were converted to the corresponding gene symbols based on the annotation information on the platform. The GSE96964 dataset contained a human osteoblast cell line (hFOB1.19) and 7 human osteosarcoma cell lines (U2OS, MTX300, HOS, MG63,143B, ZOS, and ZOSM). The GSE65071 dataset contained 15 healthy people and 20 patients with osteosarcoma. The GSE33382 dataset contained 3 human osteoblasts and 83 osteosarcoma samples. The TARGET database (https://ocg.cancer.gov/programs/target) is a database of childhood tumours that seeks to determine molecular changes in the development and progression of refractory childhood cancer using a comprehensive genomic approach. This database contains gene expression profiles and clinical information for patients with osteosarcomas. We obtained 101 osteosarcoma samples with complete clinical information from the TARGET database. Approval by the Ethics Committee was not necessary because all data were collected from publicly available databases (GEO and TARGET).

Identification of differential expression of circRNAs, miRNAs and mRNAs

The raw data were processed by background correction and normalization by using the affy package of R/Bioconductor. The Limma and edgeR packages were used to identify differentially expressed circRNAs, miRNAs and mRNAs between normal samples and tumour samples. Log [fold change (FC)] > 2 and adj. P value < 0.05 were considered to indicate a statistically significant difference.

Construction of the ceRNA network

Information about circRNAs is available on the circBase website (http://www.circbase.org/). The cancer-specific circRNA database (CSCD) (https://gb.whu.edu.cn/CSCD/) was used to predict target miRNAs. These target miRNAs were further screened by miRNAs differentially expressed in GSE65071. Then, the miRDB, TargetScan and miRTarBase databases were employed to predict specific miRNAs. Only mRNAs recognized by the 3 databases were considered to be candidate targets and were intersected with the identified differentially expressed mRNAs to identify the differentially expressed mRNAs targeted by the differentially expressed miRNAs. Finally, the circRNA-miRNA-mRNA regulatory network was constructed by screening the circRNA-miRNA and miRNA-mRNA groups. The data were visualized using Cytoscape 3.7.1.

GO and KEGG functional enrichment analysis

Database for Annotation, Visualization and Integrated Discovery (DAVID; version 6.7; http://david.ncifcrf.gov) is an online bioinformatics database that integrates biological data and analysis tools and provides gene annotation information and protein data. KEGG is a database resource for understanding advanced functions and biological systems from large-scale molecular data generated by high-throughput experimental techniques. GO is a major bioinformatics tool for annotating genes and analysing their biological processes (BPs), molecular functions (MFs) and cellular components (CCs). To determine the function of differentially expressed mRNAs in the network, analysis was performed using DAVID. An adj. P-value < 0.05 was considered to indicate a statistically significant difference.

Survival analysis

To analyse the mRNAs selected from the ceRNA network, the OS mRNA expression profile data and clinical sample data in the TARGET database were downloaded. The overall survival analysis was performed using Kaplan–Meier curves. The log-rank test was used for statistical analysis. The threshold for survival prognosis significance was a P-value < 0.05. A subnetwork was constructed based on these mRNAs.

Results

Identification of differentially expressed circRNAs, miRNAs and mRNAs

The basic information regarding the expression of circRNAs, miRNAs and mRNAs in OS and control tissues in the three microarray datasets (GSE96964, GSE65071 and GSE33382) is listed in Table 1. In the circRNA expression profile data, a total of 15 circRNAs were downregulated in OS (Fig. 2, Additional file 1: Figure S1). The basic information for the 15 circRNAs is listed in Table 2. One of the circRNAs (hsa_circ_0000144) was not found in the CSCD database. From the CSCD database, we predicted that 503 miRNAs might target these 14 circRNAs (Additional file 2: Table S1). The basic structural pattern of the 14 circRNAs is shown in Fig. 3. In the miRNA expression profile data, a total of 136 miRNAs were upregulated in OS (Additional file 3: Table S2). Twenty-four intersecting miRNAs were obtained and are shown in Additional file 4: Table S3 (Fig. 4a). A total of 324 potential target genes (Additional file 5: Table S4) were predicted by the TargetScan, miRTarBase and miRDB databases. A total of 52 mRNAs were downregulated in OS in the mRNA expression profile data (Additional file 6: Table S5). There were 52 intersecting mRNAs (Fig. 4b), as shown in Additional file 7: Table S6.

Table 1 Basic information of the 4 microarray datasets from GEO and TARGET
Fig. 2
figure 2

Heatmap of the fifteen differentially expressed circRNAs in the eight microarray datasets. By analysing the circRNA dataset, fifteen differentially expressed circRNAs were identified (|log 2 (fold change (FC))| > 2 and an adjusted P value < 0.05 were considered)

Table 2 Basic characteristics of the fifteen differentially expressed circRNAs
Fig. 3
figure 3

Structural patterns of the fourteen circRNAs. an A total of 14 circRNA structures were obtained from the circRNA website. Red represents the position where the miRNA may bind, blue represents the position where the protein may bind, and yellow represents the open reading frame

Fig. 4
figure 4

The differentially expressed miRNAs and mRNAs in the microarray datasets. a Twenty-four miRNAs were obtained from the intersection of miRNAs predicted from the website and the differentially expressed miRNAs identified in the miRNA dataset (GEO dataset). The expression levels of 24 miRNAs in the miRNA dataset are shown in a. b A total of 52 differentially expressed mRNAs were obtained by analysing the mRNA dataset and predicting the targeted binding of genes by the 24 miRNAs. The expression levels of these 52 mRNAs in the mRNA dataset are shown in b. (|log 2 (fold change (FC))| > 2, and an adjusted P value < 0.05 were the cut-off criteria). *P < 0.05, **P < 0.01, ***P < 0.001

Construction of the ceRNA network

We used 14 circRNAs, 24 miRNAs and 52 mRNAs in Cytoscape 3.7.1 to construct a circRNA-miRNA-mRNA visualization network (Fig. 5).

Fig. 5
figure 5

CeRNA network of circRNA-miRNA-mRNA interactions in OS. Red indicates circRNAs, green indicates miRNAs, and blue indicates mRNAs

Functional and pathway enrichment analysis

The results indicated that the identified mRNAs were mainly enriched in ‘response to extracellular stimulus’, ‘response to oxygen levels’, ‘response to nutrient levels’ and ‘regulation of protein serine/threonine kinase activity’ (BPs) (Fig. 6a); ‘protein kinase complex’, ‘transferase complex, transferring phosphorus-containing groups’ and ‘neuronal cell body’ (CCs) (Fig. 6b); and ‘ubiquitin protein ligase binding’ and ‘ubiquitin-like protein ligase binding’ (MFs) (Fig. 6c). KEGG pathway analysis revealed strong enrichment in the ‘PI3K-Akt signaling pathway’, ‘Proteoglycans in cancer’, ‘Cell cycle’ and ‘FoxO signaling pathway’ (Fig. 6d).

Fig. 6
figure 6

GO and KEGG functional enrichment analyses of the differentially expressed mRNAs in the ceRNA network. a The enrichment of biological processes. b The enrichment of cellular components. c The enrichment of molecular functions. d The enrichment analysis of the KEGG pathways. A P < 0.05 was considered to indicate a statistically significant difference

Survival analysis and subnetwork construction of related mRNAs

We downloaded OS mRNA expression profile data (101 samples) and clinical sample data from the TARGET database. We performed survival analysis for the 52 mRNAs in the constructed circRNA-miRNA-mRNA network. Four mRNAs (ARID5B, ELL2, PHC2, and STAT4) were significantly associated with survival prognosis in patients with OS (Fig. 7). We constructed a circRNA-miRNA-mRNA subnetwork with these four mRNAs (Fig. 8).

Fig. 7
figure 7

Survival analysis of the differentially expressed mRNAs in the ceRNA network. The TARGET database was used to analyse the survival prognosis of 52 mRNAs in osteosarcoma patients. a The ARID5B gene high expression group in osteosarcoma samples had a good prognosis compared with the low expression group (P = 0.041). b The ELL2 gene high expression group had a good prognosis compared with the low expression group (P = 0.011). c The PHC2 gene high expression group had a good prognosis compared with the low expression group (P = 0.006). d The STAT4 gene high expression group had a good prognosis compared with the low expression group (P = 0.007). (P value < 0.05 was the cut-off criterion)

Fig. 8
figure 8

Construction of the circRNA-miRNA-mRNA subnetwork based on mRNAs with survival prognostic potential. Red indicates circRNAs, green indicates miRNAs, and blue indicates mRNAs

Discussion

CircRNAs are a class of circular noncoding RNAs that have been observed in a variety of cancers [18,19,20], but their function is still unclear. Due to their special circular structure, circRNAs have high stability and great potential as new tumour markers [21]. Recent studies have found that circRNAs are involved in the development and progression of various cancers [22,23,24], such as gastric cancer, breast cancer, and colorectal cancer. Current evidence suggests that circRNAs can target miRNAs, often referred to as “miRNA sponges,” to reduce the levels of miRNAs and release their targeted inhibition of mRNAs, thereby regulating the expression of protein-coding genes [25,26,27]. These studies have found that circRNAs may play a key role in the development and progression of cancers.

OS is a common bone malignancy [1]. Current treatments for OS are mainly surgery and chemotherapy, but patient prognosis remains unsatisfactory [4]. Due to their stable structure, circRNAs may be a new strategy in the targeted therapy of OS [21]. A related study found that the circRNA LRP6 negatively regulates the levels of KLF2 and APC to promote the development of OS [28]. Dysregulated circRNA_100876 suppresses the proliferation of OS cancer cells by targeting microRNA-136 [29]. The upregulation of the circular RNA circ_0001721 predicts unfavourable prognosis in OS and facilitates cell cycle progression by sponging miR-569 and miR-599 [30]. Although these studies have found that circRNAs are dysregulated in OS and are involved in the development of cancer, the specific mechanism of circRNAs in OS remains unclear.

To investigate the role of circRNAs in OS, a circRNA-miRNA-mRNA regulatory network was constructed based on bioinformatics predictions and transcriptome data. GO enrichment analysis showed that genes in this network were mainly involved in the regulation of protein serine/threonine kinase activity, the transferase complex and ubiquitin protein ligase binding. Many studies have found that serine/threonine kinases play an important role in the occurrence and development of diseases. Insulin-like growth factor I (IGF-1) induces BNIP3 expression through a known AKT serine/threonine kinase 1 (AKT1)-mediated inhibitory phosphorylation of glycogen synthase kinase 3β (GSK-3β), which promotes cell growth and tumourigenesis [31]. Cooperative recruitment of Arl4A and Pak1 to the plasma membrane contributes to sustained Pak1 activation for cell migration [32]. Furthermore, the ubiquitination of proteins has been found to be involved in the cancer process. Suppression of the ubiquitin pathway by small molecule binding to ubiquitin enhances the doxorubicin sensitivity of cancer cells [33]. Epigenetic silencing of ubiquitin-specific protease 4 by snail1 contributes to macrophage-dependent inflammation and therapeutic resistance in lung cancer [34]. KEGG pathway analysis showed that the network was enriched in the PI3K-Akt signalling pathway, cell cycle and FoxO signalling, which have been confirmed by a large number of studies to participate in cancer progression [35,36,37,38]. These results suggest that the circRNAs in this network may play a key role in the occurrence and development of osteosarcoma. This result is worthy of further investigation.

To further study which circRNAs affect the prognosis of OS patients, we analysed the correlation between the expression levels of 52 mRNAs in this network and overall survival in OS patients based on TARGET data. Our study found that four mRNAs (ARID5B, ELL2, PHC2, and STAT4) may affect the prognosis of patients with OS. We constructed a circRNA-miRNA-mRNA subnetwork. There were 8 circRNAs (hsa_circ_0002052, hsa_circ_0088214, hsa_circ_0088227, hsa_circ_0008792, hsa_circ_0061265, hsa_circ_0002017, hsa_circ_00054442 and hsa_circ_0077765) and four miRNAs (miR-33a-3p, miR-205-3141, miR-576-3p) in the network. Studies have found that STAT4 is closely related to the prognosis of a variety of cancers, such as breast, liver, and stomach cancers [39,40,41]. By the targeting of ELL2 by microRNA-299 in glioblastoma multiforme, the upregulation of microRNA-299 is associated with a poor prognosis [42]. A previous study found that Phc2 controls the mobilization of bone marrow haematopoietic stem cells and progenitor cells by inhibiting Vcam1 expression [43]. ARID5B influences antimetabolite drug sensitivity and prognosis of acute lymphoblastic leukaemia [44]. Although we have identified four prognostic genes that have not been reported in OS, our predictions are consistent with studies in other cancers. Whether the dysregulation of these four genes will affect the prognosis of OS patients requires more in-depth research in the future.

At present, there are few reports on the mechanism of circRNAs in osteosarcoma. The novelty of this study is not only that the circRNA-miRNA-mRNA network was constructed through the GEO database but is also related to the prognosis of OS through the TARGET database. Previous studies identified OS survival prognostic genes based only on the GEO database, which contains few clinical samples of osteosarcoma. The TARGET database contains more osteosarcoma samples with clinical information, which can provide us with richer data for identifying genes related to survival prognosis. However, this study is also limited because it is mainly based on the analysis of sequencing data and may require further experimental exploration in the future.

Conclusion

We identified differentially expressed circRNAs, miRNAs and mRNAs from publicly available microarray data to construct a circRNA-related ceRNA network. In addition, four mRNAs related to prognosis were identified based on the TARGET database. With these four mRNAs, we constructed a circRNA-miRNA-mRNA subnetwork that may be associated with prognosis in OS. The findings of this study may provide new insights into the pathogenesis of OS and suggest potential therapeutic targets that warrant further investigation.

Availability of data and materials

We declare that the data and materials in this study are provided free of charge to scientists for non-commercial purposes.

References

  1. Ritter J, Bielack SS. Osteosarcoma. Ann Oncol. 2010;21(Suppl 7):vii320–5.

    Article  PubMed  Google Scholar 

  2. Geller DS, Gorlick R. Osteosarcoma: a review of diagnosis, management, and treatment strategies. Clin Adv Hematol Oncol. 2010;8(10):705–18.

    PubMed  Google Scholar 

  3. Yang L, Ge D, Chen X, Qiu J, Yin Z, Zheng S, et al. FOXP4-AS1 participates in the development and progression of osteosarcoma by downregulating LATS1 via binding to LSD1 and EZH2. Biochem Biophys Res Commun. 2018;502(4):493–500.

    Article  CAS  PubMed  Google Scholar 

  4. Link MP, Goorin AM, Miser AW, Green AA, Pratt CB, Belasco JB, et al. The effect of adjuvant chemotherapy on relapse-free survival in patients with osteosarcoma of the extremity. N Engl J Med. 1986;314(25):1600–6.

    Article  CAS  PubMed  Google Scholar 

  5. Jaffe N, Frei E 3rd, Traggis D, Bishop Y. Adjuvant methotrexate and citrovorum-factor treatment of osteogenic sarcoma. N Engl J Med. 1974;291(19):994–7.

    Article  CAS  PubMed  Google Scholar 

  6. Tian Z, Yang G, Jiang P, Zhang L, Wang J, Sun S. Long non-coding RNA Sox4 promotes proliferation and migration by activating Wnt/beta-catenin signaling pathway in osteosarcoma. Pharmazie. 2017;72(9):537–42.

    CAS  PubMed  Google Scholar 

  7. Min X, Heng H, Yu HL, Dan M, Jie C, Zeng Y, et al. Anticancer effects of 10-hydroxycamptothecin induce apoptosis of human osteosarcoma through activating caspase-3, p53 and cytochrome c pathways. Oncol Lett. 2018;15(2):2459–64.

    PubMed  Google Scholar 

  8. Ma C, Zhang Z, Cui Y, Yuan H, Wang F. Silencing FAT10 inhibits metastasis of osteosarcoma. Int J Oncol. 2016;49(2):666–74.

    Article  CAS  PubMed  Google Scholar 

  9. Marks EI, Yee NS. Molecular genetics and targeted therapeutics in biliary tract carcinoma. World J Gastroenterol. 2016;22(4):1335–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lin F, Li Z, Hua Y, Lim YP. Proteomic profiling predicts drug response to novel targeted anticancer therapeutics. Expert Rev Proteomics. 2016;13(4):411–20.

    Article  CAS  PubMed  Google Scholar 

  11. Buzdin A, Sorokin M, Garazha A, Sekacheva M, Kim E, Zhukov N, et al. Molecular pathway activation—new type of biomarkers for tumor morphology and personalized selection of target drugs. Semin Cancer Biol. 2018;53:110–24.

    Article  CAS  PubMed  Google Scholar 

  12. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, et al. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365(2):141–8.

    Article  CAS  PubMed  Google Scholar 

  13. Petkovic S, Muller S. RNA circularization strategies in vivo and in vitro. Nucleic Acids Res. 2015;43(4):2454–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Li J, Yang J, Zhou P, Le Y, Zhou C, Wang S, et al. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res. 2015;5(2):472–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhou X, Natino D, Qin Z, Wang D, Tian Z, Cai X, et al. Identification and functional characterization of circRNA-0008717 as an oncogene in osteosarcoma through sponging miR-203. Oncotarget. 2018;9(32):22288–300.

    Article  PubMed  Google Scholar 

  16. Song YZ, Li JF. Circular RNA hsa_circ_0001564 regulates osteosarcoma proliferation and apoptosis by acting miRNA sponge. Biochem Biophys Res Commun. 2018;495(3):2369–75.

    Article  CAS  PubMed  Google Scholar 

  17. Xu B, Yang T, Wang Z, Zhang Y, Liu S, Shen M. CircRNA CDR1as/miR-7 signals promote tumor growth of osteosarcoma with a potential therapeutic and diagnostic value. Cancer Manag Res. 2018;10:4871–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Chen BJ, Huang S, Janitz M. Changes in circular RNA expression patterns during human foetal brain development. Genomics. 2019;111(4):753–8.

    Article  CAS  PubMed  Google Scholar 

  19. Patop IL, Kadener S. circRNAs in Cancer. Curr Opin Genet Dev. 2018;48:121–7.

    Article  CAS  PubMed  Google Scholar 

  20. Yu T, Wang Y, Fan Y, Fang N, Wang T, Xu T, et al. CircRNAs in cancer metabolism: a review. J Hematol Oncol. 2019;12(1):90.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Jiang XM, Li ZL, Li JL, Xu Y, Leng KM, Cui YF, et al. A novel prognostic biomarker for cholangiocarcinoma: circRNA Cdr1as. Eur Rev Med Pharmacol Sci. 2018;22(2):365–71.

    PubMed  Google Scholar 

  22. Jiang F, Shen X. Current prevalence status of gastric cancer and recent studies on the roles of circular RNAs and methods used to investigate circular RNAs. Cell Mol Biol Lett. 2019;24:53.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Li Z, Chen Z, Hu G, Jiang Y. Roles of circular RNA in breast cancer: present and future. Am J Transl Res. 2019;11(7):3945–54.

    PubMed  PubMed Central  CAS  Google Scholar 

  24. Hao S, Cong L, Qu R, Liu R, Zhang G, Li Y. Emerging roles of circular RNAs in colorectal cancer. Onco Targets Ther. 2019;12:4765–77.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Chen J, Chen T, Zhu Y, Li Y, Zhang Y, Wang Y, et al. circPTN sponges miR-145-5p/miR-330-5p to promote proliferation and stemness in glioma. J Exp Clin Cancer Res. 2019;38(1):398.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Jin C, Shi L, Li Z, Liu W, Zhao B, Qiu Y, et al. Circ_0039569 promotes renal cell carcinoma growth and metastasis by regulating miR-34a-5p/CCL22. Am J Transl Res. 2019;11(8):4935–45.

    PubMed  PubMed Central  CAS  Google Scholar 

  27. Ren S, Liu J, Feng Y, Li Z, He L, Li L, et al. Knockdown of circDENND4C inhibits glycolysis, migration and invasion by up-regulating miR-200b/c in breast cancer under hypoxia. J Exp Clin Cancer Res. 2019;38(1):388.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Zheng S, Qian Z, Jiang F, Ge D, Tang J, Chen H, et al. CircRNA LRP6 promotes the development of osteosarcoma via negatively regulating KLF2 and APC levels. Am J Transl Res. 2019;11(7):4126–38.

    PubMed  PubMed Central  CAS  Google Scholar 

  29. Jin J, Chen A, Qiu W, Chen Y, Li Q, Zhou X, et al. Dysregulated circRNA_100876 suppresses proliferation of osteosarcoma cancer cells by targeting microRNA-136. J Cell Biochem. 2019;120(9):15678–87.

    Article  CAS  PubMed  Google Scholar 

  30. Li L, Guo L, Yin G, Yu G, Zhao Y, Pan Y. Upregulation of circular RNA circ_0001721 predicts unfavorable prognosis in osteosarcoma and facilitates cell progression via sponging miR-569 and miR-599. Biomed Pharmacother. 2019;109:226–32.

    Article  CAS  PubMed  Google Scholar 

  31. Riis S, Murray JB, O’Connor R. IGF-1 Signalling regulates mitochondria dynamics and turnover through a conserved GSK-3beta-Nrf2-BNIP3 pathway. Cells. 2020;9(1):147.

    Article  CAS  PubMed Central  Google Scholar 

  32. Chen KJ, Chiang TC, Yu CJ, Lee FS. Arl4A and Pak1 cooperative recruitment to plasma membrane contributes to sustained Pak1 activation for cell migration. J Cell Sci. 2020.

  33. Nguyen T, Ho M, Kim K, Yun SI, Mizar P, Easton JW, et al. Suppression of the ubiquitin pathway by small molecule binding to ubiquitin enhances doxorubicin sensitivity of the cancer cells. Molecules. 2019;24(6):1073.

    Article  CAS  PubMed Central  Google Scholar 

  34. Lai CY, Yeh DW, Lu CH, Liu YL, Chuang YC, Ruan JW, et al. Epigenetic silencing of ubiquitin specific protease 4 by Snail1 contributes to macrophage-dependent inflammation and therapeutic resistance in lung cancer. Cancers. 2020;12(1):148.

    Article  CAS  PubMed Central  Google Scholar 

  35. Lin Y, Chen Z, Zheng Y, Liu Y, Gao J, Lin S, et al. MiR-506 targets UHRF1 to inhibit colorectal cancer proliferation and invasion via the KISS1/PI3K/NF-kappaB signaling axis. Front Cell Dev Biol. 2019;7:266.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Song M, Wang X, Luo Y, Liu Z, Tan W, Ye P, et al. Cantharidin suppresses gastric cancer cell migration/invasion by inhibiting the PI3K/Akt signaling pathway via CCAT1. Chem Biol Interact. 2020;317:108939.

    Article  CAS  Google Scholar 

  37. Qian Y, Wu X, Wang H, Hou G, Han X, Song W. PAK1 silencing is synthetic lethal with CDK4/6 inhibition in gastric cancer cells via regulating PDK1 expression. Human Cell. 2020. https://doi.org/10.1007/s13577-019-00317-6

    Article  PubMed  Google Scholar 

  38. Yan Y, Huang H. Interplay among PI3K/AKT, PTEN/FOXO and AR signaling in prostate cancer. Adv Exp Med Biol. 2019;1210:319–31.

    Article  PubMed  CAS  Google Scholar 

  39. Wang S, Yu L, Shi W, Li X, Yu L. Prognostic roles of signal transducers and activators of transcription family in human breast cancer. Biosci Rep. 2018;38(6).

  40. Nishi M, Batsaikhan BE, Yoshikawa K, Higashijima J, Tokunaga T, Takasu C, et al. High STAT4 expression indicates better disease-free survival in patients with gastric cancer. Anticancer Res. 2017;37(12):6723–9.

    CAS  PubMed  Google Scholar 

  41. Wang G, Chen JH, Qiang Y, Wang DZ, Chen Z. Decreased STAT4 indicates poor prognosis and enhanced cell proliferation in hepatocellular carcinoma. World J Gastroenterol. 2015;21(13):3983–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang Q, Zhang XW, Ma YS, Lu GX, Xie RT, Yang HQ, et al. Up-regulated microRNA-299 corrected with poor prognosis of glioblastoma multiforme patients by targeting ELL2. Jpn J Clin Oncol. 2017;47(7):590–6.

    Article  PubMed  Google Scholar 

  43. Bae J, Choi SP, Isono K, Lee JY, Park SW, Choi CY, et al. Phc2 controls hematopoietic stem and progenitor cell mobilization from bone marrow by repressing Vcam1 expression. Nat Commun. 2019;10(1):3496.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Xu H, Zhao X, Bhojwani D, Shuyu E, Goodings C, Zhang H, et al. ARID5B influences antimetabolite drug sensitivity and prognosis of acute lymphoblastic leukemia. Clin Cancer Res. 2020;26(1):256–64.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank the reviewers for their constructive comments.

Funding

There is no funding to report.

Author information

Authors and Affiliations

Authors

Contributions

The contributions of all authors in the study are equivalent. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Baochuang Qi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors consent to publication.

Competing interests

The authors declare that they have no competing interests in this work.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

The expression levels of fifteen circRNAs in osteosarcoma.

Additional file 2: Table S1.

Five hundred and three miRNAs predicted to target these 14 circRNAs.

Additional file 3: Table S2.

A total of 136 miRNAs were upregulated in OS.

Additional file 4: Table S3.

Twenty-four intersecting miRNAs.

Additional file 5: Table S4.

A total of 324 potential target genes.

Additional file 6: Table S5.

A total of 52 mRNAs were downregulated in OS in the mRNA expression profile data.

Additional file 7: Table S6.

Fifty-two intersecting mRNAs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qiu, Y., Pu, C., Li, Y. et al. Construction of a circRNA-miRNA-mRNA network based on competitive endogenous RNA reveals the function of circRNAs in osteosarcoma. Cancer Cell Int 20, 48 (2020). https://doi.org/10.1186/s12935-020-1134-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-020-1134-1

Keywords