Abstract:Background: Cytosine methylation is widespread in most eukaryotic genomes and is known to play a substantial role in various regulatory pathways. Unmethylated cytosines may be converted to uracil through the addition of sodium bisulphite, allowing genome-wide quantification of cytosine methylation via high-throughput sequencing. The data thus acquired allows the discovery of methylation 'loci'; contiguous regions of methylation consistently methylated across biological replicates. The mapping of these loci allows for associations with other genomic factors to be identified, and for analyses of differential methylation to take place. Results: The segmentSeq R package is extended to identify methylation loci from high-throughput sequencing data from multiple experimental conditions. A statistical model is then developed that accounts for biological replication and variable rates of non-conversion of cytosines in each sample to compute posterior likelihoods of methylation at each locus within an empirical Bayesian framework. The same model is used as a basis for analysis of differential methylation between multiple experimental conditions with the baySeq R package. We demonstrate the capability of this method to analyse complex data sets in an analysis of data derived from multiple Dicer-like mutants in Arabidopsis. This reveals several novel behaviours at distinct sets of loci in response to loss of one or more of the Dicer-like proteins that indicate an antagonistic relationship between the Dicer-like proteins at at least some methylation loci. Finally, we show in simulation studies that this approach can be significantly more powerful in the detection of differential methylation than many existing methods in data derived from both mammalian and plant systems. Conclusions: The methods developed here make it possible to analyse high-throughput sequencing of the methylome of any given organism under a diverse set of experimental conditions. The methods are able to identify methylation loci and evaluate the likelihood that a region is truly methylated under any given experimental condition, allowing for downstream analyses that characterise differences between methylated and non-methylated regions of the genome. Futhermore, diverse patterns of differential methylation may also be characterised from these data.

HMM-Fisher: identifying differential methylation using a hidden Markov model and Fisher’s exact test

HMM-DM: identifying differentially methylated regions using a hidden Markov model

Detection of differentially methylated regions from bisulfite-seq data by hidden Markov models incorporating genome-wide methylation level distributions

Spatially Enhanced Differential RNA Methylation Analysis from Affinity-Based Sequencing Data with Hidden Markov Model.

Differential analysis of RNA methylome with improved spatial resolution

Data-Driven-Based Approach to Identifying Differentially Methylated Regions Using Modified 1D Ising Model.

Identification of Differentially Methylated Sites with Weak Methylation Effects.

Identification of Methylation States of DNA Regions for Illumina Methylation BeadChip

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

A Hybrid Method Based on Beta-Binomial Bayesian Hierarchical Model and Combination of Ranking Method to Detect Differential Methylation Regions in Bisulfite Sequencing Data

Identifying Differential Methylation in Cancer Epigenetics via a Bayesian Functional Regression Model

Qdmr: A Quantitative Method For Identification Of Differentially Methylated Regions By Entropy

Evaluation of epitranscriptome-wide N6-methyladenosine differential analysis methods

DeepH&M: Estimating single-CpG hydroxymethylation and methylation levels from enrichment and restriction enzyme sequencing methods

DRME: Count-based Differential RNA Methylation Analysis at Small Sample Size Scenario

High accuracy methylation identification tools on single molecular level for PacBio HiFi data

A Novel Computational Method For Detecting Dna Methylation Sites With Dna Sequence Information And Physicochemical Properties

DNMHMM: an Approach to Identify the Differential Nucleosome Regions in Multiple Cell Types Based on a Hidden Markov Model.

Methods for discovering genomic loci exhibiting complex patterns of differential methylation

Detecting differentially methylated mRNA from MeRIP-Seq with likelihood ratio test.

FastDMA: an Infinium HumanMethylation450 Beadchip Analyzer