Diagnostic biomarkers and potential treatment targets for Ankylosing Spondylitis

Abstract Background Ankylosing spondylitis (AS) is a chronic autoimmune disease that primarily affects the axial joints. Immune cells play a key role in the pathogenesis of AS. This study integrated bioinformatics methods with experimental validation to explore the role of natural killer (NK) cells in AS. Methods Two microarray datasets, GSE25101 and GSE73754, were selected, and the scRNA-seq data were obtained from GSE194315 and Liu’s research. Differentially expressed genes (DEGs) and functional enrichment analysis were performed respectively. weighted gene co-expression network analysis (WGCNA) was conducted to identify key modules of co-expressed genes and genes involved in NK cell function. The diagnostic value of the identified key genes were evaluated using ROC curves, logistic regression analysis, and a nomogram. real-time PCR (RT-PCR) was used to quantified the expression of genes. Statistical analysis was conducted using the R software package, and a P-value of less than 0.05 was considered statistically significant. Results Pathways enrichment analysis revealed the involvement of NK cell-mediated immune pathways and regulation of the innate immune response, indicating the crucial role of innate immunity, especially NK cells, in AS pathogenesis. The construction of a co-expression network revealed that the MElightyellow module was most relevant to the NK cell-mediated immune pathway. IL2RB, CD247, PLEKHF1, EOMES, S1PR5, FGFBP2 from the MElightyellow module were identified as key genes involved in NK cell-mediated immune response and served as potential diagnostic biomarkers for AS, with moderate to high diagnostic values based on AUC values. Further analysis using scRNA-seq profiling revealed the higher expression level of IL2RB, CD247, PLEKHF1, S1PR5, FGFBP2 in NK cells compared to that in other cell types. CD247, PLEKHF1, EOMES, S1PR5, and FGFBP2 were reduced expressed in AS patients as compare to control group verified by scRNA-seq data, CD247, EOMES, FGFBP2, IL2RB and S1PR5 were reduced expressed verified by RT-PCR, and PLEKHF1, S1PR5, and FGFBP2 was upregulated after TNF-α blocker therapy. Conclusion The study revealed the potential role of NK cells in the pathogenesis of AS and identified diagnostic biomarkers and potential treatment targets for AS.