Evolution of pigment synthesis pathways by gene and genome duplication in fish
© Braasch et al; licensee BioMed Central Ltd. 2007
Received: 20 December 2006
Accepted: 11 May 2007
Published: 11 May 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.
Pigment synthesis genes in human and teleost fish
LG 15 (sdy)
LG 13 (i-b/-1/-4/-6)
Omy, Oni and other cichlids
Eha, Gpe, Ipu, Omy, Oni
Cau, Omy, Oni, Ppr, Ssa
Cau, Ipu, Omy, Ssa
LG 11 (fdv)
Hsp, Omy, Ppr
Abu, Omy, Ssa
LG 4 (i-3)
LG 12 (b)
Omy, several cichlids
LG 18 (gol)
Hhi, Omy, Ppr
Fhe, Omy, Ssa
Omy, Pfl, Ppr, Pre, Ssa
Ipu, Omy, Pfl, Pol, Ppr, Ssa
Fhe, Sch, Ssa
Cca, Ipu, Omy, Ppr, Ssa
scaf 53, scaf 178ψ
Fhe, Hhi, Ipu, Omy, Pfl
Man, Omy, Ppr, Ssa
Cca, Omy, Ppr, Psa, Ssa, Hsp
Omy, Ppr, Ssa
LG 9 (esr)
Abu, Hhi, Omy, Ppr
Gene duplications in the melanin synthesis pathway
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 184.108.40.206) promotes the first two rate-limiting steps of melanin synthesis from tyrosine to DOPA and DOPAquinone as well as two later steps. Dct (EC 220.127.116.11) 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].
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.
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
GTP cyclohydrolase I and its feedback regulatory protein
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 18.104.22.168) and the sepiapterin reductase (Spr; EC 22.214.171.124) (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].
Enzymes of the H4biopterin regeneration pathway
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).
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).
- Fujii R: Coloration and chromatophores. The Physiology of Fishes. Edited by: Evans DH. 1993, Boca Raton, Florida: CRC Press, 535-562.Google Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Froschauer A, Braasch I, Volff JN: Fish genomes, comparative genomics and vertebrate evolution. Current Genomics. 2006, 1: 43-57. 10.2174/138920206776389766.View ArticleGoogle Scholar
- Volff JN: Genome evolution and biodiversity in teleost fish. Heredity. 2005, 94: 280-94. 10.1038/sj.hdy.6800635.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Hearing VJ, Tsukamoto K: Enzymatic control of pigmentation in mammals. FASEB J. 1991, 5: 2902-9.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Postlethwait J: The zebrafish genome: a review and msx gene case study. Genome Dyn. 2006, 2: 183-97.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ohno S: Evolution by gene duplication. 1970, New York: Springer-VerlagView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Holland PW, Garcia-Fernandez J, Williams NA, Sidow A: Gene duplications and the origins of vertebrate development. Dev. 1994, 125-33. SupplGoogle Scholar
- Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005, 3: e314-10.1371/journal.pbio.0030314.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Donoghue PCJ, Purnell MA: Genome duplication, extinction and vertebrate evolution. Trends Ecol Evol. 2005, 20: 312-319. 10.1016/j.tree.2005.04.008.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Trezise AE, Collin SP: Opsins: evolution in waiting. Curr Biol. 2005, 15: R794-6. 10.1016/j.cub.2005.09.025.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-5. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
- Swofford D: PAUP*: Phylogenetic analysis using parismony (* and other methods). 2000, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Simple radiation hybrid mapping in zebrafish using the LN54 panel. [http://mgchd1.nichd.nih.gov:8000/zfrh/beta.cgi]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.