Skip to main content
  • Research article
  • Open access
  • Published:

Association between ankylosing spondylitis and m6A methylation



N6-methyl adenosine (m6A) is the most common reversible mRNA modification in eukaryotes implicated in key roles in various biological processes. The purpose of our analysis was to examine the association of ankylosing spondylitis (AS) with m6A methylation.


We obtained 72 samples from the data set GSE73754, including 52 AS patients and 20 healthy people. We divided the samples into two groups: the experimental group and the control group, and then observed the differences of 26 m6A related genes in the two groups. We also analyzed the correlation between different m6A genes. We used a random forest tree model to screen seven m6A signature genes associated with AS to evaluate its prevalence. Next, the samples were classified according to the m6a content and differential genes. Immune analysis, gene ontology, and KEGG enrichment analyses were performed. Finally, we scored each sample with m6a and analyzed the relationship between different samples and inflammation-related factors.

Results and conclusion

In conclusion, we screened out AS-related genes and the nomogram showed that they were negatively correlated with the incidence of AS. And we found that AS may have some relationship with immunity. Our analysis results could provide further insights into the treatment of AS.


Ankylosing spondylitis (AS), a common clinical type of spondyloarthropathy, is an autoimmune disease characterized by chronic inflammation of axial joints and tendon–ligament–bone attachments. It mainly invades the spine and sacroiliac joints accompanied by multi-system damage. As the disease progresses, patients may experience limited spinal movement, stiffness, and destruction of the hip joint. In severe cases, it can lead to disability and seriously affect the patient's quality of life [1]. The pathogenesis of AS is complex with genetic factors, interactions between microorganisms and the host, differences in hormone levels, imbalance in the intestinal flora, and abnormal immune responses implicated in the process [2,3,4]. At present, it is believed that disorders of the immune system play drive the occurrence and development of AS. Several studies have shown that both innate and adaptive immune responses are involved in the pathological processes underlying AS.

Similar to DNA methylation, N6-methyl adenosine (m6A) RNA methylation modifications are dynamically reversible and can be catalyzed by methyltransferases (“writers”) or erased by demethylases (“erasers”). Reads are recognized by a set of binding proteins ("readers") with said functions [5,6,7]. Although the series of processes do not change the principle of complementary base pairing, they functionally (including RNA splicing, transport, and translation) determine the direction of RNA molecules. Thus, m6A methylation modifications are involved in various diseases, especially key in tumor progression [8, 9]. In recent years, studies have shown that m6A methylation regulators are significantly differentially expressed in several malignant tumors, and thus, the influence of abnormal m6A methylation mediated by regulatory factors on various malignant tumors has received extensive attention [10,11,12]. However, m6A methylation regulators have not been studied in AS, and we intended to elucidate the relationship between m6A methylation and AS.

In this study, we analyzed the expression of m6A methylation regulators in the samples and screened seven disease-related risk genes using the random forest tree model. We observed the relationship of immune cells with immunity in different samples. Our purpose was to examine the correlation between AS and m6A methylation; our findings are expected to provide insight into the diagnosis and treatment of AS in the future.


Data processing and downloading

The GSE73754 dataset was obtained from the Gene Expression Omnibus (GEO) database. A total of 72 samples, including 52 from AS patients and 20 from healthy individuals, were extracted. A total of 26 m6A-related genes including writers (METTL3, METTL14, METTL16, WTAP, VIRMA, ZC3H13, RBM15, RBM15B, CBLL1), readers (YTHDC1, YTHDC2, YTHDF1, YTHDF2, YTHDF3, HNRNPC, FMR1, LRPPRC, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, RBMX, ELAVL1, IGF2BP1), and erasers (FTO, ALKBH5) were screened.

Consensus cluster and immune cell infiltration analyses

