Abstract:MOTIVATION: Methylation of CpG dinucleotides is a prevalent epigenetic modification that is required for proper development in vertebrates. Genome-wide DNA methylation assays have become increasingly common, and this has enabled characterization of DNA methylation in distinct stages across differentiating cellular lineages. Changes in CpG methylation are essential to cellular differentiation; however, current methods for modeling methylation dynamics do not account for the dependency structure between precursor and dependent cell types.RESULTS: We developed a continuous-time Markov chain approach, based on the observation that changes in methylation state over tissue differentiation can be modeled similarly to DNA nucleotide changes over evolutionary time. This model explicitly takes precursor to descendant relationships into account and enables inference of CpG methylation dynamics. To illustrate our method, we analyzed a high-resolution methylation map of the differentiation of mouse stem cells into several blood cell types. Our model can successfully infer unobserved CpG methylation states from observations at the same sites in related cell types (90% correct), and this approach more accurately reconstructs missing data than imputation based on neighboring CpGs (84% correct). Additionally, the single CpG resolution of our methylation dynamics estimates enabled us to show that DNA sequence context of CpG sites is informative about methylation dynamics across tissue differentiation. Finally, we identified genomic regions with clusters of highly dynamic CpGs and present a likely functional example. Our work establishes a framework for inference and modeling that is well suited to DNA methylation data, and our success suggests that other methods for analyzing DNA nucleotide substitutions will also translate to the modeling of epigenetic phenomena.AVAILABILITY AND IMPLEMENTATION: Source code is available at www.kostkalab.net/software.

Detecting DNA Modifications from SMRT Sequencing Data by Modeling Sequence Context Dependence of Polymerase Kinetic.

Modeling Kinetic Rate Variation in Third Generation DNA Sequencing Data to Detect Putative Modifications to DNA Bases

data to detect putative modifications to DNA bases Modeling kinetic rate variation in third generation DNA sequencing

Determining RNA Natural Modifications and Nucleoside Analog-Labeled Sites by a Chemical/Enzyme-Induced Base Mutation Principle

Qdnamod: a Statistical Model-Based Tool to Reveal Intercellular Heterogeneity of DNA Modification from SMRT Sequencing Data

Context-dependent DNA polymerization effects can masquerade as DNA modification signals

SMAC: identifying DNA N -methyladenine (6mA) at the single-molecule level using SMRT CCS data

SMAC: Identifying DNA N6-methyladenine (6ma) at the Single-Molecule Level Using SMRT CCS Data

[The principle and application of the single-molecule real-time sequencing technology]

Discrimination Between Different DNA Lesions by Monitoring Single-Molecule Polymerase Stalling Kinetics During Nanopore Sequencing

Sequence Terminus Dependent PCR for Site-Specific Mutation and Modification Detection

Nanopore strand-specific mismatch enables de novo detection of bacterial DNA modifications

Modeling DNA methylation dynamics with approaches from phylogenetics

Adapting nanopore sequencing basecalling models for modification detection via incremental learning and anomaly detection

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Probabilistic Model Based Error Correction of Various Mutant Sequences Analyzed by the Single-Molecule Real-Time Sequencing

Correcting modification-mediated errors in nanopore sequencing by nucleotide demodification and reference-based correction

Detection of DNA base modifications by deep recurrent neural network on Oxford Nanopore sequencing data

Gaussian mixture model-based unsupervised nucleotide modification number detection using nanopore-sequencing readouts

De novo Identification of DNA Modifications Enabled by Genome-Guided Nanopore Signal Processing

A Hybrid HMM Approach for the Dynamics of DNA Methylation