{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,17]],"date-time":"2026-09-17T16:45:56Z","timestamp":1789663556315,"version":"4.0.1"},"reference-count":46,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2022,9,12]],"date-time":"2022-09-12T00:00:00Z","timestamp":1662940800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100012166","name":"National Key Research and Development Program of China","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100012166","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Immunol."],"abstract":"<jats:p>Ankylosing spondylitis (AS) is a common inflammatory spondyloarthritis affecting the spine and sacroiliac joint that finally results in sclerosis of the axial skeleton. Aside from human leukocyte antigen B27, transcriptomic biomarkers in blood for AS diagnosis still remain unknown. Hence, this study aimed to investigate credible AS-specific mRNA biomarkers from the whole blood of AS patients by analyzing an mRNA expression profile (GSE73754) downloaded Gene Expression Omnibus, which includes AS and healthy control blood samples. Weighted gene co-expression network analysis was performed and revealed three mRNA modules associated with AS. By performing gene set enrichment analysis, the functional annotations of these modules revealed immune biological processes that occur in AS. Several feature mRNAs were identified by analyzing the hubs of the protein-protein interaction network, which was based on the intersection between differentially expressed mRNAs and mRNA modules. A machine learning-based feature selection method, SVM-RFE, was used to further screen out 13 key feature mRNAs. After verifying by qPCR, IL17RA, Sqstm1, Picalm, Eif4e, Srrt, Lrrfip1, Synj1 and Cxcr6 were found to be significant for AS diagnosis. Among them, Cxcr6, IL17RA and Lrrfip1 were correlated with severity of AS symptoms. In conclusion, our findings provide a framework for identifying the key mRNAs in whole blood of AS that is conducive for the development of novel diagnostic markers for AS.<\/jats:p>","DOI":"10.3389\/fimmu.2022.956027","type":"journal-article","created":{"date-parts":[[2022,9,12]],"date-time":"2022-09-12T06:29:16Z","timestamp":1662964156000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":16,"title":["Identification of diagnostic mRNA biomarkers in whole blood for ankylosing spondylitis using WGCNA and machine learning feature selection"],"prefix":"10.3389","volume":"13","author":[{"given":"Yaguang","family":"Han","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yiqin","family":"Zhou","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Haobo","family":"Li","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Zhenyu","family":"Gong","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ziye","family":"Liu","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Huan","family":"Wang","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Bo","family":"Wang","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiaojian","family":"Ye","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yi","family":"Liu","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2022,9,12]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","DOI":"10.1002\/art.20990","article-title":"The challenge of diagnosis and classification in early ankylosing spondylitis: do we need new criteria","volume":"52","author":"Rudwaleit","year":"2005","journal-title":"Arthritis Rheum"},{"key":"B2","doi-asserted-by":"publisher","DOI":"10.3390\/cells8020108","article-title":"Age at disease onset and diagnosis delay in HLA-B27 negative vs Positive patients ankylosing spondylitis","volume":"23","author":"Feldtkeller","year":"2003","journal-title":"Rheum Int"},{"key":"B3","doi-asserted-by":"publisher","DOI":"10.1111\/1756-185X.14080","article-title":"Epigenetics of ankylosing spondylitis: Recent developments","volume":"24","author":"Yang","year":"2021","journal-title":"Int J Rheum Dis"},{"key":"B4","doi-asserted-by":"publisher","DOI":"10.1111\/1756-185X.12779","article-title":"The ongoing quest for biomarkers in ankylosing spondylitis","volume":"18","author":"Danve","year":"2015","journal-title":"Int J Rheum Dis"},{"key":"B5","doi-asserted-by":"publisher","first-page":"296","DOI":"10.1038\/nrrheum.2012.41","article-title":"Genetics of spondyloarthritis\u2013beyond","volume":"8","author":"Reveille","year":"2012","journal-title":"Nat Rev Rheum"},{"key":"B6","doi-asserted-by":"publisher","DOI":"10.1136\/ard.58.12.737","article-title":"Two forms of reactive arthritis","volume":"58","author":"Toivanen","year":"1999","journal-title":"Ann Rheum Dis"},{"key":"B7","doi-asserted-by":"publisher","DOI":"10.1136\/ard.2011.155044","article-title":"Human leucocyte antigen risk alleles for psoriatic arthritis among patients with psoriasis","volume":"71","author":"Eder","year":"2012","journal-title":"Ann Rheum Dis"},{"key":"B8","doi-asserted-by":"publisher","DOI":"10.1136\/ard.2007.071605","article-title":"Inflammatory biomarkers, disease activity and spinal disease measures in patients with ankylosing spondylitis after treatment with infliximab","volume":"67","author":"Visvanathan","year":"2008","journal-title":"Ann Rheum Dis"},{"key":"B9","doi-asserted-by":"publisher","DOI":"10.1002\/art.22589","article-title":"Serum matrix metalloproteinase 3 is an independent predictor of structural damage progression in patients with ankylosing spondylitis","volume":"56","author":"Maksymowych","year":"2007","journal-title":"Arthritis Rheum"},{"key":"B10","doi-asserted-by":"publisher","DOI":"10.3899\/jrheum.101128","article-title":"Serum MMP-3 level as a biomarker for monitoring and predicting response to etanercept treatment in ankylosing spondylitis","volume":"38","author":"Arends","year":"2011","journal-title":"J Rheum"},{"key":"B11","doi-asserted-by":"publisher","first-page":"41","DOI":"10.1016\/j.gene.2017.09.049","article-title":"Construction and analysis of protein-protein interaction network correlated with ankylosing spondylitis","volume":"638","author":"Kanwal","year":"2018","journal-title":"Gene"},{"key":"B12","doi-asserted-by":"publisher","first-page":"559","DOI":"10.1186\/1471-2105-9-559","article-title":"WGCNA: an r package for weighted correlation network analysis","volume":"9","author":"Langfelder","year":"2008","journal-title":"BMC Bioinf"},{"key":"B13","doi-asserted-by":"publisher","DOI":"10.3389\/fphys.2018.00399","article-title":"Identification of biomarkers associated with pathological stage and prognosis of clear cell renal cell carcinoma by Co-expression network analysis","volume":"9","author":"Chen","year":"2018","journal-title":"Front Physiol"},{"key":"B14","doi-asserted-by":"publisher","first-page":"505701","DOI":"10.1155\/2009\/505701","article-title":"Global expression patterns of three festuca species exposed to different doses of glyphosate using the affymetrix GeneChip wheat genome array","volume":"2009","author":"Cebeci","year":"2009","journal-title":"Comp Funct Genomics"},{"key":"B15","doi-asserted-by":"publisher","DOI":"10.1109\/TNB.2005.853657","article-title":"Multiple SVM-RFE for gene selection in cancer classification with expression data","volume":"4","author":"Duan","year":"2005","journal-title":"IEEE Trans Nanobioscience"},{"key":"B16","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btm036","article-title":"MSVM-RFE: extensions of SVM-RFE for multiclass gene selection on DNA microarray data","volume":"23","author":"Zhou","year":"2007","journal-title":"Bioinformatics"},{"key":"B17","doi-asserted-by":"publisher","DOI":"10.1002\/art.1780270401","article-title":"Evaluation of diagnostic criteria for ankylosing spondylitis. a proposal for modification of the new York criteria","volume":"27","author":"van der Linden","year":"1984","journal-title":"Arthritis Rheum"},{"key":"B18","doi-asserted-by":"publisher","DOI":"10.1002\/art.39464","article-title":"Sexual dimorphism in the Th17 signature of ankylosing spondylitis","volume":"68","author":"Gracey","year":"2016","journal-title":"Arthritis Rheum"},{"key":"B19","doi-asserted-by":"publisher","DOI":"10.1023\/A:1008334909089","article-title":"Imputation of missing data in industrial databases","volume":"11","author":"Lakshminarayan","year":"1999","journal-title":"Appl Intell"},{"key":"B20","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btg405","article-title":"Affy\u2013analysis of affymetrix GeneChip data at the probe level","volume":"20","author":"Gautier","year":"2004","journal-title":"Bioinformatics"},{"key":"B21","doi-asserted-by":"publisher","first-page":"1989","DOI":"10.7150\/jca.23762","article-title":"Development of a three-gene prognostic signature for hepatitis b virus associated hepatocellular carcinoma based on integrated transcriptomic analysis","volume":"9","author":"Yang","year":"2018","journal-title":"J Cancer"},{"key":"B22","doi-asserted-by":"publisher","first-page":"2344","DOI":"10.3390\/ijms20092344","article-title":"Predicting apoptosis protein subcellular locations based on the protein overlapping property matrix and tri-gram encoding","volume":"20","author":"Yang","year":"2019","journal-title":"Int J Mol Sci"},{"key":"B23","first-page":"389","article-title":"Gene selection for cancer classification using support vector machines","volume-title":"Machine Learning","author":"Guyon","year":"2002"},{"key":"B24","doi-asserted-by":"publisher","DOI":"10.2202\/1544-6115.1128","article-title":"A general framework for weighted gene co-expression network analysis","volume":"4","author":"Zhang","year":"2005","journal-title":"Stat Appl Genet Mol Biol"},{"key":"B25","doi-asserted-by":"publisher","DOI":"10.1109\/TNB.2014.2337912","article-title":"Prediction of essential proteins based on overlapping essential modules","volume":"13","author":"Zhao","year":"2014","journal-title":"IEEE Trans Nanobioscience"},{"key":"B26","doi-asserted-by":"publisher","DOI":"10.1002\/art.1780401015","article-title":"Susceptibility to ankylosing spondylitis in twins: the role of genes","volume":"40","author":"Brown","year":"1997","journal-title":"HLA Environ. Arthritis Rheum"},{"key":"B27","doi-asserted-by":"publisher","DOI":"10.3899\/jrheum.130515","article-title":"Higher expression of whole blood microRNA-21 in patients with ankylosing spondylitis associated with programmed cell death 4 mRNA expression and collagen cross-linked c-telopeptide concentration","volume":"41","author":"Huang","year":"2014","journal-title":"J Rheum"},{"key":"B28","doi-asserted-by":"publisher","DOI":"10.1136\/bmj.38954.689583.DE","article-title":"Diagnosis and management of ankylosing spondylitis","volume":"333","author":"McVeigh","year":"2006","journal-title":"BMJ"},{"key":"B29","doi-asserted-by":"publisher","DOI":"10.1111\/1756-185X.14004","article-title":"Chaperone\/scaffolding\/adaptor protein 14-3-3\u03b7 (eta): A diagnostic marker of rheumatoid arthritis","volume":"23","author":"Wei","year":"2020","journal-title":"Int J Rheum Dis"},{"key":"B30","doi-asserted-by":"publisher","first-page":"387","DOI":"10.1038\/s41584-021-00625-y","article-title":"Ankylosing spondylitis: an autoimmune or autoinflammatory disease","volume":"17","author":"Mauro","year":"2021","journal-title":"Nat Rev Rheum"},{"key":"B31","doi-asserted-by":"publisher","first-page":"5","DOI":"10.1111\/nyas.13933","article-title":"Targeting tumor necrosis factor receptors in ankylosing spondylitis","volume":"1442","author":"Lata","year":"2019","journal-title":"Ann N Y Acad Sci"},{"key":"B32","doi-asserted-by":"publisher","DOI":"10.1016\/j.biopha.2018.05.097","article-title":"Role of innate immune system in the pathogenesis of ankylosing spondylitis","volume":"105","author":"Vanaki","year":"2018","journal-title":"BioMed Pharmacother"},{"key":"B33","doi-asserted-by":"publisher","DOI":"10.1074\/jbc.M106339200","article-title":"Expression of the osteoblast differentiation factor RUNX2 (Cbfa1\/AML3\/Pebp2alpha a) is inhibited by tumor necrosis factor-alpha","volume":"277","author":"Gilbert","year":"2002","journal-title":"J Biol Chem"},{"key":"B34","doi-asserted-by":"publisher","DOI":"10.1152\/ajpendo.00534.2004","article-title":"The p55 TNF receptor mediates TNF inhibition of osteoblast differentiation independently of apoptosis","volume":"288","author":"Gilbert","year":"2005","journal-title":"Am J Physiol Endocrinol Metab"},{"key":"B35","doi-asserted-by":"publisher","first-page":"33","DOI":"10.1016\/S1043-4666(03)00106-6","article-title":"Tumor necrosis factor-alpha inhibits pre-osteoblast differentiation through its type-1 receptor","volume":"22","author":"Abbas","year":"2003","journal-title":"Cytokine"},{"key":"B36","first-page":"25","article-title":"Gene selection by sequential search wrapper approaches in microarray cancer class prediction","volume":"12","author":"Inza","year":"2002","journal-title":"J Intell  Fuzzy Syst Appl Eng Technol"},{"key":"B37","doi-asserted-by":"publisher","DOI":"10.1038\/nrd3794","article-title":"Targeting IL-17 and TH17 cells in chronic inflammation","volume":"11","author":"Miossec","year":"2012","journal-title":"Nat Rev Drug Discovery"},{"key":"B38","doi-asserted-by":"publisher","DOI":"10.1136\/annrheumdis-2019-215356","article-title":"The role of IL-17A in axial spondyloarthritis and psoriatic arthritis: recent advances and controversies","volume":"78","author":"McGonagle","year":"2019","journal-title":"Ann Rheum Dis"},{"key":"B39","doi-asserted-by":"publisher","DOI":"10.1002\/art.24968","article-title":"Identification of RGS1 as a candidate biomarker for undifferentiated spondylarthritis by genome-wide expression profiling and real-time polymerase chain reaction","volume":"60","author":"Gu","year":"2009","journal-title":"Arthritis Rheum"},{"key":"B40","doi-asserted-by":"publisher","DOI":"10.1002\/eji.201344184","article-title":"Interleukin 17 contributes to the chronicity of inflammatory diseases such as rheumatoid arthritis","volume":"44","author":"Benedetti","year":"2014","journal-title":"Eur J Immunol"},{"key":"B41","doi-asserted-by":"publisher","DOI":"10.1016\/S0140-6736(18)31946-9","article-title":"Ixekizumab, an interleukin-17A antagonist in the treatment of ankylosing spondylitis or radiographic axial spondyloarthritis in patients previously untreated with biological disease-modifying anti-rheumatic drugs (COAST-v): 16 week results of a phase 3 randomised, double-blind, active-controlled and placebo-controlled trial","volume":"392","author":"van der Heijde","year":"2018","journal-title":"Lancet"},{"key":"B42","doi-asserted-by":"publisher","first-page":"599","DOI":"10.1002\/art.40753","article-title":"Efficacy and safety of ixekizumab in the treatment of radiographic axial spondyloarthritis: Sixteen-week results from a phase III randomized, double-blind, placebo-controlled trial in patients with prior inadequate response to or intolerance of tumor necrosis factor inhibitors","volume":"71","author":"Deodhar","year":"2019","journal-title":"Arthritis Rheum"},{"key":"B43","doi-asserted-by":"publisher","first-page":"1566","DOI":"10.1136\/annrheumdis-2012-202925","article-title":"Evidence that autophagy, but not the unfolded protein response, regulates the expression of IL-23 in the gut of patients with ankylosing spondylitis and subclinical gut inflammation","volume":"73","author":"Ciccia","year":"2014","journal-title":"Ann Rheum Dis"},{"key":"B44","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pone.0097701","article-title":"The function of the chemokine receptor CXCR6 in the T cell response of mice against listeria monocytogenes","volume":"9","author":"Heesch","year":"2014","journal-title":"PloS One"},{"key":"B45","doi-asserted-by":"publisher","DOI":"10.1073\/pnas.1700109114","article-title":"Brd4 modulates the innate immune response through Mnk2-eIF4E pathway-dependent translational control of I\u03baB\u03b1","volume":"114","author":"Bao","year":"2017","journal-title":"Proc Natl Acad Sci USA"},{"key":"B46","doi-asserted-by":"publisher","first-page":"108","DOI":"10.3390\/cells8020108","article-title":"Multidisciplinary roles of LRRFIP1\/GCF2 in human biological systems and diseases","volume":"8","author":"Takimoto","year":"2019","journal-title":"Cells"}],"container-title":["Frontiers in Immunology"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fimmu.2022.956027\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,9,12]],"date-time":"2022-09-12T06:29:26Z","timestamp":1662964166000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fimmu.2022.956027\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,9,12]]},"references-count":46,"alternative-id":["10.3389\/fimmu.2022.956027"],"URL":"https:\/\/doi.org\/10.3389\/fimmu.2022.956027","relation":{},"ISSN":["1664-3224"],"issn-type":[{"value":"1664-3224","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,9,12]]},"article-number":"956027"}}