Reconstructing Statistics of Promoter Switching from Reporter Protein Population Snapshot Data

The use of fluorescent reporter proteins is an established experimental approach for dynamic quantification of gene expression over time. Yet, the observed fluorescence levels are only indirect measurements of the relevant promoter activity. At the level of population averages, reconstruction of mean activity profiles from mean fluorescence profiles has been addressed with satisfactory results. At the single cell level, however, promoter activity is generally different from cell to cell. Making sense of this variability is at the core of single-cell modelling, but complicates the reconstruction task. Here we discuss reconstruction of promoter activity statistics from time-lapse population snapshots of fluorescent reporter statistics, as obtained e.g. by flow-cytometric measurements of a dynamical gene expression experiment. After discussing the problem in the framework of stochastic modelling, we provide an estimation method based on convex optimization. We then instantiate it in the fundamental case of a single promoter switch, reflecting a typical random promoter activation or deactivation, and discuss estimation results from in silico experiments.

[1]  Peter S. Swain,et al.  The Fidelity of Dynamic Signaling by Noisy Biomolecular Networks , 2013, PLoS Comput. Biol..

[2]  John Lygeros,et al.  Grey-box techniques for the identification of a controlled gene expression model , 2014, 2014 European Control Conference (ECC).

[3]  Sheng Wu,et al.  StochKit2: software for discrete stochastic simulation of biochemical systems with events , 2011, Bioinform..

[4]  John Lygeros,et al.  Local Identification of Piecewise Deterministic Models of Genetic Networks , 2009, HSCC.

[5]  Gabriele Lillacci,et al.  The signal within the noise: efficient inference of stochastic gene regulation models using fluorescence histograms and stochastic simulations , 2013, Bioinform..

[6]  Hidde de Jong,et al.  Experimental and computational validation of models of fluorescent and luminescent reporter genes in bacteria , 2010, BMC Systems Biology.

[7]  G. Picci,et al.  Linear Stochastic Systems: A Geometric Approach to Modeling, Estimation and Identification , 2016 .

[8]  D. Gillespie.  The chemical Langevin equation , 2000 .

[9]  Frank Allgöwer,et al.  Identification of models of heterogeneous cell populations from population snapshot data , 2011, BMC Bioinformatics.

[10]  John Lygeros,et al.  Moment-Based Methods for Parameter Inference and Experiment Design for Stochastic Biochemical Reaction Networks , 2015, ACM Trans. Model. Comput. Simul..

[11]  J. Hespanha.  Modelling and analysis of stochastic hybrid systems , 2006 .

[12]  Fabian J. Theis,et al.  Reconstructing gene regulatory dynamics from high-dimensional single-cell snapshot data , 2015, Bioinform..

[13]  F. Fages,et al.  Long-term model predictive control of gene expression at the population and single-cell levels , 2012, Proceedings of the National Academy of Sciences.

[14]  G. Wahba.  Spline models for observational data , 1990 .

[15]  John Lygeros,et al.  Identification of genetic network dynamics with unate structure , 2010, Bioinform..

[16]  Brian Munsky,et al.  Listening to the noise: random fluctuations reveal gene network parameters , 2009, Molecular systems biology.

[17]  Paul J. Choi,et al.  Quantifying E. coli Proteome and Transcriptome with Single-Molecule Sensitivity in Single Cells , 2010, Science.

[18]  Andrew J. Millar,et al.  Reconstruction of transcriptional dynamics from gene reporter data using differential equations , 2008, Bioinform..

[19]  T. Elston,et al.  Stochasticity in gene expression: from theories to phenotypes , 2005, Nature Reviews Genetics.

[20]  J. Lygeros,et al.  Moment-based inference predicts bimodality in transient gene expression , 2012, Proceedings of the National Academy of Sciences.

[21]  M. Khammash,et al.  Systematic Identification of Signal-Activated Stochastic Gene Regulation , 2013, Science.

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

[23]  M. Thattai,et al.  Intrinsic noise in gene regulatory networks , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[24]  David A. Rand,et al.  Bayesian inference of biochemical kinetic parameters using the linear noise approximation , 2009, BMC Bioinformatics.

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

[26]  Valentin Zulkower,et al.  Robust reconstruction of gene expression profiles from reporter gene data using linear inversion , 2015, Bioinform..

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

[28]  D. Pincus,et al.  In silico feedback for in vivo regulation of a gene expression circuit , 2011, Nature Biotechnology.

[29]  Eugenio Cinquemani,et al.  Reconstruction of promoter activity statistics from reporter protein population snapshot data , 2015, 2015 54th IEEE Conference on Decision and Control (CDC).

[30]  J. Lygeros,et al.  Moment estimation for chemically reacting systems by extended Kalman filtering. , 2011, The Journal of chemical physics.

[31]  Linda R. Petzold,et al.  Stochastic modelling of gene regulatory networks , 2005 .

[32]  Eugenio Cinquemani,et al.  Inference of Quantitative Models of Bacterial Promoters from Time-Series Reporter Gene Data , 2015, PLoS Comput. Biol..

[33]  M. Peter,et al.  Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings , 2013, Nature Methods.

[34]  D. di Bernardo,et al.  How to infer gene networks from expression profiles , 2007, Molecular systems biology.