{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,28]],"date-time":"2026-07-28T03:15:37Z","timestamp":1785208537839,"version":"3.55.0"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1008585","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,1,22]],"date-time":"2021-01-22T00:00:00Z","timestamp":1611273600000}}],"reference-count":40,"publisher":"Public Library of Science (PLoS)","issue":"1","license":[{"start":{"date-parts":[[2021,1,11]],"date-time":"2021-01-11T00:00:00Z","timestamp":1610323200000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000051","name":"National Human Genome Research Institute","doi-asserted-by":"publisher","award":["R01 HG009937"],"award-info":[{"award-number":["R01 HG009937"]}],"id":[{"id":"10.13039\/100000051","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100008982","name":"National Science Foundation","doi-asserted-by":"publisher","award":["NSF CCF-1750472"],"award-info":[{"award-number":["NSF CCF-1750472"]}],"id":[{"id":"10.13039\/501100008982","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100008982","name":"National Science Foundation","doi-asserted-by":"publisher","award":["CNS-1763680"],"award-info":[{"award-number":["CNS-1763680"]}],"id":[{"id":"10.13039\/501100008982","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>Experimental single-cell approaches are becoming widely used for many purposes, including investigation of the dynamic behaviour of developing biological systems. Consequently, a large number of computational methods for extracting dynamic information from such data have been developed. One example is RNA velocity analysis, in which spliced and unspliced RNA abundances are jointly modeled in order to infer a \u2018direction of change\u2019 and thereby a future state for each cell in the gene expression space. Naturally, the accuracy and interpretability of the inferred RNA velocities depend crucially on the correctness of the estimated abundances. Here, we systematically compare five widely used quantification tools, in total yielding thirteen different quantification approaches, in terms of their estimates of spliced and unspliced RNA abundances in five experimental droplet scRNA-seq data sets. We show that there are substantial differences between the quantifications obtained from different tools, and identify typical genes for which such discrepancies are observed. We further show that these abundance differences propagate to the downstream analysis, and can have a large effect on estimated velocities as well as the biological interpretation. Our results highlight that abundance quantification is a crucial aspect of the RNA velocity analysis workflow, and that both the definition of the genomic features of interest and the quantification algorithm itself require careful consideration.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1008585","type":"journal-article","created":{"date-parts":[[2021,1,12]],"date-time":"2021-01-12T17:34:48Z","timestamp":1610472888000},"page":"e1008585","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":69,"title":["Preprocessing choices affect RNA velocity results for droplet scRNA-seq data"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-3833-2169","authenticated-orcid":true,"given":"Charlotte","family":"Soneson","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9798-2079","authenticated-orcid":true,"given":"Avi","family":"Srivastava","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8463-1675","authenticated-orcid":true,"given":"Rob","family":"Patro","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2269-4934","authenticated-orcid":true,"given":"Michael B.","family":"Stadler","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2021,1,11]]},"reference":[{"issue":"1","key":"pcbi.1008585.ref001","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1038\/nmeth.2801","article-title":"Method of the year 2013","volume":"11","year":"2014","journal-title":"Nat Methods"},{"issue":"1","key":"pcbi.1008585.ref002","doi-asserted-by":"crossref","first-page":"25","DOI":"10.1038\/nmeth.2769","article-title":"The promise of single-cell sequencing","volume":"11","author":"J Eberwine","year":"2014","journal-title":"Nat Methods"},{"issue":"1","key":"pcbi.1008585.ref003","doi-asserted-by":"crossref","first-page":"171","DOI":"10.1038\/nprot.2014.006","article-title":"Full-length RNA-seq from single cells using Smart-seq2","volume":"9","author":"S Picelli","year":"2014","journal-title":"Nat Protoc"},{"key":"pcbi.1008585.ref004","doi-asserted-by":"crossref","first-page":"14049","DOI":"10.1038\/ncomms14049","article-title":"Massively parallel digital transcriptional profiling of single cells","volume":"8","author":"GXY Zheng","year":"2017","journal-title":"Nat Commun"},{"key":"pcbi.1008585.ref005","doi-asserted-by":"crossref","first-page":"45","DOI":"10.1007\/978-1-4939-9240-9_4","article-title":"CEL-Seq2-Single-Cell RNA Sequencing by Multiplexed Linear Amplification","volume":"1979","author":"I Yanai","year":"2019","journal-title":"Methods Mol Biol"},{"key":"pcbi.1008585.ref006","doi-asserted-by":"crossref","first-page":"129","DOI":"10.1016\/j.copbio.2019.03.001","article-title":"Beyond bulk: a review of single cell transcriptomics methodologies and applications","volume":"58","author":"A Kulkarni","year":"2019","journal-title":"Curr Opin Biotechnol"},{"key":"pcbi.1008585.ref007","doi-asserted-by":"crossref","first-page":"70","DOI":"10.1016\/j.copbio.2019.12.005","article-title":"The emergence and promise of single-cell temporal-omics approaches","volume":"63","author":"AR Lederer","year":"2020","journal-title":"Curr Opin Biotechnol"},{"key":"pcbi.1008585.ref008","doi-asserted-by":"crossref","first-page":"547","DOI":"10.1038\/s41587-019-0071-9","article-title":"A comparison of single-cell trajectory inference methods","volume":"37","author":"W Saelens","year":"2019","journal-title":"Nat Biotechnol"},{"key":"pcbi.1008585.ref009","doi-asserted-by":"crossref","first-page":"494","DOI":"10.1038\/s41586-018-0414-6","article-title":"RNA velocity of single cells","volume":"560","author":"G La Manno","year":"2018","journal-title":"Nature"},{"key":"pcbi.1008585.ref010","doi-asserted-by":"crossref","first-page":"529","DOI":"10.1038\/msb.2011.62","article-title":"Coupled pre-mRNA and mRNA dynamics unveil operational strategies underlying transcriptional responses to stimuli","volume":"7","author":"A Zeisel","year":"2011","journal-title":"Mol Syst Biol"},{"key":"pcbi.1008585.ref011","doi-asserted-by":"crossref","first-page":"722","DOI":"10.1038\/nbt.3269","article-title":"Analysis of intronic and exonic reads in RNA-seq data characterizes transcriptional and post-transcriptional regulation","volume":"33","author":"D Gaidatzis","year":"2015","journal-title":"Nat Biotechnol"},{"key":"pcbi.1008585.ref012","doi-asserted-by":"crossref","first-page":"e89673","DOI":"10.1371\/journal.pone.0089673","article-title":"SnapShot-Seq: A method for extracting genome-wide, in Vivo mRNA dynamics from a single total RNA sample","volume":"9","author":"JM Gray","year":"2014","journal-title":"PLoS One"},{"key":"pcbi.1008585.ref013","article-title":"Generalizing RNA velocity to transient cell states through dynamical modeling","author":"V Bergen","year":"2020","journal-title":"Nat Biotechnol"},{"key":"pcbi.1008585.ref014","doi-asserted-by":"crossref","first-page":"4472","DOI":"10.1093\/bioinformatics\/btz279","article-title":"The Barcode, UMI, Set format and BUStools","volume":"35","author":"P Melsted","year":"2019","journal-title":"Bioinformatics"},{"key":"pcbi.1008585.ref015","article-title":"Modular and efficient pre-processing of single-cell RNA-seq","author":"P Melsted","year":"2019","journal-title":"bioRxiv"},{"key":"pcbi.1008585.ref016","doi-asserted-by":"crossref","first-page":"65","DOI":"10.1186\/s13059-019-1670-y","article-title":"Alevin efficiently estimates accurate gene abundances from dscRNA-seq data","volume":"20","author":"A Srivastava","year":"2019","journal-title":"Genome Biol"},{"key":"pcbi.1008585.ref017","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","article-title":"STAR: ultrafast universal RNA-seq aligner","volume":"29","author":"A Dobin","year":"2013","journal-title":"Bioinformatics"},{"issue":"3","key":"pcbi.1008585.ref018","doi-asserted-by":"crossref","first-page":"491","DOI":"10.1101\/gr.209601.116","article-title":"UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy","volume":"27","author":"T Smith","year":"2017","journal-title":"Genome Res"},{"key":"pcbi.1008585.ref019","doi-asserted-by":"crossref","first-page":"giy059","DOI":"10.1093\/gigascience\/giy059","article-title":"zUMIs\u2014A fast and flexible pipeline to process RNA sequencing data with UMIs","volume":"7","author":"S Parekh","year":"2018","journal-title":"GigaScience"},{"issue":"1","key":"pcbi.1008585.ref020","doi-asserted-by":"crossref","first-page":"78","DOI":"10.1186\/s13059-018-1449-6","article-title":"dropEst: pipeline for accurate estimation of molecular counts in droplet-based single-cell RNA-seq experiments","volume":"19","author":"V Petukhov","year":"2018","journal-title":"Genome Biol"},{"issue":"8","key":"pcbi.1008585.ref021","doi-asserted-by":"crossref","first-page":"e1006361","DOI":"10.1371\/journal.pcbi.1006361","article-title":"scPipe: A flexible R\/Bioconductor preprocessing pipeline for single-cell RNA-sequencing data","volume":"14","author":"L Tian","year":"2018","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1008585.ref022","doi-asserted-by":"crossref","first-page":"dev173849","DOI":"10.1242\/dev.173849","article-title":"Comprehensive single cell mRNA profiling reveals a detailed roadmap for pancreatic endocrinogenesis","volume":"146","author":"A Bastidas-Ponce","year":"2019","journal-title":"Development"},{"key":"pcbi.1008585.ref023","doi-asserted-by":"crossref","first-page":"290","DOI":"10.1038\/s41593-017-0056-2","article-title":"Conserved properties of dentate gyrus neurogenesis across postnatal development revealed by single-cell RNA sequencing","volume":"21","author":"H Hochgerner","year":"2018","journal-title":"Nat Neurosci"},{"key":"pcbi.1008585.ref024","doi-asserted-by":"crossref","first-page":"1650","DOI":"10.1016\/j.celrep.2018.10.026","article-title":"The Mammalian Spermatogenesis Single-Cell Transcriptome, from Spermatogonial Stem Cells to Spermatids","volume":"25","author":"BP Hermann","year":"2018","journal-title":"Cell Rep"},{"key":"pcbi.1008585.ref025","doi-asserted-by":"crossref","first-page":"4169","DOI":"10.1038\/s41467-019-12054-3","article-title":"Cell type-specific transcriptional programs in mouse prefrontal cortex during adolescence and addiction","volume":"10","author":"A Bhattacherjee","year":"2019","journal-title":"Nat Commun"},{"key":"pcbi.1008585.ref026","doi-asserted-by":"crossref","first-page":"1696","DOI":"10.1038\/s41593-019-0491-3","article-title":"Single-cell transcriptomic profiling of the aging mouse brain","volume":"22","author":"M Ximerakis","year":"2019","journal-title":"Nat Neurosci"},{"issue":"D1","key":"pcbi.1008585.ref027","doi-asserted-by":"crossref","first-page":"D766","DOI":"10.1093\/nar\/gky955","article-title":"GENCODE reference annotation for the human and mouse genomes","volume":"47","author":"A Frankish","year":"2019","journal-title":"Nucleic Acids Res"},{"key":"pcbi.1008585.ref028","unstructured":"Moses L, Pachter L. BUSpaRse: kallisto | bustools R utilities; 2019. Available from http:\/\/bioconductor.org\/packages\/BUSpaRse\/"},{"key":"pcbi.1008585.ref029","unstructured":"Pag\u00e8s H. BSgenome: Software infrastructure for efficient representation of full genomes and their SNPs; 2019. Available from http:\/\/bioconductor.org\/packages\/BSgenome\/"},{"key":"pcbi.1008585.ref030","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1003118","article-title":"Software for Computing and Annotating Genomic Ranges","volume":"9","author":"M Lawrence","year":"2013","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1008585.ref031","doi-asserted-by":"crossref","first-page":"417","DOI":"10.1038\/nmeth.4197","article-title":"Salmon provides fast and bias-aware quantification of transcript expression","volume":"14","author":"R Patro","year":"2017","journal-title":"Nat Methods"},{"issue":"1","key":"pcbi.1008585.ref032","doi-asserted-by":"crossref","first-page":"239","DOI":"10.1186\/s13059-020-02151-8","article-title":"Alignment and mapping methodology influence transcript abundance estimation","volume":"21","author":"A Srivastava","year":"2020","journal-title":"Genome Biol"},{"key":"pcbi.1008585.ref033","doi-asserted-by":"crossref","first-page":"525","DOI":"10.1038\/nbt.3519","article-title":"Near-optimal probabilistic RNA-seq quantification","volume":"34","author":"NL Bray","year":"2016","journal-title":"Nat Biotechnol"},{"key":"pcbi.1008585.ref034","unstructured":"Lun A, Risso D. SingleCellExperiment: S4 Classes for Single Cell Data; 2019. Available from http:\/\/bioconductor.org\/packages\/SingleCellExperiment\/."},{"key":"pcbi.1008585.ref035","article-title":"Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R","author":"DJ McCarthy","year":"2017","journal-title":"Bioinformatics"},{"issue":"29","key":"pcbi.1008585.ref036","doi-asserted-by":"crossref","first-page":"861","DOI":"10.21105\/joss.00861","article-title":"UMAP: Uniform Manifold Approximation and Projection","volume":"3","author":"L McInnes","year":"2018","journal-title":"J. Open Source Softw"},{"key":"pcbi.1008585.ref037","doi-asserted-by":"crossref","first-page":"841","DOI":"10.1093\/bioinformatics\/btq033","article-title":"BEDTools: a flexible suite of utilities for comparing genomic features","volume":"26","author":"AR Quinlan","year":"2010","journal-title":"Bioinformatics"},{"key":"pcbi.1008585.ref038","doi-asserted-by":"crossref","first-page":"2204","DOI":"10.1093\/bioinformatics\/btq351","article-title":"BigWig and BigBed: enabling browsing of large distributed datasets","volume":"26","author":"WJ Kent","year":"2010","journal-title":"Bioinformatics"},{"key":"pcbi.1008585.ref039","doi-asserted-by":"crossref","first-page":"335","DOI":"10.1007\/978-1-4939-3578-9_16","volume-title":"Statistical Genomics: Methods and Protocols","author":"F Hahne","year":"2016"},{"key":"pcbi.1008585.ref040","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1186\/s13059-017-1382-0","article-title":"SCANPY: large-scale single-cell gene expression data analysis","volume":"19","author":"FA Wolf","year":"2018","journal-title":"Genome Biol"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1008585","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,1,22]],"date-time":"2021-01-22T00:00:00Z","timestamp":1611273600000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008585","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,1,22]],"date-time":"2021-01-22T18:05:46Z","timestamp":1611338746000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1008585"}},"subtitle":[],"editor":[{"given":"Min","family":"Li","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2021,1,11]]},"references-count":40,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2021,1,11]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1008585","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2020.03.13.990069","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,1,11]]}}}