Spectral analysis of pair-correlation bandwidth: application to cell biology images

Calibrating tumor growth and invasion parameters with spectral spatial analysis of cancer biopsy tissues

The reaction-diffusion equation is widely used in mathematical models of cancer. The calibration of model parameters based on limited clinical data is critical to using reaction-diffusion equation simulations for reliable predictions on a per-patient basis. Here, we focus on cell-level data as routinely available from tissue biopsies used for clinical cancer diagnosis. We analyze the spatial architecture in biopsy tissues stained with multiplex immunofluorescence. We derive a two-point correlation function and the corresponding spatial power spectral distribution. We show that this data-deduced power spectral distribution can fit the power spectrum of the solution of reaction-diffusion equations that can then identify patient-specific tumor growth and invasion rates. This approach allows the measurement of patient-specific critical tumor dynamical properties from routinely available biopsy material at a single snapshot in time.

Detection and characterization of chemotaxis without cell tracking

Swarming has been observed in various biological systems from collective animal movements to immune cells. In the cellular context, swarming is driven by the secretion of chemotactic factors. Despite the critical role of chemotactic swarming, few methods to robustly identify and quantify this phenomenon exist. Here, we present a novel method for the analysis of time series of positional data generated from realizations of agent-based processes. We convert the positional data for each individual time point to a function measuring agent aggregation around a given area of interest, hence generating a functional time series. The functional time series, and a more easily visualized swarming metric of agent aggregation derived from these functions, provide useful information regarding the evolution of the underlying process over time. We extend our method to build upon the modelling of collective motility using drift-diffusion partial differential equations (PDEs). Using a functional linear model, we are able to use the functional time series to estimate the drift and diffusivity terms associated with the underlying PDE. By producing an accurate estimate for the drift coefficient, we can infer the strength and range of attraction or repulsion exerted on agents, as in chemotaxis. Our approach relies solely on using agent positional data. The spatial distribution of diffusing chemokines is not required, nor do individual agents need to be tracked over time. We demonstrate our approach using random walk simulations of chemotaxis and experiments investigating cytotoxic T cells interacting with tumouroids.

Discrete Manhattan and Chebyshev pair correlation functions in k dimensions

Pair correlation functions provide a summary statistic which quantifies the amount of spatial correlation between objects in a spatial domain. While pair correlation functions are commonly used to quantify continuous-space point processes, the on-lattice discrete case is less studied. Recent work has brought attention to the discrete case, wherein on-lattice pair correlation functions are formed by normalizing empirical pair distances against the probability distribution of random pair distances in a lattice with Manhattan and Chebyshev metrics. These distance distributions are typically derived on an ad hoc basis as required for specific applications. Here we present a generalized approach to deriving the probability distributions of pair distances in a lattice with discrete Manhattan and Chebyshev metrics, extending the Manhattan and Chebyshev pair correlation functions to lattices in k dimensions. We also quantify the variability of the Manhattan and Chebyshev pair correlation functions, which is important to understanding the reliability and confidence of the statistic.

Identifying density-dependent interactions in collective cell behaviour

Scratch assays are routinely used to study collective cell behaviour in vitro. Typical experimental protocols do not vary the initial density of cells, and typical mathematical modelling approaches describe cell motility and proliferation based on assumptions of linear diffusion and logistic growth. Jin et al. (Jin et al. 2016 J. Theor. Biol. 390, 136-145 (doi:10.1016/j.jtbi.2015.10.040)) find that the behaviour of cells in scratch assays is density-dependent, and show that standard modelling approaches cannot simultaneously describe data initiated across a range of initial densities. To address this limitation, we calibrate an individual-based model to scratch assay data across a large range of initial densities. Our model allows proliferation, motility, and a direction bias to depend on interactions between neighbouring cells. By considering a hierarchy of models where we systematically and sequentially remove interactions, we perform model selection analysis to identify the minimum interactions required for the model to simultaneously describe data across all initial densities. The calibrated model is able to match the experimental data across all densities using a single parameter distribution, and captures details about the spatial structure of cells. Our results provide strong evidence to suggest that motility is density-dependent in these experiments. On the other hand, we do not see the effect of crowding on proliferation in these experiments. These results are significant as they are precisely the opposite of the assumptions in standard continuum models, such as the Fisher-Kolmogorov equation and its generalizations.

Corrected pair correlation functions for environments with obstacles

Environments with immobile obstacles or void regions that inhibit and alter the motion of individuals within that environment are ubiquitous. Correlation in the location of individuals within such environments arises as a combination of the mechanisms governing individual behavior and the heterogeneous structure of the environment. Measures of spatial structure and correlation have been successfully implemented to elucidate the roles of the mechanisms underpinning the behavior of individuals. In particular, the pair correlation function has been used across biology, ecology, and physics to obtain quantitative insight into a variety of processes. However, naively applying standard pair correlation functions in the presence of obstacles may fail to detect correlation, or suggest false correlations, due to a reliance on a distance metric that does not account for obstacles. To overcome this problem, here we present an analytic expression for calculating a corrected pair correlation function for lattice-based domains containing obstacles. We demonstrate that this obstacle pair correlation function is necessary for isolating the correlation associated with the behavior of individuals, rather than the structure of the environment. Using simulations that mimic cell migration and proliferation we demonstrate that the obstacle pair correlation function recovers the short-range correlation known to be present in this process, independent of the heterogeneous structure of the environment. Further, we show that the analytic calculation of the obstacle pair correlation function derived here is significantly faster to implement than the corresponding numerical approach.

