Identification and validation of FPR1, FPR2, IL17RA and TLR7 as immunogenic cell death related genes in osteoarthritis

Immunogenic cell death (ICDs) has gained increasing attention for its significant clinical efficacy in various diseases. Similarly, more and more attention has been paid in the role of immune factors in the pathological process of osteoarthritis (OA). The objective of this study is to reveal the relationship between ICD-related genes and the process of OA at the gene level through bioinformatics analysis. In this study, Limma R package was applied to identify differentially expressed genes (DEG), and OA related module genes were determined by weighted gene co-expression network analysis. The ICD-related genes were extracted from a previous study. The module genes related to DEGs and ICD were overlapped. Then, hub genes were identified by a series of analyses using the Least absolute shrinkage and selection operator and random forest algorithm, the expression level and diagnostic value of hub genes were evaluated by Logistic regression. In addition, we used Spearman rank correlation analysis to clarify the relationship between hub genes and infiltrating immune cells and immune pathways. The expression levels of FPR1, FPR2, IL17RA, and TLR7 was verified in SD rat knee joint model of OA by immunohistochemistry. The expression levels of FPR1, FPR2, IL17RA, and TLR7 mRNA were detected in the IL-1β induced rat chondrocytes in qPCR experiment in vitro. Four hub genes (FPR1, FPR2, IL17RA, and TLR7) were ultimately identified as OA biomarkers associated with ICD. And knockdown of TLR7 reversed collagen II and ADAMTS-5 degradation in IL-1β-stimulated chondrocytes. This research may provide new immune related biomarkers for the diagnosis of OA and serve as a reference for disease treatment monitoring.