Cell-to-cell variability in JAK2/STAT5 pathway components and cytoplasmic volumes define survival threshold in erythroid progenitor cells

Survival or apoptosis is a binary decision in individual cells. Yet, at the cell population level, a graded increase in survival of CFU-E cells is observed upon stimulation with Erythropoietin (Epo). To identify components of JAK2/STAT5 signal transduction that contribute to the graded population response, a cell population-level model calibrated with experimental data was extended to study the behavior in single cells. The single-cell model showed that the high cell-to-cell variability in nuclear phosphorylated STAT5 is caused by variability in the amount of EpoR:JAK2 complexes and of SHP1 as well as the extent of nuclear import due to the large variance in the cytoplasmic volume of CFU-E cells. 24 to 118 pSTAT5 molecules in the nucleus for 120 min are sufficient to ensure cell survival. Thus, variability in membrane-associated processes are responsible to convert a switch-like behavior at the single-cell level to a graded population level response. Highlights Mathematical modeling enables integration of heterogeneous data Single-cell modeling captures binary decision process Multiple sources of cell-to-cell variability in erythroid progenitor cells Minimal amount of active STAT5 sufficient for survival of erythroid progenitor cells

[1]  R. Milo,et al.  Assessment of GFP Tag Position on Protein Localization and Growth Fitness in Yeast. , 2019, Journal of molecular biology.

[2]  Thomas S. Ligon,et al.  Multi-experiment nonlinear mixed effect modeling of single-cell translation kinetics after transfection , 2018, npj Systems Biology and Applications.

[3]  Jan Hasenauer,et al.  A Hierarchical, Data-Driven Approach to Modeling Single-Cell Populations Predicts Latent Causes of Cell-To-Cell Variability. , 2018, Cell systems.

[4]  Jens Timmer,et al.  Profile likelihood-based analyses of infectious disease models , 2018, Statistical methods in medical research.

[5]  Jan Hasenauer,et al.  PESTO: Parameter EStimation TOolbox , 2017, Bioinform..

[6]  G. Borthakur,et al.  Erythroleukemia-historical perspectives and recent advances in diagnosis and management. , 2017, Blood reviews.

[7]  Jens Timmer,et al.  Mathematical model of early Reelin-induced Src family kinase-mediated signaling , 2017, PloS one.

[8]  L. Goentoro,et al.  Sensing relative signal in the Tgf-β/Smad pathway , 2017, Proceedings of the National Academy of Sciences.

[9]  Melanie Boerries,et al.  Protein abundance of AKT and ERK pathway components governs cell type‐specific regulation of proliferation , 2017, Molecular systems biology.

[10]  Fabian J. Theis,et al.  Parameter estimation for dynamical systems with discrete events and logical operations , 2016, Bioinform..

[11]  Koichi Takahashi,et al.  Modeling Cellular Noise Underlying Heterogeneous Cell Responses in the Epidermal Growth Factor Signaling Pathway , 2016, PLoS Comput. Biol..

[12]  M. Selbach,et al.  Kinetic Analysis of Protein Stability Reveals Age-Dependent Degradation , 2016, Cell.

[13]  Jens Timmer,et al.  Driving the Model to Its Limit: Profile Likelihood Based Model Reduction , 2016, PloS one.

[14]  Ursula Klingmüller,et al.  Identification of Cell Type-Specific Differences in Erythropoietin Receptor Signaling in Primary Erythroid and Lung Cancer Cells , 2016, PLoS Comput. Biol..

[15]  L. Pelkmans,et al.  Passive Noise Filtering by Cellular Compartmentalization , 2016, Cell.

[16]  Fabian J Theis,et al.  Diffusion pseudotime robustly reconstructs lineage branching , 2016, Nature Methods.

[17]  Serena Liu,et al.  Faculty Opinions recommendation of Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells. , 2016 .

[18]  L. Pelkmans,et al.  Control of Transcript Variability in Single Mammalian Cells , 2015, Cell.

[19]  Fabian J. Theis,et al.  Data2Dynamics: a modeling environment tailored to parameter estimation in dynamical systems , 2015, Bioinform..

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

[21]  J. Zikherman,et al.  The role of T cell receptor signaling thresholds in guiding T cell fate decisions. , 2015, Current opinion in immunology.

[22]  Lutz Brusch,et al.  Mathematical modelling of fluid transport and its regulation at multiple scales , 2015, Biosyst..

