Construction of a diagnostic model for osteoarthritis based on transcriptomic immune-related genes

Osteoarthritis (OA) is a leading cause of disability globally, affecting over 500 million individuals worldwide. However, accurate and early diagnosis of OA is challenging to achieve. Immune-related genes play an essential role in OA development. Therefore, the objective of this study was to develop a diagnostic model for OA based on immune-related genes identified in synovial membrane.

The gene expression profile of OA were downloaded based on four datasets. The significantly differentially expressed genes (DEGs) between OA and control groups were selected. The differential immune cells were analyzed, followed by immune-related DEGs screening. WGCNA was used to screen module genes and these genes were further selected through optimization algorithm. Then, nomogram model was constructed. Chemical drug small molecule related to OA was predicted. Finally, expression levels of several key genes were validated by qRT-PCR through construction of OA rat models.

The total 656 DEGs were obtained. Eight immune cells were significantly differential between two groups, and 317 immune-related DEGs were obtained. WGCNA identified three modules. The genes in modules were significantly involved in 15 pathways, involving in 65 genes. Then 12 DEGs were screened as the final optimal combination of DEGs, such as CEBPB, CXCL1, JUND, GABARAPL2 and PDGFC. The Nomogram model was also constructed. Furthermore, the chemical small molecules, such as acetaminophen, aspirin, and caffeine were predicted. The expression levels of CEBPB, CXCL1, GABARAPL2 and PDGFC were validated in OA rat models.

A diagnostic model based on twelve immune related genes was constructed. These model genes, such as CEBPB, CXCL1, GABARAPL2, and PDGFC, may serve as diagnostic biomarkers and immunotherapeutic targets.