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

Developmental plasticity and variability in the formation of egg-spots, a pigmentation ornament in the cichlid Astatotilapia calliptera

Vertebrate pigmentation patterns are highly diverse, yet we have a limited understanding of how evolutionary changes to genetic, cellular, and developmental mechanisms generate variation. To address this, we examine the formation of a sexually-selected male ornament exhibiting inter- and intraspecific variation, the egg-spot pattern, consisting of circular yellow-orange markings on the male anal fins of haplochromine cichlid fishes. We focus on Astatotilapia calliptera, the ancestor-type species of the Malawi cichlid adaptive radiation of over 850 species. We identify a key role for iridophores in initializing egg-spot aggregations composed of iridophore-xanthophore associations. Despite adult sexual dimorphism, aggregations initially form in both males and females, with development only diverging between the sexes at later stages. Unexpectedly, we found that the timing of egg-spot initialization is plastic. The earlier individuals are socially isolated, the earlier the aggregations form, with iridophores being the cell type that responds to changes to the social environment. Furthermore, we observe apparent competitive interactions between adjacent egg-spot aggregations, which strongly suggests that egg-spot patterning results mostly from cell-autonomous cellular interactions. Together, these results demonstrate that A. calliptera egg-spot development is an exciting model for investigating pigment pattern formation at the cellular level in a system with developmental plasticity, sexual dimorphism, and intraspecific variation. As A. calliptera represents the ancestral bauplan for egg-spots, these findings provide a baseline for informed comparisons across the incredibly diverse Malawi cichlid radiation.

Soliton approximation in continuum models of leader-follower behavior

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.

The Unreasonable Effectiveness of Reaction Diffusion in Vertebrate Skin Color Patterning

In 1952, Alan Turing published the reaction-diffusion (RD) mathematical framework, laying the foundations of morphogenesis as a self-organized process emerging from physicochemical first principles. Regrettably, this approach has been widely doubted in the field of developmental biology. First, we summarize Turing's line of thoughts to alleviate the misconception that RD is an artificial mathematical construct. Second, we discuss why phenomenological RD models are particularly effective for understanding skin color patterning at the meso/macroscopic scales, without the need to parameterize the profusion of variables at lower scales. More specifically, we discuss how RD models (a) recapitulate the diversity of actual skin patterns, (b) capture the underlying dynamics of cellular interactions, (c) interact with tissue size and shape, (d) can lead to ordered sequential patterning, (e) generate cellular automaton dynamics in lizards and snakes, (f) predict actual patterns beyond their statistical features, and (g) are robust to model variations. Third, we discuss the utility of linear stability analysis and perform numerical simulations to demonstrate how deterministic RD emerges from the underlying chaotic microscopic agents.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Characterization and functional analysis of pax3 in body color transition of polychromatic Midas cichlids (Amphilophus citrinellus)

As the representative genetic and economic trait of ornamental fish, skin color has a strong impact on speciation and adaptation. However, the genetic basis of skin color pigmentation, differentiation and change is still not understood. The Midas cichlid fish with three typical body color transition stages of "black-gray‑gold" is an ideal model system for investigating the formation and change of fish body color. In this study, to investigate the regulatory role of the pair box 3 (pax3) gene in the early body color fading process of Midas cichlids, the complete cDNA sequence (3513 bp) of pax3 was successfully isolated from Midas cichlids (Amphilophus Citrinellus), and found to encode polypeptides of 491 amino acids. Expression patterns of the pax3 gene in tissues of Midas cichlids during different periods, including embryonic development and body color fading stages were detected by quantitative real-time PCR. The qRT-PCR analysis showed that pax3 was expressed in all tissues of adult fish, with a higher expression level in muscle and skin. The highest expression level in muscle tissue was significantly higher than that in other tissues (P < 0.05). During embryonic development, the expression tendency of pax3 was first increased and then decreased. In the three typical stages of early skin color fading from black to gold, pax3 expression in skin, caudal fin and scales all showed a downward trend. The expression level in the black stage was significantly higher than that in other stages (P < 0.05). Positive signal of pax3 protein was detected in the three typical skin color conversion stages, and the highest positive signal intensity was detected in the black stage, which was consistent with qRT-PCR results. After pax3 RNA interference, pax3 and the downstream genes mitf and tyr all decreased, while dct mRNA expression increased in the skin of fish. Western blotting also showed a decrease in pax3 protein concentration. Those results suggest that pax3 plays an important role in skin color formation, distribution and change in Midas cichlids through the melanogenesis pathway.

