Identification and validation of immune-related genes in osteoarthritic synovial fibroblasts

OA was generally considered as a non-inflammatory disease dominated by articular cartilage degeneration. However, the role of synovitis in OA pathogenesis has received increasing attention. Recent studies support that OA patients have a pro-inflammatory/catabolic synovial environment similar to RA patients, promoting the occurrence and development of OA. Therefore, we investigated the co-immune-related genes and pathways of OA and RA to explore whether part of the pathogenesis of RA synovitis can be used to explain OA synovitis.

Data of GSE29746 and GSE12021 were downloaded from the Gene Expression Omnibus (GEO) database. Compared with control group, differentially expressed genes (DEGs) of OA and RA groups were screened separately by R software, Venny website was used to screen co-DEGs. Metascape was used to screen the common enriched terms and pathways between OA and RA. STRING website and Cytoscape software were used to map protein-protein interaction (PPI) networks and screen co-hub genes. GSE29746 was selected as the test dataset, and GSE12021 as the validation dataset for validate the co-hub genes. The results were validated by western blotting (WB) and real-time quantitative polymerase chain reaction (qPCR) of clinical synovial samples.

We identified 573 OA-related DEGs, 148 RA-related DEGs, and 52 co-DEGs, revealing 14 common enriched terms, most of which were related to immune inflammation. IL7R was the only upregulated co-hub gene between OA and RA in the PPI network, consistent with the validation dataset. IL7R was highly expressed in clinical osteoarthritic synovial samples (P < 0.001).

Our findings suggested that IL7R is a critical co-DEG in OA and RA and confirmed the involvement of immune inflammation in disease pathogenesis. Furthermore, it confirms the role of IL7R in synovial inflammation in RA and OA synovitis and provides evidence for further investigation of OA immune inflammation.