Pax7 is required for establishment of the xanthophore lineage in zebrafish embryos

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

Molecular characterization of Pax7 and its role in melanin synthesis in Crassostrea gigas

The paired-box 7 (Pax7) is a transcription factor crucial for skin color polymorphism. However, the mechanism underlying the pigmentation associated with Pax7 in mollusks have yet to be elucidated. In this study, the cDNA sequence of Pax7 in the Pacific oyster Crassostrea gigas (CgPax7) was characterized. Phylogenetically, the identity of deduced amino acid sequence was similar to that of other mollusks and contained 463 amino acids, with conserved features of paired domain (PRD), homeobox domain (HD) and octapeptide. Gene expression analysis revealed that CgPax7 was markedly increased at D-shaped larvae stage and ubiquitously expressed in six examined tissues in adult oyster. The result of whole-mount in situ hybridization (WMISH) showed a restricted pattern of CgPax7 expression on margins of shell valves at D-shaped and umbo larvae stages. Additionally, although CgPax7 silencing had no significant effect on CgMitf expression, it significantly inhibited the expressions of CgPax7, CgTyr, CgTyrp1, CgTyrp2 and CgCdk2, genes involved in Tyr-mediated melanin synthesis. Furthermore, CgPax7 knockdown obviously decreased the tyrosinase activity. Less brown-granules at mantle edge was detected by micrographic examination and melanosomes defect was observed by transmission electron microscopy. It was demonstrated that CgPax7 play a key role in melanin synthesis by regulating Tyr-pathway in C. gigas. These findings indicated the potential framework by which mollusks pigmentation.

Single-cell resolution of MET- and EMT-like programs in osteoblasts during zebrafish fin regeneration

Zebrafish regenerate fin rays following amputation through epimorphic regeneration, a process that has been proposed to involve the epithelial-to-mesenchymal transition (EMT). We performed single-cell RNA sequencing (scRNA-seq) to elucidate osteoblastic transcriptional programs during zebrafish caudal fin regeneration. We show that osteoprogenitors are enriched with components associated with EMT and its reverse, mesenchymal-to-epithelial transition (MET), and provide evidence that the EMT markers cdh11 and twist2 are co-expressed in dedifferentiating cells at the amputation stump at 1 dpa, and in differentiating osteoblastic cells in the regenerate, the latter of which are enriched in EMT signatures. We also show that esrp1, a regulator of alternative splicing in epithelial cells that is associated with MET, is expressed in a subset of osteoprogenitors during outgrowth. This study provides a single cell resource for the study of osteoblastic cells during zebrafish fin regeneration, and supports the contribution of MET- and EMT-associated components to this process.

•

Osteoblasts express EMT/MET signatures during zebrafish fin regeneration

•

De/re-differentiating osteoblasts express cdh11, an EMT marker

•

A subset of osteoprogenitors express the MET marker esrp1

•

Our scRNA-seq data can be explored online

Genetic and correlative light and electron microscopy evidence for the unique differentiation pathway of erythrophores in brown trout skin

Based on their cell ultrastructure, two types of erythrophores in the spotted skin regions of brown trout (Salmo trutta) were previously described. To test the hypothesis regarding the origin of a new cell type following genome duplication, we analysed the gene and paralogue gene expression patterns of erythrophores in brown trout skin. In addition, the ultrastructure of both erythrophore types was precisely examined using transmission electron microscopy (TEM) and correlative light microscopy and electron microscopy (CLEM). Ultrastructural differences between the sizes of erythrophore inclusions were confirmed; however, the overlapping inclusion sizes blur the distinction between erythrophore types, which we have instead defined as cell subtypes. Nevertheless, the red spots of brown trout skin with subtype 2 erythrophores, exhibited unique gene expression patterns. Many of the upregulated genes are involved in melanogenesis or xanthophore differentiation. In addition, sox10, related to progenitor cells, was also upregulated in the red spots. The expressions of paralogues derived from two genome duplication events were also analysed. Multiple paralogues were overexpressed in the red spots compared with other skin regions, suggesting that the duplicated gene copies adopted new functions and contributed to the origin of a new cell subtype that is characteristic for red spot. Possible mechanisms regarding erythrophore origin are proposed and discussed. To the best of our knowledge, this is the first study to evaluate pigment cell types in the black and red spots of brown trout skin using the advanced CLEM approach together with gene expression profiling.

Analysis of yellow mutant rainbow trout transcriptomes at different developmental stages reveals dynamic regulation of skin pigmentation genes