In order to comprehensively analyze the relationship between m6A methylation and AS, we used the "ConsensusClusterPlus" package in R to perform a consistent cluster analysis for the patients in the dataset. According to the Cumulative Distribution Function (CDF) curve, stratification was performed based on the optimal “k” value. Principal component analysis (PCA) was used to further confirm the rationality of the consistent cluster-based grouping, and differential analysis for these subtypes was performed. The abundance of immune cells was assessed by single-sample gene set enrichment analysis (ssGSEA) to further investigate correlations, and graphs were plotted to visualize features.

Random forest model (RF) and support vector machine model (SVM)

We utilized RF and SVM to predict the onset of AS. We built the former using the "RandomForest" package, assessed the importance of each m6A regulator, and screened seven m6A regulators. For SVM, four features were captured, including maximum margin hyperplane, soft margin, kernel function, and separating hyperplane. Each data point was plotted in the n-dimensional space, and a suitable plane was found that well-distinguished the two groups. We again validate these models using inverse cumulative distribution of residuals, box plots of residuals, and receiver operating characteristic (ROC) curves.

Differential gene expression and functional enrichment between groups

Differentially expressed genes (DEGs) in different clusters were screened using the "limma" package in R, with a significance criterion of P < 0.05. The "clusterprofiler" package in R was employed to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses for the screened DEGs.

Calculation of the m6A score

We quantified subtypes by PCA for m6A-based scoring [13]. The formula for the m6A score was as follows: m6A score = PC1i, where PC1 refers to principal component 1 and i represents the level of the corresponding DEG [14].

Statistical analysis

Statistical analysis was performed using R (version 4.2.0). Linear regression analysis and Pearson correlation coefficient (r) were used to determine the correlation between gene expression patterns. Nonparametric one-way analysis of variance (ANOVA) was used to compare the variables between different groups. The comparison between two groups was made using the t-test. A value of P < 0.05 was indicative of statistical significance.


Overview of m6a regulators in AS

In the GSE73754 dataset, we divided the samples into an experimental group and a control group and observed the differences in the expression of 26 m6A-related genes between them. The expression of seven of the 26 m6a regulators differed significantly between the two groups (Fig. 1A, B). We constructed a circle diagram, and the locations of the 26 m6A-related genes on the chromosomes are shown (Fig. 1C). We also analyzed the associations among different m6A genes and found a certain linear correlation among CBLL1, RBM15B, and FTO, whereby increased expressions of CBLL1 and RBM15B corresponded to high levels of FTO (Fig. 1D, E).

Fig. 1
figure 1

m6a regulators in AS. A Heat map of differentially expressed m6a regulators between the two groups. B Boxplot of differential expression of 26 m6a regulators between the two groups. C Chromosomal location of m6a-related genes. D The relationship between FTO and CBLL1. E The relationship between FTO and RBM15B

Selection of the RF and SVM models

Next, we compared the RF and SVM models to evaluate the most suitable method for screening disease eigengenes. The residual value of the RF model was lower according to the "inverse cumulative distribution of residuals" and "boxplot of residuals", indicating that it was more suitable for screening AS characteristic genes as compared to the SVM model (Fig. 2 A, C). The area under the ROC curve of the RF model was higher than that of the SVM model, suggesting greater accuracy of the former as compared to the latter. We thus chose the RF model for subsequent analysis (Fig. 2B). Next, we screened AS signature genes using the RF model and found the point with the smallest cross-validation error in the control group, experimental group, and all samples. Subsequently, we found the number of trees corresponding to the point. Ultimately, we identified seven AS signature genes: IGFBP1、YTHDC1、METTL14、CBLL1、RBMX、HNRNPC、ALKBH5. And we scored these based on their importance (Fig. 2D,E).

Fig. 2
figure 2

Construction of the RF and SVM models. A Boxplots of residual values of RF and SVM models. B ROC curves for RF and SVM models. C Inverse cumulative distribution plot of residuals of RF and SVM models. D Genes in the RF model ranked according to their scores, suggestive of importance; the top three are RBMX, METTL14, and ALKBH5. E The influence of the number of decision trees on the error

