Identification of new overlapping and disease-specific genetic risk factors for rheumatoid arthritis and radiographic axial spondyloarthritis: a meta-analysis of three large European populations and functional characterization

IntroductionThis study conducted a meta-analysis across three large European cohorts (UKBB, FinnGen, and REPAIR), including 12,660 rheumatoid arthritis (RA) cases, 2,446 radiographic axial spondyloarthritis (r-axSpA) cases, and over 530,000 shared controls.MethodsTen independent SNPs in CARMIL1, GRM4, ITPR3, PRSS16, ZNF322, HTT, IKZF1, MANEA, and MGAM2 were analyzed, and functional characterization was performed through cytokine and protein assessments as well as eQTL analyses.ResultsTen independent SNPs were significantly associated with both RA and r-axSpA. Risk alleles included HTTrs363075A, IKZF1rs12718261A, MANEArs72920280T, and MGAM2rs73158426G, while CARMIL1rs72831267C, GRM4rs2495964G, ITPR3rs77601296A, ITPR3rs9469540T, PRSS16rs72843633T, and ZNF322rs6901425G had protective effects.

Functional analysis showed that GRM4rs2495964G was linked to decreased CCL25 levels (p = 0.00030), and ITPR3rs9469540T to reduced IL10 production after LPS stimulation (p = 1.3×10−4). The ZNF322rs6901425G allele was associated with reduced TNFB and increased TGM2 levels (p = 9.60×10−4 and p = 3.00×10−4), both involved in immune signaling and tissue remodeling.

Disease-specific associations were found in BTN2A1, BTN3A2, and H2BC11.