Yellow mutant rainbow trout (YR), an economically important aquaculture species, is popular among consumers due to its excellent meat quality and attractive appearance. Skin color is a key economic trait for YR, but little is known about the molecular mechanism of skin color development. In this study, YR skin transcriptomes were analyzed to explore temporal expression patterns of pigmentation-related genes in three different stages of skin color development. In total, 16,590, 16,682, and 5619 genes were differentially expressed between fish at 1 day post-hatching (YR1d) and YR45d, YR1d and YR90d, and YR45d and YR90d. Numerous differentially expressed genes (DEGs) associated with pigmentation were identified, and almost all of them involved in pteridine and carotenoid synthesis were significantly upregulated in YR45d and YR90d compared to YR1d, including GCH1, PTS, QDPR, CSFIR1, SLC2A11, SCARB1, DGAT2, PNPLA2, APOD, and BCO2. Interestingly, many DEGs enriched in melanin synthesis pathways were also significantly upregulated, including melanogenesis (MITF, MC1R, SLC45A2, OCA2, and GPR143), tyrosine metabolism (TYR, TYRP1, and DCT), and MAPK signaling (KITA) pathways. Using short time-series expression miner, we identified eight differential gene expression pattern profiles, and DEGs in profile 7 were associated with skin pigmentation. Protein–protein interaction network analysis showed that two modules were related to xanthophores and melanophores. In addition, 1,812,329 simple sequence repeats and 2,011,334 single-nucleotide polymorphisms were discovered. The results enhance our understanding of the molecular mechanism underlying skin pigmentation in YR, and could accelerate the molecular breeding of fish species with valuable skin color traits and will likely be highly informative for developing new therapeutic approaches to treat pigmentation disorders and melanoma.

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.

Genetic basis of orange spot formation in the guppy (Poecilia reticulata)

Background

To understand the evolutionary significance of female mate choice for colorful male ornamentation, the underlying regulatory mechanisms of such ornamentation must be understood for examining how the ornaments are associated with “male qualities” that increase the fitness or sexual attractiveness of offspring. In the guppy (Poecilia reticulata), an established model system for research on sexual selection, females prefer males possessing larger and more highly saturated orange spots as potential mates. Although previous studies have identified some chromosome regions and genes associated with orange spot formation, the regulation and involvement of these genetic elements in orange spot formation have not been elucidated. In this study, the expression patterns of genes specific to orange spots and certain color developmental stages were investigated using RNA-seq to reveal the genetic basis of orange spot formation.

Results

Comparing the gene expression levels of male guppy skin with orange spots (orange skin) with those without any color spots (dull skin) from the same individuals identified 1102 differentially expressed genes (DEGs), including 630 upregulated genes and 472 downregulated genes in the orange skin. Additionally, the gene expression levels of the whole trunk skin were compared among the three developmental stages and 2247 genes were identified as DEGs according to color development. These analyses indicated that secondary differentiation of xanthophores may affect orange spot formation.

Conclusions

The results suggested that orange spots might be formed by secondary differentiation, rather than de novo generation, of xanthophores, which is induced by Csf1 and thyroid hormone signaling pathways. Furthermore, we suggested candidate genes associated with the areas and saturation levels of orange spots, which are both believed to be important for female mate choice and independently regulated. This study provides insights into the genetic and cellular regulatory mechanisms underlying orange spot formation, which would help to elucidate how these processes are evolutionarily maintained as ornamental traits relevant to sexual selection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01942-2.

Cyclical fate restriction: a new view of neural crest cell fate specification

Neural crest cells are crucial in development, not least because of their remarkable multipotency. Early findings stimulated two hypotheses for how fate specification and commitment from fully multipotent neural crest cells might occur, progressive fate restriction (PFR) and direct fate restriction, differing in whether partially restricted intermediates were involved. Initially hotly debated, they remain unreconciled, although PFR has become favoured. However, testing of a PFR hypothesis of zebrafish pigment cell development refutes this view. We propose a novel 'cyclical fate restriction' hypothesis, based upon a more dynamic view of transcriptional states, reconciling the experimental evidence underpinning the traditional hypotheses.

Genetic compensation between Pax3 and Pax7 in zebrafish appendicular muscle formation

BACKGROUND

Migrating muscle progenitors delaminate from the somite and subsequently form muscle tissue in distant anatomical regions such as the paired appendages, or limbs. In amniotes, this process requires a signaling cascade including the transcription factor paired box 3 (Pax3).

RESULTS

In this study, we found that, unlike in mammals, pax3a/3b double mutant zebrafish develop near to normal appendicular muscle. By analyzing numerous mutant combinations of pax3a, pax3b and pax7a, and pax7b, we determined that there is a feedback system and a compensatory mechanism between Pax3 and Pax7 in this developmental process, even though Pax7 alone is not required for appendicular myogenesis. pax3a/3b/7a/7b quadruple mutant developed muscle-less pectoral fins.

