What Population Reveals about Individual Cell Identity: Single-Cell Parameter Estimation of Models of Gene Expression in Yeast

Significant cell-to-cell heterogeneity is ubiquitously observed in isogenic cell populations. Consequently, parameters of models of intracellular processes, usually fitted to population-averaged data, should rather be fitted to individual cells to obtain a population of models of similar but non-identical individuals. Here, we propose a quantitative modeling framework that attributes specific parameter values to single cells for a standard model of gene expression. We combine high quality single-cell measurements of the response of yeast cells to repeated hyperosmotic shocks and state-of-the-art statistical inference approaches for mixed-effects models to infer multidimensional parameter distributions describing the population, and then derive specific parameters for individual cells. The analysis of single-cell parameters shows that single-cell identity (e.g. gene expression dynamics, cell size, growth rate, mother-daughter relationships) is, at least partially, captured by the parameter values of gene expression models (e.g. rates of transcription, translation and degradation). Our approach shows how to use the rich information contained into longitudinal single-cell data to infer parameters that can faithfully represent single-cell identity.

[1]  James C. W. Locke,et al.  Using movies to analyse gene circuit dynamics in single cells , 2009, Nature Reviews Microbiology.

[2]  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.

[3]  J. Raser,et al.  Control of Stochasticity in Eukaryotic Gene Expression , 2004, Science.

[4]  P. Sorger,et al.  Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis , 2009, Nature.

[5]  J. Raser,et al.  Noise in Gene Expression: Origins, Consequences, and Control , 2005, Science.

[6]  A. Oudenaarden,et al.  Nature, Nurture, or Chance: Stochastic Gene Expression and Its Consequences , 2008, Cell.

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

[8]  Johan Paulsson,et al.  Separating intrinsic from extrinsic fluctuations in dynamic biological systems , 2011, Proceedings of the National Academy of Sciences.

[9]  Megan N. McClean,et al.  Signal processing by the HOG MAP kinase pathway , 2008, Proceedings of the National Academy of Sciences.

[10]  Marc Lavielle,et al.  Maximum likelihood estimation in nonlinear mixed effects models , 2005, Comput. Stat. Data Anal..

[11]  Joaquín Moreno,et al.  Specific and global regulation of mRNA stability during osmotic stress in Saccharomyces cerevisiae. , 2009, RNA.

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

[13]  Mario di Bernardo,et al.  In-Vivo Real-Time Control of Protein Expression from Endogenous and Synthetic Gene Networks , 2014, PLoS Comput. Biol..

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

[15]  C. Pesce,et al.  Regulated cell-to-cell variation in a cell-fate decision system , 2005, Nature.

[16]  Fabian J. Theis,et al.  ODE Constrained Mixture Modelling: A Method for Unraveling Subpopulation Structures and Dynamics , 2014, PLoS Comput. Biol..

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

[18]  É. Moulines,et al.  Convergence of a stochastic approximation version of the EM algorithm , 1999 .

[19]  I. Herskowitz,et al.  Unique and redundant roles for HOG MAPK pathway components as revealed by whole-genome expression analysis. , 2003, Molecular biology of the cell.

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

[21]  Gunnar Cedersund,et al.  Nonlinear mixed-effects modelling for single cell estimation: when, why, and how to use it , 2015, BMC Systems Biology.

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

[23]  António Martins,et al.  A member of the sugar transporter family, Stl1p is the glycerol/H+ symporter in Saccharomyces cerevisiae. , 2005, Molecular biology of the cell.

[24]  Mattias Goksör,et al.  A Nonlinear Mixed Effects Approach for Modeling the Cell-To-Cell Variability of Mig1 Dynamics in Yeast , 2015, PloS one.

[25]  Marc Lavielle,et al.  Mixed Effects Models for the Population Approach: Models, Tasks, Methods and Tools , 2014 .

[26]  David A. Rand,et al.  Measurement of single-cell dynamics , 2010, Nature.

[27]  Sui Huang Non-genetic heterogeneity of cells in development: more than just noise , 2009, Development.

[28]  Dirk Drasdo,et al.  Modeling Dynamics of Cell-to-Cell Variability in TRAIL-Induced Apoptosis Explains Fractional Killing and Predicts Reversible Resistance , 2014, PLoS Comput. Biol..

[29]  P. Swain,et al.  Stochastic Gene Expression in a Single Cell , 2002, Science.

[30]  Eugenio Cinquemani,et al.  Identification of biological models from single-cell data: A comparison between mixed-effects and moment-based inference , 2013, 2013 European Control Conference (ECC).

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

[32]  Mike West,et al.  Bayesian Learning from Marginal Data in Bionetwork Models , 2011, Statistical applications in genetics and molecular biology.

[33]  A. van Oudenaarden,et al.  Noise Propagation in Gene Networks , 2005, Science.

[34]  B. Snijder,et al.  Origins of regulated cell-to-cell variability , 2011, Nature Reviews Molecular Cell Biology.

[35]  Jared E. Toettcher,et al.  Light-based feedback for controlling intracellular signaling dynamics , 2011, Nature Methods.

[36]  Marc Lavielle,et al.  The use of the SAEM algorithm in MONOLIX software for estimation of population pharmacokinetic-pharmacodynamic-viral dynamics parameters of maraviroc in asymptomatic HIV subjects , 2010, Journal of Pharmacokinetics and Pharmacodynamics.