Nomogram of AS signature genes

Next, we constructed a nomogram based on the screened AS characteristic genes. According to the expression of different genes, we calculated the corresponding scores and inferred the incidence of AS according to these scores (Fig. 3A). According to the nomogram, a positive correlation between ALKBH5 and AS, while a negative correlation between RBMX, METTL14, HNRNPC, IGFBP1, YTHDC1, and CBLL1 and AS was observed. The distance between the solid and dashed lines in the calibration curve was relatively close, indicating a substantially high accuracy of the model (Fig. 3D). The decision and clinical impact curves also demonstrated that the constructed nomogram was reliable (Fig. 3 B,C).

Fig. 3
figure 3

The establishment of the nomogram. A Nomogram constructed using seven characteristic genes screened based on the RF model. B Clinical impact curve to assess model’s accuracy. C Decision curve. The farther the red curve is from the gray curve, the higher the model’s accuracy. D Calibration curve. The closer the solid and dashed lines, the more accurate the model is

m6a typing

According to differential m6a expressions across samples, we performed a genotyping analysis. At k = 2, the CDF showed the maximum value; thus, we divided all samples into two subtypes (Fig. 4A-C). Subsequent PCA demonstrated that all samples could be distinguished based on the corresponding m6a levels, suggesting that our typing was reliable (Fig. 4D). From the boxplot and gene heat map, the expression of most AS characteristic genes in cluster A was found to be relatively high, while that of the genes in cluster B was relatively low. Significant differences in the expressions of CBLL1, HNRNPC, RBMX, and ALKBH5 were found between the two clusters (Fig. 4E, F).

Fig. 4
figure 4

Typing according to m6a-related gene expression. A Consensus cluster analysis. B CDF curve for consensus clustering. C At k = 2, the CDF showed the maximum value. D PCA shows a significant difference between clusters A and B. E Boxplot showing the differences in m6a signature gene expressions between phenotypes. F Heat map of the differences in m6a signature gene expression between different phenotypes

Characteristics of m6a typing

Next, immune cell analysis based on m6a typing was performed to examine the link between m6a gene signature and immunity. Cluster A showed higher infiltration of B-, CD4 T-, CD8 T-, natural killer-, and mast cells relative to cluster B, while the latter showed a higher neutrophil abundance (Fig. 5A). We also analyzed the correlation between m6a signature genes and immune cells. As shown in the figure, RBMX exhibited the highest correlation with immune cells and was used for subsequent analysis (Fig. 5B). The proportions of immature B cells, natural killer cells, CD4T cells, and B cells in the RBMX high-expression group were higher than those in the RBMX low-expression group (Fig. 5C). We identified 38 DEGs between clusters A and B and drew a Venn diagram (Fig. 5D). For these DEGs, we performed GO and KEGG enrichment analyses (Fig. 5E, F). The processes of erythrocyte development, myeloid cell development, and erythrocyte differentiation were enriched. The results of the KEGG enrichment analysis are shown in figure.

Fig. 5
figure 5

Biological characteristics of m6a typing. A Differential analysis for immune cells based on m6a typing. B Immune cell correlation analysis for m6a typing; RBMX shows the highest correlation with immune cells. C Differential analysis for immune cells between the RBMX high- and RBMX low-expression groups. D Venn diagram of DEGs in clusters A and B. E GO enrichment analysis for m6a typing. F KEGG enrichment analysis for m6a typing

Features of DEG typing

Next, we classified the samples according to DEG levels. At k = 2, CDF showed the maximum value. We, therefore, divided the sample into two groups, whereby the between-group difference was the smallest and the within-group difference was the largest (Fig. 6 A-C). As shown in the heat map, most genes in cluster B show higher expression relative to cluster A (Fig. 6D). m6a differential analysis for genotyping was performed, and most of the m6a signature genes in cluster A showed higher expressions relative to cluster B (Fig. 6E). Finally, we analyzed immune cell proportions between the two clusters. In cluster A, the abundances of CD4 T and CD8 T cells were higher than those in cluster B, while the proportions of dendritic and delta T cells were lower than those in cluster B (Fig. 6F).

