Mathematical models of cell colonization of uniformly growing domains

During the development of vertebrate embryos, cell migrations occur on an underlying tissue domain in response to some factor, such as nutrient. Over the time scale of days in which this cell migration occurs, the underlying tissue is itself growing. Consequently cell migration and colonization is strongly affected by the tissue domain growth. Numerical solutions for a mathematical model of chemotactic migration with no domain growth can lead to travelling waves of cells with constant velocity; the addition of domain growth can lead to travelling waves with nonconstant velocity. These observations suggest a mathematical approximation to the full system equations, allowing the method of characteristics to be applied to a simplified chemotactic migration model. The evolution of the leading front of the migrating cell wave is analysed. Linear, exponential and logistic uniform domain growths are considered. Successful colonization of a growing domain depends on the competition between cell migration velocity and the velocity and form of the domain growth, as well as the initial penetration distance of the cells. In some instances the cells will never successfully colonize the growing domain. These models provide an insight into cell migration during embryonic growth, and its dependence upon the form and timing of the domain growth.     

Landscape geometry and travelling waves in the larch budmoth

Travelling waves in cyclic populations refer to temporal shifts in peak densities moving across space in a wave‐like fashion. The epicentre hypothesis states that peak densities begin in specific geographic foci and then spread into adjoining areas. Travelling waves have been confirmed in a number of population systems, begging questions about their causes. Herein we apply a newly developed statistical technique, wavelet phase analysis, to historical data to document that the travelling waves in larch budmoth (LBM) outbreaks arise from two epicentres, both located in areas with high concentrations of favourable habitat. We propose that the spatial arrangement of the landscape mosaic is responsible for initiating the travelling waves. We use a tri‐trophic model of LBM dynamics to demonstrate that landscape heterogeneity (specifically gradients in density of favourable habitat) alone, is capable of inducing waves from epicentres. Our study provides unique evidence of how landscape features can mould travelling waves.     

Travelling Waves in Hybrid Chemotaxis Models

Hybrid models of chemotaxis combine agent-based models of cells with partial differential equation models of extracellular chemical signals. In this paper, travelling wave properties of hybrid models of bacterial chemotaxis are investigated. Bacteria are modelled using an agent-based (individual-based) approach with internal dynamics describing signal transduction. In addition to the chemotactic behaviour of the bacteria, the individual-based model also includes cell proliferation and death. Cells consume the extracellular nutrient field (chemoattractant), which is modelled using a partial differential equation. Mesoscopic and macroscopic equations representing the behaviour of the hybrid model are derived and the existence of travelling wave solutions for these models is established. It is shown that cell proliferation is necessary for the existence of non-transient (stationary) travelling waves in hybrid models. Additionally, a numerical comparison between the wave speeds of the continuum models and the hybrid models shows good agreement in the case of weak chemotaxis and qualitative agreement for the strong chemotaxis case. In the case of slow cell adaptation, we detect oscillating behaviour of the wave, which cannot be explained by mean-field approximations.     

A Simple Mathematical Model Inspired by the Purkinje Cells: From Delayed Travelling Waves to Fractional Diffusion

Recently, several experiments have demonstrated the existence of fractional diffusion in the neuronal transmission occurring in the Purkinje cells, whose malfunctioning is known to be related to the lack of voluntary coordination and the appearance of tremors. Also, a classical mathematical feature is that (fractional) parabolic equations possess smoothing effects, in contrast with the case of hyperbolic equations, which typically exhibit shocks and discontinuities. In this paper, we show how a simple toy-model of a highly ramified structure, somehow inspired by that of the Purkinje cells, may produce a fractional diffusion via the superposition of travelling waves that solve a hyperbolic equation. This could suggest that the high ramification of the Purkinje cells might have provided an evolutionary advantage of “smoothing” the transmission of signals and avoiding shock propagations (at the price of slowing a bit such transmission). Although an experimental confirmation of the possibility of such evolutionary advantage goes well beyond the goals of this paper, we think that it is intriguing, as a mathematical counterpart, to consider the time fractional diffusion as arising from the superposition of delayed travelling waves in highly ramified transmission media. The case of a travelling concave parabola with sufficiently small curvature is explicitly computed. The new link that we propose between time fractional diffusion and hyperbolic equation also provides a novelty with respect to the usual paradigm relating time fractional diffusion with parabolic equations in the limit. This paper is written in such a way as to be of interest to both biologists and mathematician alike. In order to accomplish this aim, both complete explanations of the objects considered and detailed lists of references are provided.     

A diffusive SI model with Allee effect and application to FIV

