A quantitative modelling approach to zebrafish pigment pattern formation

Phenylthiourea-mediated experimental depigmentation reduces seizurogenic response of pentylenetetrazol in zebrafish larva

Zebrafish larvae exposed to chemoconvulsants show behavioral seizures and electrographic abnormalities similar to the other mammalian models, making it a potential tool in epilepsy research. During the embryonic stage, zebrafish remains transparent which enables real-time developmental detection and in-situ gene/protein expression. However, pigmentation during the larval stage restricts transparency. Phenylthiourea (1-phenyl-2-thiourea; PTU) is a commonly used pigmentation blocker that maintains larval transparency. It is widely used along with chemoconvulsants to study in situ expressions in epileptic larvae, however, its effect on seizures largely remains unknown. Therefore, in the present study, the effect of PTU-mediated depigmentation was studied on pentylenetetrazol (PTZ)-induced seizures in zebrafish larvae. After spawning, the fish embryos were subjected to standard depigmentation protocol using 0.13 mM PTU. At 7-days post fertilization seizures were induced using 8 mM PTZ. PTU exposure significantly reduced PTZ-mediated hyperactive responses indicated by decreased distance travelled and swimming velocity of the larvae. Furthermore, PTU-exposed depigmented larvae also showed an increase in the latency to the onset of PTZ-mediated clonic-like seizures. The results concluded that PTU depigmentation protocol reduces the seizurogenic response of PTZ, hence its usage for imaging zebrafish larvae must be carefully monitored to avoid erroneous results.

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.

Identifying in vivo genetic dependencies of melanocyte and melanoma development

The advent of large-scale sequencing in both development and disease has identified large numbers of candidate genes that may be linked to important phenotypes. Validating the function of these candidates in vivo is challenging, due to low efficiency and low throughput of most model systems. We have developed a rapid, scalable system for assessing the role of candidate genes using zebrafish. We generated transgenic zebrafish in which Cas9 was knocked-in to the endogenous mitfa locus, a master transcription factor of the melanocyte lineage. We used this system to identify both cell-autonomous and non-cell autonomous regulators of normal melanocyte development. We then applied this to the melanoma setting to demonstrate that loss of genes required for melanocyte survival can paradoxically promote more aggressive phenotypes, highlighting that in vitro screens can mask in vivo phenotypes. Our high-efficiency genetic approach offers a versatile tool for exploring developmental processes and disease mechanisms that can readily be applied to other cell lineages.

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.

Multiple gene-drug prediction tool reveals Rosiglitazone based treatment pathway for non-segmental vitiligo

Vitiligo is a skin disease characterized by selective loss of melanocytes, which seriously affects the appearance and causes great psychological stress to patients. In this study, we performed a comprehensive analysis of two vitiligo microarray datasets from the GEO database using bioinformatics tools to identify 297 up-regulated mRNAs and 186 down-regulated mRNAs, revealing important roles for pathways related to melanin synthesis, tyrosine metabolism, and inflammatory factors, such as "PPAR signaling pathway", "tyrosine metabolism", "nonalcoholic fatty liver disease (NAFLD) pathway", "melanogenesis", and "IL-17 signaling pathway". Combining the Search Tool for Interacting Chemicals (STITCH) database 5.0 and the drug-gene interaction database 3.0 (DGIdb), we identified that the PPAR-γ agonist rosiglitazone may promote melanin synthesis via EDNRB. Next, we investigated the mechanism of rosiglitazone and PPAR-γ pathway in promoting melanin production. Consistent with the results of bioinformatics analysis, the expression levels of PPAR-γ, EDNRB, and TYR were significantly reduced in human non-segmental vitiligo skin along with the reduction of MITF, a key gene for epidermal melanogenesis. Meanwhile, rosiglitazone increased melanin synthesis capacity in melanocytes and zebrafish by activating PPAR-γ and upregulating TYR, TYRP-1, and TYRP-2. Conversely, treatment of melanocytes with the PPAR-γ antagonist GW resulted in inhibition of melanin synthesis and expression of melanin-related factors. At the same time, simultaneous treatment of rosiglitazone with GW reversed the inhibitory effect of GW on melanin synthesis. In this study, we identified that rosiglitazone, an important insulin sensitizer, promotes melanin synthesis in melanocytes by increasing PPAR-γ activity and upregulating the expression levels of EDNRB and TYR. These findings may provide new ideas for exploring the pathogenesis and potential therapeutic targets of non-segmental vitiligo.

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.

