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]     

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.     

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.     

The fate of the duplicated androgen receptor in fishes: a late neofunctionalization event?

Background  

Based on the observation of an increased number of paralogous genes in teleost fishes compared with other vertebrates and on the conserved synteny between duplicated copies, it has been shown that a whole genome duplication (WGD) occurred during the evolution of Actinopterygian fish. Comparative phylogenetic dating of this duplication event suggests that it occurred early on, specifically in teleosts. It has been proposed that this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish, notably by allowing the sub- or neo-functionalization of many duplicated genes.     

Gene Duplications and Evolution of Vertebrate Voltage-Gated Sodium Channels

Voltage-gated sodium channels underlie action potential generation in excitable tissue. To establish the evolutionary mechanisms that shaped the vertebrate sodium channel α-subunit (SCNA) gene family and their encoded Na_v1 proteins, we identified all SCNA genes in several teleost species. Molecular cloning revealed that teleosts have eight SCNA genes, compared to ten in another vertebrate lineage, mammals. Prior phylogenetic analyses have indicated that the genomes of both teleosts and tetrapods contain four monophyletic groups of SCNA genes, and that tandem duplications expanded the number of genes in two of the four mammalian groups. However, the number of genes in each group varies between teleosts and tetrapods, suggesting different evolutionary histories in the two vertebrate lineages. Our findings from phylogenetic analysis and chromosomal mapping of Danio rerio genes indicate that tandem duplications are an unlikely mechanism for generation of the extant teleost SCNA genes. Instead, analyses of other closely mapped genes in D. rerio as well as of SCNA genes from several teleost species all support the hypothesis that a whole-genome duplication was involved in expansion of the SCNA gene family in teleosts. Interestingly, despite their different evolutionary histories, mRNA analyses demonstrated a conservation of expression patterns for SCNA orthologues in teleosts and tetrapods, suggesting functional conservation. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Axel Meyer]     

Teleost fish with specific genome duplication as unique models of vertebrate evolution

Whole-genome duplication (WGD) is believed to be one of the major evolutionary events that shaped the genome organization of vertebrates. Here, we review recent research on vertebrate genome evolution, specifically on WGD and its consequences for gene and genome evolution in teleost fishes. Recent genome analyses confirmed that all vertebrates experienced two rounds of WGD early in their evolution, and that teleosts experienced a subsequent additional third-round (3R)-WGD. The 3R-WGD was estimated to have occurred 320–400 million years ago in a teleost ancestor, but after its divergence from a common ancestor with living non-teleost actinopterygians (Bichir, Sturgeon, Bowfin, and Gar) based on the analyses of teleost-specific duplicate genes. This 3R-WGD was confirmed by synteny analysis and ancestral karyotype inference using the genome sequences of Tetraodon and medaka. Most of the tetrapods, on the other hand, have not experienced an additional WGD; however, they have experienced repeated chromosomal rearrangements throughout the whole genome. Therefore, different types of chromosomal events have characterized the genomes of teleosts and tetrapods, respectively. The 3R-WGD is useful to investigate the consequences of WGD because it is an evolutionarily recent WGD and thus teleost genomes retain many more WGD-derived duplicates and “traces” of their evolution. In addition, the remarkable morphological, physiological, and ecological diversity of teleosts may facilitate understanding of macrophenotypic evolution on the basis of genetic/genomic information. We highlight the teleosts with 3R-WGD as unique models for future studies on ecology and evolution taking advantage of emerging genomics technologies and systems biology environments.     

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.     

The evolution of teleost pigmentation and the fish‐specific genome duplication

Teleost fishes have evolved a unique complexity and diversity of pigmentation and colour patterning that is unmatched among vertebrates. Teleost colouration is mediated by five different major types of neural‐crest derived pigment cells, while tetrapods have a smaller repertoire of such chromatophores. The genetic basis of teleost colouration has been mainly uncovered by the cloning of pigmentation genes in mutants of zebrafish Danio rerio and medaka Oryzias latipes. Many of these teleost pigmentation genes were already known as key players in mammalian pigmentation, suggesting partial conservation of the corresponding developmental programme among vertebrates. Strikingly, teleost fishes have additional copies of many pigmentation genes compared with tetrapods, mainly as a result of a whole‐genome duplication that occurred 320–350 million years ago at the base of the teleost lineage, the so‐called fish‐specific genome duplication. Furthermore, teleosts have retained several duplicated pigmentation genes from earlier rounds of genome duplication in the vertebrate lineage, which were lost in other vertebrate groups. It was hypothesized that divergent evolution of such duplicated genes may have played an important role in pigmentation diversity and complexity in teleost fishes, which therefore not only provide important insights into the evolution of the vertebrate pigmentary system but also allow us to study the significance of genome duplications for vertebrate biodiversity.     

