A two-phase thin-film model for cell-induced gel contraction incorporating osmotic effects.

We present a mathematical model of an experiment in which cells are cultured within a gel, which in turn floats freely within a liquid nutrient medium. Traction forces exerted by the cells on the gel cause it to contract over time, giving a measure of the strength of these forces. Building upon our previous work (Reoch et al. in J Math Biol 84(5):31, 2022), we exploit the fact that the gels used frequently have a thin geometry to obtain a reduced model for the behaviour of a thin, two-dimensional cell-seeded gel. We find that steady-state solutions of the reduced model require the cell density and volume fraction of polymer in the gel to be spatially uniform, while the gel height may vary spatially. If we further assume that all three of these variables are initially spatially uniform, this continues for all time and the thin film model can be further reduced to solving a single, non-linear ODE for gel height as a function of time. The thin film model is further investigated for both spatially-uniform and varying initial conditions, using a combination of analytical techniques and numerical simulations. We show that a number of qualitatively different behaviours are possible, depending on the composition of the gel (i.e., the chemical potentials) and the strength of the cell traction forces. However, unlike in the earlier one-dimensional model, we do not observe cases where the gel oscillates between swelling and contraction. For the case of initially uniform cell and gel density, our model predicts that the relative change in the gels' height and length are equal, which justifies an assumption previously used in the work of Stevenson et al. (Biophys J 99(1):19-28, 2010). Conversely, however, even for non-uniform initial conditions, we do not observe cases where the length of the gel changes whilst its height remains constant, which have been reported in another model of osmotic swelling by Trinschek et al. (AIMS Mater Sci 3(3):1138-1159, 2016; Phys Rev Lett 119:078003, 2017).

A mathematical model for nutrient-limited uniaxial growth of a compressible tissue.

We consider the uniaxial growth of a tissue or colony of cells, where a nutrient (or some other chemical) required for cell proliferation is supplied at one end, and is consumed by the cells. An example would be the growth of a cylindrical yeast colony in the experiments described by Vulin et al. (2014). We develop a reaction-diffusion model of this scenario which couples nutrient concentration and cell density on a growing domain. A novel element of our model is that the tissue is assumed to be compressible. We define replicative regions, where cells have sufficient nutrient to proliferate, and quiescent regions, where the nutrient level is insufficient for this to occur. We also define pathlines, which allow us to track individual cell paths within the tissue. We begin our investigation of the model by considering an incompressible tissue where cell density is constant before exploring the solution space of the full compressible model. In a large part of the parameter space, the incompressible and compressible models give qualitatively similar results for both the nutrient concentration and cell pathlines, with the key distinction being the variation in density in the compressible case. In particular, the replicative region is located at the base of the tissue, where nutrient is supplied, and nutrient concentration decreases monotonically with distance from the nutrient source. However, for a highly-compressible tissue with small nutrient consumption rate, we observe a counter-intuitive scenario where the nutrient concentration is not necessarily monotonically decreasing, and there can be two replicative regions. For parameter values given in the paper by Vulin et al. (2014), the incompressible model slightly overestimates the colony length compared to experimental observations; this suggests the colony may be somewhat compressible. Both incompressible and compressible models predict that, for these parameter values, cell proliferation is ultimately confined to a small region close to the colony base.

Modelling foraging competition between solitarious and gregarious organisms in increasingly heterogeneous environments.