Assessing the exploratory profile of two zebrafish populations: influence of anxiety-like phenotypes and independent trials on homebase-related parameters and exploration

Anxiety is a protective behavior when animals face aversive conditions. The open field test (OFT) is used to assess the spatio-temporal dynamics of exploration, in which both homebase formation and recognition of environmental cues may reflect habituation to unfamiliar conditions. Because emotional- and affective-like states influence exploration patterns and mnemonic aspects, we aimed to verify whether the exploratory behaviors of two zebrafish populations showing distinct baselines of anxiety differ in two OFT sessions. Firstly, we assessed the baseline anxiety-like responses of short fin (SF) and leopard (LEO) populations using the novel tank test (NTT) and light-dark test (LDT) in 6-min trials. Fish were later tested in two consecutive days in the OFT, in which the spatial occupancy and exploratory profile were analyzed for 30min. In general, LEO showed pronounced diving behavior and scototaxis in the NTT and LDT, respectively, in which an "anxiety index" corroborated their exacerbated anxiety-like behavior. In the OFT, the SF population spent less time to establish the homebase in the 1^st trial, while only LEO showed a markedly reduction in the latency to homebase formation in the 2^nd trial. Both locomotion and homebase-related activities were decreased in the 2^nd trial, in which animals also revealed increased occupancy in the center area of the apparatus. Moreover, we verified a significant percentage of homebase conservation for both populations, while only SF showed reduced the number of trips and increased the average length of trips. Principal component analyses revealed that distinct factors accounted for total variances between trials for each population tested. While homebase exploration was reduced in the 2^nd trial for SF, an increased occupancy in the center area and hypolocomotion were the main factors that contribute to the effects observed in LEO during re-exposure to the OFT. In conclusion, our novel data support the homebase conservation in zebrafish subjected to independent OFT sessions, as well as corroborate a population-dependent effect on specific behavioral parameters related to exploration.

Fish-specific N-terminal domain sequence in Connexin 39.4 plays an important role in zebrafish stripe formation by regulating the opening and closing of gap junctions and hemichannels

BACKGROUND

Connexin 39.4 (Cx39.4) is involved in zebrafish (Danio rerio) skin patterning; when mutated, zebrafish exhibit a wavy stripe/labyrinth pattern instead of stripes. Cx39.4 is unique in that it has two additional serine/arginine (SR) residues, Ser2 and Arg3, at positions 2 and 3. Here, I investigated the role of these SR residues in Cx39.4 function.

METHODS

To examine the SR residues in Cx39.4, mutants of the SR residues were generated. Voltage-clamp recordings were performed using Xenopus oocytes to characterize the channel properties of the mutants. Transgenic zebrafish expressing each mutant were generated, and the effects of each mutation on fish skin patterning were evaluated.

RESULTS

The Cx39.4R3K mutant showed essentially the same properties as the wild-type (Cx39.4WT) in both electrophysiological analyses, leading to transgenic, complete phenotype rescue. Both the Cx39.4R3A mutant and deletion mutant of SR residues (Cx39.4delSR) showed a faster decay of gap junction activity and abnormal hemichannel activity, resulting in wide stripes and interstripes that indicate instability. Although the Cx39.4R3D mutant showed no channel activity in gap junctions or hemichannels, it caused unstable phenotypes in the transgene, namely a completely rescued phenotype in some individuals and loss of melanophores in others.

CONCLUSIONS

The SR residues in the NT domain of Cx39.4 are critical for the regulation of channel function, which appears to determine skin patterning.

GENERAL SIGNIFICANCE

These results elucidate the roles of the two SR residues unique to the NT domain of Cx39.4 in its channel function, which is important for zebrafish stripe pattern formation.

