Objective

Previous studies have shown that differential DNA methylation is associated with SLE susceptibility. How DNA methylation affects the clinical heterogeneity of SLE has not been fully defined. We conducted this study to identify differentially methylated CpG sites associated with nephritis among women with SLE.

Methods

The methylation status of 428 229 CpG sites across the genome was characterised for peripheral blood cells from 322 women of European descent with SLE, 80 of whom had lupus nephritis, using the Illumina HumanMethylation450 BeadChip. Multivariable linear regression adjusting for population substructure and leucocyte cell proportions was used to identify differentially methylated sites associated with lupus nephritis. The influence of genetic variation on methylation status was investigated using data from a genome-wide association study of SLE. Pathway analyses were used to identify biological processes associated with lupus nephritis.

Results

We identified differential methylation of 19 sites in 18 genomic regions that was associated with nephritis among patients with SLE (false discovery rate q<0.05). Associations for four sites in HIF3A, IFI44 and PRR4 were replicated when examining methylation data derived from CD4+ T cells collected from an independent set of patients with SLE. These associations were not driven by genetic variation within or around the genomic regions. In addition, genes associated with lupus nephritis in a prior genome-wide association study were not differentially methylated in this epigenome-wide study. Pathway analysis indicated that biological processes involving type 1 interferon responses and the development of the immune system were associated with nephritis in patients with SLE.

Conclusions

Differential methylation of genes regulating the response to tissue hypoxia and interferon-mediated signalling is associated with lupus nephritis among women with SLE. These findings have not been identified in genetic studies of lupus nephritis, suggesting that epigenome-wide association studies can help identify the genomic differences that underlie the clinical heterogeneity of SLE.

Objectives

The Xq28 region containing IRAK1 and MECP2 has been identified as a risk locus for systemic lupus erythematosus (SLE) in previous genetic association studies. However, due to the strong linkage disequilibrium between IRAK1 and MECP2, it remains unclear which gene is affected by the underlying causal variant(s) conferring risk of SLE.

Methods

We fine-mapped ≥136 SNPs in a ∼227 kb region on Xq28, containing IRAK1, MECP2 and seven adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15 783 case-control subjects derived from four different ancestral groups.

Results

Multiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at p<5×10(-8) with consistent association in subjects with African ancestry. Of these, six SNPs located in the TMEM187-IRAK1-MECP2 region captured the underlying causal variant(s) residing in a common risk haplotype shared by all four ancestral groups. Among them, rs1059702 best explained the Xq28 association signals in conditional testings and exhibited the strongest p value in transancestral meta-analysis (p(meta )= 1.3×10(-27), OR=1.43), and thus was considered to be the most likely causal variant. The risk allele of rs1059702 results in the amino acid substitution S196F in IRAK1 and had previously been shown to increase NF-κB activity in vitro. We also found that the homozygous risk genotype of rs1059702 was associated with lower mRNA levels of MECP2, but not IRAK1, in SLE patients (p=0.0012) and healthy controls (p=0.0064).

Conclusions

These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.