Evolution of follistatin in teleosts revealed through phylogenetic, genomic and expression analyses

Follistatin (Fst) inhibits transforming growth factor-β (TGF-B) proteins and is a known regulator of amniote myogenesis. Here, we used phylogenetic, genomic and experimental approaches to study its evolution in teleosts. Phylogenetic analyses suggested that one fst gene (fst1) is common to euteleosts, but a second gene (fst2) is conserved specifically within the Ostariophysi. Zebrafish fst1/2 respectively appear on chromosomes 5 and 10 in two genomic regions, each with conserved synteny to a single region in tetrapods. Interestingly, other teleosts have two corresponding chromosomal regions with a similar repertoire of paralogues. Phylogenetic reconstruction clustered these gene duplicates into two sister clades branching from tetrapod sequences. We suggest that an ancestral fst-containing chromosome was duplicated during the teleost whole genome duplication, but that fst2 was lost in lineages external to the Ostariophysi. We show that Fst1 of teleosts/mammals has evolved under strong purifying selection, but the N-terminal of Fst2 may have evolved under positive selection. Furthermore, the tissue-specific expression of zebrafish fst2 was restricted to fewer tissues compared to its paralogue and the single fst1 orthologue of Atlantic salmon (Salmo salar). Zebrafish fst1/2 may have subfunctionalized relative to non-duplicated vertebrate lineages, as several regions in the fst promoter of tetrapods were conserved with one paralogue, but not both. Finally, we examined the embryonic expression of fst1 in a teleost outside the Ostariophysi (Atlantic salmon). During segmentation, fst1 was expressed in the anterior somite compartment but was excluded from muscle progenitors that strongly expressed myogenic regulatory factors (MRFs). Later, fst1 was expressed in myogenic progenitors of the pectoral fin buds and also within the pax7 ⁺ cell layer external to the myotome. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The duplication of the <Emphasis Type="Italic">Hox</Emphasis> gene clusters in teleost fishes

Higher teleost fishes, including zebrafish and fugu, have duplicated their Hox genes relative to the gene inventory of other gnathostome lineages. The most widely accepted theory contends that the duplicate Hox clusters orginated synchronously during a single genome duplication event in the early history of ray-finned fishes. In this contribution we collect and re-evaluate all publicly available sequence information. In particular, we show that the short Hox gene fragments from published PCR surveys of the killifish Fundulus heteroclitus, the medaka Oryzias latipes and the goldfish Carassius auratus can be used to determine with little ambiguity not only their paralog group but also their membership in a particular cluster. Together with a survey of the genomic sequence data from the pufferfish Tetraodon nigroviridis we show that at least percomorpha, and possibly all eutelosts, share a system of 7 or 8 orthologous Hox gene clusters. There is little doubt about the orthology of the two teleost duplicates of the HoxA and HoxB clusters. A careful analysis of both the coding sequence of Hox genes and of conserved non-coding sequences provides additional support for the “duplication early” hypothesis that the Hox clusters in teleosts are derived from eight ancestral clusters by means of subsequent gene loss; the data remain ambiguous, however, in particular for the HoxC clusters. Assuming the “duplication early” hypothesis we use the new evidence on the Hox gene complements to determine the phylogenetic positions of gene-loss events in the wake of the cluster duplication. Surprisingly, we find that the resolution of redundancy seems to be a slow process that is still ongoing. A few suggestions on which additional sequence data would be most informative for resolving the history of the teleostean Hox genes are discussed. Supplemental material is available at http://www.bioinf.uni-leipzig.de/Publications/SUPPLEMENTS/04-006/.     

Loss of the Birt–Hogg–Dubé gene product folliculin induces longevity in a hypoxia‐inducible factor–dependent manner

