Promoter methylation of the MGAT3 and BACH2 genes correlates with the composition of the immunoglobulin G glycome in inflammatory bowel disease

RESEARCH Open Access Promoter methylation of the MGAT3 and BACH2 genes correlates with the composition of the immunoglobulin G glycome in inflammatory bowel disease Marija Klasić 1† , Dora Markulin 1† , Aleksandar Vojta 1 , Ivana Samaržija 1 , Ivan Biruš 1 , Paula Dobrinić 1 , Nicholas T. Ventham 2 , Irena Trbojević-Akmačić 3 , Mirna Šimurina 3 , Jerko Štambuk 3 , Genadij Razdorov 3 , Nicholas A. Kennedy 2,5 , Jack Satsangi 2,10 , Ana M. Dias 6 , Salome Pinho 6 , Vito Annese 7 , Anna Latiano 8 , Renata D’Inca 9 , IBD consortium, Gordan Lauc 3,4 and Vlatka Zoldoš 1* Abstract Background: Many genomeand epigenome-wide association studies (GWAS and EWAS) and studies of promoter methylation of candidate genes for inflammatory bowel disease (IBD) have demonstrated significant associations between genetic and epigenetic changes and IBD. Independent GWA studies have identified genetic variants in the BACH2,IL6ST,LAMB1,IKZF1,andMGAT3 loci to be associated with both IBD and immunoglobulin G (IgG) glycosylation. Methods: Using bisulfite pyrosequencing, we analyzed CpG methylation in promoter regions of these five genes from peripheral blood of several hundred IBD patients and healthy controls (HCs) from two independent cohorts, respectively. Results: We found significant differences in the methylation levels in the MGAT3 and BACH2 genes between both Crohn’s disease and ulcerative colitis when compared to HC. The same pattern of methylation changes was identified for both genes in CD19 + B cells isolated from the whole blood of a subset of the IBD patients. A correlation analysis was performed between the MGAT3 and BACH2 promoter methylation and individual IgG glycans, measured in the same individuals of the two large cohorts. MGAT3 promoter methylation correlated significantly with galactosylation, sialylation, and bisecting GlcNAc on IgG of the same patients, suggesting that activity of the GnT-III enzyme, encoded by this gene, might be altered in IBD. The correlations between the BACH2 promoter methylation and IgG glycans were less obvious, since BACH2 is not a glycosyltransferase and therefore may affect IgG glycosylation only indirectly. Conclusions: Our results suggest that epigenetic deregulation of key glycosylation genes might lead to an increase in pro-inflammatory properties of IgG in IBD through a decrease in galactosylation and sialylation and an increase of bisecting GlcNAc on digalactosylated glycan structures. Finally, we showed that CpG methylation in the promoter of the MGAT3 gene is altered in CD3 + T cells isolated from inflamed mucosa of patients with ulcerative colitis from a third smaller cohort, for which biopsies were available, suggesting a functional role of this glyco-gene in IBD pathogenesis. * Correspondence: [email protected] † Marija Klasićand Dora Markulin contributed equally to this work. 1 Department of Biology, Division of Molecular Biology, Faculty of Science, University of Zagreb, Horvatovac 102a, 10000 Zagreb, Croatia Full list of author information is available at the end of the article © The Author(s). 2018 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Klasićet al. Clinical Epigenetics (2018) 10:75 https://doi.org/10.1186/s13148-018-0507-y Background Inflammatory bowel disease (IBD) is a chronic intestinal inflammatory condition classified in two major forms— Crohn’s disease (CD) and ulcerative colitis (UC)—which exhibit etiologically and clinically distinct features. Nowadays, IBD affects 2.5–3 million people in Europe and causes considerable morbidity [1]. Despite numerous clinical, genetic, and other experimental studies, our understanding of IBD development and progression remains incomplete. It is generally accepted that IBD represents an aberrant immune response to gut microbiota in genetically susceptible individuals [2]. Genome-wide association studies (GWAS) have identified over 200 genetic susceptibility loci, the majority of which were associated with both forms of IBD in genome-wide meta-analysis [3–7]. However, common genetic variants account only for 8.2 and 13.1% heritability of UC and CD, respectively [7]. Interactionofanindividual’s gut microbiome, immune system, genetic background, and environmental factors, such as smoking, diet, drugs, and physical activity [2,8–10], makes IBD a complex etiopathogenic entity. The challenge is therefore to identify additional factors involved in the development and progression of this disease, especially given its rapidly increasing incidence. It is probable that epigenetics play a key role in the interactions between environmental, microbial, and genetic factors that participate in IBD development and progression. These include DNA methylation and histone modifications, as well as some other epigenetic mechanisms [11–13]; for a review, see [14,15]. DNA methylation remains the most studied epigenetic modification, readily assayed in a large number of individuals/samples. Hypermethylation of gene promoters is generally associated with gene silencing, while promoter hypomethylation is associated with gene activation [16]. Environmentally changed DNA methylation pattern may contribute to the development of many complex diseases by mediating the interplay between external and internal factors and the gene expression [17–21]. There are also data to suggest that the aforementioned environmental modifiers of IBD can also affect DNA methylation [17–19,22]. Epigenetic component of IBD has been addressed in many studies, mostly by whole genome methylation analysis performed on peripheral blood mononuclear cells (PBMCs) or mucosal tissue, revealing regions differentially methylated between the diseaseandhealthystate,aswellasbetweenCDand UC [11–13,23–26]. The majority of eukaryotic proteins are modified by addition of complex oligosaccharides (glycans) through the process of glycosylation. Therefore, glycans are an integral part of nearly all membrane and secreted proteins, including components of the immune system [27]. Aberrant protein glycosylation is implicated in virtually every human complex disease, including inflammation [28–31]. Previous studies have suggested that N-glycosylation of secreted and membrane proteins might be regulated epigenetically and that aberrant glycosylation profiles in disease can arise through aberrant epigenetics [32–38]. A comprehensive review about the role of protein glycosylation in IBD has been given recently [39]. N-glycosylation of serum-circulating proteins (such as the acute phase proteins; immunoglobulin G, IgG; and immunoglobulin A, IgA) or whole plasma N-glycome (i.e., N-glycans present on all plasma proteins) has been the focus of IBD biomarker discovery [36,40–43]. In addition, our partners from IBD consortium and others established that altered glycosylation of IgG, which is a key effector of the humoral immune system, has a role in balancing inflammation at the systemic level [42–46]. GWA studies indicated associations of IBD with several loci involved in protein glycosylation [47,48]. More recently, the first GWAS of IgG glycosylation identified 16 loci specifically associated with changes in IgG glycosylation [49]. Interestingly, five of these loci showed pleiotropy with IBD: MGAT3, a glyco-gene encoding for a glycosyltransferase, GnT-III; LAMB1, a member of transmembrane glycoprotein family of extracellular matrix; the IL6ST, a signal transducer shared by many cytokines; IKZF1; and BACH2, transcription factors involved in B cell differentiation, activation, and maturation. Only the MGAT3 is a classical glyco-gene with a known function in IgG glycosylation, while the exact functional roles for other four GWAS hits in IgG glycosylation or IBD remain unknown. In this study, we investigated promoter methylation differences in these five genes, associated with both IBD and IgG glycosylation, in peripheral whole blood of several hundred IBD patients from two independent cohorts. We also correlated promoter methylation data with IgG glycosylation data analyzed previously for the same IBD patients by our partners from the IBD consortium [43,46,50]. Peripheral blood was used for DNA methylation analysis and serum or plasma was used for glycan analysis, since one of our goals was the search for potential IBD biomarkers. As peripheral whole blood is a heterogeneous cell mixture with specific methylation pattern for each of the cell types [51], we also analyzed promoter methylation of our candidate genes in CD19 + B cells and CD3 + T cells isolated from peripheral blood mononuclear cells (PBMCs). B cells were of our particular interest since these cells produce IgG on their membrane and are precursors of plasma cells which secrete IgG. We have further explored if aberrant promoter methylation recorded in peripheral whole blood of IBD Klasićet al. Clinical Epigenetics (2018) 10:75 Page 2 of 14 patients can be a proxy for epigenetic events occurring in the inflamed mucosa. To address this question, we analyzed DNA methylation from PBMCs, CD3 + T cells isolated from PBMCs, and CD3 + T cells isolated from inflamed colonic mucosa of UC patients from the third smaller cohort, for which biopsies were available. Methods Patient selection and ethics Patients were recruited prospectively from Edinburgh, UK, and Florence, Italy, as a part of the IBD-BIOM project. The recruitment of patients from Edinburgh has been described elsewhere [13,43]. Briefly, we recruited IBD patients prospectively as close as possible to the date of diagnosis from gastroenterology outpatient and endoscopy appointments between 2012 and 2015. We recruited symptomatic controls from gastroenterology clinics during the same period. In these individuals, we had excluded IBD and other organic bowel pathology following biochemical and/or endoscopic investigations. We recruited a further healthy volunteer cohort with no gastrointestinal symptoms. IBD patients were stratified by disease type (ulcerative colitis, UC, and Crohn’s disease, CD). Detailed genetic, phenotypic, and other data regarding IBD cases are given in Additional file 1: Tables S1 and S2. Florence cohort was collected through the network of the Italian Group for IBD (IG-IBD) since the beginning of 2001 and first described in 2005 [1] following an internal validation of phenotyping. Subsequently, longitudinal update has been performed on a yearly basis. Ethical approvals were obtained from Tayside Committee on Medical Ethics B, and all patients and controls provided written, informed consent (LREC 06/S1101/16, LREC 2000/4/192). Florence recruitment details IBD patients were prospectively recruited as close as possible to the date of diagnosis from gastroenterology outpatient and endoscopy appointments between years 2012 and 2015 in different tertiary referral centers in San Giovanni Rotondo, Rome, Rozzano (Milan), Padua, and Florence, Italy. Symptomatic controls were recruited in the same centers (gastroenterology clinics) during the same period. In these individuals, IBD and other organic bowel pathology were excluded by biochemical and/or endoscopic investigations. IBD patients were stratified by disease type (ulcerative colitis, UC, and Crohn’s disease, CD). Samples were obtained with the same methodology (see further) and centrally collected at San Giovanni Rotondo, Italy. Sample collection We collected whole blood at the time of patient recruitment into 9-ml serum Z-clot activator tubes (Greiner), allowed them to clot at 4 °C for 60 min, and then centrifuged at 2500×gfor 15 min. The serum was aliquoted off and stored at −80 °C until further analysis. A subset of patients and controls recruited in Edinburgh (Additional file 1: Table S3) underwent immunomagnetic cell separation to obtain CD19 + B cells. The methods have previously been detailed elsewhere [13]. Venepuncture using 9-ml K3 EDTA vacuette (Greiner) tubes was performed to obtain between of 18 and 36 ml of EDTA-buffered blood. An initial Ficoll (Ficoll-Paque, GE Healthcare, Bucks, UK) density gradient centrifugation was performed to obtain peripheral blood mononuclear cells. Cells labelled with antibody-coated microbeads (human CD8 + and CD19 + microbeads, 20 μlper1×10 7 cells) were immunomagnetic separated using the autoMACs Pro cell separator (Miltenyi, Germany). CD19 + separations were performed following an initial CD8 + depletion step. Nucleic acids were extracted using AllPrep (Qiagen, Hilden, Germany) according to the manufacturer’s guidance and stored at −80 °C. Colonic biopsies from controls and UC patients with inactive and active form of disease were mechanically dissociated to prepare single-cell suspensions using Hanks’balanced salt solution modified medium, without calcium chloride and magnesium sulfate (HBSS) (Sigma), with penicillin/streptomycin and gentamicin. PBMCs were obtained by density gradient centrifugation using Lymphoprep. CD3 + T cells (from biopsies and blood) were magnetically sorted by using the EasySep™Human T Cell Enrichment Kit (STEMCELL) following the manufacturer’s instructions. Following cell isolation, DNA extraction was performed using the Invisorb Spin Tissue Mini Kit (Stratec Molecular) following the manufacturer’s instructions. DNA methylation analysis We analyzed promoter methylation of the candidate genes in the DNA from whole blood, as well as from the separated CD19 + B cells. In addition, for the MGAT3, which is a glycosyltransferase with direct and known function in IgG glycosylation [52], we analyzed promoter methylation—in DNA from PBMCs, CD3 + T cells isolated from PBMCs, and CD3 + T cells isolated from the colonic mucosa of healthy controls and UC patients (classified according to active and inactive form of the disease) of the third independent smaller subcohort collected by the Gastroenterology Department of Centro Hospitalar do Porto-Hospital de Santo António, Portugal (Additional file 1: Table S4). All specimens were subjected to histological examination and classification. All participants gave informed consent about all clinical procedures, and research protocols were approved by the ethics committee of CHP/HSA, Portugal (233/ 12(179-DEFI/177-CES). Klasićet al. Clinical Epigenetics (2018) 10:75 Page 3 of 14 For DNA methylation analysis, 500 ng of DNA from whole blood was bisulfite converted using EZ-96 DNA Methylation Gold kit (Zymo Research, Freiburg, Germany), and 100 ng of DNA from CD19 + B cells, PBMCs, and T cells was converted using EZ DNA Methylation Gold kit (Zymo Research, Freiburg, Germany) according to the manufacturer’s protocol. Two to six pyrosequencing assays were developed for promoter regions of each of the five candidate genes (BACH2,MGAT3,IL6ST,IKZF1, and LAMB1). The selection of analyzed CpG sites was random for assays 2–5 of the MGAT3 gene. CpG sites within the MGAT3 assay 1 were selected based on the GEO (Gene Expression Omnibus) database where methylation data were obtained using Illumina HumanMethylation450 BeadChip v1.1 technology. For the BACH2 gene, assays were selected based on location of differentially methylated CpGs in different cell lines tested by ENCODE project, using Illumina HumanMethylation450 BeadChip v1.1 technology (a newer version, the Infinium MethylationEPIC 850K was not available at the time). We used traditional bisulfite-based protocols which cannot discriminate between 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) as oxidative bisulfite (oxBS-450K) method can [53]. However, recent studies have shown that global DNA hydroxymethylation is very low in blood cells [54,55]. Furthermore, hydroxymethylation is significantly depleted from promotors and CpG islands, while enriched in the gene bodies [53,56]. Based on the estimated statistical power, we did initial screening on 60 patients for each pyrosequencing assay, after which we excluded those genes (pyrosequencing assays) that did not show any statistically significant differences between IBD patients and healthy controls. Pyrosequencing assays for LAMB1,IL6ST, and IKZF1 are shown in Additional file 2: Figure S1. We continued to analyze promoter methylation only in the BACH2 and MGAT3 genes. Specific regions were amplified using PyroMark PCR kit (Qiagen, Hilden, Germany). The cycling conditions for the BACH2 gene were as follows: initial polymerase activation step for 15 min at 95 °C followed by 50 cycles of 30 s denaturation at 95 °C, primer annealing for 30 s at primer-specific temperatures (Additional file 1: Table S5), and 30 s at 72 °C, with final extension at 72 °C for 10 min. The cycling protocol used for amplification of the MGAT3 gene fragments was described previously [35], with the annealing temperature adjusted to 55 °C for the fragment 1 performed on DNA from CD19 + B cells. For quantitative measurement of DNA methylation level at specific CpG sites, PCR-amplified bisulfite-converted DNA was sequenced using the PyroMark Q24 Advanced pyrosequencing system (Qiagen) according to the manufacturer’s recommendations. Sequences of PCR and pyrosequencing primers for the BACH2 and the MGAT3 genes are listed in Additional file 1: Table S5. EpiTect PCR Control DNA Set (methylated and unmethylated bisulfite-converted human DNA, Qiagen) was used as a control for PCR and pyrosequencing reactions. Statistical analysis The nonparametric Mann-Whitney Utest was used to compare the methylation status of CpG sites encompassed by the pyrosequencing assays in the MGAT3 and BACH2 genes between the two independent groups: HC compared to each of CD or UC. Significance threshold was set at p< 0.05 with additional Bonferroni correction for multiple testing. Given that age was our primary concern as a potential confounder, we visualized the age in the three groups (CD, UC, and HC) for the samples included in each analysis as violin plots (Additional file 3: Figure S2) and assured there was no significant difference between the age groups (p> 0.05) using the Mann-Whitney Utest. This was done to assure the validity and strengthen the rationale for the selection of statistical methods. ForthedataoftheMGAT3 promoter methylation from PBMCs, CD3 + T cells isolated from blood, and CD3 + T cells isolated from inflamed colonic mucosa (the Porto cohort), the Mann-Whitney Utest was applied with Bonferroni correction accounting for 15 CpG sites. Glycan analysis Glycans present on IgG were analyzed from serum of over 1000 IBD (UC and CD) patients and healthy controls in the Edinburgh cohort using ultra performance liquid chromatography (UPLC) [43,50]. In the Florence cohort, plasma samples of 3500 IBD patients and healthy controls was used for analysis of IgG glycopeptides by liquid chromatography coupled to mass spectrometry (LC-MS) [46]. The data for IgG glycosylation analysis were used in this work for correlation analysis with promoter methylation data of MGAT3 and BACH2 genes, with matching samples from the very same patients and healthy controls. Isolation of IgG from blood plasma IgG has been isolated from blood plasma by affinity chromatography using CIM Protein G 96-well plate (BIA Separations, Ajdovščina, Slovenia) and vacuum manifold (Pall Corporation, Port Washington, NY, USA) as previously described [57,58]. In short, plasma samples (50–90 μl) were diluted with 1 × PBS, pH 7.4 in the ratio 1:7. All samples were filtered through 0.45 and 0.2-μm AcroPrep GHP filter plates (Pall Corporation) using vacuum manifold and immediately applied to preconditioned Protein G plate. After washing of the Klasićet al. Clinical Epigenetics (2018) 10:75 Page 4 of 14 Protein G plate, IgG was eluted with 0.1 mol L −1 formic acid and immediately neutralized with ammonium bicarbonate to pH 7.0. Protein G plate was regenerated and stored at 4 °C. IgG glycosylation analysis using ultra-performance liquid chromatography N-glycans from isolated IgG in the Edinburgh cohort were released with PNGase F after drying 300 μl of each IgG elution fraction, labeled with 2-aminobenzamide and excess of regents removed by clean-up using hydrophilic interaction liquid chromatography solid phase extraction (HILIC-SPE). Fluorescently labeled and purified N-glycans were separated by HILIC-UPLC using Acquity UPLC instrument (Waters, Milford, MA, USA) as previously described [43]. Samples were separated into 24 peaks [57], and the amount of N-glycans in each chromatographic peak was expressed as a percentage of total integrated area (% area). IgG glycosylation analysis using liquid chromatography coupled to mass spectrometry In the Florence cohort, Fc-specific IgG glycopeptides were analyzed after IgG purification, overnight trypsin digestion at 37 °C, and reverse-phase purification on Chromabond C18 beads using vacuum manifold as described [46,59]. Samples were analyzed using nanoliquid chromatography coupled to mass spectrometry (nanoLC-MS), on a nanoACQUITY UPLC system (Waters, Milford Massachusetts, USA) coupled to quadrupole-TOF-MS (Compact; Bruker Daltonics, Bremen, Germany) equipped with a sheath-flow ESI sprayer (capillary electrophoresis ESI-MS sprayer; Agilent Technologies, Santa Clara, USA) as previously described [46]. The nanoACQUITY UPLC system and the Bruker Compact Q-TOF-MS were operated under HyStar software version 3.2. Data was processed as described previously [46,60]. This resulted in the extraction of 16 IgG1, 16 IgG2/3, and 11 IgG4 glycoforms. The tryptic Fc-glycopeptides for IgG2 and IgG3 subclasses have identical peptide moieties in the Caucasian population and are therefore not distinguishable with this methodology. Annotation of the spectra was done based on accurate mass according to the relevant literature [40,57]. Correlation analysis Methylation data for the BACH2 (assay 2) and the MGAT3 (assays 1 and 2) genes (obtained for the two large cohorts) were filtered according to the peak quality by rejecting peaks marked as “failed”by the pyrosequencing software. Average methylation across all assayed CpG sites was calculated for each pyrosequencing assay in each cohort. Methylation results for individual patients were matched with their corresponding glycan profiles. Sizes of datasets and patient classes obtained after including complete records (i.e., both methylation data and glycan profiles present) are shown in Additional file 1: Table S6. Individual glycan structures were represented as relative abundances and batch-corrected. Percentage of structures with bisecting N-acetylglucosamine was calculated for each cohort as a derived trait at this point. Glycan structures identified by each method were translated to Oxford notation, and only the 13 structures present in both the Edinburgh and Florence datasets were considered for correlation. We used IgG1 data from the Florence cohort, as this isoform was the most abundant. Three additional derived traits were calculated: ratios of FA2B to FA2, FA2BG1 to FA2G1, and FA2BG2 to FA2G2. Pearson correlation between CpG methylation data and 17 glycan features (13 structures and 4 derived traits) was calculated. Significance threshold was set at p< 0.05 with additional Bonferroni correction for 17-fold multiple testing. Methylation of assayed CpG sites in promoters of the BACH2 and MGAT3 genes was correlated with measured glycan structures. Pearson correlation coefficient along with the associated pvalue was calculated between average CpG methylation (for all genes/assays) and each measured IgG glycan structure. Calculation was done on pairwise complete observations. Only correlations with the pvalue below 0.01 were considered further. Next, correlation coefficients for all CpG assays were calculated, which was used to rank glycan structures according to regulation by the assayed region. Glycan structures with the strongest correlation (either positive or negative) to CpG methylation were then used to explain regulatory effects. All calculations and data visualizations were done in R language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria). Visualization of correlations was done using the R package “corrplot.” Results Promoter methylation of the candidate genes in whole blood and B cells of IBD patients In order to assess the level of methylation in CpG islands of the five candidate genes (BACH2,MGAT3, IKZF1,LAMB1,andIL6ST), associated with both IBD and IgG glycosylation by GWAS, we developed several pyrosequencing assays for each of the genes (Fig. 1and Additional file 2:FigureS1).Weperformed initial screening of the pyrosequencing assays on 60 patients. Overall cytosine methylation levels were very low for LAMB1 (average value per group < 8%), for IL6ST (< 3.5%), and for IKZF1 (< 4%) in the assayed portion of their promoters; therefore, we Klasićet al. Clinical Epigenetics (2018) 10:75 Page 5 of 14 could not identify differential methylation. We then excluded these genes from further analysis. MGAT3 and BACH2 promoter methylation was analyzed in several hundred IBD patients and healthy controls from two independent cohorts (Additional file 1: Tables S1 and S2). In these genes, we analyzed methylation level at 47 CpGs covered by five pyrosequencing assays in the BACH2 gene: 21 CpG sites were in the promoter region, 1 CpG site was in first exon, and 25 CpG sites were located in the first intron of the gene. A total of 32 CpG sites, covered by five pyrosequencing assays, was analyzed for MGAT3: 18 CpG sites were located in the promoter region and 14 CpG sites in the first intron (Fig. 1). Most of those CpG sites were located within CpG islands of the both genes. We found differential CpG methylation between IBD patients and HC within the assay A2, located at 213–368 bp upstream (relative to the gene orientation) of the TSS in the BACH2 promoter and within the assays A1 and A2, located in the CpG island 1 of the MGAT3 gene. The same pattern of differential methylation at these CpG sites was observed in whole blood of patients and HC from two large independent cohorts (Fig. 2). CpG methylation level was generally low (up to 20%) in the assayed portion of the BACH2 promoter; however, significant differences between HC and CD methylation level were recorded at CpG sites 4, 5, 6, and 8 (Fig. 2a). For the assayed portion of the MGAT3 promoter, general methylation level was high, with all CpG sites showing a reproducibly significant difference between HC andbothCDandUC.CpGsites2,13,and15 showed significant differences only for CD but not for UC. Direction of change was different for the two genes—differentially methylated CpG sites within the BACH2 promoter were hypomethylated, while those for the MGAT3 gene were hypermethylated in disease compared to healthy individuals. These results were confirmed on CD19 + Bcells isolated from peripheral whole blood of the independent, smaller patient sample from the Edinburgh cohort (67 samples). The CpG sites 1–5 and 12–13 in the MGAT3 promoter were differentially methylated between CD and HC. Only CpG site 5 within the assay A2 of the BACH2 gene showed change in the methylation level between HC and CD in CD19 + B cells (Fig. 2b). There were no differences in the methylation level of the same CpG sites within assayed fragments of the BACH2 and MGAT3 genes between UC and HC. It is worth noting that the same pattern of CpG methylation differences was observed in PBMCs of the IBD patients and HC from both large independent cohorts, and most of the CpG sites within the assayed portion of the MGAT3 promoter were also differentially methylated in CD19 + B cells from the subset of IBD patients from the Edinburgh cohort (Fig. 2a,b). Fig. 1 Positions of the BACH2 and MGAT3 genes in the human genome and relative positions of the fragments analyzed for methylation level (pyrosequencing assays) within these genes. For each pyrosequencing assay (A1–A5), the region amplified by PCR is shown. Positions of the genes on the chromosomes are shown using chromosome models (red vertical lines). Coordinates are relative to the hg19 human genome assembly. The genes are displayed in the direction corresponding to their reading frames. Annotations (CpG islands and pyrosequencing assays) are to scale. TSS transcription start site Klasićet al. Clinical Epigenetics (2018) 10:75 Page 6 of 14 a b Fig. 2 (See legend on next page.) Klasićet al. Clinical Epigenetics (2018) 10:75 Page 7 of 14 Promoter methylation of the MGAT3 gene in CD3 + T cells from PBMCs and inflamed colonic mucosa of UC patients We included in our investigation biopsy samples of UC patients from an independent cohort from the Gastroenterology Department of Centro Hospitalar do Porto-Hospital de Santo António, Portugal. Given the technical challenges in obtaining DNA and RNA from a small number of purified cells from inflamed colonic mucosa, a subset of patients with active and inactive phase of UC was selected for methylation analysis from three sources: (1) PBMCs, (2) CD3 + T cells isolated from PBMCs, and (3) CD3 + T cells isolated from colonic mucosa (see also Additional file 1: Table S4). Inter-individual variation of MGAT3 methylation level measured from PBMCs and from CD3 + T cells isolated from PBMCs was quite large—it varied from 47 to 94% and from 26 to 90%, respectively. Therefore, we could not find any difference in CpG methylation level between UC patients and HC in assayed fragments of the MGAT3 promoter, neither in PBMCs nor in CD3 + T cells isolated from PBMCs. However, we recorded a total of 7 (out of 15) differentially methylated CpG sites in CD3 + T cells isolated from colonic mucosa of UC patients with active disease compared with HC (Fig. 3). Overall, the methylation level of CpGs within assayed fragments of the MGAT3 promoter was high in CD3 + T cells from healthy colonic mucosa (between 77 and 98%). When compared to inflamed mucosa of UC patients with active phase of the disease, the same CpG sites were hypomethylated, with the highest difference at the CpG position 10 (13.24%; p= 5.08 × 10 −5 ; Fig. 3). In inactive UC, no significant differences could be found after Bonferroni correction for multiple testing. It is worth noting that the methylation pattern in CD3 + T cells isolated from inflamed colonic mucosa differed from the methylation patterns in PBMCs and for CD3 + T cells isolated from PBMCs. The latter two were very similar and had much lower methylation levels than that measured for CD3 + T cells from inflamed colonic mucosa (Fig. 3). Also, MGAT3 methylation level was increased in UC compared with HC when measured from PBMCs or CD3 + T cells isolated from PBMCs (hypermethylation), while it was decreased (hypomethylation) when measured from CD3 + T cells isolated from inflamed colonic mucosa in comparison with healthy mucosa. Correlation between the MGAT3 and BACH2 promoter methylation and IgG glycosylation There were statistically significant correlations that replicated across assays and cohorts between the MGAT3 promoter methylation and glycan structures FA2, FA2G2, FA2BG2, and FA2G2S1, as well as the derived trait of the ratio of FA2B to FA2 (Fig. 4a). All correlations except with FA2 were negative. No reproducible significant correlations could be found between BACH2 promoter methylation and the glycan structures (Fig. 4a). In order to infer a mechanistic pathway of the observed correlations, we mapped them to the glycan biosynthesis pathways (Fig. 4b). The ratio of bisecting glycans to FA2 was taken as an indicator of MGAT3 (GnT-III) activity. This interpretation allowed us to infer lower GnT-III enzymatic activity when the promoter of the MGAT3 gene was methylated. Increase in MGAT3 promoter methylation correlated with a decrease in certain galactosylated and sialylated structures (Fig. 4b). In addition to the decreased levels of bisecting GlcNAc on non-galactosylated glycans (B/FA2), the most significant effect of the MGAT3 promoter methylation on IgG glycome composition was a decrease of IgG galactosylation. Discussion Results from this study strongly indicate that the MGAT3 and BACH2 genes play an important role in IBD pathogenesis and suggest a possible disease pathway mediated by the pro-inflammatory properties of IgG antibodies acquired by alterations in Fc glycosylation. Our recent study, performed on a large cohort of over 1000 IBD patients, reported a significant difference in IgG glycome composition in both UC and CD compared to healthy controls [43,46]. We found a decrease in quantity of galacosylated glycans in both CD and UC, as well as a decrease in sialylated glycans and an increase of bisecting GlcNAc on digalactosylated glycan structures on IgG in CD. Indeed, alternative N-glycosylation of an IgG molecule influences its function—pro-inflammatory and anti-inflammatory activity depends on the (See figure on previous page.) Fig. 2 Box plot of CpG methylation in peripheral whole blood for the BACH2 and MGAT3 genes in the Edinburgh and Florence cohorts and in B cells from a subset of patients from Edinburgh cohort. Groups were compared using the Mann-Whitney Utest with significance threshold of p= 0.05, corrected for multiple testing using the Bonferroni method. aMethylation levels were generally low in the assayed portion of the BACH2 gene promoter, with significant differences between HC and CD methylation at CpG sites 4, 5, 6, and 8 (replicated in both cohorts). For the MGAT3 gene, general methylation level was high, with all CpG sites showing a reproducibly significant difference between HC and both CD and UC, except for CpG sites 2, 13, and 15 for which reproducible significant differences were found only between HC and CD. bIn B cells, isolated from PBMCs of a subset of the patients from the Edinburg cohort, differential methylation was found at the CpG position 5 of the BACH2 gene (assay 2) between HC and CD, while for the MGAT3 gene, differentially methylated were CpG sites 1–5, 12, and 13 between HC and CD. CD Crohn’s disease, UC ulcerative colitis, HC healthy controls Klasićet al. Clinical Epigenetics (2018) 10:75 Page 8 of 14 glycans added on the Cy2 domain of its Fc region [29]. These glycans are of a biantennary complex type with or without bisecting GlcNAc, core fucose, galactose, and sialic acid residues [61]. Recently, this was confirmed in a large multi-centric study of IgG glycome in IBD [46]. Therefore, glycan changes observed on IgG in peripheral blood of UC and CD patients are obviously associated with increased inflammatory potential of IgG, suggesting functional relevance of IgG glycosylation for IBD. Here, we propose a possible mechanism underlying the aberrant IgG glycosylation pattern observed in IBD [43,46]. Out of five candidate genes analyzed in this work, the MGAT3, a glycosyltransferase which participates in synthesis of IgG glycans, and the BACH2,a transcription factor and a master regulator of a network of genes relevant for B cell integrity [62,63], showed differential methylation in peripheral blood of both CD and UC patients when compared to healthy individuals. Even though we identified changes in methylation level for both UC and CD compared to HC, the differences were more pronounced for CD. This is concordant with other studies that explored either whole genome methylation or promoter methylation of candidate genes in IBD [11]. The extent of the change in IgG glycome composition was also consistently higher in CD than UC compared to HC [43,46]. The protein encoded by the MGAT3 gene (N-acetylglucosaminyltransferase III, GnT-III) is responsible for significant functional alteration of glycans on the Fc region of an IgG antibody. The GnT-III adds N-acetylglucosamine (GlcNAc) on β1,4-linked mannose in the three-mannose core of N-glycans, producing bisecting GlcNAc structures. In the same CD patients, who showed changed MGAT3 promoter methylation level in peripheral blood cells, a significant increase in the percentage of bisecting GlcNAc on glycans of circulating IgG antibodies was recorded, too. The association of the MGAT3 with both IgG N-glycosylation [49] and Crohn’s disease [4,5] suggests that N-glycans with bisecting GlcNAc could be involved in CD pathogenesis through functional effect on IgG antibody. Correlations between BACH2 and MGAT3 promoter methylation and glycan structures have given further insight into the changes of IgG glycosylation pattern mediated by those two genes (Fig. 4b). The MGAT3 a b c Fig. 3 Box plot of CpG methylation level in the MGAT3 gene promoter (assays A1 and A2) analyzed from PBMCs (a), CD3 + T cells isolated from PBMCs (b), and CD3 + T cells isolated from inflamed colonic mucosa (c) from the independent cohort of Porto. Changes between UC patients with active disease and HC were statistically significant only in CD3 + T cells isolated from inflamed colonic mucosa at CpG positions 3 and 7–12 (p< 0.05 after Bonferroni correction for 15 hypotheses). PBMC peripheral blood mononuclear cells, UC ulcerative colitis, HC healthy controls Klasićet al. Clinical Epigenetics (2018) 10:75 Page 9 of 14