Aldh2 is a lineage-specific metabolic gatekeeper in melanocyte stem cells

Melanocyte stem cells (McSCs) in zebrafish serve as an on-demand source of melanocytes during growth and regeneration, but metabolic programs associated with their activation and regenerative processes are not well known. Here, using live imaging coupled with scRNA-sequencing, we discovered that, during regeneration, quiescent McSCs activate a dormant embryonic neural crest transcriptional program followed by an aldehyde dehydrogenase (Aldh) 2 metabolic switch to generate progeny. Unexpectedly, although ALDH2 is well known for its aldehyde-clearing mechanisms, we find that, in regenerating McSCs, Aldh2 activity is required to generate formate – the one-carbon (1C) building block for nucleotide biosynthesis – through formaldehyde metabolism. Consequently, we find that disrupting the 1C cycle with low doses of methotrexate causes melanocyte regeneration defects. In the absence of Aldh2, we find that purines are the metabolic end product sufficient for activated McSCs to generate progeny. Together, our work reveals McSCs undergo a two-step cell state transition during regeneration, and that the reaction products of Aldh2 enzymes have tissue-specific stem cell functions that meet metabolic demands in regeneration.

Summary: In zebrafish melanocyte regeneration, quiescent McSCs respond by re-expressing a neural crest identity, followed by an Aldh2-dependent metabolic switch to generate progeny.

Oca2 targeting using CRISPR/Cas9 in the Malawi cichlid Astatotilapia calliptera

Identifying genetic loci underlying trait variation provides insights into the mechanisms of diversification, but demonstrating causality and characterizing the role of genetic loci requires testing candidate gene function, often in non-model species. Here we establish CRISPR/Cas9 editing in Astatotilapia calliptera, a generalist cichlid of the remarkably diverse Lake Malawi radiation. By targeting the gene oca2 required for melanin synthesis in other vertebrate species, we show efficient editing and germline transmission. Gene edits include indels in the coding region, probably a result of non-homologous end joining, and a large deletion in the 3' untranslated region due to homology-directed repair. We find that oca2 knock-out A. calliptera lack melanin, which may be useful for developmental imaging in embryos and studying colour pattern formation in adults. As A. calliptera resembles the presumed generalist ancestor of the Lake Malawi cichlids radiation, establishing genome editing in this species will facilitate investigating speciation, adaptation and trait diversification in this textbook radiation.

Isolation of Chromatophores from Brown Trout (Salmo trutta) Skin

As new techniques for single-cell analysis are evolving, there is a growing need for streamlined and flexible protocols for single-cell suspension preparation and single-cell isolation. Based on previous research on the model species zebrafish (Danio rerio), a novel protocol for the isolation of chromatophores from brown trout (Salmo trutta) skin was developed and specific gene expression in chromatophore types evaluated with qPCR. This is the first report on a protocol for the isolation of chromatophores from the skin of salmonids. The method is simple and fast (taking about 90 min), requires no special equipment and is in line with the high ethical standards required in animal research, as no animals need to be sacrificed. In addition, skin biopsies can be taken from the same fish at multiple time points, allowing monitoring of chromatophore differentiation and their involvement in the formation of pigmentation patterns. The protocol should be effective for the vast majority of salmonid species due to similarities in skin morphology.

Tfap2b specifies an embryonic melanocyte stem cell that retains adult multifate potential

Melanocytes, the pigment-producing cells, are replenished from multiple stem cell niches in adult tissue. Although pigmentation traits are known risk factors for melanoma, we know little about melanocyte stem cell (McSC) populations other than hair follicle McSCs and lack key lineage markers with which to identify McSCs and study their function. Here we find that Tfap2b and a select set of target genes specify an McSC population at the dorsal root ganglia in zebrafish. Functionally, Tfap2b is required for only a few late-stage embryonic melanocytes, and is essential for McSC-dependent melanocyte regeneration. Fate mapping data reveal that tfap2b⁺ McSCs have multifate potential, and are the cells of origin for large patches of adult melanocytes, two other pigment cell types (iridophores and xanthophores), and nerve-associated cells. Hence, Tfap2b confers McSC identity in early development, distinguishing McSCs from other neural crest and pigment cell lineages, and retains multifate potential in the adult zebrafish.