[23]  Fabian J Theis,et al.  Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells , 2015, Nature Biotechnology.

[24]  Ursula Klingmüller,et al.  Identification of isoform-specific dynamics in phosphorylation-dependent STAT5 dimerization by quantitative mass spectrometry and mathematical modeling. , 2014, Journal of proteome research.

[25]  Paul C. Bressloff,et al.  Stochastic Processes in Cell Biology , 2014, Interdisciplinary Applied Mathematics.

[26]  B. Kuster,et al.  Mass-spectrometry-based draft of the human proteome , 2014, Nature.

[27]  Daniel Hidalgo,et al.  Erythropoiesis: from molecular pathways to system properties. , 2014, Advances in experimental medicine and biology.

[28]  Christian Widmer,et al.  Determinants of robustness in spindle assembly checkpoint signalling , 2013, Nature Cell Biology.

[29]  Fabian J Theis,et al.  Lessons Learned from Quantitative Dynamical Modeling in Systems Biology , 2013, PloS one.

[30]  A. Coulon,et al.  Eukaryotic transcriptional dynamics: from single molecules to cell populations , 2013, Nature Reviews Genetics.

[31]  Uwe D. Hanebeck,et al.  Efficient deterministic dirac mixture approximation of Gaussian distributions , 2013, 2013 American Control Conference.

[32]  Bruce Tidor,et al.  Combined Model of Intrinsic and Extrinsic Variability for Computational Network Design with Application to Synthetic Biology , 2013, PLoS Comput. Biol..

[33]  Ursula Klingmüller,et al.  Cellular ERK phospho-form profiles with conserved preference for a switch-like pattern. , 2013, Journal of proteome research.

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

[35]  L. Hennighausen,et al.  Sequential activation of genetic programs in mouse mammary epithelium during pregnancy depends on STAT5A/B concentration , 2012, Nucleic acids research.

[36]  Karen Sachs,et al.  Multiplexed mass cytometry profiling of cellular states perturbed by small-molecule regulators , 2012, Nature Biotechnology.

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

[38]  Lucas Pelkmans,et al.  Using Cell-to-Cell Variability—A New Era in Molecular Biology , 2012, Science.

[39]  Wayne A Moore,et al.  Update for the logicle data scale including operational code implementations , 2012, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[40]  Pradeep S Rajendran,et al.  Single-cell dissection of transcriptional heterogeneity in human colon tumors , 2011, Nature Biotechnology.

[41]  J. Timmer,et al.  Division of labor by dual feedback regulators controls JAK2/STAT5 signaling over broad ligand range , 2011, Molecular systems biology.

[42]  Eva Balsa-Canto,et al.  AMIGO, a toolbox for advanced model identification in systems biology using global optimization , 2011, Bioinform..

[43]  M. Selbach,et al.  Global quantification of mammalian gene expression control , 2011, Nature.

[44]  A Raue,et al.  Identifiability and observability analysis for experimental design in nonlinear dynamical models. , 2010, Chaos.

[45]  Y. Liu,et al.  A Key Commitment Step in Erythropoiesis Is Synchronized with the Cell Cycle Clock through Mutual Inhibition between PU.1 and S-Phase Progression , 2010, PLoS biology.

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

[47]  T. Maiwald,et al.  Materials and Methods SOM Text Figs. S1 to S16 References Materials and Methods , 2022 .

[48]  Claire V. Harper,et al.  Population robustness arising from cellular heterogeneity , 2010, Proceedings of the National Academy of Sciences.

[49]  Lea Goentoro,et al.  Evidence that fold-change, and not absolute level, of beta-catenin dictates Wnt signaling. , 2009, Molecular cell.

[50]  Uri Alon,et al.  Dynamics and variability of ERK2 response to EGF in individual living cells. , 2009, Molecular cell.

[51]  Sabrina L Spencer,et al.  Non-genetic Cell-to-cell Variability and the Consequences for Pharmacology This Review Comes from a Themed Issue on Omics Edited the Distribution of Protein Abundance and Resulting Variability in Phenotype Measuring Cell-to-cell Variation , 2022 .

[52]  P. Liberali,et al.  Population context determines cell-to-cell variability in endocytosis and virus infection , 2009, Nature.

[53]  Sarala M. Wimalaratne,et al.  The Systems Biology Graphical Notation , 2009, Nature Biotechnology.

[54]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

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

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