A minimal reaction-diffusion model for the spatiotemporal spread of an infectious disease is considered. The model is motivated by the Feline Immunodeficiency Virus (FIV) which causes AIDS in cat populations. Because the infected period is long compared with the lifespan, the model incorporates the host population growth. Two different types are considered: logistic growth and growth with a strong Allee effect. In the model with logistic growth, the introduced disease propagates in form of a travelling infection wave with a constant asymptotic rate of spread. In the model with Allee effect the spatiotemporal dynamics are more complicated and the disease has considerable impact on the host population spread. Most importantly, there are waves of extinction, which arise when the disease is introduced in the wake of the invading host population. These waves of extinction destabilize locally stable endemic coexistence states. Moreover, spatially restricted epidemics are possible as well as travelling infection pulses that correspond either to fatal epidemics with succeeding host population extinction or to epidemics with recovery of the host population. Generally, the Allee effect induces minimum viable population sizes and critical spatial lengths of the initial distribution. The local stability analysis yields bistability and the phenomenon of transient epidemics within the regime of disease-induced extinction. Sustained oscillations do not exist.     

Travelling wave patterns in a model of the spinal pattern generator using spiking neurons

The aim of this study is to produce travelling waves in a planar net of artificial spiking neurons. Provided that the parameters of the waves – frequency, wavelength and orientation – can be sufficiently controlled, such a network can serve as a model of the spinal pattern generator for swimming and terrestrial quadruped locomotion. A previous implementation using non-spiking, sigmoid neurons lacked the physiological plausibility that can only be attained using more realistic spiking neurons. Simulations were conducted using three types of spiking neuronal models. First, leaky integrate-and-fire neurons were used. Second, we introduced a phenomenological bursting neuron. And third, a canonical model neuron was implemented which could reproduce the full dynamics of the Hodgkin–Huxley neuron. The conditions necessary to produce appropriate travelling waves corresponded largely to the known anatomy and physiology of the spinal cord. Especially important features for the generation of travelling waves were the topology of the local connections – so-called off-centre connectivity – the availability of dynamic synapses and, to some extent, the availability of bursting cell types. The latter were necessary to produce stable waves at the low frequencies observed in quadruped locomotion. In general, the phenomenon of travelling waves was very robust and largely independent of the network parameters and emulated cell types.     

Gene-culture waves of advance

Major genetic and cultural changes may have been coupled during hominid evolution. Since hominids have had a wide geographical distribution for about one million years, any mutant gene or cultural innovation that became established had to spread from its origin. A pair of nonlinear diffusion equations is derived which models the propagation of a mutant gene and a cultural innovation. Both are assumed to originate in the same locality along a linear habitat. The mutant gene and its allele are semidominant, and the two cultural choices are transmitted according to what I call the logistic attraction-repulsion model. The genes influence cultural choice, and the two interact to determine fitness. Of particular interest is the case in which mutant gene and cultural innovation are mutually dependent, neither being able to spread without the other. Each equation of the pair is similar in form to Fisher's equation, with a linear function of the other dependent variable replacing the constant coefficient in the reaction term. The partial differential equations are solved numerically to obtain the asymptotic speeds. Their form also suggests an heuristic argument which has proved useful, but I have been unable to obtain any analytic results. The waves of the system are shown to be of two types, synchronous and asynchronous. When genes and culture are mutually dependent, synchronous travelling waves can exist. However, their existence is dependent on initial conditions, and the speed of propagation is slow.     

The geometry and motion of reaction-diffusion waves on closed two-dimensional manifolds

Chemical or biological systems modelled by reaction diffusion (R.D.) equations which support simple one-dimensional travelling waves (oscillatory or otherwise) may be expected to produce intricate two or three-dimensional spatial patterns, either stationary or subject to certain motion. Such structures have been observed experimentally. Asymptotic considerations applied to a general class of such systems lead to fundamental restrictions on the existence and geometrical form of possible structures. As a consequence of the geometrical setting, it is a straightforward matter to consider the propagation of waves on closed two-dimensional manifolds. We derive a fundamental equation for R.D. wave propagation on surfaces and discuss its significance. We consider the existence and propagation of rotationally symmetric and double spiral waves on the sphere and on the torus. On leave of absence from: Department of Mathematics, Glasgow College of Technology, Cowcaddens Road, Glasgow G4 0BA, Scotland, UK     

