Modeling DNA methylation profiles through a dynamic equilibrium between methylation and demethylation

DNA methylation is a heritable epigenetic mark that plays a key role in regulating gene expression. Mathematical modeling has been extensively applied to unravel the regulatory mechanisms of this process. In this study, we aimed to investigate DNA methylation by performing a high-depth analysis of particular loci, and by subsequent modeling of the experimental results. In particular, we performed an in-deep DNA methylation profiling of two genomic loci surrounding the transcription start site of the D-Aspartate Oxidase and the D-Serine Oxidase genes in different samples (n = 51). We found evidence of cell-to-cell differences in DNA methylation status. However, these cell differences were maintained between different individuals, which indeed showed very similar DNA methylation profiles. Therefore, we hypothesized that the observed pattern of DNA methylation was the result of a dynamic balance between DNA methylation and demethylation, and that this balance was identical between individuals. We hence developed a simple mathematical model to test this hypothesis. Our model reliably captured the characteristics of the experimental data, suggesting that DNA methylation and demethylation work together in determining the methylation state of a locus. Furthermore, our model suggested that the methylation status of neighboring cytosines plays an important role in this balance.