Microarray profiling of circular RNAs in steroid-associated osteonecrosis of the femoral head: Observational study

The aim of this study was to elucidate the molecular mechanisms and to identify the differential expression of circular RNAs (circRNAs) for steroid-associated osteonecrosis of the femoral head (SONFH) using bioinformatics analysis. circRNA microarray was performed with 3 SONFH tissues and the adjacent normal tissues, and differentially expressed circRNA were identified by limma package in R. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery database. In addition, a differentially expressed genes (DEG)-associated circRNA/microRNA (miRNA) interaction was predicted by combination of TargetScan and miRanda, and the circRNA/miRNA interaction network generated by the cytoscape software. A total of 647 differentially expressed circRNAs, including 433 upregulated and 214 downregulated circRNA were identified. The most enriched GO terms for upregulated and downregulated circRNA were extracellular matrix organization and leukocyte activation in biological process; extracellular matrix and spindle pole in cellular component; integrin binding and ATP binding in molecular function, and KEGG pathway enrichment analyses showed that the upregulated and downregulated circRNA were strongly associated with Protein digestion and absorption and Cell cycle. Moreover, a total of 212 differentially expressed messenger RNAs (mRNAs), including 113 upregulated and 99 downregulated genes were identified. In addition, from the analysis of miRNA, long noncoding RNAs, mRNA, and circRNA networks, we found that hsa_circ_0008136 and hsa_circ_0074758 were respectively the upregulated and downregulated circRNA with highest degrees. The identified circRNA and mRNA could be implicated in the progression of human SONFH. The findings could lead to a better understanding of SONFH pathogenesis.