{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,2,22]],"date-time":"2025-02-22T00:45:07Z","timestamp":1740185107438,"version":"3.37.3"},"reference-count":24,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2018,8,23]],"date-time":"2018-08-23T00:00:00Z","timestamp":1534982400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100013429","name":"Delaware INBRE","doi-asserted-by":"publisher","award":["NIGMS P20GM103446"],"award-info":[{"award-number":["NIGMS P20GM103446"]}],"id":[{"id":"10.13039\/100013429","id-type":"DOI","asserted-by":"publisher"}]},{"name":"State of Delaware and the Delaware Biotechnology Institute"},{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["IOS1545780"],"award-info":[{"award-number":["IOS1545780"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"NIH","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,3,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Several methods have been proposed for the paired RNA-seq analysis. However, many of them do not consider the heterogeneity in treatment effect among pairs that can naturally arise in real data. In addition, it has been reported in literature that the false discovery rate (FDR) control of some popular methods has been problematic. In this paper, we present a full hierarchical Bayesian model for the paired RNA-seq count data that accounts for variation of treatment effects among pairs and controls the FDR through the posterior expected FDR.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>Our simulation studies show that most competing methods can have highly inflated FDR for small to moderate sample sizes while PairedFB is able to control FDR close to the nominal levels. Furthermore, PairedFB has overall better performance in ranking true differentially expressed genes (DEGs) on the top than others, especially when the sample size gets bigger or when the heterogeneity level of treatment effects is high. In addition, PairedFB can be applied to identify the biologically significant DEGs with controlled FDR. The real data analysis also indicates PairedFB tends to find more biologically relevant genes even when the sample size is small. PairedFB is also shown to be robust with respect to the model misspecification in terms of its relative performance compared to others.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>Software to implement this method (PairedFB) can be downloaded at: https:\/\/sites.google.com\/a\/udel.edu\/qiujing\/publication.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Supplementary information<\/jats:title>\n                  <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty731","type":"journal-article","created":{"date-parts":[[2018,8,22]],"date-time":"2018-08-22T19:19:20Z","timestamp":1534965560000},"page":"787-797","source":"Crossref","is-referenced-by-count":5,"title":["PairedFB: a full hierarchical Bayesian model for paired RNA-seq data with heterogeneous treatment effects"],"prefix":"10.1093","volume":"35","author":[{"given":"Yuanyuan","family":"Bian","sequence":"first","affiliation":[{"name":"Department of Statistics, University of Missouri, Columbia, MO, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Chong","family":"He","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of Missouri, Columbia, MO, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jie","family":"Hou","sequence":"additional","affiliation":[{"name":"Department of Electrical Engineering and Computer Science, University of Missouri, Columbia, MO, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jianlin","family":"Cheng","sequence":"additional","affiliation":[{"name":"Department of Electrical Engineering and Computer Science, University of Missouri, Columbia, MO, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jing","family":"Qiu","sequence":"additional","affiliation":[{"name":"Department of Applied Economics and Statistics, University of Delaware, Newark, DE, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2018,8,23]]},"reference":[{"key":"2023020108345686100_bty731-B1","first-page":"555","volume-title":"Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables","author":"Abramowitz","year":"1964"},{"key":"2023020108345686100_bty731-B2","doi-asserted-by":"crossref","first-page":"446","DOI":"10.1111\/j.1541-0420.2005.00499.x","article-title":"Bayesian covariance selection in generalized linear mixed models","volume":"62","author":"Cai","year":"2006","journal-title":"Biometrics"},{"key":"2023020108345686100_bty731-B3","doi-asserted-by":"crossref","first-page":"110.","DOI":"10.1186\/1471-2105-14-110","article-title":"Differential expression analysis for paired RNA-seq data","volume":"14","author":"Chung","year":"2013","journal-title":"BMC Bioinformatics"},{"key":"2023020108345686100_bty731-B4","doi-asserted-by":"crossref","first-page":"87","DOI":"10.1515\/sagmb-2015-0011","article-title":"What if we ignore the random effects when analyzing RNA-seq data in a multifactor experiment","volume":"15","author":"Cui","year":"2016","journal-title":"Stat. Appl. Genet. Mol. Biol"},{"key":"2023020108345686100_bty731-B5","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v040.i08","article-title":"Rcpp: seamless R and C++ integration","volume":"40","author":"Eddelbuettel","year":"2011","journal-title":"J. Stat. Softw"},{"key":"2023020108345686100_bty731-B6","doi-asserted-by":"crossref","first-page":"e0170632.","DOI":"10.1371\/journal.pone.0170632","article-title":"A comparison of RNA-seq results from paired formalin-fixed paraffin-embedded and fresh-frozen glioblastoma tissue samples","volume":"12","author":"Esteve-Codina","year":"2017","journal-title":"PLoS One"},{"volume-title":"GNU Scientific Library Reference Manual.","year":"2009","author":"Gough","key":"2023020108345686100_bty731-B7"},{"key":"2023020108345686100_bty731-B8","doi-asserted-by":"crossref","first-page":"12335.","DOI":"10.1038\/srep12335","article-title":"Robust gene expression and mutation analyses of RNA-sequencing of formalin-fixed diagnostic tumor samples","volume":"5","author":"Graw","year":"2015","journal-title":"Sci. Rep"},{"key":"2023020108345686100_bty731-B9","doi-asserted-by":"crossref","first-page":"S2.","DOI":"10.1186\/1471-2164-14-S8-S2","article-title":"Evaluation of read count based RNAseq analysis methods","volume":"14","author":"Guo","year":"2013","journal-title":"BMC Genomics"},{"key":"2023020108345686100_bty731-B10","doi-asserted-by":"crossref","first-page":"135.","DOI":"10.1186\/1471-2105-14-135","article-title":"Empirical Bayesian analysis of paired high-throughput sequencing data with a beta-binomial distribution","volume":"14","author":"Hardcastle","year":"2013","journal-title":"BMC Bioinformatics"},{"key":"2023020108345686100_bty731-B11","doi-asserted-by":"crossref","first-page":"R29.","DOI":"10.1186\/gb-2014-15-2-r29","article-title":"Voom: precision weights unlock linear model analysis tools for RNA-seq read counts","volume":"15","author":"Law","year":"2014","journal-title":"Genome Biol"},{"key":"2023020108345686100_bty731-B12","doi-asserted-by":"crossref","first-page":"555","DOI":"10.1007\/s13253-015-0227-0","article-title":"A semi-parametric Bayesian approach for differential expression analysis of RNA-seq data","volume":"20","author":"Liu","year":"2015","journal-title":"J. Agric. Biol. Environ. Stat"},{"key":"2023020108345686100_bty731-B13","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol"},{"key":"2023020108345686100_bty731-B14","doi-asserted-by":"crossref","first-page":"4288","DOI":"10.1093\/nar\/gks042","article-title":"Differential expression analysis of multifactor RNA-seq experiments with respect to biological variation","volume":"40","author":"McCarthy","year":"2012","journal-title":"Nucleic Acids Res"},{"key":"2023020108345686100_bty731-B15","doi-asserted-by":"crossref","first-page":"155","DOI":"10.1093\/biostatistics\/5.2.155","article-title":"Detecting differential gene expression with a semiparametric hierarchical mixture method","volume":"5","author":"Newton","year":"2004","journal-title":"Biostatistics"},{"key":"2023020108345686100_bty731-B16","first-page":"7","article-title":"CODA: convergence diagnosis and output analysis for mcmc","volume":"6","author":"Plummer","year":"2006","journal-title":"R News"},{"key":"2023020108345686100_bty731-B17","doi-asserted-by":"crossref","first-page":"R25.","DOI":"10.1186\/gb-2010-11-3-r25","article-title":"A scaling normalization method for differential expression analysis of RNA-seq data","volume":"11","author":"Robinson","year":"2010","journal-title":"Genome Biol"},{"key":"2023020108345686100_bty731-B18","doi-asserted-by":"crossref","first-page":"2881","DOI":"10.1093\/bioinformatics\/btm453","article-title":"Moderated statistical tests for assessing differences in tag abundance","volume":"23","author":"Robinson","year":"2007","journal-title":"Bioinformatics"},{"key":"2023020108345686100_bty731-B19","doi-asserted-by":"crossref","first-page":"2587","DOI":"10.1214\/10-AOS792","article-title":"Bayes and empirical-bayes multiplicity adjustment in the variable-selection problem","volume":"38","author":"Scott","year":"2010","journal-title":"Ann. Stat"},{"key":"2023020108345686100_bty731-B20","doi-asserted-by":"crossref","first-page":"91.","DOI":"10.1186\/1471-2105-14-91","article-title":"A comparison of methods for differential expression analysis of RNA-seq data","volume":"14","author":"Soneson","year":"2013","journal-title":"BMC Bioinformatics"},{"key":"2023020108345686100_bty731-B21","doi-asserted-by":"crossref","first-page":"15545","DOI":"10.1073\/pnas.0506580102","article-title":"Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles","volume":"102","author":"Subramanian","year":"2005","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023020108345686100_bty731-B22","doi-asserted-by":"crossref","first-page":"e9317.","DOI":"10.1371\/journal.pone.0009317","article-title":"Tumor transcriptome sequencing reveals allelic expression imbalances associated with copy number alterations","volume":"5","author":"Tuch","year":"2010","journal-title":"PLoS One"},{"key":"2023020108345686100_bty731-B23","doi-asserted-by":"crossref","first-page":"113","DOI":"10.1093\/biostatistics\/kxs031","article-title":"Bayesian analysis of RNA sequencing data by estimating multiple shrinkage priors","volume":"14","author":"Van De Wiel","year":"2013","journal-title":"Biostatistics"},{"key":"2023020108345686100_bty731-B24","doi-asserted-by":"crossref","first-page":"e3215.","DOI":"10.1371\/journal.pone.0003215","article-title":"The evolving transcriptome of head and neck squamous cell carcinoma: a systematic review","volume":"3","author":"Yu","year":"2008","journal-title":"PLoS One"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/5\/787\/48965892\/bioinformatics_35_5_787.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/5\/787\/48965892\/bioinformatics_35_5_787.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,2,1]],"date-time":"2023-02-01T19:40:17Z","timestamp":1675280417000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/5\/787\/5078474"}},"subtitle":[],"editor":[{"given":"Bonnie","family":"Berger","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2018,8,23]]},"references-count":24,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2019,3,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty731","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"type":"print","value":"1367-4803"},{"type":"electronic","value":"1367-4811"}],"subject":[],"published-other":{"date-parts":[[2019,3,1]]},"published":{"date-parts":[[2018,8,23]]}}}