CONCLUSIONS

We found that Pax3 and Pax7 are redundantly required during appendicular myogenesis in zebrafish, where Pax7 is able to activate the same developmental programs as Pax3 in the premigratory progenitor cells.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which multiple species share the same conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which vivid colour and pattern are produced in a Müllerian mimicry complex of poison frogs. We did this by first assembling a high-quality de novo genome assembly for the mimic poison frog Ranitomeya imitator. This assembled genome is 6.8 Gbp in size, with a contig N50 of 300 Kbp R. imitator and two colour morphs from both Ranitomeya fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes that are differentially expressed throughout development, many of them related to melanocyte development, melanin synthesis, iridophore development and guanine synthesis. Polytypic differences within species may be the result of differences in expression and/or timing of expression, whereas convergence for colour pattern between species does not appear to be due to the same changes in gene expression. In addition, we identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species. Finally, we hypothesize that genes in the keratin family are important for producing different structural colours within these frogs.

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.

An atlas of neural crest lineages along the posterior developing zebrafish at single-cell resolution

Neural crest cells (NCCs) are vertebrate stem cells that give rise to various cell types throughout the developing body in early life. Here, we utilized single-cell transcriptomic analyses to delineate NCC-derivatives along the posterior developing vertebrate, zebrafish, during the late embryonic to early larval stage, a period when NCCs are actively differentiating into distinct cellular lineages. We identified several major NCC/NCC-derived cell-types including mesenchyme, neural crest, neural, neuronal, glial, and pigment, from which we resolved over three dozen cellular subtypes. We dissected gene expression signatures of pigment progenitors delineating into chromatophore lineages, mesenchyme cells, and enteric NCCs transforming into enteric neurons. Global analysis of NCC derivatives revealed they were demarcated by combinatorial hox gene codes, with distinct profiles within neuronal cells. From these analyses, we present a comprehensive cell-type atlas that can be utilized as a valuable resource for further mechanistic and evolutionary investigations of NCC differentiation.

The MITF paralog tfec is required in neural crest development for fate specification of the iridophore lineage from a multipotent pigment cell progenitor

Understanding how fate specification of distinct cell-types from multipotent progenitors occurs is a fundamental question in embryology. Neural crest stem cells (NCSCs) generate extraordinarily diverse derivatives, including multiple neural, skeletogenic and pigment cell fates. Key transcription factors and extracellular signals specifying NCSC lineages remain to be identified, and we have only a little idea of how and when they function together to control fate. Zebrafish have three neural crest-derived pigment cell types, black melanocytes, light-reflecting iridophores and yellow xanthophores, which offer a powerful model for studying the molecular and cellular mechanisms of fate segregation. Mitfa has been identified as the master regulator of melanocyte fate. Here, we show that an Mitf-related transcription factor, Tfec, functions as master regulator of the iridophore fate. Surprisingly, our phenotypic analysis of tfec mutants demonstrates that Tfec also functions in the initial specification of all three pigment cell-types, although the melanocyte and xanthophore lineages recover later. We show that Mitfa represses tfec expression, revealing a likely mechanism contributing to the decision between melanocyte and iridophore fate. Our data are consistent with the long-standing proposal of a tripotent progenitor restricted to pigment cell fates. Moreover, we investigate activation, maintenance and function of tfec in multipotent NCSCs, demonstrating for the first time its role in the gene regulatory network forming and maintaining early neural crest cells. In summary, we build on our previous work to characterise the gene regulatory network governing iridophore development, establishing Tfec as the master regulator driving iridophore specification from multipotent progenitors, while shedding light on possible cellular mechanisms of progressive fate restriction.

An improved de novo assembling and polishing of Solea senegalensis transcriptome shed light on retinoic acid signalling in larvae