Signaling through the hypoxia‐inducible factor hif‐1 controls longevity, metabolism, and stress resistance in Caenorhabditis elegans. Hypoxia‐inducible factor (HIF) protein levels are regulated through an evolutionarily conserved ubiquitin ligase complex. Mutations in the VHL gene, encoding a core component of this complex, cause a multitumor syndrome and renal cell carcinoma in humans. In the nematode, deficiency in vhl‐1 promotes longevity mediated through HIF‐1 stabilization. However, this longevity assurance pathway is not yet understood. Here, we identify folliculin (FLCN) as a novel interactor of the hif‐1/vhl‐1 longevity pathway. FLCN mutations cause Birt–Hogg–Dubé syndrome in humans, another tumor syndrome with renal tumorigenesis reminiscent of the VHL disease. Loss of the C. elegans ortholog of FLCN F22D3.2 significantly increased lifespan and enhanced stress resistance in a hif‐1‐dependent manner. F22D3.2, vhl‐1, and hif‐1 control longevity by a mechanism distinct from insulin‐like signaling. Daf‐16 deficiency did not abrogate the increase in lifespan mediated by flcn‐1. These findings define FLCN as a player in HIF‐dependent longevity signaling and connect organismal aging, stress resistance, and regulation of longevity with the formation of renal cell carcinoma.     

Evolutionary divergence of the APETALA1 and CAULIFLOWER proteins

Abstract APETALA1 (AP1) and CAULIFLOWER (CAL) are a pair of paralogous genes that were generated through the pre‐Brassicaceae whole‐genome duplication event. AP1 and CAL have both partially redundant and unique functions. Previous studies have shown that the K and C regions of their proteins are essential for the functional divergence. However, which differences in these regions are the major contributors and how the differences were accumulated remain unknown. In the present study, we compared the sequences of the two proteins and identified five gaps and 55 amino acid replacements between them. Investigation of genomic sequences further indicated that the differences in the proteins were caused by non‐synonymous substitutions and changes in exon–intron structures. Reconstruction of three‐dimensional structures revealed that the sequence divergence of AP1 and CAL has resulted in differences between the two in terms of the number, length, position and orientation of α‐helices, especially in the K and C regions. Comparisons of sequences and three‐dimensional structures of ancestral proteins with AP1 and CAL suggest that the ancestral AP1 protein experienced fewer changes, whereas the ancestral CAL protein accumulated more changes shortly after gene duplication, relative to their common ancestor. Thereafter, AP1‐like proteins experienced few mutations, whereas CAL‐like proteins were not conserved until the diversification of the Brassicaceae lineage I. This indicates that AP1‐ and CAL‐like proteins evolved asymmetrically after gene duplication. These findings provide new insights into the functional divergence of AP1 and CAL genes.     

Phosphorylation network rewiring by gene duplication

Elucidating how complex regulatory networks have assembled during evolution requires a detailed understanding of the evolutionary dynamics that follow gene duplication events, including changes in post‐translational modifications. We compared the phosphorylation profiles of paralogous proteins in the budding yeast Saccharomyces cerevisiae to that of a species that diverged from the budding yeast before the duplication of those genes. We found that 100 million years of post‐duplication divergence are sufficient for the majority of phosphorylation sites to be lost or gained in one paralog or the other, with a strong bias toward losses. However, some losses may be partly compensated for by the evolution of other phosphosites, as paralogous proteins tend to preserve similar numbers of phosphosites over time. We also found that up to 50% of kinase–substrate relationships may have been rewired during this period. Our results suggest that after gene duplication, proteins tend to subfunctionalize at the level of post‐translational regulation and that even when phosphosites are preserved, there is a turnover of the kinases that phosphorylate them.     

Dietary variations within a family of ambush predators (Platycephalidae) occupying different habitats and environments in the same geographical region