Fig. 6
figure 6

Typing according to DEG levels. A Consensus cluster analysis. B CDF curve for consensus cluster analysis. C At k = 2, the CDF shows the maximum value. D The heat map of DEGs between clusters A and B. E Boxplots showing differences in m6a signature gene expression between different types. F Differences in immune cells between different types

m6a score-based signature and inflammatory factors

PCA is an efficient statistical analysis technique to quantify m6a signatures. Based on the results of our analysis, m6a cluster A showed a lower score than m6a cluster B, while there were no differences between the scores of gene clusters A and B (Fig. 7A, B). In the Sankey diagram, a certain similarity between the results of m6a typing and genotyping was observed (Fig. 7C). Inflammatory factors and immunity are closely related. Finally, we analyzed the relationship between different types and inflammatory factors. The data showed no differences in the expression of inflammatory factors among subtypes (Fig. 7D, E).

Fig. 7
figure 7

m6a scores for different types. A Differences in m6a scores of genotypes. B Differences in m6a scores for m6a types. C Sankey diagram for different types and m6a scores. D Differences in the expression of inflammatory factors based on m6a typing. E Differences in the expression of inflammatory factors based on gene typing


AS is a chronic autoimmune disease [15], the etiology of which remains clear. Currently, it is believed that AS is related to genetic, infection, immune, and endocrine factors [16,17,18]. AS involves the spine, sacroiliac joints, hip joints, and other bones, and can also include tendons and ligaments around the joints, resulting in ossification and ankylosis of these structures [19, 20]. The disease causes physical dysfunction in patients and simultaneous serious psychological problems, including loss of confidence in life, resulting in a heavy burden to the family and society [21, 22]. As an important modification in the transcriptome, m6a methylation affects the stability, splicing, and translation of mRNA, thus playing an important role in post-transcriptional regulation. Previous studies show that regulation of m6A methylation modification is closely related to the occurrence and development of tumors via tumor-promoting effects through several tumor-related molecules or signaling pathways [23]. Methylation of m6A can affect the progression of osteosarcoma by affecting tumor cell proliferation, apoptosis, invasion ability, and migration ability [24]. And m6A methylation can have a certain impact on osteoporosis by regulating osteoblast differentiation and regulating osteoclast bone resorption function [25, 26]. However, the relationship between m6a methylation and AS is unknown. In this study, we used the RF model to construct a nomogram to estimate the incidence of AS; we hypothesized that m6a methylation may be related to AS.