Generation of a Transgenic Zebrafish Line for In Vivo Assessment of Hepatic Apoptosis

Hepatic apoptosis is involved in a variety of pathophysiologic conditions in the liver, including hepatitis, steatosis, and drug-induced liver injury. The development of easy-to-perform and reliable in vivo assays would thus greatly enhance the efforts to understand liver diseases and identify associated genes and potential drugs. In this study, we developed a transgenic zebrafish line that was suitable for the assessment of caspase 3 activity in the liver by using in vivo fluorescence imaging. The larvae of transgenic zebrafish dominantly expressed Casper3GR in the liver under control of the promoter of the phosphoenolpyruvate carboxykinase 1 gene. Casper3GR is composed of two fluorescent proteins, tagGFP and tagRFP, which are connected via a peptide linker that can be cleaved by activated caspase 3. Under tagGFP excitation conditions in zebrafish that were exposed to the well-characterized hepatotoxicant isoniazid, we detected increased and decreased fluorescence associated with tagGFP and tagRFP, respectively. This result suggests that isoniazid activates caspase 3 in the zebrafish liver, which digests the linker between tagGFP and tagRFP, resulting in a reduction in the Förster resonance energy transfer to tagRFP upon tagGFP excitation. We also detected isoniazid-induced inhibition of caspase 3 activity in zebrafish that were treated with the hepatoprotectants ursodeoxycholic acid and obeticholic acid. The transgenic zebrafish that were developed in this study could be a powerful tool for identifying both hepatotoxic and hepatoprotective drugs, as well as for analyzing the effects of the genes of interest to hepatic apoptosis.

Whole-organism 3D quantitative characterization of zebrafish melanin by silver deposition micro-CT

We previously described X-ray histotomography, a high-resolution, non-destructive form of X-ray microtomography (micro-CT) imaging customized for three-dimensional (3D), digital histology, allowing quantitative, volumetric tissue and organismal phenotyping (Ding et al., 2019). Here, we have combined micro-CT with a novel application of ionic silver staining to characterize melanin distribution in whole zebrafish larvae. The resulting images enabled whole-body, computational analyses of regional melanin content and morphology. Normalized micro-CT reconstructions of silver-stained fish consistently reproduced pigment patterns seen by light microscopy, and further allowed direct quantitative comparisons of melanin content across wild-type and mutant samples, including subtle phenotypes not previously noticed. Silver staining of melanin for micro-CT provides proof-of-principle for whole-body, 3D computational phenomic analysis of a specific cell type at cellular resolution, with potential applications in other model organisms and melanocytic neoplasms. Advances such as this in whole-organism, high-resolution phenotyping provide superior context for studying the phenotypic effects of genetic, disease, and environmental variables.

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.

Nile Tilapia: A Model for Studying Teleost Color Patterns

The diverse color patterns of cichlid fishes play an important role in mate choice and speciation. Here we develop the Nile tilapia (Oreochromis niloticus) as a model system for studying the developmental genetics of cichlid color patterns. We identified 4 types of pigment cells: melanophores, xanthophores, iridophores and erythrophores, and characterized their first appearance in wild-type fish. We mutated 25 genes involved in melanogenesis, pteridine metabolism, and the carotenoid absorption and cleavage pathways. Among the 25 mutated genes, 13 genes had a phenotype in both the F0 and F2 generations. None of F1 heterozygotes had phenotype. By comparing the color pattern of our mutants with that of red tilapia (Oreochromis spp), a natural mutant produced during hybridization of tilapia species, we found that the pigmentation of the body and eye is controlled by different genes. Previously studied genes like mitf, kita/kitlga, pmel, tyrb, hps4, gch2, csf1ra, pax7b, and bco2b were proved to be of great significance for color patterning in tilapia. Our results suggested that tilapia, a fish with 4 types of pigment cells and a vertically barred wild-type color pattern, together with various natural and artificially induced color gene mutants, can serve as an excellent model system for study color patterning in vertebrates.

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.

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.