Thin-film lubrication model for biofilm expansion under strong adhesion.

Understanding microbial biofilm growth is important to public health because biofilms are a leading cause of persistent clinical infections. In this paper, we develop a thin-film model for microbial biofilm growth on a solid substratum to which it adheres strongly. We model biofilms as two-phase viscous fluid mixtures of living cells and extracellular fluid. The model explicitly tracks the movement, depletion, and uptake of nutrients and incorporates cell proliferation via a nutrient-dependent source term. Notably, our thin-film reduction is two dimensional and includes the vertical dependence of cell volume fraction. Numerical solutions show that this vertical dependence is weak for biologically feasible parameters, reinforcing results from previous models in which this dependence was neglected. We exploit this weak dependence by writing and solving a simplified one-dimensional model that is computationally more efficient than the full model. We use both the one- and two-dimensional models to predict how model parameters affect expansion speed and biofilm thickness. This analysis reveals that expansion speed depends on cell proliferation, nutrient availability, cell-cell adhesion on the upper surface, and slip on the biofilm-substratum interface. Our numerical solutions provide a means to qualitatively distinguish between the extensional flow and lubrication regimes, and quantitative predictions that can be tested in future experiments.

Modelling uniaxial non-uniform yeast colony growth: Comparing an agent-based model and continuum approximations.

Biological experiments have shown that yeast can be restricted to grow in a uniaxial direction, vertically upwards from an agar plate to form a colony. The growth occurs as a consequence of cell proliferation driven by a nutrient supply at the base of the colony, and the height of the colony has been observed to increase linearly with time. Within the colony the nutrient concentration is non-constant and yeast cells throughout the colony will therefore not have equal access to nutrient, resulting in non-uniform growth. In this work, an agent based model is developed to predict the microscopic spatial distribution of labelled cells within the colony when the probability of cell proliferation can vary in space and time. We also describe a method for determining the average trajectories or pathlines of labelled cells within a colony growing in a uniaxial direction, enabling us to connect the microscopic and macroscopic behaviours of the system. We present results for six cases, which involve different assumptions for the presence or absence of a quiescent region (where no cell proliferation occurs), the size of the proliferative region, and the spatial variation of proliferation rates within the proliferative region. These six cases are designed to provide qualitative insight into likely growth scenarios whilst remaining amenable to analysis. We compare our macroscopic results to experimental observations of uniaxial colony growth for two cases where only a fixed number of cells at the base of the colony can proliferate. The model predicts that the height of the colony will increase linearly with time in both these cases, which is consistent with experimental observations. However, our model shows how different functional forms for the spatial dependence of the proliferation rate can be distinguished by tracking the pathlines of cells at different positions in the colony. More generally, our methodology can be applied to other biological systems exhibiting uniaxial growth, providing a framework for classifying or determining regions of uniform and non-uniform growth.

A rigid body framework for multicellular modeling.

Off-lattice models are a well-established approach in multicellular modeling, where cells are represented as points that are free to move in space. The representation of cells as point objects is useful in a wide range of settings, particularly when large populations are involved; however, a purely point-based representation is not naturally equipped to deal with objects that have length, such as cell boundaries or external membranes. Here we introduce an off-lattice modeling framework that exploits rigid body mechanics to represent objects using a collection of conjoined one-dimensional edges in a viscosity-dominated system. This framework can be used to represent cells as free moving polygons, to allow epithelial layers to smoothly interact with themselves, to model rod-shaped cells such as bacteria and to robustly represent membranes. We demonstrate that this approach offers solutions to the problems that limit the scope of current off-lattice multicellular models.

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.

A thin-film extensional flow model for biofilm expansion by sliding motility.

In the presence of glycoproteins, bacterial and yeast biofilms are hypothesized to expand by sliding motility. This involves a sheet of cells spreading as a unit, facilitated by cell proliferation and weak adhesion to the substratum. In this paper, we derive an extensional flow model for biofilm expansion by sliding motility to test this hypothesis. We model the biofilm as a two-phase (living cells and an extracellular matrix) viscous fluid mixture, and model nutrient depletion and uptake from the substratum. Applying the thin-film approximation simplifies the model, and reduces it to one-dimensional axisymmetric form. Comparison with Saccharomyces cerevisiae mat formation experiments reveals good agreement between experimental expansion speed and numerical solutions to the model with O ( 1 ) parameters estimated from experiments. This confirms that sliding motility is a possible mechanism for yeast biofilm expansion. Having established the biological relevance of the model, we then demonstrate how the model parameters affect expansion speed, enabling us to predict biofilm expansion for different experimental conditions. Finally, we show that our model can explain the ridge formation observed in some biofilms. This is especially true if surface tension is low, as hypothesized for sliding motility.