The present and future of Turing models in developmental biology

The Turing model (or reaction-diffusion model), first published in 1952, is a mathematical model that can account for autonomy in the morphogenesis of organisms. Although initially controversial, the model has gradually gained wider acceptance among experimental embryologists due to the accumulation of experimental data to support it. More recently, this model and others based on it have been used not only to explain biological phenomena conceptually but also as working hypotheses for molecular-level experiments and as internal components of more-complex 3D models. In this Spotlight, I will provide a personal perspective from an experimental biologist on some of the recent developments of the Turing model.

Zebrafish airinemes optimize their shape between ballistic and diffusive search

In addition to diffusive signals, cells in tissue also communicate via long, thin cellular protrusions, such as airinemes in zebrafish. Before establishing communication, cellular protrusions must find their target cell. Here we demonstrate that the shapes of airinemes in zebrafish are consistent with a finite persistent random walk model. The probability of contacting the target cell is maximized for a balance between ballistic search (straight) and diffusive search (highly curved, random). We find that the curvature of airinemes in zebrafish, extracted from live cell microscopy, is approximately the same value as the optimum in the simple persistent random walk model. We also explore the ability of the target cell to infer direction of the airineme's source, finding that there is a theoretical trade-off between search optimality and directional information. This provides a framework to characterize the shape, and performance objectives, of non-canonical cellular protrusions in general.

Studies of Turing pattern formation in zebrafish skin

Skin patterns are the first example of the existence of Turing patterns in living organisms. Extensive research on zebrafish, a model organism with stripes on its skin, has revealed the principles of pattern formation at the molecular and cellular levels. Surprisingly, although the networks of cell-cell interactions have been observed to satisfy the 'short-range activation and long-range inhibition' prerequisites for Turing pattern formation, numerous individual reactions were not envisioned based on the classical reaction-diffusion model. For example, in real skin, it is not an alteration in concentrations of chemicals, but autonomous migration and proliferation of pigment cells that establish patterns, and cell-cell interactions are mediated via direct contact through cell protrusions. Therefore, the classical reaction-diffusion mechanism cannot be used as it is for modelling skin pattern formation. Various studies are underway to adapt mathematical models to the experimental findings on research into skin patterns, and the purpose of this review is to organize and present them. These novel theoretical methods could be applied to autonomous pattern formation phenomena other than skin patterns. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Pigment pattern morphospace of Danio fishes: evolutionary diversification and mutational effects

Molecular and cellular mechanisms underlying variation in adult form remain largely unknown. Adult pigment patterns of fishes in the genus Danio, which includes zebrafish, Danio rerio, consist of horizontal stripes, vertical bars, spots and uniform patterns, and provide an outstanding opportunity to identify causes of species level variation in a neural crest derived trait. Understanding pigment pattern variation requires quantitative approaches to assess phenotypes, yet such methods have been mostly lacking for pigment patterns. We introduce metrics derived from information theory that describe patterns and pattern variation in Danio fishes. We find that these metrics used singly and in multivariate combinations are suitable for distinguishing general pattern types, and can reveal even subtle phenotypic differences attributable to mutations. Our study provides new tools for analyzing pigment pattern in Danio and potentially other groups, and sets the stage for future analyses of pattern morphospace and its mechanistic underpinnings.

