Abstract:Deoxyribonucleic acid (DNA) methylation is an epigenetic alteration crucial for regulating stress responses. Identifying large-scale DNA methylation at single nucleotide resolution is made possible by whole genome bisulfite sequencing. An essential task following the generation of bisulfite sequencing data is to detect differentially methylated cytosines (DMCs) among treatments. Most statistical methods for DMC detection do not consider the dependency of methylation patterns across the genome, thus possibly inflating type I error. Furthermore, small sample sizes and weak methylation effects among different phenotype categories make it difficult for these statistical methods to accurately detect DMCs. To address these issues, the wavelet-based functional mixed model (WFMM) was introduced to detect DMCs. To further examine the performance of WFMM in detecting weak differential methylation events, we used both simulated and empirical data and compare WFMM performance to a popular DMC detection tool methylKit. Analyses of simulated data that replicated the effects of the herbicide glyphosate on DNA methylation in Arabidopsis thaliana show that WFMM results in higher sensitivity and specificity in detecting DMCs compared to methylKit, especially when the methylation differences among phenotype groups are small. Moreover, the performance of WFMM is robust with respect to small sample sizes, making it particularly attractive considering the current high costs of bisulfite sequencing. Analysis of empirical Arabidopsis thaliana data under varying glyphosate dosages, and the analysis of monozygotic (MZ) twins who have different pain sensitivities-both datasets have weak methylation effects of <1%-show that WFMM can identify more relevant DMCs related to the phenotype of interest than methylKit. Differentially methylated regions (DMRs) are genomic regions with different DNA methylation status across biological samples. DMRs and DMCs are essentially the same concepts, with the only difference being how methylation information across the genome is summarized. If methylation levels are determined by grouping neighboring cytosine sites, then they are DMRs; if methylation levels are calculated based on single cytosines, they are DMCs.

DMseg: a Python algorithm for de novo detection of differentially or variably methylated regions

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

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate

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

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

Identification of Differentially Methylated Sites with Weak Methylation Effects.

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

Identification of Functionally Methylated Regions Based on Discriminant Analysis Through Integrating Methylation and Gene Expression Data

A comprehensive evaluation of computational tools to identify differential methylation regions using RRBS data

A Comprehensive Evaluation of Alignment Software for Reduced Representation Bisulfite Sequencing Data

RankDMG: a differentially methylated gene analysis method for three special scenarios

Evaluation of epitranscriptome-wide N6-methyladenosine differential analysis methods

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

Detecting Differentially Methylated Loci for Illumina Array Methylation Data Based on Human Ovarian Cancer Data.

Detection of differentially methylated CpG sites between tumor samples with uneven tumor purities.

Metrics for Evaluating Differentially Methylated Region Sets Predicted from BS-seq Data.

Abstract P390: E-value over P-value and Its Adjustments: Improving Dna Methylation Analysis with an Application to Hypertension Studies

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

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

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

Detecting differentially methylated regions in bisulfite sequencing data using quasi-binomial mixed models with smooth covariate effect estimates