QUANTIFYING INTRINSIC AND EXTRINSIC NOISE IN GENE TRANSCRIPTION USING THE LINEAR NOISE APPROXIMATION: AN APPLICATION TO SINGLE CELL DATA

被引:27
作者
Finkenstadt, Barbel [1 ]
Woodcock, Dan J. [2 ]
Komorowski, Michal [3 ]
Harper, Claire V. [4 ]
Davis, Julian R. E. [5 ]
White, Mike R. H. [4 ]
Rand, David A. [2 ]
机构
[1] Univ Warwick, Dept Stat, Coventry CV4 7AL, W Midlands, England
[2] Univ Warwick, Syst Biol Ctr, Coventry CV4 7AL, W Midlands, England
[3] Polish Acad Sci, Div Modelling Biol & Med, PL-02106 Warsaw, Poland
[4] Univ Manchester, Fac Life Sci, Manchester M13 9PT, Lancs, England
[5] Univ Manchester, Fac Med & Human Sci, Manchester M13 9PT, Lancs, England
基金
英国生物技术与生命科学研究理事会;
关键词
Linear noise approximation; kinetic parameter estimation; intrinsic and extrinsic noise; state space model and Kalman filter; Bayesian hierarchical modeling; BAYESIAN-INFERENCE; STOCHASTIC SIMULATION; MODELS; FLUCTUATIONS;
D O I
10.1214/13-AOAS669
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
A central challenge in computational modeling of dynamic biological systems is parameter inference from experimental time course measurements. However, one would not only like to infer kinetic parameters but also study their variability from cell to cell. Here we focus on the case where single-cell fluorescent protein imaging time series data are available for a population of cells. Based on van Kampen's linear noise approximation, we derive a dynamic state space model for molecular populations which is then extended to a hierarchical model. This model has potential to address the sources of variability relevant to single-cell data, namely, intrinsic noise due to the stochastic nature of the birth and death processes involved in reactions and extrinsic noise arising from the cell-to-cell variation of kinetic parameters. In order to infer such a model from experimental data, one must also quantify the measurement process where one has to allow for nonmeasurable molecular species as well as measurement noise of unknown level and variance. The availability of multiple single-cell time series data here provides a unique testbed to fit such a model and quantify these different sources of variation from experimental data.
引用
收藏
页码:1960 / 1982
页数:23
相关论文
共 36 条
  • [1] [Anonymous], 2006, MARKOV CHAIN MONTE C
  • [2] [Anonymous], 2011, Stochastic modelling for systems biology
  • [3] UNDERSTANDING THE METROPOLIS-HASTINGS ALGORITHM
    CHIB, S
    GREENBERG, E
    [J]. AMERICAN STATISTICIAN, 1995, 49 (04) : 327 - 335
  • [4] Acceleration of the Multiple-Try Metropolis algorithm using antithetic and stratified sampling
    Craiu, Radu V.
    Lemieux, Christiane
    [J]. STATISTICS AND COMPUTING, 2007, 17 (02) : 109 - 120
  • [5] Numerical techniques for maximum likelihood estimation of continuous-time diffusion processes
    Durham, GB
    Gallant, AR
    [J]. JOURNAL OF BUSINESS & ECONOMIC STATISTICS, 2002, 20 (03) : 297 - 316
  • [6] Likelihood inference for discretely observed nonlinear diffusions
    Elerian, O
    Chib, S
    Shephard, N
    [J]. ECONOMETRICA, 2001, 69 (04) : 959 - 993
  • [7] Fast evaluation of fluctuations in biochemical networks with the linear noise approximation
    Elf, J
    Ehrenberg, M
    [J]. GENOME RESEARCH, 2003, 13 (11) : 2475 - 2484
  • [8] Stochastic gene expression in a single cell
    Elowitz, MB
    Levine, AJ
    Siggia, ED
    Swain, PS
    [J]. SCIENCE, 2002, 297 (5584) : 1183 - 1186
  • [9] FEARNHEAD P., 2013, INFERENCE REACTION N
  • [10] Reconstruction of transcriptional dynamics from gene reporter data using differential equations
    Finkenstadt, Barbel
    Heron, Elizabeth A.
    Komorowski, Michal
    Edwards, Kieron
    Tang, Sanyi
    Harper, Claire V.
    Davis, Julian R. E.
    White, Michael R. H.
    Millar, Andrew J.
    Rand, David A.
    [J]. BIOINFORMATICS, 2008, 24 (24) : 2901 - 2907