- Research article
- Open Access
Evolution of pigment synthesis pathways by gene and genome duplication in fish
BMC Evolutionary Biology volume 7, Article number: 74 (2007)
Coloration and color patterning belong to the most diverse phenotypic traits in animals. Particularly, teleost fishes possess more pigment cell types than any other group of vertebrates. As the result of an ancient fish-specific genome duplication (FSGD), teleost genomes might contain more copies of genes involved in pigment cell development than tetrapods. No systematic genomic inventory allowing to test this hypothesis has been drawn up so far for pigmentation genes in fish, and almost nothing is known about the evolution of these genes in different fish lineages.
Using a comparative genomic approach including phylogenetic reconstructions and synteny analyses, we have studied two major pigment synthesis pathways in teleost fish, the melanin and the pteridine pathways, with respect to different types of gene duplication. Genes encoding three of the four enzymes involved in the synthesis of melanin from tyrosine have been retained as duplicates after the FSGD. In the pteridine pathway, two cases of duplicated genes originating from the FSGD as well as several lineage-specific gene duplications were observed. In both pathways, genes encoding the rate-limiting enzymes, tyrosinase and GTP-cyclohydrolase I (GchI), have additional paralogs in teleosts compared to tetrapods, which have been generated by different modes of duplication. We have also observed a previously unrecognized diversity of gchI genes in vertebrates. In addition, we have found evidence for divergent resolution of duplicated pigmentation genes, i.e., differential gene loss in divergent teleost lineages, particularly in the tyrosinase gene family.
Mainly due to the FSGD, teleost fishes apparently have a greater repertoire of pigment synthesis genes than any other vertebrate group. Our results support an important role of the FSGD and other types of duplication in the evolution of pigmentation in fish.
Coloration and color patterning of skin, scales, feathers, and hair belong to the most diverse phenotypic traits in vertebrates and have a plethora of functions such as camouflage, warning or threatening of predators, and species recognition [1, 2]. Coloration is the result of diverse pigments synthesized by pigment cells or chromatophores, which are derived from the neural crest. There are noticeable differences in the number of chromatophore types among vertebrate groups. Mammals and birds possess only the brown to black melanocytes, while amphibians and reptiles additionally have the yellow to red xantho-/erythrophores and the reflecting iridophores. In teleost fish, up to five different pigment cell types have been identified, with white leucophores and blue cyanophores in addition to the aforementioned cell types (reviewed in ). Some pigment cell types in teleosts are even further partitioned into distinct sublineages that are under different genetic control [3, 4].
The genetic basis of pigment cell development and differentiation is largely conserved between mammals and teleosts. Many genes such as Sox10, Mitf, Kit and Ednrb, some of them first identified through the cloning of coat color mutations in mice, have subsequently been found to be involved in pigmentation in teleost fish as well [5–8]. Other genes with functions in pigmentation like slc24a5 were identified first in teleosts and later on in mammals . However, an important difference between teleost fish and tetrapods has recently emerged from several studies on particular fish species. For some single copy pigmentation genes of tetrapods, two paralogous genes are present in teleost genomes, possibly as the result of a fish-specific whole-genome duplication (FSGD) that occurred ~250 to 350 million years ago (mya) in a common ancestor of teleosts (reviewed in [10–13]). Examples of such duplicated genes include sox10a and sox10b , mitfa and mitfb [15, 16], kita and kitb , csf1ra and csf1rb  and pomca and pomcb , for which at least one of the duplicates has been shown to participate in pigment cell development in fish. These genes encode transcription factors (sox10, mitf), signaling molecules (pomc) or cell-surface receptors (kit, csf1r) and are involved in neural crest specification (sox10) or commitment of pigment cell precursors to a particular chromatophore fate (mitf, kit: melanophores; csf1r: xanthophores).
A major step in chromatophore differentiation is the biosynthesis of the pigment displayed by the respective type of pigment cells. Although there are sporadic reports of duplicated genes for pigment synthesis enzymes in specific teleost lineages [20–22], no systematic genomic analysis has been performed so far to determine the complete set of duplicated pigmentation genes in fish and to better understand how pigment synthesis pathways as a whole have been affected by the FSGD.
In the present studies, we have analyzed genes involved in the biosynthesis of the dark pigment melanin, which is produced by melanophores [23, 24], and of the pteridine pigments synthesized in xanthophores (reviewed in ). We find that the FSGD had a deep impact on the melanin synthesis pathway, with three out of four enzyme-encoding genes being duplicates in teleosts. The pteridine synthesis pathway has been affected to a lesser degree by the FSGD, with two of nine enzymes represented by two teleost-specific paralogs. Several cases of lineage-specific duplication were also observed in the pteridine pathway. In both pathways, genes encoding the rate-limiting enzymes are duplicated in teleosts compared to tetrapods, with different modes of duplication being involved.
Phylogenetic analyses and synteny studies were applied to the reconstruction of the evolution of genes involved in pigmentation pathways in ray-finned fish. The combination of these complementary approaches is particularly helpful to deduce the evolutionary history of gene families following gene and genome duplication events (reviewed in ). Most genes duplicated during large-scale and in particular whole genome duplications are part of so-called paralogons, i.e., chromosomal blocks of two or more duplicated (paralogous) genes ( and references therein). The detection of such paralogons helps to draw conclusions on the origin of duplicated genes when phylogenetic signals are not informative. In the present study, sequence information from the genome assemblies of five teleost species – two pufferfish species (Tetraodon nigroviridis, Takifugu rubripes), medaka (Oryzias latipes), three-spined stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) – and from diverse teleost expressed sequence tag (EST) resources have been compared to data from tetrapods (human, mouse, chicken, frog) and different invertebrate outgroups (urochordates, sea urchin, fruitfly, nematodes). The results of these surveys are summarized in Table 1; accession numbers are given in the Additional files [see Additional files 1, 2].
Gene duplications in the melanin synthesis pathway
Melanin, the dark pigment of melanophores, is synthesized from tyrosine within a specialized organelle, the melanosome [23, 24]. In mammals and birds, two types of melanin are produced, the black to brown eumelanin and the lighter pheomelanin, but in teleost fishes only eumelanin has been observed . The eumelanin synthesis pathway is presented in Figure 1. Disruption of melanogenesis leads to reduced pigmentation intensity culminating in complete albinism .
Tyrosinase gene family
Vertebrate melanin synthesis involves the members of the tyrosinase gene family: tyrosinase (tyr), tyrosinase-related protein 1 (tyrp1) and dopachrome tautomerase (dct; also known as tyrosinase-related protein 2) [23, 24]. Tyrosinase (EC 22.214.171.124) promotes the first two rate-limiting steps of melanin synthesis from tyrosine to DOPA and DOPAquinone as well as two later steps. Dct (EC 126.96.36.199) converts DOPAchrome to DHICA, and Tyrp1 is involved in the formation of indole-5,6-quinone carboxylic acid from DHICA (Figure 1). During the early evolution of the chordate lineage, an ancestral tyrosinase gene was duplicated before the divergence of urochordates and vertebrates leading to tyrosinase and a tyrosinase-related gene. The latter one was subsequently duplicated in the vertebrate lineage giving rise to tyrp1 and dct [28, 29].
The phylogeny of the tyrosinase gene family in vertebrates is presented in Figure 2a. For tyr and tyrp1, duplications were observed in teleosts, while dct is present as a single copy in teleosts just like in tetrapods. As shown previously, tyr was duplicated in teleosts during the course of the FSGD after divergence from the more basal actinopterygian lineages of sturgeon (Acipenser baerii) and gar (Lepisosteus platyrhynchus) leading to two tyr copies in Takifugu and cichlids . Consistently, two tyr paralogs were found in pufferfishes, stickleback, medaka and rainbow trout. In the zebrafish, however, only the tyra paralog was detected, suggesting that tyrb has been lost. The two tyr paralogs previously reported for the tetraploid rainbow trout (Oncorhynchus mykiss)  are – according to the present phylogeny – the result of the FSGD rather than of another round of genome duplication in the ancestor of salmonid fishes [30, 31]. Furthermore, synteny data also provide strong evidence for the duplication of tyr during the course of the FSGD. In human, TYR is found on chromosome 11q14-q21 and the (co-)orthologs of nearby genes are also found in close proximity to the tyr paralogs in pufferfishes, stickleback, medaka and zebrafish (Figure 2b). Most of these genes are present on one of the two paralogous chromosomal blocks (paralogons) in teleost genomes, but two genes, grm5 and rab38, are also present in duplicate. The presence of these paralogons in diverse teleost genomes is consistent with the duplication of tyr during the course of the FSGD. This is also in agreement with previous studies showing that the tyr gene-containing chromosomes 13 and 14 from medaka as well as chromosomes 15 and 10 from zebrafish (Figure 2b) are derived from the same protochromosome [32–34].
Furthermore, our analysis demonstrated that tyr is not the only gene found from the melanin pathway to be duplicated in fish. Two paralogs of tyrp1 were identified in medaka, zebrafish, stickleback and fathead minnow (Pimephales promelas), while only one complete tyrp1 paralog (tyrp1a) was detected in pufferfishes. In Tetraodon, additionally a region in scaffold 13631 with partial but significant sequence similarity to tyrp1b was found. However, some splice sites of this sequence are degenerated and the putative coding sequence contains a stop codon. We confirmed the presence of this stop codon by sequencing of genomic DNA [GenBank: EF183530], thus excluding the possibility of a sequencing error. Hence, this sequence represents most likely a tyrp1b pseudogene.
In the phylogeny of the entire tyrosinase gene family based on protein sequences (Figure 2a), the tree topology is not consistent with a duplication of tyrp1 during the FSGD. In contrast, a separate maximum likelihood phylogeny of vertebrate tyrp1 genes based on nucleotide sequences suggests the duplication of tyrp1b during the course of the FSGD [see Additional file 3]. This was also confirmed by synteny data (Figure 2c), which are generally considered as more reliable than molecular phylogenies to reconstruct large-scale duplication history : the region of human chromosome 9p23 containing TYRP1 is syntenic to two tyrp1-containing paralogons in medaka, stickleback and zebrafish. Accordingly, the respective medaka chromosomes 1 and 18 have been shown to contain large duplicated segments having been formed from a same protochromosome by the FSGD .
In zebrafish, the previously described paralog tyrp1b is found on chromosome 1 in the present genome assembly (Zv6) and was mapped to the corresponding linkage group (LG) 1 . The newly found tyrp1a paralog is found in Zv6 on chromosome 11, but was not genetically mapped so far. As a paralogous relationship between zebrafish chromosomes 1 and 11 has not been reported so far and since there are frequent discrepancies between mapping data of zebrafish genes and their chromosomal assignment in current genome assemblies ( and own observations), we mapped tyrp1a using the radiation hybrid panel LN54 . The tyrp1a gene was assigned not to LG 11 (as expected from Zv6 genome assembly analysis) but to LG 7 at a distance of 0.00 cR from marker Z21714 with a LOD score of 27.4. However, a paralogous LG1–LG7 relationship has also not been reported for zebrafish so far. These data suggest the presence of a newly identified paralogon in the zebrafish genome. Kitb, another pigmentation gene duplicate that has its origin in the FSGD , is found 3' of tyrp1b (Figure 2c). However, this gene is not part of the tyrp1 paralogon, as kita is found on LG 20 and not on LG 7 in zebrafish and the human KIT is found on chromosome 4 and not on chromosome 9.
A later step in melanin synthesis from indole-5,6-quinone carboxylic acid to eumelanin is presumably catalyzed by the Silver protein (also known as Pmel17) [37, 38] (Figure 1). A recent study identified two paralogs of silver (silv) in the zebrafish . We reconstructed the evolution of the silv gene in the vertebrate lineage and additionally found two paralogs of silv in both pufferfishes as well as in medaka and stickleback, with a tree topology consistent with a duplication in the course of the FSGD (Figure 3a). This conclusion is once again strongly supported by synteny data (Figure 3b): human SILV is found on chromosome 12q13-q14, a region that is highly syntenic to Tetraodon chromosomes 9 and 11, medaka chromosomes 5 and 7 and zebrafish chromosomes 11 and 23, which contain the silv paralogons. As all these chromosomes evolved from a same protochromosome through duplication during the FSGD [32–34, 39], it is most likely that the silv duplicates in teleosts stem from this event.
Three non-related genes, oca2, aim1 and slc24a5, encode for transporter proteins residing in the melanosomal membrane and being essential for melanin synthesis (Figure 1). Loss-of-functions mutations in these genes lead to reduced melanin pigmentation in teleosts and mammals [9, 40–42]. In contrast to tyr, tyrp1 and silv, all three transporters are encoded by a single gene in the teleost species analyzed [see Additional file 4].
Gene duplications in the pteridine synthesis pathway
The yellow to reddish pteridine pigments of xanthophores are synthesized from GTP through the pteridine synthesis pathway (Figure 4) (reviewed in ). The pteridine pathway is composed of three component pathways. The first one leads to the de novo synthesis of tetrahydrobiopterin (H4biopterin) from GTP. H4biopterin is a cofactor for neurotransmitter synthesis and tyrosinase activity in melanophores . The second component is the regeneration pathway of oxidized H4biopterin after it has acted as cofactor . The third pathway shares several steps with the first one and leads to the formation of the yellow pigments, sepiapterin and its derivatives, as well as probably to the reddish drosopterin, which is also found in teleost fishes, especially in poecilids .
GTP cyclohydrolase I and its feedback regulatory protein
The first, rate limiting step in pteridine synthesis is catalyzed by the GTP cyclohydrolase I (GchI; EC 188.8.131.52) (Figure 4). GchI expression is an initial step for melanophore and xanthophore differentiation due to its involvement in the different component pathways . Only one gchI gene has been reported so far in vertebrates, including the teleosts rainbow trout  and zebrafish . The present survey for gchI sequences in vertebrates, however, revealed an unexpected diversity of gchI genes in teleosts and amphibians (Figure 5). While a single gchI gene was found in mammals and birds, two gchI genes were identified in frog, medaka, zebrafish, and fathead minnow, and even three genes are present in pufferfishes and stickleback. Taking conserved syntenies into account, three groups of gchI genes became obvious.
The first group, termed gchIa, is phylogenetically well defined and includes the single gchI gene from mammals and birds and one copy from frog and teleosts including the known rainbow trout gene. All vertebrate gchIa orthologs are found within a region of conserved synteny. In rainbow trout, even two copies of gchIa are found, which might be the result of the salmonid autotetraploidization event that occurred 25 to 100 mya [30, 31]. Paralogs of other pigment synthesis enzymes in trout and/or salmon (silv, gchfr, dhprb, pam; Figure 6,Table 1, [see Additional files 1, 2]) might also have resulted from this salmonid-specific genome duplication.
The second group of gchI genes, gchIb, consists of the second gchI gene from frog, the previously known zebrafish gene and further teleost genes. The orthology of gchIb genes is well supported by conserved syntenies between frog and teleosts. The third group, gchIc, has been found so far only in pufferfishes, stickleback and the gilthead seabream (Sparus aurata) and is also phylogenetically and syntenically well defined.
How are the three gchI groups related to each other? We confirmed by RHP mapping the chromosomal allocations of gchIa and gchIb in the zebrafish genome assembly on chromosomes 17 and 12, respectively. GchIa was assigned to LG 17 at a distance of 0.00 cR from marker fc19b04 with a LOD score of 18.4. GchIb was mapped to LG 12 at a distance of 0.00 cR from marker fc18g04 with a LOD score of 15.3. As LG 17 and LG 12 do not seem to have evolved by protochromosome duplication during the FSGD  and due to the presence of gchIb in amphibians, the duplication that led to gchIa and gchIb seems to be older than the split between ray-finned fishes and tetrapods. Both gchIa and gchIb genes are found in proximity to members of the socs gene family: gchIa is linked to socs4 in all vertebrates examined, gchIb to socs5 in frog and zebrafish (Figure 5a). Socs4 and socs5 are the closest related members within the socs gene family . Therefore it seems most likely that gchIa/b and socs4/5 precursors were duplicated together, possibly during one of the two earlier rounds of genome duplication having taken place during the early evolution of the vertebrate lineage (1R or 2R) [30, 47, 48]. Socs5 is also found in mammals and birds within a syntenic region that resembles the gchIb region of frog and teleosts (Figure 5b) suggesting that gchIb was lost secondarily in these lineages. The human regions containing GCHIA/SOCS4 and SOCS5 are found on chromosomes 14 (Figure 5a) and 2 (Figure 5b), respectively, which were shown to contain many paralogous genes that arose during the 1R/2R genome duplications .
The origin of gchIc found in pufferfishes, stickleback and gilthead seabream remains unclear. Genes surrounding gchIc are not related to those of the other gchI regions and the human orthologs of these genes are found on chromosome 19q13. The corresponding chromosomal region on medaka chromosome 4 seems to be highly conserved in gene order with pufferfishes and stickleback, but a large gap is found in the medaka genome assembly between pik3r2 and ankrd47 (not shown). Thus, gchIc might also be present in medaka but absent from the current genome assembly. However, no EST or shotgun trace sequence from medaka was found that could represent gchIc. In zebrafish, a less conserved chromosomal block is found on chromosome 2. If gchIc genes arose as a paralog of gchIa or gchIb as result of the FSGD, one would expect to find gchIc on other chromosomes (Tetraodon: chr 14 or 3; medaka: chr 24 or 8; zebrafish: chr 20 or 3) [32–34]. Thus, there is no evidence that gchIc has been formed during the FSGD and its relationships to the other gchI groups remain elusive. It might be possible that gchIc arose by a lineage-specific gene duplication or that it is also a remnant of earlier rounds of genome duplication in vertebrates that has been maintained only in some teleost lineages.
The GchI enzymatic activity is regulated by the H4biopterin-dependent GTP cyclohydrolase I feedback regulatory protein (Gchfr) . In most teleost species, a single gchfr gene was found. In contrast, two gchfr genes were identified in rainbow trout and Atlantic salmon (Salmo salar) (Figure 6). The phylogeny suggests duplication of gchfr in the common ancestor of these salmonid fishes, which fits well the salmonid autotetraploidization event [30, 31].
6-pyruvoyltetrahydropterin synthase and sepiapterin reductase
Subsequent steps of pteridine synthesis are catalyzed by the 6-pyruvoyltetrahydropterin synthase (Pts; EC 184.108.40.206) and the sepiapterin reductase (Spr; EC 220.127.116.11) (Figure 4). In the guppy, pts expression correlates with the presence of xanthophore-based yellow color patterns . A single pts gene was found in all vertebrates analyzed including teleosts [see Additional file 5a].
Sepiapterin reductase (Spr) catalyzes the next three steps towards H4biopterin (as well as a step in the third component pathway; see below). We could identify only one spr gene in tetrapods as well as in medaka and Tetraodon. In zebrafish, stickleback and Takifugu, however, two spr genes were found (Figure 7). Although the relationships between spr genes are phylogenetically not fully resolved, the genomic regions containing the duplicated teleost spr genes are syntenic to each other as well as to human chromosome 2, where SPR is found (Figure 7b). Both teleost spr genes are in close proximity to paralogs of smyd1. Additional duplicated genes are also found in these chromosomal regions. Furthermore, zebrafish chromosomes 5 and 8, which contain spra and sprb, respectively, evolved from the same pre-FSGD protochromosome [33, 34]. It is therefore most likely that the spr duplicates found in some teleost lineages are remnants of the FSGD.
Enzymes of the H4biopterin regeneration pathway
The H4biopterin regeneration pathway involves the enzymes Pcbd (pterin-4 alpha-carbinolamine dehydratase/dimerization cofactor of hepatocyte nuclear factor 1 alpha (TCF1); EC 18.104.22.168) and Dhpr (Dihydropteridine reductase; EC 22.214.171.124) . Pcbd is a bifunctional protein having a function as a dimerizing co-factor of the HNF1 homeobox transcription factors in addition to its enzymatic activity. A transcriptional target of the Pcbd/HNF1 complex is the tyrosinase promoter, pointing to a particular importance for melanophore differentiation in human . In contrast to invertebrates, two pcbd genes, pcbd1 and pcbd2, were found in tetrapods and teleosts (Figure 8a), which is consistent with a previous analysis of this gene family . Thus, pcbd was duplicated early in the vertebrate lineage, possibly during one of the two rounds of genome duplication (1R/2R). Two copies of pcbd1 were identified in Takifugu. While one copy (located in scaffold 53) has the usual structure with four coding exons, the other copy (scaffold 178) consists of a single exon (Figure 8b). A polyA sequence downstream of the stop codon and target site duplications indicative of a retrotransposition event were detected (not shown). Furthermore, a mutation in pcbd1 from scaffold 178 turned codon 31 into a premature stop codon (Figure 8b). Thus, this pcbd1 copy most likely represents a "processed" pseudogene generated by retrotransposition of a pcbd1 mRNA.
Several types of duplications have affected dhpr genes in teleosts (Figure 9). First of all, while a single dhpr gene is found in tetrapods as well as in pufferfishes, medaka and stickleback, two dhpr copies are present in salmon, trout and fathead minnow. The zebrafish genome contains even three copies of dhpr. Two major clades of dhpr genes become apparent in teleosts through phylogenetic analysis. The first clade, termed dhpra, consists of most of the teleost sequences including pufferfishes, medaka, stickleback, the zebrafish sequence found on chromosome 14 and one copy each of salmon, trout and fathead minnow. The second clade, dhprb, contains the second copy from salmon, trout and fathead minnow as well as the two zebrafish sequences found on chromosome 1 in scaffolds 64 and 67, respectively. These two zebrafish paralogs on chromosome 1 obviously arose by duplication in the family Cyprinidae and encode putative proteins that show 82% sequence similarity. Both copies of dhprb, termed dhprba (scaffold 64) and dhprbb (scaffold 67), could be amplified by PCR and sequenced from zebrafish cDNA with paralog-specific primer sets [GenBank: EF183528, EF183529], excluding genome assembly artifacts. In the zebrafish genome assembly, additional but partial sequences of dhprba and dhprbb are present in scaffold 63 and scaffold 67, respectively. These sequences were not included in further analyses.
The analysis of the dhpr-containing regions in vertebrate genomes revealed that the two main clades, dhpra and dhprb, might originate from the FSGD (Figure 9b): Genes surrounding the human DHPR gene on chromosome 4 are found in vicinity to the teleost dhpra gene (Tetraodon: chr 20; medaka: chr 10: zebrafish: chr 14) as well as on another chromosome (Tetraodon: chr 18; medaka: chr 1; zebrafish: chr 1). All these chromosomes evolved by duplication of the same protochromosome during the course of the FSGD . Later on, dhprb was further duplicated in the lineage leading to zebrafish probably through intrachromosomal gene duplication. This led to the formation of dhprba and dhprbb on chromosome 1, where they are separated by approximately 1 Mb.
Enzymes involved in pteridine pigment synthesis
The third component pathway that leads to the formation of the yellow pteridine sepiapterin and its derivatives branches off from the first component pathway by hypothetical enzymatic reactions (Figure 4). Subsequent reactions require Spr (see above) and Xod/Xdh (xanthine oxidase/xanthine dehydrogenase; EC 126.96.36.199/EC 188.8.131.52) . As in tetrapods, Xod/Xdh is represented by a single gene in teleost genomes [see Additional file 5b].
The biosynthetic pathway for the reddish drosopterin has not been elucidated yet in vertebrates and only one enzyme of the pathway in Drosophila, clot, has been characterized at the molecular level . A single thioredoxin-like 5 gene, the vertebrate ortholog of Drosophila clot, is found in tetrapods and teleosts as well [see Additional file 5c].
Finally, the switch between the H4biopterin and sepiapterin synthesis might be regulated by PAM (protein associated with Myc), which is affected in the zebrafish esrom mutant that has reduced yellow pigmentation . The pam gene is single-copy in teleosts and tetrapods [see Additional file 5d].
Duplication of pigmentation genes: molecular mechanisms and evolutionary fates
In the present study, we have analyzed the two major pigment synthesis pathways in vertebrates, the melanin and the pteridine pathways, with respect to gene and genome duplications particularly within the teleost lineage. Seventeen vertebrate pigmentation genes were analyzed and various modes of duplication were observed. On the one hand, different rounds of genome duplication have expanded several pigment gene families. Five clear cases of FSGD-based duplications (tyr, tyrp1, silv, spr, dhpr) were found (29%). Other duplications might be the result of earlier rounds of genome duplication (1R/2R) [30, 47] in the vertebrate stem lineage (gchIa/b, pcbd1/2). In addition, gene duplications generated by the recent salmonid-specific autotetraploidization [30, 31] could be also detected (silv, gchIa, gchfr, dhprb, pam). On the other hand, lineage-specific local gene duplications were also identified: the duplication of dhprb in the zebrafish, the duplication by retrotransposition of pcbd1 in Takifugu and possibly the occurrence of gchIc in a common ancestor of pufferfishes, stickleback and perciforms. Although the majority of duplicated genes in vertebrate genomes were created by whole genome duplications , lineage-specific duplications of pigmentation genes, which have also been found for the urochordate Ciona intestinalis , seem to be a common theme in chordate evolution. In conclusion, teleost fishes have a greater potential repertoire of pigment synthesis genes than all other vertebrate groups. However, entirely duplicated synthesis pathways are not observed, and the function of both paralogs in pigmentation pathways remains generally to be demonstrated.
The impact of genome duplications on entire metabolic pathways in the vertebrate lineage has been studied so far only for the glycolysis . Based on a similar approach to that used in the present study, the authors showed that none of the three rounds of genome duplication in the vertebrate lineage (1R/2R/FSGD) led to a completely duplicated glycolytic pathway in extant genomes. In total, 46% of the glycolytic enzymes in vertebrates were duplicated in teleosts due to the FSGD (11 of 24 enzymes) . Here, 75% (3/4) of the melanogenic enzymes and 22% (2/9) of the enzymes from the pteridine pathway were found to be duplicated during the FSGD. Although the value for melanogenesis seems to be elevated in comparison to pteridine synthesis and glycolysis, all differences between the three pathways (glycolysis, melanogenesis and pteridine pathway) are statistically not significant (χ2-test, p > 0,05).
Generally, three different fates of duplicated genes are observed (reviewed by ). In most cases, one duplicate gets lost due to functional redundancy. This process of non-functionalization was estimated to have occurred in a range of 76% of FSGD duplicates in zebrafish  and 76 to 85% in the pufferfish lineage [33, 39, 58]. Here, for pigmentation genes 71% (12/17) were found to be reduced from two to one copy in teleosts after the FSGD but before the split of Ostariophysii (zebrafish) and Neoteleostei (medaka, stickleback and pufferfishes). Including lineage-specific losses, the ratio of non-functionalization for pigment synthesis genes is 82% (14/17) for pufferfish and medaka and 76% (13/17) for stickleback and zebrafish suggesting that pigment synthesis genes do not deviate from the global trend. Two other fates of gene duplicates might lead to the retention of both copies within a genome. Either one copy obtains a new function (neo-functionalization) or the original gene functions are divided between the two duplicates (sub-functionalization). Recently, it was shown that combinations of both mechanisms are possible (sub-neo-functionalization) . Asymmetric evolution, which might be an indicator for neo-functionalization, has been observed for many duplicated genes in teleosts [58, 60, 61] including pigmentation genes . Neo-functionalization of duplicated enzymes can lead, for example, to the evolution of new substrate specificities  or even of entirely new functions not associated with the enzymatic property . Subfunctionalization of duplicated enzymes might occur at the level of gene expression leading, e.g., to tissue-specific expression ( and references therein) or at the protein level, when a duplicate becomes specialized for a certain substrate . Whether and how functional divergence of duplicated pigment synthesis enzymes has occurred in teleosts will be an important focus of future studies.
Evolution of the melanin synthesis pathway
The melanin synthesis pathway involves four enzymes. Three of them were found to be duplicated in teleosts as result of the FSGD. In the tyrosinase gene family, FSGD-duplication was observed for tyr and tyrp1, while dct was present as a single copy gene in all lineages analyzed. However, the retention of tyrosinase gene family members after the FSGD is variable between the different lineages (Table 1). Tyra was lost in the zebrafish and tyrp1b in the pufferfishes, while medaka and stickleback have retained both copies of tyr and tyrp1. Thus, the tyrosinase gene family is a good example for divergent resolution, i.e., differential loss of gene duplicates in divergent lineages, a mechanism that might facilitate speciation [13, 65–67].
Mutational disruption of melanin synthesis at different steps of the pathway might lead to diverse forms of albinism . Tyrosinase is the first, rate-limiting enzyme of melanogenesis. In the zebrafish, loss-of-function in the single tyr gene, tyra, leads to an albino phenotype . In the medaka, several albino mutants were identified that are also affected in the tyra paralog . Our data provide evidence for the presence of tyrb in the medaka but the functions of this paralog in teleosts remain unresolved. No tyrb mutant is available at present in fish. The fact that some tyra mutations in the medaka lead to a complete albino phenotype  suggests that tyrb cannot substitute for tyra. This is in agreement with functional studies of the two tyr duplicates in the rainbow trout : simultaneous morpholino knock-down of both paralogs reduces pigmentation in the eye and the skin to the same amount as knocking-down tyr paralogs separately. Since knock-down of tyrb gene function in the rainbow trout leads to reduced pigmentation in the eye and the skin , tyrb seems to be involved in melanin synthesis too. Tyrosinase is involved in several steps of melanogenesis (Figure 1), and it is therefore possible that teleosts tyr paralogs might have become subfunctionalized and specialized for individual steps of the pathway.
There is so far no evidence supporting the functional divergence of tyrp1 paralogs in fish. Mutation of tyrp1 in mammals leads to reduced pigmentation . No tyrp1 mutant has been identified in teleosts until today, possibly due to a functional redundancy of tyrp1 duplicates. Interestingly, in the present study a putative regulator of Tyrp1 function was also found to be duplicated in teleosts as result of the FSGD: tyr duplicates in teleosts are genetically linked to duplicates of rab38 (Figure 2b). Rab38 is thought to play a role in sorting Tyrp1 to the melanosome in mice .
The duplication of the silver gene has been previously described in the zebrafish . Our study shows that this duplication is indeed the result of the FSGD and that silver has also been retained in duplicate in pufferfishes, medaka and stickleback. In zebrafish, silva is expressed in melanophores and the retinal pigment epithelium (RPE) of the eye, while silvb expression is restricted to the RPE . The expression of silv paralogs is similar to the expression of duplicated mitf transcription factor genes . In mammals, Silv transcription is dependent on Mitf [71, 72]. It will be highly interesting to investigate the differential regulation of silv paralogs by Mitf duplicates in different teleost lineages.
Due to the limited knowledge of gene functions it remains elusive at present, whether there is a correlation between excess of genes involved in melanin synthesis and the vast diversity of coloration in fish. Functional experiments on the divergence of pigmentation gene duplicates are currently carried out in our laboratory to elucidate this question.
Evolution of the pteridine synthesis pathway
The pteridine synthesis pathway has been less affected by the FSGD than the melanin pathway, but several cases of lineage-specific duplication were observed.
GchI is the first and rate-limiting enzyme of pteridine synthesis. In this analysis, we have observed an unforeseen diversity of gchI genes in vertebrates. We could identify two clades of gchI genes, gchIa and gchIb, which most likely arose through genome duplication during early vertebrate evolution, as well as a third clade of unresolved origin, gchIc, which is only found in some teleost species. The GchI enzyme is required at the initial step of the synthesis of both H4biopterin and pteridine pigments (Figure 4). GchIa has been found in all vertebrate lineages and is therefore most likely involved in H4bioterin formation. GchIb is only found in those lineages that possess xanthophores: teleost fishes and amphibians. Furthermore, gchIb from zebrafish, which is a paralog of the mammalian gchIa (and not its ortholog as previously thought), is expressed in the xanthophore lineage (but also in melanophores and neurons) . We therefore propose that gchIb plays a major role in the synthesis of pteridine pigments of xanthophores and that it was lost secondarily in mammals and birds concomitantly to the loss of xanthophores in these lineages. Functional studies in teleosts and amphibians will be necessary to test this hypothesis.
Spr is involved in both the de novo synthesis of H4bioterin and the production of pteridine pigment after the split between both component pathways (Figure 4). Interestingly, the spr gene is found to be duplicated as result of the FSGD in zebrafish, stickleback and Tetraodon. It might be possible that each of the spr paralogs has become specialized for one component pathway, but expression data for duplicated teleost spr genes are not available at present. Sprb paralogs might have been lost quite recently in medaka as well as in Tetraodon after its split from Takifugu. This is a good example for the former observation that anciently duplicated genes still can be lost after millions of years .
Finally, dhpr in zebrafish illustrates how different evolutionary scenarios can progressively shape pigmentation gene families. After the duplication of dhpr in the FSGD, both dhpr paralogs were retained in Ostariophysii (zebrafish, fathead minnow) until dhprb was further duplicated in the zebrafish lineage, while dhprb was apparently lost from pufferfishes, medaka and stickleback.
Evolution by genome duplication: the pigmentary system
The evolutionary significance of whole genome duplications is still widely debated. The two presumed rounds of genome duplication early in the vertebrate lineage (1R/2R) have been linked to an increase in phenotypic complexity and to the evolution of vertebrate-specific traits such as the neural crest [30, 47]. Several authors have suggested that the divergent evolution of duplicates generated by the FSGD might be involved in species diversity in teleost fishes, which represent approximately 50% of all vertebrate species (reviewed in [10–13]). However, these hypotheses have been questioned based on the fossil record . In addition, a reduced probability of extinction in teleost fishes compared to other vertebrates probably due to the FSGD has been proposed, since mutational robustness, increased genetic variation, and increased tolerance to environmental conditions could be by-products of genome duplication .
With regard to the pigmentary system, it has been previously suggested that the FSGD had a major importance for the evolution of pigmentation genes in teleost fish . The present study puts further evidence in this direction by showing that pigment synthesis pathways (and the melanin synthesis pathway in particular) have been affected by the FSGD. Interestingly, the genetic repertoire for color perception, i.e., the opsin gene family, has also been expanded by duplications in the teleost lineage [75, 76]. It remains to be elucidated whether the diversity and complexity of coloration observed in teleost fishes compared to other vertebrate groups are causally linked to the expansion of pigmentation gene families as result of the FSGD. This FSGD might also have provided the genetic raw material for the diversity of coloration within teleosts since species-specific sequence evolution of duplicated genes is a common mechanism in this group . Furthermore, our study points out lineage-specific patterns of loss and retention of duplicated pigmentation genes in teleosts. Divergent resolution of duplicated genes might facilitate speciation events [13, 65–67].
The present study shows that teleost fishes have a greater repertoire of pigment synthesis genes than any other vertebrate group mainly due to the fish-specific genome duplication but also as result of other types of gene duplications. Thus, pigmentation genes from teleosts offer an excellent opportunity to study the effects of gene and genome duplication on gene regulatory, protein-protein interaction and metabolic networks (e.g., specification of chromatophore fates, receptor-ligand interactions and pigment synthesis, respectively) and their connections. Future studies on functional divergence of duplicated pigmentation genes will reveal important insights into the significance of gene and genome duplication for the evolution of vertebrate phenotypes.
Sequence database surveys
Nucleotide sequences of pigmentation genes from ray-finned fishes were identified using BLAST searches against GenBank (nr and EST databases), the current genome assemblies and Trace Archives at Ensembl  of zebrafish (Zv 6), medaka (version 1), Tetraodon (version 7), Fugu (version 4.0), stickleback (BROAD S1) as well as TIGR gene indices  of cichlids (Astatotilapia burtoni, Haplochromis chilotes and Haplochromis sp. 'red tail sheller'), catfish (Ictalurus punctatus), killifish (Fundulus heteroclitus), rainbow trout and salmon. Usually the human gene was used as query sequence. If necessary, coding sequences were annotated manually from genome assemblies based on sequence homology to other species. In some cases species-specific EST clusters were assembled. Similarly, sequences from human, mouse, chicken, frog, ascidian (Ciona intestinalis), sea urchin (Strongylocentrotus purpuratus), fruitfly (Drosophila melanogaster) and nematode (Caenorhabditis elegans) were obtained from GenBank or Ensembl under inclusion of information given in ref. 29.
Sequence alignments and phylogenetic reconstructions
All nucleotide sequences obtained from BLAST searches were loaded into BioEdit , translated into proteins, and aligned using ClustalW  as implemented in BioEdit. Alignments were carefully checked and ambiguously aligned regions were removed prior to phylogeny analyses. Identical sequences were removed.
Larger draft neighbour-joining trees were obtained with MEGA3 . Based on these trees, outgroups for final phylogenies were chosen. These were either the closest human paralog to the gene under investigation (in case of larger vertebrate gene families) or invertebrate orthologs.
Final protein maximum likelihood phylogenies were computed with PHYML [82, 83] with 100 bootstrap replicates. Models of protein evolution and parameter values were determined with ProtTest . PAUP  was used to obtain neighbor-joining bootstrap values of 10,000 replicates.
Syntenic relationships between human and teleosts genomes within the chromosomal regions containing the gene of interest were inferred using the Reciprocal Blast Hit method .
Sequences of 15–20 genes surrounding the human ortholog were used as initial queries for BLAST searches against the five teleost genome assemblies at Ensembl , followed by reciprocal BLAST searches of the best hits against human and other teleost genomes.
Radiation hybrid panel mapping
The zebrafish radiation hybrid panel LN54  was used according to the supplier's instructions (Marc Ekker, University of Ottawa) to map tyrp1a, gchIa and gchIb. The following primer sets were used: Dre-tyrp1a-ex1F: 5'-ATGTTTGGACTTTATGGA GC-3', Dre-tyrp1a-ex1R: 5'-GTCAAACCCGCTGTAGTTC-3' (annealing temperature TA: 56°C); Dre-gchIA-ex1F: 5'-AAGAAACTGACGGAGCGATC-3', Dre-gchIA-ex1R: 5'-TCTCCTGGTATCCCTTGGTG-3' TA: 56°C); Dre-gchIB-ex1F: 5'-CAATGGCAAAATCGTCACAG-3', Dre-gchIB-ex1R: 5'-TGGTCTCGTGGTATC CCTTAG-3' (TA: 52°C). The obtained RHVECTORs were submitted to the LN54 radiation hybrid map website  to get chromosomal positions.
Sequencing of Tetraodon tyrp1b pseudogene and zebrafish dhprb genes
The Tetraodon tyrpb1 pseudogene was amplified from genomic DNA and sequenced using primers Tni-ps-trp1b-F1 (5'- AACCTGGACACAAAGCCTCAC-3') and Tni-ps-trp1b-R1 (5'-ATGGTAGGAGAGAGCACGCAC-3') (TA: 62°C).
Zebrafish (strain WüAB) total RNA was extracted from various adult tissues using the TRIzol Reagent (Invitrogen, Karlsruhe, Germany). cDNA was synthesized with the RevertAid TM First Strand cDNA Synthesis Kit (Fermentas Life Science, St. Leon-Rot, Germany) and pooled. Dhprb sequences were amplified from the cDNA pool using paralog specific primer sets: Dre-dhprba-ex1F: 5'-CTCGTGAAGACAGAATGGCAG-3', Dre-dhprba-ex7R: 5' TGCTTTCTCCAGTCGTCCAC-3' (TA: 60°C); Dre-dhprbb-ex1F: 5'-AGCGAAGTAAAGAAAGTGATTG-3', Dre-dhprbb-ex7R: 5'-TAGGGGTAG CCACTGTTCTG-3' (TA: 58°C). PCR products were cloned using the TA Cloning Kit (Invitrogen, Karlsruhe, Germany) and subsequently sequenced. Sequencing was performed with a CEQ 2000XL system (Beckman Coulter, Krefeld, Germany).
Fujii R: Coloration and chromatophores. The Physiology of Fishes. Edited by: Evans DH. 1993, Boca Raton, Florida: CRC Press, 535-562.
Bagnara JT: Comparative anatomy and physiology of pigment cells in nonmammalian tissues. The Pigmentary System: Physiology and Pathophysiology. Edited by: Nordlund JJ, Boissy RE, Hearing VJ, King RA, Ortonne JP. 1998, New York: Oxford University Press, 9-40.
Mellgren EM, Johnson SL: The evolution of morphological complexity in zebrafish stripes. Trends Genet. 2002, 18: 128-34. 10.1016/S0168-9525(01)02614-2.
Odenthal J, Rossnagel K, Haffter P, Kelsh RN, Vogelsang E, Brand M, van Eeden FJ, Furutani-Seiki M, Granato M, Hammerschmidt M, Heisenberg CP, Jiang YJ, Kane DA, Mullins MC, Nusslein-Volhard C: Mutations affecting xanthophore pigmentation in the zebrafish, Danio rerio. Development. 1996, 123: 391-8.
Dutton KA, Pauliny A, Lopes SS, Elworthy S, Carney TJ, Rauch J, Geisler R, Haffter P, Kelsh RN: Zebrafish colourless encodes sox10 and specifies non-ectomesenchymal neural crest fates. Development. 2001, 128: 4113-25.
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.
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.
Parichy DM, Mellgren EM, Rawls JF, Lopes SS, Kelsh RN, Johnson SL: 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. 10.1006/dbio.2000.9899.
Lamason RL, Mohideen MA, Mest JR, Wong AC, Norton HL, Aros MC, Jurynec MJ, Mao X, Humphreville VR, Humbert JE, Sinha S, Moore JL, Jagadeeswaran P, Zhao W, Ning G, Makalowska I, McKeigue PM, O'Donnell D, Kittles R, Parra EJ, Mangini NJ, Grunwald DJ, Shriver MD, Canfield VA, Cheng KC: SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science. 2005, 310: 1782-6. 10.1126/science.1116238.
Postlethwait J, Amores A, Cresko W, Singer A, Yan YL: Subfunction partitioning, the teleost radiation and the annotation of the human genome. Trends Genet. 2004, 20: 481-90. 10.1016/j.tig.2004.08.001.
Meyer A, Van de Peer Y: From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays. 2005, 27: 937-45. 10.1002/bies.20293.
Froschauer A, Braasch I, Volff JN: Fish genomes, comparative genomics and vertebrate evolution. Current Genomics. 2006, 1: 43-57. 10.2174/138920206776389766.
Volff JN: Genome evolution and biodiversity in teleost fish. Heredity. 2005, 94: 280-94. 10.1038/sj.hdy.6800635.
Lang M, Miyake T, Braasch I, Tinnemore D, Siegel N, Salzburger W, Amemiya CT, Meyer A: A BAC library of the East African haplochromine cichlid fish Astatotilapia burtoni. J Exp Zoolog B Mol Dev Evol. 2006, 306: 35-44. 10.1002/jez.b.21068.
Lister JA, Close J, Raible DW: Duplicate mitf genes in zebrafish: complementary expression and conservation of melanogenic potential. Dev Biol. 2001, 237: 333-44. 10.1006/dbio.2001.0379.
Altschmied J, Delfgaauw J, Wilde B, Duschl J, Bouneau L, Volff JN, Schartl M: Subfunctionalization of duplicate mitf genes associated with differential degeneration of alternative exons in fish. Genetics. 2002, 161: 259-67.
Mellgren EM, Johnson SL: kitb, a second zebrafish ortholog of mouse Kit. Dev Genes Evol . 2005, 215 (9): 470-77. 10.1007/s00427-005-0001-3.
Braasch I, Salzburger W, Meyer A: Asymmetric evolution in two fish-specifically duplicated receptor tyrosine kinase paralogons involved in teleost coloration. Mol Biol Evol. 2006, 23: 1192-202. 10.1093/molbev/msk003.
de Souza FS, Bumaschny VF, Low MJ, Rubinstein M: Subfunctionalization of expression and peptide domains following the ancient duplication of the proopiomelanocortin gene in teleost fishes. Mol Biol Evol. 2005, 22: 2417-27. 10.1093/molbev/msi236.
Hoegg S, Brinkmann H, Taylor JS, Meyer A: Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish. J Mol Evol. 2004, 59: 190-203. 10.1007/s00239-004-2613-z.
Schonthaler HB, Lampert JM, von Lintig J, Schwarz H, Geisler R, Neuhauss SC: A mutation in the silver gene leads to defects in melanosome biogenesis and alterations in the visual system in the zebrafish mutant fading vision. Dev Biol. 2005, 284: 421-36.
Boonanuntanasarn S, Yoshizaki G, Iwai K, Takeuchi T: Molecular cloning, gene expression in albino mutants and gene knockdown studies of tyrosinase mRNA in rainbow trout. Pigment Cell Res. 2004, 17: 413-21. 10.1111/j.1600-0749.2004.00166.x.
Hearing VJ, Tsukamoto K: Enzymatic control of pigmentation in mammals. FASEB J. 1991, 5: 2902-9.
del Marmol V, Beermann F: Tyrosinase and related proteins in mammalian pigmentation. FEBS Lett. 1996, 381: 165-8. 10.1016/0014-5793(96)00109-3.
Ziegler I: The pteridine pathway in zebrafish: regulation and specification during the determination of neural crest cell-fate. Pigment Cell Res. 2003, 16: 172-82. 10.1034/j.1600-0749.2003.00044.x.
Postlethwait J: The zebrafish genome: a review and msx gene case study. Genome Dyn. 2006, 2: 183-97.
Spritz RA, Chiang PW, Oiso N, Alkhateeb A: Human and mouse disorders of pigmentation. Curr Opin Genet Dev. 2003, 13: 284-9. 10.1016/S0959-437X(03)00059-5.
Sato S, Toyoda R, Katsuyama Y, Saiga H, Numakunai T, Ikeo K, Gojobori T, Yajima I, Yamamoto H: Structure and developmental expression of the ascidian TRP gene: insights into the evolution of pigment cell-specific gene expression. Dev Dyn. 1999, 215: 225-37. 10.1002/(SICI)1097-0177(199907)215:3<225::AID-AJA5>3.0.CO;2-S.
Takeuchi K, Satou Y, Yamamoto H, Satoh N: A genome-wide survey of genes for enzymes involved in pigment synthesis in an ascidian, Ciona intestinalis. Zoolog Sci. 2005, 22: 723-34. 10.2108/zsj.22.723.
Ohno S: Evolution by gene duplication. 1970, New York: Springer-Verlag
Allendorf FW, Thorgaard GH: Tetraploidy and the evolution of salmonid fishes. Evolutionary genetics of fish. Edited by: Turner JB. 1984, New York: Plenum Press, 1-53.
Naruse K, Tanaka M, Mita K, Shima A, Postlethwait J, Mitani H: A medaka gene map: the trace of ancestral vertebrate proto-chromosomes revealed by comparative gene mapping. Genome Res. 2004, 14: 820-8. 10.1101/gr.2004004.
Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS: The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res. 2005, 15: 1307-14. 10.1101/gr.4134305.
Kohn M, Hogel J, Vogel W, Minich P, Kehrer-Sawatzki H, Graves JA, Hameister H: Reconstruction of a 450-My-old ancestral vertebrate protokaryotype. Trends Genet. 2006, 22: 203-10. 10.1016/j.tig.2006.02.008.
Guryev V, Koudijs MJ, Berezikov E, Johnson SL, Plasterk RH, van Eeden FJ, Cuppen E: Genetic variation in the zebrafish. Genome Res. 2006, 16: 491-7. 10.1101/gr.4791006.
Hukriede N, Fisher D, Epstein J, Joly L, Tellis P, Zhou Y, Barbazuk B, Cox K, Fenton-Noriega L, Hersey C, Miles J, Sheng X, Song A, Waterman R, Johnson SL, Dawid IB, Chevrette M, Zon LI, McPherson J, Ekker M: The LN54 radiation hybrid map of zebrafish expressed sequences. Genome Res. 2001, 11: 2127-32. 10.1101/gr.210601.
Chakraborty AK, Platt JT, Kim KK, Kwon BS, Bennett DC, Pawelek JM: Polymerization of 5,6-dihydroxyindole-2-carboxylic acid to melanin by the pmel 17/silver locus protein. Eur J Biochem. 1996, 236: 180-8. 10.1111/j.1432-1033.1996.t01-1-00180.x.
Theos AC, Truschel ST, Raposo G, Marks MS: The Silver locus product Pmel17/gp100/Silv/ME20: controversial in name and in function. Pigment Cell Res. 2005, 18: 322-36. 10.1111/j.1600-0749.2005.00269.x.
Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biemont C, Skalli Z, Cattolico L, Poulain J, De Berardinis V, Cruaud C, Duprat S, Brottier P, Coutanceau JP, Gouzy J, Parra G, Lardier G, Chapple C, McKernan KJ, McEwan P, Bosak S, Kellis M, Volff JN, Guigo R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechavi M, Laudet V, Schachter V, Quetier F, Saurin W, Scarpelli C, Wincker P, Lander ES, Weissenbach J, Roest Crollius H: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-57. 10.1038/nature03025.
Fukamachi S, Asakawa S, Wakamatsu Y, Shimizu N, Mitani H, Shima A: Conserved function of medaka pink-eyed dilution in melanin synthesis and its divergent transcriptional regulation in gonads among vertebrates. Genetics. 2004, 168: 1519-27. 10.1534/genetics.104.030494.
Protas ME, Hersey C, Kochanek D, Zhou Y, Wilkens H, Jeffery WR, Zon LI, Borowsky R, Tabin CJ: Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism. Nat Genet. 2006, 38: 107-11. 10.1038/ng1700.
Fukamachi S, Shimada A, Shima A: Mutations in the gene encoding B, a novel transporter protein, reduce melanin content in medaka. Nat Genet. 2001, 28: 381-5. 10.1038/ng584.
Henze M, Rempeters G, Anders F: Pteridines in the skin of xiphophorine fish (Poeciliidae). Comp Biochem Physiol B. 1977, 56: 35-46. 10.1016/0305-0491(77)90219-X.
Maier J, Witter K, Gutlich M, Ziegler I, Werner T, Ninnemann H: Homology cloning of GTP-cyclohydrolase I from various unrelated eukaryotes by reverse-transcription polymerase chain reaction using a general set of degenerate primers. Biochem Biophys Res Commun. 1995, 212: 705-11. 10.1006/bbrc.1995.2026.
Pelletier I, Bally-Cuif L, Ziegler I: Cloning and developmental expression of zebrafish GTP cyclohydrolase I. Mech Dev. 2001, 109: 99-103. 10.1016/S0925-4773(01)00516-0.
Hilton DJ, Richardson RT, Alexander WS, Viney EM, Willson TA, Sprigg NS, Starr R, Nicholson SE, Metcalf D, Nicola NA: Twenty proteins containing a C-terminal SOCS box form five structural classes. Proc Natl Acad Sci USA. 1998, 95: 114-9. 10.1073/pnas.95.1.114.
Holland PW, Garcia-Fernandez J, Williams NA, Sidow A: Gene duplications and the origins of vertebrate development. Dev. 1994, 125-33. Suppl
Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005, 3: e314-10.1371/journal.pbio.0030314.
Chavan B, Gillbro JM, Rokos H, Schallreuter KU: GTP cyclohydrolase feedback regulatory protein controls cofactor 6-tetrahydrobiopterin synthesis in the cytosol and in the nucleus of epidermal keratinocytes and melanocytes. J Invest Dermatol. 2006, 126: 2481-9. 10.1038/sj.jid.5700425.
Ben J, Lim TM, Phang VP, Chan WK: Cloning and tissue expression of 6-pyruvoyl tetrahydropterin synthase and xanthine dehydrogenase from Poecilia reticulata. Mar Biotechnol (NY). 2003, 5: 568-78. 10.1007/s10126-002-0121-y.
Schallreuter KU, Kothari S, Hasse S, Kauser S, Lindsey NJ, Gibbons NC, Hibberts N, Wood JM: In situ and in vitro evidence for DCoH/HNF-1 alpha transcription of tyrosinase in human skin melanocytes. Biochem Biophys Res Commun. 2003, 301: 610-6. 10.1016/S0006-291X(02)03076-0.
Rose RB, Pullen KE, Bayle JH, Crabtree GR, Alber T: Biochemical and structural basis for partially redundant enzymatic and transcriptional functions of DCoH and DCoH2. Biochemistry. 2004, 43: 7345-55. 10.1021/bi049620t.
Giordano E, Peluso I, Rendina R, Digilio A, Furia M: The clot gene of Drosophila melanogaster encodes a conserved member of the thioredoxin-like protein superfamily. Mol Genet Genomics. 2003, 268: 692-7.
Le Guyader S, Maier J, Jesuthasan S: Esrom, an ortholog of PAM (protein associated with c-myc), regulates pteridine synthesis in the zebrafish. Dev Biol. 2005, 277: 378-86. 10.1016/j.ydbio.2004.09.029.
Blomme T, Vandepoele K, De Bodt S, Simillion C, Maere S, Van de Peer Y: The gain and loss of genes during 600 million years of vertebrate evolution. Genome Biol. 2006, 7: R43.1-12. 10.1186/gb-2006-7-5-r43.
Steinke D, Hoegg S, Brinkmann H, Meyer A: Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates. BMC Biol. 2006, 4: 16-10.1186/1741-7007-4-16.
Postlethwait JH, Woods IG, Ngo-Hazelett P, Yan YL, Kelly PD, Chu F, Huang H, Hill-Force A, Talbot WS: Zebrafish comparative genomics and the origins of vertebrate chromosomes. Genome Res. 2000, 10: 1890-902. 10.1101/gr.164800.
Brunet FG, Crollius HR, Paris M, Aury JM, Gibert P, Jaillon O, Laudet V, Robinson-Rechavi M: Gene loss and evolutionary rates following whole-genome duplication in teleost fishes. Mol Biol Evol. 2006, 23: 1808-16. 10.1093/molbev/msl049.
He X, Zhang J: Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics. 2005, 169: 1157-64. 10.1534/genetics.104.037051.
Van de Peer Y, Taylor JS, Braasch I, Meyer A: The ghost of selection past: rates of evolution and functional divergence of anciently duplicated genes. J Mol Evol. 2001, 53: 436-46. 10.1007/s002390010233.
Steinke D, Salzburger W, Braasch I, Meyer A: Many genes in fish have species-specific asymmetric rates of molecular evolution. BMC Genomics. 2006, 7: 20-10.1186/1471-2164-7-20.
Rodriguez-Trelles F, Tarrio R, Ayala FJ: Convergent neofunctionalization by positive Darwinian selection after ancient recurrent duplications of the xanthine dehydrogenase gene. Proc Natl Acad Sci USA. 2003, 100: 13413-7. 10.1073/pnas.1835646100.
Zhang J, Rosenberg HF, Nei M: Positive Darwinian selection after gene duplication in primate ribonuclease genes. Proc Natl Acad Sci USA. 1998, 95: 3708-13. 10.1073/pnas.95.7.3708.
Zhang J, Zhang YP, Rosenberg HF: Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat Genet. 2002, 30: 411-5. 10.1038/ng852.
Taylor JS, Van de Peer Y, Meyer A: Genome duplication, divergent resolution and speciation. Trends Genet. 2001, 17: 299-301. 10.1016/S0168-9525(01)02318-6.
Scannell DR, Byrne KP, Gordon JL, Wong S, Wolfe K: Multiple rounds of speciation associated with reciprocal gene loss in polyploid yeasts. Nature. 2006, 440: 341-45. 10.1038/nature04562.
Aury JM, Jaillon O, Duret L, Noel B, Jubin C, Porcel BM, Segurens B, Daubin V, Anthouard V, Aiach N, Arnaiz O, Billaut A, Beisson J, Blanc I, Bouhouche K, Camara F, Duharcourt S, Guigo R, Gogendeau D, Katinka M, Keller AM, Kissmehl R, Klotz C, Koll F, Le Mouel A, Lepere G, Malinsky S, Nowacki M, Nowak JK, Plattner H, Poulain J, Ruiz F, Serrano V, Zagulski M, Dessen P, Betermier M, Weissenbach J, Scarpelli C, Schachter V, Sperling L, Meyer E, Cohen J, Wincker P: Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia. Nature. 2006, 444: 171-178. 10.1038/nature05230.
Page-McCaw PS, Chung SC, Muto A, Roeser T, Staub W, Finger-Baier KC, Korenbrot JI, Baier H: Retinal network adaptation to bright light requires tyrosinase. Nat Neurosci. 2004, 7: 1329-36. 10.1038/nn1344.
Iida A, Inagaki H, Suzuki M, Wakamatsu Y, Hori H, Koga A: The tyrosinase gene of the i(b) albino mutant of the medaka fish carries a transposable element insertion in the promoter region. Pigment Cell Res. 2004, 17: 158-64. 10.1046/j.1600-0749.2003.00122.x.
Loftus SK, Larson DM, Baxter LL, Antonellis A, Chen Y, Wu X, Jiang Y, Bittner M, Hammer JA, Pavan WJ: Mutation of melanosome protein RAB38 in chocolate mice. Proc Natl Acad Sci USA. 2002, 99: 4471-6. 10.1073/pnas.072087599.
Baxter LL, Pavan WJ: Pmel17 expression is Mitf-dependent and reveals cranial melanoblast migration during murine development. Gene Expr Patterns. 2003, 3: 703-7. 10.1016/j.modgep.2003.07.002.
Du J, Miller AJ, Widlund HR, Horstmann MA, Ramaswamy S, Fisher DE: MLANA/MART1 and SILV/PMEL17/GP100 are transcriptionally regulated by MITF in melanocytes and melanoma. Am J Pathol. 2003, 163: 333-43.
Donoghue PCJ, Purnell MA: Genome duplication, extinction and vertebrate evolution. Trends Ecol Evol. 2005, 20: 312-319. 10.1016/j.tree.2005.04.008.
Crow KD, Wagner GP: Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. What is the role of genome duplication in the evolution of complexity and diversity?. Mol Biol Evol. 2006, 23: 887-92. 10.1093/molbev/msj083.
Trezise AE, Collin SP: Opsins: evolution in waiting. Curr Biol. 2005, 15: R794-6. 10.1016/j.cub.2005.09.025.
Matsumoto Y, Fukamachi S, Mitani H, Kawamura S: Functional characterization of visual opsin repertoire in Medaka (Oryzias latipes). Gene. 2006, 371: 268-78. 10.1016/j.gene.2005.12.005.
Ensembl Genome Browser. [http://www.ensembl.org/]
TIGR Gene Indices. [http://www.tigr.org/tdb/tgi/]
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-80. 10.1093/nar/22.22.4673.
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-63. 10.1093/bib/5.2.150.
Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online–a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33: W557-9. 10.1093/nar/gki352.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-5. 10.1093/bioinformatics/bti263.
Swofford D: PAUP*: Phylogenetic analysis using parismony (* and other methods). 2000, Sunderland, Massachusetts: Sinauer Associates
Simple radiation hybrid mapping in zebrafish using the LN54 panel. [http://mgchd1.nichd.nih.gov:8000/zfrh/beta.cgi]
We would like to thank all members of the BioFuture Group and Agnès Dettaï, Kathrin Lampert, Christopher Untucht and Michael Fackelmann for helpful support and discussions, as well as Dirk Steinke and Simone Hoegg for sharing results prior to publication. The zebrafish radiation hybrid panel LN54 was kindly provided by Marc Ekker. Our work is funded by grants from the German Federal Ministry of Education and Research (BMBF, Biofuture program, to JNV), the German Research Society (DFG, to MS and JNV), the Association pour la Recherche sur le Cancer (ARC, to JNV), and the Centre National de la Recherche Scientifique (CNRS, to JNV).
IB participated in the design of the study, carried out bioinformatic analyses and molecular studies and wrote the manuscript. MS was involved in the design of the study and helped to draft the manuscript. JNV participated in the design of the study, helped with data analyses and to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Nucleotide accession numbers of melanin synthesis genes. GenBank accession numbers, Ensembl accession numbers or TIGR EST clusters (TC) used for phylogenetic analyses are given. EST denotes manually assembled EST clusters. Partial sequences that were not included in final phylogenetic trees are indicated by #, pseudogenes by ψ. See Table 1 for species abbreviations. (PDF 153 KB)
Additional File 2: Nucleotide accession numbers of pteridine synthesis genes. GenBank accession numbers, Ensembl accession numbers or TIGR EST clusters (TC) used for phylogenetic analyses are given. EST denotes manually assembled EST clusters. scaf: scaffold, ctg: contig of Ensembl genome assembly. Partial sequences that were not included in final phylogenetic trees are indicated by #, pseudogenes by ψ. See Table 1 for species abbreviations. (PDF 179 KB)
Additional File 3: Nucleotide phylogeny of tyrp1 genes in vertebrates. Maximum-likelihood phylogeny of tyrp1 genes based on a 1681 nucleotide alignment (GTR+I+G model; parameter values estimated from the dataset). The tree is rooted with human TYR and TYRP1 genes. Numbers at the branches denote bootstrap values (maximum likelihood/neighbor joining) above 50%. The topology of the tree is consistent with the duplication of tyrp1 during the FSGD. (PDF 144 KB)
Additional File 4: Molecular phylogeny of melanosomal transporters: Oca2, Aim1, Slc24a5. Maximum-likelihood phylogeny of (a) Oca2 (854 AA; rooted with Oca2 from sea urchin), (b) Aim1 (701 AA; rooted with human SLC45A1), and (c) Slc24a5 (651 AA; rooted with human SLC24A3/4 proteins). Numbers at the branches denote bootstrap values (maximum likelihood/neighbor joining) above 50%. No duplications were observed in teleosts. (PDF 175 KB)
Additional File 5: Molecular phylogeny of pteridine synthesis enzymes: Pts, Xod/Xdh, Tnxl5, Pam. Maximum-likelihood phylogeny of (a) Pts (158 AA; rooted with Pts from Ciona), (b) Xod/Xdh (1453 AA; rooted with human AOX1; an Aox1 sequence from guppy (Poecilia reticulata) is wrongly annotated in GenBank as Xod/Xdh), (c) Tnxl5 (130 AA; rooted with Clot from Drosophila), and (d) Pam (4932 AA; rooted with Pam from Drosophila). Numbers at the branches denote bootstrap values (maximum likelihood/neighbor joining) above 50%. No duplications were observed in teleosts. (PDF 183 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Braasch, I., Schartl, M. & Volff, JN. Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol Biol 7, 74 (2007). https://doi.org/10.1186/1471-2148-7-74
- Genome Duplication
- Fathead Minnow
- Melanin Synthesis
- Vertebrate Lineage