Ankylosing spondylitis (AS) is a chronic immune-mediated inflammatory disease affecting the axial skeleton, characterized by progressive structural damage and functional impairment. Although biologic therapies targeting tumor necrosis factor and interleukin-17 have improved clinical outcomes, a substantial proportion of patients fail to achieve sustained disease control. Emerging evidence suggests that metabolic alterations may contribute to AS pathogenesis; however, systematic characterization of metabolism-related biomarkers and their regulatory networks remains limited, and the interplay between metabolic dysfunction and immune dysregulation in AS is poorly understood. Two whole-blood GEO datasets (GSE25101, GSE73754; n = 104) were integrated as the primary analytical cohort. A third dataset (GSE11886, n = 18; monocyte-derived macrophages) was included for exploratory cross-tissue analysis. Differential expression analysis identified 847 DEGs, which were refined to 16 metabolism-related genes through weighted gene co-expression network analysis (WGCNA) and GeneCards database filtering. Eleven machine learning algorithms with 5-fold cross-validation were applied to construct diagnostic models and identify hub genes. Validation analyses included immune cell infiltration estimation using CIBERSORT, metabolic pathway activity assessment via ssGSEA, single-cell transcriptomics from GSE268839, functional enrichment through GSEA/GSVA, and chromosomal localization analysis. A competing endogenous RNA (ceRNA) regulatory network was constructed to map post-transcriptional regulation. Natural compounds from 66 AS-treating traditional Chinese medicines were screened against hub genes using deep learning-based binding prediction. Multiple machine learning algorithms achieved comparable cross-validated performance (CV AUC range 0.741–0.836; top five models: 0.805–0.836) using the six hub genes (MFN2, SLC27A3, RHOB, SMG7, AKR1B1, LCOR) identified through SHAP-based feature importance analysis of the PLS model. Leave-one-dataset-out validation between the two whole-blood cohorts showed that all algorithms exceeded an AUC of 0.77 in Round 1 (validate: GSE73754, n = 72; best AUC 0.861), while Round 2 (validate: GSE25101, n = 32) yielded more modest performance (best AUC, 0.715) reflecting the smaller validation sample. Exploratory application to GSE11886 (macrophage-derived samples) showed near-chance performance, consistent with the tissue-source discrepancy. AS patients exhibited significant downregulation of oxidative phosphorylation, TCA cycle, and glycolysis pathways (p < 0.01), accompanied by elevated glutathione metabolism (p < 0.001). Immune cell deconvolution revealed reduced CD8+ T cell proportions correlating with MFN2 downregulation, and increased neutrophil frequencies correlating with SLC27A3 upregulation. Exploratory single-cell analysis indicated that RHOB expression was relatively enriched in border-associated macrophages and fibroblasts, while AKR1B1 was more prominently expressed in vascular endothelial cells and plasmacytoid dendritic cells. The ceRNA network identified 21 miRNAs and 65 lncRNAs forming 86 regulatory interactions, with four key regulatory axes (SATB1-AS1/miR-539-5p/LCOR, FAM95B1/miR-223-3p/RHOB, LINC01106/miR-106a-5p/MFN2, AATBC/miR-185-5p/SMG7) predicted to regulate hub gene expression. Compound screening identified betaine, pyruvic acid, citric acid, etc., as top-ranking candidates, with MFN2 showing the highest binding capacity among hub genes. This study provides an integrative framework linking metabolic reprogramming with immune dysfunction in AS. The six-gene diagnostic signature showed preliminary discriminatory ability in the available datasets, while the ceRNA regulatory network and natural compound screening results prioritize candidate regulatory pathways and compounds for future validation. These findings advance our understanding of AS pathogenesis and may guide future biomarker development and targeted intervention strategies.
Dan et al. (Mon,) studied this question.