Genome-wide impact of hydrogen peroxide on maintenance DNA methylation in replicating cells
Epigenetics & Chromatin volume 14, Article number: 17 (2021)
Environmental factors, such as oxidative stress, have the potential to modify the epigenetic landscape of cells. We have previously shown that DNA methyltransferase (DNMT) activity can be inhibited by sublethal doses of hydrogen peroxide (H2O2). However, site-specific changes in DNA methylation and the reversibility of any changes have not been explored. Using bead chip array technology, differential methylation was assessed in Jurkat T-lymphoma cells following exposure to H2O2.
Sublethal H2O2 exposure was associated with an initial genome-wide decrease in DNA methylation in replicating cells, which was largely corrected 72 h later. However, some alterations were conserved through subsequent cycles of cell division. Significant changes to the variability of DNA methylation were also observed both globally and at the site-specific level.
This research indicates that increased exposure to H2O2 can result in long-term alterations to DNA methylation patterns, providing a mechanism for environmental factors to have prolonged impact on gene expression.
Epigenetic modification of chromatin provides a mechanistic basis through which environmental stimuli can modulate gene expression. Epigenetic marks, such as methylation of cytosine and histone modifications, are utilized during development to generate a wide range of cellular phenotypes from the same DNA code. While the majority of epigenetic changes in adult organisms are due to the continued differentiation of stem cells, they can also be modified during ageing [1, 2] and by interactions with environmental stimuli such as cigarette smoke, pollutants, and alcohol consumption [3,4,5,6,7,8,9]. However, there is a lack of mechanistic studies that investigate how environmental factors impact DNA methylation.
Several potential environmental modifiers are associated with increased oxidative stress, including aging and chronic inflammation [10, 11]. Mitochondrial dysfunction is linked to increased superoxide and hydrogen peroxide (H2O2) production in a range of pathologies . Ageing, smoking, radiation and xenobiotic metabolism can also result in increased H2O2 in biological systems [13,14,15,16,17,18,19]. Immune cells such as neutrophils and macrophages generate a range of reactive oxygen species following NADPH oxidase (NOX) activation, many of which are cell permeable, and contribute to oxidative stress during inflammation [11, 20, 21]. H2O2 and myeloperoxidase-derived oxidants produced by immune cells can oxidize redox-sensitive thiol proteins and methionine [22,23,24]. Importantly, methionine is a precursor of the methyl donor S-adenosyl-methionine (SAM).
We have previously reported that H2O2 and chloramines generated by phagocytic immune cells are able to inhibit DNA methyltransferase (DNMT) activity during DNA replication . These results suggested that this inhibition was due to the oxidation of DNMT1’s catalytic site cysteine and occurred at concentrations that did not have major impacts on cell proliferation. DNMTs copy DNA methylation to the nascent DNA strand during DNA replication, a process which ensures that the correct pattern of methylation is maintained throughout cellular division. However, in a cell culture system, only the chloramines had a significant impact on global DNA methylation because they also depleted SAM .
In this study, we have revisited the impact of H2O2 on maintenance DNA methylation in cultured cells by using a more sensitive genomic approach that individually monitors 850,000 CpG sites across the human genome. We have also investigated the reversible nature of any modifications analysing changes seen at early and late timepoints. We reveal that sublethal doses of H2O2 caused a significant disturbance to DNA methylation patterns, increasing overall variability and site-specific methylation in the 1–2 h immediately following exposure to this oxidant. Many of the changes were lost 72 h after exposure, but not all. Collectively, these findings describe a potential mechanism for various environmental factors and pathologies that promote oxidative stress to impact the methylation pattern of proliferating cells.
We propose that oxidative stress is able to inhibit the maintenance methylation that occurs following DNA replication in proliferating cells. Cultured Jurkat T-lymphoma cells were synchronized with a thymidine block to maximize the number undergoing DNA synthesis at the time of exposure to oxidative stress. Flow cytometry was used to monitor cell cycle transitions and determine the period of maximal DNA replication during the S-phase (Fig. 1a). The asynchronous population of Jurkat cells shows that just before cell division the majority of cells were in the G0/G1 phase, ~ 20% were in S-phase and ~ 14% of cells were in G2/M phase. The period of 2–5 h post-release was when the greatest number of cells entered the S-phase (Fig. 1b), and the time we chose to introduce H2O2.
H2O2 is consumed by Jurkat cells and media constituents within ~ 30 min of exposure [26, 27]. Cells were therefore exposed to two doses of H2O2, at 2 and 3 h post-release, to maximize the potential for inhibition of maintenance methylation during DNA replication. To determine the highest concentration of H2O2 that could be used without significant cell death, synchronized cells were exposed to a range of concentrations (Fig. 2a). Treatment with 15 µM H2O2 resulted in a large decrease in viability after 24 h, while 5 µM and 10 µM H2O2 treatments only had a small effect. With 10 µM H2O2, which was selected for subsequent experiments, there was no significant difference in viability between the treatment and control groups at 4 or 72 h (p = 0.25) (Fig. 2b). Assessment of cell growth indicated that only the control cells doubled in cell density after 24 h, whereas the H2O2-treated cells had significantly lower percentage growth (p = 0.02) (Fig. 2c). Cells were diluted to a concentration of 0.25 × 106 cells/ mL at 24 h, and the decrease in fold growth for the treatment samples compared to control was still observed at 48 and 72 h, although not statistically significant (p = 0.1 and 0.6). This slowing of growth rate was confirmed by the analysis of CFSE dilution by flow cytometry (Fig. 3). Collectively this data shows that both control and treated cells underwent at least one round of division by the 72-h time point.
DNA methylation analyses
DNA methylation profiles were assessed using the Illumina Infinium MethylationEPIC 850K array at 4 h and 72 h after the addition of cytidine to release cells from thymidine block, which was 2 and 70 h from the first H2O2 exposure, respectively. This technology assesses DNA methylation at over 850,000 CpG sites across the human genome, with nucleotide-level resolution.
Principal component analysis (PCA) of DNA methylation and hierarchical clustering
Multidimensional scaling was used to assess the top contributing factors in data variability, which corresponded with treatment, and to a lesser extent time (Fig. 4a). Overall, there was a distinct separation between samples from the control and H2O2 treatment groups. Samples from the control group formed a relatively tight cluster with less spatial separation, and as expected there was no separation between time points for the control samples. The H2O2 treatment samples from the 72-h time point more closely represented the control samples than the H2O2 treatment samples from control at the 4-h time point.
Hierarchical clustering was performed on a distance matrix calculated using the β-values for all probes, regardless of significance. Except for one sample, which incorrectly grouped within the treatment, this approach successfully differentiated between the control and H2O2 treatment groups (Fig. 4b).
Relative average DNA methylation change
At the 4-h time point, the mean methylation level of the H2O2 treatment group was 3.33% less than the control group (p = 0.013, t(5) = 3.72) (Fig. 5a). However, at 72 h this significant difference was no longer observed (p = 0.10, t(4) = 2.2). There was no significant difference between the 4-h and 72-h time points within each group, although the H2O2 treatment group demonstrated a minor increase and had substantially reduced variability (Fig. 5a).
Plotting the individual probe-wise comparisons of significance versus effect size change confirmed there was a substantial decrease in methylation after H2O2 treatment, with several probes from the 4-h time point demonstrating a significant log2FC greater than ~ 1. There was a comparative reduction in the effect size and significance at the 72-h time point, which indicates that the H2O2 treatment and control samples were more similar to each other than at the 4-h time point (Fig. 5b, c). Together these observations suggested that on a global scale, there were significant decreases in DNA methylation in the treatment group compared to the control group at 4 h, which were largely restored by 72 h (Fig. 5).
Analysis of variability in DNA methylation after H2O2 treatment
Principal component analysis and global methylation change suggested that there was a higher level of heterogeneity in the DNA methylation profiles after H2O2 treatment. To investigate the individual probes that contributed towards this observation, we assessed the differential variability at individual CpG positions (DVPs) using a general linear model that accounted for biological replicates. At the 4-h time point, 1035 DVPs demonstrated a significant increase in variability in the H2O2 treatment group compared to the control, and 286 demonstrated a significant decrease in variability. At the 72-h time point only 41 DVPs demonstrated an increase in variability and 115 demonstrated a decrease. Forty-six DVPs were common between the two time points, however, only 18% (7 DVPs) showed a consistent direction of change (Fig. 6).
Differentially methylated CpGs at the 4-h time point
To identify changes that associated with H2O2 treatment we assessed DNA methylation levels at each individual probe using the “topTreat” algorithm, within the R package, Limma. To summarize the results, each probe was mapped to the corresponding genomic region and unadjusted significance was plotted across the whole genome as a Manhattan plot (Fig. 7). This analysis demonstrated that the significant changes were distributed across the genome, and did not appear to be concentrated at specific genomic locations. There were 1162 probes that demonstrated a significant decrease in methylation using a log2FC (M-values) cutoff of 1, and 89 probes that demonstrated a significant increase in methylation (adj. p < 0.05). Increasing the log2FC cutoff to 1.2 identified 87 probes that demonstrated a significant decrease in methylation and three probes that demonstrated a significant increase (Additional file 1: Table S1). The top most significant probes had an initial methylation value of 20–25% in the control group, and a value of 10–15% in the H2O2 treatment group (Fig. 8a). Therefore, the raw data for the top six most significant probes that did not conform to this general trend are presented in Fig. 8 (including the three probes that demonstrated increased methylation). The top 20 significantly differentially methylated positions with a log2FC > 1.2 in order of largest fold change are presented in Table 1.
Differentially methylated CpGs at the 72-h time point
At the 72-h time point, 19 probes demonstrated a significant decrease in methylation and one probe demonstrated a significant increase in methylation (adj. p < 0.05), when assessed using a log2FC cutoff of 1 (Additional file 1: Table S2). The magnitude of these changes was less than 10%, therefore there were no significant probes detected at a log2FC cutoff of 1.2 at the 72-h time point.
Seventeen of the 20 significant probes detected at the 72-h time point were also detected in the 4-h time point when a log2FC of > 1 was used (Table 2) and the raw data from a subset of these probes are presented in Fig. 9. Two probes were also detected in the 4-h time point when a log2FC > 1.2 was used (cg05577994 mapping to gene: EMX2OS at chr10:119,254,968, adj. p < 0.05 and cg18363176 mapping to gene: FLI1 at chr11:128,606,110, adj. p < 0.04). Both probes had a smaller log2FC at 72 h than at 4 h. However, all other significant probes showed a substantially larger log2FC at 72 h than at 4 h (Table 2).
Top differentially methylated gene regions
We next investigated if multiple probes mapping to the same genomic loci demonstrated consistent methylation changes. We limited this analysis to differentially methylated regions (DMRs) that contained five or more adjacent significant CpGs, with an average change in methylation larger than 5%. Under these parameters, we detected 70 significant DMRs at the 4-h time point and five significant DMRs at the 72-h time point (Additional file 1: Tables S3 and S4). Two significant DMRs were consistent between the 4-h and 72-h time points (Table 3 and Fig. 10).
The genes corresponding to the top most significant CpG sites at 4 h, were assessed for significant biological relevance by pathway analysis within the GO and KEGG databases. Neither comparison identified significantly enriched pathways, after correction for multiple testing (Additional file 1: Table S5).
We used methylation arrays to investigate if exposure to H2O2 during DNA replication could interfere with the transfer of DNA methylation patterns in proliferating cells. H2O2 exposure during the S-phase of DNA replication disrupted DNA methylation of the nascent strand in proliferating cells 2 h after the first exposure to H2O2, with most but not all of these changes reversed by 72 h. By 72 h, all cells had undergone at least one round of division, indicating that these changes were representative of the epigenome of daughter cells. Sublethal H2O2 is known to alter gene expression in the short-term through impacting the activity of redox-sensitive transcription factors. Our results provide a mechanism by which H2O2 may have long-term impacts on gene expression and cell function.
Because aliquots of the treatment and control samples were split from the same source immediately prior to treatment, methylome changes relating to cell cycle were controlled for in the study design. The design of these experiments was such that DNA methylation was assessed 1–2 h after treatment, thus comparing the efficiency of nascent DNA methylation between the control and treatment samples. DNA methylation was then assessed again at 72 h to establish if the observed changes were corrected, or were inherited into the epigenome of subsequent cell generations. Cell viability and live cell counts indicated that approximately 72 h was required before cellular replication in the treatment group resembled that of the control group. This is likely to reflect the length of time required for repair of DNA damage, and re-establishment of the epigenome, and therefore DNA methylation was not assessed at the 24- and 48-h time points. An unsupervised assessment of data variability indicated that the treatment samples from the 4-h time point grouped distinctively from all other samples. There was no noticeable separation between the control samples, which more closely resembled the 72-h treatment samples.
Analysis of the average relative DNA methylation levels between groups indicated that H2O2 exposure corresponds with a significant decrease of ~ 3.3% in average methylation between treatment and control at the 4-h time point, however, there was no significant change between the 72-h time points. This change was driven by the combined effect of many probes with a log2FC approximately equal to one, however, at 72 h the effect sizes (and therefore significance) were substantially reduced. Hierarchical clustering was relatively successful at distinguishing the control samples from the treatment samples, and also identified that samples from the same biological replicate were more closely related, an effect which was subsequently controlled for in the statistical design.
These findings build on our previous work where we showed that while H2O2 inhibited DNMT, there was no detectable change in global methylation levels . The main experimental difference our previous work was that H2O2 exposure occurred at the time of cytidine release. We modified our previous protocol by treating the cells twice with 10 µM H2O2 during the S-phase of DNA replication, at 2 and 3 h after the replication block was released, in order to maximize H2O2 exposure during the time period where DNMTs are most active . Furthermore, while our previous study utilized a mass spectrometry approach which analysed the methylation of incorporated heavy-labelled deoxycytidine in the nascent strand, the current study used array technology that increased both the magnitude of the effect size and sensitivity, enabling detection of very small changes that were not previously possible using traditional methods. Our current work also highlighted the importance of nucleotide resolution, which enabled the detection of both increased and decreased methylation.
Intriguingly, there was a substantial increase in differential variability at both global and site-specific measurements. Differential variability of methylation at for the 4-h time point was significantly increased for the H2O2 treated samples compared to controls. Theoretically, his could either be due to increased variability at specific sites or random sites across the genome. We found that 1035 probes had significantly increased variability associated with H2O2 treatment at 4 h, which was reduced (41 probes) by 72 h. Therefore, site-specific increases in variability were contributing to the global variability in our system. It is important to note that there was also low consistency in the direction of change between the two time points, where 82% of the probes that were significant at both time points showed contrasting results.
The role of stochasticity in biological systems has received increasing attention in multiple contexts . Whether cell-intrinsic or extrinsic, stochasticity can impact on a wide range of cellular processes including epigenetics and cell differentiation. Multicellular organisms have a remarkable capacity to produce consistent phenotypic outcomes even when challenged by variable conditions [30, 31]. This “buffering capacity” of normal cells is dysregulated in cancer and can be unmasked by treatment with a toxic stress such as chemotherapy. In practical terms, epigenetic heterogeneity in response to an environmental stress enables a range of transcriptional states across a population, some of which might confer resistance. This has been observed in acute myeloid leukaemia where epigenetic heterogeneity has been correlated with poorer outcomes and a shorter time to relapse . Our finding that treatment with H2O2 increased the variability of methylation at the early time point is consistent with this and potentially points to a mechanism that allows a proportion of cancer cells across a population to adapt to oxidative stress.
We also investigated significant differentially methylated sites that demonstrated the largest fold change after H2O2 treatment. At the 4-h time point, 1162 probes showed a significant reduction in DNA methylation (log2FC > 1) in the treatment samples compared to control samples. In addition to this, 90 probes had an effect size larger than log2FC of 1.2, which suggests that in these instances the methylation level of the parental DNA strand may have been influenced . With the exception of two probes (mapping to EMX2OS and FLI1), these 90 probes did not show a significant change at the 72-h time point. This suggests that the large changes observed at the 4-h time point were eventually corrected, or occurred in cells that eventually died. The majority of the significant probes at the 72-h time had a smaller effect size at 4 h (log2FC ~ 1), which suggests that they represent a delayed effect. Although no significant enrichment of gene pathways were observed, several of the top most differentially methylated genes are disrupted in pathologies such as diabetes and cancer [34,35,36,37,38,39,40,41,42,43,44,45].
Two gene regions corresponding to the SPATA16 gene and LOC100506102 (also overlapping with the small nuclear RNA: snoU13.410–201, and the miRNA: RP11-322J23.1–001) demonstrated significant differential methylation at both 4 h and 72 h. The potential results of these changes are unclear; however, the SPATA16 gene is related with sperm count and motility. Extensive research has suggested that reactive oxygen species are a contributing factor to male infertility , with oxidative stress identified as a common pathology in approximately half of all infertile men . Deregulation of snoRNA has been observed in several different cancer types including chronic lymphocytic leukaemia , hepatocellular carcinoma , colorectal cancer  and endometrial cancer , and also metabolic stress disorder . Their roles as diagnostic and prognostic biomarkers have also been extensively studied [50, 53, 54].
Oxidants can directly damage and modify DNA, leading to mutations and genetic instability, but they can also cause interference in epigenetic pathways. One of the most common products of cellular oxidative stress is oxidized guanine or 8-hydroxy-2′-deoxyguanosine (8-OHdG) [55,56,57]. Increased frequencies of 8-OHdG in hepatocellular carcinoma patient DNA have been associated with increased oxidative stress, and in vitro studies have revealed that H2O2 exposure (50—250 µM) increased the occurrence of 8-OHdG and corresponded with altered methylation patterns in promoter regions of tumour suppressor genes . When 8-OHdG is formed in close proximity to CpG sites, it has been shown to interfere with DNMT1 activity by impeding its access to the nascent strand, resulting in impaired maintenance methylation . There is also evidence that 8-OHdG, through the recruitment of 8-oxoguanine DNA glycosylase-1, activates demethylation by ten-eleven translocation methylcytosine dioxygenase 1 (TET1) . However, one study has suggested that the cumulative occurrence of 8-OHdG is low , and no reports have demonstrated specific enrichment at CpG sites. Therefore, 8-OHdG formation due to H2O2 exposure, particularly at low concentrations, would be unlikely to account for the large genome-wide decreases in methylation observed in our study. Another potential mechanism for regulation of methylation levels by oxidants is decreased TET activity. TET enzymes are sensitive to redox regulation, as their iron-dependent catalytic regions require ascorbate for optimal activity [62,63,64,65]. Oxidants such as H2O2, chloramines and superoxide, ubiquitous in an environment of oxidative stress, have been shown to deplete ascorbate levels, and could contribute to reduced TET activity [66, 67].
Alternatively, since there have been several reports that H2O2 can directly inhibit DNMT1 activity, it is more likely that the transient inhibition of DNMT1 due to H2O2 exposure is largely responsible for the decreased methylomic changes observed [25, 68]. However, the methylation changes from the 4-h time point were spread across the genome, with no strong evidence for enrichment of specific gene pathways, and yet sequence-specific effects on DNA methylation were observed across four independent experiments. This result would not be anticipated if H2O2 were acting through inhibition of DNMT1 activity alone. One possibility is that since many of the changes were consistent across multiple adjacent probes, the effects of H2O2 treatment relate to the positioning of DNMT1 at the time of exposure. The cells were released from synchronization at the same time, therefore they would progress through the cell cycle at a similar rate. This may indicate that there are undiscovered mechanisms that combine with inhibition of DNMT activity to yield sequence-specific alterations in DNA methylation.
Most of the large epigenetic changes reported in the literature are associated with long-term exposure to environmental stressors such as smoking and chronic inflammatory disease states [5, 69,70,71,72], it was therefore unexpected to observe a large number of methylation changes with a relatively short oxidant exposure period. However, our findings have a substantial effect size, and the oxidant concentration used in this analysis is similar to in vivo concentrations observed during chronic infection [73,74,75,76]. Although it is unlikely that a large population of cells would exist in S-phase at any one time point in vivo, in the situation of prolonged inflammation, a large proportion of cells will be exposed to a range of oxidants during cellular replication. Therefore, it is conceivable that H2O2-induced inhibition of DNMT could also occur in an asynchronous population of replicating cells in an inflammatory environment. Furthermore, although the current study was limited by the utilization of a single cell type, DNMT inhibition due to H2O2 treatment has been reported in other cell lines , and therefore methylation changes as a result of oxidative stress should be futher explored in other cell lines as well as primary cells.
These data provide supporting evidence for a mechanistic link between environmental oxidative stress and DNA methylation. We propose that oxidative stress can influence DNMT1 activity, causing longer term changes in cytosine methylation patterns. Although it is accepted that environmental factors can influence the epigenome, few studies have investigated the pathways through which this can occur and the effect in subsequent cell generations. These data demonstrate that exposure to oxidative stress during cellular replication corresponds with substantial decreases in DNA methylation, and while the majority of these changes are subsequently corrected, a smaller number of sites remain demethylated. Furthermore, exposure to H2O2 resulted in increased variability of methylation in a site-specific manner. This is of great interest for cancer biology where metabolic reprogramming and the tumour microenvironment combine to create a milieu of oxidative stress that can promote altered epigenetic profiles and modify cell behaviour of both neoplastic and normal cells.
The results of this research indicate that the human Jurkat T-lymphoma epigenome is susceptible to sustained changes in DNA methylation after exposure to oxidative stress during cellular replication. These studies need to be expanded to other proliferating cells. H2O2 exposure was associated with genome-wide decreases in DNA methylation, which were largely corrected during subsequent cycles of cell division. The magnitude of the observed changes generally decrease with time, however, a subset of changes that displayed a relatively small initial effect size at 4 h demonstrated an increased effect size at 72 h. These results suggest that oxidant-induced changes in DNA methylation can be inherited by subsequent cell generations. Furthermore, we show that treatment with H2O2 affects variability of methylation in a global and site-specific manner.
Cell culture reagents
All experiments were performed using E6.1 Jurkat CD4 + /CD8- T-cell lymphoma cell line. Cells were cultured in RPMI 1640 + 2 mM glutamine + 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (all from Gibco, Life Technologies), and incubated in 5% CO2 at 37 °C. Hydrogen peroxide was diluted into PBS and its concentration determined spectrophotometrically (ε240 = 43.6 M−1 cm−1).
Cellular replication of E6.1 Jurkat cells was synchronized using a single thymidine block protocol to arrest cells between late G1 and early S-phase. Briefly, cells (1 million cells/ml) were treated with thymidine (2′-deoxythymidine, Sigma) at a final concentration of 1 mM and incubated for 18 h. To promote release into the S-phase, the block was alleviated by two media changes and addition of cytidine (2′-deoxycytidine, Sigma) to a final concentration of 50 µM. Cell cycle progression was assessed by flow cytometry every two hours over a 12-h period by analysing cellular DNA content following cell fixing in ethanol and staining with propidium iodide (PI) (Sigma). This approach was used to detect cells in each of the three major phases of the cycle (G0/G1, S and G2/M), and also allowed for detection of apoptotic cells with minimal DNA content.
Cell culture procedure
Each of four replicates originated from frozen cells of differing passage lineages, and the following procedure was performed for each replicate, on independent days. Cellular synchronization was performed on 200 × 106 cells in 200 mL as described above. Immediately after release, the cells were split into two equal volumes at 0.5 × 106/ mL for treatment and control samples. At two hours post-release, cells in media were treated with H2O2 at a final concentration of 10 µM. This was repeated one hour later, meaning the cells received two H2O2 bolus doses totaling 20 µM. Cells were split to a concentration of 250,000 cells/ mL after 24 h. Samples containing approximately five million cells, were collected at four hours and 72-h post-release for DNA extraction. Flow cytometry for cell viability (using PI), cell cycle analyses, and live cell counts were performed at pre-block, pre-release, and 4, 24, 48 and 72-h post-release time points.
Cell viability, proliferation and cell cycle progression
The percentage of viable cells at each major time point was assessed by the exclusion of PI. Cell cycle analyses were performed by fixation with ethanol followed by incubation with PI for 30 min. Cell proliferation was observed by labelling cells with carboxyfluorescein diacetate succinimidyl ester (CFSE; Invitrogen, New Zealand), as described previously .
DNA extraction and bisulfite conversion
DNA was extracted from cell samples using a GeneJet DNA extraction kit (ThermoFisher Scientific, USA) according to the manufacturer’s instructions. Sodium bisulfite treatment of genomic DNA was performed on 1000 ng of DNA using the Zymo Research EZ DNA methylation kit, according to the manufacturer’s recommended protocol for use on the Illumina Infinium MethylationEPIC 850K array. Cells harvested at 4 h and 72 h after cell cycle release were assessed for DNA methylation profiles using the Illumina methylationEPIC bead chip array. Methylated and non-methylated values were determined using the minfi pipeline, within R. To control for DNA methylation changes that occur naturally during cellular replication, the block-release protocol was performed on all cells together.
DNA methylation profiles of the control samples were subtracted from the treatment samples. Data analyses were performed using the Minfi and Limma Bioconductor software packages within the R statistical program (www.R-project.org). All packages, programs and pipelines used in these analyses are freely available, and the workflow was based upon published scripts . Any additional scripts are available from the authors on request. All samples passed quality control together using Subset-quantile Within Array Normalization (SWAN)  and variance-stabilizing transformation [80, 81]. Probes containing detection p-values > 0.05 for 1% or more samples, were excluded from further analysis. Because the cell samples were of the same sex, all chromosomes were retained, however, probes identified as having polymorphic hybridizing potential and homology to common SNPs  were removed. The final data frame contained 809,911 probes available for analysis. Multidimensional scaling of the top 1000 methylation values was performed using pairwise distance method for gene selection, and normalized, filtered methylation values. Hierarchical clustering was performed on a “minkowski” distance matrix calculated using the β-values for all probes, regardless of significance .
Average relative DNA methylation changes
Average relative DNA methylation changes were calculated using mean β-values from all probes for each sample in the four groups. The absolute mean differences between control and treatment groups at each time point were calculated using paired t-tests.
Differentially variable positions (DVPs) were identified using the DiffVar algorithm  within the missMethyl package . This analysis used a similar linear model as below, however, biological replicates were incorporated into the statistical design as covariates.
Identification of differentially methylated CpGs
The methylation status of each probe was calculated using normalized probe signals represented as methylation values (M-values) and β-values. M-values were generated within Minfi as the log2 ratio of the signal intensities of methylated probe divided by the unmethylated probe. Unless stated otherwise, statistical analyses were performed using the M-values. β-values (average DNA methylation level for each probe) were used for data visualization and range from 0 (unmethylated) to 1 (methylated). β-values were generated by dividing the methylated probe signal with the sum of the methylated and unmethylated probe signals . Normalized β-values and M-values were manually assessed for fit (Additional file 1: Figure S1).
Differentially methylated positions which correlated with treatment were identified using a linear regression model within the Limma package, with adjustment for multiple testing. This was calculated by comparing the methylation measurements of the control samples with the treatment samples, for each of the two time points (4-h post-release, and 72-h post-release). Samples that originated from the same cell passage were treated as a biological replicate and the correlation coefficient between these samples was incorporated into the statistical design using Limma’s duppcorr function. The top most significant, differentially methylated CpG positions were identified by log fold change (logFC) weighted functions using Limma’s topTreat algorithm. Adjustment for multiple correction was performed using the “Benjamini, Hochberg” method within Limma and Quantile–Quantile plots for the linear regression model were assessed at each time point using the observed against expected –log10(p-value) (Additional file 1: Fig. S2). The observed data show a minor inflation of p-values smaller than expected by chance. Therefore, we substantially reduced our criteria for significance by excluding probes that displayed less than a 10% mean difference in methylation between treatment and control samples. This also removed significant changes observed at sites that displayed methylation values less than 10% or greater than 90%, as changes within these ranges are less likely to be of biological relevance.
Pathway analysis was performed by comparison with the Kyoto Encyclopedia of Genes and Genomes (KEGG)  and gene ontology databases (GO)  using the missMethyl R package , with correction for probe bias.
Differentially methylated regions
Differentially methylated regions were interrogated within the Minfi package using the statistical package DMRcate . A methylation differential cutoff of 10 was used, and unless stated otherwise significance was determined using a false discovery rate (FDR) of 0.05 in conjunction with a p-value cutoff of 0.05.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files. Original scripts used for data analysis are available from the corresponding authors upon request.
Carboxyfluorescein diacetate succinimidyl ester
Principal component analysis
Differential variability at individual CpG positions
Differentially methylated regions
Gene ontology database
Kyoto Encyclopedia of Genes and Genomes
Bork S, et al. DNA methylation pattern changes upon long-term culture and aging of human mesenchymal stromal cells. Aging Cell. 2010;9:54–63.
Jiang S, Guo Y. Epigenetic Clock: DNA Methylation in Aging. Stem Cells Int. 2020;2020:1047896.
Anderson OS, Sant KE, Dolinoy DC. Nutrition and epigenetics: an interplay of dietary methyl donors, one-carbon metabolism and DNA methylation. J Nutr Biochem. 2012;23:853–9.
Rönn T, et al. A six months exercise intervention influences the genome-wide DNA methylation pattern in human adipose tissue. PLOS Genet. 2013;9:1–16.
Wan ES, et al. Cigarette smoking behaviors and time since quitting are associated with differential DNA methylation across the human genome. Hum Mol Genet. 2012;21:3073–82.
Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003;33:245–54.
Breitling LP, Yang R, Korn B, Burwinkel B, Brenner H. Tobacco-smoking-related differential DNA methylation: 27K discovery and replication. Am J Hum Genet. 2011;88:450–7.
Breton CV, et al. Prenatal tobacco smoke exposure affects global and gene-specific DNA methylation. Am J Respir Crit Care Med. 2009;180:462–7.
Waterland RA, Jirtle RL. Transposable elements: targets for early nutritional effects on epigenetic gene regulation. Mol Cell Biol. 2003;23:5293–300.
Beckman KB, Ames BN. The free radical theory of aging matures. Physiol Rev. 1998;78:547–81.
Federico A, Morgillo F, Tuccillo C, Ciardiello F, Loguercio C. Chronic inflammation and oxidative stress in human carcinogenesis. Int J Cancer. 2007;121:2381–6.
Chouchani ET, et al. Ischaemic accumulation of succinate controls reperfusion injury through mitochondrial ROS. Nature. 2014;515:431–5.
Harman D. The Biologic Clock: The Mitochondria? J Am Geriatr Soc. 1972;20:145–7.
Sohal RS, Sohal BH. Hydrogen peroxide release by mitochondria increases during aging. Mech Ageing Dev. 1991;57:187–202.
Kamceva G, et al. Cigarette Smoking and Oxidative Stress in Patients with Coronary Artery Disease. Open access Maced J Med Sci. 2016;4:636–40.
Church DF, Pryor WA. Free-radical chemistry of cigarette smoke and its toxicological implications. Environ Health Perspect. 1985;64:111–26.
Dayal D, Martin SM, Limoli CL, Spitz DR. Hydrogen peroxide mediates the radiation-induced mutator phenotype in mammalian cells. Biochem J. 2008;413:185–91.
Bhattacharyya S, Sinha K, Sil PC. Cytochrome P450s: mechanisms and biological implications in drug metabolism and its interaction with oxidative stress. Curr Drug Metab. 2014;15:719–42.
Gómez-Tabales J, García-Martín E, Agúndez JAG, Gutierrez-Merino C. Modulation of CYP2C9 activity and hydrogen peroxide production by cytochrome b 5. Sci Rep. 2020;10:1–13.
Winterbourn CC, Kettle AJ, Hampton MB. Reactive Oxygen Species and Neutrophil Function. Annu Rev Biochem. 2016;85:765–92.
Winterbourn CC. Reconciling the chemistry and biology of reactive oxygen species. Nat Chem Biol. 2008;4:278–86.
Grisham MB, Jefferson MM, Melton DF, Thomas EL. Chlorination of endogenous amines by isolated neutrophils Ammonia-dependent bactericidal, cytotoxic, and cytolytic activities of the chloramines. J Biol Chem. 1984;259:10404–13.
Midwinter RG, Cheah FC, Moskovitz J, Vissers MC, Winterbourn CC. IκB is a sensitive target for oxidation by cell-permeable chloramines: Inhibition of NF-κB activity by glycine chloramine through methionine oxidation. Biochem J. 2006;396:71–8.
Laporte A, Lortz S, Schaal C, Lenzen S, Elsner M. Hydrogen peroxide permeability of cellular membranes in insulin-producing cells. Biochim Biophys Acta Biomembr. 2020;1862:183096.
OConnor, K. M., Das, A. B., Winterbourn, C. C. & Hampton, M. B. Inhibition of DNA methylation in proliferating human lymphoma cells by immune cell oxidants. J. Biol. Chem. (2020).
Hampton MB, Orrenius S. Dual regulation of caspase activity by hydrogen peroxide: implications for apoptosis. FEBS Lett. 1997;414:552–6.
Hua Long, L. & Halliwell, B. Oxidation and Generation of Hydrogen Peroxide by Thiol Compounds in Commonly Used Cell Culture Media. Biochem. Biophys. Res. Commun. 286, 991–994 (2001).
Desjobert C, et al. Combined analysis of DNA methylation and cell cycle in cancer cells. Epigenetics. 2015;10:82–91.
Symmons O, Raj A. What’s luck got to do with it: single cells, multiple fates, and biological nondeterminism. Mol Cell. 2016;62:788–802.
Waddington CH. The epigenotype. Int J Epidemiol. 2012;41:10–3.
Pujadas E, Feinberg AP. Regulated noise in the epigenetic landscape of development and disease. Cell. 2012;148:1123–31.
Li S, et al. Distinct evolution and dynamics of epigenetic and genetic heterogeneity in acute myeloid leukemia. Nat Med. 2016;22:792–9.
Valinluck V, Sowers LC. Endogenous cytosine damage products alter the site selectivity of human DNA maintenance methyltransferase DNMT1. Cancer Res. 2007;67:946–50.
Shan Z, et al. Interactions between zinc transporter-8 gene (SLC30A8) and plasma zinc concentrations for impaired glucose regulation and type 2 diabetes. Diabetes. 2014;63:1796–803.
Lee YH, et al. Association between polymorphisms in SLC30A8, HHEX, CDKN2A/B, IGF2BP2, FTO, WFS1, CDKAL1, KCNQ1 and type 2 diabetes in the Korean population. J Hum Genet. 2008;53:991–8.
Xing M. Oxidative stress: a new risk factor for thyroid cancer. Endocr Relat Cancer. 2012;19:C7–11.
Addis L, et al. Microdeletions of ELP4 are associated with language impairment, autism spectrum disorder, and mental retardation. Hum Mutat. 2015;36:842–50.
Rutter GA, Chimienti F. SLC30A8 mutations in type 2 diabetes. Diabetologia. 2015;58:31–6.
Lau W, Andrew T, Maniatis N. High-resolution genetic maps identify multiple type 2 diabetes loci at regulatory hotspots in African Americans and Europeans. Am J Hum Genet. 2017;100:803–16.
Chi GC, et al. Long-term outdoor air pollution and DNA methylation in circulating monocytes: Results from the Multi-Ethnic Study of Atherosclerosis (MESA). Environ Heal A Glob Access Sci Source. 2016;15:1–12.
van der Vaart, H., Postma, D. S., Timens, W. & Ten Hacken, N. H. T. Acute effects of cigarette smoke on inflammation and oxidative stress: a review. Thorax 59, 713 LP – 721 (2004).
Seman NA, Mohamud WNW, Östenson C-G, Brismar K, Gu HF. Increased DNA methylation of the SLC30A8 gene promoter is associated with type 2 diabetes in a Malay population. Clin Epigenetics. 2015;7:30.
Voisin S, et al. Many obesity-associated SNPs strongly associate with DNA methylation changes at proximal promoters and enhancers. Genome Med. 2015;7:103.
Nikoshkov A, et al. Epigenetic DNA methylation in the promoters of the Igf1 receptor and insulin receptor genes in db/db mice. Epigenetics. 2011;6:405–9.
Meyer KF, et al. The fetal programming effect of prenatal smoking on Igf1r and Igf1 methylation is organ- and sex-specific. Epigenetics. 2017;12:1076–91.
Dam AHDM, et al. Homozygous mutation in SPATA16 is associated with male infertility in human globozoospermia. Am J Hum Genet. 2007;81:813–20.
Tremellen K. Oxidative stress and male infertility—a clinical perspective. Hum Reprod Update. 2008;14:243–58.
Ronchetti D, et al. Small nucleolar RNAs as new biomarkers in chronic lymphocytic leukemia. BMC Med Genomics. 2013;6:27.
Xu G, et al. Small nucleolar RNA 113–1 suppresses tumorigenesis in hepatocellular carcinoma. Mol Cancer. 2014;13:1–14.
Okugawa Y, et al. Clinical significance of SNORA42 as an oncogene and a prognostic biomarker in colorectal cancer. Gut. 2017;66:107–17.
Ravo M, et al. Small non-coding RNA deregulation in endometrial carcinogenesis. Oncotarget. 2015;6:4677–91.
Michel CI, et al. Small nucleolar RNAs U32a, U33, and U35a are critical mediators of metabolic stress. Cell Metab. 2011;14:33–44.
Liao J, et al. Small nucleolar RNA signatures as biomarkers for non-small-cell lung cancer. Mol Cancer. 2010;9:198.
Krishnan P, et al. Profiling of small nucleolar RNAs by next generation sequencing: Potential new players for breast cancer prognosis. PLoS ONE. 2016;11:1–17.
Cheng KC, Cahill DS, Kasai H, Nishimura S, Loeb LA. 8-Hydroxyguanine, an abundant form of oxidative DNA damage, causes G → T and A → C substitutions. J Biol Chem. 1992;267:166–72.
Shigenaga MK, Gimeno CJ, Ames BN. Urinary 8-hydroxy-2’-deoxyguanosine as a biological marker of in vivo oxidative DNA damage. Proc Natl Acad Sci U S A. 1989;86:9697–701.
Kuchino Y, et al. Misreading of DNA templates containing 8-hydroxydeoxyguanosine at the modified base and at adjacent residues. Nature. 1987;327:77–9.
Nishida N, et al. Reactive oxygen species induce epigenetic instability through the formation of 8-hydroxydeoxyguanosine in human hepatocarcinogenesis. Dig Dis. 2013;31:459–66.
Turk PW, Laayoun A, Smith SS, Weitzman SA. DNA adduct 8-hydroxyl-2’-deoxyguanosine (8-hydroxyguanine) affects function of human DNA methyltransferase. Carcinogenesis. 1995;16:1253–5.
Zhou X, et al. OGG1 is essential in oxidative stress induced DNA demethylation. Cell Signal. 2016;28:1163–71.
Fraga CG, Shigenaga MK, Park JW, Degan P, Ames BN. Oxidative damage to DNA during aging: 8-Hydroxy-2’-deoxyguanosine in rat organ DNA and urine. Proc Natl Acad Sci U S A. 1990;87:4533–7.
Yin R, et al. Ascorbic acid enhances tet-mediated 5-methylcytosine oxidation and promotes DNA demethylation in mammals. J Am Chem Soc. 2013;135:10396–403.
Blaschke K, et al. Vitamin C induces Tet-dependent DNA demethylation and a blastocyst-like state in ES cells. Nature. 2013;500:222–6.
Vissers, M. C. M. & Das, A. B. Ascorbate as an enzyme cofactor. in VNew Biochemical and Functional Insights (eds. Chen", ["Qi & Vissers"], "Margreet C M) 28 (CRC Press, 2020).
Das, A. B., Smith-Díaz, C. C. & Vissers, M. C. M. Emerging epigenetic therapeutics for myeloid leukemia: modulating demethylase activity with ascorbate. Haematologica Online ahe, (2020).
Parker A, Cuddihy SL, Son TG, Vissers MCM, Winterbourn CC. Roles of superoxide and myeloperoxidase in ascorbate oxidation in stimulated neutrophils and H 2O 2-treated HL60 cells. Free Radic Biol Med. 2011;51:1399–405.
Thomas EL, Learn DB, Jefferson MM, Weatherred W. Superoxide-dependent oxidation of extracellular reducing agents by isolated neutrophils. J Biol Chem. 1988;263:2178–86.
Kilgore JA, et al. Identification of DNMT1 selective antagonists using a novel scintillation proximity assay. J Biol Chem. 2013;288:19673–84.
Shenker NS, et al. Epigenome-wide association study in the European Prospective Investigation into Cancer and Nutrition (EPIC-Turin) identifies novel genetic loci associated with smoking. Hum Mol Genet. 2012;22:843–51.
Guida F, et al. Dynamics of smoking-induced genome-wide methylation changes with time since smoking cessation. Hum Mol Genet. 2015;24:2349–59.
Glória L, et al. DNA hypomethylation and proliferative activity are increased in the rectal mucosa of patients with long-standing ulcerative colitis. Cancer. 1996;78:2300–6.
Chandra A, Senapati S, Roy S, Chatterjee G, Chatterjee R. Epigenome-wide DNA methylation regulates cardinal pathological features of psoriasis. Clin Epigenetics. 2018;10:108.
Roberts CK, et al. A high-fat, refined-carbohydrate diet induces endothelial dysfunction and oxidant/antioxidant imbalance and depresses NOS protein expression. J Appl Physiol. 2005;98:203–10.
Tsukimori K, Yoshitomi T, Morokuma S, Fukushima K, Wake N. Serum uric acid levels correlate with plasma hydrogen peroxide and protein carbonyl levels in preeclampsia. Am J Hypertens. 2008;21:1343–6.
Forman HJ, Bernardo A, Davies KJA. What is the concentration of hydrogen peroxide in blood and plasma? Arch Biochem Biophys. 2016;603:48–53.
Kalyanaraman B, Sohnle PG. Generation of free radical intermediates from foreign compounds by neutrophil-derived oxidants. J Clin Invest. 1985;75:1618–22.
Hock BD, et al. Effect of activated human polymorphonuclear leucocytes on T lymphocyte proliferation and viability. Immunology. 2012;137:249–58.
Maksimovic J, Phipson B, Oshlack A. A cross-package Bioconductor workflow for analysing methylation array data. F1000Research. 2017;5:1–51.
Maksimovic J, Gordon L, Oshlack A. SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol. 2012;13:1–12.
Shi, W., Oshlack, A. & Smyth, G. K. Optimizing the noise versus bias trade-off for Illumina whole genome expression BeadChips. Nucleic Acids Res. 38, (2010).
Lin SM, Du P, Huber W, Kibbe WA. Model-based variance-stabilizing transformation for Illumina microarray data. Nucleic Acids Res. 2008;36:1–9.
Pidsley R, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17:1–17.
Stevens AJ, et al. Methylomic changes in response to micronutrient supplementation and MTHFR genotype. Epigenomics. 2018;10:1201–14.
Phipson B, Oshlack A. DiffVar: a new method for detecting differential variability with application to methylation in cancer and aging. Genome Biol. 2014;15:465.
Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32:286–8.
Du P, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28:27–30.
Carlson, M., Falcon, S., Pages, H. & Li, N. GO.db: A set of annotation maps describing the entire Gene Ontology. R package version 2.9.0 (2007).
Peters TJ, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015;8:6.
We are grateful to Diane Kenwright (University of Otago, Wellington, New Zealand), for providing supporting resources. Dr. John Pearson for advice with the bioinformatic analysis.
This work was supported by the Health Research Council of New Zealand, Lotteries Health Research Commission of New Zealand, the Canterbury Medical Research Foundation of New Zealand and the Postgraduate Tassell Scholarship for Cancer Research.
Ethics approval and consent to participate
Ethics approval for this study was obtained from the University of Otago Ethics committee, reference URA/06/12/083/AM06.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Seddon, A.R., Liau, Y., Pace, P.E. et al. Genome-wide impact of hydrogen peroxide on maintenance DNA methylation in replicating cells. Epigenetics & Chromatin 14, 17 (2021). https://doi.org/10.1186/s13072-021-00388-6
- DNA methylation
- DNA methyltransferase
- Oxidative stress
- Hydrogen peroxide
- Redox signalling