{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,28]],"date-time":"2026-02-28T10:08:01Z","timestamp":1772273281485,"version":"3.50.1"},"reference-count":42,"publisher":"Oxford University Press (OUP)","issue":"17","license":[{"start":{"date-parts":[[2018,9,9]],"date-time":"2018-09-09T00:00:00Z","timestamp":1536451200000},"content-version":"vor","delay-in-days":8,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100000266","name":"EPSRC","doi-asserted-by":"publisher","award":["1265509"],"award-info":[{"award-number":["1265509"]}],"id":[{"id":"10.13039\/501100000266","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000268","name":"BBSRC","doi-asserted-by":"publisher","award":["BB\/L006340\/1"],"award-info":[{"award-number":["BB\/L006340\/1"]}],"id":[{"id":"10.13039\/501100000268","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000268","name":"BBSRC","doi-asserted-by":"publisher","award":["BB\/M017982\/1"],"award-info":[{"award-number":["BB\/M017982\/1"]}],"id":[{"id":"10.13039\/501100000268","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100000265","name":"Medical Research Council","doi-asserted-by":"publisher","award":["MR\/M013170\/1"],"award-info":[{"award-number":["MR\/M013170\/1"]}],"id":[{"id":"10.13039\/501100000265","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2018,9,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Transcription in single cells is an inherently stochastic process as mRNA levels vary greatly between cells, even for genetically identical cells under the same experimental and environmental conditions. We present a stochastic two-state switch model for the population of mRNA molecules in single cells where genes stochastically alternate between a more active ON state and a less active OFF state. We prove that the stationary solution of such a model can be written as a mixture of a Poisson and a Poisson-beta probability distribution. This finding facilitates inference for single cell expression data, observed at a single time point, from flow cytometry experiments such as FACS or fluorescence in situ hybridization (FISH) as it allows one to sample directly from the equilibrium distribution of the mRNA population. We hence propose a Bayesian inferential methodology using a pseudo-marginal approach and a recent approximation to integrate over unobserved states associated with measurement error.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>We provide a general inferential framework which can be widely used to study transcription in single cells from the kind of data arising in flow cytometry experiments. The approach allows us to separate between the intrinsic stochasticity of the molecular dynamics and the measurement noise. The methodology is tested in simulation studies and results are obtained for experimental multiple single cell expression data from FISH flow cytometry experiments.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>All analyses were implemented in R. Source code and the experimental data are available at https:\/\/github.com\/SimoneTiberi\/Bayesian-inference-on-stochastic-gene-transcription-from-flow-cytometry-data.<\/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\/bty568","type":"journal-article","created":{"date-parts":[[2018,7,6]],"date-time":"2018-07-06T01:20:56Z","timestamp":1530840056000},"page":"i647-i655","source":"Crossref","is-referenced-by-count":30,"title":["Bayesian inference on stochastic gene transcription from flow cytometry data"],"prefix":"10.1093","volume":"34","author":[{"given":"Simone","family":"Tiberi","sequence":"first","affiliation":[{"name":"Institute of Molecular Life Sciences, University of Z\u00fcrich, Z\u00fcrich, Switzerland"},{"name":"Swiss Institue of Bioinformatics, University of Z\u00fcrich, Z\u00fcrich, Switzerland"},{"name":"Department of Statistics, University of Warwick, Coventry, UK"}]},{"given":"Mark","family":"Walsh","sequence":"additional","affiliation":[{"name":"School of Life Sciences, University of Warwick, Coventry, UK"}]},{"given":"Massimo","family":"Cavallaro","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of Warwick, Coventry, UK"},{"name":"School of Life Sciences, University of Warwick, Coventry, UK"}]},{"given":"Daniel","family":"Hebenstreit","sequence":"additional","affiliation":[{"name":"School of Life Sciences, University of Warwick, Coventry, UK"}]},{"given":"B\u00e4rbel","family":"Finkenst\u00e4dt","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of Warwick, Coventry, UK"}]}],"member":"286","published-online":{"date-parts":[[2018,9,8]]},"reference":[{"key":"2023061313512807200_bty568-B1","doi-asserted-by":"crossref","first-page":"697","DOI":"10.1214\/07-AOS574","article-title":"The pseudo-marginal approach for efficient Monte Carlo computations","volume":"37","author":"Andrieu","year":"2009","journal-title":"Ann. Stat"},{"key":"2023061313512807200_bty568-B2","doi-asserted-by":"crossref","first-page":"1139","DOI":"10.1093\/genetics\/164.3.1139","article-title":"Estimation of population growth or decline in genetically monitored populations","volume":"164","author":"Beaumont","year":"2003","journal-title":"Genetics"},{"key":"2023061313512807200_bty568-B3","doi-asserted-by":"crossref","first-page":"271","DOI":"10.1016\/j.molcel.2007.11.035","article-title":"A 5\u2019 splice site enhances the recruitment of basal transcription initiation factors in vivo","volume":"29","author":"Damgaard","year":"2008","journal-title":"Mol. Cell"},{"key":"2023061313512807200_bty568-B4","doi-asserted-by":"crossref","first-page":"17454","DOI":"10.1073\/pnas.1213530109","article-title":"Transcriptional burst frequency and burst size are equally modulated across the human genome","volume":"109","author":"Dar","year":"2012","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023061313512807200_bty568-B5","doi-asserted-by":"crossref","first-page":"20160833.","DOI":"10.1098\/rsif.2016.0833","article-title":"Stochastic models of gene transcription with upstream drives: exact solution and sample path characterization","volume":"14","author":"Dattani","year":"2017","journal-title":"J. R. Soc. Interface"},{"key":"2023061313512807200_bty568-B6","doi-asserted-by":"crossref","first-page":"120","DOI":"10.1063\/1.1750549","article-title":"Statistical fluctuations in autocatalytic reactions","volume":"8","author":"Delbr\u00fcck","year":"1940","journal-title":"J. Chem. Phys"},{"key":"2023061313512807200_bty568-B7","doi-asserted-by":"crossref","first-page":"646","DOI":"10.1016\/j.cels.2017.10.013","article-title":"Asymmetry between activation and deactivation during a transcriptional pulse","volume":"5","author":"Dunham","year":"2017","journal-title":"Cell Syst"},{"key":"2023061313512807200_bty568-B8","doi-asserted-by":"crossref","first-page":"e08494","DOI":"10.7554\/eLife.08494","article-title":"Spatially coordinated dynamic gene transcription in living pituitary tissue","volume":"5","author":"Featherstone","year":"2016","journal-title":"Elife"},{"key":"2023061313512807200_bty568-B9","doi-asserted-by":"crossref","first-page":"1960","DOI":"10.1214\/13-AOAS669","article-title":"Quantifying intrinsic and extrinsic noise in gene transcription using the linear noise approximation: an application to single cell data","volume":"7","author":"Finkenst\u00e4dt","year":"2013","journal-title":"Ann. Appl. Stat"},{"key":"2023061313512807200_bty568-B10","author":"Flegal","year":"2017"},{"key":"2023061313512807200_bty568-B11","doi-asserted-by":"crossref","DOI":"10.1201\/9781482296426","volume-title":"Markov Chain Monte Carlo Stochastic Simulation for Bayesian Inference","author":"Gamerman","year":"2006","edition":"2"},{"key":"2023061313512807200_bty568-B12","doi-asserted-by":"crossref","first-page":"1025","DOI":"10.1016\/j.cell.2005.09.031","article-title":"Real-time kinetics of gene activity in individual bacteria","volume":"123","author":"Golding","year":"2005","journal-title":"Cell"},{"key":"2023061313512807200_bty568-B13","doi-asserted-by":"crossref","first-page":"223","DOI":"10.2307\/3318737","article-title":"An adaptive metropolis algorithm","volume":"7","author":"Haario","year":"2001","journal-title":"Bernoulli"},{"key":"2023061313512807200_bty568-B14","doi-asserted-by":"crossref","first-page":"e1000607.","DOI":"10.1371\/journal.pbio.1000607","article-title":"Dynamic analysis of stochastic transcription cycles","volume":"9","author":"Harper","year":"2011","journal-title":"PLoS Biol"},{"key":"2023061313512807200_bty568-B15","doi-asserted-by":"crossref","first-page":"97","DOI":"10.1093\/biomet\/57.1.97","article-title":"Monte Carlo sampling methods using Markov chains and their applications","volume":"57","author":"Hastings","year":"1970","journal-title":"Biometrika"},{"key":"2023061313512807200_bty568-B16","doi-asserted-by":"crossref","first-page":"333","DOI":"10.1016\/j.tig.2013.04.001","article-title":"Are gene loops the cause of transcriptional noise?","volume":"29","author":"Hebenstreit","year":"2013","journal-title":"Trends Genet"},{"key":"2023061313512807200_bty568-B17","doi-asserted-by":"crossref","first-page":"497","DOI":"10.1038\/msb.2011.28","article-title":"Rna sequencing reveals two major classes of gene expression levels in metazoan cells","volume":"7","author":"Hebenstreit","year":"2011","journal-title":"Mol. Syst. Biol"},{"key":"2023061313512807200_bty568-B18","doi-asserted-by":"crossref","first-page":"233","DOI":"10.1145\/358598.358630","article-title":"A spectral method for confidence interval generation and run length control in simulations","volume":"24","author":"Heidelberger","year":"1981","journal-title":"Commun. ACM"},{"key":"2023061313512807200_bty568-B19","doi-asserted-by":"crossref","first-page":"1109","DOI":"10.1287\/opre.31.6.1109","article-title":"Simulation run length control in the presence of an initial transient","volume":"31","author":"Heidelberger","year":"1983","journal-title":"Oper. Res"},{"key":"2023061313512807200_bty568-B20","doi-asserted-by":"crossref","first-page":"655","DOI":"10.1093\/biostatistics\/kxv010","article-title":"A stochastic transcriptional switch model for single cell imaging data","volume":"16","author":"Hey","year":"2015","journal-title":"Biostatistics"},{"key":"2023061313512807200_bty568-B21","doi-asserted-by":"crossref","DOI":"10.1002\/0471715816","volume-title":"Univariate Discrete Distributions","author":"Johnson","year":"2005"},{"key":"2023061313512807200_bty568-B22","doi-asserted-by":"crossref","first-page":"451","DOI":"10.1038\/nrg1615","article-title":"Stochasticity in gene expression: from theories to phenotypes","volume":"6","author":"Kaern","year":"2005","journal-title":"Nat. Rev. Genet"},{"key":"2023061313512807200_bty568-B23","doi-asserted-by":"crossref","first-page":"3116","DOI":"10.1016\/S0006-3495(01)75949-8","article-title":"Stochasticity in transcriptional regulation: origins, consequences, and mathematical representations","volume":"81","author":"Kepler","year":"2001","journal-title":"Biophys. J"},{"key":"2023061313512807200_bty568-B24","doi-asserted-by":"crossref","first-page":"R7.","DOI":"10.1186\/gb-2013-14-1-r7","article-title":"Inferring the kinetics of stochastic gene expression from single-cell RNA-sequencing data","volume":"14","author":"Kim","year":"2013","journal-title":"Genome Biol"},{"key":"2023061313512807200_bty568-B25","doi-asserted-by":"crossref","first-page":"1087","DOI":"10.1063\/1.1699114","article-title":"Equation of state calculations by fast computing machines","volume":"21","author":"Metropolis","year":"1953","journal-title":"J. Chem. Phys"},{"key":"2023061313512807200_bty568-B26","doi-asserted-by":"crossref","first-page":"335","DOI":"10.1080\/01621459.1949.10483310","article-title":"The Monte Carlo method","volume":"44","author":"Metropolis","year":"1949","journal-title":"J. Am. Stat. Assoc"},{"key":"2023061313512807200_bty568-B27","doi-asserted-by":"crossref","first-page":"20563","DOI":"10.1073\/pnas.1312310110","article-title":"Stimulus-induced modulation of transcriptional bursting in a single mammalian gene","volume":"110","author":"Molina","year":"2013","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023061313512807200_bty568-B28","doi-asserted-by":"crossref","first-page":"179","DOI":"10.1007\/s002110100285","article-title":"Computing the confluent hypergeometric function, M (a, b, x)","volume":"90","author":"Muller","year":"2001","journal-title":"Numerische Mathematik"},{"key":"2023061313512807200_bty568-B29","doi-asserted-by":"crossref","first-page":"183","DOI":"10.1126\/science.1216379","article-title":"Using gene expression noise to understand gene regulation","volume":"336","author":"Munsky","year":"2012","journal-title":"Science"},{"key":"2023061313512807200_bty568-B30","doi-asserted-by":"crossref","first-page":"157","DOI":"10.1016\/j.plrev.2005.03.003","article-title":"Models of stochastic gene expression","volume":"2","author":"Paulsson","year":"2005","journal-title":"Phys. Life Rev"},{"key":"2023061313512807200_bty568-B31","doi-asserted-by":"crossref","first-page":"222","DOI":"10.1006\/tpbi.1995.1027","article-title":"Markovian modeling of gene-product synthesis","volume":"48","author":"Peccoud","year":"1995","journal-title":"Theor. Popul. Biol"},{"key":"2023061313512807200_bty568-B32","author":"Plummer","year":"2016"},{"key":"2023061313512807200_bty568-B33","volume-title":"R: A Language and Environment for Statistical Computing","author":"R Core Team","year":"2016"},{"key":"2023061313512807200_bty568-B34","doi-asserted-by":"crossref","first-page":"e309.","DOI":"10.1371\/journal.pbio.0040309","article-title":"Stochastic mRNA synthesis in mammalian cells","volume":"4","author":"Raj","year":"2006","journal-title":"PLoS Biol"},{"key":"2023061313512807200_bty568-B35","doi-asserted-by":"crossref","first-page":"e1000704","DOI":"10.1371\/journal.pcbi.1000704","article-title":"Effects of transcriptional pausing on gene expression dynamics","volume":"6","author":"Rajala","year":"2010","journal-title":"PLoS Comput. Biol"},{"key":"2023061313512807200_bty568-B36","doi-asserted-by":"crossref","first-page":"17256","DOI":"10.1073\/pnas.0803850105","article-title":"Analytical distributions for stochastic gene expression","volume":"105","author":"Shahrezaei","year":"2008","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023061313512807200_bty568-B37","author":"Singh","year":"2013"},{"key":"2023061313512807200_bty568-B38","doi-asserted-by":"crossref","first-page":"554","DOI":"10.1038\/ng.821","article-title":"General properties of transcriptional time series in Escherichia coli","volume":"43","author":"So","year":"2011","journal-title":"Nat. Genet"},{"key":"2023061313512807200_bty568-B39","doi-asserted-by":"crossref","first-page":"472","DOI":"10.1126\/science.1198817","article-title":"Mammalian genes are transcribed with widely different bursting kinetics","volume":"332","author":"Suter","year":"2011","journal-title":"Science"},{"key":"2023061313512807200_bty568-B40","doi-asserted-by":"crossref","first-page":"6994","DOI":"10.1073\/pnas.1400049111","article-title":"Phenotypic switching in gene regulatory networks","volume":"111","author":"Thomas","year":"2014","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023061313512807200_bty568-B41","doi-asserted-by":"crossref","first-page":"748","DOI":"10.1038\/nbt.2642","article-title":"Single-cell gene expression analysis reveals genetic associations masked in whole-tissue experiments","volume":"31","author":"Wills","year":"2013","journal-title":"Nat. Biotechnol"},{"key":"2023061313512807200_bty568-B42","doi-asserted-by":"crossref","first-page":"e1003161","DOI":"10.1371\/journal.pcbi.1003161","article-title":"Cell-cycle dependence of transcription dominates noise in gene expression","volume":"9","author":"Zopf","year":"2013","journal-title":"PLoS Comput. Biol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/17\/i647\/50582556\/bioinformatics_34_17_i647.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/17\/i647\/50582556\/bioinformatics_34_17_i647.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,6,13]],"date-time":"2023-06-13T13:53:41Z","timestamp":1686664421000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/34\/17\/i647\/5093237"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2018,9,1]]},"references-count":42,"journal-issue":{"issue":"17","published-print":{"date-parts":[[2018,9,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty568","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2018,9,1]]},"published":{"date-parts":[[2018,9,1]]}}}