Pair correlation functions for identifying spatial correlation in discrete domains

Identifying and quantifying spatial correlation are important aspects of studying the collective behavior of multiagent systems. Pair correlation functions (PCFs) are powerful statistical tools that can provide qualitative and quantitative information about correlation between pairs of agents. Despite the numerous PCFs defined for off-lattice domains, only a few recent studies have considered a PCF for discrete domains. Our work extends the study of spatial correlation in discrete domains by defining a new set of PCFs using two natural and intuitive definitions of distance for a square lattice: the taxicab and uniform metric. We show how these PCFs improve upon previous attempts and compare between the quantitative data acquired. We also extend our definitions of the PCF to other types of regular tessellation that have not been studied before, including hexagonal, triangular, and cuboidal. Finally, we provide a comprehensive PCF for any tessellation and metric, allowing investigation of spatial correlation in irregular lattices for which recognizing correlation is less intuitive.

Quantifying the dominant growth mechanisms of dimorphic yeast using a lattice-based model

A mathematical model is presented for the growth of yeast that incorporates both dimorphic behaviour and nutrient diffusion. The budding patterns observed in the standard and pseudohyphal growth modes are represented by a bias in the direction of cell proliferation. A set of spatial indices is developed to quantify the morphology and compare the relative importance of the directional bias to nutrient concentration and diffusivity on colony shape. It is found that there are three different growth modes: uniform growth, diffusion-limited growth (DLG) and an intermediate region in which the bias determines the morphology. The dimorphic transition due to nutrient limitation is investigated by relating the directional bias to the nutrient concentration, and this is shown to replicate the behaviour observed in vivo Comparisons are made with experimental data, from which it is found that the model captures many of the observed features. Both DLG and pseudohyphal growth are found to be capable of generating observed experimental morphologies.

Variable species densities are induced by volume exclusion interactions upon domain growth

In this work we study the effect of domain growth on spatial correlations in agent populations containing multiple species. This is important as heterogenous cell populations are ubiquitous during the embryonic development of many species. We have previously shown that the long-term behavior of an agent population depends on the way in which domain growth is implemented. We extend this work to show that, depending on the way in which domain growth is implemented, different species dominate in multispecies simulations. Continuum approximations of the lattice-based model that ignore spatial correlations cannot capture this behavior, while those that explicitly account for spatial correlations can. The results presented here show that the precise mechanism of domain growth can determine the long-term behavior of multispecies populations and, in certain circumstances, establish spatially varying species densities.

Spatial moment dynamics for collective cell movement incorporating a neighbour-dependent directional bias

The ability of cells to undergo collective movement plays a fundamental role in tissue repair, development and cancer. Interactions occurring at the level of individual cells may lead to the development of spatial structure which will affect the dynamics of migrating cells at a population level. Models that try to predict population-level behaviour often take a mean-field approach, which assumes that individuals interact with one another in proportion to their average density and ignores the presence of any small-scale spatial structure. In this work, we develop a lattice-free individual-based model (IBM) that uses random walk theory to model the stochastic interactions occurring at the scale of individual migrating cells. We incorporate a mechanism for local directional bias such that an individual's direction of movement is dependent on the degree of cell crowding in its neighbourhood. As an alternative to the mean-field approach, we also employ spatial moment theory to develop a population-level model which accounts for spatial structure and predicts how these individual-level interactions propagate to the scale of the whole population. The IBM is used to derive an equation for dynamics of the second spatial moment (the average density of pairs of cells) which incorporates the neighbour-dependent directional bias, and we solve this numerically for a spatially homogeneous case.

Spatial structure arising from neighbour-dependent bias in collective cell movement

Mathematical models of collective cell movement often neglect the effects of spatial structure, such as clustering, on the population dynamics. Typically, they assume that individuals interact with one another in proportion to their average density (the mean-field assumption) which means that cell-cell interactions occurring over short spatial ranges are not accounted for. However, in vitro cell culture studies have shown that spatial correlations can play an important role in determining collective behaviour. Here, we take a combined experimental and modelling approach to explore how individual-level interactions give rise to spatial structure in a moving cell population. Using imaging data from in vitro experiments, we quantify the extent of spatial structure in a population of 3T3 fibroblast cells. To understand how this spatial structure arises, we develop a lattice-free individual-based model (IBM) and simulate cell movement in two spatial dimensions. Our model allows an individual's direction of movement to be affected by interactions with other cells in its neighbourhood, providing insights into how directional bias generates spatial structure. We consider how this behaviour scales up to the population level by using the IBM to derive a continuum description in terms of the dynamics of spatial moments. In particular, we account for spatial correlations between cells by considering dynamics of the second spatial moment (the average density of pairs of cells). Our numerical results suggest that the moment dynamics description can provide a good approximation to averaged simulation results from the underlying IBM. Using our in vitro data, we estimate parameters for the model and show that it can generate similar spatial structure to that observed in a 3T3 fibroblast cell population.