Epigenome-wide methylation haplotype association analysis identified HLA-DRB1, HLA-DRB5 and HLA-DQB1 as risk factors for rheumatoid arthritis

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.