Uncovering the effect of RNA polymerase steric interactions on gene expression noise: analytical distributions of nascent and mature RNA numbers

The telegraph model is the standard model of stochastic gene expression, which can be solved exactly to obtain the distribution of mature RNA numbers per cell. A modification of this model also leads to an analytical distribution of the nascent RNA numbers. These solutions are routinely used for the analysis of single-cell data, including the inference of transcriptional parameters. However, these models neglect important mechanistic features of transcription elongation, such as the stochastic movement of RNA polymerases and their steric interactions. Here we construct a model of gene expression describing promoter switching between inactive and active states, binding of RNA polymerases in the active state, their stochastic movement including steric interactions along the gene, and their unbinding leading to a mature transcript that subsequently decays. We derive the steady-state distributions of the nascent and mature RNA numbers in two important limiting cases: constitutive expression and slow promoter switching. We show that RNA fluctuations are suppressed by steric interactions between RNA polymerases, and that this suppression can even lead to sub-Poissonian fluctuations; these effects are most pronounced for nascent RNA and less prominent for mature RNA, since the latter is not a direct sensor of transcription. We find a relationship between the parameters of our microscopic mechanistic model and those of the standard models that ensures excellent consistency in their prediction of the first and second RNA number moments over vast regions of parameter space, encompassing slow, intermediate, and rapid promoter switching, provided the RNA number distributions are Poissonian or super-Poissonian. Furthermore, we identify the limitations of inference from mature RNA data, specifically showing that it cannot differentiate between highly distinct RNA polymerase traffic patterns on a gene.

[1]  Abhyudai Singh,et al.  The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian , 2023, bioRxiv.

[2]  T. Lenstra,et al.  Quantifying how post-transcriptional noise and gene copy number variation bias transcriptional parameter inference from mRNA distributions , 2022, bioRxiv.

[3]  J. Szavits-Nossan,et al.  Steady-state distributions of nascent RNA for general initiation mechanisms , 2022, bioRxiv.

[4]  J. Szavits-Nossan,et al.  Mean-field theory accurately captures the variation of copy number distributions across the mRNA life cycle. , 2022, Physical review. E.

[5]  R. Grima,et al.  Distinguishing between models of mammalian gene expression: telegraph-like models versus mechanistic models , 2021, bioRxiv.

[6]  W. Du,et al.  Neural network aided approximation and parameter inference of non-Markovian models of gene expression , 2021, Nature Communications.

[7]  R. Grima,et al.  Frequency domain analysis of fluctuations of mRNA and protein copy numbers within a cell lineage: theory and experimental validation , 2020, bioRxiv.

[8]  M. Khammash,et al.  Frequency spectra and the color of cellular noise , 2020, Nature Communications.

[9]  B. Waclaw,et al.  Current-density relation in the exclusion process with dynamic obstacles. , 2020, Physical review. E.

[10]  Tatiana Filatova,et al.  Statistics of Nascent and Mature RNA Fluctuations in a Stochastic Model of Transcriptional Initiation, Elongation, Pausing, and Termination , 2020, Bulletin of Mathematical Biology.

[11]  E. Petfalski,et al.  Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates , 2020, bioRxiv.

[12]  J. Szavits-Nossan,et al.  Dynamics of ribosomes in mRNA translation under steady- and nonsteady-state conditions. , 2020, Physical review. E.

[13]  Ramon Grima,et al.  Analytical distributions for detailed models of stochastic gene expression in eukaryotic cells , 2020, Proceedings of the National Academy of Sciences.

[14]  P. Hrabák Time-headway distribution for random-sequential-update TASEP with periodic and open boundaries , 2020 .

[15]  S. Karthika,et al.  Totally asymmetric simple exclusion process with resetting , 2020, Journal of Physics A: Mathematical and Theoretical.

[16]  R. Grima,et al.  Small protein number effects in stochastic models of autoregulated bursty gene expression. , 2020, The Journal of chemical physics.

[17]  T. Lenstra,et al.  Live‐cell imaging reveals the interplay between transcription factors, nucleosomes, and bursting , 2019, The EMBO journal.

[18]  R. Brewster,et al.  Probing mechanisms of transcription elongation through cell-to-cell variability of RNA polymerase , 2019, bioRxiv.

[19]  Lucy Ham,et al.  Extrinsic noise and heavy-tailed laws in gene expression , 2019, bioRxiv.