Early development and quorum sensing in bacterial biofilms

 We develop mathematical models to examine the formation, growth and quorum sensing activity of bacterial biofilms. The growth aspects of the model are based on the assumption of a continuum of bacterial cells whose growth generates movement, within the developing biofilm, described by a velocity field. A model proposed in Ward et al. (2001) to describe quorum sensing, a process by which bacteria monitor their own population density by the use of quorum sensing molecules (QSMs), is coupled with the growth model. The resulting system of nonlinear partial differential equations is solved numerically, revealing results which are qualitatively consistent with experimental ones. Analytical solutions derived by assuming uniform initial conditions demonstrate that, for large time, a biofilm grows algebraically with time; criteria for linear growth of the biofilm biomass, consistent with experimental data, are established. The analysis reveals, for a biologically realistic limit, the existence of a bifurcation between non-active and active quorum sensing in the biofilm. The model also predicts that travelling waves of quorum sensing behaviour can occur within a certain time frame; while the travelling wave analysis reveals a range of possible travelling wave speeds, numerical solutions suggest that the minimum wave speed, determined by linearisation, is realised for a wide class of initial conditions. Received: 10 February 2002 / Revised version: 29 October 2002 / Published online: 19 March 2003 Key words or phrases: Bacterial biofilm – Quorum sensing – Mathematical modelling – Numerical solution – Asymptotic analysis – Travelling wave analysis     

Homoclinic solutions in mechanical systems with small dissipation. Application to DNA dynamics

 This paper analyzes the problem of persistence of homoclinic solutions to perturbed systems of second order ODE's. These systems arise from PDE's, when considering solutions in the form of travelling waves. It is shown that homoclinic solutions persist in the presence of dissipation. Dissipation can be balanced by nonautonomous terms of compact support, which are controlled by a single parameter. This result is applied to prove the existence of torsional pulse-like travelling waves propagating along a nonelastic DNA molecule. In this case the energy is added to system by advancing the RNA polymerase. Received: 17 August 2000 / Revised version: 26 October 2001 / Published online: 14 March 2002     

Travelling wave solutions of diffusive Lotka-Volterra equations

We establish the existence of travelling wave solutions for two reaction diffusion systems based on the Lotka-Volterra model for predator and prey interactions. For simplicity, we consider only 1 space dimension. The waves are of transition front type, analogous to the travelling wave solutions discussed by Fisher and Kolmogorov et al. for a scalar reaction diffusion equation. The waves discussed here are not necessarily monotone. For any speed c there is a travelling wave solution of transition front type. For one of the systems discussed here, there is a distinguished speed c^* dividing the waves into two types, waves of speed c < c^* being one type, waves of speed c ? c^* being of the other type. We present numerical evidence that for this system the wave of speed c^* is stable, and that c^* is an asymptotic speed of propagation in some sense. For the other system, waves of all speeds are in some sense stable. The proof of existence uses a shooting argument and a Lyapunov function. We also discuss some possible biological implications of the existence of these waves.     

Travelling waves of a delayed SIR epidemic model with nonlinear incidence rate and spatial diffusion

This paper is concerned with the existence of travelling waves to a SIR epidemic model with nonlinear incidence rate, spatial diffusion and time delay. By analyzing the corresponding characteristic equations, the local stability of a disease-free steady state and an endemic steady state to this system under homogeneous Neumann boundary conditions is discussed. By using the cross iteration method and the Schauder's fixed point theorem, we reduce the existence of travelling waves to the existence of a pair of upper-lower solutions. By constructing a pair of upper-lower solutions, we derive the existence of a travelling wave connecting the disease-free steady state and the endemic steady state. Numerical simulations are carried out to illustrate the main results.     

From periodic travelling waves to travelling fronts in the spike-diffuse-spike model of dendritic waves

In the vertebrate brain excitatory synaptic contacts typically occur on the tiny evaginations of neuron dendritic surface known as dendritic spines. There is clear evidence that spine heads are endowed with voltage-dependent excitable channels and that action potentials invade spines. Computational models are being increasingly used to gain insight into the functional significance of a spine with an excitable membrane. The spike-diffuse-spike (SDS) model is one such model that admits to a relatively straightforward mathematical analysis. In this paper we demonstrate that not only can the SDS model support solitary travelling pulses, already observed numerically in more detailed biophysical models, but that it has periodic travelling wave solutions. The exact mathematical treatment of periodic travelling waves in the SDS model is used, within a kinematic framework, to predict the existence of connections between two periodic spike trains of different interspike interval. The associated wave front in the sequence of interspike intervals travels with a constant velocity without degradation of shape, and might therefore be used for the robust encoding of information.     

Comparison of definitions of the lag phase and the exponential phase in bacterial growth

M.H. ZWIETERING, F.M. ROMBOUTS AND K. VAN 'T RIET. 1992. Different definitions of the lag time and of the duration of the exponential phase can be used to calculate these quantities from growth models. The conventional definitions were compared with newly proposed definitions. It appeared to be possible to derive values for the lag time and the duration of the exponential phase from the growth models, and differences between the various definitions could be quantified. All the different values can be calculated from the growth parameters μ _m , and a. Therefore, it appeared to be unnecessary to use complicated mathematical equations: simple equations were adequate. For the Gompertz model the conventional definition of the lag time did not differ appreciably from the newly proposed definition. The end-point of the exponential phase and thus the duration of the exponential phase differed considerably for the two definitions. For the logistic model the two definitions lead to considerable differences for all quantities. It is recommended that the conventional definition is used for calculating the lag time. For the duration of the exponential phase it is recommended that the new definition is used. The value can be calculated, however, directly from the conventional growth parameters.     

