Comparative Analysis of <Emphasis Type="Italic">Period</Emphasis> Genes in Teleost Fish Genomes

Period (Per) is a canonical circadian clock gene. The fruit fly, an invertebrate, has one per gene, while the human, a tetrapod vertebrate, has three Per genes. Per1, Per2, and Per3 of the tetrapods were generated from two rounds of ancient genome duplications from the ancestral chordate Per gene. Searching for five teleost fish genomes in a combination of phylogenetic, splicing site, and syntenic analyses revealed that zebrafish have two per1 genes, per1a and per1b, one per2, and one per3; medaka, fugu, and tetraodon each have two per2 genes, per2a and per2b, one per1, and one per3; sticklebacks also have per2a, per2b, and one per1 but lack per3; and per1a/per1b in zebrafish and per2a/per2b in madaka, fugu, tetraodon, and stickleback are ancient duplicates. While the dN/dS ratios of the five fish per duplicates are all <1, suggesting that they likely have been subject to purifying selection, the Tajima relative rate test showed that zebrafish per1a/per1b and fugu and medaka per2a/per2b have asymmetric evolutionary rates, implicating that one of these duplicates might have been under positive selection or relaxed functional constraint. Further, in situ hybridization showed that zebrafish per1a and per1b clearly have distinct patterns of temporal and spatial expression. These results support the notion that extra copies of teleost per genes were generated from the fish-specific genome duplication, and divergent resolution after the duplication resulted in retention of different per duplicates in different fish, most of which have diverged significantly.     

Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative genomics of duplicate γ-glutamyl transferase genes in teleosts: medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), green spotted pufferfish (Tetraodon nigroviridis), fugu (Takifugu rubripes), and zebrafish (Danio rerio)

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we examine multiple paralogous genes of γ-glutamyl transferase (GGT) in several distantly related teleost species including medaka, stickleback, green spotted pufferfish, fugu, and zebrafish. Through mining genome databases, we have identified multiple GGT orthologs. Duplicate (paralogous) GGT sequences for GGT1 (GGT1 a and b), GGTL1 (GGTL1 a and b), and GGTL3 (GGTL3 a and b) were identified for each species. Phylogenetic analysis suggests that GGTs are ancient proteins conserved across most metazoan phyla and those paralogous GGTs in teleosts likely arose from the serial 3R genome duplication events. A third GGTL1 gene (GGTL1c) was found in green spotted pufferfish; however, this gene is not present in medaka, stickleback, or fugu. Similarly, one or both paralogs of GGTL3 appear to have been lost in green spotted pufferfish, fugu, and zebrafish. Syntenic relationships were highly maintained between duplicated teleost chromosomes, among teleosts and across ray-finned (Actinopterygii) and lobe-finned (Sarcopterygii) species. To assess subfunction partitioning, six medaka GGT genes were cloned and assessed for developmental and tissue-specific expression. On the basis of these data, we propose a modification of the "duplication-degeneration-complementation" model of subfunction partitioning where quantitative differences rather than absolute differences in gene expression are observed between gene paralogs. Our results demonstrate that multiple GGT genes have been retained within teleost genomes. Questions remain, however, regarding the functional roles of multiple GGTs in these species.     

Comparative genomics provides evidence for an ancient genome duplication event in fish