Modeling Uniaxial Nonuniform Cell Proliferation.

Growth in biological systems occurs as a consequence of cell proliferation fueled by a nutrient supply. In general, the nutrient gradient of the system will be nonconstant, resulting in biased cell proliferation. We develop a uniaxial discrete cellular automaton with biased cell proliferation using a probability distribution which reflects the nutrient gradient of the system. An explicit probability mass function for the displacement of any tracked cell under the cellular automaton model is derived and verified against averaged simulation results; this displacement distribution has applications in predicting cell trajectories and evolution of expected site occupancies.

TAMMiCol: Tool for analysis of the morphology of microbial colonies.

Many microbes are studied by examining colony morphology via two-dimensional top-down images. The quantification of such images typically requires each pixel to be labelled as belonging to either the colony or background, producing a binary image. While this may be achieved manually for a single colony, this process is infeasible for large datasets containing thousands of images. The software Tool for Analysis of the Morphology of Microbial Colonies (TAMMiCol) has been developed to efficiently and automatically convert colony images to binary. TAMMiCol exploits the structure of the images to choose a thresholding tolerance and produce a binary image of the colony. The images produced are shown to compare favourably with images processed manually, while TAMMiCol is shown to outperform standard segmentation methods. Multiple images may be imported together for batch processing, while the binary data may be exported as a CSV or MATLAB MAT file for quantification, or analysed using statistics built into the software. Using the in-built statistics, it is found that images produced by TAMMiCol yield values close to those computed from binary images processed manually. Analysis of a new large dataset using TAMMiCol shows that colonies of Saccharomyces cerevisiae reach a maximum level of filamentous growth once the concentration of ammonium sulfate is reduced to 200 µM. TAMMiCol is accessed through a graphical user interface, making it easy to use for those without specialist knowledge of image processing, statistical methods or coding.

Characterizing the shape patterns of dimorphic yeast pseudohyphae.

Pseudohyphal growth of the dimorphic yeast Saccharomyces cerevisiae is analysed using two-dimensional top-down binary images. The colony morphology is characterized using clustered shape primitives (CSPs), which are learned automatically from the data and thus do not require a list of predefined features or a priori knowledge of the shape. The power of CSPs is demonstrated through the classification of pseudohyphal yeast colonies known to produce different morphologies. The classifier categorizes the yeast colonies considered with an accuracy of 0.969 and standard deviation 0.041, demonstrating that CSPs capture differences in morphology, while CSPs are found to provide greater discriminatory power than spatial indices previously used to quantify pseudohyphal growth. The analysis demonstrates that CSPs provide a promising avenue for analysing morphology in high-throughput assays.

Diffusion-Limited Growth of Microbial Colonies.

The emergence of diffusion-limited growth (DLG) within a microbial colony on a solid substrate is studied using a combination of mathematical modelling and experiments. Using an agent-based model of the interaction between microbial cells and a diffusing nutrient, it is shown that growth directed towards a nutrient source may be used as an indicator that DLG is influencing the colony morphology. A continuous reaction-diffusion model for microbial growth is employed to identify the parameter regime in which DLG is expected to arise. Comparisons between the model and experimental data are used to argue that the bacterium Bacillus subtilis can undergo DLG, while the yeast Saccharomyces cerevisiae cannot, and thus the non-uniform growth exhibited by this yeast must be caused by the pseudohyphal growth mode rather than limited nutrient availability. Experiments testing directly for DLG features in yeast colonies are used to confirm this hypothesis.

Nutrient-limited growth with non-linear cell diffusion as a mechanism for floral pattern formation in yeast biofilms.