Summary: A multidimensional morphospace for pigment patterns yields quantitative insights into the evolution and genetics of diverse pigment patterns across zebrafish and related species.

Genetics of white color and iridophoroma in "Lemon Frost" leopard geckos

The squamates (lizards and snakes) are close relatives of birds and mammals, with more than 10,000 described species that display extensive variation in a number of important biological traits, including coloration, venom production, and regeneration. Due to a lack of genomic tools, few genetic studies in squamates have been carried out. The leopard gecko, Eublepharis macularius, is a popular companion animal, and displays a variety of coloration patterns. We took advantage of a large breeding colony and used linkage analysis, synteny, and homozygosity mapping to investigate a spontaneous semi-dominant mutation, “Lemon Frost”, that produces white coloration and causes skin tumors (iridophoroma). We localized the mutation to a single locus which contains a strong candidate gene, SPINT1, a tumor suppressor implicated in human skin cutaneous melanoma (SKCM) and over-proliferation of epithelial cells in mice and zebrafish. Our work establishes the leopard gecko as a tractable genetic system and suggests that a tumor suppressor in melanocytes in humans can also suppress tumor development in iridophores in lizards.

The squamates (lizards and snakes) comprise a diverse group of reptiles, with more than 10,000 described species that display extensive variation in a number of important biological traits, including coloration. In this manuscript, we used quantitative genetics and genomics to map the mutation underlying white coloration in the Lemon Frost morph of the common leopard gecko, Eublepharis macularius. Lemon Frost geckos have increased white body coloration with brightened yellow and orange areas. This morph also displays a high incidence of iridophoroma, a tumor of white-colored cells. We obtained phenotype information and DNA samples from geckos in a large breeding colony and used genome sequencing and genetic linkage analysis to localize the Lemon Frost mutation to a single locus. This locus contains a strong candidate gene, SPINT1, a tumor suppressor implicated in human skin cutaneous melanoma. Together with other recent advances, our work brings reptiles into the modern genetics era.

Mathematical modeling of the interaction between yolk utilization and fish growth in zebrafish, Danio rerio

Optimal embryonic development plays a major role in the health of an individual beyond the developmental stage. Nutritional perturbation during development is associated with cardiovascular and metabolic disease later in life. With both nutritional uptake and overall growth being risk factors for eventual health, it is necessary to understand not only the behavior of the processes during development but also their interactions. In this study, we used differential equations, image analyses, curve fittings, parameter estimation and laboratory experiments to quantify the rate of yolk absorption and its effect on early development of a vertebrate model (Danio rerio). Findings from this study establish a nonlinear functional relationship between nutrient absorption and early fish growth. We found that the rate of change in fish length and yolk utilization is logistic, that is the yolk decays rapidly for a period of time before leveling out. An interesting finding from this study is that yolk utilization reaches its maximum at 84 h post-fertilization. We validated our mathematical models against experimental observations, making them powerful tools for replication and future simulations.

A complex genetic architecture in zebrafish relatives Danio quagga and D. kyathit underlies development of stripes and spots

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons.

Pigment patterns of fishes are diverse and function in a wide range of behaviors. Common pattern themes include stripes and spots, exemplified by the closely related minnows Danio quagga and D. kyathit, respectively. We show that these patterns arise late in development owing to alterations in the development and arrangements of pigment cells. In the closely related model organism zebrafish (D. rerio) single genes can switch the pattern from stripes to spots. Yet, we show that pattern differences between D. quagga and D. kyathit have a more complex genetic basis, depending on multiple genes and interactions between these genes. Our findings illustrate the importance of characterizing naturally occurring genetic variants, in addition to laboratory induced mutations, for a more complete understanding of pigment pattern development and evolution.

Reaction-diffusion in a growing 3D domain of skin scales generates a discrete cellular automaton

