{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2024,8,20]],"date-time":"2024-08-20T15:13:07Z","timestamp":1724166787737},"reference-count":26,"publisher":"Oxford University Press (OUP)","issue":"2","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2014,1,15]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:p>Motivation: With the advent of high-throughput sequencing technology, bisulfite\u2013sequencing-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.<\/jats:p><jats:p>Results: In 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 &amp;gt;0.97. These results suggest the effectiveness of our proposed procedures.<\/jats:p><jats:p>Availability and implementation: The proposed procedures are implemented in R and are available at http:\/\/www.stat.purdue.edu\/\u223ccheng70\/code.html.<\/jats:p><jats:p>Contact: \u00a0cheng70@purdue.edu<\/jats:p><jats:p>Supplementary information: \u00a0Supplementary data are is available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btt674","type":"journal-article","created":{"date-parts":[[2013,11,23]],"date-time":"2013-11-23T03:29:55Z","timestamp":1385177395000},"page":"172-179","source":"Crossref","is-referenced-by-count":14,"title":["A classification approach for DNA methylation profiling with bisulfite next-generation sequencing data"],"prefix":"10.1093","volume":"30","author":[{"given":"Longjie","family":"Cheng","sequence":"first","affiliation":[{"name":"Department of Statistics, Purdue University, West Lafayette, IN 47906, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yu","family":"Zhu","sequence":"additional","affiliation":[{"name":"Department of Statistics, Purdue University, West Lafayette, IN 47906, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2013,11,21]]},"reference":[{"key":"2023012710392670700_btt674-B1","doi-asserted-by":"crossref","first-page":"286","DOI":"10.1080\/01621459.1985.10478110","article-title":"Estimation of allocation rates in a cluster analysis context","volume":"80","author":"Basford","year":"1985","journal-title":"J. Am. Stat. Assoc."},{"key":"2023012710392670700_btt674-B2","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1111\/j.2517-6161.1995.tb02031.x","article-title":"Controlling the false discovery rate: a practical and powerful approach to multiple testing","volume":"57","author":"Benjamini","year":"1995","journal-title":"J. R. Stat. Soc. Series B Stat. Methodol."},{"key":"2023012710392670700_btt674-B3","doi-asserted-by":"crossref","first-page":"1106","DOI":"10.1038\/nbt.1681","article-title":"Quantitative comparison of genome-wide DNA methylation mapping technologies","volume":"28","author":"Bock","year":"2010","journal-title":"Nat. Biotechnol."},{"key":"2023012710392670700_btt674-B4","doi-asserted-by":"crossref","first-page":"1715","DOI":"10.1093\/bioinformatics\/btp312","article-title":"Hierarchical Hidden Markov Model with application to joint analysis of ChIP-chip and ChIP-seq data","volume":"25","author":"Choi","year":"2009","journal-title":"Bioinformatics"},{"key":"2023012710392670700_btt674-B5","doi-asserted-by":"crossref","first-page":"1351","DOI":"10.1214\/009053606000001460","article-title":"Size, power and false discovery rates","volume":"35","author":"Efron","year":"2007","journal-title":"Ann. Stat."},{"key":"2023012710392670700_btt674-B6","first-page":"143","article-title":"A modified false discovery rate multiple-comparisons procedure for discrete data, applied to human immunodeficiency virus genetics","volume":"54","author":"Gilbert","year":"2005","journal-title":"Appl. Stat."},{"key":"2023012710392670700_btt674-B7","doi-asserted-by":"crossref","first-page":"133","DOI":"10.1038\/nmeth.1414","article-title":"Genome-scale DNA methylation mapping of clinical samples at single-nucleotide resolution","volume":"7","author":"Gu","year":"2010","journal-title":"Nat. Methods"},{"key":"2023012710392670700_btt674-B8","doi-asserted-by":"crossref","first-page":"1097","DOI":"10.1038\/nbt.1682","article-title":"Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications","volume":"28","author":"Harris","year":"2010","journal-title":"Nat. Biotechnol."},{"key":"2023012710392670700_btt674-B9","doi-asserted-by":"crossref","first-page":"43","DOI":"10.1142\/9789814329804_0003","article-title":"A false discovery rate procedure for categorical data","volume-title":"Recent Advances in Biostatistics: False Discovery Rates, Survival Analysis, and Related Topics","author":"Heyse","year":"2011"},{"key":"2023012710392670700_btt674-B10","doi-asserted-by":"crossref","first-page":"145","DOI":"10.1038\/nmeth.1828","article-title":"DNA methylome analysis using short bisulfite sequencing data","volume":"9","author":"Krueger","year":"2012","journal-title":"Nat. Methods"},{"key":"2023012710392670700_btt674-B11","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1038\/nrg2732","article-title":"Principles and challenges of genomewide DNA methylation analysis","volume":"11","author":"Laird","year":"2010","journal-title":"Nat. Rev. Genet."},{"key":"2023012710392670700_btt674-B12","doi-asserted-by":"crossref","first-page":"315","DOI":"10.1038\/nature08514","article-title":"Human DNA methylomes at base resolution show widespread epigenomic differences","volume":"462","author":"Lister","year":"2009","journal-title":"Nature"},{"key":"2023012710392670700_btt674-B13","doi-asserted-by":"crossref","first-page":"68","DOI":"10.1038\/nature09798","article-title":"Hotspots of aberrant epigenomic reprogramming in human induced pluripotent stem cells","volume":"471","author":"Lister","year":"2011","journal-title":"Nature"},{"key":"2023012710392670700_btt674-B14","doi-asserted-by":"crossref","DOI":"10.1002\/0471721182","volume-title":"Finite Mixture Models","author":"McLachlan","year":"2000"},{"key":"2023012710392670700_btt674-B15","doi-asserted-by":"crossref","first-page":"98","DOI":"10.1182\/blood-2012-03-420273","article-title":"Transgene-free iPSCs generated from small volume peripheral blood nonmobilized CD34+ cells","volume":"121","author":"Merling","year":"2013","journal-title":"Blood"},{"key":"2023012710392670700_btt674-B16","doi-asserted-by":"crossref","first-page":"369","DOI":"10.1186\/1471-2105-11-369","article-title":"HPeak: an HMM-based algorithm for defining read-enriched regions in ChIP-Seq data","volume":"11","author":"Qin","year":"2010","journal-title":"BMC Bioinformatics"},{"key":"2023012710392670700_btt674-B17","doi-asserted-by":"crossref","first-page":"597","DOI":"10.1038\/nrg1655","article-title":"DNA methylation and human disease","volume":"6","author":"Robertson","year":"2005","journal-title":"Nat. Rev. Genet."},{"key":"2023012710392670700_btt674-B18","doi-asserted-by":"crossref","first-page":"226","DOI":"10.1016\/j.ymeth.2009.05.003","article-title":"High-throughput bisulfite sequencing in mammalian genomes","volume":"48","author":"Smith","year":"2009","journal-title":"Methods"},{"key":"2023012710392670700_btt674-B19","doi-asserted-by":"crossref","first-page":"479","DOI":"10.1111\/1467-9868.00346","article-title":"A direct approach to false discovery rates","volume":"64","author":"Storey","year":"2002","journal-title":"J. R. Statist. Soc. Series B Stat. Methodol."},{"key":"2023012710392670700_btt674-B20","doi-asserted-by":"crossref","first-page":"2013","DOI":"10.1214\/aos\/1074290335","article-title":"The positive false discovery rate: a Bayesian interpretation and the q-value","volume":"31","author":"Storey","year":"2003","journal-title":"Ann. Stat."},{"key":"2023012710392670700_btt674-B21","doi-asserted-by":"crossref","first-page":"9440","DOI":"10.1073\/pnas.1530509100","article-title":"Statistical significance for genome-wide studies","volume":"100","author":"Storey","year":"2003","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2023012710392670700_btt674-B22","doi-asserted-by":"crossref","first-page":"901","DOI":"10.1198\/016214507000000545","article-title":"Oracle and adaptive compound decision rules for false discovery rate control","volume":"102","author":"Sun","year":"2007","journal-title":"J. Am. Stat. Assoc."},{"key":"2023012710392670700_btt674-B23","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1214\/aos\/1176346060","article-title":"On the convergence properties of the EM algorithm","volume":"11","author":"Wu","year":"1983","journal-title":"Ann. Stat."},{"key":"2023012710392670700_btt674-B24","doi-asserted-by":"crossref","first-page":"e21034","DOI":"10.1371\/journal.pone.0021034","article-title":"Statistical quantification of methylation levels by next-generation sequencing","volume":"6","author":"Wu","year":"2011","journal-title":"PLoS One"},{"key":"2023012710392670700_btt674-B25","doi-asserted-by":"crossref","first-page":"873","DOI":"10.1093\/bioinformatics\/btq057","article-title":"Fast and SNP-tolerant detection of complex variants and splicing in short reads","volume":"26","author":"Wu","year":"2010","journal-title":"Bioinformatics"},{"key":"2023012710392670700_btt674-B26","article-title":"Read-mapping using personalized diploid reference genome for RNA sequencing data reduced bias for detecting allele-specific expression","volume-title":"BIBMW 2012 Workshop on Data-Mining of Next Generation Sequencing","author":"Yuan","year":"2013"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/30\/2\/172\/48916047\/bioinformatics_30_2_172.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/30\/2\/172\/48916047\/bioinformatics_30_2_172.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,5,20]],"date-time":"2024-05-20T07:47:26Z","timestamp":1716191246000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/30\/2\/172\/227206"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2013,11,21]]},"references-count":26,"journal-issue":{"issue":"2","published-print":{"date-parts":[[2014,1,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btt674","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2014,1,15]]},"published":{"date-parts":[[2013,11,21]]}}}