Senegalese sole is an economically important flatfish species in aquaculture and an attractive model to decipher the molecular mechanisms governing the severe transformations occurring during metamorphosis, where retinoic acid seems to play a key role in tissue remodeling. In this study, a robust sole transcriptome was envisaged by reducing the number of assembled libraries (27 out of 111 available), fine-tuning a new automated and reproducible set of workflows for de novo assembling based on several assemblers, and removing low confidence transcripts after mapping onto a sole female genome draft. From a total of 96 resulting assemblies, two "raw" transcriptomes, one containing only Illumina reads and another with Illumina and GS-FLX reads, were selected to provide SOLSEv5.0, the most informative transcriptome with low redundancy and devoid of most single-exon transcripts. It included both Illumina and GS-FLX reads and consisted of 51,348 transcripts of which 22,684 code for 17,429 different proteins described in databases, where 9527 were predicted as complete proteins. SOLSEv5.0 was used as reference for the study of retinoic acid (RA) signalling in sole larvae using drug treatments (DEAB, a RA synthesis blocker, and TTNPB, a RA-receptor agonist) for 24 and 48 h. Differential expression and functional interpretation were facilitated by an updated version of DEGenes Hunter. Acute exposure of both drugs triggered an intense, specific and transient response at 24 h but with hardly observable differences after 48 h at least in the DEAB treatments. Activation of RA signalling by TTNPB specifically increased the expression of genes in pathways related to RA degradation, retinol storage, carotenoid metabolism, homeostatic response and visual cycle, and also modified the expression of transcripts related to morphogenesis and collagen fibril organisation. In contrast, DEAB mainly decreased genes related to retinal production, impairing phototransduction signalling in the retina. A total of 755 transcripts mainly related to lipid metabolism, lipid transport and lipid homeostasis were altered in response to both treatments, indicating non-specific drug responses associated with intestinal absorption. These results indicate that a new assembling and transcript sieving were both necessary to provide a reliable transcriptome to identify the many aspects of RA action during sole development that are of relevance for sole aquaculture.

Absence of Desmin in Myofibers of the Zebrafish Extraocular Muscles

Purpose

To study the medial rectus (MR) muscle of zebrafish (Daniorerio) with respect to the pattern of distribution of desmin and its correlation to distinct types of myofibers and motor endplates.

Methods

The MRs of zebrafish were examined using confocal microscopy in whole-mount longitudinal specimens and in cross sections processed for immunohistochemistry with antibodies against desmin, myosin heavy chain isoforms, and innervation markers. Desmin patterns were correlated to major myofiber type and type of innervation. A total of 1382 myofibers in nine MR muscles were analyzed.

Results

Four distinct desmin immunolabeling patterns were found in the zebrafish MRs. Approximately a third of all slow myofibers lacked desmin, representing 8.5% of the total myofiber population. The adult zebrafish MR muscle displayed en grappe, en plaque, and multiterminal en plaque neuromuscular junctions (NMJs) with intricate patterns of desmin immunolabeling.

Conclusions

The MRs of zebrafish showed important similarities with the human extraocular muscles with regard to the pattern of desmin distribution and presence of the major types of NMJs and can be regarded as an adequate model to further study the role of desmin and the implications of heterogeneity in cytoskeletal protein composition.

Translational Relevance

The establishment of a zebrafish model to study the cytoskeleton in muscles that are particularly resistant to muscle disease opens new avenues to understand human myopathies and muscle dystrophies and may provide clues to new therapies.

Evaluation of hypopigmentation in embryonic zebrafish induced by emerging disinfection byproduct, 3, 5-di-I-tyrosylalanine

Halogenated dipeptides, 3, 5-di-I-tyrosylalanine (DIYA), have been identified as novel disinfection byproducts (DBPs), following chloramination of authentic water. However, little is known about their toxicity. Zebrafish embryos were used to assess the toxicity of novel iodinated DBPs (I-DBPs). Although DIYA did not exhibit high acute toxicity to embryonic zebrafish (LC₅₀ > 2 mM), it significantly inhibited pigmentation of melanophores and xanthophores on head, trunk and tail at 500 μM as determined by photographic analysis. Whereas N-phenylthiourea (PTU) as a pigment inhibitor did not inhibit development of yellow pigments. Colorimetric detection of melanin further confirmed these results. Quantitative real time polymerase chain reaction (qRT-PCR) measurements indicated that genes (dct, slc24a5, tyr, tyrp1a, tyrp1b, silva) associated with the melanogenesis pathway were dramatically down-regulated following exposure to 500 μM DIYA. In addition, enzymatic activity of tyrosinase (TYR) decreased, also demonstrating that the underlying mechanism of hypopigmentation was attributed to the disruption of melanogenesis pathway. Transcription levels of xanthophore genes (gch2, bnc2, csf1a, csf1b, pax7a and pax7b) were also monitored by qRT-PCR assay. DIYA exposure up-regulated expression of gch2 and bnc2, but not csf1 and pax7. Tested DIYA analogues, brominated tyrosine was unlikely to inhibit pigmentation, indicating that the iodine substitution and dipeptides structure are of important structural feature for the inhibition of pigmentation. In this study, we observed that DIYA inhibited melanogenesis related genes, which might contribute to pigmentation defects. Moreover, as an emerging I-DBPs, the developmental toxicity of aromatic dipeptides should be further studied.

