How does cellular contact affect differentiation mediated pattern formation?

Linking discrete and continuous models of cell birth and migration

Self-organization of individuals within large collectives occurs throughout biology. Mathematical models can help elucidate the individual-level mechanisms behind these dynamics, but analytical tractability often comes at the cost of biological intuition. Discrete models provide straightforward interpretations by tracking each individual yet can be computationally expensive. Alternatively, continuous models supply a large-scale perspective by representing the 'effective' dynamics of infinite agents, but their results are often difficult to translate into experimentally relevant insights. We address this challenge by quantitatively linking spatio-temporal dynamics of continuous models and individual-based data in settings with biologically realistic, time-varying cell numbers. Specifically, we introduce and fit scaling parameters in continuous models to account for discrepancies that can arise from low cell numbers and localized interactions. We illustrate our approach on an example motivated by zebrafish-skin pattern formation, in which we create a continuous framework describing the movement and proliferation of a single cell population by upscaling rules from a discrete model. Our resulting continuous models accurately depict ensemble average agent-based solutions when migration or proliferation act alone. Interestingly, the same parameters are not optimal when both processes act simultaneously, highlighting a rich difference in how combining migration and proliferation affects discrete and continuous dynamics.

A quantitative modelling approach to zebrafish pigment pattern formation

Pattern formation is a key aspect of development. Adult zebrafish exhibit a striking striped pattern generated through the self-organisation of three different chromatophores. Numerous investigations have revealed a multitude of individual cell-cell interactions important for this self-organisation, but it has remained unclear whether these known biological rules were sufficient to explain pattern formation. To test this, we present an individual-based mathematical model incorporating all the important cell-types and known interactions. The model qualitatively and quantitatively reproduces wild type and mutant pigment pattern development. We use it to resolve a number of outstanding biological uncertainties, including the roles of domain growth and the initial iridophore stripe, and to generate hypotheses about the functions of leopard. We conclude that our rule-set is sufficient to recapitulate wild-type and mutant patterns. Our work now leads the way for further in silico exploration of the developmental and evolutionary implications of this pigment patterning system.

Modeling Stripe Formation on Growing Zebrafish Tailfins

As zebrafish develop, black and gold stripes form across their skin due to the interactions of brightly colored pigment cells. These characteristic patterns emerge on the growing fish body, as well as on the anal and caudal fins. While wild-type stripes form parallel to a horizontal marker on the body, patterns on the tailfin gradually extend distally outward. Interestingly, several mutations lead to altered body patterns without affecting fin stripes. Through an exploratory modeling approach, our goal is to help better understand these differences between body and fin patterns. By adapting a prior agent-based model of cell interactions on the fish body, we present an in silico study of stripe development on tailfins. Our main result is a demonstration that two cell types can produce stripes on the caudal fin. We highlight several ways that bone rays, growth, and the body-fin interface may be involved in patterning, and we raise questions for future work related to pattern robustness.

Iridophores as a source of robustness in zebrafish stripes and variability in Danio patterns

Zebrafish (Danio rerio) feature black and yellow stripes, while related Danios display different patterns. All these patterns form due to the interactions of pigment cells, which self-organize on the fish skin. Until recently, research focused on two cell types (melanophores and xanthophores), but newer work has uncovered the leading role of a third type, iridophores: by carefully orchestrated transitions in form, iridophores instruct the other cells, but little is known about what drives their form changes. Here we address this question from a mathematical perspective: we develop a model (based on known interactions between the original two cell types) that allows us to assess potential iridophore behavior. We identify a set of mechanisms governing iridophore form that is consistent across a range of empirical data. Our model also suggests that the complex cues iridophores receive may act as a key source of redundancy, enabling both robust patterning and variability within Danio.

Modelling stripe formation in zebrafish: an agent-based approach

Zebrafish have distinctive black stripes and yellow interstripes that form owing to the interaction of different pigment cells. We present a two-population agent-based model for the development and regeneration of these stripes and interstripes informed by recent experimental results. Our model describes stripe pattern formation, laser ablation and mutations. We find that fish growth shortens the necessary scale for long-range interactions and that iridophores, a third type of pigment cell, help align stripes and interstripes.

Aggregation patterns from nonlocal interactions: Discrete stochastic and continuum modeling

Conservation equations governed by a nonlocal interaction potential generate aggregates from an initial uniform distribution of particles. We address the evolution and formation of these aggregating steady states when the interaction potential has both attractive and repulsive singularities. Currently, no existence theory for such potentials is available. We develop and compare two complementary solution methods, a continuous pseudoinverse method and a discrete stochastic lattice approach, and formally show a connection between the two. Interesting aggregation patterns involving multiple peaks for a simple doubly singular attractive-repulsive potential are determined. For a swarming Morse potential, characteristic slow-fast dynamics in the scaled inverse energy is observed in the evolution to steady state in both the continuous and discrete approaches. The discrete approach is found to be remarkably robust to modifications in movement rules, related to the potential function. The comparable evolution dynamics and steady states of the discrete model with the continuum model suggest that the discrete stochastic approach is a promising way of probing aggregation patterns arising from two- and three-dimensional nonlocal interaction conservation equations.