Microarray Expression Profile and Functional Analysis of Circular RNAs in Osteosarcoma

Abstract

Background/Aims: Osteosarcoma (OS) is the most common primary malignant bone tumor in children and adolescents. However, the molecular mechanisms regulating osteosarcoma tumorigenesis and progression are still poorly understood. Circular RNAs (circRNAs) have been identified as microRNA sponges and are involved in many important biological processes. This study aims to investigate the global changes in the expression pattern of circRNAs in osteosarcoma and provide a comprehensive understanding of differentially expressed circRNAs. Methods: Microarray based circRNA expression was determined in osteosarcoma cell lines and compared with hFOB1.19, which was used as the normal control. We confirmed the microarray data by real time-qPCR in both osteosarcoma cell lines and tissues. The circRNA/microRNA/mRNA interaction network was predicted using bioinformatics. Gene Ontology analysis and 4 annotation tools for pathway analysis (KEGG, Biocarta, PANTHER and Reactome) were used to predict the functions of differentially expressed circRNAs. Results: We revealed a number of differentially expressed circRNAs and 12 of them were confirmed, which suggests a potential role of circRNAs in OS. Among these differentially expressed circRNAs, hsa_circRNA_103801 was up-regulated in both osteosarcoma cell lines and tissues, while hsa_circRNA_104980 was down-regulated. The most likely potential target miRNAs for hsa_circRNA_103801 include hsa-miR-370-3p, hsa-miR-338-3p and hsa-miR-877-3p, while the most potential target miRNAs of hsa_circRNA_104980 consist of hsa-miR-1298-3p and hsa-miR-660-3p. Functional analysis found that hsa_circRNA_103801 was involved in pathways in cancer, such as the HIF-1, VEGF and angiogenesis pathway, the Rap1 signaling pathway and the PI3K-Akt signaling pathway, while hsa_circRNA_104980 was related to some pathways such as the tight junction pathway. Conclusions: This study has identified the comprehensive expression profile of circRNAs in osteosarcoma for the first time. And the ceRNA network prediction and bioinformatics functional analysis could provide a comprehensive understanding of hsa_circRNA_103801 and hsa_circRNA_104980, which may be involved in the initiation and progression of osteosarcoma. The present study indicates that circRNAs may play important roles in osteosarcoma and thus serve as biomarkers of osteosarcoma diagnosis and treatment.