{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,23]],"date-time":"2026-08-23T15:44:16Z","timestamp":1787499856212,"version":"build-2736575974"},"reference-count":43,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2021,12,15]],"date-time":"2021-12-15T00:00:00Z","timestamp":1639526400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"name":"Computational Systems Biology Graduate School","award":["GRK1772"],"award-info":[{"award-number":["GRK1772"]}]},{"DOI":"10.13039\/501100001659","name":"Deutsche Forschungsgemeinschaft","doi-asserted-by":"crossref","id":[{"id":"10.13039\/501100001659","id-type":"DOI","asserted-by":"crossref"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2022,2,7]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Single-cell RNA sequencing determines RNA copy numbers per cell for a given gene. However, technical noise poses the question how observed distributions (output) are connected to their cellular distributions (input).<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We model a single-cell RNA sequencing setup consisting of PCR amplification and sequencing, and derive probability distribution functions for the output distribution given an input distribution. We provide copy number distributions arising from single transcripts during PCR amplification with exact expressions for mean and variance. We prove that the coefficient of variation of the output of sequencing is always larger than that of the input distribution. Experimental data reveals the variance and mean of the input distribution to obey characteristic relations, which we specifically determine for a HeLa dataset. We can calculate as many moments of the input distribution as are known of the output distribution (up to all). This, in principle, completely determines the input from the output distribution.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and Implementation<\/jats:title>\n                    <jats:p>Source code freely available at https:\/\/github.com\/danielschw188\/InputOutputSCRNASeq.<\/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\/btab841","type":"journal-article","created":{"date-parts":[[2021,12,12]],"date-time":"2021-12-12T15:06:19Z","timestamp":1639321579000},"page":"1336-1343","source":"Crossref","is-referenced-by-count":2,"title":["On the relation between input and output distributions of scRNA-seq experiments"],"prefix":"10.1093","volume":"38","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-2825-3352","authenticated-orcid":false,"given":"Daniel","family":"Schwabe","sequence":"first","affiliation":[{"name":"Mathematical Cell Physiology, Max Delbr\u00fcck Center for Molecular Medicine in the Helmholtz Association , 13125 Berlin, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7137-1114","authenticated-orcid":false,"given":"Martin","family":"Falcke","sequence":"additional","affiliation":[{"name":"Mathematical Cell Physiology, Max Delbr\u00fcck Center for Molecular Medicine in the Helmholtz Association , 13125 Berlin, Germany"},{"name":"Department of Physics, Humboldt University Berlin , 12489 Berlin, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2021,12,15]]},"reference":[{"key":"2023020108573601700_btab841-B1","doi-asserted-by":"crossref","first-page":"R106","DOI":"10.1186\/gb-2010-11-10-r106","article-title":"Differential expression analysis for sequence count data","volume":"11","author":"Anders","year":"2010","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B2","doi-asserted-by":"crossref","first-page":"584","DOI":"10.1038\/nmeth.4263","article-title":"SCnorm: robust normalization of single-cell RNA-seq data","volume":"14","author":"Bacher","year":"2017","journal-title":"Nat. Methods"},{"key":"2023020108573601700_btab841-B3","doi-asserted-by":"crossref","first-page":"1388","DOI":"10.1101\/gr.3820805","article-title":"Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels","volume":"15","author":"Bengtsson","year":"2005","journal-title":"Genome Res"},{"key":"2023020108573601700_btab841-B4","doi-asserted-by":"crossref","first-page":"4996","DOI":"10.1016\/j.ces.2010.05.046","article-title":"Efficiency of the polymerase chain reaction","volume":"65","author":"Booth","year":"2010","journal-title":"Chem. Eng. Sci"},{"key":"2023020108573601700_btab841-B5","doi-asserted-by":"crossref","first-page":"1008","DOI":"10.1038\/s41587-021-00875-x","article-title":"Bayesian inference of gene expression states from single-cell RNA-seq data","volume":"39","author":"Breda","year":"2021","journal-title":"Nat. Biotechnol"},{"key":"2023020108573601700_btab841-B6","doi-asserted-by":"crossref","first-page":"411","DOI":"10.1038\/nbt.4096","article-title":"Integrating single-cell transcriptomic data across different conditions, technologies, and species","volume":"36","author":"Butler","year":"2018","journal-title":"Nat. Biotechnol"},{"key":"2023020108573601700_btab841-B7","doi-asserted-by":"crossref","first-page":"183","DOI":"10.1186\/s13059-020-02103-2","article-title":"Bayesian model selection reveals biological origins of zero inflation in single-cell transcriptomics","volume":"21","author":"Choi","year":"2020","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B8","doi-asserted-by":"crossref","first-page":"1183","DOI":"10.1126\/science.1070919","article-title":"Stochastic gene expression in a single cell","volume":"297","author":"Elowitz","year":"2002","journal-title":"Science"},{"key":"2023020108573601700_btab841-B9","doi-asserted-by":"crossref","first-page":"390","DOI":"10.1038\/s41467-018-07931-2","article-title":"Single-cell RNA-seq denoising using a deep count autoencoder","volume":"10","author":"Eraslan","year":"2019","journal-title":"Nat. Commun"},{"key":"2023020108573601700_btab841-B10","doi-asserted-by":"crossref","first-page":"296","DOI":"10.1186\/s13059-019-1874-1","article-title":"Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression","volume":"20","author":"Hafemeister","year":"2019","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B11","doi-asserted-by":"crossref","first-page":"421","DOI":"10.1038\/nbt.4091","article-title":"Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors","volume":"36","author":"Haghverdi","year":"2018","journal-title":"Nat. Biotechnol"},{"key":"2023020108573601700_btab841-B12","doi-asserted-by":"crossref","first-page":"422","DOI":"10.1186\/1471-2105-11-422","article-title":"baySeq: empirical Bayesian methods for identifying differential expression in sequence count data","volume":"11","author":"Hardcastle","year":"2010","journal-title":"BMC Bioinformatics"},{"key":"2023020108573601700_btab841-B13","volume-title":"The Theory of Branching Process","author":"Harris","year":"1964"},{"key":"2023020108573601700_btab841-B14","doi-asserted-by":"crossref","first-page":"562","DOI":"10.1093\/biostatistics\/kxx053","article-title":"Missing data and technical variability in single-cell RNA-sequencing experiments","volume":"19","author":"Hicks","year":"2018","journal-title":"Biostatistics"},{"key":"2023020108573601700_btab841-B15","doi-asserted-by":"crossref","first-page":"218","DOI":"10.1186\/s13059-020-02132-x","article-title":"A systematic evaluation of single-cell RNA-sequencing imputation methods","volume":"21","author":"Hou","year":"2020","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B16","doi-asserted-by":"crossref","first-page":"539","DOI":"10.1038\/s41592-018-0033-z","article-title":"SAVER: gene expression recovery for single-cell RNA sequencing","volume":"15","author":"Huang","year":"2018","journal-title":"Nat. Methods"},{"key":"2023020108573601700_btab841-B17","volume-title":"Stochastic Processes in Physics and Chemistry. North-Holland Personal Library","author":"Kampen","year":"2007","edition":"3rd edn"},{"key":"2023020108573601700_btab841-B18","doi-asserted-by":"crossref","first-page":"9253","DOI":"10.1073\/pnas.86.23.9253","article-title":"Fidelity of DNA polymerases in DNA amplification","volume":"86","author":"Keohavong","year":"1989","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023020108573601700_btab841-B19","doi-asserted-by":"crossref","first-page":"414","DOI":"10.1038\/335414a0","article-title":"Amplification and analysis of DNA sequences in single human sperm and diploid cells","volume":"335","author":"Li","year":"1988","journal-title":"Nature"},{"key":"2023020108573601700_btab841-B20","doi-asserted-by":"crossref","first-page":"997","DOI":"10.1038\/s41467-018-03405-7","article-title":"An accurate and robust imputation method scimpute for single-cell RNA-seq data","volume":"9","author":"Li","year":"2018","journal-title":"Nat. Commun"},{"key":"2023020108573601700_btab841-B21","doi-asserted-by":"crossref","first-page":"1053","DOI":"10.1038\/s41592-018-0229-2","article-title":"Deep generative modeling for single-cell transcriptomics","volume":"15","author":"Lopez","year":"2018","journal-title":"Nat. Methods"},{"key":"2023020108573601700_btab841-B22","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":"2023020108573601700_btab841-B23","doi-asserted-by":"crossref","first-page":"75","DOI":"10.1186\/s13059-016-0947-7","article-title":"Pooling across cells to normalize single-cell RNA sequencing data with many zero counts","volume":"17","author":"Lun","year":"2016","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B24","doi-asserted-by":"crossref","first-page":"181","DOI":"10.1007\/s11357-021-00339-9","article-title":"Cell-to-cell variation in gene expression and the aging process","volume":"43","author":"Mendenhall","year":"2021","journal-title":"Geroscience"},{"key":"2023020108573601700_btab841-B25","doi-asserted-by":"crossref","first-page":"14","DOI":"10.3390\/cells9010014","article-title":"Single-cell expression variability implies cell function","volume":"9","author":"Osorio","year":"2019","journal-title":"Cells"},{"key":"2023020108573601700_btab841-B26","doi-asserted-by":"crossref","first-page":"101","DOI":"10.1016\/S0006-3495(96)79205-6","article-title":"Theoretical uncertainty of measurements using quantitative polymerase chain reaction","volume":"71","author":"Peccoud","year":"1996","journal-title":"Biophys. J"},{"key":"2023020108573601700_btab841-B27","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1038\/s41467-017-02554-5","article-title":"A general and flexible method for signal extraction from single-cell RNA-seq data","volume":"9","author":"Risso","year":"2018","journal-title":"Nat. Commun"},{"key":"2023020108573601700_btab841-B28","doi-asserted-by":"crossref","first-page":"969","DOI":"10.3109\/1040841X.2015.1122571","article-title":"Gene expression variability in clonal populations: causes and consequences","volume":"42","author":"Roberfroid","year":"2016","journal-title":"Crit. Rev. Microbiol"},{"key":"2023020108573601700_btab841-B29","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"edgeR: a Bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2023020108573601700_btab841-B30","doi-asserted-by":"crossref","first-page":"1350","DOI":"10.1126\/science.2999980","article-title":"Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia","volume":"230","author":"Saiki","year":"1985","journal-title":"Science"},{"key":"2023020108573601700_btab841-B31","doi-asserted-by":"crossref","first-page":"770","DOI":"10.1038\/s41588-021-00873-4","article-title":"Separating measurement and expression models clarifies confusion in single-cell RNA sequencing analysis","volume":"53","author":"Sarkar","year":"2021","journal-title":"Nat. Genet"},{"key":"2023020108573601700_btab841-B32","doi-asserted-by":"crossref","first-page":"e9946","DOI":"10.15252\/msb.20209946","article-title":"The transcriptome dynamics of single cells during the cell cycle","volume":"16","author":"Schwabe","year":"2020","journal-title":"Mol. Syst. Biol"},{"key":"2023020108573601700_btab841-B33","doi-asserted-by":"crossref","first-page":"037111","DOI":"10.1063\/1.3184537","article-title":"From puffs to global ca2+ signals: how molecular properties shape global signals","volume":"19","author":"Skupin","year":"2009","journal-title":"Chaos"},{"key":"2023020108573601700_btab841-B34","doi-asserted-by":"crossref","first-page":"12947","DOI":"10.1073\/pnas.93.23.12947","article-title":"Efficiency of DNA replication in the polymerase chain reaction","volume":"93","author":"Stolovitzky","year":"1996","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023020108573601700_btab841-B35","doi-asserted-by":"crossref","first-page":"533","DOI":"10.1093\/nar\/gkz1134","article-title":"Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises","volume":"48","author":"Sun","year":"2020","journal-title":"Nucleic Acids Res"},{"key":"2023020108573601700_btab841-B36","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1016\/j.bdq.2015.01.005","article-title":"How good is a PCR efficiency estimate: recommendations for precise and robust qPCR efficiency assessments","volume":"3","author":"Svec","year":"2015","journal-title":"Biomol. Detect. Quantif"},{"key":"2023020108573601700_btab841-B37","doi-asserted-by":"crossref","first-page":"147","DOI":"10.1038\/s41587-019-0379-5","article-title":"Droplet scRNA-seq is not zero-inflated","volume":"38","author":"Svensson","year":"2020","journal-title":"Nat. Biotechnol"},{"key":"2023020108573601700_btab841-B38","doi-asserted-by":"crossref","first-page":"ra59","DOI":"10.1126\/scisignal.2005237","article-title":"Reliable encoding of stimulus intensities within random sequences of intracellular ca2+ spikes","volume":"7","author":"Thurley","year":"2014","journal-title":"Sci. Signal"},{"key":"2023020108573601700_btab841-B39","doi-asserted-by":"crossref","first-page":"12","DOI":"10.1186\/s13059-019-1850-9","article-title":"A benchmark of batch-effect correction methods for single-cell RNA sequencing data","volume":"21","author":"Tran","year":"2020","journal-title":"Genome Biol"},{"key":"2023020108573601700_btab841-B40","doi-asserted-by":"crossref","first-page":"026601","DOI":"10.1088\/0034-4885\/77\/2\/026601","article-title":"Noise in biology","volume":"77","author":"Tsimring","year":"2014","journal-title":"Rep. Prog. Phys"},{"key":"2023020108573601700_btab841-B41","doi-asserted-by":"crossref","first-page":"39921","DOI":"10.1038\/srep39921","article-title":"Batch effects and the effective design of single-cell gene expression studies","volume":"7","author":"Tung","year":"2017","journal-title":"Sci. Rep"},{"key":"2023020108573601700_btab841-B42","doi-asserted-by":"crossref","first-page":"716","DOI":"10.1016\/j.cell.2018.05.061","article-title":"Recovering gene interactions from single-cell data using data diffusion","volume":"174","author":"van Dijk","year":"2018","journal-title":"Cell"},{"key":"2023020108573601700_btab841-B43","doi-asserted-by":"crossref","first-page":"49","DOI":"10.1089\/cmb.1995.2.49","article-title":"Modeling the polymerase chain reaction","volume":"2","author":"Weiss","year":"1995","journal-title":"J. Comput. Biol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btab841\/41940739\/btab841.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/5\/1336\/49009304\/btab841.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/5\/1336\/49009304\/btab841.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,2,1]],"date-time":"2023-02-01T15:28:44Z","timestamp":1675265324000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/38\/5\/1336\/6462187"}},"subtitle":[],"editor":[{"given":"Peter","family":"Robinson","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2021,12,15]]},"references-count":43,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2022,2,7]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btab841","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2021.10.08.463628","asserted-by":"object"}]},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2022,3,1]]},"published":{"date-parts":[[2021,12,15]]}}}