We previously showed that the adult ocellated lizard skin colour pattern is effectively generated by a stochastic cellular automaton (CA) of skin scales. We additionally suggested that the canonical continuous 2D reaction-diffusion (RD) process of colour pattern development is transformed into this discrete CA by reduced diffusion coefficients at the borders of scales (justified by the corresponding thinning of the skin). Here, we use RD numerical simulations in 3D on realistic lizard skin geometries and demonstrate that skin thickness variation on its own is sufficient to cause scale-by-scale coloration and CA dynamics during RD patterning. In addition, we show that this phenomenon is robust to RD model variation. Finally, using dimensionality-reduction approaches on large networks of skin scales, we show that animal growth affects the scale-colour flipping dynamics by causing a substantial decrease of the relative length scale of the labyrinthine colour pattern of the lizard skin.

Evolution of pigment cells and patterns: recent insights from teleost fishes

Skin pigment patterns of vertebrates are stunningly diverse, and nowhere more so than in teleost fishes. Several species, including relatives of zebrafish, recently evolved cichlid fishes of East Africa, clownfishes, deep sea fishes, and others are providing insights into pigment pattern evolution. This overview describes recent advances in understanding periodic patterns, like stripes and spots, the loss of patterns, and the role of cell-type diversification in generating pigmentation phenotypes. Advances in this area are being facilitated by the application of modern methods of gene editing, genomics, computational analysis, and other approaches to non-traditional model organisms having interesting pigmentary phenotypes. Several topics worthy of future attention are outlined as well.

Reproduction of patterns in melanocytic proliferations by agent-based simulation and geometric modeling

Spatio-temporal patterns of melanocytic proliferations observed in vivo are important for diagnosis but the mechanisms that produce them are poorly understood. Here we present an agent-based model for simulating the emergence of the main biologic patterns found in melanocytic proliferations. Our model portrays the extracellular matrix of the dermo-epidermal junction as a two-dimensional manifold and we simulate cellular migration in terms of geometric translations driven by adhesive, repulsive and random forces. Abstracted cellular functions and melanocyte-matrix interactions are modeled as stochastic events. For identification and validation we use visual renderings of simulated cell populations in a horizontal perspective that reproduce growth patterns observed in vivo by sequential dermatoscopy and corresponding vertical views that reproduce the arrangement of melanocytes observed in histopathologic sections. Our results show that a balanced interplay of proliferation and migration produces the typical reticular pattern of nevi, whereas the globular pattern involves additional cellular mechanisms. We further demonstrate that slight variations in the three basic cellular properties proliferation, migration, and adhesion are sufficient to produce a large variety of morphological appearances of nevi. We anticipate our model to be a starting point for the reproduction of more complex scenarios that will help to establish functional connections between abstracted microscopic behavior and macroscopic patterns in all types of melanocytic proliferations including melanoma.

Zebrafish Melanophores Suggest Novel Functions of Cell Chirality in Tissue Formation

Several types of cells show left–right asymmetric behavior, unidirectional rotation, or spiral movements. For example, neutrophil-like differentiated HL60 (dHL60) cells show leftward bias in response to chemoattractant. Neurons extend neurites, creating a clockwise spiral. Platelet cells shows unidirectional spiral arrangements of actin fibers. In the microfabricated culture environment, groups of C2C12 cells (mouse myoblast cell line) were autonomously aligned in a counter-clockwise spiral pattern, and isolated C2C12 cells showed unidirectional spiral pattern of the actin skeleton. This biased directionality suggested that these cells have inherent cell chirality. In addition to these cells, we recently found that melanophores of zebrafish also have an intrinsic cellular chirality that was shown by their counter-clockwise self-rotation. Although this cell chirality is obvious, the function of the cell chirality is still unclear. In this review, we compare the cell chirality of melanophores of zebrafish with other cell chirality and consider the function of cell chirality in morphogenesis.

Role of the Connexin C-terminus in skin pattern formation of Zebrafish

Background

Zebrafish display a striped skin pattern on their body; two types of connexins, namely, Connexin39.4 (Cx39.4) and Connexin41.8 (Cx41.8), are involved in stripe pattern formation. Herein, we investigated the role of the C-terminal (CT) domains of Cx39.4 and Cx41.8 in vivo and in vitro.

