Nonlinear dynamics and pattern bifurcations in a model for vegetation stripes in semi-arid environments

In many semi-arid environments, vegetation is self-organised into spatial patterns. The most striking examples of this are on gentle slopes, where striped patterns are typical, running parallel to the contours. Previously, Klausmeier [1999. Regular and irregular patterns in semiarid vegetation. Science 284, 1826-1828.] has proposed a model for vegetation stripes based on competition for water. Here, we present a detailed study of the patterned solutions in the full nonlinear model, using numerical bifurcation analysis of both the pattern odes and the model pdes. We show that patterns exist for a wide range of rainfall levels, and in particular for much lower rainfall than have been considered by previous authors. Moreover, we show that for many rainfall levels, patterns with a variety of different wavelengths are stable, with mode selection dependent on initial conditions. This raises the possibility of hysteresis, and in numerical solutions of the model we show that pattern selection depends on rainfall history in a relatively simple way.     

On a simplified model for pattern formation in honey bee colonies

We present a simplified version of a previously presented model (Camazine et al. (1990)) that generates the characteristic pattern of honey, pollen and brood which develops on combs in honey bee colonies. We demonstrate that the formation of a band of pollen surrounding the brood area is dependent on the assumed form of the honey and pollen removal terms, and that a significant pollen band arises as the parameter controlling the rate of pollen input passes through a bifurcation value. The persistence of the pollen band after a temporary increase in pollen input can be predicted from the model. We also determine conditions on the parameters which ensure the accumulation of honey in the periphery and demonstrate that, although there is an important qualitative difference between the simplified and complete models, an analysis of the simplified version helps us understand many biological aspects of the more complex complete model. Corresponding author     

Applications of a model for scale-invariant pattern formation in developing systems

A fundamental problem in developmental biology concerns the proportioning of the developing tissue of a morphallactic system into different cell types in a way that is independent of the overall size of the tissue. The two main models for positional information in pattern formation, the source-sink models and the Turing reaction-diffusion models, have shortcomings that limit their applicability. In a previous paper, we described a model that can produce perfectly scale-invariant spatial patterns and analyzed some of its mathematical properties. In the present paper, we demonstrate some of the shortcomings of the standard reaction-diffusion models and discuss the applicability of our model to developmental systems.     

Travelling waves and pattern formation in a model for fungal development

 Under a variety of conditions, the hyphal density within the expanding outer edge of growing fungal mycelia can be spatially heterogeneous or nearly uniform. We conduct an analysis of a system of reaction-diffusion equations used to model the growth of fungal mycelia and the subsequent development of macroscopic patterns produced by differing hyphal and hence biomass densities. Both local and global results are obtained using analytical and numerical techniques. The emphasis is on qualitative results, including the effects of changes in parameter values on the structure of the solution set. Received 22 November 1995; received in revised form 17 May 1996     

Travelling waves in a tissue interaction model for skin pattern formation

Tissue interaction plays a major role in many morphogenetic processes, particularly those associated with skin organ primordia. We examine travelling wave solutions in a tissue interaction model for skin pattern formation which is firmly based on the known biology. From a phase space analysis we conjecture the existence of travelling waves with specific wave speeds. Subsequently, analytical approximations to the wave profiles are derived using perturbation methods. We then show numerically that such travelling wave solutions do exist and that they are in good agreement with our analytical results. Finally, the biological implications of our analysis are discussed.     

A new model for biological pattern formation

Various non-equilibrium growth models have been used to explore the development of morphology in biological systems. Here we review a class of biological growth models which exhibit fractal structures and discuss the relationship of these models to a variety of other phenomena.     

Numerical investigation of spatial pattern in a vegetation model with feedback function

The vegetative cover in semi-arid lands typically occurs as patches of individual species more or less separated from one another by bare ground. Klausmeier [1999. Regular and irregular patterns in semiarid vegetation. Science 284 (5421), 1826-1828] reported that the vegetation striped patterns can grow lying along the contours of gentle slopes. He has proposed a model of vegetation stripes based on competition for water. In this paper, our main aim is to study the positive feedback effects between the water and biomass on the vegetation spatial pattern formation within a nonsaturated soil, which arises from the suction of water by the roots and processes of water resource redistribution. According to the dispersion relation formula, we discuss the changes of the wavelength, wave speed, as well as the conditions of the spatial pattern formation. Our numerical results show that trees are more sensitive than grasses to the positive feedback function to format the spatial heterogenous pattern, and the stronger positive feedback increases the parameters region where vegetation bands occur, which indicates that the positive feedback raises the possibility of shift from green to desert states in semi-arid areas for the long term. Our numerical results also show that the positive feedback can increase the migration velocity of the vegetation stripes.     

