Causal relationship between the immune cells and ankylosing spondylitis: univariable, bidirectional, and multivariable Mendelian randomization

Background: Ankylosing spondylitis (AS) is an autoimmune disease that affects millions of individuals. Immune cells have been recognized as having a crucial role in the pathogenesis of AS. However, their relationship has not been fully explored. Methods: We chose to employ Mendelian randomization (MR) to investigate the potential correlation between immune cells and AS. We sourced the data on immune cells from the latest genome-wide association studies (GWASs). We obtained data on AS from the FinnGen consortium. Our comprehensive univariable MR analysis covered 731 immune cells to explore its potential causal relationship with AS. The primary analysis method was inverse-variance weighted (IVW). Additionally, we used Cochran's Q test and the MR-Egger intercept test to assess the presence of pleiotropy and heterogeneity. We examined whether our results could be influenced by individual single-nucleotide polymorphisms (SNPs) using the leave-one-out test. We conducted a bidirectional MR to investigate the reverse relationship. We also applied multivariable MR to decrease the potential influence between the immune cells. Results: Overall, our univariable MR analysis revealed eight immune cells associated with AS. Among these, four immune cells contributed to an increased risk of AS, while four immune cells were identified as protective factors for AS. However, the Bonferroni test confirmed only one risk factor and one protective factor with a significance level of p < 6.84E-05. CD8 on effector memory CD8(+) T cell could increase the risk of AS (p: 1.2302E-05, OR: 2.9871, 95%CI: 1.8289-4.8786). HLA DR on CD33(dim) HLA DR+ CD11b(+) could decrease the risk of AS (p: 1.2301E-06, OR: 0.5446, 95%CI: 0.4260-0.6962). We also identified a bidirectional relationship between CD4 on CD39(+) activated CD4 regulatory T cells and AS utilizing the bidirectional MR. To address potential confounding among immune cells, we employed multivariable MR analysis, which revealed that only one immune cell had an independent effect on AS. HLA DR on CD33(dim) HLA DR+ CD11b(+) could decrease the risk of AS (p: 2.113E-06, OR: 0.0.5423, 95%CI: 0.4210-0.6983). Our findings were consistently stable and reliable. Conclusions: Our findings indicated a potential link between immune cells and AS, which could provide a new idea for future research. Nevertheless, the specific underlying mechanisms require further exploration.