Abstract:MOTIVATIONWith the advent of high-throughput sequencing technology, bisulfite-sequencing-based DNA methylation profiling methods have emerged as the most promising approaches due to their single-base resolution and genome-wide coverage. However, statistical analysis methods for analyzing this type of methylation data are not well developed. Although the most widely used proportion-based estimation method is simple and intuitive, it is not statistically adequate in dealing with the various sources of noise in bisulfite-sequencing data. Furthermore, it is not biologically satisfactory in applications that require binary methylation status calls.RESULTSIn this article, we use a mixture of binomial model to characterize bisulfite-sequencing data, and based on the model, we propose to use a classification-based procedure, called the methylation status calling (MSC) procedure, to make binary methylation status calls. The MSC procedure is optimal in terms of maximizing the overall correct allocation rate, and the false discovery rate (FDR) and false non-discovery rate (FNDR) of MSC can be estimated. To control FDR at any given level, we further develop an FDR-controlled MSC procedure, which combines a local FDR-based adaptive procedure with the MSC procedure. Both simulation study and real data application are carried out to examine the performance of the proposed procedures. It is shown in our simulation study that the estimates of FDR and FNDR of the MSC procedure are appropriate. Simulation study also demonstrates that the FDR-controlled MSC procedure is valid in controlling FDR at a prespecified level and is more powerful than the individual binomial testing procedure. In the real data application, the MSC procedure exhibits an estimated FDR of 0.1426 and an estimated FNDR of 0.0067. The overall correct allocation rate is >0.97. These results suggest the effectiveness of our proposed procedures.AVAILABILITY AND IMPLEMENTATIONThe proposed procedures are implemented in R and are available at http://www.stat.purdue.edu/*cheng70/code.html.

HBI: a hierarchical Bayesian interaction model to estimate cell-type-specific methylation quantitative trait loci incorporating priors from cell-sorted bisulfite sequencing data

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

Comprehensive mechanistic characterization of mQTLs in an Asian population

IMAGE: High-Powered Detection of Genetic Effects on DNA Methylation Using Integrated Methylation QTL Mapping and Allele-Specific Analysis

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

A unified framework for cell-type-specific eQTLs prioritization by integrating bulk and scRNA-seq data

Mapping QTLs with Epistatic Effects and QTL×environment Interactions by Mixed Linear Model Approaches

CHIT: an allele-specific method for testing the association between molecular quantitative traits and phenotype–genotype interaction

BayMeth: improved DNA methylation quantification for affinity capture sequencing data using a flexible Bayesian approach

Bayesian estimation of cell type–specific gene expression with prior derived from single-cell data

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Efficient and accurate mixed model association tool for single-cell eQTL analysis

Bayesian Modeling of Haplotype Effects in Multiparent Populations.

Methylation QTLs Are Associated with Coordinated Changes in Transcription Factor Binding, Histone Modifications, and Gene Expression Levels

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

A Classification Approach for DNA Methylation Profiling with Bisulfite Next-Generation Sequencing Data

Genetic control of DNA methylation is largely shared across European and East Asian populations

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

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

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