Density dependent phase polyphenism is the exhibiting of two or more distinct phenotypes from a single genotype depending on local population density. The most well known insect to exhibit this phenomenon is the locust, with whom the profound effect on behaviour leads to the classification of the two phases; solitarious, where locusts actively avoid other locusts, and gregarious, where locusts are strongly attracted to other locusts. It has been shown that food distributions at both small and large scales have an effect on the process of gregarisation. While gregarisation offers advantages, such as greater predator avoidance, the relationship between phase polyphenism and potential foraging benefits is still not fully understood. In this paper, we explore the effect of gregarisation on foraging within increasingly heterogeneous environments using a partial differential equation model. We first consider a single two dimensional simulation of a spatially heterogeneous environment to understand the mechanics of gregarious/solitarious foraging. We then look at the steady state foraging advantage (measured as the ratio of per-capita contact with food) in environments ranging from homogeneous to very spatially heterogeneous. Finally, we perform a parameter sensitivity analysis to find which model parameters have the greatest effect on foraging advantage. We find that during the aggregation stage, prior to the onset of marching (which we do not model here), in increasingly heterogeneous food environments it is better to be gregarious than solitarious. In addition, we find that this is intrinsic to the gregarious/solitarious behavioural dynamic as it occurs almost regardless of the model parameters. That is to say, it doesn't matter how fast the organisms disperse or how strong their long range interactions as long as there is the solitarious/gregarious behaviour the gregarious foraging advantage will exist.

A mathematical model for cell-induced gel contraction incorporating osmotic effects.