[57]  Uwe D. Hanebeck,et al.  Localized Cumulative Distributions and a multivariate generalization of the Cramér-von Mises distance , 2008, 2008 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems.

[58]  Jens Timmer,et al.  An error model for protein quantification , 2007, Bioinform..

[59]  O. Rath,et al.  MAP kinase signalling pathways in cancer , 2007, Oncogene.

[60]  Yoh Iwasa,et al.  Robustness of the signal transduction system of the mammalian JAK/STAT pathway and dimerization steps. , 2006, Journal of theoretical biology.

[61]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..

[62]  Leonore A Herzenberg,et al.  Interpreting flow cytometry data: a guide for the perplexed , 2006, Nature Immunology.

[63]  Mario Roederer,et al.  A new “Logicle” display method avoids deceptive effects of logarithmic scaling for low signals and compensated data , 2006, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[64]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[65]  D. Barber,et al.  Turning cells red: signal transduction mediated by erythropoietin. , 2005, Trends in cell biology.

[66]  P. Chattopadhyay,et al.  Seventeen-colour flow cytometry: unravelling the immune system , 2004, Nature Reviews Immunology.

[67]  Rudolph van der Merwe,et al.  Sigma-point kalman filters for probabilistic inference in dynamic state-space models , 2004 .

[68]  B. Neel,et al.  The 'Shp'ing news: SH2 domain-containing tyrosine phosphatases in cell signaling. , 2003, Trends in biochemical sciences.

[69]  J. Timmer,et al.  Identification of nucleocytoplasmic cycling as a remote sensor in cellular signaling by databased modeling , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[70]  H. Lodish,et al.  Ineffective erythropoiesis in Stat5a(-/-)5b(-/-) mice due to decreased survival of early erythroblasts. , 2001, Blood.

[71]  M. H. Nguyen,et al.  A common epitope is shared by activated signal transducer and activator of transcription-5 (STAT5) and the phosphorylated erythropoietin receptor: implications for the docking model of STAT activation. , 2001, Blood.

[72]  A. Yoshimura,et al.  CIS3/SOCS-3 Suppresses Erythropoietin (EPO) Signaling by Binding the EPO Receptor and JAK2* , 2000, The Journal of Biological Chemistry.

[73]  H. Lodish,et al.  Fetal Anemia and Apoptosis of Red Cell Progenitors in Stat5a−/−5b−/− Mice A Direct Role for Stat5 in Bcl-XL Induction , 1999, Cell.

[74]  Douglas M. Bates,et al.  Unconstrained parametrizations for variance-covariance matrices , 1996, Stat. Comput..

[75]  U Klingmüller,et al.  Multiple tyrosine residues in the cytosolic domain of the erythropoietin receptor promote activation of STAT5. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[76]  B. Groner,et al.  Identification of tyrosine residues within the intracellular domain of the erythropoietin receptor crucial for STAT5 activation. , 1996, The EMBO journal.

[77]  Rudolf Jaenisch,et al.  Generation of committed erythroid BFU-E and CFU-E progenitors does not require erythropoietin or the erythropoietin receptor , 1995, Cell.

[78]  N. Copeland,et al.  A novel cytokine‐inducible gene CIS encodes an SH2‐containing protein that binds to tyrosine‐phosphorylated interleukin 3 and erythropoietin receptors. , 1995, The EMBO journal.

[79]  D. Levy,et al.  Prolactin, growth hormone, erythropoietin and granulocyte‐macrophage colony stimulating factor induce MGF‐Stat5 DNA binding activity. , 1995, The EMBO journal.

[80]  M. Koury,et al.  Survival or death of individual proerythroblasts results from differing erythropoietin sensitivities: a mechanism for controlled rates of erythrocyte production. , 1993, Blood.

[81]  M. Potter,et al.  History of the BALB/c family. , 1985, Current topics in microbiology and immunology.

[82]  W. Nijhof,et al.  Cell Kinetic Behaviour of A Synchronized Population of Erythroid Precursor Cells in vitro , 1984, Cell and tissue kinetics.

[83]  W. Nijhof,et al.  Isolation and characterization of the erythroid progenitor cell: CFU-E , 1983, The Journal of cell biology.

[84]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[85]  N. Macdonald,et al.  Time delay in simple chemostat models , 1976, Biotechnology and bioengineering.

[86]  P. Marks,et al.  The effect of erythropoietin on colonial growth of erythroid precursor cells in vitro. , 1974, Proceedings of the National Academy of Sciences of the United States of America.