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

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.

Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

Ultrastructural analysis of throat dermal tissue and chromatophore components in the threespine stickleback (Gasterosteus aculeatus)

The threespine stickleback (Gasterosteus aculeatus) is an important model for studying the evolution of nuptial coloration, but histological analyses of color are largely lacking. Previous analyses of one nuptial coloration trait, orange-red coloration along the body, have indicated carotenoids are the main pigment producing this color. In addition, recent gene expression studies found variation in the correlates of throat coloration between the sexes and between populations, raising the possibility of variation in the mechanisms underlying superficially similar coloration. We used transmission electron microscopy (TEM) to investigate the histological correlates of color in the throat dermal tissue of threespine stickleback from Western North America, within and between sexes, populations, and ecotypes. Ultrastructural analysis revealed carotenoid-containing erythrophores to be the main chromatophore component associated with orange-red coloration in both males and females across populations. In individuals where some darkening of the throat tissue was present, with no obvious orange-red coloration, erythrophores were not detected. Melanophore presence was more population-specific in expression, including being the only chromatophore component detected in a population of darker fish. We found no dermal chromatophore units within colorless throat tissue. This work confirms the importance of carotenoids and the erythrophore in producing orange-red coloration across sexes, as well as melanin within the melanophore in producing darkened coloration, but does not reveal broad histological differences among populations with similar coloration.

Myron Gordon Award Lecture 2023: Painting the neural crest: How studying pigment cells illuminates neural crest cell biology