There are approximately 25 000 species in the division Teleostei and most are believed to have arisen during a relatively short period of time ca. 200 Myr ago. The discovery of 'extra' Hox gene clusters in zebrafish (Danio rerio), medaka (Oryzias latipes), and pufferfish (Fugu rubripes), has led to the hypothesis that genome duplication provided the genetic raw material necessary for the teleost radiation. We identified 27 groups of orthologous genes which included one gene from man, mouse and chicken, one or two genes from tetraploid Xenopus and two genes from zebrafish. A genome duplication in the ancestor of teleost fishes is the most parsimonious explanation for the observations that for 15 of these genes, the two zebrafish orthologues are sister sequences in phylogenies that otherwise match the expected organismal tree, the zebrafish gene pairs appear to have been formed at approximately the same time, and are unlinked. Phylogenies of nine genes differ a little from the tree predicted by the fish-specific genome duplication hypothesis: one tree shows a sister sequence relationship for the zebrafish genes but differs slightly from the expected organismal tree and in eight trees, one zebrafish gene is the sister sequence to a clade which includes the second zebrafish gene and orthologues from Xenopus, chicken, mouse and man. For these nine gene trees, deviations from the predictions of the fish-specific genome duplication hypothesis are poorly supported. The two zebrafish orthologues for each of the three remaining genes are tightly linked and are, therefore, unlikely to have been formed during a genome duplication event. We estimated that the unlinked duplicated zebrafish genes are between 300 and 450 Myr. Thus, genome duplication could have provided the genetic raw material for teleost radiation. Alternatively, the loss of different duplicates in different populations (i.e. 'divergent resolution') may have promoted speciation in ancient teleost populations.     

Consequences of hoxb1 duplication in teleost fish

Vertebrate evolution is characterized by gene and genome duplication events. There is strong evidence that a whole-genome duplication occurred in the lineage leading to the teleost fishes. We have focused on the teleost hoxb1 duplicate genes as a paradigm to investigate the consequences of gene duplication. Previous analysis of the duplicated zebrafish hoxb1 genes suggested they have subfunctionalized. The combined expression pattern of the two zebrafish hoxb1 genes recapitulates the expression pattern of the single Hoxb1 gene of tetrapods, possibly due to degenerative changes in complementary cis-regulatory elements of the duplicates. Here we have tested the hypothesis that all teleost duplicates had a similar fate post duplication, by examining hoxb1 genes in medaka and striped bass. Consistent with this theory, we found that the ancestral Hoxb1 expression pattern is subdivided between duplicate genes in a largely similar fashion in zebrafish, medaka, and striped bass. Further, our analysis of hoxb1 genes reveals that sequence changes in cis-regulatory regions may underlie subfunctionalization in all teleosts, although the specific changes vary between species. It was previously shown that zebrafish hoxb1 duplicates have also evolved different functional capacities. We used misexpression to compare the functions of hoxb1 duplicates from zebrafish, medaka and striped bass. Unexpectedly, we found that some biochemical properties, which were paralog specific in zebrafish, are conserved in both duplicates of other species. This work suggests that the fate of duplicate genes varies across the teleost group.     

