- Open Access
Bar, stripe and spot development in sand-dwelling cichlids from Lake Malawi
EvoDevo volume 10, Article number: 18 (2019)
Melanic patterns such as horizontal stripes, vertical bars and spots are common among teleost fishes and often serve roles in camouflage or mimicry. Extensive research in the zebrafish model has shown that the development of horizontal stripes depends on complex cellular interactions between melanophores, xanthophores and iridophores. Little is known about the development of horizontal stripes in other teleosts, and even less is known about bar or spot development. Here, we compare chromatophore composition and development of stripes, bars and spots in two cichlid species of sand-dwellers from Lake Malawi—Copadichromis azureus and Dimidiochromis compressiceps.
(1) In D. compressiceps, stripes are made of dense melanophores underlaid by xanthophores and overlaid by iridophores. Melanophores and xanthophores are either loose or absent in interstripes, and iridophores are dense. In C. azureus, spots and bars are composed of a chromatophore arrangement similar to that of stripes but are separated by interbars where density of melanophores and xanthophores is only slightly lower than in stripes and iridophore density appears slightly greater. (2) Stripe, bar and spot chromatophores appear in the skin at metamorphosis. Stripe melanophores directly differentiate along horizontal myosepta into the adult pattern. In contrast, bar number and position are dynamic throughout development. As body length increases, new bars appear between old ones or by splitting of old ones through new melanophore appearance, not migration. Xanthophore and iridophore distributions follow melanophore patterns. (3) Metamorphic pigmentation arises in cichlids in a fashion similar to that described in zebrafish: melanophore progenitors derived from the medial route of neural crest migration migrate from the vicinity of the neural tube to the skin during metamorphosis.
The three pigment cell types forming stripes, bars and spots arise in the skin at metamorphosis. Stripes develop by differentiation of melanophores along horizontal myosepta, while bars do not develop along patent anatomical boundaries and increase in number in relation with body size. We propose that metamorphic melanophore differentiation and migratory arrest upon arrival to the skin lead to stripe formation, while bar formation must be supported by extensive migration of undifferentiated melanophores in the skin.
The evolutionary diversification of body coloration is a fascinating subject because body color is at the center of multiple selective pressures, such as light protection, female mate choice and predation. As such, coloration likely plays an important role in the creation and maintenance of biodiversity, yet the genetic and developmental changes required for the evolution of new patterns with adaptive value remain largely unknown.
Various mechanisms underlie the diversity of colors and color patterns observed among vertebrates. The melanocyte is the only pigment cell type, or chromatophore, found in birds and mammals . Melanocytes may produce different kinds of melanin pigment to produce colors ranging from yellowish (pheomelanin) to black (eumelanin). Other colors seen in birds and mammals result from food-derived pigments interacting with skin extracellular matrix [2, 3]. A greater diversity of chromatophore types underlies body coloration in reptiles, amphibians and fish, with seven types described in fish: xanthophores, erythrophores, iridophores, leucophores, melanophores, cyanophores and erythro-iridophores [4,5,6,7].
Vertebrate chromatophores are derived from an embryonic population of pluripotent and highly migratory cells called neural crest (NC) cells. These cells emerge from the dorsal neural tube toward the end of neurulation and migrate throughout the embryo to give rise to the skeletal and connective tissues of the face, neurons and glia of the peripheral nervous system and chromatophores, among many other derivatives . Trunk NC cells migrate through two distinct routes previously thought to be associated with distinct fates. Chromatophore progenitors migrate along a “lateral” route between surface ectoderm and somite, while neuronal and glial progenitors migrate along a “medial” route between somite and neural tube [9,10,11]. Recent evidence has shown that chromatophores also arise from NC cells that travel along the medial route in zebrafish and give rise to melanophores and iridophores during metamorphosis or regeneration [12, 13]. It is unknown if this mechanism is conserved in other teleosts or other vertebrate clades.
The highly diverse patterns of chordate body colors result from the particular distribution of pigments and/or chromatophores and in many instances of their particular combinations and overlapping patterns in the skin. The numerous studies focusing on zebrafish horizontal stripe formation have made the developmental and genetic basis of this complex pattern the best understood of all animals. Three chromatophore types, black melanophores, yellow xanthophores and silver/blue iridophores make up the zebrafish stripe pattern. Dark stripes are composed of all three chromatophore types with melanophores providing darkness, while iridophores provide iridescence and blue hue or green hue when combined with xanthophores. Light interstripes are made of xanthophores overlaid on iridophores [14, 15]. This adult pattern arises gradually during metamorphosis, spanning the period of 3 to 6/7 weeks post-fertilization [16, 17]. A key finding in zebrafish stripe development is that the distribution of each chromatophore depends on interactions with other chromatophore types. Melanophores prevent the spreading of dense iridophores into dark stripe territory, iridophores and xanthophores promote melanophore numbers, yet xanthophores also prevent spreading of iridophores into dark stripe regions, maintain the integrity of melanophore stripes and prevent spreading of melanophores in light stripe regions [18,19,20,21]. It is thus crucial to consider all chromatophore types when attempting to understand body color development in other organisms. While horizontal stripes are commonly found in other teleost species, the development of other color pattern elements, such as vertical bars and spots, remains poorly studied. It is unknown if the cell–cell interactions described in danios exist in other genera and if so, how they generate complex body coloration. In addition, only a handful of studies have identified the types of genetic and developmental changes driving the evolution of such complex patterns [22,23,24].
The great diversity of color patterns seen in the cichlid species flocks of East African rift lakes makes them great models to elucidate the developmental basis of previously unexplored patterns. The young ages of the species assemblies found in lakes Victoria (500,000 years) and Malawi (4 million years; ) make them amenable to genetic mapping by allowing the production of fertile hybrids in laboratory conditions. Moreover, the presence of highly similar yet independently evolved color patterns in species from these lakes offers great opportunities to elucidate the mechanisms underlying convergent evolution  as recently demonstrated by Kratochwill et al. .
Here, we describe and compare the development of the distinct color patterns of two cichlid species from Lake Malawi, Copadichromis azureus and Dimidiochromis compressiceps, in order to identify the developmental differences underlying their divergent color patterns. The color pattern of C. azureus includes three dark spots and a series of vertical dark bars. The D. compressiceps pattern includes three horizontal dark stripes. Similar to zebrafish, we find that three pigment cell types underlie the color patterns of both species, melanophores, xanthophores and iridophores, with different distributions in D. compressiceps and C. azureus adults. The time of first appearance of each chromatophore type is similar in D. compressiceps and C. azureus. Development of the adult pattern begins during larval development, essentially defining the onset of metamorphosis. Xanthophores appear first in both species, followed by melanophores, and finally iridophores. D. compressiceps stripes lay directly on top of horizontal myotome boundaries, while C. azureus vertical stripes do not clearly correlate with patent anatomical boundaries. We demonstrate by in situ hybridization that migration of NC cells takes place between 42 h post-fertilization (hpf) and 54hpf in both species along the lateral and medial routes and that, similar to zebrafish, pigment progenitors migrate from the vicinity of the neural tube to the skin during metamorphosis.
1. Adult pigmentation in C. azureus and D. compressiceps
The C. azureus body pigment pattern consists of the superimposition of two patterns commonly found in Lake Malawi, but also lakes Victoria and Tanganyika [26,27,28]: (1) an alternating series of dark bars and light interbars, and (2) three black spots (S1–3; Fig. 1a). Young adults have nine bars (b1–9) and eight interbars (i1–8; Fig. 1a), while more bars are observed on larger individuals. A series of five melanic patches variably aligned with vertical bars are present at the base of the dorsal fin. Individuals can modify the contrast between bars and interbars by changing bar darkness. A silver sheen covers the entire body with the exception of the spots, suggesting the widespread presence of iridophores. This pattern is seen in both males and females, although the silver sheen becomes metallic blue when males become dominant. In order to identify (1) what types of chromatophores are present in C. azureus and (2) how their particular arrangement creates the C. azureus pattern, we observed the skin of live individuals at various optical magnifications. Because pigment cells overlap one another in the skin, live specimens were photographed before and after l-adrenaline treatment to induce melanin granule aggregation in order to obtain a clear relationship between the abundance of particular chromatophores and the larger patterns observed on adult specimens: bars, interbars and spots. Three pigment cell types were identified: melanophores (dark brown/black), xanthophores (yellow-orange) and iridophores (silvery-white; Fig. 1b–g). All three chromatophore types are found in both the dark bars and the light interbars (Fig. 1d, e), as well as the spots (Fig. 1f, g). Melanophores were 23% more abundant in bars (99.4 cell/mm2) than interbars (81 cell/mm2; Fig. 1h), and cell size was larger in bars (Fig. 1d). Xanthophore abundance was 11% greater in bars (61.3 cell/mm2) than interbars (55.3 cell/mm2; Fig. 1d–e, h). Iridophores were generally overlaid on melanophores in spots, bars and interbars (Fig. 1d, f), where they seemed to provide greater cover although they could not be counted. Melanophore and xanthophore abundance was similar in spots as in bars (89.7 cell/mm2 and 68.3 cell/mm2; Fig. 1g, h), although melanophores were darker in spots (Fig. 1i).
The D. compressiceps body pigment pattern is also found on numerous cichlids from lakes Malawi, Tanganyika and Victoria and includes three horizontal dark stripes separated by light silvery interstripes (Fig. 1j, k; [27, 28]). An additional dark stripe runs along the dorsal midline (Fig. 1j, k). Dark stripes tend to disappear in dominant males as a metallic blue sheen develops over the entire body. As in C. azureus, melanophores, xanthophores and iridophores were identified by microscopy on live specimens before and after l-adrenaline treatment (Fig. 1k–o). All dark stripes contain the three chromatophores, while dorsal interstripes contain iridophores and melanophores, and ventral interstripes contain iridophores alone (Fig. 1k–o). Differences in melanophore density differentiate the stripes: melanophore density is the highest in the midline stripe (SML), intermediate in the dorsolateral (SDL) and dorsal (SD) stripes and lowest in the ventrolateral stripe (SVL), resulting in a discontinuous appearance (Fig. 1j–p).
2. Pre-metamorphic pigmentation development in D. compressiceps and C. azureus
In order to understand how each species-specific pigment pattern develops, live specimens from both species were imaged daily throughout embryonic and larval development, which ends after roughly 2 weeks of development when the yolk becomes fully depleted. Rates of development are extremely similar in D. compressiceps and C. azureus, and very similar to that described for the Nile tilapia . Melanophores are the first pigment cells to differentiate in both species starting at 42hpf on the yolk (data not shown), followed by a 4-day period with minimal change to this pigmentation pattern (Fig. 2a–i, m), except in the eye where epithelial melanization and iridophores were first observed at 4dpf (Fig. 2a, e) and 6dpf (Fig. 2c, g), respectively.
A new wave of skin pigmentation started at 6dpf (Fig. 2c, g), a day before overt morphological changes indicating the start of metamorphosis in both species such as separation of dorsal, caudal and anal fins and appearance of caudal fin ray elements (CFRE) (Fig. 3). Melanophores and xanthophores become visible on the dorsal head at 6dpf (Fig. 2c, g), followed by a progressive expansion of xanthophores alone along the dorsal trunk and tail at 7dpf (Fig. 4a–h). Xanthophore expansion proceeded ventrally over the next 3 days, reaching the lateral midline by 8dpf/2CFRE (Fig. 4e–h) and the ventral midline by 10dpf/5CFRE in both species (Fig. 4i–l). Following this full body coverage, xanthophore distribution became concentrated under bar melanophores in C. azureus (Fig. 4m, n) or stripe melanophores in D. compressiceps (Fig. 4o, p) by the end of the larval stage (9CFRE).
Unlike xanthophores, metamorphic melanophores in the trunk (anterior to anus) and/or tail (posterior to anus) first appear deep in larvae of both species at 8dpf/2CFRE—on the dorsal side of the neural tube (Fig. 2i, m). These deep melanophores appear along the entire length of the neural tube in C. azureus, while they are restricted to the tail region in D. compressiceps at this stage. Skin melanophores are first observed the next day (9dpf/3CFRE) in both species—at the base of the dorsal fin, and at different regions along the lateral midline (Figs. 2j, n and 5a, b).
3. Metamorphic pigment pattern development in C. azureus
Three distinct patches (P) that prefigure the adult spots, but also the midparts of presumptive bars b4, b7 and b9 are observed in the skin at 9dpf/3CFRE along the lateral midline of C. azureus: P1 on the trunk, P2 on the tail, and P3 on the caudal fin peduncle (Figs. 2j and 5a). New melanophore appearance within and around these three patches leads to increased melanophore density as well as progressive expansion of each patch until 12dpf/6CFRE (Figs. 2j–l, q and 5a, c). By 13dpf/7CFRE new melanophore appearance leads to the formation of 2 new patches: P4 between P1 and P2, and P5 between P2 and P3 (Figs. 2r and 5e). Interestingly, appearance of these new melanophore patches takes place between differentiating iridophore streaks (Figs. 2r, s and 5e). At 15dpf/8CFRE, a new melanophore patch, P6, appears anterior to P1 (Fig. 2s). The progressive definition of C. azureus bars and interbars occurs concomitantly with the appearance of a series of dark and light patches along the dorsal midline. Similar to lateral midline patches, the first pigment cells appearing at the dorsal midline are melanophores concentrated in three melanophore patches (D1–3) at the base of the dorsal fin at 11dpf/5CFRE (Figs. 2l and 5c). New melanophores increase the cell density around these three centers, and iridophores start appearing in the two regions separating the patches by 6CFRE/12dpf (Fig. 2q). By 13dpf/7CFRE, D1 becomes resegmented into D1a and D1b, and D2 into D2a and D2b by new iridophore patches (Figs. 2r and 5e). One more melanophore dorsal patch (D4) flanked by iridophores also appears on the caudal fin peduncle (Figs. 2r and 5e). Some, but not all of these dorsal patches are dorso-ventrally aligned with the presumptive bars developing during the same period over C. azureus flanks.
In order to determine whether lateral patch melanophores migrate ventrally from dorsal patches as previously proposed in Mesoamerican cichlids  and more generally what cell behaviors underlie stripe morphogenesis, 4 C. azureus larvae were imaged daily from 10dpf/4CFRE to 18dpf/9CFRE at high magnification to follow individual melanophore behavior during bar formation (Fig. 6). Most melanophores did not move after appearing in the skin (Fig. 6a–i), and few mitoses were observed. Bars are thus shaped by the addition of new melanophores in C. azureus.
4. Metamorphic pigment pattern development in D. compressiceps
In contrast to C. azureus, the posterior half of the midlateral stripe in D. compressiceps first appears at 9dpf/3CFRE as a continuous stripe of melanophores along the lateral midline in the tail region (posterior to anus; Figs. 2n, 5b). Melanophores are also observed at the base of the dorsal fin at this stage, prefiguring the dorsal stripe (Fig. 5b). Melanophore density increases each day by the addition of new melanophores to posterior SML and SD (Figs. 2o, p, u–x and 5b, d, f, h). Melanophores forming the anterior half of SML appear dorsal to the yolk at 12dpf/6CFRE (Fig. 2u), along with those forming the dorsal (SD) and ventrolateral stripes (SVL). New melanophore additions through later stages increase cell density in each stripe without changing the pattern prefiguring the adult stripes (Fig. 2u–x).
Interstripe iridophores appear dorsal and ventral to the midlateral stripe at 13dpf/7CFRE (Figs. 2v and 5f), and new iridophore additions through later stages give rise to interstripes X1DL and X1VL (Figs. 2x and 5h). The adult D. compressiceps pigment pattern is essentially complete by the end of the larval period (17dpf/9CFRE) (Figs. 2x and 5h).
5. Timing of chromatophore progenitor development
Based on the observed onset of adult pigmentation at metamorphosis in D. compressiceps and C. azureus, we next asked what are the embryonic origins of metamorphic chromatophores in these cichlids. In zebrafish, metamorphic melanophore and iridophore progenitors arise from stem cell pools associated with dorsal root ganglia, which are derived from ventral/medial route NC cells and replace embryonic/larval pigment cells [12, 13]. These progenitors migrate to the skin along sensory axons traveling between myotomes . In contrast, larval xanthophores are derived from dorsal/lateral route NC cells and remain in the skin at metamorphosis and into adulthood .
In order to determine when chromatophore progenitors arise and migrate in D. compressiceps and C. azureus, we conducted in situ hybridization to detect the expression of NC- and chromatophore-specific markers at various developmental stages in D. compressiceps and C. azureus embryos. FoxD3 and sox10 label premigratory and migratory NC cells, respectively [30,31,32,33]. Kit-l, mitfa, tyr and tyr-p label melanophore progenitors [34,35,36]; csf1ra, gch and xdh label xanthophore progenitors [18, 37, 38]; pnp4 and bmp4 label iridophore progenitors . In addition, ednrb1a labels all three chromatophore types in zebrafish , as well as migratory NC cells in mammals .
Gene expression patterns were first examined at 36, 42, 48, 54 and 60hpf based on previously described timing of cranial NC migration in another cichlid, the Nile Tilapia . Sox10 and foxD3 expression was detected at 36 and 42hpf in C. azureus (Fig. 7a–d) and D. compressiceps (data not shown). At 36hpf, foxD3 is expressed in the hindbrain and in premigratory cranial NC cells, as well as two streams of NC cells just beginning to migrate bilaterally, which we interpret as the mandibular NC stream (directly posterior to the eye) and the post-otic NC stream (posterior to otic vesicle) (Fig. 7a). Sox10 expression was also detected in the hindbrain and cranial NC streams, albeit at lower intensity (Fig. 7b).
At 42hpf, foxd3 expression was still detected in all three cranial NC streams, as well as in two distinct cell populations in the anterior trunk. One population was found ventrally between the neural tube and the first two somites, which we interpret as NC cells migrating along the medial/ventral route. A second population consisted of individual cells on top of the neural tube, which we interpret as NC cells migrating along the lateral route (Fig. 7c). At this and later stages, laterally migrating NC cells were present further posteriorly than medial/ventrally migrating NC cells, consistent with an earlier onset of lateral migration during the anterior–posterior sequence of trunk NC migration, as described in other vertebrates, including zebrafish . foxD3 was also detected in the medial parts of somites 3–7, as previously described in zebrafish . Sox10 expression at these stages was restricted to the pre- and post-otic streams of cranial NC but was also detected in medially migrating NC at the levels of somites 1–2, and in laterally migrating NC cells further posteriorly, similar to foxd3 (Fig. 7d).
Trunk expression of foxD3 and sox10 at subsequent stages (48, 54hpf) was detected in successively more posterior streams of NC cells. A notable difference between the medial and laterally migrating populations was the duration of foxd3 and sox10 expression: expression of both was maintained in NC streams migrating medially but downregulated in streams migrating laterally (Fig. 7e, f, h–k). Ednrb1a expression was also detected at 48 and 54hpf in both species. In the trunk, ednrb1a expression labeled medial NC streams similar to foxd3 and sox10 (Fig. 7g, l, m), and also lateral NC streams at 54hpf (Fig. 7l). Expression of foxd3, sox10 or ednrb1a decreased by 60hpf and was no longer detected at 66 or 72hpf, and none of the other markers listed above produced detectable signals from 36 to 72hpf.
The undifferentiated progenitors of zebrafish metamorphic melanophores and iridophores reside at dorsal root ganglia until metamorphosis, when they migrate along nerve axons and between myotomes to the skin . In order to determine whether similar pigment progenitors exist in cichlids, we conducted in situ hybridization on cryosections at premetamorphosis (96hpf) and metamorphosis stages (11dpf/5CFRE) using ednrb1a as a pan-chromatophore marker and mitfa as a melanophore marker. Expression of both genes was first detected at 11dpf/5CFRE (Fig. 7n, o) along the neural tube, between myotomes, and in the skin, consistent with a migration of chromatophore progenitors from the vicinity of the neural tube to the skin during metamorphosis. No difference in mitfa or ednrb1a expression was detected between D. compressiceps and C. azureus at these stages.
Chromatophore distribution relative to anatomical landmarks
How do spatial cues instruct pigment cells to differentiate in particular regions of the embryo or larva? Previous findings in zebrafish have shown that metamorphic melanophore progenitors (melanoblasts) travel from dorsal root ganglia to the skin along axons between myotomes . Here, we provide evidence for a similar mechanism in cichlids. This suggests a simple mechanism for horizontal stripe formation along horizontal myosepta: melanoblasts migrating to the skin along horizontal myosepta could give rise to horizontal stripes by differentiating in the skin upon arrival. In contrast, additional melanoblast migration would be required in the skin for vertical bar formation.
Dimidiochromis compressiceps dark stripes are clearly overlaid on horizontal myosepta: SML melanophores differentiate along the lateral midline myoseptum, stripe SD along the dorsal midline myoseptum, and stripes SDL and SVL along the dorsal and ventral lateral myosepta, respectively (Fig. 8b, d, f, h). In contrast, C. azureus vertical bars or spots arise from a progressively greater number of lateral melanophore patches that progressively elongate dorso-ventrally from new melanophores (Fig. 8a, c, e, g), without clear relationships with underlying myosepta. A major developmental difference resulting in the different color patterns of C. azureus and D. compressiceps may thus be in the migratory behavior of melanoblasts in the skin (Fig. 8i). Accordingly, the D. compressiceps pattern may result from the differentiation of melanoblasts into melanophores with limited migratory potential once they reach the skin, while the C. azureus pattern would require the migration of melanoblasts in the skin before differentiation into melanophores (Fig. 8i), since differentiated melanophores are distributed throughout the skin regardless of myoseptum position.
Developmental regulation of bar number in old-world and new-world cichlids
Vertical bars similar to those found in C. azureus are found in many cichlid species from both Africa and the Americas. Bar ontogeny was described in 40 species of neotropical cichlids of the subfamily Cichlasomatinae with various numbers of bars at adults stages . Interestingly, the authors found that 9 vertical bars (including 2 in the head) develop first in all species examined, followed by species-specific bar divisions or fusions. Similar to C. azureus, bars form progressively starting with an initial pattern of lateral and dorsal midline melanophore clusters in Cichlasomatinae cichlid larvae. Dorsal and/or ventral migration of melanophores away from the clusters was reported to lead to the formation of bars. Our observation of individual melanophores in the skin of individual C. azureus larvae through several days of development did not show such migration of differentiated melanophores. Instead, supported exit points of metamorphic melanoblasts at horizontal myosepta suggest dorsal and ventral migration of melanophores through the skin in the undifferentiated state to result in the vertical bar pattern of C. azureus.
The timing of metamorphosis was also recently experimentally shown to modulate the number of vertical bars in a Mesoamerican cichlid, the convict cichlid Amatitlania nigrofasciata . Consistent with a perturbation of metamorphosis , it was found that the timing and level of thyroid signaling activity affect (1) the timing of larval to juvenile coloration, and  the number of vertical bars. Individuals under elevated thyroid hormone levels from embryo to adult stages start pigmentation metamorphosis earlier than untreated individuals and develop 7 bars instead of 9, while individuals with lower thyroid hormone levels start pigmentation metamorphosis later and end up with 12 bars.
How does the timing of metamorphosis affect the number of bars? A common developmental pattern reported in barred Mesoamerican cichlids and here in C. azureus is the resolution of a few broad melanophore patches into numerous thin bars from larval to adult stages. In C. azureus new bars appear in growing larvae by new melanophore appearance in the skin—not rearrangement of already present melanophores. New bars form between old ones that have become spaced apart through body growth. Body growth also correlates with bar splitting—it creates new space between bar melanophores, which is filled by new melanophores within daughter bars but remains vacant between daughter bars. These observations suggest that increase in physical space along the anterior–posterior axis is a positive regulator of bar number. The increase in bar number associated with delayed metamorphosis in the Mesoamerican cichlid A. nigrofasciata  also suggests that physical space at the time of first bar appearance is also a critical regulator of bar number. A. nigrofasciata individuals with delayed metamorphosis developed roughly the same number of bars as C. azureus at juvenile stages; however, these bars were thinner and thus more widely spaced, perhaps leaving more room for later-appearing bars.
Genetic determinants of cichlid color patterns
The ontogeny of zebrafish body coloration is by far the best understood of all vertebrates because of the numerous mutant screens that have led to the identification of genes required for melanophore, xanthophore and iridophore development, but also for the interactions among and between these cell types and with their environment [47, 48]. However, this understanding is largely restricted to the development of horizontal stripes. The diversity of cichlid color patterns in closely related species and their amenability to genetic mapping makes them a great model system to explore the developmental and genetic underpinnings of pigment pattern evolution in vertebrates. Research efforts have so far been far greater to identify the genetic factors regulating cichlid pigmentation than their developmental basis. Although most naturally occurring trait variants are under the control of multiple genes, as illustrated in Albertson et al. , two instances of monogenic traits were discovered in East African cichlids. Cis-regulatory variants affecting pax7a expression were found to be associated with the development of the orange-blotch camouflage pattern in Mbuna cichlids of Lake Malawi [50, 51]. More recently, a cis-regulatory variant affecting agrp2 expression was shown to determine the presence/absence of the horizontal dark stripe found in lake Victoria cichlids Haplochromis sauvagei and Pundamilia nyererei . These authors also showed that difference in agrp2 expression reliably predicted the presence/absence of dark stripe in cichlids from other East African Lakes, thus demonstrating a clear case of convergence under the control of a single gene .
Here, we identified potential developmental differences underlying pigment pattern divergence in D. compressiceps and C. azureus. Identification of the loci regulating these developmental differences will constitute the next step toward a broader understanding of the developmental evolution of body coloration in vertebrates.
The three same chromatophore types make up D. compressiceps and C. azureus body coloration: melanophores, iridophores and xanthophores. In D. compressiceps, stripes are made of dense melanophores underlaid by xanthophores and overlaid by iridophores. Melanophores and xanthophores are either loose or absent in interstripes, and iridophores are dense. In C. azureus, spots and bars are composed of a chromatophore arrangement similar to that of stripes but are separated by interbars where density of melanophores and xanthophores is only slightly lower than in stripes and iridophore density appears slightly greater.
Differences in pigmentation between C. azureus and D. compressiceps arise at metamorphosis when melanophores appear in the skin as arrangements prefiguring stripes, bars or spots. In D. compressiceps, stripe melanophores differentiate along horizontal myosepta and the adult pattern is essentially in place by the end of the larval period. In contrast, the number of bars developing in C. azureus larvae is smaller than in adults. Bars appear as a series of lateral melanophore patches that progressively elongate dorso-ventrally by new melanophore additions. New bars appear either between old ones, or through splitting of old ones. The place of new melanophore appearance, not migration, determines new bar formation. Increase in xanthophore density in developing bars and stripes follows that of melanophores. Iridophores appear in interstripes after stripe formation, while their appearance in interbars is concomitant with bar formation.
Analysis of neural crest migration by ISH indicates that melanophore progenitors migrate from the vicinity of the neural tube to the skin along myosepta at metamorphosis, as previously described in zebrafish . Time-lapse analysis of bar development during metamorphosis shows that melanophores do not migrate after differentiation. We thus propose that metamorphic melanophore differentiation and migratory arrest upon arrival to the skin lead to stripe formation, while bar formation must be supported by extensive dorsal and/or ventral migration of undifferentiated melanophores in the skin.
All animals were reared and euthanized following protocols approved by the IACUC at the University of North Carolina, Wilmington. D. compressiceps and C. azureus brood stocks were purchased from a pet store and maintained separately at 28 °C in 55- or 75-gallon tanks. Adults were fed Northfin Food Cichlid Formula 3 mm sinking pellets. Fertilized eggs were collected from mouth-brooding females a few hours after fertilization. Embryos were then staged under a microscope, reared in cichlid egg tumblers and imaged daily until the end of the larval period (yolk depletion). Embryos and larvae were staged following Fujimura and Okada, 2007 .
For imaging of live specimens, individuals were anesthetized in tank water using Tricaine (MP Biomedicals) following protocols approved by the IACUC at the University of North Carolina, Wilmington and immobilized on a drop of 3% methyl cellulose (Alfa Aesar). l-Adrenaline treatment (Alfa Aesar) was carried out under the same conditions at a concentration of 4.5 mg/mL. Chromatophore density averages reported in Fig. 1h, p were obtained from regions representative of each pattern in C. azureus (n = 3) or D. compressiceps (n = 3). To measure and compare melanophore darkness between bars and spots (Fig. 1i), images were  converted to greyscale and  inverted in Photoshop. Integrated density was then measured in Fiji, with darker regions displaying higher integrated density values than lighter regions. For time-lapse imaging of C. azureus larvae, 4 specimens were imaged daily from 10dpf/4CFRE to 18dpf/9CFRE at high magnification to follow individual melanophore behavior during bar formation. Each specimen was anesthetized daily as described above, imaged and returned to its tank. Daily image capture was determined to be adequate to follow new melanophore appearance after multiple-daily image captures showed barely noticeable changes to the melanophore pattern. Melanophore migratory activity was assessed by recording  the position of melanophores relative to stable landmarks (myosepta) and  the relative positions of closest neighbors over the entire duration of the time-lapse imaging. The replacement of a single melanophore by 2 melanophores in successive time points was interpreted as a mitosis. All observations and images were made using a dissecting microscope Leica M165FC equipped with Planapo 1.6X M-series objective and Prior L200 fluorescent light source. Image capture was made with a Leica DFC7000T camera controlled by the LAS X software. Chromatophore counts were performed in FIJI. Image processing was performed in Adobe Photoshop.
In situ hybridization
Embryos and larvae were fixed and processed for whole mount in situ hybridization as described in Le Pabic et al. : embryos were fixed in 4% PFA for 3 days at 4 °C, rinsed in PBT and transferred to 100% methanol at − 20 °C for storage. For cryosections, embryos and larvae were processed as described in Albertson et al. . In situ hybridization probes were generated from RT-PCR amplified regions of C. azureus foxd3, sox10, kit-l, mitfa, tyr, tyr-p, csf1ra, gch, xdh, pnp4, bmp4, ednrb1a cDNAs using the primers listed in Table 1. mRNA was extracted from deyolked 54hpf C. azureus embryos. In situ hybridization was conducted using antisense probes prepared from pGem-t-easy vectors (Promega) containing cDNA fragments corresponding to C. azureus foxd3, sox10, kit-l, mitfa, tyr, tyr-p, csf1ra, gch, xdh, pnp4, bmp4 and ednrb1a.
Availability of data and materials
The datasets used and/or analyzed during this study are available from the corresponding author on reasonable request.
Caro T. The colours of extant mammals. Semin Cell Dev Biol. 2013;24:542–52.
Roulin A, Ducrest AL. Genetics of colouration in birds. Semin Cell Dev Biol. 2013;24:594–608.
Kelsh RN, Harris ML, Colanesi S, Erickson CA. Stripes and belly-spots—a review of pigment cell morphogenesis in vertebrates. Semin Cell Dev Biol. 2009;20:90–104.
Olsson M, Stuart-Fox D, Ballen C. Genetics and evolution of colour patterns in reptiles. Semin Cell Dev Biol. 2013;24:529–41.
Grether GF, Kolluru GR, Nersissian K. Individual colour patches as multicomponent signals. Biol Rev. 2004;79:583–610.
Fujii R. The regulation of motile activity in fish chromatophores. Pigment Cell Res. 2000;13:300–19.
Goda M, et al. Integumental reddish-violet coloration owing to novel dichromatic chromatophores in the teleost fish, Pseudochromis diadema. Pigment Cell Melanoma Res. 2011;24:614–7.
Le Douarin NM, Kalcheim C. The neural crest. 2nd ed. Cambridge: Cambridge University Press; 1999.
Erickson CA, Weston JA. An SEM analysis of neural crest migration in the mouse. J Embryol Exp Morphol. 1983;74:97–118.
Raible DW, et al. Segregation and early dispersal of neural crest cells in the embryonic zebrafish. Dev Dyn. 1992;195:29–42.
Weston JA, Butler SL. Temporal factors affecting localization of neural crest cells in chicken embryo. Dev Biol. 1966;14:246.
Budi EH, Patterson LB, Parichy DM. Post-embryonic nerve-associated precursors to adult pigment cells: genetic requirements and dynamics of morphogenesis and differentiation. Plos Genet. 2011;7:e1002044.
Dooley CM, Mongera A, Walderich B, Nusslein-Volhard C. On the embryonic origin of adult melanophores: the role of ErbB and Kit signalling in establishing melanophore stem cells in zebrafish. Development. 2013;140:1003–13.
Hirata M, Nakamura K, Kanemaru T, Shibata Y, Kondo S. Pigment cell organization in the hypodermis of zebrafish. Dev Dyn. 2003;227:497–503.
Hirata M, Nakamura K, Kondo S. Pigment cell distributions in different tissues of the zebrafish, with special reference to the striped pigment pattern. Dev Dyn. 2005;234:293–300.
Singh AP, Schach U, Nusslein-Volhard C. Proliferation, dispersal and patterned aggregation of iridophores in the skin prefigure striped colouration of zebrafish. Nat Cell Biol. 2014;16:604.
Mahalwar P, Walderich B, Singh AP, Nusslein-Volhard C. Local reorganization of xanthophores fine-tunes and colors the striped pattern of zebrafish. Science. 2014;345:1362–4.
Maderspacher F, Nusslein-Volhard C. Formation of the adult pigment pattern in zebrafish requires leopard and obelix dependent cell interactions. Development. 2003;130:3447–57.
Frohnhofer HG, Krauss J, Maischein HM, Nusslein-Volhard C. Iridophores and their interactions with other chromatophores are required for stripe formation in zebrafish. Development. 2013;140:2997–3007.
Patterson LB, Parichy DM. Interactions with iridophores and the tissue environment required for patterning melanophores and xanthophores during zebrafish adult pigment stripe formation. Plos Genet. 2013;9:e1003561.
Nakamasu A, Takahashi G, Kanbe A, Kondo S. Interactions between zebrafish pigment cells responsible for the generation of Turing patterns. Proc Natl Acad Sci USA. 2009;106:8429–34.
Manceau M, Domingues VS, Mallarino R, Hoekstra HE. The developmental role of agouti in color pattern evolution. Science. 2011;331:1062–5.
Kratochwil CF, et al. Agouti-related peptide 2 facilitates convergent evolution of stripe patterns across cichlid fish radiations. Science. 2018;362:457.
Spiewak JE, et al. Evolution of endothelin signaling and diversification of adult pigment pattern in Danio fishes. Plos Genet. 2018;14:e1007538.
Delvaux D. Musée royal de l’Afrique Centrale (Tervuren), Département de Géologie et Minéralogie, rapport annuel, 1995, vol. 1993; 1994. p. 99–108.
Kocher TD, Conroy JA, McKaye KR, Stauffer JR. Similar morphologies of cichlid fish in lakes Tanganyika and Malawi are due to convergence. Mol Phylogenet Evol. 1993;2:158–65.
Barel CDN, van Oijen MJP, Witte F, Witte-Maas ELM. An introduction to the taxonomy and morphology of the haplochromine cichlids from Lake Victoria. Neth J Zool. 1977;27:381–9.
Seehausen O, Mayhew PJ, Van Alphen JJM. Evolution of colour patterns in East African cichlid fish. J Evol Biol. 1999;12:514–34.
Rican O, Musilova Z, Muska M, Novak J. Development of coloration patterns in neotropical cichlids (Teleostei: Cichlidae: Cichlasomatinae). Folia Zool. 2005;54:1–46.
Kelsh RN, Dutton K, Medlin J, Eisen JS. Expression of zebrafish fkd6 in neural crest-derived glia. Mech Dev. 2000;93:161–4.
Lister JA, et al. Zebrafish Foxd3 is required for development of a subset of neural crest derivatives. Dev Biol. 2006;290:92–104.
Stewart RA, et al. Zebrafish foxd3 is selectively required for neural crest specification, migration and survival. Dev Biol. 2006;292:174–88.
Dutton KA, et al. Zebrafish colourless encodes sox10 and specifies non-ectomesenchymal neural crest fates. Development. 2001;128:4113–25.
Parichy DM, Rawls JF, Pratt SJ, Whitfield TT, Johnson SL. Zebrafish sparse corresponds to an orthologue of c-kit and is required for the morphogenesis of a subpopulation of melanocytes, but is not essential for hematopoiesis or primordial germ cell development. Development. 1999;126:3425–36.
Lister JA, Robertson CP, Lepage T, Johnson SL, Raible DW. Nacre encodes a zebrafish microphthalmia-related protein that regulates neural-crest-derived pigment cell fate. Development. 1999;126:3757–67.
Johnson SL, Nguyen AN, Lister JA. Mitfa is required at multiple stages of melanocyte differentiation but not to establish the melanocyte stem cell. Dev Biol. 2011;350:405–13.
Parichy DM, Ransom DG, Paw B, Zon LI, Johnson SL. An orthologue of the kit-related gene fms is required for development of neural crest-derived xanthophores and a subpopulation of adult melanocytes in the zebrafish, Danio rerio. Development. 2000;127:3031–44.
Parichy DM, Turner JM. Temporal and cellular requirements for Fms signaling during zebrafish adult pigment pattern development. Development. 2003;130:817–33.
Lang MR, Patterson LB, Gordon TN, Johnson SL, Parichy DM. Basonuclin-2 requirements for zebrafish adult pigment pattern development and female fertility. Plos Genet. 2009;5:e1000744.
Parichy DM, et al. Mutational analysis of endothelin receptor b1 (rose) during neural crest and pigment pattern development in the zebrafish Danio rerio. Dev Biol. 2000;227:294–306.
Lee HO, Levorse JM, Shin MK. The endothelin receptor-B is required for the migration of neural crest-derived melanocyte and enteric neuron precursors. Dev Biol. 2003;259:162–75.
Le Pabic P, Stellwag EJ, Scemama J-L. Embryonic development and skeletogenesis of the pharyngeal jaw apparatus in the cichlid Nile tilapia (Oreochromis niloticus). Anat Rec. 2009;292:1780–800.
Raible DW, Eisen JS. Specification and commitment in zebrafish trunk neural crest development. Dev Biol. 1994;163:536.
Odenthal J, Nusslein-Volhard C. Fork head domain genes in zebrafish. Dev Genes Evol. 1998;208:245–58.
Prazdnikov DV, Shkil FN. Experimental evidence of the role of heterochrony in evolution of the Mesoamerican cichlids pigment patterns. Evol Dev. 2019;21:3–15.
McMenamin SK, Parichy DM. Metamorphosis in teleosts. Anim Metamorph. 2013;103:127–65.
Singh AP, Nusslein-Volhard C. Zebrafish stripes as a model for vertebrate colour pattern formation. Curr Biol. 2015;25:R81–92.
Parichy DM, Spiewak JE. Origins of adult pigmentation: diversity in pigment stem cell lineages and implications for pattern evolution. Pigment Cell Melanoma Res. 2015;28:31–50.
Albertson RC, et al. Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes. Mol Ecol. 2014;23:5135–50.
Roberts RB, Ser JR, Kocher TD. Sexual conflict resolved by invasion of a novel sex determiner in Lake Malawi cichlid fishes. Science. 2009;326:998–1001.
Roberts RB, Moore EC, Kocher TD. An allelic series at pax7a is associated with colour polymorphism diversity in Lake Malawi cichlid fish. Mol Ecol. 2017;26:2625–39.
Fujimura K, Okada N. Development of the embryo, larva and early juvenile of Nile tilapia Oreochromis niloticus (Pisces: Cichlidae). Developmental staging system. Dev Growth Differ. 2007;49:301–24.
Le Pabic P, Stellwag EJ, Brothers SN, Scemama J-L. Comparative analysis of Hox paralog group 2 gene expression during Nile tilapia (Oreochromis niloticus) embryonic development. Dev Genes Evol. 2007;217:749–58.
Albertson RC, et al. Molecular pedomorphism underlies craniofacial skeletal evolution in Antarctic notothenioid fishes. BMC Evol Biol. 2010;10:4.
Research reported in this study was supported by NIH/NIDCR R01DE13828 to PL and TFS.
Ethics approval and consent to participate
All animals were reared and euthanized following protocols approved by the IACUC at the University of North Carolina, Wilmington.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.