[20]  R. Sandberg,et al.  Genomic encoding of transcriptional burst kinetics , 2019, Nature.

[21]  B. Waclaw,et al.  Quantitative modelling predicts the impact of DNA methylation on RNA polymerase II traffic , 2019, Proceedings of the National Academy of Sciences.

[22]  Sandeep Choubey Nascent RNA kinetics: Transient and steady state behavior of models of transcription. , 2018, Physical review. E.

[23]  S. Legewie,et al.  Estrogen‐dependent control and cell‐to‐cell variability of transcriptional bursting , 2018, Molecular systems biology.

[24]  Julia Zeitlinger,et al.  Paused RNA polymerase II inhibits new transcriptional initiation , 2017, Nature Genetics.

[25]  Lorenza Vitale,et al.  GeneBase 1.1: a tool to summarize data from NCBI gene datasets and its application to an update of human gene statistics , 2016, Database J. Biol. Databases Curation.

[26]  Ido Golding,et al.  Stochastic Kinetics of Nascent RNA. , 2016, Physical review letters.

[27]  Charles G. Morgan,et al.  A Mechanistic Model for Cooperative Behavior of Co-transcribing RNA Polymerases , 2016, PLoS Comput. Biol..

[28]  Samuel M. D. Oliveira,et al.  Dissecting the stochastic transcription initiation process in live Escherichia coli , 2016, DNA research : an international journal for rapid publication of reports on genes and genomes.

[29]  R. Grima,et al.  Molecular finite-size effects in stochastic models of equilibrium chemical systems. , 2015, The Journal of chemical physics.

[30]  Yu Rim Lim,et al.  Quantitative Understanding of Probabilistic Behavior of Living Cells Operated by Vibrant Intracellular Networks , 2015 .

[31]  A. Raj,et al.  Single mammalian cells compensate for differences in cellular volume and DNA copy number through independent global transcriptional mechanisms. , 2015, Molecular cell.

[32]  S. Itzkovitz,et al.  Bursty gene expression in the intact mammalian liver. , 2015, Molecular cell.

[33]  Ryan A. Kellogg,et al.  Noise Facilitates Transcriptional Control under Dynamic Inputs , 2015, Cell.

[34]  Niraj Kumar,et al.  Exact distributions for stochastic gene expression models with bursting and feedback. , 2014, Physical review letters.

[35]  Brian Munsky,et al.  Transcription Factors Modulate c-Fos Transcriptional Bursts , 2014, Cell reports.

[36]  J. Lis,et al.  Genome-wide dynamics of Pol II elongation and its interplay with promoter proximal pausing, chromatin, and exons , 2014, eLife.

[37]  N. Popović,et al.  Phenotypic switching in gene regulatory networks , 2014, Proceedings of the National Academy of Sciences.

[38]  Sandeep Choubey,et al.  Deciphering Transcriptional Dynamics In Vivo by Counting Nascent RNA Molecules , 2013, PLoS Comput. Biol..

[39]  S. Klumpp,et al.  Backtracking dynamics of RNA polymerase: pausing and error correction , 2013, Journal of physics. Condensed matter : an Institute of Physics journal.

[40]  J. Marioni,et al.  Inferring the kinetics of stochastic gene expression from single-cell RNA-sequencing data , 2013, Genome Biology.

[41]  D. Chowdhury Stochastic mechano-chemical kinetics of molecular motors: A multidisciplinary enterprise from a physicist’s perspective , 2012, 1207.6070.

[42]  Tianshou Zhou,et al.  Analytical Results for a Multistate Gene Model , 2012, SIAM J. Appl. Math..

[43]  R. Lipowsky,et al.  Translation by Ribosomes with mRNA Degradation: Exclusion Processes on Aging Tracks , 2011 .

[44]  Nacho Molina,et al.  Mammalian Genes Are Transcribed with Widely Different Bursting Kinetics , 2011, Science.

[45]  R. Segev,et al.  GENERAL PROPERTIES OF THE TRANSCRIPTIONAL TIME-SERIES IN ESCHERICHIA COLI , 2011, Nature Genetics.

[46]  A. Oudenaarden,et al.  Cellular Decision Making and Biological Noise: From Microbes to Mammals , 2011, Cell.

[47]  Stefan Klumpp,et al.  Pausing and Backtracking in Transcription Under Dense Traffic Conditions , 2011 .