Phylogenetic Timing of the Fish-Specific Genome Duplication Correlates with the Diversification of Teleost Fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages (Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes (fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.[Reviewing Editor: Martin Kreitman]     

Duplication and divergence of fgf8 functions in teleost development and evolution

Fibroblast growth factors play critical roles in many aspects of embryo patterning that are conserved across broad phylogenetic distances. To help understand the evolution of fibroblast growth factor functions, we identified members of the Fgf8/17/18-subfamily in the three-spine stickleback Gasterosteus aculeatus, and investigated their evolutionary relationships and expression patterns. We found that fgf17b is the ortholog of tetrapod Fgf17, whereas the teleost genes called fgf8 and fgf17a are duplicates of the tetrapod gene Fgf8, and thus should be called fgf8a and fgf8b. Phylogenetic analysis supports the view that the Fgf8/17/18-subfamily expanded during the ray-fin fish genome duplication. In situ hybridization experiments showed that stickleback fgf8 duplicates exhibited common and unique expression patterns, indicating that tissue specialization followed the gene duplication event. Moreover, direct comparison of stickleback and zebrafish embryonic expression patterns of fgf8 co-orthologs suggested lineage-specific independent subfunction partitioning and the acquisition or the loss of ortholog functions. In tetrapods, Fgf8 plays an important role in the apical ectodermal ridge of the developing pectoral appendage. Surprisingly, differences in the expression of fgf8a in the apical ectodermal ridge of the pectoral fin bud in zebrafish and stickleback, coupled with the role of fgf16 and fgf24 in teleost pectoral appendage show that different Fgf genes may play similar roles in limb development in various vertebrates.     

The "fish-specific" Hox cluster duplication is coincident with the origin of teleosts

The Hox gene complement of zebrafish, medaka, and fugu differs from that of other gnathostome vertebrates. These fishes have seven to eight Hox clusters compared to the four Hox clusters described in sarcopterygians and shark. The clusters in different teleost lineages are orthologous, implying that a "fish-specific" Hox cluster duplication has occurred in the stem lineage leading to the most recent common ancestor of zebrafish and fugu. The timing of this event, however, is unknown. To address this question, we sequenced four Hox genes from taxa representing basal actinopterygian and teleost lineages and compared them to known sequences from shark, coelacanth, zebrafish, and other teleosts. The resulting gene genealogies suggest that the fish-specific Hox cluster duplication occurred coincident with the origin of crown group teleosts. In addition, we obtained evidence for an independent Hox cluster duplication in the sturgeon lineage (Acipenseriformes). Finally, results from HoxA11 suggest that duplicated Hox genes have experienced diversifying selection immediately after the duplication event. Taken together, these results support the notion that the duplicated Hox genes of teleosts were causally relevant to adaptive evolution during the initial teleost radiation.     

Diverse forms of guanylyl cyclases in medaka fish -- their genomic structure and phylogenetic relationships to those in vertebrates and invertebrates

Fish species such as medaka fish, fugu, and zebrafish contain more guanylyl cyclases (GCs) than do mammals. These GCs can be divided into two types: soluble GCs and membrane GCs. The latter are further divided into four subfamilies: (i) natriuretic peptide receptors, (ii) STa/guanylin receptors, (iii) sensory-organ-specific membrane GCs, and (iv) orphan receptors. Phylogenetic analyses of medaka fish GCs, along with those of fugu and zebrafish, suggest that medaka fish is a much closer relative to fugu than to zebrafish. Analyses of nucleotide data available on a web site (http://www.ncbi. nlm.nih.gov/) of GCs from a range of organisms from bacteria to vertebrates suggest that gene duplication, and possibly chromosomal duplication, play important roles in the divergence of GCs. In particular, the membrane GC genes were generated by chromosomal duplication before the divergence of tetrapods and teleosts.     

Evolutionary Patterns and Selective Pressures of Odorant/Pheromone Receptor Gene Families in Teleost Fishes

Background

Teleost fishes do not have a vomeronasal organ (VNO), and their vomeronasal receptors (V1Rs, V2Rs) are expressed in the main olfactory epithelium (MOE), as are odorant receptors (ORs) and trace amine-associated receptors (TAARs). In this study, to obtain insights into the functional distinction among the four chemosensory receptor families in teleost fishes, their evolutionary patterns were examined in zebrafish, medaka, stickleback, fugu, and spotted green pufferfish.

Methodology/Principal Findings

Phylogenetic analysis revealed that many lineage-specific gene gains and losses occurred in the teleost fish TAARs, whereas only a few gene gains and losses have taken place in the teleost fish vomeronasal receptors. In addition, synonymous and nonsynonymous nucleotide substitution rate ratios (K_A/K_S) in TAARs tended to be higher than those in ORs and V2Rs.

Conclusions/Significance

Frequent gene gains/losses and high K_A/K_S in teleost TAARs suggest that receptors in this family are used for detecting some species-specific chemicals such as pheromones. Conversely, conserved repertoires of V1R and V2R families in teleost fishes may imply that receptors in these families perceive common odorants for teleosts, such as amino acids. Teleost ORs showed intermediate evolutionary pattern between TAARs and vomeronasal receptors. Many teleost ORs seem to be used for common odorants, but some ORs may have evolved to recognize lineage-specific odors.     

Are all fishes ancient polyploids?

Euteleost fishes seem to have more copies of many genes than their tetrapod relatives. Three different mechanisms could explain the origin of these 'extra' fish genes. The duplicates may have been produced during a fish-specific genome duplication event. A second explanation is an increased rate of independent gene duplications in fish. A third possibility is that after gene or genome duplication events in the common ancestor of fish and tetrapods, the latter lost more genes. These three hypotheses have been tested by phylogenetic tree reconstruction. Phylogenetic analyses of sequences from human, mouse, chicken, frog (Xenopus laevis), zebrafish (Danio rerio) and pufferfish (Takifugu rubripes) suggest that ray-finned fishes are likely to have undergone a whole genome duplication event between 200 and 450 million years ago. We also comment here on the evolutionary consequences of this ancient genome duplication.     

Identification and analysis of additional copies of the platelet-derived growth factor receptor and colony stimulating factor 1 receptor genes in fugu

The receptors for the platelet-derived growth factor (PDGFRalpha and PDGFRbeta) belong to a subfamily of protein tyrosine kinase receptors that also includes kit and the colony stimulating factor-1 receptor (CSF1R). In mammals, the genes encoding PDGFRalpha and PDGFRbeta are tandemly linked to the kit and CSF1R genes, respectively. Based on the structural similarity and genomic organization of these four genes, it has been suggested that they arose from an ancestral protein tyrosine kinase receptor gene by two rounds of duplication. We have previously cloned the PDGFRbeta and CSF1R genes from the pufferfish, Fugu rubripes, and shown that they are tandemly linked like the mammalian genes [Genome Res. 6 (1996) 1185]. We have now cloned two additional members of this gene family, fPDGFRbeta2 and fCSF1R2 from the fugu and shown that these two genes are also tandemly linked. This indicates that the PDGFRbeta-CSF1R locus has been duplicated in the lineage leading to fugu. The fugu fPDGFRbeta2 and fCSF1R2 genes contain three and one extra introns, respectively, compared with other members of this family. Polymerase chain reaction cloning of a conserved region of PDGFRbeta gene from other ray-finned fishes identified two copies in the zebrafish (order Cypriniformes) and sunfish (order Tetraodontiformes). These results are discussed in the context of the proposed teleost lineage-specific whole genome duplication hypothesis.     

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.     

Duplication of phospholipase C-delta gene family in fish genomes

Fishes possess more genes than other vertebrates, possibly because of a genome duplication event during the evolution of the teleost (ray-finned) fish lineage. To further explore this idea, we cloned five genes encoding phosphoinositide-specific phospholipase C-delta (PLC-delta), designated respectively PoPLC-deltas, from olive flounder (Paralichthys olivaceus), and we performed phylogenetic analysis and sequence comparison to compare our putative gene products (PoPLC-deltas) with the sequences of known human PLC isoforms. The deduced amino acid sequences shared high sequence identity with human PLC-delta1, -delta3, and -delta4 isozymes and exhibited similar primary structures. In phylogenetic analysis of PoPLC-deltas with PLC-deltas of five teleost fishes (zebrafish, stickleback, medaka, Tetraodon, and Takifugu), three tetrapods (human, chicken, and frog), and two tunicates (sea squirt and pacific sea squirt), whose putative sequences of PLC-delta are available in Ensembl genome browser, the result also indicated that the two paralogous genes corresponding to each PLC-delta isoform originated from fish-specific genome duplication prior to the divergence of teleost fish. Our analyses suggest that an ancestral PLC-delta gene underwent three rounds of genome duplication during the evolution of vertebrates, leading to the six genes of three PLC-delta isoforms in teleost fish.     

Phylogenetic and evolutionary relationships and developmental expression patterns of the zebrafish <Emphasis Type="Italic">twist</Emphasis> gene family

Four members of the twist gene family (twist1a, 1b, 2, and 3) are found in the zebrafish, and they are thought to have arisen through three rounds of gene duplication, two of which occurred prior to the tetrapod-fish split. Phylogenetic analysis groups most of the vertebrate Twist1 peptides into clade I, except for the Twist1b proteins of the acanthopterygian fish (medaka, pufferfish, stickleback), which clustered within clade III. Paralogies and orthologies among the zebrafish, medaka, and human twist genes were determined using comparative synteny analysis of the chromosomal regions flanking these genes. Comparative nucleotide substitution analyses also revealed a faster rate of nucleotide mutation/substitution in the acanthopterygian twist1b compared to the zebrafish twist1b, thus accounting for their anomalous phylogenetic clustering. We also observed minimal expression overlap among the four twist genes, suggesting that despite their significant peptide similarity, their regulatory controls have diverged considerably, with minimal functional redundancy between them. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative analysis of her genes during fish somitogenesis suggests a mouse/chick-like mode of oscillation in medaka

Somitogenesis is the key developmental step, which divides the vertebrate body axis into segmentally repeated structures. It requires an intricate process of pre-patterning, which is driven by an oscillator mechanism consisting of the Delta–Notch pathway and various hairy- and Enhancer of split-related (her) genes. The subset of her genes, which are necessary to set up the segmentation clock, reveal a complex scenario of interactions. To understand which her genes are essential core players in this process, we compared the expression patterns of somitogenesis-relevant her genes in zebrafish and medaka (Oryzias latipes). Most of the respective medaka genes (Ol-her) are duplicated like what has been shown for zebrafish (Dr-her) and pufferfish genes (Fr-her). However, zebrafish genes show some additional copies and significant differences in expression patterns. For the paralogues Dr-her1 and Dr-her11, only one copy exists in the medaka (Ol-her1/11), which combines the expression patterns found for both zebrafish genes. In contrast to Dr-her5, the medaka orthologue appears to play a role in somitogenesis because it is expressed in the presomitic mesoderm (PSM). PSM expression also suggests a role for both Ol-her13 genes, homologues of mouse Hes6 (mHes6), in this process, which would be consistent with a conserved mHes6 homologue gear in the segmentation clock exclusively in lower vertebrates. Members of the mHes5 homologue group seem to be involved in somite formation in all vertebrates (e.g. Dr- and Ol-her12), although different paralogues are additionally recruited in zebrafish (e.g. Dr-her15) and medaka (e.g. Ol-her4). We found that the linkage between duplicates is strongly conserved between pufferfish and medaka and less well conserved in zebrafish. Nevertheless, linkage and orientation of several her duplicates are identical in all three species. Therefore, small-scale duplications must have happened before whole genome duplication occurred in a fish ancestor. Expression of multiple stripes in the intermediate PSM, characteristic for the zebrafish orthologues, is absent in all somitogenesis-related her genes of the medaka. In fact, the expression mode of Ol-her1/11 and Ol-her5 indicates dynamism similar to the hairy clock genes in chicken and mouse. This suggests that Danio rerio shows a rather derived clock mode when compared to other fish species and amniotes or that, alternatively, the clock mode evolved independently in zebrafish, medaka and mouse or chicken.An erratum to this article can be found at     

Evolution of trace amine associated receptor (TAAR) gene family in vertebrates: lineage-specific expansions and degradations of a second class of vertebrate chemosensory receptors expressed in the olfactory epithelium

The trace amine-associated receptors (TAARs) form a specific family of G protein-coupled receptors in vertebrates. TAARs were initially considered neurotransmitter receptors, but recent study showed that mouse TAARs function as chemosensory receptors in the olfactory epithelium. To clarify the evolutionary dynamics of the TAAR gene family in vertebrates, near-complete repertoires of TAAR genes and pseudogenes were identified from the genomic assemblies of 4 teleost fishes (zebrafish, fugu, stickleback, and medaka), western clawed frogs, chickens, 3 mammals (humans, mice, and opossum), and sea lampreys. Database searches revealed that fishes had many putatively functional TAAR genes (13-109 genes), whereas relatively small numbers of TAAR genes (3-22 genes) were identified in tetrapods. Phylogenetic analysis of these genes indicated that the TAAR gene family was subdivided into 5 subfamilies that diverged before the divergence of ray-finned fishes and tetrapods. In tetrapods, virtually all TAAR genes were located in 1 specific region of their genomes as a gene cluster; however, in fishes, TAAR genes were scattered throughout more than 2 genomic locations. This possibly reflects a whole-genome duplication that occurred in the common ancestor of ray-finned fishes. Expression analysis of zebrafish and stickleback TAAR genes revealed that many TAARs in these fishes were expressed in the olfactory organ, suggesting the relatively high importance of TAARs as chemosensory receptors in fishes. A possible evolutionary history of the vertebrate TAAR gene family was inferred from the phylogenetic and comparative genomic analyses.