Stuart T. Johnston

How much information can be obtained from tracking the position of the leading edge in a scratch assay

Moving cell fronts are an essential feature of wound healing, development and disease. The rate at which a cell front moves is driven, in part, by the cell motility, quantified in terms of the cell diffusivity D, and the cell proliferation rate λ. Scratch assays are a commonly reported procedure used to investigate the motion of cell fronts where an initial cell monolayer is scratched, and the motion of the front is monitored over a short period of time, often less than 24 h. The simplest way of quantifying a scratch assay is to monitor the progression of the leading edge. Use of leading edge data is very convenient because, unlike other methods, it is non-destructive and does not require labelling, tracking or counting individual cells among the population. In this work, we study short-time leading edge data in a scratch assay using a discrete mathematical model and automated image analysis with the aim of investigating whether such data allow us to reliably identify D and λ. Using a naive calibration approach where we simply scan the relevant region of the (D, λ) parameter space, we show that there are many choices of D and λ for which our model produces indistinguishable short-time leading edge data. Therefore, without due care, it is impossible to estimate D and λ from this kind of data. To address this, we present a modified approach accounting for the fact that cell motility occurs over a much shorter time scale than proliferation. Using this information, we divide the duration of the experiment into two periods, and we estimate D using data from the first period, whereas we estimate λ using data from the second period. We confirm the accuracy of our approach using in silico data and a new set of in vitro data, which shows that our method recovers estimates of D and λ that are consistent with previously reported values except that that our approach is fast, inexpensive, non-destructive and avoids the need for cell labelling and cell counting.

Estimating cell diffusivity and cell proliferation rate by interpreting IncuCyte ZOOM™ assay data using the Fisher-Kolmogorov model

BackgroundStandard methods for quantifying IncuCyte ZOOM™ assays involve measurements that quantify how rapidly the initially-vacant area becomes re-colonised with cells as a function of time. Unfortunately, these measurements give no insight into the details of the cellular-level mechanisms acting to close the initially-vacant area. We provide an alternative method enabling us to quantify the role of cell motility and cell proliferation separately. To achieve this we calibrate standard data available from IncuCyte ZOOM™ images to the solution of the Fisher-Kolmogorov model.ResultsThe Fisher-Kolmogorov model is a reaction-diffusion equation that has been used to describe collective cell spreading driven by cell migration, characterised by a cell diffusivity, D, and carrying capacity limited proliferation with proliferation rate, λ, and carrying capacity density, K. By analysing temporal changes in cell density in several subregions located well-behind the initial position of the leading edge we estimate λ and K. Given these estimates, we then apply automatic leading edge detection algorithms to the images produced by the IncuCyte ZOOM™ assay and match this data with a numerical solution of the Fisher-Kolmogorov equation to provide an estimate of D. We demonstrate this method by applying it to interpret a suite of IncuCyte ZOOM™ assays using PC-3 prostate cancer cells and obtain estimates of D, λ and K. Comparing estimates of D, λ and K for a control assay with estimates of D, λ and K for assays where epidermal growth factor (EGF) is applied in varying concentrations confirms that EGF enhances the rate of scratch closure and that this stimulation is driven by an increase in D and λ, whereas K is relatively unaffected by EGF.ConclusionsOur approach for estimating D, λ and K from an IncuCyte ZOOM™ assay provides more detail about cellular-level behaviour than standard methods for analysing these assays. In particular, our approach can be used to quantify the balance of cell migration and cell proliferation and, as we demonstrate, allow us to quantify how the addition of growth factors affects these processes individually.

Modelling the movement of interacting cell populations: a moment dynamics approach.

Mathematical models describing the movement of multiple interacting subpopulations are relevant to many biological and ecological processes. Standard mean-field partial differential equation descriptions of these processes suffer from the limitation that they implicitly neglect to incorporate the impact of spatial correlations and clustering. To overcome this, we derive a moment dynamics description of a discrete stochastic process which describes the spreading of distinct interacting subpopulations. In particular, we motivate our model by mimicking the geometry of two typical cell biology experiments. Comparing the performance of the moment dynamics model with a traditional mean-field model confirms that the moment dynamics approach always outperforms the traditional mean-field approach. To provide more general insight we summarise the performance of the moment dynamics model and the traditional mean-field model over a wide range of parameter regimes. These results help distinguish between those situations where spatial correlation effects are sufficiently strong, such that a moment dynamics model is required, from other situations where spatial correlation effects are sufficiently weak, such that a traditional mean-field model is adequate.

Quantifying the effect of experimental design choices for in vitro scratch assays.

Scratch assays are often used to investigate potential drug treatments for chronic wounds and cancer. Interpreting these experiments with a mathematical model allows us to estimate the cell diffusivity, D, and the cell proliferation rate, λ. However, the influence of the experimental design on the estimates of D and λ is unclear. Here we apply an approximate Bayesian computation (ABC) parameter inference method, which produces a posterior distribution of D and λ, to new sets of synthetic data, generated from an idealised mathematical model, and experimental data for a non-adhesive mesenchymal population of fibroblast cells. The posterior distribution allows us to quantify the amount of information obtained about D and λ. We investigate two types of scratch assay, as well as varying the number and timing of the experimental observations captured. Our results show that a scrape assay, involving one cell front, provides more precise estimates of D and λ, and is more computationally efficient to interpret than a wound assay, with two opposingly directed cell fronts. We find that recording two observations, after making the initial observation, is sufficient to estimate D and λ, and that the final observation time should correspond to the time taken for the cell front to move across the field of view. These results provide guidance for estimating D and λ, while simultaneously minimising the time and cost associated with performing and interpreting the experiment.

Co-operation, Competition and Crowding: A Discrete Framework Linking Allee Kinetics, Nonlinear Diffusion, Shocks and Sharp-Fronted Travelling Waves

Invasion processes are ubiquitous throughout cell biology and ecology. During invasion, individuals can become isolated from the bulk population and behave differently. We present a discrete, exclusion-based description of the birth, death and movement of individuals. The model distinguishes between individuals that are part of, or are isolated from, the bulk population by imposing different rates of birth, death and movement. This enables the simulation of various co-operative or competitive mechanisms, where there is either a positive or negative benefit associated with being part of the bulk population, respectively. The mean-field approximation of the discrete process gives rise to 22 different classes of partial differential equation, which can include Allee kinetics and nonlinear diffusion. Here we examine the ability of each class of partial differential equation to support travelling wave solutions and interpret the long time behaviour in terms of the individual-level parameters. For the first time we show that the strong Allee effect and nonlinear diffusion can result in shock-fronted travelling waves. We also demonstrate how differences in group and individual motility rates can influence the persistence of a population and provide conditions for the successful invasion of a population.

Filling the gaps: A robust description of adhesive birth-death-movement processes.

Existing continuum descriptions of discrete adhesive birth-death-movement processes provide accurate predictions of the average discrete behavior for limited parameter regimes. Here we present an alternative continuum description in terms of the dynamics of groups of contiguous occupied and vacant lattice sites. Our method provides more accurate predictions, is valid in parameter regimes that could not be described by previous continuum descriptions, and provides information about the spatial clustering of occupied sites. Furthermore, we present a simple analytic approximation of the spatial clustering of occupied sites at late time, when the system reaches its steady-state configuration.