The aim of this study was to compare nonrandom associations between physically adjacent single methylation polymorphism loci among rheumatoid arthritis (RA) and normal subjects for investigating RA-risk methylation haplotypes (meplotype). With 354 ACPA-positive RA patients and 335 normal controls selected from a case-control study based on Swedish population, we conducted the first RA epigenome-wide meplotype association study using our software EWAS2.0, mainly including (i) converted the β value to methylation genotype (menotype) data, (ii) identified methylation disequilibrium (MD) block, (iii) calculated frequent of each meplotypes in MD block and performed case-control association test and (iv) screened for RA-risk meplotypes by odd ratio (OR) and p-values. Ultimately, 545 meplotypes on 334 MD blocks were identified significantly associated with RA (p-value < .05). These meplotypes were mapped to 329 candidate genes related to RA. Subsequently, combined with gene optimization, eight RA-risk meplotypes were identified on three risk genes: HLA-DRB1, HLA-DRB5 and HLA-DQB1. Our results reported the relationship between DNA methylation pattern on HLA-DQB1 and the risk of RA for the first time, demonstrating the co-demethylation of 'cg22984282' and 'cg13423887' on HLA-DQB1 gene (meplotype UU, p-value = 2.90E - 6, OR = 1.68, 95% CI = [1.35, 2.10]) may increase the risk of RA. Our results demonstrates the potential of methylation haplotype analysis to identify RA-related genes from a new perspective and its applicability to the study of other disease.
Keywords: DNA methylation; HLA-DQB1 antigen; HLA-DRB1 chains; HLA-DRB5 chains; disease susceptibility; rheumatoid arthritis.
© 2023 John Wiley & Sons Ltd.