This study has determined the extents and basis for variations in the composition of the prey ingested by the abundant species of a family highly adapted for ambush predation, i.e. Platycephalidae, in a region (south‐western Australia) where that family is found in different habitats and environments. Dietary data were thus collected for Leviprora inops and Platycephalus laevigatus from seagrass in marine embayments and for Platycephalus westraliae from over sand in an estuary. These were then collated with those recorded previously for Platycephalus speculator from over sand and in seagrass in an estuary and for Platycephalus longispinis from over sand in coastal marine waters. While crustaceans and teleosts together dominated the diet of all five species, their percentage volumetric dietary contributions varied greatly, with those of crustaceans ranging from 7% for L. inops to 65% for P. speculator and those of teleosts ranging from 29% for P. longispinis to 91% for L. inops. For analyses, the data were separated into two sets. The first comprised the 17 dietary categories of invertebrates and all identified and unidentified teleosts collectively, while the second consisted of the 23 identified teleost families, both of which were subjected to permutational analysis of variance (PERMANOVA), analysis of similarities (ANOSIM) and a new (two‐way) version of the RELATE procedure. The diets of three species changed seasonally, when using invertebrate dietary categories and teleosts collectively, but with only one species, when employing identified teleost families, probably reflecting a greater tendency for invertebrate than teleost prey abundance to change during the year. On the basis of dietary data for invertebrate taxa + teleosts collectively, the diets of three of the five species changed serially with body size, with a fourth species feeding, throughout life, predominantly on the carid Palaemonetes australis. Based on identified teleost families, the diets of the three species that fed predominantly on teleosts underwent serial size‐related changes. Although L. inops and the co‐occurring P. laevigatus both consume large volumes of teleosts, the former ingests larger, less demersal and more mobile prey, e.g. the labrids Haletta semifasciata and Neoodax balteatus, than the latter, e.g. the scorpaenid Gymnapistes marmoratus, reflecting the possession by L. inops of a far longer head and larger buccal cavity. Circumstantial evidence suggests that the large differences in the volumes of crustaceans and teleosts consumed by each platycephalid species are related to differences in the relative availability of these prey in the different habitats or environments of each species.     

Investigating ancient duplication events in the Arabidopsis genome

The complete genomic analysis of Arabidopsis thaliana has shown that a major fraction of the genome consists of paralogous genes that probably originated through one or more ancient large-scale gene or genome duplication events. However, the number and timing of these duplications still remains unclear, and several different hypotheses have been put forward recently. Here, we reanalyzed duplicated blocks found in the Arabidopsis genome described previously and determined their date of divergence based on silent substitution estimations between the paralogous genes and, where possible, by phylogenetic reconstruction. We show that methods based on averaging protein distances of heterogeneous classes of duplicated genes lead to unreliable conclusions and that a large fraction of blocks duplicated much more recently than assumed previously. We found clear evidence for one large-scale gene or even complete genome duplication event somewhere between 70 to 90 million years ago. Traces pointing to a much older (probably more than 200 million years) large-scale gene duplication event could be detected. However, for now it is impossible to conclude whether these old duplicates are the result of one or more large-scale gene duplication events. abbreviations dA, fraction of amino acid substitutions; Kn, number of nonsynonymous substitutions per nonsynonymous site; Ks, number of synonymous substitutions per synonymous site; MYA, million years ago     

Comparative analysis of teleost fish genomes reveals preservation of different ancient clock duplicates in different fishes

Clock (Circadian locomotor output cycle kaput) was the first vertebrate circadian clock gene identified in a mouse forward genetics mutagenesis screen. It encodes a bHLH-PAS protein that is highly conserved throughout evolution. Tetrapods also have the second Clock gene, Clock2 or Npas2 (Neuronal PAS domain protein 2). Conversely, the fruit fly, an invertebrate, has only one clock gene. Interrogation of the five teleost fish genome databases revealed that the zebrafish and the Japanese pufferfish (fugu) each have three clock genes, whereas the green spotted pufferfish (tetraodon), the Japanese medaka fish and the three-spine stickleback each have two clock genes. Phylogenetic and splice site analyses indicated that zebrafish and fugu each have two clock1 genes, clock1a and clock1b and one clock2; tetraodon also have clock1a and clock1b but do not have clock2; and medaka and stickleback each have clock1b and one clock2. Genome neighborhood analysis further showed that clock1a/clock1b in zebrafish, fugu and tetraodon is an ancient duplicate. While the dN/dS ratios of these three fish clock duplicates are all <1, indicating that purifying selection has acted upon them; the Tajima relative rate test showed that all three fish clock duplicates have asymmetric evolutionary rates, implicating that one of these duplicates have been under positive selection or relaxed functional constraint. These results support the view that teleost fish clock genes were generated from an ancient genome-wide duplication, and differential gene loss after the duplication resulted in retention of different ancient duplicates in different teleost fishes, which could have contributed to the evolution of the distinct fish circadian clock mechanisms.