Abstract
BACKGROUND: The pathogenesis of ankylosing spondylitis (AS) is still not clear, and immune-related genes have not been systematically explored in AS. The purpose of this paper was to identify the potential early biomarkers most related to immunity in AS and develop a predictive disease risk model with bioinformatic methods and the Gene Expression Omnibus database (GEO) to improve diagnostic and therapeutic efficiency. METHODS: To identify differentially expressed genes and create a gene coexpression network between AS and healthy samples, we downloaded the AS-related datasets GSE25101 and GSE73754 from the GEO database and employed weighted gene coexpression network analysis (WGCNA). We used the GSVA, GSEABase, limma, ggpubr, and reshape2 packages to score immune data and investigated the links between immune cells and immunological functions by using single-sample gene set enrichment analysis (ssGSEA). The value of the core gene set and constructed model for early AS diagnosis was investigated by using receiver operating characteristic (ROC) curve analysis. RESULTS: Biological function and immune score analyses identified central genes related to immunity, key immune cells, key related pathways, gene modules, and the coexpression network in AS. Granulysin (GNLY), Granulysin (GZMK), CX3CR1, IL2RB, dysferlin (DYSF), and S100A12 may participate in AS development through NK cells, CD8(+) T cells, Th1 cells, and other immune cells and represent potential biomarkers for the early diagnosis of AS occurrence and progression. Furthermore, the T cell coinhibitory pathway may be involved in AS pathogenesis. CONCLUSION: The AS disease risk model constructed based on immune-related genes can guide clinical diagnosis and treatment and may help in the development of personalized immunotherapy.