[48]  P. Hrabák,et al.  Inter-particle gap distribution and spectral rigidity of the totally asymmetric simple exclusion process with open boundaries , 2010, 1011.0196.

[49]  Tao Jia,et al.  Applications of Little's Law to stochastic models of gene expression. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[50]  J. Hespanha,et al.  Optimal feedback strength for noise suppression in autoregulatory gene networks. , 2009, Biophysical journal.

[51]  G. Schütz,et al.  RNA polymerase motors: dwell time distribution, velocity and dynamical phases , 2009, 0904.2625.

[52]  F. Hayot,et al.  Stochasticity of gene products from transcriptional pulsing. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[53]  F. Bruggeman,et al.  Elongation dynamics shape bursty transcription and translation , 2009, Proceedings of the National Academy of Sciences.

[54]  Terence Hwa,et al.  Stochasticity and traffic jams in the transcription of ribosomal RNA: Intriguing role of termination and antitermination , 2008, Proceedings of the National Academy of Sciences.

[55]  D. Larson,et al.  Single-RNA counting reveals alternative modes of gene expression in yeast , 2008, Nature Structural &Molecular Biology.

[56]  Vahid Shahrezaei,et al.  Analytical distributions for stochastic gene expression , 2008, Proceedings of the National Academy of Sciences.

[57]  D. Chowdhury,et al.  Transcriptional bursts: A unified model of machines and mechanisms , 2008, 0804.1227.

[58]  N. Cohen,et al.  Fluctuations, pauses, and backtracking in DNA transcription. , 2008, Biophysical journal.

[59]  Debashish Chowdhury,et al.  Interacting RNA polymerase motors on a DNA track: effects of traffic congestion and intrinsic noise on RNA synthesis. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[60]  Nir Friedman,et al.  Linking stochastic dynamics to population distribution: an analytical framework of gene expression. , 2006, Physical review letters.

[61]  J R Yates,et al.  RNA polymerase II elongation factors Spt4p and Spt5p play roles in transcription elongation by RNA polymerase I and rRNA processing , 2006, Proceedings of the National Academy of Sciences.

[62]  F. Essler,et al.  Bethe ansatz solution of the asymmetric exclusion process with open boundaries. , 2005, Physical review letters.

[63]  M. Nomura,et al.  Histones are required for transcription of yeast rRNA genes by RNA polymerase I. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[64]  Johan Paulsson,et al.  Models of stochastic gene expression , 2005 .

[65]  M. Nomura,et al.  Tor pathway regulates Rrn3p-dependent recruitment of yeast RNA polymerase I to the promoter but does not participate in alteration of the number of active genes. , 2003, Molecular biology of the cell.

[66]  C. Rao,et al.  Control, exploitation and tolerance of intracellular noise , 2002, Nature.

[67]  Arnab Majumdar,et al.  Distribution of time-headways in a particle-hopping model of vehicular traffic , 1998 .

[68]  J. Peccoud,et al.  Markovian Modeling of Gene-Product Synthesis , 1995 .

[69]  B. Derrida,et al.  Exact solution of a 1d asymmetric exclusion model using a matrix formulation , 1993 .

[70]  E. Domany,et al.  Phase transitions in an exactly soluble one-dimensional exclusion process , 1993, cond-mat/9303038.

[71]  Bernard Derrida,et al.  Exact correlation functions in an asymmetric exclusion model with open boundaries , 1993 .

[72]  Eytan Domany,et al.  An exact solution of a one-dimensional asymmetric exclusion model with open boundaries , 1992 .

[73]  O. Miller,et al.  Transcription mapping of the Escherichia coli chromosome by electron microscopy , 1989, Journal of bacteriology.

[74]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[75]  A. Pipkin,et al.  Kinetics of biopolymerization on nucleic acid templates , 1968, Biopolymers.

[76]  L. Takács,et al.  On a coincidence problem concerning telephone traffic , 1958 .

[77]  D. Kendall Stochastic Processes Occurring in the Theory of Queues and their Analysis by the Method of the Imbedded Markov Chain , 1953 .

[78]  John D. C. Little,et al.  A PROOF FOR THE QUEUING FORMULA : , 2015 .

[79]  M. R. Evanst,et al.  A : Mathematical and General Exact solution of a 1 D asymmetric exclusion model using a matrix formulation , 2002 .

[80]  J D Littler,et al.  A PROOF OF THE QUEUING FORMULA , 1961 .