An Efficient Moments-Based Inference Method for Within-Host Bacterial Infection Dynamics

Abstract

AbstractOver the last ten years, isogenic tagging (IT) has revolutionised the study of bacterial infection dynamics in laboratory animal models. However, quantitative analysis of IT data has been hindered by the piecemeal development of relevant statistical models. The most promising approach relies on stochastic Markovian models of bacterial population dynamics within and among organs. Here we present an efficient numerical method to fit such stochastic dynamic models to in vivo experimental IT data. A common approach to statistical inference with stochastic dynamic models relies on producing large numbers of simulations, but this remains a slow and inefficient method for all but simple problems, especially when tracking bacteria in multiple locations simultaneously. Instead, we derive and solve the systems of ordinary differential equations for the two lower-order moments of the stochastic variables (mean, variance and covariance). For any given model structure, and assuming linear dynamic rates, we demonstrate how the model parameters can be efficiently and accurately estimated by divergence minimisation. We then apply our method to an experimental dataset and compare the estimates and goodness-of-fit to those obtained by maximum likelihood estimation. While both sets of parameter estimates had overlapping confidence regions, the new method produced lower values for the division and death rates of bacteria: these improved the goodness-of-fit at the second time point at the expense of that of the first time point. This flexible framework can easily be applied to a range of experimental systems. Its computational efficiency paves the way for model comparison and optimal experimental design.Author SummaryRecent advancements in technology have meant that microbiologists are producing vast amounts of experimental data. However, statistical methods by which we can analyse that data, draw informative inference, and test relevant hypotheses, are much needed. Here, we present a new, efficient inference tool for estimating parameters of stochastic models, with a particular focus on models of within-host bacterial dynamics. The method relies on matching the two lower-order moments of the experimental data (i.e., mean, variance and covariance), to the moments from the mathematical model. The method is verified, and particular choices justified, through a number of simulation studies. We then use this method to estimate models that have been previously estimated using a “gold-standard” maximum likelihood procedure.List of symbols:A: number of animalsT: number of tagged strainsn: number of organsNi: number of bacteria in organ imij: migration rate from organ i to organ jki: killing rate in organ iri: replication rate in organ iτi: observation time iA, B, C: matricesλ: vector of transition ratesB: Number of bootstrap samplesθ∗: MDE parameter estimateAbbreviations:ABC: approximate Bayesian computationIT: isogenic taggingLV: live vaccineMARE: mean absolute relative errorMDE: minimum divergence estimateMLE: maximum likelihood estimateqPCR: quantitative polymerase chain reactionWITS: wildtype isogenic tagged strain