Global DNA methylation patters in peripheral blood and tissues from IgG4-RD patients and controls
DNA methylation profile of B cells and CD4+ T cells from 10 HCs and 10 IgG4-RD patients was analysed. DNA methylation of salivary gland tissues from 4 IgG4-RD patients and controls was also assayed. The workflow was shown in Figure 1A. The overall DNA methylation beta value displayed similar distribution patterns between two groups (Figure 1B). 44 sites were hypomethlated and 166 sites were hypermethylated in B cells and 260 sites were hypomethylated and 112 sites were hypermethylated in CD4+ T cells from IgG4-RD patients compared to healthy controls (Figure 1C). A total of 36945 sites were hypomethylated and 78370 sites were hypermethylated in salivary gland tissues of IgG4-RD patients compared with controls (Figure 1C).
Next, genes and DMPs (Differential Methylation Probes) of salivary gland tissues were obtained and compared with those from B cells and CD4+ T cells in IgG4-RD patients. As presented in Figure 2A, 15 DMPs were common to the three cell types. To further identify the methylated status of common genes in three cell types, the delta beta value were calculated. We found that DPM2 (cg21181453), IQCK (cg10266221), ABCC13 (cg05699681, cg04985582) were hypermethylated and MBP (cg18455083) was hypomethylated in B cells, CD4+ T cells and tissue of IgG4-RD patients compared with controls (Figure 2B). Although TMEM108 (cg18478854) and LOC441897 (cb18346412) were hypermethylated in both B cells and CD4+ T cells, they were found hypomethylated in tissues of IgG4-RD patients compared with controls (Figure 2B).
DNA methylation status in B cells from IgG4-RD patients
Nexted, we analysed CD19+ B cells from 10 IgG4-RD patients and 10 healthy controls and identified 166 hypermethylated and 44 hypomethylated probes (Figure 3A). Among them, the majority of DMPs were present within intergenic regions (IGR) or within gene body (Figure 3A). The top associated hypomethylated and hypermethylated probes with annotated genes are shown in table S2 and table S3. IQCK, UMODL1, AHDC1, LARS2, RASA3, USP16, BTBD11, CLECL1, LOC441897 and ABCC13 genes showed hypermethylated, TMEM9B, CLIC6, MIR492, ARNTL, TNN, MBP, CDK15, PIGL, HLA-DRB1 and MYO1D showed hypomethylated at CpG sites in B cells from IgG4-RD patients.
GO analysis of the DMPs revealed enrichment for multiple categories . Among those, interferon-gamma-mediated signaling pathway (related genes: TRIM68, HLA-DRB1, HLA-DQB2 and IFNGR2), regulation of canonical Wnt signaling pathway (related genes: STK4, BICC1, KANK1, ARNTL, TNN and PSMD5) were shown to be enriched (Figure 3B).
DNA methylation status in CD4+ T cells from IgG4-RD patients
Furthermore, we examined CD4+ B cells from 10 IgG4-RD patients and 10 healthy controls and identified 112 hypermethylated and 260 hypomethylated probes (Figure 3C). Similar with B cells, the majority of DMPs of CD4+ T cells were also present within intergenic regions (IGR) or within gene body (Figure 3C). The top associated hypomethylated and hypermethylated probes with annotated genes are shown in table S4 and table S5. MAST4, AHDC1, BICC1, USP16, IQCK, LOC727677, MYO1D, LOC441897, CAV2 and TGFBR2 genes showed hypermethylated, TMEM9B, CLIC6, CCS, PAWR, FIP1L1, ARHGAP15, ARHGAP26, PKD1L1, HLA-DRB1 and HLA-DQB2 showed hypomethylated at CpG sites in CD4+ T cells from IgG4-RD patients.
GO analysis of the DMPs revealed enrichment for multiple categories . Among those, cellular response to oxidative stress (related genes: FUT8, ARNTL, CFLAR, PTPRK, ETS1, TRAP1, FOXO1, OXR1, PAWR, PRR5L, CCS, PYROXD1 and MIR21), T cell activation (related genes: CLECL1, TGFBR2, ZMIZ1, EOMES, CD3E, TIGIT, CLEC7A, PAWR, KIF13B, CD300A, FZD7 and MIR21) and Wnt-activated receptor activity (related genes: LRP5 and FZD7) were shown to be enriched (Figure 3D)..
Identified DMPs and DMRs in salivary gland tissue of IgG4-RD patients
Furthermore, we studies the DNA methylation status in the affected organ, the salivary gland of the IgG4-RD. A total of 78370 hypermethylated and 36945 hypomethylated probes were identified in tissues of IgG4-RD patients compared with healthy controls (Figure 4A-C). Among them, the majority of DMPs were also present within intergenic regions (IGR) or within gene body (Figure 4C). The top 25 associated hypomethylated and hypermethylated probes with annotated genes are shown in table S6 and table S7. LOH12CR1, PFKP, ARSB, SNX29, TBCD, STARD13, GALNT10, EIF4E3, VGLL4, TGFBR3, GLG1, RBM34, EPAS1, ITPKB, TBS22D1, SLC44A2, TRIOBP, ACAD11, KCNAB2, PDE4D, RERE, ARID1B, UAP1, LGALS3 and LOC102724933 genes showed hypermethylated, CD69, RUNX1, TMEM229B, JAZF1, SLC7A6, FAM69A, EVL, CYTIP, WIPF1, TNFSF8, LOC102724, NCKAP1L, SETBP1, KIAA0748, PPP1R16B, GRAP2, CASP8, ARSG, ITGB2-AS1, MGAT1, CD37, MAP2K2, TNFAIP8L2 and RCSD1 showed hypomethylated at CpG sites in tissues from IgG4-RD patients.
KEGG pathway analysis identified “Calcium signaling pathway” ,“cAMP signaling pathway” ,“MAPK signaling pathway”, “Oxytocin signaling pathway”, “Rap1 signaling pathway”, “Hippo signaling pathway”, “VEGF signaling pathway” and “cGMP-PKG signaling pathway” in the differentially methylated genes in tissues from IgG4-RD patients (Figure 4D). GO analysis of the DMPs revealed enrichment for positive regulation of cell adhesion, T cell activation (Figure 4E).
Finally, we characterized the most significant DNA Methylation Regions (DMRs) in salivary gland of the IgG4-RD patients (table S8). To better characterize genes that were detected in the DMRs, GO and KEGG pathway analyses were performed. KEGG pathway analyses identified genes involved in the cAMP signaling pathway, Calcium signaling pathway, PI3K-Akt signaling pathway, PPAR signaling pathway, Rap1 signaling pathway and Ras signaling pathway (Figure 5A). GO analyses revealed genes were most enriched in extracelluar matrix organization and ion channer activity (Figure 5B).