Flow-distributed oscillation, flow-velocity modulation, and resonance

Morphogen-directed cell fate boundaries: slow passage through bifurcation and the role of folded saddles

One of the fundamental mechanisms in embryogenesis is the process by which cells differentiate and create tissues and structures important for functioning as a multicellular organism. Morphogenesis involves diffusive process of chemical signalling involving morphogens that pre-pattern the tissue. These morphogens influence cell fate through a highly nonlinear process of transcriptional signalling. In this paper, we consider this multiscale process in an idealised model for a growing domain. We focus on intracellular processes that lead to robust differentiation into two cell lineages through interaction of a single morphogen species with a cell fate variable that undergoes a bifurcation from monostability to bistability. In particular, we investigate conditions that result in successful and robust pattern formation into two well-separated domains, as well as conditions where this fails and produces a pinned boundary wave where only one part of the domain grows. We show that successful and unsuccessful patterning scenarios can be characterised in terms of presence or absence of a folded saddle singularity for a system with two slow variables and one fast variable; this models the interaction of slow morphogen diffusion, slow parameter drift through bifurcation and fast transcription dynamics. We illustrate how this approach can successfully model acquisition of three cell fates to produce three-domain "French flag" patterning, as well as for a more realistic model of the cell fate dynamics in terms of two mutually inhibiting transcription factors.

Flow-induced transitions in bistable systems

We studied transitions between spatiotemporal patterns that can be induced in a spatially extended nonlinear chemical system by a unidirectional flow in combination with constant inflow concentrations. Three different scenarios were investigated. (i) Under conditions where the system exhibited two stable fixed points, the propagation direction of trigger fronts could be reversed, so that domains of the less stable fixed point invaded the system. (ii) For bistability between a stable fixed point and a limit cycle we observed that above a critical flow velocity, the unstable focus at the center of the limit cycle could be stabilized. Increasing the flow speed further, a regime of damped flow-distributed oscillations was found and, depending on the boundary values at the inflow, finally the stable fixed point dominated. Similarly, also in the case of spatiotemporal chaos (iii), the unstable steady state could be stabilized and was replaced by the stable fixed point with increasing flow velocity. We finally outline a linear stability analysis that can explain part of our findings.

Harmonic resonant excitation of flow-distributed oscillation waves and Turing patterns driven at a growing boundary

We perform numerical studies of a reaction-diffusion system that is both Turing and Hopf unstable, and that grows by addition at a moving boundary (which is equivalent by a Galilean transformation to a reaction-diffusion-advection system with a fixed boundary and a uniform flow). We model the conditions of a recent set of experiments which used a temporally varying illumination in the boundary region to control the formation of patterns in the bulk of the photosensitive medium. The frequency of the illumination variations can select patterns from among the competing instabilities of the medium. In the usual case, the waves that are excited have frequencies (as measured at a constant distance from the upstream boundary) matching the driving frequency. In contrast to the usual case, we find that both Turing patterns and flow-distributed oscillation waves can be excited by forcing at subharmonic multiples of the wave frequencies. The final waves (with frequencies at integer multiples of the driving frequency) are formed by a process in which transient wave fronts break up and reconnect. We find ratios of response to driving frequency as high as 10.

Selection of flow-distributed oscillation and Turing patterns by boundary forcing in a linearly growing, oscillating medium

We studied the response of a linearly growing domain of the oscillatory chemical chlorine dioxide-iodide-malonic acid (CDIMA) medium to periodic forcing at its growth boundary. The medium is Hopf-, as well as Turing-unstable and the system is convectively unstable. The results confirm numerical predictions that two distinct modes of pattern can be excited by controlling the driving frequency at the boundary, a flow-distributed-oscillation (FDO) mode of traveling waves at low values of the forcing frequency f , and a mode of stationary Turing patterns at high values of f . The wavelengths and phase velocities of the experimental patterns were compared quantitatively with results from dynamical simulations and with predictions from linear dispersion relations. The results for the FDO waves agreed well with these predictions, and obeyed the kinematic relations expected for phase waves with frequencies selected by the boundary driving frequency. Turing patterns were also generated within the predicted range of forcing frequencies, but these developed into two-dimensional structures which are not fully accounted for by the one-dimensional numerical and analytical models. The Turing patterns excited by boundary forcing persist when the forcing is removed, demonstrating the bistability of the unforced, constant size medium. Dynamical simulations at perturbation frequencies other than those of the experiments showed that in certain ranges of forcing frequency, FDO waves become unstable, breaking up into harmonic waves of different frequency and wavelength and phase velocity.