Being red, blue and green: the genetic basis of coloration differences in the strawberry poison frog (Oophaga pumilio)

Background

Animal coloration is usually an adaptive attribute, under strong local selection pressures and often diversified among species or populations. The strawberry poison frog (Oophaga pumilio) shows an impressive array of color morphs across its distribution in Central America. Here we quantify gene expression and genetic variation to identify candidate genes involved in generating divergence in coloration between populations of red, green and blue O. pumilio from the Bocas del Toro archipelago in Panama.

Results

We generated a high quality non-redundant reference transcriptome by mapping the products of genome-guided and de novo transcriptome assemblies onto a re-scaffolded draft genome of O. pumilio. We then measured gene expression in individuals of the three color phenotypes and identified color-associated candidate genes by comparing differential expression results against a list of a priori gene sets for five different functional categories of coloration – pteridine synthesis, carotenoid synthesis, melanin synthesis, iridophore pathways (structural coloration), and chromatophore development. We found 68 candidate coloration loci with significant expression differences among the color phenotypes. Notable upregulated examples include pteridine synthesis genes spr, xdh and pts (in red and green frogs); carotenoid metabolism genes bco2 (in blue frogs), scarb1 (in red frogs), and guanine metabolism gene psat1 (in blue frogs). We detected significantly higher expression of the pteridine synthesis gene set in red and green frogs versus blue frogs. In addition to gene expression differences, we identified 370 outlier SNPs on 162 annotated genes showing signatures of diversifying selection, including eight pigmentation-associated genes.

Conclusions

Gene expression in the skin of the three populations of frogs with differing coloration is highly divergent. The strong signal of differential expression in pteridine genes is consistent with a major role of these genes in generating the coloration differences among the three morphs. However, the finding of differentially expressed genes across pathways and functional categories suggests that multiple mechanisms are responsible for the coloration differences, likely involving both pigmentary and structural coloration. In addition to regulatory differences, we found potential evidence of differential selection acting at the protein sequence level in several color-associated loci, which could contribute to the color polymorphism.

Comparative transcriptomics reveals candidate carotenoid color genes in an East African cichlid fish

BACKGROUND

Carotenoids contribute significantly to animal body coloration, including the spectacular color pattern diversity among fishes. Fish, as other animals, derive carotenoids from their diet. Following uptake, transport and metabolic conversion, carotenoids allocated to body coloration are deposited in the chromatophore cells of the integument. The genes involved in these processes are largely unknown. Using RNA-Sequencing, we tested for differential gene expression between carotenoid-colored and white skin regions of a cichlid fish, Tropheus duboisi "Maswa", to identify genes associated with carotenoid-based integumentary coloration. To control for positional gene expression differences that were independent of the presence/absence of carotenoid coloration, we conducted the same analyses in a closely related population, in which both body regions are white.

RESULTS

A larger number of genes (n = 50) showed higher expression in the yellow compared to the white skin tissue than vice versa (n = 9). Of particular interest was the elevated expression level of bco2a in the white skin samples, as the enzyme encoded by this gene catalyzes the cleavage of carotenoids into colorless derivatives. The set of genes with higher expression levels in the yellow region included genes involved in xanthophore formation (e.g., pax7 and sox10), intracellular pigment mobilization (e.g., tubb, vim, kif5b), as well as uptake (e.g., scarb1) and storage (e.g., plin6) of carotenoids, and metabolic conversion of lipids and retinoids (e.g., dgat2, pnpla2, akr1b1, dhrs). Triglyceride concentrations were similar in the yellow and white skin regions. Extracts of integumentary carotenoids contained zeaxanthin, lutein and beta-cryptoxanthin as well as unidentified carotenoid structures.

CONCLUSION

Our results suggest a role of carotenoid cleavage by Bco2 in fish integumentary coloration, analogous to previous findings in birds. The elevated expression of genes in carotenoid-rich skin regions with functions in retinol and lipid metabolism supports hypotheses concerning analogies and shared mechanisms between these metabolic pathways. Overlaps in the sets of differentially expressed genes (including dgat2, bscl2, faxdc2 and retsatl) between the present study and previous, comparable studies in other fish species provide useful hints to potential carotenoid color candidate genes.

Regulation of zebrafish melanocyte development by ligand-dependent BMP signaling

