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.

A family of mixture models for beta valued DNA methylation data

A novel family of beta mixture models for the differential analysis of DNA methylation data: An application to prostate cancer

A Nonparametric Bayesian Approach for Clustering Bisulfate-Based DNA Methylation Profiles

HMM-DM: identifying differentially methylated regions using a 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

BDBB: A Novel Beta-distribution-based Biclustering Algorithm for Revealing Local Co-methylation Patterns in Epi-transcriptome Profiling Data

Clustering DNA methylation expressions using nonparametric beta mixture model

BBM: A novel beta-binomial-distribution-based biclustering algorithm for mining m 6 A co-methylation patterns

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

A Variational Bayes Beta Mixture Model for Feature Selection in Dna Methylation Studies

MBMM: Moment Estimating Beta Mixture Model-based Clustering Algorithm for M6a Co-methylation Module Mining

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

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

Identification of Differentially Methylated Sites with Weak Methylation Effects.

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

Empirical Bayes Model Comparisons for Differential Methylation Analysis

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

A Bayesian hierarchical model to estimate DNA methylation conservation in colorectal tumors

A Bayesian framework for genome-wide inference of DNA methylation levels

Classification of Large DNA Methylation Datasets for Identifying Cancer Drivers

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