Cellular retinoic acid-binding proteins are essential for hindbrain patterning and signal robustness in zebrafish.

The vitamin A derivative retinoic acid (RA) is a morphogen that patterns the anterior-posterior axis of the vertebrate hindbrain. Cellular retinoic acid-binding proteins (Crabps) transport RA within cells to both its nuclear receptors (RARs) and degrading enzymes (Cyp26s). However, mice lacking Crabps are viable, suggesting that Crabp functions are redundant with those of other fatty acid-binding proteins. Here we show that Crabps in zebrafish are essential for posterior patterning of the hindbrain and that they provide a key feedback mechanism that makes signaling robust as they are able to compensate for changes in RA production. Of the four zebrafish Crabps, Crabp2a is uniquely RA inducible and depletion or overexpression of Crabp2a makes embryos hypersensitive to exogenous RA. Computational models confirm that Crabp2a improves robustness within a narrow concentration range that optimizes a 'robustness index', integrating spatial information along the RA morphogen gradient. Exploration of signaling parameters in our models suggests that the ability of Crabp2a to transport RA to Cyp26 enzymes for degradation is a major factor in promoting robustness. These results demonstrate a previously unrecognized requirement for Crabps in RA signaling and hindbrain development, as well as a novel mechanism for stabilizing morphogen gradients despite genetic or environmental fluctuations in morphogen availability.

Red queen dynamics in specific predator-prey systems.

The dynamics of a predator-prey system are studied, with a comparison of discrete and continuous strategy spaces. For a [Formula: see text] system, the average strategies used in the discrete and continuous case are shown to be the same. It is further shown that the inclusion of constant prey switching in the discrete case can have a stabilising effect and reduce the number of available predator types through extinction.

Noise drives sharpening of gene expression boundaries in the zebrafish hindbrain.

Morphogens provide positional information for spatial patterns of gene expression during development. However, stochastic effects such as local fluctuations in morphogen concentration and noise in signal transduction make it difficult for cells to respond to their positions accurately enough to generate sharp boundaries between gene expression domains. During development of rhombomeres in the zebrafish hindbrain, the morphogen retinoic acid (RA) induces expression of hoxb1a in rhombomere 4 (r4) and krox20 in r3 and r5. Fluorescent in situ hybridization reveals rough edges around these gene expression domains, in which cells co-express hoxb1a and krox20 on either side of the boundary, and these sharpen within a few hours. Computational analysis of spatial stochastic models shows, surprisingly, that noise in hoxb1a/krox20 expression actually promotes sharpening of boundaries between adjacent segments. In particular, fluctuations in RA initially induce a rough boundary that requires noise in hoxb1a/krox20 expression to sharpen. This finding suggests a novel noise attenuation mechanism that relies on intracellular noise to induce switching and coordinate cellular decisions during developmental patterning.

Cell proliferation and oxygen diffusion in a vascularising scaffold.

The supply of oxygen to proliferating cells within a scaffold is a key factor for the successful building of new tissue in soft tissue engineering applications. A recent in vivo model, where an arteriovenous loop is placed in a scaffold, allows a vascularising network to form within a scaffold, establishing an oxygen source within, rather than external, to the scaffold. A one-dimensional model of oxygen concentration, cell proliferation and cell migration inside such a vascularising scaffold is developed and investigated. In addition, a vascularisation model is presented, which supports a vascularisation front which moves at a constant speed. The effects of vascular growth, homogenous and heterogenous seeding, diffusion of cells and critical hypoxic oxygen concentration are considered. For homogenous seeding, a relationship between the speed of the vascular front and a parameter defining the rate of oxygen diffusion relative to the rate of oxygen consumption determines whether a hypoxic region exists at some time. In particular, an estimate of the length of time that a fixed point in the scaffold will remain under hypoxic conditions is determined. For heterogenous seeding, a Fisher-like travelling wave of cells is established behind the vascular front. These findings provide a fundamental understanding of the important interplay between the parameters and allows for a theoretical assessment of a seeding strategy in a vascularising scaffold.

Travelling waves of attached and detached cells in a wound-healing cell migration assay.

During a wound-healing cell migration assay experiment, cells are observed to detach and undergo mitosis before reattaching as a pair of cells on the substrate. During experiments with mice 3T3 fibroblasts, cell detachment can be confined to the wavefront region or it can occur over the whole wave profile. A multi-species continuum approach to modelling a wound-healing assay is taken to account for this phenomenon. The first cell population is composed of attached motile cells, while the second population is composed of detached immotile cells undergoing mitosis and returning to the migrating population after successful division. The first model describes cell division occurring only in the wavefront region, while a second model describes cell division over the whole of the scrape wound. The first model reverts to the Fisher equation when the reattachment rate relative to the detachment rate is large, while the second case does not revert to the Fisher equation in any limit. The models yield travelling wave solutions. The minimum wave speed is slower than the minimum Fisher wave speed and its dependence on the cell detachment and reattachment rate parameters is analysed. Approximate travelling wave profiles of the two cell populations are determined asymptotically under certain parameter regimes.

Multi-scale modeling of a wound-healing cell migration assay.

A continuum model and a discrete model are developed to capture the population-scale and cell-scale behavior in a wound-healing cell migration assay created from a scrape wound in a confluent cell monolayer. During wound closure, the cell population forms a sustained traveling wave, with close contact between cells behind the wavefront. Cells exhibit contact inhibition of migration and contact-limited proliferation. The continuum model includes the two dominant mechanisms and characteristics of cell migration and proliferation, using a cell diffusivity function that decreases with cell density and a logistic proliferative growth term. The discrete model arises naturally from the continuum model. Individual cells are simulated as continuous-time random walkers with nearest-neighbor transitions, together with a birth/death process. The migration and proliferation parameters are determined by analysing individual mice 3T3 fibroblast cell trajectories obtained during the development of a confluent cell monolayer and in a wound healing assay. The population-scale model successfully predicts the shape and speed of the traveling wave, while the discrete model is also successful in capturing the contact inhibition of migration effects.

T cell development in the thymus: from periodic seeding to constant output.

T cell development occurs in the thymus throughout life. Recent experimental findings show that the seeding of the thymus by multi-potent stem cells from the bone marrow is periodic rather than continuous, as previously assumed. However it is well known that the output rate of cells from the thymus is relatively constant. A quantitative model is used to verify the current hypotheses regarding T cell development in the steady state mouse thymus. The results show that the thymus could be at a periodic steady state with out-of-phase thymocyte populations. Experiments to examine possible periodic fluctuations in the thymus are proposed and methods for further analysis are outlined.

Modelling directional guidance and motility regulation in cell migration.

Although cell migration is an essential process in development, how cells reach their final destination is not well understood. Secreted molecules are known to have a migratory effect, but it remains unclear whether such molecules act as directional guidance cues or as motility regulators. There is potential to use signalling molecules in new medical therapies, so it is important to identify the exact role these molecules play. This paper focuses on distinguishing between inhibitory and repulsive effects produced by signalling molecules, based on recent experiments examining the effect of Slit, a secreted protein, on the migration of neurons from the brain. The primary role of Slit, whether it is an inhibitor or repellent of neurons, is in dispute. We present population-level continuum models and recast these in terms of transition probabilities governing individual cells. Various cell-sensing strategies are considered within this framework. The models are applied to the neuronal migration experiments. To resolve the particular role of Slit, simulations of the models characterising different cell-sensing strategies are compared at the population and individual cell level, providing two complementary perspectives on the system. Difficulties and limitations in deducing cell migration rules from time-lapse imaging are discussed.