Extending the Mathematical Palette for Developmental Pattern Formation: Piebaldism.

Here, we present a theoretical investigation with potential insights on developmental mechanisms. Three biological factors, consisting of two diffusing factors and a cell-autonomous immobile transcription factor are combined with different feedback mechanisms. This results in four different situations or fur patterns. Two of them reproduce classical Turing patterns: (1) regularly spaced spots, (2) labyrinth patterns or straight lines with an initial slope in the activation of the transcription factor. The third situation does not lead to patterns, but results in different homogeneous color tones. Finally, the fourth one sheds new light on the possible mechanisms leading to the formation of piebald patterns exemplified by the random patterns on the fur of some cows' strains and Dalmatian dogs. Piebaldism is usually manifested as white areas of fur, hair, or skin due to the absence of pigment-producing cells in those regions. The distribution of the white and colored zones does not reflect the classical Turing patterns. We demonstrate that these piebald patterns are of transient nature, developing from random initial conditions and relying on a system's bistability. We show numerically that the presence of a cell-autonomous factor not only expands the range of reaction diffusion parameters in which a pattern may arise, but also extends the pattern-forming abilities of the reaction-diffusion equations.

The feasibility of equilibria in large ecosystems: A primary but neglected concept in the complexity-stability debate.

The consensus that complexity begets stability in ecosystems was challenged in the seventies, a result recently extended to ecologically-inspired networks. The approaches assume the existence of a feasible equilibrium, i.e. with positive abundances. However, this key assumption has not been tested. We provide analytical results complemented by simulations which show that equilibrium feasibility vanishes in species rich systems. This result leaves us in the uncomfortable situation in which the existence of a feasible equilibrium assumed in local stability criteria is far from granted. We extend our analyses by changing interaction structure and intensity, and find that feasibility and stability is warranted irrespective of species richness with weak interactions. Interestingly, we find that the dynamical behaviour of ecologically inspired architectures is very different and richer than that of unstructured systems. Our results suggest that a general understanding of ecosystem dynamics requires focusing on the interplay between interaction strength and network architecture.

The Effect of Gap Junctional Coupling on the Spatiotemporal Patterns of Ca2+ Signals and the Harmonization of Ca2+-Related Cellular Responses.

Calcium ions (Ca2+) are important mediators of a great variety of cellular activities e.g. in response to an agonist activation of a receptor. The magnitude of a cellular response is often encoded by frequency modulation of Ca2+ oscillations and correlated with the stimulation intensity. The stimulation intensity highly depends on the sensitivity of a cell to a certain agonist. In some cases, it is essential that neighboring cells produce a similar and synchronized response to an agonist despite their different sensitivity. In order to decipher the presumed function of Ca2+ waves spreading among connecting cells, a mathematical model was developed. This model allows to numerically modifying the connectivity probability between neighboring cells, the permeability of gap junctions and the individual sensitivity of cells to an agonist. Here, we show numerically that strong gap junctional coupling between neighbors ensures an equilibrated response to agonist stimulation via formation of Ca2+ phase waves, i.e. a less sensitive neighbor will produce the same or similar Ca2+ signal as its highly sensitive neighbor. The most sensitive cells within an ensemble are the wave initiator cells. The Ca2+ wave in the cytoplasm is driven by a sensitization wave front in the endoplasmic reticulum. The wave velocity is proportional to the cellular sensitivity and to the strength of the coupling. The waves can form different patterns including circular rings and spirals. The observed pattern depends on the strength of noise, gap junctional permeability and the connectivity probability between neighboring cells. Our simulations reveal that one highly sensitive region gradually takes the lead within the entire noisy system by generating directed circular phase waves originating from this region.