Waves of crown-of-thorns starfish outbreaks —where do they come from?

To explain the southward propagating waves of crown-of-thorns outbreaks on the Great Barrier Reef, it has been proposed that the northern region acts as a seed area for these waves. In this paper we study whether the highly variable current pattern of this region facilitates the start of local outbreaks, making the northern region the seed area or pace-maker of travelling waves of outbreaks. To this end we construct an artificial reef system, which resembles the Great Barrier Reef by having a northern area with a random current pattern, and a central and southern region with a (southward) biased current pattern. Travelling waves running from north to south are observed in the system both with and without a deviant current pattern in the northern area. This demonstrates that travelling waves are an emergent property of the system. The result indicates that the observed waves of crown-of-thorns outbreaks on the Great Barrier Reef do not necessarily imply that there is a seed area in the north.     

Generation of periodic waves by landscape features in cyclic predator-prey systems

The vast majority of models for spatial dynamics of natural populations assume a homogeneous physical environment. However, in practice, dispersing organisms may encounter landscape features that significantly inhibit their movement. We use mathematical modelling to investigate the effect of such landscape features on cyclic predator-prey populations. We show that when appropriate boundary conditions are applied at the edge of the obstacle, a pattern of periodic travelling waves develops, moving out and away from the obstacle. Depending on the assumptions of the model, these waves can take the form of roughly circular 'target patterns' or spirals. This is, to our knowledge, a new mechanism for periodic-wave generation in ecological systems and our results suggest that it may apply quite generally not only to cyclic predator-prey interactions, but also to populations that oscillate for other reasons. In particular, we suggest that it may provide an explanation for the observed pattern of travelling waves in the densities of field voles (Microtus agrestis) in Kielder Forest (Scotland-England border) and of red grouse (Lagopus lagopus scoticus) on Kerloch Moor (northeast Scotland), which in both cases move orthogonally to any large-scale obstacles to movement. Moreover, given that such obstacles to movement are the rule rather than the exception in real-world environments, our results suggest that complex spatio-temporal patterns such as periodic travelling waves are likely to be much more common in the natural world than has previously been assumed.     

Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications

A mathematical model of the spatio-temporal dynamics of a two host, two parasitoid system is presented. There is a coupling of the four species through parasitism of both hosts by one of the parasitoids. The model comprises a system of four reaction-diffusion equations. The underlying system of ordinary differential equations, modelling the host-parasitoid population dynamics, has a unique positive steady state and is shown to be capable of undergoing Hopf bifurcations, leading to limit cycle kinetics which give rise to oscillatory temporal dynamics. The stability of the positive steady state has a fundamental impact on the spatio-temporal dynamics: stable travelling waves of parasitoid invasion exhibit increasingly irregular periodic travelling wave behaviour when key parameter values are increased beyond their Hopf bifurcation point. These irregular periodic travelling waves give rise to heterogeneous spatio-temporal patterns of host and parasitoid abundance. The generation of heterogeneous patterns has ecological implications and the concepts of temporary host refuge and niche formation are considered.     

Periodic travelling waves in cyclic predator–prey systems

Predation is an established cause of cycling in prey species. Here, the ability of predation to explain periodic travelling waves in prey populations, which have recently been found in a number of spatiotemporal field studies, is examined. The nature of periodic waves in these systems, and the way in which they can be generated by the invasion of predators into a prey population is discussed. A theoretical calculation that predicts, as a function of two parameter ratios, whether such an invasion will lead to a stable periodic travelling wave that would be observed in practice is presented ‐ the alternative outcome is spatiotemporal chaos. The calculation also predicts quantitative details of the periodic waves, such as speed and amplitude. The results give new insights into the types of predator‐prey systems in which one would expect to see periodic travelling waves following an invasion by predators.     

Spatially mismatched trophic dynamics: cyclically outbreaking geometrids and their larval parasitoids

For trophic interactions to generate population cycles and complex spatio-temporal patterns, like travelling waves, the spatial dynamics must be matched across trophic levels. Here, we propose a spatial methodological approach for detecting such spatial match–mismatch and apply it to geometrid moths and their larval parasitoids in northern Norway, where outbreak cycles and travelling waves occur. We found clear evidence of spatial mismatch, suggesting that the spatially patterned moth cycles in this system are probably ruled by trophic interactions involving other agents than larval parasitoids.