Methods

To investigate the role of CT domains in vivo, we established transgenic zebrafish lines expressing the CT-domain-modified connexin series in pigmented cells and observed skin patterns in fish. To investigate the role of the CT domains in vitro, we expressed the CT-domain modified connexin series in Neuro-2a (N2a) cells and calculated the plaque formation frequency.

Results

The overexpression of Cx39.4 lacking a CT domain produced skin patterns similar to that produced by full-length Cx39.4 in the cx39.4 ^-/- mutant and in cx39.4 and cx41.8 double-knockout mutant zebrafish. Fluorescence-protein-fused CT-domain-modified Cx39.4 formed gap junction plaques between N2a cells. The overexpression of CT-truncated Cx41.8 rescued the mutant phenotype in the cx41.8 ^-/- mutant but did not function in the double knockout zebrafish. Fluorescence-protein-fused CT-truncated Cx41.8 hardly formed plaques between N2a cells without Cx39.4 but formed gap junction plaques when co-expressed with Cx39.4.

Conclusions

The CT domain of Cx39.4 is not required for protein function, at least in the pigment cells of zebrafish. However, the need for the CT domain of Cx41.8 depends on Cx39.4 expression.

General significance

These results provide evidence for the interactions between Cx39.4 and Cx41.8 in pigment cells of zebrafish and suggest that at least one connexin must have a CT domain.

In situ differentiation of iridophore crystallotypes underlies zebrafish stripe patterning

Skin color patterns are ubiquitous in nature, impact social behavior, predator avoidance, and protection from ultraviolet irradiation. A leading model system for vertebrate skin patterning is the zebrafish; its alternating blue stripes and yellow interstripes depend on light-reflecting cells called iridophores. It was suggested that the zebrafish’s color pattern arises from a single type of iridophore migrating differentially to stripes and interstripes. However, here we find that iridophores do not migrate between stripes and interstripes but instead differentiate and proliferate in-place, based on their micro-environment. RNA-sequencing analysis further reveals that stripe and interstripe iridophores have different transcriptomic states, while cryogenic-scanning-electron-microscopy and micro-X-ray diffraction identify different crystal-arrays architectures, indicating that stripe and interstripe iridophores are different cell types. Based on these results, we present an alternative model of skin patterning in zebrafish in which distinct iridophore crystallotypes containing specialized, physiologically responsive, organelles arise in stripe and interstripe by in-situ differentiation.

The skin of zebrafish is patterned by alternating blue stripes and yellow interstripes which arises from guanine crystal-containing cells called iridophores that reflect light. Here the authors track iridophores and see that they do not migrate between stripes and interstripes, but instead differentiate and proliferate in place based on their micro-environment.

Evolution of the potassium channel gene Kcnj13 underlies colour pattern diversification in Danio fish

The genetic basis of morphological variation provides a major topic in evolutionary developmental biology. Fish of the genus Danio display colour patterns ranging from horizontal stripes, to vertical bars or spots. Stripe formation in zebrafish, Danio rerio, is a self-organizing process based on cell-contact mediated interactions between three types of chromatophores with a leading role of iridophores. Here we investigate genes known to regulate chromatophore interactions in zebrafish that might have evolved to produce a pattern of vertical bars in its sibling species, Danio aesculapii. Mutant D. aesculapii indicate a lower complexity in chromatophore interactions and a minor role of iridophores in patterning. Reciprocal hemizygosity tests identify the potassium channel gene obelix/Kcnj13 as evolved between the two species. Complementation tests suggest evolutionary change through divergence in Kcnj13 function in two additional Danio species. Thus, our results point towards repeated and independent evolution of this gene during colour pattern diversification.

Tracking collective cell motion by topological data analysis