Preventing terminal differentiation is important in the development and progression of many cancers including melanoma. Recent identification of the BMP ligand GDF6 as a novel melanoma oncogene showed GDF6-activated BMP signaling suppresses differentiation of melanoma cells. Previous studies have identified roles for GDF6 orthologs during early embryonic and neural crest development, but have not identified direct regulation of melanocyte development by GDF6. Here, we investigate the BMP ligand gdf6a, a zebrafish ortholog of human GDF6, during the development of melanocytes from the neural crest. We establish that the loss of gdf6a or inhibition of BMP signaling during neural crest development disrupts normal pigment cell development, leading to an increase in the number of melanocytes and a corresponding decrease in iridophores, another neural crest-derived pigment cell type in zebrafish. This shift occurs as pigment cells arise from the neural crest and depends on mitfa, an ortholog of MITF, a key regulator of melanocyte development that is also targeted by oncogenic BMP signaling. Together, these results indicate that the oncogenic role ligand-dependent BMP signaling plays in suppressing differentiation in melanoma is a reiteration of its physiological roles during melanocyte development.

Zebrafish Pigment Pattern Formation: Insights into the Development and Evolution of Adult Form

Vertebrate pigment patterns are diverse and fascinating adult traits that allow animals to recognize conspecifics, attract mates, and avoid predators. Pigment patterns in fish are among the most amenable traits for studying the cellular basis of adult form, as the cells that produce diverse patterns are readily visible in the skin during development. The genetic basis of pigment pattern development has been most studied in the zebrafish, Danio rerio. Zebrafish adults have alternating dark and light horizontal stripes, resulting from the precise arrangement of three main classes of pigment cells: black melanophores, yellow xanthophores, and iridescent iridophores. The coordination of adult pigment cell lineage specification and differentiation with specific cellular interactions and morphogenetic behaviors is necessary for stripe development. Besides providing a nice example of pattern formation responsible for an adult trait of zebrafish, stripe-forming mechanisms also provide a conceptual framework for posing testable hypotheses about pattern diversification more broadly. Here, we summarize what is known about lineages and molecular interactions required for pattern formation in zebrafish, we review some of what is known about pattern diversification in Danio, and we speculate on how patterns in more distant teleosts may have evolved to produce a stunningly diverse array of patterns in nature.

White LED Light Exposure Inhibits the Development and Xanthophore Pigmentation of Zebrafish Embryo

Circadian rhythm in all living organisms is disturbed continuously by artificial light sources and artificial lighting has become a hazard for public health. Circadian rhythm of melatonin maintains high levels of melatonin during the night and low levels during the day. N-acetyltransferase (arylalkylamine N-acetyltransferase, AANAT) is one of the four enzymes required for melatonin synthesis and mtnr1ba is a melatonin receptor-encoding mRNA that is expressed widely in the embryonic brain. Pax7 has important roles during neural crest development and especially xanthophore pigmentation. Due to its diurnal nature, zebrafish provide a special opportunity for research on circadian rhythms that are regulated by melatonin. Here in this study, we showed that when compared with the white light control group, white LED light exposure resulted in loss of yellow pigmentation, decreased body length and locomotor activity, oxidant-antioxidant imbalance and decreased expressions of aanat2, mtnr1ba, and pax7 in zebrafish embryos. Histological analysis of this group revealed disorganization of the spaces among photoreceptor cells, decreased total retinal thickness and photoreceptor cell layer thickness compared with the control group. Artificial lighting pollution has the potential to become an important risk factor for different diseases including cancer especially for industrialized countries, therefore, more studies should be performed and necessary regulations should be made regarding this risk factor.

TRAF6 function as a novel co-regulator of Wnt3a target genes in prostate cancer

Background

Tumour necrosis factor receptor associated factor 6 (TRAF6) promotes inflammation in response to various cytokines. Aberrant Wnt3a signals promotes cancer progression through accumulation of β-Catenin. Here we investigated a potential role for TRAF6 in Wnt signaling.

Methods

TRAF6 expression was silenced by siRNA in human prostate cancer (PC3U) and human colorectal SW480 cells and by CRISPR/Cas9 in zebrafish. Several biochemical methods and analyses of mutant phenotype in zebrafish were used to analyse the function of TRAF6 in Wnt signaling.

Findings

Wnt3a-treatment promoted binding of TRAF6 to the Wnt co-receptors LRP5/LRP6 in PC3U and LNCaP cells in vitro. TRAF6 positively regulated mRNA expression of β-Catenin and subsequent activation of Wnt target genes in PC3U cells. Wnt3a-induced invasion of PC3U and SW480 cells were significantly reduced when TRAF6 was silenced by siRNA. Database analysis revealed a correlation between TRAF6 mRNA and Wnt target genes in patients with prostate cancer, and high expression of LRP5, TRAF6 and c-Myc correlated with poor prognosis. By using CRISPR/Cas9 to silence TRAF6 in zebrafish, we confirm TRAF6 as a key molecule in Wnt3a signaling for expression of Wnt target genes.

