Abstract
Analyze the relationship between genetic variations in the MTHFR gene at SNPs (rs1801131 and rs1801133) and the therapy outcomes for Iraqi patients with rheumatoid arthritis (RA). The study was conducted on a cohort of 95 RA Iraqi patients. Based on their treatment response, the cohort was divided into two groups: the responder (47 patients) and the nonresponder (48 patients), identified after at least three months of methotrexate (MTX) treatment. A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique was employed to analyze the MTHFR variations, specifically at rs1801133 and rs1801131. Overall, rs1801131 followed both codominant and dominate models, in which in the codominant model, GG [OR (95% CI) 0.11 (0.022-0.553)] and TG [OR (95% CI) 0.106 (0.021-0.528)] predict responders compared to the TT genotype; meanwhile, for the dominate model, the presence of both GG and TG genotypes [OR (95% CI) 0.108 (0.023-0.507)] together predict responders compared to the TT genotype. The A(rs1801133)G(rs1801131) haplotype was significantly associated with responders [OR (95% CI): 0.388 (0.208-0.723)], while the G(rs1801133)T(rs1801131) haplotype was associated marginally with nonresponders [OR (95% CI) 1.980 (0.965-4.064)]. In the final multivariate analysis, GG/TG(rs1801131) genotypes were independently related to responders after adjustment for patients, disease, and treatment characteristics, while TT(rs1801131) genotypes were associated with nonresponders. The Iraqi RA patients showed genetic polymorphism in MTHFR gene rs1801131 with T carrier allele associated with nonresponders to MTX therapy. The rs1801131 followed both codominant and dominant models. The G-carried allele for rs1801131 showed an independent association with responder to MTX therapy after adjustment for patients, disease, and treatment characteristics.