By modifying and calibrating an active vertex model to experiments, we have simulated numerically a confluent cellular monolayer spreading on an empty space and the collision of two monolayers of different cells in an antagonistic migration assay. Cells are subject to inertial forces and to active forces that try to align their velocities with those of neighboring ones. In agreement with experiments in the literature, the spreading test exhibits formation of fingers in the moving interfaces, there appear swirls in the velocity field, and the polar order parameter and the correlation and swirl lengths increase with time. Numerical simulations show that cells inside the tissue have smaller area than those at the interface, which has been observed in recent experiments. In the antagonistic migration assay, a population of fluidlike Ras cells invades a population of wild type solidlike cells having shape parameters above and below the geometric critical value, respectively. Cell mixing or segregation depends on the junction tensions between different cells. We reproduce the experimentally observed antagonistic migration assays by assuming that a fraction of cells favor mixing, the others segregation, and that these cells are randomly distributed in space. To characterize and compare the structure of interfaces between cell types or of interfaces of spreading cellular monolayers in an automatic manner, we apply topological data analysis to experimental data and to results of our numerical simulations. We use time series of data generated by numerical simulations to automatically group, track and classify the advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These techniques of topological data analysis are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant for tissue engineering, biological effects of materials, tissue and organ regeneration.

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.

Linking genotype, cell behavior, and phenotype: multidisciplinary perspectives with a basis in zebrafish patterns

Zebrafish are characterized by dark and light stripes, but mutants display a rich variety of altered patterns. These patterns arise from the interactions of brightly colored pigment cells, making zebrafish a self-organization problem. The diversity of patterns present in zebrafish and other emerging fish models provides an excellent system for elucidating how genes, cell behavior, and visible animal characteristics are related. With the goal of highlighting how experimental and mathematical approaches can be used to link these scales, I overview current descriptions of zebrafish patterning, describe advances in the understanding of the mechanisms underlying cell communication, and discuss new work that moves beyond zebrafish to explore patterning in evolutionary relatives.

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.

Topological data analysis of zebrafish patterns

Self-organized pattern behavior is ubiquitous throughout nature, from fish schooling to collective cell dynamics during organism development. Qualitatively these patterns display impressive consistency, yet variability inevitably exists within pattern-forming systems on both microscopic and macroscopic scales. Quantifying variability and measuring pattern features can inform the underlying agent interactions and allow for predictive analyses. Nevertheless, current methods for analyzing patterns that arise from collective behavior capture only macroscopic features or rely on either manual inspection or smoothing algorithms that lose the underlying agent-based nature of the data. Here we introduce methods based on topological data analysis and interpretable machine learning for quantifying both agent-level features and global pattern attributes on a large scale. Because the zebrafish is a model organism for skin pattern formation, we focus specifically on analyzing its skin patterns as a means of illustrating our approach. Using a recent agent-based model, we simulate thousands of wild-type and mutant zebrafish patterns and apply our methodology to better understand pattern variability in zebrafish. Our methodology is able to quantify the differential impact of stochasticity in cell interactions on wild-type and mutant patterns, and we use our methods to predict stripe and spot statistics as a function of varying cellular communication. Our work provides an approach to automatically quantifying biological patterns and analyzing agent-based dynamics so that we can now answer critical questions in pattern formation at a much larger scale.

Two MYB Proteins in a Self-Organizing Activator-Inhibitor System Produce Spotted Pigmentation Patterns

Many organisms exhibit visually striking spotted or striped pigmentation patterns. Developmental models predict that such spatial patterns can form when a local autocatalytic feedback loop and a long-range inhibitory feedback loop interact. At its simplest, this self-organizing network only requires one self-activating activator that also activates a repressor, which inhibits the activator and diffuses to neighboring cells. However, the molecular activators and inhibitors fully fitting this versatile model remain elusive in pigmentation systems. Here, we characterize an R2R3-MYB activator and an R3-MYB repressor in monkeyflowers (Mimulus). Through experimental perturbation and mathematical modeling, we demonstrate that the properties of these two proteins correspond to an activator-inhibitor pair in a two-component, reaction-diffusion system, explaining the formation of dispersed anthocyanin spots in monkeyflower petals. Notably, disrupting this pattern impacts pollinator visitation. Thus, subtle changes in simple activator-inhibitor systems are likely essential contributors to the evolution of the remarkable diversity of pigmentation patterns in flowers.

