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.

Exploring T Cell and NK Cell Involvement in Ankylosing Spondylitis Through Single‐Cell Sequencing

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

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

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

Diagnostic biomarkers and potential treatment targets for Ankylosing Spondylitis

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

The shared role of neutrophils in ankylosing spondylitis and ulcerative colitis

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

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

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

Integrated Proteomics and Single-Cell Mass Cytometry Analysis Dissects the Immune Landscape of Ankylosing Spondylitis.

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

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

Screening for Key Genes and Transcription Factors in Ankylosing Spondylitis by RNA-Seq

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

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

TCR Repertoire and CDR3 Motif Analyses Depict the Role of Αβ T Cells in Ankylosing Spondylitis.

Single‐cell analysis reveals innate immunity dynamics in ankylosing spondylitis

Identification of key genes in Ankylosing spondylitis.

Whole-blood Gene Expression Profiling in Ankylosing Spondylitis Identifies Novel Candidate Genes That May Contribute to the Inflammatory and Tissue-Destructive Disease Aspects.