It has been 30 (!!) years since I began working on zebrafish pigment cells, as a postdoc in the laboratory of Prof. Christiane Nüsslein-Volhard. There, I participated in the first large-scale mutagenesis screen in zebrafish, focusing on pigment cell mutant phenotypes. The isolation of colourless, shady, parade and choker mutants allowed us (as a postdoc in Prof. Judith Eisen's laboratory, and then in my own laboratory at the University of Bath since 1997) to pursue my ambition to address long-standing problems in the neural crest field. Thus, we have studied how neural crest cells choose individual fates, resulting in our recent proposal of a new, and potentially unifying, model which we call Cyclical Fate Restriction, as well as addressing how pigment cell patterns are generated. A key feature of our work in the last 10 years has been the use of mathematical modelling approaches to clarify our biological models and to refine our interpretations. None of this would have been possible without a hugely talented group of laboratory members and other collaborators from around the world-it has been, and I am sure will continue to be, a pleasure and privilege to work with you all!

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

kcnj13 regulates pigment cell shapes in zebrafish and has diverged by cis-regulatory evolution between Danio species

Teleost fish of the genus Danio are excellent models to study the genetic and cellular bases of pigment pattern variation in vertebrates. The two sister species Danio rerio and Danio aesculapii show divergent patterns of horizontal stripes and vertical bars that are partly caused by the divergence of the potassium channel gene kcnj13. Here, we show that kcnj13 is required only in melanophores for interactions with xanthophores and iridophores, which cause location-specific pigment cell shapes and thereby influence colour pattern and contrast in D. rerio. Cis-regulatory rather than protein coding changes underlie kcnj13 divergence between the two Danio species. Our results suggest that homotypic and heterotypic interactions between the pigment cells and their shapes diverged between species by quantitative changes in kcnj13 expression during pigment pattern diversification.

Craniofacial diversity across Danionins and the effects of TH status on craniofacial morphology of two Danio species

The model zebrafish ( Danio rerio ) belongs to the Danioninae subfamily with a range of informative phenotypes. However, the craniofacial diversity across the subfamily is not fully described. To better understand craniofacial phenotypes across Danioninae we used microCT and 3D geometric morphometrics to capture skull shapes from nine species. The Danio species examined showed largely similar skull shapes, although D. aesculapii , the sister species to D. rerio showed a unique morphology. Two non- Danio species examined, Chela dadiburjori and Devario aequipinnatus showed distinct skull morphologies unique from those of other species examined. Thyroid hormone regulates skeletal development and remodeling, and we asked if changes in developmental thyroid hormone metabolism could underlie some of the craniofacial diversity across Danioninae. We reared two Danio species under altered thyroid profiles, finding that hypothyroid individuals from both species showed corresponding morphological shifts in skull shape. Hypothyroid Danios showed skull morphologies closer to that of Chela and unlike any of the examined wild-type Danio species. We provide an examination of the evolved craniofacial diversity across Danioninae, and demonstrate that alterations to thyroid hormone have the capacity to create unique skull phenotypes.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Distinct genetic origins of eumelanin intensity and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. In order to understand genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

Effects of Natural and Synthetic Astaxanthin on Growth, Body Color, and Transcriptome and Metabolome Profiles in the Leopard Coral Grouper (Plectropomus leopardus)

Natural and synthetic astaxanthin can promote pigmentation in fish. In this study, the effects of dietary astaxanthin on growth and pigmentation were evaluated in leopard coral grouper (Plectropomus leopardus). Fish were assigned to three groups: 0% astaxanthin (C), 0.02% natural astaxanthin (HP), and 0.02% synthetic astaxanthin (AS). Brightness (L*) was not influenced by astaxanthin. However, redness (a*) and yellowness (b*) were significantly higher for fish fed astaxanthin-containing diets than fish fed control diets and were significantly higher in the HP group than in the AS group. In a transcriptome analysis, 466, 33, and 32 differentially expressed genes (DEGs) were identified between C and HP, C and AS, and AS and HP, including various pigmentation-related genes. DEGs were enriched for carotenoid deposition and other pathways related to skin color. A metabolome analysis revealed 377, 249, and 179 differential metabolites (DMs) between C and HP, C and AS, and AS and HP, respectively. In conclusion, natural astaxanthin has a better coloration effect on P. leopardus, which is more suitable as a red colorant in aquaculture. These results improve our understanding of the effects of natural and synthetic astaxanthin on red color formation in fish.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

A mechanism for red coloration in vertebrates

Red coloration is a salient feature of the natural world. Many vertebrates produce red color by converting dietary yellow carotenoids into red ketocarotenoids via an unknown mechanism. Here, we show that two enzymes, cytochrome P450 2J19 (CYP2J19) and 3-hydroxybutyrate dehydrogenase 1-like (BDH1L), are sufficient to catalyze this conversion. In birds, both enzymes are expressed at the sites of ketocarotenoid biosynthesis (feather follicles and red cone photoreceptors), and genetic evidence implicates these enzymes in yellow/red color variation in feathers. In fish, the homologs of CYP2J19 and BDH1L are required for ketocarotenoid production, and we show that these enzymes are sufficient to produce ketocarotenoids in cell culture and when ectopically expressed in fish skin. Finally, we demonstrate that the red-cone-enriched tetratricopeptide repeat protein 39B (TTC39B) enhances ketocarotenoid production when co-expressed with CYP2J19 and BDH1L. The discovery of this mechanism of ketocarotenoid biosynthesis has major implications for understanding the evolution of color diversity in vertebrates.

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.

Characterization of Collagen from Three Genetic Lines (Gray, Red and F1) of Oreochromis niloticus (Tilapia) Skin in Young and Old Adults

From tilapia (Oreochromis niloticus) farming, the by-products have been identified as a source of collagen that could be used for the development of dermocosmetics or pharmaceutical products. However, the characteristics of collagen related to a specific strain or culture must be well defined prior to its application. Collagen was extracted from the skin of three strains of tilapia: red YY males (YY: two Y-type sex chromosomes), XX gray females, and the F1: offspring of crossing red YY males with XX gray females; at different ages in the adult phase, using acetic acid and pepsin enzyme. The characteristics of acid-soluble collagen (ASC) and pepsin-soluble collagen (PSC) were shown by SDS-PAGE band profiles to be similar to bovine collagen type I (SIGMA), the PSC of gray tilapia being more fragile to temperature changes, consistent with the results of fractional viscosity. The characteristics of the F1 progeny were prioritized for being a commercially productive and sustainable source for the extraction of collagen, and the ASC form, being the one with the greatest stability and advantage over PSC, of importance to our investigations, leads to a controlled digestion as in the case of peptide induction, and also in the development of natural products in the pharmaceutical and/or dermocosmetic industry. Evaluations of the triple helix structure by FT-IR, X-ray diffraction and UV-visible spectroscopy give similar results between the strains: red, gray, and F1, and between ages in the adult form F1 (15, 24, and 36 months of age). Consequently, the skin of tilapia in adult form is recommended sustainably for up to 24 months of age where the collagen is obtained with the use of acetic acid without enzymatic treatment.

Analysis of pigment cell composition, pigment content, tyrosinase content and activity of three kinds of loaches Misgurnus anguillicaudatus from Poyang Lake

Pigment cell composition, pigment content, tyrosinase content and activity analysis were investigated on three kinds of loaches Misgurnus anguillicaudatus: big blackspot loaches (BBL), small blackspot loaches (SBL) and non-blackspot loaches (NBL), from Poyang Lake. Results showed that there were three types of skin pigment cells, namely melanophores, xanthophores and iridophores. Melanophores in dorsum were more than those in abdomen. Melanophore cytosomes in BBL were larger than those in SBL and NBL, and melanosomes were the largest in stage four. The melanophores in dorsal skin of SBL or NBL were small cell bodies, spindle-like and in chain distribution. There was an extremely significant difference in melanin content in BBL between the dorsum and abdomen (P < 0.01). There were no significant differences in melanin abdominal content, lutein and carotenoid contents among three kinds of loaches (P > 0.05). In dorsal skin, tyrosinase content was the highest in BBL, and it was significantly lower in NBL than in BBL and SBL (P < 0.01). This study reveals the differences in pigment and tyrosinase content in three kinds of loaches and provides a theoretical basis for further study of the mechanism of black spot formation.

Cell Fate Decisions in the Neural Crest, from Pigment Cell to Neural Development

The neural crest shows an astonishing multipotency, generating multiple neural derivatives, but also pigment cells, skeletogenic and other cell types. The question of how this process is controlled has been the subject of an ongoing debate for more than 35 years. Based upon new observations of zebrafish pigment cell development, we have recently proposed a novel, dynamic model that we believe goes some way to resolving the controversy. Here, we will firstly summarize the traditional models and the conflicts between them, before outlining our novel model. We will also examine our recent dynamic modelling studies, looking at how these reveal behaviors compatible with the biology proposed. We will then outline some of the implications of our model, looking at how it might modify our views of the processes of fate specification, differentiation, and commitment.

A highly contiguous nuclear genome assembly of the mandarinfish Synchiropus splendidus (Syngnathiformes: Callionymidae)

The fish order Syngnathiformes has been referred to as a collection of misfit fishes, comprising commercially important fish such as red mullets as well as the highly diverse seahorses, pipefishes, and seadragons-the well-known family Syngnathidae, with their unique adaptations including male pregnancy. Another ornate member of this order is the species mandarinfish. No less than two types of chromatophores have been discovered in the spectacularly colored mandarinfish: the cyanophore (producing blue color) and the dichromatic cyano-erythrophore (producing blue and red). The phylogenetic position of mandarinfish in Syngnathiformes, and their promise of additional genetic discoveries beyond the chromatophores, made mandarinfish an appealing target for whole-genome sequencing. We used linked sequences to create synthetic long reads, producing a highly contiguous genome assembly for the mandarinfish. The genome assembly comprises 483 Mbp (longest scaffold 29 Mbp), has an N50 of 12 Mbp, and an L50 of 14 scaffolds. The assembly completeness is also high, with 92.6% complete, 4.4% fragmented, and 2.9% missing out of 4584 BUSCO genes found in ray-finned fishes. Outside the family Syngnathidae, the mandarinfish represents one of the most contiguous syngnathiform genome assemblies to date. The mandarinfish genomic resource will likely serve as a high-quality outgroup to syngnathid fish, and furthermore for research on the genomic underpinnings of the evolution of novel pigmentation.

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.

Development and genetics of red coloration in the zebrafish relative Danio albolineatus

Animal pigment patterns play important roles in behavior and, in many species, red coloration serves as an honest signal of individual quality in mate choice. Among Danio fishes, some species develop erythrophores, pigment cells that contain red ketocarotenoids, whereas other species, like zebrafish (D. rerio) only have yellow xanthophores. Here, we use pearl danio (D. albolineatus) to assess the developmental origin of erythrophores and their mechanisms of differentiation. We show that erythrophores in the fin of D. albolineatus share a common progenitor with xanthophores and maintain plasticity in cell fate even after differentiation. We further identify the predominant ketocarotenoids that confer red coloration to erythrophores and use reverse genetics to pinpoint genes required for the differentiation and maintenance of these cells. Our analyses are a first step toward defining the mechanisms underlying the development of erythrophore-mediated red coloration in Danio and reveal striking parallels with the mechanism of red coloration in birds.

Immunoglobulin superfamily receptor Junctional adhesion molecule 3 (Jam3) requirement for melanophore survival and patterning during formation of zebrafish stripes

Adhesive interactions are essential for tissue patterning and morphogenesis yet difficult to study owing to functional redundancies across genes and gene families. A useful system in which to dissect roles for cell adhesion and adhesion-dependent signaling is the pattern formed by pigment cells in skin of adult zebrafish, in which stripes represent the arrangement of neural crest derived melanophores, cells homologous to melanocytes. In a forward genetic screen for adult pattern defects, we isolated the pissarro (psr) mutant, having a variegated phenotype of spots, as well as defects in adult fin and lens. We show that psr corresponds to junctional adhesion protein 3b (jam3b) encoding a zebrafish orthologue of the two immunoglobulin-like domain receptor JAM3 (JAM-C), known for roles in adhesion and signaling in other developing tissues, and for promoting metastatic behavior of human and murine melanoma cells. We found that zebrafish jam3b is expressed post-embryonically in a variety of cells including melanophores, and that jam3b mutants have defects in melanophore survival. Jam3b supported aggregation of cells in vitro and was required autonomously by melanophores for an adherent phenotype in vivo. Genetic analyses further indicated both overlapping and non-overlapping functions with the related receptor, Immunoglobulin superfamily 11 (Igsf11) and Kit receptor tyrosine kinase. These findings suggest a model for Jam3b function in zebrafish melanophores and hint at the complexity of adhesive interactions underlying pattern formation.

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.