Long-term Temporal Stability of Peripheral Blood DNA Methylation Profiles in Patients With Inflammatory Bowel Disease

Background & Aims There is great current interest in the potential application of DNA methylation alterations in peripheral blood leukocytes (PBLs) as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). However, the intra-individual stability of PBL methylation in IBD has not been characterized. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array. Methods We followed a cohort of 46 adult patients with IBD (36 Crohn’s disease [CD], 10 ulcerative colitis [UC]; median age, 44 years; interquartile range [IQR] 27–56 years; 50% female) that received standard care follow-up at the Amsterdam University Medical Centers. Paired PBL samples were collected at 2 time points with a median of 7 years (range, 2–9 years) in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively. Results Around 60% of all EPIC array loci presented poor intra-individual stability (ICC <0.50); 78.114 (≈9%) showed good (ICC, 0.75–0.89), and 41.274 (≈5%) showed excellent (ICC ≥0.90) stability, between both measured time points. Focusing on previously identified consistently differentially methylated positions indicated that 22 CD-, 11 UC-, and 24 IBD-associated loci demonstrated high stability (ICC ≥0.75) over time; of these, we observed a marked stability of CpG loci associated to the HLA genes. Conclusions Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBD-associated epigenetic biomarkers with increased potential for independent validation. These data also have potential implications in understanding disease pathogenesis.


BACKGROUND & AIMS:
There is great current interest in the potential application of DNA methylation alterations in peripheral blood leukocytes (PBLs) as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). However, the intra-individual stability of PBL methylation in IBD has not been characterized. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array.