Biological tissues are composed of cells surrounded by the extracellular matrix (ECM). The ECM can be thought of as a fibrous polymer network, acting as a natural scaffolding to provide mechanical support to the cells. Reciprocal mechanical and chemical interactions between the cells and the ECM are crucial in regulating the development of tissues and maintaining their functionality. Hence, to maintain in vivo-like behaviour when cells are cultured in vitro, they are often seeded in a gel, which aims to mimic the ECM. In this paper, we present a mathematical model that incorporates cell-gel interactions together with osmotic pressure to study the mechanical behaviour of biological gels. In particular, we consider an experiment where cells are seeded within a gel, which gradually compacts due to forces exerted on it by the cells. Adopting a one-dimensional Cartesian geometry for simplicity, we use a combination of analytical techniques and numerical simulations to investigate how cell traction forces interact with osmotic effects (which can lead to either gel swelling or contraction depending on the gel's composition). Our results show that a number of qualitatively different behaviours are possible, depending on the composition of the gel (i.e. its chemical potentials) and the strength of the cell traction forces. A novel prediction of our model is that there are cases where the gel oscillates between swelling and contraction; to our knowledge, this behaviour has not been reported in experiments. We also consider how physical parameters like drag and viscosity affect the manner in which the gel evolves.

Correction: Modelling locust foraging: How and why food affects group formation.

[This corrects the article DOI: 10.1371/journal.pcbi.1008353.].

Modelling locust foraging: How and why food affects group formation.

Locusts are short horned grasshoppers that exhibit two behaviour types depending on their local population density. These are: solitarious, where they will actively avoid other locusts, and gregarious where they will seek them out. It is in this gregarious state that locusts can form massive and destructive flying swarms or plagues. However, these swarms are usually preceded by the aggregation of juvenile wingless locust nymphs. In this paper we attempt to understand how the distribution of food resources affect the group formation process. We do this by introducing a multi-population partial differential equation model that includes non-local locust interactions, local locust and food interactions, and gregarisation. Our results suggest that, food acts to increase the maximum density of locust groups, lowers the percentage of the population that needs to be gregarious for group formation, and decreases both the required density of locusts and time for group formation around an optimal food width. Finally, by looking at foraging efficiency within the numerical experiments we find that there exists a foraging advantage to being gregarious.

Influences of transversely isotropic rheology and translational diffusion on the stability of active suspensions.

Suspensions of self-motile, elongated particles are a topic of significant current interest, exemplifying a form of 'active matter'. Examples include self-propelling bacteria, algae and sperm, and artificial swimmers. Ericksen's model of a transversely isotropic fluid (Ericksen 1960 Colloid Polym. Sci.173, 117-122 (doi:10.1007/bf01502416)) treats suspensions of non-motile particles as a continuum with an evolving preferred direction; this model describes fibrous materials as diverse as extracellular matrix, textile tufts and plant cell walls. Director-dependent effects are incorporated through a modified stress tensor with four viscosity-like parameters. By making fundamental connections with recent models for active suspensions, we propose a modification to Ericksen's model, mainly the inclusion of self-motility; this can be considered the simplest description of an oriented suspension including transversely isotropic effects. Motivated by the fact that transversely isotropic fluids exhibit modified flow stability, we conduct a linear stability analysis of two distinct cases, aligned and isotropic suspensions of elongated active particles. Novel aspects include the anisotropic rheology and translational diffusion. In general, anisotropic effects increase the instability of small perturbations, while translational diffusion stabilizes a range of wave-directions and, in some cases, a finite range of wavenumbers, thus emphasizing that both anisotropy and translational diffusion can have important effects in these systems.

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.

Understanding interactions between populations: Individual based modelling and quantification using pair correlation functions.

Understanding the underlying mechanisms that produce the huge variety of swarming and aggregation patterns in animals and cells is fundamental in ecology, developmental biology, and regenerative medicine, to name but a few examples. Depending upon the nature of the interactions between individuals (cells or animals), a variety of different large-scale spatial patterns can be observed in their distribution; examples include cell aggregates, stripes of different coloured skin cells, etc. For the case where all individuals are of the same type (i.e., all interactions are alike), a considerable literature already exists on how the collective organisation depends on the inter-individual interactions. Here, we focus on the less studied case where there are two different types of individuals present. Whilst a number of continuum models of this scenario exist, it can be difficult to compare these models to experimental data, since real cells and animals are discrete. In order to overcome this problem, we develop an agent-based model to simulate some archetypal mechanisms involving attraction and repulsion. However, with this approach (as with experiments), each realisation of the model is different, due to stochastic effects. In order to make useful comparisons between simulations and experimental data, we need to identify the robust features of the spatial distributions of the two species which persist over many realisations of the model (for example, the size of aggregates, degree of segregation or intermixing of the two species). In some cases, it is possible to do this by simple visual inspection. In others, the features of the pattern are not so clear to the unaided eye. In this paper, we introduce a pair correlation function (PCF), which allows us to analyse multi-species spatial distributions quantitatively. We show how the differing strengths of inter-individual attraction and repulsion between species give rise to different spatial patterns, and how the PCF can be used to quantify these differences, even when it might be impossible to recognise them visually.

Pattern formation in multiphase models of chemotactic cell aggregation.

We develop a continuum model for the aggregation of cells cultured in a nutrient-rich medium in a culture well. We consider a 2D geometry, representing a vertical slice through the culture well, and assume that the cell layer depth is small compared with the typical lengthscale of the culture well. We adopt a continuum mechanics approach, treating the cells and culture medium as a two-phase mixture. Specifically, the cells and culture medium are treated as fluids. Additionally, the cell phase can generate forces in response to environmental cues, which include the concentration of a chemoattractant that is produced by the cells within the culture medium. The model leads to a system of coupled nonlinear partial differential equations for the volume fraction and velocity of the cell phase, the culture medium pressure and the chemoattractant concentration, which must be solved subject to appropriate boundary and initial conditions. To gain insight into the system, we consider two model reductions, appropriate when the cell layer depth is thin compared to the typical length scale of the culture well: a (simple) 1D and a (more involved) thin-film extensional flow reduction. By investigating the resulting systems of equations analytically and numerically, we identify conditions under which small amplitude perturbations to a homogeneous steady state (corresponding to a spatially uniform cell distribution) can lead to a spatially varying steady state (pattern formation). Our analysis reveals that the simpler 1D reduction has the same qualitative features as the thin-film extensional flow reduction in the linear and weakly nonlinear regimes, motivating the use of the simpler 1D modelling approach when a qualitative understanding of the system is required. However, the thin-film extensional flow reduction may be more appropriate when detailed quantitative agreement between modelling predictions and experimental data is desired. Furthermore, full numerical simulations of the two model reductions in regions of parameter space when the system is not close to marginal stability reveal significant differences in the evolution of the volume fraction and velocity of the cell phase, and chemoattractant concentration.

Identifying the necrotic zone boundary in tumour spheroids with pair-correlation functions.

Automatic identification of the necrotic zone boundary is important in the assessment of treatments on in vitro tumour spheroids. This has been difficult especially when the difference in cell density between the necrotic and viable zones of a tumour spheroid is small. To help overcome this problem, we developed novel one-dimensional pair-correlation functions (PCFs) to provide quantitative estimates of the radial distance of the necrotic zone boundary from the centre of a tumour spheroid. We validate our approach on synthetic tumour spheroids in which the position of the necrotic zone boundary is known a priori It is then applied to nine real tumour spheroids imaged with light sheet-based fluorescence microscopy. PCF estimates of the necrotic zone boundary are compared with those of a human expert and an existing standard computational method.

A structure-function analysis of the left ventricle.

This study presents a structure-function analysis of the mammalian left ventricle and examines the performance of the cardiac capillary network, mitochondria, and myofibrils at rest and during simulated heavy exercise. Left ventricular external mechanical work rate was calculated from cardiac output and systemic mean arterial blood pressure in resting sheep (Ovis aries; n = 4) and goats (Capra hircus; n = 4) under mild sedation, followed by perfusion-fixation of the left ventricle and quantification of the cardiac capillary-tissue geometry and cardiomyocyte ultrastructure. The investigation was then extended to heavy exercise by increasing cardiac work according to published hemodynamics of sheep and goats performing sustained treadmill exercise. Left ventricular work rate averaged 0.017 W/cm3 of tissue at rest and was estimated to increase to ∼0.060 W/cm3 during heavy exercise. According to an oxygen transport model we applied to the left ventricular tissue, we predicted that oxygen consumption increases from 195 nmol O2·s-1·cm-3 of tissue at rest to ∼600 nmol O2·s-1·cm-3 during heavy exercise, which is within 90% of the oxygen demand rate and consistent with work remaining predominantly aerobic. Mitochondria represent 21-22% of cardiomyocyte volume and consume oxygen at a rate of 1,150 nmol O2·s-1·cm-3 of mitochondria at rest and ∼3,600 nmol O2·s-1·cm-3 during heavy exercise, which is within 80% of maximum in vitro rates and consistent with mitochondria operating near their functional limits. Myofibrils represent 65-66% of cardiomyocyte volume, and according to a Laplacian model of the left ventricular chamber, generate peak fiber tensions in the range of 50 to 70 kPa at rest and during heavy exercise, which is less than maximum tension of isolated cardiac tissue (120-140 kPa) and is explained by an apparent reserve capacity for tension development built into the left ventricle.

An investigation of the influence of extracellular matrix anisotropy and cell-matrix interactions on tissue architecture.

Mechanical interactions between cells and the fibrous extracellular matrix (ECM) in which they reside play a key role in tissue development. Mechanical cues from the environment (such as stress, strain and fibre orientation) regulate a range of cell behaviours, including proliferation, differentiation and motility. In turn, the ECM structure is affected by cells exerting forces on the matrix which result in deformation and fibre realignment. In this paper we develop a mathematical model to investigate this mechanical feedback between cells and the ECM. We consider a three-phase mixture of collagen, culture medium and cells, and formulate a system of partial differential equations which represents conservation of mass and momentum for each phase. This modelling framework takes into account the anisotropic mechanical properties of the collagen gel arising from its fibrous microstructure. We also propose a cell-collagen interaction force which depends upon fibre orientation and collagen density. We use a combination of numerical and analytical techniques to study the influence of cell-ECM interactions on pattern formation in tissues. Our results illustrate the wide range of structures which may be formed, and how those that emerge depend upon the importance of cell-ECM interactions.

On the derivation of approximations to cellular automata models and the assumption of independence.

Cellular automata are discrete agent-based models, generally used in cell-based applications. There is much interest in obtaining continuum models that describe the mean behaviour of the agents in these models. Previously, continuum models have been derived for agents undergoing motility and proliferation processes, however, these models only hold under restricted conditions. In order to narrow down the reason for these restrictions, we explore three possible sources of error in deriving the model. These sources are the choice of limiting arguments, the use of a discrete-time model as opposed to a continuous-time model and the assumption of independence between the state of sites. We present a rigorous analysis in order to gain a greater understanding of the significance of these three issues. By finding a limiting regime that accurately approximates the conservation equation for the cellular automata, we are able to conclude that the inaccuracy between our approximation and the cellular automata is completely based on the assumption of independence.

Distinguishing between mechanisms of cell aggregation using pair-correlation functions.

Many cell types form clumps or aggregates when cultured in vitro through a variety of mechanisms including rapid cell proliferation, chemotaxis, or direct cell-to-cell contact. In this paper we develop an agent-based model to explore the formation of aggregates in cultures where cells are initially distributed uniformly, at random, on a two-dimensional substrate. Our model includes unbiased random cell motion, together with two mechanisms which can produce cell aggregates: (i) rapid cell proliferation and (ii) a biased cell motility mechanism where cells can sense other cells within a finite range, and will tend to move towards areas with higher numbers of cells. We then introduce a pair-correlation function which allows us to quantify aspects of the spatial patterns produced by our agent-based model. In particular, these pair-correlation functions are able to detect differences between domains populated uniformly at random (i.e. at the exclusion complete spatial randomness (ECSR) state) and those where the proliferation and biased motion rules have been employed - even when such differences are not obvious to the naked eye. The pair-correlation function can also detect the emergence of a characteristic inter-aggregate distance which occurs when the biased motion mechanism is dominant, and is not observed when cell proliferation is the main mechanism of aggregate formation. This suggests that applying the pair-correlation function to experimental images of cell aggregates may provide information about the mechanism associated with observed aggregates. As a proof of concept, we perform such analysis for images of cancer cell aggregates, which are known to be associated with rapid proliferation. The results of our analysis are consistent with the predictions of the proliferation-based simulations, which supports the potential usefulness of pair correlation functions for providing insight into the mechanisms of aggregate formation.

Non-local models for the formation of hepatocyte-stellate cell aggregates.

Liver cell aggregates may be grown in vitro by co-culturing hepatocytes with stellate cells. This method results in more rapid aggregation than hepatocyte-only culture, and appears to enhance cell viability and the expression of markers of liver-specific functions. We consider the early stages of aggregate formation, and develop a new mathematical model to investigate two alternative hypotheses (based on evidence in the experimental literature) for the role of stellate cells in promoting aggregate formation. Under Hypothesis 1, each population produces a chemical signal which affects the other, and enhanced aggregation is due to chemotaxis. Hypothesis 2 asserts that the interaction between the two cell types is by direct physical contact: the stellates extend long cellular processes which pull the hepatocytes into the aggregates. Under both hypotheses, hepatocytes are attracted to a chemical they themselves produce, and the cells can experience repulsive forces due to overcrowding. We formulate non-local (integro-partial differential) equations to describe the densities of cells, which are coupled to reaction-diffusion equations for the chemical concentrations. The behaviour of the model under each hypothesis is studied using a combination of linear stability analysis and numerical simulations. Our results show how the initial rate of aggregation depends upon the cell seeding ratio, and how the distribution of cells within aggregates depends on the relative strengths of attraction and repulsion between the cell types. Guided by our results, we suggest experiments which could be performed to distinguish between the two hypotheses.

A mathematical model of liver cell aggregation in vitro.

The behavior of mammalian cells within three-dimensional structures is an area of intense biological research and underpins the efforts of tissue engineers to regenerate human tissues for clinical applications. In the particular case of hepatocytes (liver cells), the formation of spheroidal multicellular aggregates has been shown to improve cell viability and functionality compared to traditional monolayer culture techniques. We propose a simple mathematical model for the early stages of this aggregation process, when cell clusters form on the surface of the extracellular matrix (ECM) layer on which they are seeded. We focus on interactions between the cells and the viscoelastic ECM substrate. Governing equations for the cells, culture medium, and ECM are derived using the principles of mass and momentum balance. The model is then reduced to a system of four partial differential equations, which are investigated analytically and numerically. The model predicts that provided cells are seeded at a suitable density, aggregates with clearly defined boundaries and a spatially uniform cell density on the interior will form. While the mechanical properties of the ECM do not appear to have a significant effect, strong cell-ECM interactions can inhibit, or possibly prevent, the formation of aggregates. The paper concludes with a discussion of our key findings and suggestions for future work.