Travelling lipid domains in a dynamic model for protein-induced pattern formation in biomembranes

Cell membranes are composed of a mixture of lipids. Many biological processes require the formation of spatial domains in the lipid distribution of the plasma membrane. We have developed a mathematical model that describes the dynamic spatial distribution of acidic lipids in response to the presence of GMC proteins and regulating enzymes. The model encompasses diffusion of lipids and GMC proteins, electrostatic attraction between acidic lipids and GMC proteins as well as the kinetics of membrane attachment/detachment of GMC proteins. If the lipid-protein interaction is strong enough, phase separation occurs in the membrane as a result of free energy minimization and protein/lipid domains are formed. The picture is changed if a constant activity of enzymes is included into the model. We chose the myristoyl-electrostatic switch as a regulatory module. It consists of a protein kinase C that phosphorylates and removes the GMC proteins from the membrane and a phosphatase that dephosphorylates the proteins and enables them to rebind to the membrane. For sufficiently high enzymatic activity, the phase separation is replaced by travelling domains of acidic lipids and proteins. The latter active process is typical for nonequilibrium systems. It allows for a faster restructuring and polarization of the membrane since it acts on a larger length scale than the passive phase separation. The travelling domains can be pinned by spatial gradients in the activity; thus the membrane is able to detect spatial clues and can adapt its polarity dynamically to changes in the environment.     

A cellular oscillator model for periodic pattern formation

In this paper, we present a model for pattern formation in developing organisms that is based on cellular oscillators (CO). An oscillatory process within cells serves as a developmental clock whose period is tightly regulated by cell autonomous or non-autonomous mechanisms. A spatial pattern is generated as a result of an initial temporal ordering of the cell oscillators freezing into spatial order as the clocks slow down and stop at different times or phases in their cycles. We apply a CO model to vertebrate somitogenesis and show that we can reproduce the dynamics of periodic gene expression patterns observed in the pre-somitic mesoderm. We also show how varying somite lengths can be generated with the CO model. We then discuss the model in view of experimental evidence and its relevance to other instances of biological pattern formation, showing its versatility as a pattern generator.     

Taxis-driven pattern formation in a predator-prey model with group defense

We consider a reaction-diffusion(-taxis) predator-prey system with group defense in the prey. Taxis-driven instability can occur if the group defense influences the taxis rate (Wang et al., 2017). We elaborate that this mechanism is indeed possible but biologically unlikely to be responsible for pattern formation in such a system. Conversely, we show that patterns in excitable media such as spatiotemporal Sierpinski gasket patterns occur in the reaction-diffusion model as well as in the reaction-diffusion-taxis model. If group defense leads to a dome-shaped functional response, these patterns can have a rescue effect on the predator population in an invasion scenario. Preytaxis with prey repulsion at high prey densities can intensify this mechanism leading to taxis-induced persistence. In particular, taxis can increase parameter regimes of successful invasions and decrease minimum introduction areas necessary for a successful invasion. Last, we consider the mean period of the irregular oscillations. As a result of the underlying mechanism of the patterns, this period is two orders of magnitude smaller than the period in the nonspatial system. Counter-intuitively, faster-moving predators lead to lower oscillation periods and eventually to extinction of the predator population. The study does not only provide valuable insights on theoretical spatially explicit predator-prey models with group defense but also comparisons of ecological data with model simulations.     

A model of morphogenetic pattern formation

A model for the morphogenetic movement of surfaces composed of cellular monolayers is proposed. The cells are presumed joined at their lateral surfaces. An otherwise unspecified substance called a "morphogen" is introduced which is the agent of change in the individual cell (or cell-like region). The distribution of these cellular deformations define a surface (the middle surface, through the middle of the cell heights) via equations given for the Gauss and Mean curvatures of the surface defined at each point. The Gauss curvature as a function of the morphogen level determines the metric of the surface "g(u, v)" in conformal co-ordinates u, v. A unique equation for the morphogen distribution over the survace is presented which has the property of size invariance, that is, the model "regulates" without need of further arguments. The two resulting coupled equations for the metric and the morphogen, eqns (4) and (2), both non-linear equations, are to be solved self-consistently, once the individual cell deformation as a function of morphogen is given. The surface geometry determines the morphogen distribution, and the morphogen distribution in turn affects the surface geometry. Extension of the model to two or more morphogens is straightforward, and the key property of "regulation" or size invariance of the model is retained. Numerical integration of the two coupled equations is carried out in the case of axial symmetry, and the results presented by the case that individual cells deform by changing the ratio of their apical to basal areas, as well as their heights. Gastrulation in small regulating holoblastic eggs (e.g. starfish, sea urchin and amphioxus) is discussed in light of the present model.     