METHODS:
We followed a cohort of 46 adult patients with IBD (36 Crohn's disease [CD], 10 ulcerative colitis [UC]; median age, 44 years; interquartile range [IQR] 27-56 years; 50% female) that received standard care follow-up at the Amsterdam University Medical Centers. Paired PBL samples were collected at 2 time points with a median of 7 years (range, 2-9 years) in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively.
CONCLUSIONS: Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBDassociated epigenetic biomarkers with increased potential for independent validation. These data also have potential implications in understanding disease pathogenesis. C rohn's disease (CD) and ulcerative colitis (UC) are chronic relapsing and remitting inflammatory bowel diseases (IBDs) characterized by a wide variety of phenotypic manifestations. 1 Although the etiology of IBD remains unknown, it is thought to arise as a result of a complex interplay between the host and microbial composition, triggered by environmental factors, such as tobacco smoking or diet. [2][3][4] Accordingly, much effort has been invested in understanding the interaction between host and environment, which is thought to be mediated by the epigenome. 5 The epigenome represents the set of mitotically heritable modifications that can affect gene transcriptions without altering the primary DNA sequence. 6 DNA methylation, one of the most studied epigenetic mechanisms, involves the attachment of methyl groups to cytosine-phosphate-guanine (CpG) nucleotide sequences on the DNA. This covalent attachment is mitotically heritable and can, under certain conditions, regulate gene expression, thereby altering cellular behavior. 7 Over the past decade, multiple epigenome-wide association studies (EWAS) have sought to characterize, classify, and predict IBD and its various phenotypes using DNA methylation. [8][9][10][11][12][13] However, most EWAS in IBD to date have been cross-sectional in design, reporting aberrant DNA methylation signatures in peripheral blood leukocytes (PBLs) and/or mucosal tissue, 9,11,12,[14][15][16][17][18][19][20] with only a single longitudinal study in mucosal tissue 12 and PBLs. 21 Previous literature has shown that the intra-individual variability of DNA methylation is most prominent during the early stages of life, which gradually diminishes and presents a more stable phenotype after 5 years of life. 22,23 Nonetheless, the influence of aging on genome-wide DNA methylation has been well-described in monozygotic twins 24,25 and unrelated healthy populations, 26,27 demonstrating a global decrease in methylation as individuals age, as well as site-specific increases in methylation in CpG-rich areas, both of which are thought to result from dynamic external and internal environmental changes. [28][29][30] As epigenetics, and thus DNA methylation, is cell-type specific, observed differences found in heterogeneous populations such as PBLs or tissue might reflect differences in the cellular composition. [31][32][33] Nonetheless, age-related differences were found in more homogeneous populations, such as purified T-cells and monocytes. 15,34,35 Despite the strong effects of age on DNA methylation, a high correlation between baseline and follow-up methylation data in pediatric IBD mucosal tissue has been observed. 12 In contrast, IBD-associated differences in blood have shown to largely revert back to patterns observed in non-IBD controls during follow-up as the result of treatment and normalization of C-reactive protein (CRP). 21 It is noteworthy that these studies focused only on a subset of IBD-associated CpGs, and did not report on the long-term stability of all CpG probes located on the Illumina Human-Methylation EPIC BeadChip array. Although temporal stability and intra-individual variability in PBL-derived DNA methylation has been investigated in adult healthy individuals [36][37][38][39] and patients with systemic lupus erythematosus, 40 no such study has been conducted in patients with IBD. There is widespread interest in the application of epigenetic markers in personalization of treatment. 41 If epigenetic biomarkers are to be used as pathognomonic for IBD or its (sub)phenotypes, the features of interest would need to remain stable throughout the duration of the disease and, hence, over time, without being affected by various internal/external exposures. In addition, for biomarker development, loci that are time-stable reduce the number of false-positive findings, thereby increasing the probability of independent replication. Furthermore, selection of time-stable epigenetic biomarkers would help overcome current practical barriers in sample collection at specific time frames, thereby facilitating the use of larger samples sizes with similar phenotypes needed to enhance predictive power. We therefore sought to identify CpG positions that present stable DNA methylation in PBLs obtained from a well-generalizable cohort of adult patients with IBD with a 7-year median time between collected DNA samples.

Patient Demographics
A total of 46 adult patients with IBD (36 CD, 10 UC) with a median age of 44 years (interquartile range [IQR], 27-56 years) and median disease duration of 12 years (IQR, 7-21 years) were included. Gender, surgical history, disease location, and disease behavior were balanced within this cohort (Table 1). Notably, 32 patients (69.6%) were previously treated with an anti-tumor necrosis factor, prior to T1 sampling. Between T1 and T2 during regular IBD care follow-up, 10 patients (21.7%) underwent IBD-related surgery, 24 patients (52.2%) were treated with vedolizumab, and 14 patients (30.4%) were treated with ustekinumab, reflecting the tertiary referral population seen at the Amsterdam University Medical Centers (UMC). No significant differences in median CRP (P ¼ .97) or leukocyte count (P ¼ .85) between T1 and T2 were observed (Figure 1, B-C).

Time-associated Differential Methylation Expectedly Associates With Age-related CpGs
We first investigated the differences in methylation between both time points, identifying 194,391 (z23%) differentially methylated positions (DMPs) when comparing T1 and T2 at a false-discovery rate (FDR)-adjusted P-value of < .05 ( Figure 2, A-C), which we termed time-associated DMPs. As our sample of interest was derived from peripheral blood, we investigated whether differences in the cellular composition were observable. Comparing the predicted blood cell composition yielded significant increase of the B-(P ¼ .017) and CD4þ T-cells (P ¼ .013) over time, whereas the neutrophils present a significant decrease (P ¼ .02) (Figure 2, D).
Expectedly, the time-associated differences were enriched for age-associated CpGs, which have been defined as the "epigenetic clocks" from Horvath, 42 Hannum, 43 Levine, 44 and Knight 45 (Figure 2, E). Furthermore, for these specific epigenetic clock CpGs, we observed a general hypomethylated pattern at T2 relative to T1 CpG sites ( Figure 2, F), suggesting that the observed differences in DNA methylation are enriched for age-related differences. Functional enrichment analyses of the time-associated DMPs displayed several cancer-associated pathways ( Figure 3).

Time-invariant, Stable Methylated Probes are Enriched in Genes Involved in Cell Adhesion
To identify CpGs that were consistently methylated at both time points, we performed intra-class correlation (ICC) analysis, which indicated that the majority of the CpGs (517.576 probes or around 60%) present poor intraindividual stability over time (ICC <0.50) ( Table 2). Conversely, 119.388 CpGs (z14%) displayed a statistically significant high ICC (!0.75), which we termed stably methylated positions (SMPs). Expectedly, CpGs with high ICC values typically presented less difference in mean methylation ( Figure 4). We reasoned that probes that were associated with sites known to harbor genetic variants, both intentional and unintentional, 46 should present the highest stability, as the genome of an individual typically does not change over the course of 7 years. Indeed, splitting the data by modality suggested that the CpG sites associated with known germline variants, namely those that were included   for quality control purposes, presented high (>0.9) ICC values ( Figure 5, A and Table 2). Of the SMPs with ICC values over 0.9, 15.766 SMPs (around 2%) presented no indication that they bind predicted or potential genetic variants, which we classified as hyper-stable methylated positions (HSMPs) ( Table 3 and Figure 5

Stability Analysis of Previous IBD-associated DMPs, HLA, and IBD-susceptibility genes
We next investigated whether previously reported IBDassociated DMPs were found to be invariant over time. To do so, we evaluated ICC values of 255 CD-associated, 103 UC-associated, and 221 IBD-associated consistent DMPs identified in our systematic review and meta-analysis on CD-, UC-, and IBD-associated differential methylation, 48 which included a total of 552 samples (177 CD, 132 UC, and 243 HC) from 4 different EWAS. 14,[17][18][19] Focusing on the stability of these DMPs in this cohort, we show that the majority (151 or 59.2% CD-associated, 73 or 70.9% UCassociated, and 156 or 70.6% IBD-associated) present poor to moderate stability, indicating that the methylation status of these DMPs are affected by age or other exposures over time (Figure 8, A; Supplementary Table 1). Nonetheless, 22 CD-associated (12.4%), 11 UC-associated (8.3%), and 24 IBD-associated (9.9%) loci show good to excellent stability over time, providing evidence that these CD-, UC-, or IBD-associated DMPs are unaffected by aging or the exposures over time (Figure 8, A; Supplementary Table 1).
Among the many IBD-associated DMPs, we specifically zoomed in on VMP1 (cg12054453 and cg16936953) as well as RPS6KA2 (cg17501210), as they were identified in multiple IBD-EWAS, 17,18,21 as well as shown found to be among the most significant IBD-associated DMPs in our metaanalysis. 48 We observed moderate consistency over time, with noticeable overall hypermethylation at T2 relative to T1 for the aforementioned 3 CpGs (Figure 8, B).
In addition to our meta-analysis, we also interrogated the CpGs that were CD-associated but CRP-independent, as reported on by Somineni et al. 21 Interrogation thereof using our cohort revealed that cg25112191 (RORC), cg13707793 (CXXC5), cg21049840 (GPR183), cg06460200 (GPR183), cg06366627 (DIDO1), and cg15860510 (ESP8L3) presented poor consistency (ICC <0.5); cg04570316 (GMNN), cg02240291 (SMARCD3), cg00092736 (ESPNL), and cg00092736 (ESPNL) presented moderate consistency (ICC, Note: An overview of the ICC values classified using the system presented by Koo and Li. 47 QC GV ¼ Quality control probes that bind genetic variant only. Annotated GV ¼ Methylation probes that are annotated to bind genetic variants at the CpG of interest. Predicted GV ¼ Methylation probes that are annotated to bind genetic variants that were annotated by Gaphunter to be caused by genetic variants. Methylation ¼ Methylation probes for which we have no evidence that they bind genetic variants. GV, Genetic variants; ICC, interclass correlation. We next interrogated the stability of all CpG loci associated with several well-known GWAS-identified IBD risk genes involved in IBD pathogenesis, namely ATG16L1, NOD2, IL23R, CARD9, FUT2, TYK2, and TNFSF15, 49,50 as well as specific IBD-associated major histocompatibility complex encoding HLA genes previously reported on in GWAS studies, namely HLA-DRB1, HLA-DQB1, HLA-DQA1, HLA-DPA1, HLA-DPB1, HLA-A, HLA-B, and HLA-C. 49,[51][52][53][54] Comparing all IBD risk genes, we noticed that the HLA genes presented the highest stability, all of which had a median ICC score over 0.5, whereas the majority of CpGs that annotate to non-HLA IBD risk genes had poor ICC values (<0.5) (Figure 9, A). Nonetheless, for each of these non-HLA IBD risk genes, we identified highly stable methylated positions, several of which located to transcription start sites or first exons ( Figure 10 and Supplementary  Table 2), implicating potential regulatory function.
As DNA methylation measured using deamination technologies cannot distinguish DNA methylation from genetic variants located at the CpG of interest, 46 we investigated whether such technical artefacts were found among the HLA SMPs by interrogating the dbSNP (v151) database for catalogued variants, as well as by investigating for a typical clustered methylation signal when probes hybridize with genetic variants (GVs) using Gaphunter. 55 Notably, most of the high ICC values were found for CpGs that presented some type of clustering typical of GVs but were not necessarily catalogued in dbSNP (v151) (Figure 9, B). In addition, HLA class II genes (HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, and HLA-DRB1) appeared to have a larger proportion of highly stable probes compared with HLA class I genes (HLA-A, HLA-B, and HLA-C) (Figure 9, B). Besides technical artefacts, DNA methylation itself can be affected by GVs that occur in the vicinity. 56 As such, we cross-referenced our observations with a previous large-scale genotyping study of the HLA region in both patients with CD and UC. 52 We indeed found multiple probes within the vicinity (<1 Kb) of CD-or UCassociated HLA alleles, many of which did not appear to have an annotated dbSNP identifier, indicating that the observed differences were either unannotated GVs, or CpGs whose methylation status is strictly controlled by neighboring IBD-associated GVs. Notably, several highly stable probes (ICC !0.75) found to be annotated to HLA-C, HLA-DPB1, and HLA-DPA1 were located far away (>1 Kb) from any of the IBDassociated GVs, did not associate with catalogued GVs, nor were identified as a potential GV by Gaphunter (Figure 9, C-E; Figure 11; and Supplementary Table 3), suggesting by default strong methylation stability over time.

Discussion
Biomarker research often involves samples taken prior to or within a strictly pre-defined timeframe to the outcome of interest to mitigate the number of additional variables. In this study, we performed long-term longitudinal stability analyses of the PBL DNA methylome obtained from a cohort of adult patients with IBD (36 CD and 10 UC) that were collected at 2 time points separated by a median of 7 years reflective of a real-life tertiary referral population.  CGID, Illumina CpG identifier; Coordinate, Genomic coordinate of the CpG on the human genome (build hg19); DM, mean difference in percentage methylation; ICC, intraclass correclation coefficient; Gene, associated gene as well as the location in the gene; P-value ICC , P-value associated with the intraclass correlation coefficient; Padjusted ICC , Benjamini-Hochberg-adjusted P-value associated with the intraclass correlation coefficient; P-value DM , P-value associated with the mean difference in percentage methylation; P-adjusted DM , Benjamini-Hochberg-adjusted P-value associated with the mean difference in percentage methylation.
Our observations indicate that the majority of all loci (z60%) measured on the Illumina HumanMethylation EPIC BeadChip array present notable intra-individual variation in methylation over time (ICC <0.5), which is enriched for age-associated CpGs. Nonetheless, not all time-associated DMPs in our cohort were previously reported as age-associated CpGs. Although we observed no significant differences in CRP and leukocyte count between both time points, other external or environmental exposures, such as smoking, dietary alterations, therapy failure or switch, IBD-related surgery, or disease progression might have altered methylation status contributing to the observed time-associated differences. By contrast, 119.388 (z14%) and 41.274 (z5%) loci presented a highly stable pattern across both time points, with ICC values !0.75 and !0.90, respectively. Such loci retained their degree of methylation even after the aforementioned known IBD-associated and unknown external exposures as well as differences in cellular composition, suggesting time stability. Previous studies investigating DNA methylation stability in PBLs of healthy adults using both the Illumina Human-Methylation EPIC BeadChip array, as well as its predecessor the Illumina HumanMethylation 450k BeadChip array, presented similar observations. In these studies, 16.9% to 23% of the CpG loci presented a moderate/good (ICC, 0.50-0.79), whereas 8.3% to 12.9% of the CpG loci presented a good/ excellent (ICC !0.8) stability over a span of 1 to 6 years. 37,39 Focusing on IBD, when interrogating 255 CD-associated, 103 UC-associated, and 221 IBD-associated DMPs identified in our own meta-analysis 48 of 4 IBD EWAS, 14,[17][18][19] we observed that the majority presented poor to moderate stability, suggesting that the aforementioned IBD-associated loci might also be affected by exposures over time that might or might not be related with IBD. Although interesting, such time-variant probes should be interpreted with care when used as predictive biomarkers, given their association with exposures that occurred during both time points. To that end, our data could be used as a resource to preselect time-invariant CpG loci before independent validation when performing IBDassociated EWAS, 57 thereby increasing the potential to identify replicable predictive biomarkers better reflecting the underlying biology of IBD. In addition, such an approach would enable a larger pool of samples to be used as samples need not to be obtained within the same age range when performing DNA methylation studies on IBD and its phenotypes.
When specifically interrogating the IBD-associated probes cg12054453 (VMP1), cg16936953 (VMP1), and cg17501210 (RPS6KA2), we find moderate consistency over time with a noticeable hypermethylation at T2 compared with T1. Similar analyses performed on the CD-associated yet CRP-independent probes reported by Somineni et al 21 showed that the majority of these CpGs did not present long-term stability in our cohort. Differences between our observations and that of Somineni might be attributable to differences between adult and paediatric cohorts or might simply reflect non-inflammatory changes in methylation that occur over time. By contrast, cg09171692 (RORC) and cg02055816 (SMARCD3) presented high ICC values (0.76 and 0.87, respectively), indicating good stability in our cohort, irrespective of CRP or noninflammatory exposures. Notably, both genes have previously been associated with IBD in multiple studies 21,58,59 and therefore, show promise as stable IBD-associated loci.
Given the complex, multifactorial nature of IBD, 2 focussing on the interplay between genetic variation and DNA methylation rather than single gene mutations alone might prove more useful in understanding its molecular etiology. Previous EWAS of mucosal tissue 60 and peripheral blood 61 both demonstrated differential methylation between IBD and controls for well-known GWAS-identified IBD-associated risk genes, suggesting differential methylation of key risk genes to affect disease susceptibility. Our observations corroborate this hypothesis, showing highly stable methylation for particular CpG loci within these risk genes. Interestingly, several of these HSMPs were located in or near to the transcription start sites, potentially regulating gene transcription by maintaining the aberrant phenotype. 18 There has been extensive interest in (epi)genetic alterations of the highly polymorphic HLA region related to IBD pathogenesis, most consistently reported for HLA class II genes involved in the presentation of bacterial antigens to CD4þ Tcells. 14,18,49,[51][52][53][54]62,63 Specifically, genetic variation of classical HLA genes has been suggested to play a role in the aberrant response to the dysbiotic microbiome observed in IBD, 51 with particular impacts for the response to biological treatment [64][65][66] and the formation of anti-drug antibodies. 54 However, translation of the results into clinical practice has proven to be difficult due to the high number of polymorphisms of HLA ab heterodimers and strong linkage disequilibrium. 51 In our study, we observed multiple HSMPs in HLA genes, suggesting that the DNA methylation profile of these genes is very stable over time. Our results corroborate with previous array data showing highly significant correlations between CpG loci on several HLA class II genes of neonates compared with toddlers (r ¼ 0.83) and adults (r ¼ 0.88) with type 1 diabetes. 67 Although further interrogation of these HSMPs indicates that multiple CpGs might be actual genetic variants, we also find multiple HSMPs that are not genetic variants. Nonetheless, several of such epigenetic HSMPs do occur within the vicinity of known IBD-associated HLA-variants, providing evidence that particular HLA-alleles might impart a strong, stabilizing effect on the epigenome.

Strengths and Limitations
To our knowledge, we are the first to assess the stability of the DNA methylome obtained from PBLs of patients with IBD with a median 7-year follow-up period in a real-life disease exposure setting. This study is explorative in nature, using a moderate sample size without prior power calculation. Nonetheless, we note that studies of similar design have been conducted with a similar sample size. 37,40,57 Although we can be reasonably confident in identifying the time-invariant aspect of the SMPs, we cannot fully eliminate the possibility that the SMPs would remain stable in a more diverse IBD cohort, as the typical markers of inflammation (CRP and leukocyte count) were hardly different between both time points.

Conclusion
We observe considerable variability in DNA methylation measurements taken from PBL at 2 different time points separated by a median of 7 years. By contrast, around 14% of all CpG loci could be considered highly stable even after IBD-specific exposures during the 2 points. Focusing on these CpG loci during biomarker discovery might result in the identification of biologically relevant and more robust , ICC values of all CpG loci of class I and II HLA genes. The potential GV representing probes that were annotated with a genetic variant, predicted GV representing probes that presented a methylation signal typically found when driven by a GV, and methylation representing probes for which we have no evidence that they hybridize with any GV. Visualizations of the ICC values of all Illumina CpGs annotated to HLA-C (C), HLA-DPB1 (D), and HLA-DPA1 (E), relative to their position on each gene and grouped as potential GV (pink), predicted GV (green), or methylation (blue). Dots below represent known genetic variants as reported by Goyette53 for CD vs healthy controls (pink) and UC vs healthy controls (turquoise).
IBD-associated epigenetic biomarkers with an increased probability of independent replication.

Patient Selection
We performed a single-center, longitudinal EWAS, where we collected PBL samples from adult patients with IBD at the Amsterdam UMC. The interval between the time of sampling ranged from 2 to 9 years with a median of 7 (Figure 1, A). All included patients were historically diagnosed with either CD or UC on the basis of a combination of clinical symptoms and endoscopic inflammation as confirmed by histology per the current guidelines. 68,69 In addition, all patients received standard care followup. No additional inclusion or exclusion criteria were used as the goal was to collect a cohort of patients with IBD that reflected the overall IBD population at the Amsterdam UMC. This study was approved by the medical ethics committee of the academic medical hospital (METC NL24572.018.08 and NL53989.018.15), and written informed consent was obtained from all subjects prior to sampling.

Sample Collection and DNA Methylation Analysis
Whole peripheral blood samples were collected in a 6-mL EDTA tube and stored at À80 ºC until further processing. Genomic DNA was isolated using the QIAsymphony, whereupon the quantity of the DNA was assessed using the FLUOstar OMEGA and quality of the high-molecular weight DNA on a 0.8% agarose gel. Genomic DNA was bisulfite converted using the Zymo EZ DNA Methylation kit, randomized per plate to limit batch effects, and analyzed on the Illumina Human-Methylation EPIC BeadChip array at the Core Facility Genomics, Amsterdam UMC, Amsterdam, the Netherlands.

Statistical Analysis of Clinical Data
Baseline characteristics of all included patients were summarized using descriptive statistics. Categorical variables are presented as percentages and continuous variables as median annotated with the IQR. Differences in CRP and leukocyte count levels between T1 and T2 were calculated using the Wilcoxon signed ranks test. Analyses of clinical data were performed in IBM SPSS statistics version 26 and methylation analyses in the R statistical environment version 4.2.1.

Time-dependent DNA Methylation Data Analyses
For differential methylation analyses, raw DNA methylation data were imported into the R statistical environment using the Bioconductor minfi 70 package (version 1.36), whereupon the raw signal intensities were normalized using functional normalization 71 and converted into methylation ratios. Differential methylation analyses was performed using limma 72 (version 3.46) and eBayes 73 regressing against time point (T2 vs T1), gender, smoking behavior, disease, and blood cell distribution. Statistical significance was defined as an FDR-adjusted P-value < .05. In addition to identifying time-associated differences in methylation, we also investigated differences in methylation associated with CRP and leukocyte count. Blood cell estimations were performed using the IDOL predictor CpGs as reference. 74 Time-associated DMPs were investigated for their association with age by performing gene set enrichment analyses using the age-associated CpGs reported by Horvath, 42 Hannum, 43 Levine, 44 and Knight. 45 Visualizations were generated using ggplot2 75 (version 3.3.5) and gghighlight (version 0.3.2).

Time-independent DNA Methylation Data Analyses
For the time stability analyses, raw DNA methylation data were imported using ewastools to retain the 89 quality control probes that bind GVs. Methylation probes that might bind GVs were identified on the basis of the minfi-provided annotation files, which we termed as potential GVs. Additional GV-binding probes were estimated using the Gaphunter tool 55 as implemented in minfi, which we termed the predicted GVs. Moreover, as opposed to the differential methylation analysis, for the time stability analyses, we did not perform normalization nor did we correct for any other potential confounders (eg, gender, smoking behavior, disease, and blood cell distribution) to identify truly stable signals. Stability analyses were conducted using ICC analyses, where ICC estimates and their 95% confidence intervals were calculated using the irr package implemented in R. Specifically, a 2-way mixed, single measures, consistency analysis was performed. 76 Visualizations were generated using ggplot2 75 (version 3.3.5) and gghighlight (version 0.3.2).

Gene Ontology Enrichment and Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment Analyses
Functional enrichment analyses genes annotated to both stable and unstable methylated probes was performed using GOmeth 77 as implemented in missMethyl. 33 The gene ontology terms were grouped according to biological process, cellular component, and molecular function, and an FDR corrected P-value below .05 indicates a statistical significant difference. , and HLA-DRB1 (E) relative to their position on each gene and grouped as potential GV (pink), predicted GV (green), or methylation (blue). Dots below represent known genetic variants as reported by Goyette53 for CD vs healthy controls (pink) and UC vs healthy controls (turquoise).