Circular RNA expression profile of articular chondrocytes in an IL-1β-induced mouse model of osteoarthritis

Abstract

Osteoarthritis (OA) is a widely prevalent degenerative joint disease characterized by articular cartilage degradation and joint inflammation. The pathogenesis of OA remains unclear, leading to a lack of effective treatment. Previous studies have reported that circular RNAs (circRNAs) are involved in the development of various diseases. However, the function of circRNAs and their roles in OA is largely unknown. Therefore, we aimed to investigate changes in circRNA expression and predict their functions in OA by using bioinformatics analysis. An OA model was established in mouse articular chondrocytes (MACs) treated by interleukin-1β (IL-1β), and then the circRNA profile was screened by Next Generation Sequencing. By comparing circRNA expression in IL-1β- treated MACs and normal controls, differentially expressed circRNAs were identified during OA pathogenesis, and differential expression levels of selected circRNAs were validated by qRT-PCR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were employed to predict the functions of these circRNAs. Because circRNAs can act as “miRNA sponges”, we also constructed a circRNA-miRNA network to predict their functions. A total of 255 circRNAs were found to be differentially expressed in IL-1β-treated MACs (p≤0.05; fold-change≥2) from the expression of the normal controls. Among them, 119 circRNAs were significantly up-regulated, and the other 136 were down-regulated. Seven circRNAs were randomly selected to verify the reliability of these profiles by quantitative qRT-PCR. After obtaining the parental genes of differentially expressed circRNA, the top 30 enrichment GO entries and KEGG pathways were annotated. Then, two significantly differentially expressed circRNAs (mmu-circRNA-30365 and mmu-circRNA-36866) were identified and selected for further analysis, meanwhile a circRNA-miRNA regulation network was created and the top five most likely functional-related target miRNAs of the circRNAs were collected. Although the exact mechanisms and biological functions of these circRNAs in the development of OA need further exploration, our findings do suggest that the differentially expressed circRNAs were involved in the pathogenesis of OA. Thus, our study brings us closer to understanding the pathogenic mechanisms and finding new molecular targets for the clinical treatment of osteoarthritis.