Building a Morphogen Gradient without Diffusion in a Growing Tissue

In many developmental systems, spatial pattern arises from morphogen gradients, which provide positional information for cells to determine their fate. Typically, diffusion is thought to be the mechanism responsible for building a morphogen gradient. An alternative mechanism is investigated here. Using mathematical modeling, we demonstrate how a non-diffusive morphogen concentration gradient can develop in axially growing tissue systems, where growth is due to cell proliferation only. Two distinct cases are considered: in the first, all cell proliferation occurs in a localized zone where active transcription of a morphogen-producing gene occurs, and in the second, cell proliferation is uniformly distributed throughout the tissue, occurring in both the active transcription zone and beyond. A cell containing morphogen mRNA produces the morphogen protein, hence any gradient in mRNA transcripts translates into a corresponding morphogen protein gradient. Proliferation-driven growth gives rise to both advection (the transport term) and dilution (a reaction term). These two key mechanisms determine the resultant mRNA transcript distribution. Using the full range of uniform initial conditions, we show that advection and dilution due to cell proliferation are, in general, sufficient for morphogen gradient formation for both types of axially growing systems. In particular, mRNA transcript degradation is not necessary for gradient formation; it is only necessary with localized proliferation for one special value of the initial concentration. Furthermore, the morphogen concentration decreases with distance away from the transcription zone, except in the case of localized proliferation with the initial concentration sufficiently large, when the concentration can either increase with distance from the transcription zone or sustain a local minimum. In both localized and uniformly distributed proliferation, in order for a concentration gradient to form across the whole domain, transcription must occur in a zone equal to the initial domain size; otherwise, it will only form across part of the tissue

Is the skin an excitable medium? Pattern formation in erythema gyratum repens

Erythema gyratum repens (EGR) is a rare, inflammatory dermatosis of unknown aetiology. The morphology of the eruption is striking and displays rapidly evolving circinate and gyrate bands of erythematous and scaly skin. Although the aetiology of the pattern is unknown, it has previously been noted that the eruption shares morphologic features with the patterns of spatio-temporal chemical concentration profiles observed in the Belusov-Zhabotinski (BZ) reaction. Yet this morphologic correspondence has not been investigated further. Here we apply a simple non-linear reaction -diffusion model, previously used to describe the BZ reaction, as a template for pattern formation in EGR, and show how the mechanism may provide a biochemical basis for many of the dynamic and morphologic features of the rash. These results are supported by the results of a cellular automaton simulation approximating the dynamics of oscillatory chemical systems-the Hodgepodge machine-where the spatio-temporal patterns developed show astonishing similarities to the morphology of EGR

Random walks in nonuniform environments with local dynamic interactions

We consider a class of lattice random walk models in which the random walker is initially confined to a finite connected set of allowed sites but has the opportunity to enlarge this set by colliding with its boundaries, each such collision having a given probability of breaking through. The model is motivated by an analogy to cell motility in tissue, where motile cells have the ability to remodel extracellular matrix, but is presented here as a generic model for stochastic erosion. For the one-dimensional case, we report some exact analytic results, some mean-field type analytic approximate results and simulations. We compute exactly the mean and variance of the time taken to enlarge the interval from a single site to a given size. The problem of determining the statistics of the interval length and the walker's position at a given time is more difficult and we report several interesting observations from simulations. Our simulations include the case in which the initial interval length is random and the case in which the initial state of the lattice is a random mixture of allowed and forbidden sites, with the walker placed at random on an allowed site. To illustrate the extension of these ideas to higher-dimensional systems, we consider the erosion of the simple cubic lattice commencing from a single site and report simulations of measures of cluster size and shape and the mean-square displacement of the walker

Growth patterns in Onychophora (velvet worms): lack of a localised posterior proliferation zone

