Systematic Analysis of Transcriptomic Profile of Chondrocytes in Osteoarthritic Knee Using Next-Generation Sequencing and Bioinformatics

The phenotypic change of chondrocytes and the interplay between cartilage and subchondral bone in osteoarthritis (OA) has received much attention. Structural changes with nerve ingrowth and vascular penetration within OA cartilage may contribute to arthritic joint pain. The aim of this study was to identify differentially expressed genes and potential miRNA regulations in OA knee chondrocytes through next-generation sequencing and bioinformatics analysis. Results suggested the involvement of SMAD family member 3 (SMAD3) and Wnt family member 5A (WNT5A) in the growth of blood vessels and cell aggregation, representing features of cartilage damage in OA. Additionally, 26 dysregulated genes with potential miRNA-mRNA interactions were identified in OA knee chondrocytes. Myristoylated alanine rich protein kinase C substrate (MARCKS), epiregulin (EREG), leucine rich repeat containing 15 (LRRC15), and phosphodiesterase 3A (PDE3A) expression patterns were similar among related OA cartilage, subchondral bone and synovial tissue arrays in Gene Expression Omnibus database. The Ingenuity Pathway Analysis identified MARCKS to be associated with the outgrowth of neurite, and novel miRNA regulations were proposed to play critical roles in the pathogenesis of the altered OA knee joint microenvironment. The current findings suggest new perspectives in studying novel genes potentially contributing to arthritic joint pain in knee OA, which may assist in finding new targets for OA treatment.