Self-Organization in Pattern Formation

Self-organization is pervasive in development, from symmetry breaking in the early embryo to tissue patterning and morphogenesis. For a few model systems, the underlying molecular and cellular processes are now sufficiently characterized that mathematical models can be confronted with experiments, to explore the dynamics of pattern formation. Here, we review selected systems, ranging from cyanobacteria to mammals, where different forms of cell-cell communication, acting alone or together with positional cues, drive the patterning of cell fates, highlighting the insights that even very simple models can provide as well as the challenges on the path to a predictive understanding of development.

Neural crest development: insights from the zebrafish

Our understanding of the neural crest, a key vertebrate innovation, is built upon studies of multiple model organisms. Early research on neural crest cells (NCCs) was dominated by analyses of accessible amphibian and avian embryos, with mouse genetics providing complementary insights in more recent years. The zebrafish model is a relative newcomer to the field, yet it offers unparalleled advantages for the study of NCCs. Specifically, zebrafish provide powerful genetic and transgenic tools, coupled with rapidly developing transparent embryos that are ideal for high-resolution real-time imaging of the dynamic process of neural crest development. While the broad principles of neural crest development are largely conserved across vertebrate species, there are critical differences in anatomy, morphogenesis, and genetics that must be considered before information from one model is extrapolated to another. Here, our goal is to provide the reader with a helpful primer specific to neural crest development in the zebrafish model. We focus largely on the earliest events-specification, delamination, and migration-discussing what is known about zebrafish NCC development and how it differs from NCC development in non-teleost species, as well as highlighting current gaps in knowledge.

The identification of genes involved in the evolution of color patterns in fish

The genetic basis of morphological variation, both within and between species, provides a major topic in evolutionary biology. Teleost fish produce most elaborate color patterns, and among the more than 20 000 species a number have been chosen for more detailed analyses because they are suitable to study particular aspects of color pattern evolution. In several fish species, color variants and pattern variants have been collected, transcriptome analyses have been carried out, and the recent advent of gene editing tools, such as CRISPR/Cas9, has allowed the production of mutants. Covering mostly the literature from the last three years, we discuss the cellular basis of coloration and the identification of loci involved in color pattern differences between sister species in cichlids and Danio species, in which cis-regulatory changes seem to prevail.

The minimal gap-junction network among melanophores and xanthophores required for stripe-pattern formation in zebrafish

Connexin39.4 (Cx39.4) and Connexin41.8 (Cx41.8), two gap-junction proteins expressed in both melanophores and xanthophores, are critical for the intercellular communication among pigment cells that is necessary for generating the stripe pigment pattern of zebrafish. We previously characterized the gap-junction properties of Cx39.4 and Cx41.8, but how these proteins contribute to stripe formation remains unclear; this is because distinct types of connexins potentially form heteromeric gap junctions, which precludes accurate elucidation of individual connexin functions in vivo. Here, by arranging Cx39.4 and Cx41.8 expression in pigment cells, we identified the simplest gap-junction network required for stripe generation: Cx39.4 expression in melanophores is required but expression in xanthophores is not necessary for stripe patterning, whereas Cx41.8 expression in xanthophores is sufficient for the patterning, and Cx41.8 expression in melanophores might stabilize the stripes. Moreover, patch-clamp recordings revealed that Cx39.4 gap junctions exhibit spermidine-dependent rectification property. Our results suggest that Cx39.4 facilitates the critical cell-cell interactions between melanophores and xanthophores that mediate a unidirectional activation-signal transfer from xanthophores to melanophores, which is essential for melanophore survival.