We screened seven m6a characteristic genes using the RF model, namely ALKBH5, RBMX, METTL14, HNRNPC, IGFBP1, YTHDC1, and CBLL1. Among them, three genes, RBMX, METTL14, and ALKBH5, showed the highest scores, suggestive of their importance. RBMX, a heterogeneous ribonucleoprotein, is mainly localized in the nucleus and consists of 391 amino acid residues with an RNA-binding domain at the amino terminus of the protein [27, 28]. This structural feature reveals the physical association of RBMX with RNA, and the binding region is required for RNA processing and metabolism [29]. Therefore, the level of RBMX expression and its biological activity may affect the overall gene expression in cells to maintain normal cellular homeostasis [21]. Recent studies have implicated RBMX protein expression in apoptosis [30]. Previous studies show that all X-chromosome-located RBMX genes are significantly up-regulated in breast cancer and correlate positively with the expression of the pro-apoptotic gene, BAX [31]. Based on these previous findings, RBMX plays an important role in apoptosis or programmed cell death, which may be related to AS. Recently, m6A modification mediated by METTL14 has been implicated in the occurrence and development of tumors, and METTL14 plays a role in promoting cancer. For example, m6A methylation levels are significantly higher in approximately 70% of pancreatic cancer samples, and dysregulation of METTL14 can affect m6A levels in pancreatic cancer cells [32]. The expression of METTL14 is also up-regulated in breast cancer tissues, thus promoting breast cancer cell migration and invasion by regulating the expression of hsa-miR-146a-5p [33]. ALKBH5 demethylation plays an important role in various life processes, especially its complex role in tumors, which is worthy of further discussion. ALKBH5 can stabilize the expression of the target gene FOXM1 mRNA by reducing the levels of m6A, thereby inducing the occurrence of glioblastoma [33]. LKBH5 can also enhance the stability of BCL-2 mRNA through demethylation and promote the proliferation and invasion of ovarian cancer cells [34]. Thus, RBMX, METTL14, and ALKBH5, the three m6a characteristic genes in this study, may be related to AS. The nomogram constructed showed a negative correlation of these genes with the incidence of AS. Immune correlation analysis showed the association of these three genes with immunity, especially RBMX, which may underlie their importance in AS.

AS is an autoimmune disease. Several inflammatory T cells, monocytes, and macrophages infiltrate bones, joints, and synovial tissues of AS patients, causing progressive inflammation. Significantly increased levels of alpha globulin, gamma globulin, immunoglobulin IgG, IgA, IgM, serum complement C3, and complement C4 expressions, and Thl/Th2 and Thl7/Treg ratios have been reported in patients with AS [35, 36]. The pathogenesis of AS is related to the abnormal function of B cells. The proportion of CD19 + B cells in AS patients with active phase or peripheral joint involvement was higher than that in patients with stable phase or axial joint involvement, and the proportion of CD19 + B cells in active AS patients after 12 weeks of etanercept treatment did not decrease [37]. Another study showed that the levels of expression of BAFF and its receptor, BR3, in peripheral blood leukocytes of AS patients were significantly up-regulated, and the two were significantly positively correlated; the abundance of CD19 + CD38 + antibody-secreting cells also increased [38]. Our data suggested that both m6a cluster A and gene cluster A had greater immune infiltration relative to m6a cluster B and gene cluster B, which may be due to AS prevalence.