Previous experiments have shown that mature yeast mat biofilms develop a floral morphology, characterised by the formation of petal-like structures. In this work, we investigate the hypothesis that nutrient-limited growth is the mechanism by which these floral patterns form. To do this, we use a combination of experiments and mathematical analysis. In mat formation experiments of the yeast species Saccharomyces cerevisiae, we observe that mats expand radially at a roughly constant speed, and eventually undergo a transition from circular to floral morphology. To determine the extent to which nutrient-limited growth can explain these features, we adopt a previously proposed mathematical model for yeast growth. The model consists of a coupled system of reaction-diffusion equations for the yeast cell density and nutrient concentration, with a non-linear, degenerate diffusion term for cell spread. Using geometric singular perturbation theory and numerics, we show that the model admits travelling wave solutions in one dimension, which enables us to infer the diffusion ratio from experimental data. We then use a linear stability analysis to show that two-dimensional planar travelling wave solutions for feasible experimental parameters are linearly unstable to non-planar perturbations. This provides a potential mechanism by which petals can form, and allows us to predict the characteristic petal width. There is good agreement between these predictions, numerical solutions to the model, and experimental data. We therefore conclude that the non-linear cell diffusion mechanism provides a possible explanation for pattern formation in yeast mat biofilms, without the need to invoke other mechanisms such as flow of extracellular fluid, cell adhesion, or changes to cellular shape or behaviour.

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.

Differential Clonal Expansion in an Invading Cell Population: Clonal Advantage or Dumb Luck?

In neoplastic cell growth, clones and subclones are variable both in size and mutational spectrum. The largest of these clones are believed to represent those cells with mutations that make them the most "fit," in a Darwinian sense, for expansion in their microenvironment. Thus, the degree of quantitative clonal expansion is regarded as being determined by innate qualitative differences between the cells that originate each clone. Here, using a combination of mathematical modelling and clonal labelling experiments applied to the developmental model system of the forming enteric nervous system, we describe how cells which are qualitatively identical may consistently produce clones of dramatically different sizes: most clones are very small while a few clones we term "superstars" contribute most of the cells to the final population. The basis of this is minor stochastic variations ("luck") in the timing and direction of movement and proliferation of individual cells, which builds a local advantage for daughter cells that is cumulative. This has potentially important consequences. In cancers, especially before strongly selective cytotoxic therapy, the assumption that the largest clones must be the cells with deterministic proliferative ability may not always hold true. In development, the gradual loss of clonal diversity as "superstars" take over the population may erode the resilience of the system to somatic mutations, which may have occurred early in clonal growth.

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.

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

Images from cell biology experiments often indicate the presence of cell clustering, which can provide insight into the mechanisms driving the collective cell behaviour. Pair-correlation functions provide quantitative information about the presence, or absence, of clustering in a spatial distribution of cells. This is because the pair-correlation function describes the ratio of the abundance of pairs of cells, separated by a particular distance, relative to a randomly distributed reference population. Pair-correlation functions are often presented as a kernel density estimate where the frequency of pairs of objects are grouped using a particular bandwidth (or bin width), Δ>0. The choice of bandwidth has a dramatic impact: choosing Δ too large produces a pair-correlation function that contains insufficient information, whereas choosing Δ too small produces a pair-correlation signal dominated by fluctuations. Presently, there is little guidance available regarding how to make an objective choice of Δ. We present a new technique to choose Δ by analysing the power spectrum of the discrete Fourier transform of the pair-correlation function. Using synthetic simulation data, we confirm that our approach allows us to objectively choose Δ such that the appropriately binned pair-correlation function captures known features in uniform and clustered synthetic images. We also apply our technique to images from two different cell biology assays. The first assay corresponds to an approximately uniform distribution of cells, while the second assay involves a time series of images of a cell population which forms aggregates over time. The appropriately binned pair-correlation function allows us to make quantitative inferences about the average aggregate size, as well as quantifying how the average aggregate size changes with time.

Incomplete penetrance: The role of stochasticity in developmental cell colonization.

