{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,11]],"date-time":"2026-03-11T23:50:31Z","timestamp":1773273031310,"version":"3.50.1"},"reference-count":30,"publisher":"Oxford University Press (OUP)","issue":"24","license":[{"start":{"date-parts":[[2016,10,2]],"date-time":"2016-10-02T00:00:00Z","timestamp":1475366400000},"content-version":"vor","delay-in-days":403,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2015,12,15]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Motivation: Assigning RNA-seq reads to their transcript of origin is a fundamental task in transcript expression estimation. Where ambiguities in assignments exist due to transcripts sharing sequence, e.g. alternative isoforms or alleles, the problem can be solved through probabilistic inference. Bayesian methods have been shown to provide accurate transcript abundance estimates compared with competing methods. However, exact Bayesian inference is intractable and approximate methods such as Markov chain Monte Carlo and Variational Bayes (VB) are typically used. While providing a high degree of accuracy and modelling flexibility, standard implementations can be prohibitively slow for large datasets and complex transcriptome annotations.<\/jats:p>\n               <jats:p>Results: We propose a novel approximate inference scheme based on VB and apply it to an existing model of transcript expression inference from RNA-seq data. Recent advances in VB algorithmics are used to improve the convergence of the algorithm beyond the standard Variational Bayes Expectation Maximization algorithm. We apply our algorithm to simulated and biological datasets, demonstrating a significant increase in speed with only very small loss in accuracy of expression level estimation. We carry out a comparative study against seven popular alternative methods and demonstrate that our new algorithm provides excellent accuracy and inter-replicate consistency while remaining competitive in computation time.<\/jats:p>\n               <jats:p>Availability and implementation: The methods were implemented in R and C++, and are available as part of the BitSeq project at github.com\/BitSeq. The method is also available through the BitSeq Bioconductor package. The source code to reproduce all simulation results can be accessed via github.com\/BitSeq\/BitSeqVB_benchmarking.<\/jats:p>\n               <jats:p>Contact: \u00a0james.hensman@sheffield.ac.uk or panagiotis.papastamoulis@manchester.ac.uk or Magnus.Rattray@manchester.ac.uk<\/jats:p>\n               <jats:p>Supplementary information: \u00a0Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btv483","type":"journal-article","created":{"date-parts":[[2015,8,28]],"date-time":"2015-08-28T00:18:54Z","timestamp":1440721134000},"page":"3881-3889","source":"Crossref","is-referenced-by-count":33,"title":["Fast and accurate approximate inference of transcript expression from RNA-seq data"],"prefix":"10.1093","volume":"31","author":[{"given":"James","family":"Hensman","sequence":"first","affiliation":[{"name":"1 Sheffield Institute for Translational Neuroscience (SITraN), Sheffield, UK,"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Panagiotis","family":"Papastamoulis","sequence":"additional","affiliation":[{"name":"2 Faculty of Life Sciences,"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Peter","family":"Glaus","sequence":"additional","affiliation":[{"name":"3 School of Computer Science, The University of Manchester, Manchester, UK and"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Antti","family":"Honkela","sequence":"additional","affiliation":[{"name":"4 Helsinki Institute for Information Technology (HIIT), Department of Computer Science, University of Helsinki, Helsinki, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Magnus","family":"Rattray","sequence":"additional","affiliation":[{"name":"2 Faculty of Life Sciences,"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2015,8,26]]},"reference":[{"key":"2023051307183814000_btv483-B1","doi-asserted-by":"crossref","first-page":"251","DOI":"10.1162\/089976698300017746","article-title":"Natural gradient works efficiently in learning","volume":"10","author":"Amari","year":"1998","journal-title":"Neural Comput."},{"key":"2023051307183814000_btv483-B2","volume-title":"Pattern Recognition and Machine Learning","author":"Bishop","year":"2006"},{"key":"2023051307183814000_btv483-B3","article-title":"Near-optimal RNA-Seq quantification","volume":"arXiv","author":"Bray","year":"2015","journal-title":"arXiv (q-bio.QM)"},{"key":"2023051307183814000_btv483-B4","doi-asserted-by":"crossref","first-page":"D48","DOI":"10.1093\/nar\/gks1236","article-title":"Ensembl 2013","volume":"41","author":"Flicek","year":"2013","journal-title":"Nucleic Acids Res."},{"key":"2023051307183814000_btv483-B5","doi-asserted-by":"crossref","DOI":"10.1201\/9780429258480","volume-title":"Bayesian Data Analysis","author":"Gelman","year":"2003","edition":"2nd edn"},{"key":"2023051307183814000_btv483-B6","doi-asserted-by":"crossref","first-page":"1721","DOI":"10.1093\/bioinformatics\/bts260","article-title":"Identifying differentially expressed transcripts from RNA-seq data with biological variation","volume":"28","author":"Glaus","year":"2012","journal-title":"Bioinformatics"},{"key":"2023051307183814000_btv483-B7","article-title":"Fast variational inference in the conjugate exponential family","author":"Hensman","year":"2012","journal-title":"Adv. Neural Inf. Process. Syst. (NIPS)."},{"key":"2023051307183814000_btv483-B8","doi-asserted-by":"crossref","first-page":"383","DOI":"10.1109\/TPAMI.2014.2318711","article-title":"Fast nonparametric clustering of structured time-series","volume":"37","author":"Hensman","year":"2015","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell."},{"key":"2023051307183814000_btv483-B9","first-page":"3235","article-title":"Approximate Riemannian conjugate gradient learning for fixed-form variational Bayes","volume":"11","author":"Honkela","year":"2010","journal-title":"J. Mach. Learn. Res."},{"key":"2023051307183814000_btv483-B10","doi-asserted-by":"crossref","first-page":"1026","DOI":"10.1093\/bioinformatics\/btp113","article-title":"Statistical inferences for isoform expression in RNA-seq","volume":"25","author":"Jiang","year":"2009","journal-title":"Bioinformatics"},{"key":"2023051307183814000_btv483-B11","doi-asserted-by":"crossref","first-page":"1009","DOI":"10.1038\/nmeth.1528","article-title":"Analysis and design of RNA sequencing experiments for identifying isoform regulation","volume":"7","author":"Katz","year":"2010","journal-title":"Nat. Methods"},{"key":"2023051307183814000_btv483-B12","doi-asserted-by":"crossref","first-page":"996","DOI":"10.1101\/gr.229102","article-title":"The human genome browser at UCSC","volume":"6","author":"Kent","year":"2002","journal-title":"Genome Res."},{"key":"2023051307183814000_btv483-B13","doi-asserted-by":"crossref","first-page":"357","DOI":"10.1038\/nmeth.1923","article-title":"Fast gapped-read alignment with Bowtie 2","volume":"9","author":"Langmead","year":"2012","journal-title":"Nat. Methods"},{"key":"2023051307183814000_btv483-B14","doi-asserted-by":"crossref","first-page":"323","DOI":"10.1186\/1471-2105-12-323","article-title":"RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome","volume":"12","author":"Li","year":"2011","journal-title":"BMC Bioinformatics"},{"key":"2023051307183814000_btv483-B15","doi-asserted-by":"crossref","first-page":"493","DOI":"10.1093\/bioinformatics\/btp692","article-title":"RNA-Seq gene expression estimation with read mapping uncertainty","volume":"26","author":"Li","year":"2010","journal-title":"Bioinformatics"},{"key":"2023051307183814000_btv483-B16","doi-asserted-by":"crossref","first-page":"621","DOI":"10.1038\/nmeth.1226","article-title":"Mapping and quantifying mammalian transcriptomes by RNA-Seq","volume":"5","author":"Mortazavi","year":"2008","journal-title":"Nat. Methods"},{"key":"2023051307183814000_btv483-B17","doi-asserted-by":"crossref","first-page":"2292","DOI":"10.1093\/bioinformatics\/btt381","article-title":"TIGAR: transcript isoform abundance estimation method with gapped alignment of RNA-Seq data by variational Bayesian inference","volume":"18","author":"Nariai","year":"2013","journal-title":"Bioinformatics"},{"key":"2023051307183814000_btv483-B18","doi-asserted-by":"crossref","DOI":"10.1186\/1471-2164-15-S10-S5","article-title":"TIGAR2: sensitive and accurate estimation of transcript isoform expression with longer RNA-Seq reads","volume":"15","author":"Nariai","year":"2014","journal-title":"BMC Genomics"},{"key":"2023051307183814000_btv483-B19","doi-asserted-by":"crossref","first-page":"203","DOI":"10.1515\/sagmb-2013-0054","article-title":"Improved variational Bayes inference for transcript expression estimation","volume":"13","author":"Papastamoulis","year":"2014","journal-title":"Stat. Appl. Genet. Mol. Biol."},{"key":"2023051307183814000_btv483-B20","doi-asserted-by":"crossref","first-page":"97","DOI":"10.1016\/j.csda.2014.07.005","article-title":"On the estimation of mixtures of Poisson regression models with large number of components","volume":"93","author":"Papastamoulis","year":"2014","journal-title":"Comput. Stat. Data Anal."},{"key":"2023051307183814000_btv483-B21","doi-asserted-by":"crossref","first-page":"462","DOI":"10.1038\/nbt.2862","article-title":"Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms","volume":"32","author":"Patro","year":"2014","journal-title":"Nat. Biotechnol."},{"key":"2023051307183814000_btv483-B22","doi-asserted-by":"crossref","first-page":"71","DOI":"10.1038\/nmeth.2251","article-title":"Streaming fragment assignment for real-time analysis of sequencing experiments","volume":"10","author":"Roberts","year":"2013","journal-title":"Nat. Methods"},{"key":"2023051307183814000_btv483-B23","doi-asserted-by":"crossref","first-page":"309","DOI":"10.1214\/13-AOAS687","article-title":"Quantifying alternative splicing from paired-end RNA-sequencing data","volume":"8","author":"Rossell","year":"2014","journal-title":"Ann. Appl. Stat."},{"key":"2023051307183814000_btv483-B24","doi-asserted-by":"crossref","first-page":"903","DOI":"10.1038\/nbt.2957","article-title":"A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the sequencing quality control consortium","volume":"32","author":"SEQC\/MAQC-III Consortium","year":"2014","journal-title":"Nat. Biotechnol."},{"key":"2023051307183814000_btv483-B25","doi-asserted-by":"crossref","first-page":"320","DOI":"10.1186\/1471-2105-14-320","article-title":"Design of RNA splicing analysis null models for post hoc filtering of Drosophila head RNA-Seq data with the splicing analysis kit (Spanki)","volume":"14","author":"Sturgill","year":"2013","journal-title":"BMC Bioinformatics"},{"key":"2023051307183814000_btv483-B26","doi-asserted-by":"crossref","first-page":"1616","DOI":"10.1101\/gr.134445.111","article-title":"Deep sequencing of subcellular RNA fractions shows splicing to be predominantly co-transcriptional in the human genome but inefficient for lncRNAs","volume":"22","author":"Tilgner","year":"2012","journal-title":"Genome Res."},{"key":"2023051307183814000_btv483-B27","doi-asserted-by":"crossref","first-page":"516","DOI":"10.1038\/nbt.1621","article-title":"Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation","volume":"28","author":"Trapnell","year":"2010","journal-title":"Nat. Biotechnol."},{"key":"2023051307183814000_btv483-B28","doi-asserted-by":"crossref","first-page":"46","DOI":"10.1038\/nbt.2450","article-title":"Differential analysis of gene regulation at transcript resolution with RNA-seq","volume":"31","author":"Trapnell","year":"2013","journal-title":"Nat. Biotechnol."},{"key":"2023051307183814000_btv483-B29","doi-asserted-by":"crossref","first-page":"R13","DOI":"10.1186\/gb-2011-12-2-r13","article-title":"Haplotype and isoform specific expression estimation using multi-mapping RNA-seq reads","volume":"12","author":"Turro","year":"2011","journal-title":"Genome Biol."},{"key":"2023051307183814000_btv483-B30","doi-asserted-by":"crossref","first-page":"3150","DOI":"10.1093\/nar\/gkl396","article-title":"An expectation-maximization algorithm for probabilistic reconstructions of full-length isoforms from splice graphs","volume":"34","author":"Xing","year":"2006","journal-title":"Nucleic Acids Res."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/24\/3881\/50306921\/bioinformatics_31_24_3881.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/24\/3881\/50306921\/bioinformatics_31_24_3881.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,5,13]],"date-time":"2023-05-13T07:20:37Z","timestamp":1683962437000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/31\/24\/3881\/196984"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2015,8,26]]},"references-count":30,"journal-issue":{"issue":"24","published-print":{"date-parts":[[2015,12,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btv483","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2015,12,15]]},"published":{"date-parts":[[2015,8,26]]}}}