Interpretation

We identify TRAF6 as an important component in Wnt3a signaling to promote activation of Wnt target genes, a finding important for understanding mechanisms driving prostate cancer progression.

Fund

KAW 2012.0090, CAN 2017/544, Swedish Medical Research Council (2016-02513), Prostatacancerförbundet, Konung Gustaf V:s Frimurarestiftelse and Cancerforskningsfonden Norrland. The funders did not play a role in manuscript design, data collection, data analysis, interpretation nor writing of the manuscript.

•

TRAF6 positively regulated mRNA expression of b-Catenin and subsequent activation of Wnt target genes in prostate cancer cells in vitro.

•

High expression of LRP5, TRAF6 and c-Myc correlated with poor prognosis for patients with prostate cancer.

Thyroid hormone regulates distinct paths to maturation in pigment cell lineages

Thyroid hormone (TH) regulates diverse developmental events and can drive disparate cellular outcomes. In zebrafish, TH has opposite effects on neural crest derived pigment cells of the adult stripe pattern, limiting melanophore population expansion, yet increasing yellow/orange xanthophore numbers. To learn how TH elicits seemingly opposite responses in cells having a common embryological origin, we analyzed individual transcriptomes from thousands of neural crest-derived cells, reconstructed developmental trajectories, identified pigment cell-lineage specific responses to TH, and assessed roles for TH receptors. We show that TH promotes maturation of both cell types but in distinct ways. In melanophores, TH drives terminal differentiation, limiting final cell numbers. In xanthophores, TH promotes accumulation of orange carotenoids, making the cells visible. TH receptors act primarily to repress these programs when TH is limiting. Our findings show how a single endocrine factor integrates very different cellular activities during the generation of adult form.

Dynamic regulation of mRNA and miRNA associated with the developmental stages of skin pigmentation in Japanese ornamental carp

The Japanese ornamental carp (Cyprinus carpio var. Koi) is famous for multifarious colors and patterns, making it commonly culture and trade across the world. Although functional genes and inheritance of color traits have been commonly studied, seldom attentions were focused on the genetic regulation during the developmental process of pigmentation. To better understand the mechanism of skin color development, we observed the morphogenesis of pigment cells during the post-embryonic stages and analysed the temporal expression pattern of mRNAs/miRNAs profiles in four distinct developmental stages. 59 and 103 differentially expressed genes/miRNAs (DEGs/DEMs) associated with pigmentation and skin were identified, including pax7, mitf, tyr, tyrp1, etc., and the highest DEGs were detected at 11 days post hatching (dph). In addition, the functional characteristics of mRNAs/miRNAs associated with pteridine and carotenoid pathway were also examined. Furthermore, 65 miRNA-mRNA interaction pairs related to pigmentation, pteridines and carotenoids metabolism were detected between different stages. Interestingly, the largest pairs appeared in the transition from 11 dph to 48 dph, which had the similar trend with DEGs further manifesting the importance of 11 dph. This study produced a comprehensive programme of DEGs/DEMs during color development, which will provide resources to understand the regulation mechanism in color formation. The understanding of genetic basis in color formation might promote the production and breeding of the Koi carp.

The axolotl genome and the evolution of key tissue formation regulators

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl (Ambystoma mexicanum) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3. However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3^-/- and Pax7^-/- mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies.

Generation and characterization of Kctd15 mutations in zebrafish

Potassium channel tetramerization domain containing 15 (Kctd15) was previously found to have a role in early neural crest (NC) patterning, specifically delimiting the region where NC markers are expressed via repression of transcription factor AP-2a and inhibition of Wnt signaling. We used transcription activator-like effector nucleases (TALENs) to generate null mutations in zebrafish kctd15a and kctd15b paralogs to study the in vivo role of Kctd15. We found that while deletions producing frame-shift mutations in each paralog showed no apparent phenotype, kctd15a/b double mutant zebrafish are smaller in size and show several phenotypes including some affecting the NC, such as expansion of the early NC domain, increased pigmentation, and craniofacial defects. Both melanophore and xanthophore pigment cell numbers and early markers are up-regulated in the double mutants. While we find no embryonic craniofacial defects, adult mutants have a deformed maxillary segment and missing barbels. By confocal imaging of mutant larval brains we found that the torus lateralis (TLa), a region implicated in gustatory networks in other fish, is absent. Ablation of this brain tissue in wild type larvae mimics some aspects of the mutant growth phenotype. Thus kctd15 mutants show deficits in the development of both neural crest derivatives, and specific regions within the central nervous system, leading to a strong reduction in normal growth rates.