Cell sorting by differential cell motility: a model for pattern formation in Dictyostelium

In the slug stage of the cellular slime mold Dictyostelium discoideum, prespore cells and four types of prestalk cells show a well-defined spatial distribution in a migrating slug. We have developed a continuous mathematical model for the distribution pattern of these cell types based on the balance of force in individual cells. In the model, cell types are assumed to have different properties in cell motility, i.e. different motive force, the rate of resistance against cell movement, and diffusion coefficient. Analysis of the stationary solution of the model shows that combination of these parameters and slug speed determines the three-dimensional shape of a slug and cell distribution pattern within it. Based on experimental data of slug motive force and velocity measurements, appropriate sets of parameters were chosen so that the cell-type distribution at stationary state matches the distribution in real slugs. With these parameters, we performed numerical calculation of the model in two-dimensional space using a moving particle method. The results reproduced many of the basic features of slug morphogenesis, i.e. cell sorting, translocation of the prestalk region, elongation of the slug, and its steady migration.     

Searching for a model for use in vegetation analysis

Summary Indirect gradient analysis methods require an explicit vegetation model which must be based on direct gradient analysis studies. Various vegetation models are reviewed. Field evidence for the models is discussed. Experimental studies of species response to environmental gradients are reviewed and discussed. Three types of gradient are recognized as important for development of models: indirect environmental gradients where the environmental factor has no direct physiological influence on plant growth e.g. elevation; direct environmental gradients where the factor has a direct physiological effect on growth but is not an essential resource, e.g. pH; resource gradients where the factor is an essential resource for plant growth. The behaviour of the ecological carrying capacity and the role of competition along such gradients are shown to be important for developing vegetation models.     

A model for budding in hydra: pattern formation in concentric rings

Current models of pattern formation in Hydra propose head-and foot-specific morphogens to control the development of the body ends and along the body length axis. In addition, these morphogens are proposed to control a cellular parameter (positional value, source density) which changes gradually along the axis. This gradient determines the tissue polarity and the regional capacity to form a head and a foot, respectively, in transplantation experiments. The current models are very successful in explaining regeneration and transplantation experiments. However, some results obtained render problems, in particular budding, the asexual way of reproduction is not understood. Here an alternative model is presented to overcome these problems. A primary system of interactions controls the positional values. At certain positional values secondary systems become active which initiate the local formation of e.g. mouth, tentacles, and basal disc. (i) A system of autocatalysis and lateral inhibition is suggested to exist as proposed by Gierer and Meinhardt (Kybernetik 12 (1972) 30). (ii) The activator is neither a head nor a foot activator but rather causes an increase of the positional value. (iii) On the other hand, a generation of the activator leads to its loss from cells and therewith to a (local) decrease of the positional value. (iv) An inhibitor is proposed to exist which antagonizes an increase of the positional value. External conditions like the gradient of positional values in the surroundings and interactions with other sites of morphogen production decide whether at a certain site of activator generation the positional value will increase (head formation), decrease (foot formation) or increase in the centre and decrease in the periphery thereby forming concentric rings (bud formation). Computer-simulation experiments show basic features of budding, regeneration and transplantation.     

A model for rouleaux pattern formation of red blood cells

Human red blood cells (RBCs) in a solution form rouleaux patterns under various conditions. The degree of rouleaux formation depends on, for example, the concentration and molecular weight of added large molecules. We present a two-dimensional discrete cellular space model in which an RBC is represented by a rectangle and differential adhesion is assumed among the longer (a-site), the shorter (b-site) sides of the rectangle and the solvent. The total sum of the adhesion energy is assumed to guide the step-by-step change of the model cell configuration and also define absolutely stable patterns. We compare the set of absolutely stable patterns and cell aggregate patterns for both actual and computer-simulated cases to obtain the basic validity of our framework. Then we proceed to assess the effects of added high polymers to the adhesion parameters. We first note that under suitable conditions, decrease in a-site-solvent affinity is necessary to have complex patterns rather than increase of a-a affinity. The hypothesis that addition of high polymers reduce the a-site-solvent affinity is concomitant with a newly proposed osmotic stress theory. The parameter fitting results for the experimental phase change curves can also be interpreted as supporting more the new theory than existing traditional explanations.