Analysis of CFSE time-series data using division-, age- and label-structured population models

MOTIVATION In vitro and in vivo cell proliferation is often studied using the dye carboxyfluorescein succinimidyl ester (CFSE). The CFSE time-series data provide information about the proliferation history of populations of cells. While the experimental procedures are well established and widely used, the analysis of CFSE time-series data is still challenging. Many available analysis tools do not account for cell age and employ optimization methods that are inefficient (or even unreliable). RESULTS We present a new model-based analysis method for CFSE time-series data. This method uses a flexible description of proliferating cell populations, namely, a division-, age- and label-structured population model. Efficient maximum likelihood and Bayesian estimation algorithms are introduced to infer the model parameters and their uncertainties. These methods exploit the forward sensitivity equations of the underlying partial differential equation model for efficient and accurate gradient calculation, thereby improving computational efficiency and reliability compared with alternative approaches and accelerating uncertainty analysis. The performance of the method is assessed by studying a dataset for immune cell proliferation. This revealed the importance of different factors on the proliferation rates of individual cells. Among others, the predominate effect of cell age on the division rate is found, which was not revealed by available computational methods. AVAILABILITY AND IMPLEMENTATION The MATLAB source code implementing the models and algorithms is available from http://janhasenauer.github.io/ShAPE-DALSP/Contact: jan.hasenauer@helmholtz-muenchen.de SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  A. Hoffmann,et al.  A multi-scale approach reveals that NF-κB cRel enforces a B-cell decision to divide , 2015, Molecular systems biology.

[2]  John F Markham,et al.  Measuring lymphocyte proliferation, survival and differentiation using CFSE time-series data , 2007, Nature Protocols.

[3]  A. B. Lyons,et al.  Determination of lymphocyte division by flow cytometry. , 1994, Journal of immunological methods.

[4]  Jan Hasenauer,et al.  Modeling and parameter estimation for heterogeneous cell populations , 2013 .

[5]  D. Roose,et al.  Computational analysis of CFSE proliferation assay , 2006, Journal of mathematical biology.

[6]  H. T. Banks,et al.  Quantifying CFSE label decay in flow cytometry data , 2013, Appl. Math. Lett..

[7]  J. Smith,et al.  Do cells cycle? , 1973, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Alan S Perelson,et al.  Quantifying T lymphocyte turnover. , 2013, Journal of theoretical biology.

[9]  F. Allgöwer,et al.  Analysis and Simulation of Division- and Label-Structured Population Models , 2012, Bulletin of mathematical biology.

[10]  L. Fenton The Sum of Log-Normal Probability Distributions in Scatter Transmission Systems , 1960 .

[11]  Chyung-Ru Wang,et al.  Helper T cell differentiation is controlled by the cell cycle. , 1998, Immunity.

[12]  G. Bocharov,et al.  Asymmetry of Cell Division in CFSE-Based Lymphocyte Proliferation Analysis , 2013, Front. Immunol..

[13]  Karyn L. Sutton,et al.  A new model for the estimation of cell proliferation dynamics using CFSE data. , 2011, Journal of immunological methods.

[14]  Karel Svoboda,et al.  From cudgel to scalpel: toward precise neural control with optogenetics , 2011, Nature Methods.

[15]  Heikki Haario,et al.  DRAM: Efficient adaptive MCMC , 2006, Stat. Comput..

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

[18]  Dirk Roose,et al.  Numerical modelling of label-structured cell population growth using CFSE distribution data , 2007, Theoretical Biology and Medical Modelling.

[19]  Timm Schroeder,et al.  Long-term single-cell imaging of mammalian stem cells , 2011, Nature Methods.

[20]  Luís N. Vicente,et al.  A particle swarm pattern search method for bound constrained global optimization , 2007, J. Glob. Optim..

[21]  Stephen P. Brooks,et al.  Assessing Convergence of Markov Chain Monte Carlo Algorithms , 2007 .

[22]  H. Banks,et al.  A novel statistical analysis and interpretation of flow cytometry data , 2013, Journal of biological dynamics.

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

[24]  Andreas Meyerhans,et al.  Analysis of variability in estimates of cell proliferation parameters for cyton-based models using CFSE-based flow cytometry data , 2015 .

[25]  Dirk Roose,et al.  Estimation of Cell Proliferation Dynamics Using CFSE Data , 2011, Bulletin of mathematical biology.

[26]  Laurent Pujo-Menjouet,et al.  Analysis of cell kinetics using a cell division marker: mathematical modeling of experimental data. , 2003, Biophysical journal.

[27]  Burkhard Ludewig,et al.  Mathematical models for CFSE labelled lymphocyte dynamics: asymmetry and time-lag in division , 2014, Journal of mathematical biology.

[28]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[29]  D. Kapraun Cell Proliferation Models, CFSE-Based Flow Cytometry Data, and Quantification of Uncertainty. , 2014 .

[30]  H. Banks,et al.  Evaluating the importance of mitotic asymmetry in cyton-based models for FSE-based flow cytometry data , 2015 .

[31]  Alan S. Perelson,et al.  Estimating Lymphocyte Division and Death Rates from CFSE Data , 2006, Bulletin of mathematical biology.

[32]  P. Hodgkin,et al.  A model of immune regulation as a consequence of randomized lymphocyte division and death times , 2007, Proceedings of the National Academy of Sciences.

[33]  Gareth O. Roberts,et al.  Convergence assessment techniques for Markov chain Monte Carlo , 1998, Stat. Comput..

[34]  Fabian J Theis,et al.  High-dimensional Bayesian parameter estimation: case study for a model of JAK2/STAT5 signaling. , 2013, Mathematical biosciences.

[35]  Jhagvaral Hasbold,et al.  Activation-Induced B Cell Fates Are Selected by Intracellular Stochastic Competition , 2012, Science.