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.

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

Addressing dispersion in mis‐measured multivariate binomial outcomes: A novel statistical approach for detecting differentially methylated regions in bisulfite sequencing data

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

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

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

A Bayesian Approach for Analysis of Whole-Genome Bisulfite Sequencing Data Identifies Disease-Associated Changes in DNA Methylation

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation

A Statistical Method for Detecting Differentially Methylated Gene Functional Regions Based on Case-Control Studies

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

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

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

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

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

Identification of Differentially Methylated Sites with Weak Methylation Effects.

MethylCal: Bayesian calibration of methylation levels

Estimation of the methylation pattern distribution from deep sequencing data

A novel method for detecting association between DNA methylation and diseases using spatial information

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

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

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate

A New Approach to Decode DNA Methylome and Genomic Variants Simultaneously from Double Strand Bisulfite Sequencing.