Cell colonization during embryonic development involves cells migrating and proliferating over growing tissues. Unsuccessful colonization, resulting from genetic causes, can result in various birth defects. However not all individuals with the same mutation show the disease. This is termed incomplete penetrance, and it even extends to discordancy in monozygotic (identical) twins. A one-dimensional agent-based model of cell migration and proliferation within a growing tissue is presented, where the position of every cell is recorded at any time. We develop a new model that approximates this agent-based process - rather than requiring the precise configuration of cells within the tissue, the new model records the total number of cells, the position of the most advanced cell, and then invokes an approximation for how the cells are distributed. The probability mass function (PMF) for the most advanced cell is obtained for both the agent-based model and its approximation. The two PMFs compare extremely well, but using the approximation is computationally faster. Success or failure of colonization is probabilistic. For example for sufficiently high proliferation rate the colonization is assured. However, if the proliferation rate is sufficiently low, there will be a lower, say 50%, chance of success. These results provide insights into the puzzle of incomplete penetrance of a disease phenotype, especially in monozygotic twins. Indeed, stochastic cell behavior (amplified by disease-causing mutations) within the colonization process may play a key role in incomplete penetrance, rather than differences in genes, their expression or environmental conditions.

Quantifying two-dimensional filamentous and invasive growth spatial patterns in yeast colonies.

The top-view, two-dimensional spatial patterning of non-uniform growth in a Saccharomyces cerevisiae yeast colony is considered. Experimental images are processed to obtain data sets that provide spatial information on the cell-area that is occupied by the colony. A method is developed that allows for the analysis of the spatial distribution with three metrics. The growth of the colony is quantified in both the radial direction from the centre of the colony and in the angular direction in a prescribed outer region of the colony. It is shown that during the period of 100-200 hours from the start of the growth of the colony there is an increasing amount of non-uniform growth. The statistical framework outlined in this work provides a platform for comparative quantitative assays of strain-specific mechanisms, with potential implementation in inferencing algorithms used for parameter-rate estimation.

Approximating spatially exclusive invasion processes.

A number of biological processes, such as invasive plant species and cell migration, are composed of two key mechanisms: motility and reproduction. Due to the spatially exclusive interacting behavior of these processes a cellular automata (CA) model is specified to simulate a one-dimensional invasion process. Three (independence, Poisson, and 2D-Markov chain) approximations are considered that attempt to capture the average behavior of the CA. We show that our 2D-Markov chain approximation accurately predicts the state of the CA for a wide range of motility and reproduction rates.

Interpreting scratch assays using pair density dynamics and approximate Bayesian computation.

Quantifying the impact of biochemical compounds on collective cell spreading is an essential element of drug design, with various applications including developing treatments for chronic wounds and cancer. Scratch assays are a technically simple and inexpensive method used to study collective cell spreading; however, most previous interpretations of scratch assays are qualitative and do not provide estimates of the cell diffusivity, D, or the cell proliferation rate, λ. Estimating D and λ is important for investigating the efficacy of a potential treatment and provides insight into the mechanism through which the potential treatment acts. While a few methods for estimating D and λ have been proposed, these previous methods lead to point estimates of D and λ, and provide no insight into the uncertainty in these estimates. Here, we compare various types of information that can be extracted from images of a scratch assay, and quantify D and λ using discrete computational simulations and approximate Bayesian computation. We show that it is possible to robustly recover estimates of D and λ from synthetic data, as well as a new set of experimental data. For the first time, our approach also provides a method to estimate the uncertainty in our estimates of D and λ. We anticipate that our approach can be generalized to deal with more realistic experimental scenarios in which we are interested in estimating D and λ, as well as additional relevant parameters such as the strength of cell-to-cell adhesion or the strength of cell-to-substrate adhesion.

Assessing the role of spatial correlations during collective cell spreading.

Spreading cell fronts are essential features of development, repair and disease processes. Many mathematical models used to describe the motion of cell fronts, such as Fisher's equation, invoke a mean-field assumption which implies that there is no spatial structure, such as cell clustering, present. Here, we examine the presence of spatial structure using a combination of in vitro circular barrier assays, discrete random walk simulations and pair correlation functions. In particular, we analyse discrete simulation data using pair correlation functions to show that spatial structure can form in a spreading population of cells either through sufficiently strong cell-to-cell adhesion or sufficiently rapid cell proliferation. We analyse images from a circular barrier assay describing the spreading of a population of MM127 melanoma cells using the same pair correlation functions. Our results indicate that the spreading melanoma cell populations remain very close to spatially uniform, suggesting that the strength of cell-to-cell adhesion and the rate of cell proliferation are both sufficiently small so as not to induce any spatial patterning in the spreading populations.