Satellite-like cells contribute to pax7-dependent skeletal muscle repair in adult zebrafish

Satellite cells, also known as muscle stem cells, are responsible for skeletal muscle growth and repair in mammals. Pax7 and Pax3 transcription factors are established satellite cell markers required for muscle development and regeneration, and there is great interest in identifying additional factors that regulate satellite cell proliferation, differentiation, and/or skeletal muscle regeneration. Due to the powerful regenerative capacity of many zebrafish tissues, even in adults, we are exploring the regenerative potential of adult zebrafish skeletal muscle. Here, we show that adult zebrafish skeletal muscle contains cells similar to mammalian satellite cells. Adult zebrafish satellite-like cells have dense heterochromatin, express Pax7 and Pax3, proliferate in response to injury, and show peak myogenic responses 4-5 days post-injury (dpi). Furthermore, using a pax7a-driven GFP reporter, we present evidence implicating satellite-like cells as a possible source of new muscle. In lieu of central nucleation, which distinguishes regenerating myofibers in mammals, we describe several characteristics that robustly identify newly-forming myofibers from surrounding fibers in injured adult zebrafish muscle. These characteristics include partially overlapping expression in satellite-like cells and regenerating myofibers of two RNA-binding proteins Rbfox2 and Rbfoxl1, known to regulate embryonic muscle development and function. Finally, by analyzing pax7a; pax7b double mutant zebrafish, we show that Pax7 is required for adult skeletal muscle repair, as it is in the mouse.

Mechanisms of silver_nanoparticles induced hypopigmentation in embryonic zebrafish

Silver_nanoparticles (AgNPs) have been reported to inhibit specification of erythroid cells and to induce spinal cord deformities and cardiac arrhythmia in vertebrates, but have not been implicated in development of neural crest (NC) and pigment cells in an in vivo model yet. In current study, down-regulated expressions of NC genes pax7 and foxd3, melanophore genes mitfa and dct, and xanthophore gene gch2 in AgNPs-exposed embryos were revealed by microarray, qRT-PCR and whole-mount in situ hybridization (WISH). Then, the down-regulated expressions of melanophore genes mitfa and dct but not xanthophore gene gch2 in AgNPs-exposed embryos were found to be recovered by melanogenesis agonists palmitic acid and dibutyryl cyclic AMP (dbcAMP). Finally, Ag⁺ chelating and AgNPs coating compound l-cysteine was found to neutralize AgNPs-induced hypopigmentation in AgNPs-exposed embryos, and to recover the down-regulated expressions of both dct and gch2 to nearly normal level in embryos, suggesting that AgNPs-releasing Ag⁺ might mediate their biological effects on zebrafish pigmentation mostly. This study was firstly to unveil that AgNPs might specifically act up-stream of mitfa and pax7 genes to suppress specification and differentiation of melanophore and xanthophore lineages respectively by their releasing Ag⁺ during vertebrate embryogenesis.

An allelic series at pax7a is associated with colour polymorphism diversity in Lake Malawi cichlid fish

Despite long-standing interest in the evolution and maintenance of discrete phenotypic polymorphisms, the molecular genetic basis of such polymorphism in the wild is largely unknown. Female sex-associated blotched colour polymorphisms found in cichlids of Lake Malawi, East Africa, represent a highly successful polymorphic phenotype, found and maintained in four genera across the geographic expanse of the lake. Previously, we identified an association with an allelic variant of the paired-box transcription factor gene pax7a and blotched colour morphs in Lake Malawi cichlid fishes. Although a diverse range of blotched phenotypes are present in Lake Malawi cichlid species, they all appeared to result from an allele of pax7a that produces increased levels of transcript. Here, we examine the developmental and genetic basis of variation among blotched morphs. First, we confirm that pax7a-associated blotch morphs result primarily from modulation of melanophore development and survival. From laboratory crosses and natural population studies, we identify at least three alleles of pax7a associated with discrete subtypes of blotched morphs, in addition to the ancestral pax7a allele. Genotypes at pax7a support initial evolution of a novel pax7a allele to produce the blotched class of morphs, followed by subsequent evolution of that pax7a blotched allele to produce additional alleles associated with discrete colour morphs. Variant alleles of pax7a produce different levels of pax7a transcript, correlating with pigmentation phenotype at the cellular level. This naturally selected allelic series should serve as a case study for understanding the molecular genetic control of pax7a expression and the evolution of sex-associated alleles.