<p>Abstract</p> <p>Background</p> <p>During embryonic development of segmented animals, body segments are thought to arise from the so-called "posterior growth zone" and the occurrence of this "zone" has been used to support the homology of segmentation between arthropods, annelids, and vertebrates. However, the term "posterior growth zone" is used ambiguously in the literature, mostly referring to a region of increased proliferation at the posterior end of the embryo. To determine whether such a localised posterior proliferation zone is an ancestral feature of Panarthropoda (Onychophora + Tardigrada + Arthropoda), we examined cell division patterns in embryos of Onychophora.</p> <p>Results</p> <p>Using in vivo incorporation of the DNA replication marker BrdU (5-bromo-2'-deoxyuridine) and anti-phospho-histone H3 immunolabelling, we found that a localised posterior region of proliferating cells does not occur at any developmental stage in onychophoran embryos. This contrasts with a localised pattern of cell divisions at the posterior end of annelid embryos, which we used as a positive control. Based on our data, we present a mathematical model, which challenges the paradigm that a localised posterior proliferation zone is necessary for segment patterning in short germ developing arthropods.</p> <p>Conclusions</p> <p>Our findings suggest that a posterior proliferation zone was absent in the last common ancestor of Onychophora and Arthropoda. By comparing our data from Onychophora with those from annelids, arthropods, and chordates, we suggest that the occurrence of a "posterior growth zone" currently cannot be used to support the homology of segmentation between these three animal groups.</p

Generalized index for spatial data sets as a measure of complete spatial randomness

Spatial data sets, generated from a wide range of physical systems can be analyzed by counting the number of objects in a set of bins. Previous work has been limited to equal-sized bins, which are inappropriate for some domains (e.g., circular). We consider a nonequal size bin configuration whereby overlapping or nonoverlapping bins cover the domain. A generalized index, defined in terms of a variance between bin counts, is developed to indicate whether or not a spatial data set, generated from exclusion or nonexclusion processes, is at the complete spatial randomness (CSR) state. Limiting values of the index are determined. Using examples, we investigate trends in the generalized index as a function of density and compare the results with those using equal size bins. The smallest bin size must be much larger than the mean size of the objects. We can determine whether a spatial data set is at the CSR state or not by comparing the values of a generalized index for different bin configurations—the values will be approximately the same if the data is at the CSR state, while the values will differ if the data set is not at the CSR state. In general, the generalized index is lower than the limiting value of the index, since objects do not have access to the entire region due to blocking by other objects. These methods are applied to two applications: (i) spatial data sets generated from a cellular automata model of cell aggregation in the enteric nervous system and (ii) a known plant data distribution.Emily J. Hackett-Jones, Kale J. Davies, Benjamin J. Binder, and Kerry A. Landma

A Multicellular Model of Intestinal Crypt Buckling and Fission

Crypt fission is an in vivo tissue deformation process that is involved in both intestinal homeostasis and colorectal tumourigenesis. Despite its importance, the mechanics underlying crypt fission are currently poorly understood. Recent experimental development of organoids, organ-like buds cultured from crypt stem cells in vitro, has shown promise in shedding light on crypt fission. Drawing inspiration from observations of organoid growth and fission in vivo, we develop a computational model of a deformable epithelial tissue layer. Results from in silico experiments show the stiffness of cells and the proportions of cell subpopulations affect the nature of deformation in the epithelial layer. In particular, we find that increasing the proportion of stiffer cells in the layer increases the likelihood of crypt fission occurring. This is in agreement with and helps explain recent experimental work

Computational Modelling of Tissue-Engineered Cartilage Constructs

Cartilage is a fundamental tissue to ensure proper motion between bones and damping of mechanical loads. This tissue often suffers damage and has limited healing capacity due to its avascularity. In order to replace surgery and replacement of joints by metal implants, tissue engineered cartilage is seen as an attractive alternative. These tissues are obtained by seeding chondrocytes or mesenchymal stem cells in scaffolds and are given certain stimuli to improve establishment of mechanical properties similar to the native cartilage. However, tissues with ideal mechanical properties were not obtained yet. Computational models of tissue engineered cartilage growth and remodelling are invaluable to interpret and predict the effects of experimental designs. The current model contribution in the field will be presented in this chapter, with a focus on the response to mechanical stimulation, and the development of fully coupled modelling approaches incorporating simultaneously solute transport and uptake, cell growth, production of extracellular matrix and remodelling of mechanical properties.publishe

A continuum model for the development of tissue-engineered cartilage around a chondrocyte

The limited ability of cartilage to repair when damaged has led to the investigation of tissue engineering as a method for reconstructing cartilage. We propose a continuum multispecies model for the development of cartilage around a single chondrocyte. As in healthy cartilage, the model predicts a balance between synthesis, transport, binding and decay of matrix components. Two mechanisms are investigated for the transport of soluble matrix components: diffusion and advection, caused by displacement of the scaffold medium. Numerical results indicate that a parameter defined by the ratio of the flux of soluble components out of the chondrocyte and its diffusive flux determines which of these mechanisms is dominant. We investigate the diffusion-dominated and advection-dominated limiting cases using perturbation analysis. Using parameter values from the literature, our modelling results suggest that both diffusion and advection are significant mechanisms in developing cartilage. Moreover, in this parameter regime, results are particularly sensitive to parameter values. These two observations could explain differences observed experimentally between various scaffold media. Modelling results are also used to predict the minimum chondrocyte seeding density required to produce functional cartilage. <br/