Abstract:Background: Ankylosing spondylitis (AS) is a progressive rheumatic disease and studies reveal that the immune system is critical for the pathogenesis of AS. In the present study, various bioinformatics analysis methods were comprehensively applied, designed to identify potential key genes and inflammation states of AS. Methods: The transcriptome profiles of GSE25101 and GSE73754 obtained from the Gene Expression Omnibus (GEO) database were merged for subsequent analyses. The differentially expressed genes (DEGs) were identified using the Bioconductor package Limma and threshold values. Functional enrichment and pathway enrichment analyses were performed using the clusterProfiler package and Gene Set Enrichment Analysis (GSEA). Next, protein-protein interaction (PPI) network of the identified DEGs was constructed by the online database, the Search Tool for the Retrieval of Interacting Genes (STRING), visualization and analysis were performed through Cytoscape software. Subsequently, we applied CIBERSORT algorithm to identify subpopulation proportions of immune cells in peripheral blood samples. Finally, we validated the hub genes with the GSE18781 dataset. Samples were collected from patients to validate gene and protein expression using qRT-PCR and ELISA. Results: A total of 334 DEGs were identified, including 182 upregulated and 152 downregulated DEGs, between AS patients and normal human controls, which were primarily involved in immune response, autophagy, and natural killer cell-mediated cytotoxicity. The most prominent module and candidate biomarkers were identified from the PPI network. Biomarkers were selected for validation and their expressions were significantly decreased in peripheral blood samples which was consistent with transcriptome sequencing results. Nine genes with AUC > 0.70 were considered to be AS hub genes for ROC curve analysis, including GZMA, GZMK, PRF1, GNLY, NKG7, KLRB1, KLRD1, IL2RB and CD247. Furthermore, CIBERSORT results suggest that AS contained a higher proportion of CD8+ T cells, naive CD4+ T cells, neutrophils, and lower levels of gamma delta T cells compared with the normal controls. Conclusion: In this study, we identified DEGs combined with their closely related biological functions and propose that granule-associated proteins and immune infiltration maybe involved in the progression of ankylosing spondylitis. These validated hub genes may provide new perspectives for understanding the molecular mechanisms of ankylosing spondylitis.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Peripheral Mononuclear Cells in Patients With Ankylosing Spondylitis

Identification of the key genes and long non-coding RNAs in ankylosing spondylitis using RNA sequencing (vol 43, pg 1179, 2019)

Gene-regulatory network analysis of ankylosing spondylitis with a single-cell chromatin accessible assay

RNA Sequencing for Gene Expression Profiles in Peripheral Blood Mononuclear Cells with Ankylosing Spondylitis RNA

Diagnostic biomarkers and potential treatment targets for Ankylosing Spondylitis

Identification of immune related cells and crucial genes in the peripheral blood of ankylosing spondylitis by integrated bioinformatics analysis

Integrating multiomics and Single-Cell communication analysis to uncover Ankylosing spondylitis mechanisms

Dysregulated Non-Coding RNA Expression in T Cells from Patients with Ankylosing Spondylitis Contributes to Its Immunopathogenesis

Single‐cell RNA‐seq reveals altered NK cell subsets and reduced levels of cytotoxic molecules in patients with ankylosing spondylitis

Single Cell Transcriptome and Surface Epitope Analysis of Ankylosing Spondylitis Facilitates Disease Classification by Machine Learning

The shared role of neutrophils in ankylosing spondylitis and ulcerative colitis

Transcriptome-Wide Association Studies and Integration Analysis of mRNA Expression Profiles Identify Candidate Genes and Pathways Associated With Ankylosing Spondylitis

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Deciphering Metabolic Pathways and Protein-Protein Interaction Networks in Ankylosing Spondylitis through Single-Cell RNA Sequencing

CD96 Downregulation Promotes the Immune Response of CD4 T Cells and Associates with Ankylosing Spondylitis.

RNA sequencing and bioinformatics analysis of differentially expressed genes in the peripheral serum of ankylosing spondylitis patients

Upregulated of ANXA3, SORL1, and Neutrophils May Be Key Factors in the Progressionof Ankylosing Spondylitis

Identification of necroptosis‐related genes in ankylosing spondylitis by bioinformatics and experimental validation

Single-Cell Analysis of Patients with Axial Spondyloarthritis After Anti-TNFα Treatment: Experimental Data and Review of the Literature

Single-Cell Analysis of Patients with Axial Spondyloarthritis after Anti-TNF Treatment: Experimental Data and Review of the Literature

Identification of potential biomarkers for ankylosing spondylitis based on bioinformatics analysis