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.

Modelling ecological communities as if they were DNA

eDNAPlus: A unifying modelling framework for DNA-based biodiversity monitoring

A dynamic ancestral graph model and GPU-based simulation of a community based on metagenomic sampling

Modeling Community Dynamics Through Environmental Effects, Species Interactions and Movement

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Identifying Latent Dynamic Components in Biological Systems

Assembly processes inferred from eDNA surveys of a pond metacommunity are consistent with known species ecologies

Phylogeny structures species' interactions in experimental ecological communities

Deciphering microbial interactions and detecting keystone species with co-occurrence networks

Challenges in detecting ecological interactions using sedimentary ancient DNA data

Learning microbial interaction networks from metagenomic count data

A rapid phylogeny-based method for accurate community profiling of large-scale metabarcoding datasets

Inferring the Evolutionary Model of Community-Structuring Traits with Convolutional Kitchen Sinks

OxDNA to study species interactions

Modeling DNA methylation dynamics with approaches from phylogenetics

Modelling survival and allele complementation in the evolution of genomes with polymorphic loci

The Community Simulator: A Python package for microbial ecology

Novel community data in ecology -- properties and prospects

A Bayesian approach to inferring the phylogenetic structure of communities from metagenomic data

Modelling phylogeny in 16S rRNA gene sequencing datasets using string kernels

A Quantitative Test of Population Genetics Using Spatio-Genetic Patterns in Bacterial Colonies