In conclusion, we assessed the relationship between m6a methylation and AS, identified the characteristic genes of m6a, and discussed their relationship with immunity. Our findings may provide strategies for the treatment of AS in the future.


  1. Navid D, Mahdi M, Saeed A, Ahmadreza J. HLA-B*27 subtypes and their implications in the pathogenesis of ankylosing spondylitis. Gene. 2018;670:15–21.

    Article  Google Scholar 

  2. Zhang L, Zhang YJ, Chen J, et al. The association of HLA-B27 and Klebsiella pneumoniae in ankylosing spondylitis: a systematic review. Microb Pathog. 2018;117:49.

    Article  CAS  PubMed  Google Scholar 

  3. Tapia-Serrano R, Jimenez-Balderas FJ, Murrieta S, Bravo-Gatica C, Mintz G. Testicular function in active ankylosing spondylitis. Therapeutic response to human chorionic gonadotrophin. J Rheumatol. 1991;18(6):841.

    CAS  PubMed  Google Scholar 

  4. Manetti R. Natural killer cell stimulatory factor (interleukin 12 [IL-12]) induces T helper type 1 (Th1)-specific immune responses and inhibits the development of IL-4-producing Th cells. J Exp Med. 1993;177(4):1199–204.

    Article  CAS  PubMed  Google Scholar 

  5. Tan F, Zhao M, Xiong F, Wang Y, Guo C. N6-methyladenosine-dependent signalling in cancer progression and insights into cancer therapies. J Exp Clin Cancer Res. 2021;40(1):526.

    Article  Google Scholar 

  6. Li Y, Xiao J, Bai J, et al. Molecular characterization and clinical relevance of m6A regulators across 33 cancer types. Mol Cancer. 2019;6:18.

    Article  Google Scholar 

  7. Huang J, Dong X, Gong Z, et al. Solution structure of the RNA recognition domain of METTL3-METTL14 N(6)-methyladenosine methyltransferase. Protein Cell. 2019;10(4):272–84.

    Article  CAS  PubMed  Google Scholar 

  8. Zeng C, Huang W, Li Y, Weng H. Roles of METTL3 in cancer: mechanisms and therapeutic targeting. J Hematol Oncol. 2020;13(1):117.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Wang X, Feng J, Xue Y, et al. Structural basis of N(6)-adenosine methylation by the METTL3-METTL14 complex. Nature. 2016;534(7608):575–8.

    Article  CAS  PubMed  Google Scholar 

  10. Ping XL, Sun BF, Wang L, et al. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Nature. 2014;2:42.

    Google Scholar 

  11. Patil DP, et al. m(6)A RNA methylation promotes XIST-mediated transcriptional repression. Nature. 2016;537:369–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Yue Y, Liu J, Cui X. VIRMA mediates preferential mA mRNA methylation in 3 UTR and near stop codon and associates with alternative polyadenylation. 2018.

  13. Hnzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-Seq data. BMC Bioinform. 2013;14(1):7–7.

    Article  Google Scholar 

  14. Wang H, Zhou L. Random survival forest with space extensions for censored data. Artif Intell Med. 2017;5:69.

    Google Scholar 

  15. Ranganathan VG, et al. Pathogenesis of ankylosing spondylitis-recent advances and future directions. Science. 2010;5:963.

    Google Scholar 

  16. Goodman SM, Figgie M. Lower extremity arthroplasty in patients with inflammatory arthritis: preoperative and perioperative management. J Am Acad Orthop Surg. 2013;21(6):355–63.

    PubMed  Google Scholar 

  17. Moshirfar A, Rand FF, Kebaish KM. Fixed sagittal plane deformity: causes, prevention, and treatment options. Semin Spine Surg. 2011;23(2):135–41.

    Article  Google Scholar 

  18. Charles YP, Ntilikina Y, Collinet A, Steib JP. Combined percutaneous and open instrumentation for thoracolumbar kyphosis correction by two-level pedicle subtraction osteotomy in ankylosing spondylitis. Eur J Orthop Surg Traumatol. 2020;30(5):939–47.

    Article  PubMed  Google Scholar 

  19. Joshi AB. Total hip arthroplasty in ankylosing spondylitis. J Arthroplasty. 2002;17(4):427–33.

    Article  PubMed  Google Scholar 

  20. Guan M, Wang J, Zhao L, Xiao J, Li Z, Shi Z. Management of hip involvement in ankylosing spondylitis. Clin Rheumatol. 2013;32(8):1115–20.

    Article  PubMed  Google Scholar 

  21. Rudwaleit M, Landewe R, Van D, et al. The development of Assessment of SpondyloArthritis international Society classification criteria for axial spondyloarthritis (part I): classification of paper patients by expert opinion including uncertainty appraisal. Ann Rheum Dis. 2009;68(6):770–6.

    Article  CAS  PubMed  Google Scholar 

  22. Qian BP, Jiang J, Qiu Y, Wang B, Yu Y, Zhu ZZ. Radiographical predictors for postoperative sagittal imbalance in patients with thoracolumbar kyphosis secondary to ankylosing spondylitis after lumbar pedicle subtraction osteotomy. Spine. 2013;38(26):E1669.

    Article  PubMed  Google Scholar 

  23. Huang W, Chen TQ, Fang K, Zeng ZC, Ye H, Chen YQ. N6-methyladenosine methyltransferases: functions, regulation, and clinical potential. J Hematol Oncol. 2021;14(1):117.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wu Y, Wang Z, Shen J, et al. The role of m6A methylation in osteosarcoma biological processes and its potential clinical value. Hum Genomics. 2022;16(1):12.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wang W, Qiao SC, Wu XB, et al. Circ_0008542 in osteoblast exosomes promotes osteoclast-induced bone resorption through m6A methylation. Cell Death Dis. 2021;12(7):628.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Huang M, Xu S, Liu L, et al. m6A methylation regulates osteoblastic differentiation and bone remodeling. Front Cell Dev Biol. 2021;9:783322.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Shin KH, Kim RH, Kang MK, et al. p53 promotes the fidelity of DNA end-joining activity by, in part, enhancing the expression of heterogeneous nuclear ribonucleoprotein G. DNA Repair. 2007;6(6):830–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Shin KH, Mo KK, Park NH. Heterogeneous nuclear ribonucleoprotein G, nitric oxide, and oral carcinogenesis. Nitric Oxide. 2008;19(2):125–32.

    Article  CAS  PubMed  Google Scholar 

  29. Shin K-H. Heterogeneous nuclear ribonucleoprotein g shows tumor suppressive effect against oral squamous cell carcinoma cells. Clin Cancer Res Off J Am Assoc Cancer Res. 2006;12(10):3222–8.

    Article  CAS  Google Scholar 

  30. Sutherland LC, Rintala-Maki ND, White RD, Morin CD. RNA binding motif (RBM) proteins: a novel family of apoptosis modulators? J Cell Biochem. 2005;94(1):63.

    Article  Google Scholar 

  31. Martínez-Arribas F, Agudo D, Pollán M, et al. Positive correlation between the expression of X-chromosome RBM genes (RBMX, RBM3, RBM10) and the proapoptotic Bax gene in human breast cancer. Science. 2006;97(6):1275–82.

    Google Scholar 

  32. Wang M, Liu J, Zhao Y, He R, Qin R. Upregulation of METTL14 mediates the elevation of PERP mRNA N6 adenosine methylation promoting the growth and metastasis of pancreatic cancer. Mol Cancer. 2020;19(1):248.

    Article  CAS  Google Scholar 

  33. Yi D, Wang R, Shi X, Xu L, Yilihamu Y, Sang J. METTL14 promotes the migration and invasion of breast cancer cells by modulating N6-methyladenosine and hsa-miR-146a-5p expression. Oncol Rep. 2020;43(5):1375–86.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Zhu H, Gan X, Jiang X, Diao S, Wu H, Hu J. ALKBH5 inhibited autophagy of epithelial ovarian cancer through miR-7 and BCL-2. J Exp Clin Cancer Res. 2019;38(1):21.

    Article  Google Scholar 

  35. Wang JX, et al. Circulating levels of Th1 and Th2 chemokines in patients with ankylosing spondylitis. Cytokine. 2016;2:46.

    Google Scholar 

  36. Wang C, Liao Q, Yihe HU, Zhong DA. T lymphocyte subset imbalances in patients contribute to ankylosing spondylitis. Exp Ther Med. 2015;9(1):250–6.

    Article  CAS  PubMed  Google Scholar 

  37. Lin Q, Jie-Ruo GU, Tian-Wang LI, et al. Value of the peripheral blood B-cells subsets in patients with ankylosing spondylitis. 中华医学杂志 (英文版) . 2009;122(015):1784–1789.

  38. Chen M, Zhang L, Ren Y, Zhang K, Li S. Defective function of CD24(+)CD38(+) regulatory B cells in ankylosing spondylitis. DNA Cell Biol. 2015;35(2):69.

    Google Scholar 

Download references



Author information

Authors and Affiliations



YW designed this topic. ZL analyzed the data and plotted.

Corresponding author

Correspondence to Yansong Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Luan, Z., Wang, Y. Association between ankylosing spondylitis and m6A methylation. J Orthop Surg Res 18, 757 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: