Kinetic foundation of the zero-inflated negative binomial model for single-cell RNA sequencing data

Abstract

Single-cell RNA sequencing data have complex features such as dropout events, over-dispersion, and high-magnitude outliers, resulting in complicated probability distributions of mRNA abundances that are statistically characterized in terms of a zero-inflated negative binomial (ZINB) model. Here we provide a mesoscopic kinetic foundation of the widely used ZINB model based on the biochemical reaction kinetics underlying transcription. Using multiscale modeling and simplification techniques, we show that the ZINB distribution of mRNA abundance and the phenomenon of transcriptional bursting naturally emerge from a three-state stochastic transcription model. We further reveal a nontrivial quantitative relation between dropout events and transcriptional bursting, which provides novel insights into how and to what extent the burst size and burst frequency could reduce the dropout rate. Three different biophysical origins of over-dispersion are also clarified at the single-cell level.

19 pages, 5 figures

References in corpus (5)

• Multiple-scale stochastic processes: decimation, averaging and beyond
• Spectral solutions to stochastic models of gene expression with bursts and regulation
• Simplification of Markov chains with infinite state space and the mathematical theory of random gene expression bursts
• Coarse graining of master equations with fast and slow states
• Emergent Lévy behavior in single-cell stochastic gene expression