Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

Yukuto Sato12*, Yasuyuki Hashiguchi1 and Mutsumi Nishida1*

Author Affiliations

1 Division of Molecular Marine Biology, Ocean Research Institute, The University of Tokyo, 1-15-1 Minamidai, Nakano-ku, Tokyo 164-8639, Japan

2 Division of Population Genetics, National Institute of Genetics, 1111 Yata, Mishima, Shizuoka 411-8540, Japan

For all author emails, please log on.

BMC Evolutionary Biology 2009, 9:127  doi:10.1186/1471-2148-9-127

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/9/127


Received:13 November 2008
Accepted:5 June 2009
Published:5 June 2009

© 2009 Sato et al; licensee 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.

Abstract

Background

Recent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.

Results

We applied a state-of-the-art method of maximum-likelihood phylogenetic inference and conserved synteny analyses to each of 130 genes involved in the above biological systems of human. These analyses identified 116 orthologous gene groups between teleosts and tetrapods, and 45 pairs of 3R-WGD-derived duplicate genes among them. This suggests that more than half [(45×2)/(116+45)] = 56.5%) of the loci, probably more than ten thousand genes, present in a common ancestor of the four teleosts were still duplicated after the 3R-WGD. The estimated temporal pattern of gene loss suggested that, after the 3R-WGD, many (71/116) of the duplicated genes were rapidly lost during the initial 75 million years (MY), whereas on average more than half (27.3/45) of the duplicated genes remaining in the ancestor of the four teleosts (45/116) have persisted for about 275 MY. The 3R-WGD-derived duplicates that have persisted for a long evolutionary periods of time had significantly larger number of interacting partners and longer length of protein coding sequence, implying that they tend to be more multifunctional than the singletons after the 3R-WGD.

Conclusion

We have shown firstly the temporal pattern of gene loss process after 3R-WGD on the basis of teleost phylogeny and divergence time frameworks. The 3R-WGD-derived duplicates have not undergone constant exponential decay, suggesting that selection favoured the long-term persistence of a subset of duplicates that tend to be multi-functional. On the basis of these results obtained from the analysis of 116 orthologous gene groups, we propose that more than ten thousand of 3R-WGD-derived duplicates have experienced lineage-specific evolution, that is, the differential sub-/neo-functionalization or secondary loss between lineages, and contributed to teleost diversity.

Background

Whole-genome duplication (WGD) is believed to be one of several major evolutionary events that shaped the genomes of eukaryotes from yeasts and plants to vertebrates [1-6]. WGD, which generates dozens of thousands of duplicate genes, is usually followed by massive gene loss or the acquisition of new roles for the duplicated genes (i.e., sub-/neo-functionalization [7,8]). The evolution of these duplicate genes should occur independently among the lineages that diverged after the WGD. Such an event may contribute to the emergence of genomic differences among lineages that have experienced a WGD, and, therefore increase the genetic and phenotypic diversity among the organisms in that group.

Jawed vertebrates, which have the most complex body plan and behavioral characteristics, are thought to have experienced two rounds of (1R- and 2R-) WGD events early in their evolution, and teleost fishes experienced one more WGD (3R-WGD [6,9-12]). This notion is supported by data from several recent genomic analyses [4,13-21]. Since these events (1R-, 2R-, and 3R-WGD) occurred in ancestors shared by major phylogenetic groups such as tetrapods and teleosts, they may have been important for the formation of vertebrate-specific genomic features.

The teleost-specific 3R-WGD may be crucial to understanding certain aspects of teleost diversity; in addition, the event is very interesting given its time of occurrence. The 1R- and 2R-WGDs, which contributed to the formation of jawed-vertebrate genomes, including those of mammals, occurred before the split between tetrapods and teleosts; thus, they are very old events [5,14,15,17]. Probably because of their oldness, the 1R- and 2R-WGDs seem to have left few traces in the duplicated genes or tree topology of the gene families within the current genomes [4,22,23]. On the other hand, the 3R-WGD is estimated to have occurred in an ancestor of teleosts but after the divergence of teleosts and tetrapods [24-26]. Thus, it is the relatively recent WGD shared by a large vertebrate group, i.e., teleosts. Therefore, we can expect that teleost genomes contain many more WGD-derived duplicate genes and their 'traces' of evolution than tetrapod genomes. In addition, in teleosts, whole-genome sequence data from multiple species [16,20,27], reliable phylogenetic frameworks, and estimated divergence times between lineages [28-33] are now available. Thus, we can systematically analyze the evolution of WGD-derived duplicate genes and their lineage-specific actions by focusing on teleosts, which experienced 3R-WGD. This will provide valuable insights into the evolution of vertebrates through WGDs.

Previous studies have analyzed WGD-derived duplicate genes mainly by pairwise comparisons using the genomes of phylogenetically distant species such as human and pufferfishes [15-17,19,20]. Studies focusing on a particular gene family have also been performed (e.g., [34,35]). Studies on the medaka and green spotted puffer fish (Tetraodon) genomes suggested that they contain about 2,000 pairs of duplicate genes derived from the 3R-WGD [20]. The analyses based on pairwise comparison, however, are insufficient to address the detailed evolutionary process, e.g., the temporal- and lineage-specific manner of gene loss/persistence, after the WGD. These concerns are still not fully resolved by studies based on whole genome data from a few teleost species [36,37], and they can be resolved only by ancestral-state inferences by using multiple genome data, and phylogenetic frameworks and divergence time estimations. Thus, it is unknown how duplicated genes have been maintained or lost temporarily through lineage diversification after the WGD and how many of them contribute to the current genomes of each species.

To address these concerns, we estimated the evolutionary processes experienced by the duplicated genes after the 3R-WGD by analyzing the whole-genome data of four teleosts (zebrafish, stickleback, medaka, and Tetraodon) based on teleost phylogeny and their divergence time estimates. In our comparative analyses, we focused on signal transduction pathways involved in learning, memory, and sensory perception, which may have played crucial roles in vertebrate evolution, as well as energy metabolism, which is common to eukaryotes. As representatives of each pathway, we selected four molecular interaction networks: long-term potentiation of synaptic transmission (LTP), taste transduction (TT), olfactory transduction (OT), and tricarboxylic acid cycle (TCA). First, we sought duplicate genes derived from the 3R-WGD by intensive phylogenetic and conserved synteny analyses of the genes involved in the above networks. Next, we estimated the temporal pattern of duplicate gene loss/persistence after the 3R-WGD on the basis of teleost phylogeny and their divergence time framework [28,30,31,33]. The 3R-WGD-derived duplicates that have persisted over the course of evolution were analyzed in terms of their function, locations in the network topologies, numbers of interacting partners, and total length of coding sequences.

Results

Orthologous gene groups between tetrapods and teleosts

In this study, protein-coding genes involved in LTP, TT, OT, and TCA were analyzed to search for duplicate genes derived from the 3R-WGD, and their subsequent loss or retention was inferred. According to the KEGG pathway database [38], LTP, TT, OT, and TCA (the network schemes are based on knowledge in the human) comprise 67, 24, 32, and 27 human loci, respectively (olfactory receptor [OR] and taste receptor type 2 [T2R] genes were excluded; see Methods). Among these genes, three, eight, and five were repeatedly involved in LTP and TT; LTP and OT; and LTP, TT, and OT, respectively. After removing these overlaps, 130 human genes [see Additional file 1: Table S1] were subjected to a comparative genomic analysis with other animal genomes.

Additional file 1. Supplementary tables. This PDF file includes supplementary tables S1 and S2.

Format: PDF Size: 368KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Putative orthologs of these 130 human genes in Tetraodon, stickleback, medaka, zebrafish, chicken, clawed frog, ascidian, and fruit fly were obtained through BLASTN searches against their genomes using human protein-coding sequences as queries (for details, see Methods). The obtained genes were subjected to a series of phylogenetic analyses, including preliminary and secondary neighbor-joining (NJ) analyses and final maximum likelihood (ML) analysis. Based on the ML trees [see Additional file 2: Fig. S1–S63], the genes were classified into four categories (Figure 1; for details, see Methods). As a result, 119 of the 130 human genes were found to have orthologs in teleost genomes [see Additional file 2: Fig. S1–S63], while no clear orthologs were identified for the remaining 11 genes (PPP3CB; PPP3R2; CALM1, 2, 3, and 6; PRKCG; GUCA1C; and CLCA1, 2, and 4) [see Additional file 2: Fig. S14, S15, S16, S25, S45, and S47]. Among the 119 human genes, 3 were duplicated specifically in human [see Additional file 2: Fig. S7 and S53B]. Accordingly, we identified a total of 116 orthologous gene groups between tetrapods and teleosts, corresponding to the 130 human genes examined (Table 1).

Additional file 2. Supplementary figures. This PDF file includes supplementary figures S1 – S68.

Format: PDF Size: 983KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 1. Schematic view of the four categories of orthologous gene groups. (A) Gene groups with a 1:1 orthologous relationship between tetrapods and teleost fishes. (B) Gene groups that were duplicated specifically within the human or tetrapod lineage. (C) Gene groups that were duplicated through 3R-WGD. (D) Gene groups with no clear orthologous relationship between tetrapods and teleosts. Arrows denote gene duplication events. Triangles denote the orthologous gene clades derived from tetrapods, teleosts, or outgroup organisms (Ciona and Drosophila). Shaded triangles denote orthologous gene clades containing a human networks-related gene.

Table 1. Number of orthologous relationships identified between human and teleosts.

Persistence of the duplicated genes after the 3R-WGD

Systematic analysis of the gene phylogenies described above successfully identified the duplicated genes derived from the 3R-WGD. These identifications were ascertained based on previously published information related to doubly conserved synteny in the medaka genome [20] shown in [Additional file 2: Fig. S1–S63] and the genome-to-genome conserved synteny information presented in [Additional file 2: Fig. S64–S68]. Duplicated genes derived from the 3R-WGD were identified in at least one species in 45 of the 116 orthologous gene groups (Table 2). This indicates that the genome of the common ancestor of the four teleosts (zebrafish, medaka, stickleback, and Tetraodon) contained at least 161 (116+45) loci belonging to the 116 orthologous gene groups, meaning that 56.5% [(45×2)/161] of the loci in the ancestor were duplicate genes derived from the 3R-WGD.

Table 2. Number of network-related orthologous gene groups in which gene duplication by the 3R-WGD was detected.

Figure 2 shows a comprehensive graphical representation of the results of our data-mining, phylogenetic, and synteny analyses. Circles indicate the presence of particular genes in the genome of each species, providing information related to lineage-specific loss or the persistence of the 3R-WGD-derived duplicates within the four teleost genomes. For instance, the zebrafish genome is missing one gene for PPP1CA in its LTP network, whereas the other fishes possess both members of the 3R-WGD-derived pair. Based on the lineage-specific presence or absence of the duplicates, we counted the number of 3R-WGD-derived genes and the total number of loci in the four teleost genomes (Table 3). The average of these numbers across the four teleosts shows that, among 133 loci still present on average in a teleost genome, average 55 (55/133 = 41.4%) remain duplicated from the 3R-WGD in at least one teleost genome (not in all four teleost genomes).

thumbnailFigure 2. Schematic representation of the results of comparative genomics-based data mining and maximum likelihood phylogenetic inferences for each gene. The magenta and green circles denote the presence of genes belonging to the orthologous gene groups in which 3R-WGD-derived duplicates were detected. Gray circles indicate the presence of singleton genes. Open circles denote loci that were specifically duplicated within the species. Question marks indicate the presence of partial gene sequences that are phylogenetically unsorted. Stars denote those gene groups involved in more than one network. Abbreviations: LPT, long-term potentiation; TT, taste transduction; OT, olfactory transduction; TCA, tricarboxylic acid cycle; 3R-WGD, third-round whole genome duplication.

Table 3. Total number of network-related loci and 3R-WGD-derived duplicated loci identified in each teleost genome.

A few genes were found to have expanded dramatically by repeated duplication independent of WGD. For example, in the OT system, phosphodiesterase 1C (PDE1C) genes were duplicated specifically in the stickleback genome [see Additional file 2: Fig. S44], implying the existence of multiple bypass circuits that regulate the cyclic AMP (cAMP) pathway in the stickleback OT system. Other examples include the calmodulin (CaM or CALM) genes involved in LTP and OT. The CaM gene family included diverse members before the split of tetrapods and teleosts, and the differential members have persisted in the tetrapod and teleost genomes, respectively [see Additional file 2: Fig. S16]. Interestingly, the amino acid sequences of CALM1 and 2 in tetrapods and their closest CALMs in teleosts are identical in their alignable region.

Temporal loss or persistence of the duplicated genes after the 3R-WGD

Based on our data concerning the presence or absence of 3R-WGD-derived duplicate loci (Figure 2), we estimated the number of gene loss events on the basis of the teleost phylogeny using parsimony (for details, see Methods). The inferred gene loss events were then assigned to the phylogenetic tree with branch lengths proportional to the estimated divergence time derived from a molecular clock analysis of mitochondrial genome sequences [33] (Figure 3A).

thumbnailFigure 3. Inferred temporal process of gene loss after 3R-WGD in a teleost ancestor. (A) The estimated numbers of gene loss events in the teleost phylogeny was parsimoniously drawn from the presence or absence of loci belonging to the 116 orthologous gene groups. The branch lengths are proportional to the estimated divergence time among the lineages [33]. Gray bars indicate 95% confidence intervals for them. (B) Approximation of the number of gene loss events based on a neutral model of the loss-of-function of duplicated genes (αe-2μt) [8,65,66]. Gray and black lines show the approximation from 3R-WGD and point c, and all data points, respectively. (C) Proposed phase transition during duplicate gene loss. All data points except the 3R-WGD were approximated by a linear equation. (D) Approximation of the number of gene loss events to the upper and lower extremes of the 95% confidence intervals of estimated occurrence time of the 3R-WGD [15,17,19] and divergence times of teleosts [33]. The approximations to the upper and lower extremes are shown by black and gray lines, respectively.

The obtained temporal pattern of gene loss (Figure 3B) indicates that many (71/116 genes; see Table 2) of the duplicated genes were rapidly lost following the 3R-WGD before the split of zebrafish (node c), whereas on average more than half (27.3/45 genes; see Table 3 "Grand total": 54.5/2 = 27.3) of the duplicated genes remaining in the ancestor of the four teleosts (45/116 genes) have persisted for about 275 million years (MY). Least-squares fitting of a neutral model of the loss of gene function to the data points for the 3R-WGD and the common ancestor yielded an exponential decay curve with a slope of -0.0123 (Figure 3B gray line). On the other hand, approximation to all data points yielded a moderate curve with slope = -0.0034 (Figure 3B black line), showing that the tempo of gene loss was not constant over the course of evolution after the 3R-WGD (Figure 3C). This result was essentially unchanged when we used upper and lower extremes of the 95% confidence intervals of the estimated times of 3R-WGD and lineage divergences in the approximation (relevant slopes were -0.0106 and -0.0030, and -0.0135 and -0.0037 when upper and lower extreme values were used, respectively; see Figure 3D).

Properties of the gene groups that retained 3R-WGD-derived duplicates

We characterized the 3R-WGD-derived duplicate genes on the basis of several properties of the encoded proteins. In terms of their locations in the network topologies of LTP, TT, OT, and TCA, we were not able to find a local concentration of 3R-WGD-derived duplicates (e.g., their localization at upstream or downstream of the pathways; Figure 4). In addition, although most of the proteins included in the four systems are of different types, we found no significant difference in the frequency of 3R-WGD-derived duplicates (described in Table 2) among the four networks (χ2 = 1.7266, d.f. = 3, P = 0.6310).

thumbnailFigure 4. Distributions of the orthologous gene groups containing duplicate genes that arose through 3R-WGD (shown as gray boxes) in the network diagrams. Panels A-D: long-term potentiation (LTP), taste transduction (TT), olfactory transduction (OT), and tricarboxylic acid cycle (TCA), respectively. White boxes denote orthologous gene groups with no 3R-WGD-derived duplicates (1:1 orthologous relationship between tetrapods and teleosts). Dashed boxes denote gene groups in which no clear orthologous relationship was identified between tetrapods and teleosts. The network schemes of the LTP, TT, OT, and TCA are based on knowledge in human summarized in KEGG pathway database [38].

We next compared the function, total length, and number of interacting partners of the proteins in the networks. There was no significant difference in the frequency of each type of function (e.g., enzymes, G proteins, ion channels, phosphorylation enzymes, and receptors) between the 3R-WGD-derived duplicate genes and 1:1 orthologous genes (χ2 = 4.3984, d.f. = 5, P = 0.4936; Figure 5A). However, there was a significant difference in the length of the protein-coding sequences. The protein peptides encoded by the 3R-WGD-derived duplicate genes tended to be long (>1000 amino acids) rather than short (<200 amino acids; χ2 = 8.5044, d.f. = 2, P = 0.0142; Figure 5B; human protein were used as the hypothetical ancestral state). Additionally, the 3R-WGD-derived duplicates were enriched in longer-sized glutamate receptors (GRM1, GRM5, GRIN2A, GRIN2B, and GRIN2C) compared to the 1:1 orthologous genes, although this trend was not statistically significant (χ2 = 2.2011, d.f. = 2, P = 0.3327). The average number of interacting partners was also significantly higher for the 3R-WGD-derived duplicate genes (n1 = 71, n2 = 45, Welch's t = 2.0203, P = 0.0470) [see Additional file 1: Table S2]. Within the 3R-WGD-derived duplicates, there were no significant differences in the length of protein-coding sequences and the number of interacting partners among the genes that remain duplicated in one, two, three, and four teleost genomes (see Table 2; data not shown).

thumbnailFigure 5. Characteristics of the orthologous gene groups in which duplicated genes generated by the 3R-WGD were detected (black bar) compared to the groups in which 3R-WGD-derived duplicate genes were not detected (white bar). (A) Frequencies of six types of protein function. (B) Frequency of each of three classes of protein-coding region length (# of amino acids).

Discussion

Persistence of the duplicated genes after the WGD event

Careful analysis of more than 100 orthologous gene groups in the tetrapod and teleost genomes suggested that the genome of a common ancestor of four teleosts (zebrafish, medaka, stickleback, and Tetraodon) contained a high proportion (90/161 = 56.5%) of duplicated genes derived from the teleost-specific 3R-WGD (Table 2). This estimate was produced by extensive data mining of whole-genome databases and reliable ML phylogenetic analyses of the genes involved in LTP, TT, OT, and TCA. Although these genes represent small fraction of the whole genome, and such transcriptional, signal-transduction, and metabolic networks may have different biases for duplicate gene preservation [21,37], we found no significant difference in the frequency of 3R-WGD-derived duplicates among the four networks (see Results and Table 2). This may imply that there is no remarkable difference in retention rate of the 3R-WGD-derived duplicates between subcellular networks that would be important in vertebrate evolution, and metabolism networks that are common to eukaryotes. A possible reason for this is the relatively smaller population size of vertebrates, in which the neutral, random genetic drift is more dominant factor in their duplicate gene loss/retention [8]. In the monocellular eukaryotes, which have much larger population sizes, on the other hand, natural selection would play a more dominant role in duplicate gene loss/retention. Actually, in a ciliate Paramecium study, it is suggested that gene dosage constraints (a form of natural selection) after WGD favoured co-retention of the genes within same metabolic pathways or protein complexes [39].

If the retention ratio of 56.5% represents a trend in the whole genome, the duplicate genes generated by the 3R-WGD may have occupied more than half of the protein-coding loci in a common ancestor of the four teleosts. This may even underestimate the true proportion of the 3R-WGD-derived duplicates in the common ancestor, since some duplicates may have been lost independently in all 4 species. In any case, such a large number of duplicate loci, probably more than ten thousand genes, may have undergone lineage-specific loss or differential functionalization among teleosts after their diversification. Such lineage-specific gene evolution would characterize the genome of each teleost lineage, and consequently, contribute to various diversity among modern teleosts, which include more than 26,000 species, 500 families, and 40 orders [40].

Our results also suggests that many of the 3R-WGD-derived duplicate genes persist in modern teleost genomes. The number of the duplicates identified (Table 3) suggests that they comprise on average 41.4% of the protein-coding loci in modern teleosts. This value is higher than those reported in previous studies based on pairwise comparisons of the human and teleost genomes, i.e., 15.2% [(2134×2)/28005] to 20.0% [(2009×2)/20131] [20,41]. These estimations are based on analysis of much larger number of genes compared to those of this study, however, these might be underestimations for the retention rate of 3R-WGD-derived duplicates due to the limited sensitivity of semi-automatic BLAST-based analyses. In our study, we simultaneously analyzed multiple teleost genomes using state-of-the-art methods for phylogenetic inference. Thus, we may have more effectively identified the 3R-WGD-derived duplicates, but instead, our analysis covers limited fraction of the genomes (116 orthologous genes; about 0.5% of all protein-coding genes).

Although teleost genomes may contain large number of 3R-WGD-derived duplicate genes as suggested above, it is reported that teleost and tetrapod genomes have similar number of protein-coding genes in general (mammals: 22,000 on average in human, mouse, dog, and cow; teleosts: 23,000 on average in Tetraodon, stickleback, medaka, and zebrafish) [41]. This may be explained by both an increase of gene number in tetrapods through repeated chromosomal rearrangements and local gene duplications [16,20], and secondary loss of 1R- and 2R-WGD-derived duplicates in teleosts [37]. These issues should be addressed within the broader context of 1R/2R/3R-WGDs in early vertebrate evolution based on the emerging genome sequence data of shark (elephantfish), lamprey, lancelet, and sea urchin in the future.

Evolution of the duplicate genes derived from the 3R-WGD

The inferred temporal pattern of loss/persistence of the 3R-WGD-derived duplicates (Figure 3) provides new insights into the long-term evolution of duplicate genes after WGD. Whereas many (71/116 genes) of the duplicated genes are estimated to have been lost during the initial 75 MY before the divergence of Otocephala (the group including zebrafish; Figure 3B gray line), on average more than half (27.3/45 genes) of the genes survived until the divergence of Otocephala (45/116 genes) persisted for 275 MY (Figure 3B black line). Such a temporal pattern of duplicate gene maintenance cannot be solely explained by a complete neutral loss of gene function, which would be expected to lead to an exponential drop in gene number over time (Figure 3B gray line). Instead, it appears that the genes which persisted for 275 MY have been maintained by selection. This evolutionary scenario is essentially unchanged if we use both upper and lower extremes of 95% confidence intervals for estimated times of the 3R-WGD and lineage divergences (see Figure 3D). Our results, however, are based on a subset of genes in the genomes of only four teleosts and an approximate divergence times of teleosts with large error intervals of 50–100 MYA [33]. In addition, our results might be affected by incompleteness of the draft genome data of teleosts, especially that of zebrafish. Therefore, extrapolating our model proposed here to genome evolution in 26,000 teleost species over 300 MYA should be verified in further studies with an increased number of sample genes and genomes.

The long-term maintenance of the duplicates may be attributable to sub-/neo-functionalization of their gene functions [8,42-44] in general, although advantageous effect of the increase of gene dosage by duplication may also contribute to the duplicate gene persistence [45]. The model of subfunctionalization predicts that, initially after the WGD, many of the duplicated genes were rapidly lost due to the neutral accumulation of degenerative mutations leading to a loss of gene function (neutral phase; Figure 3C); however, a subset of duplicate genes were maintained through subfunctionalization caused by degenerative mutations leading to the loss of a subset of gene functions. Such subfunctionalized duplicates can persist for extremely long (evolutionary) periods of time by selection (selection phase; Figure 3C) and thus afford opportunities for the evolution of new genes with novel functions or more adaptive properties (sub-neofunctionalization [26,43,44]). The proposed time course of gene loss/persistence after the 3R-WGD (Figure 3) seems to be well-matched with the model of sub-neofunctionalization.

The above view is supported by the present results of several analyses related to protein multifunctionality. Because the initial subfunctionalization occurs primarily through degenerative mutations leading to the loss of a subset of gene functions, it is theoretically possible that the probability of subfunctionalization depends on the degree of multifunctionality of the genes or their encoded proteins [8]. This prediction is supported by the observation that the genes included among the 3R-WGD-derived duplicates are significantly longer in terms of their coding sequence (longer sequences likely contain more protein domains and motifs; Figure 5B) and have a significantly larger total number of interacting partners [see Additional file 1: Table S2]. These results, which are compatible with the theoretical prediction of the subfunctionalization model, imply that the 3R-WGD-derived duplicates that remain in modern teleosts have maintained initially by subfunctionalization.

The timing of the occurrences of sub-neofunctionalization is important for interpreting our results and understanding long-term genome evolution after WGDs. If the sub-neofunctionalization have completed mainly in the initial phase after the 3R-WGD but before the divergence of the four teleosts (see Figure 3), this would have contribute little to the differential gene functions among teleost lineages. However, if the neofunctionalization following subfunctionalization occurred mainly after the divergence of the four teleosts, such lineage-specific sub-neofunctionalization of many 3R-WGD-derived duplicates would have contributed to various aspects of teleost diversity. We propose that the latter is the case in general, because neofunctionalization depends on rare beneficial mutations [8] and, thus, this process will proceed gradually over long evolutionary periods of time [26]. These issues, however, should be addressed by further analyses on the spatio-temporal patterns of gene expression and/or repertoires of protein function of the respective 3R-WGD-derived duplicates in the future.

Genome and gene duplication and teleost evolution

The present study clearly shows that a few gene groups among those analyzed here have expanded dramatically by repeated gene duplications independent of WGD. One representative of them is the well known OR (olfactory receptor) gene family. In vertebrates, this family is constituted by about ten subfamilies [46], and the different subfamily has expanded dramatically in each genome of tetrapods and teleosts, respectively [46,47]. Furthermore, in the OT system, PDE1C genes were found to have been duplicated specifically in stickleback genome [see Additional file 2: Fig. S44]. PDE is a key enzyme in the cyclic AMP pathway, as it regulates the localization, duration, and amplitude of intracellular cAMP signaling [48]. Sticklebacks, therefore, may have multiple bypass circuits to regulate the cAMP pathway possibly in the OT system via sub-/neo-functionalization of these highly duplicated PDE1C genes. These multiple occurrences of the PDE1C genes may be a genetic basis for some adaptive traits of sticklebacks. Recent analysis on the stickleback PDE1C genes has favoured the notion that these multiple duplicates have been retained, not by sub-/neo-functionalization in their gene functions, but rather through their effect of increased gene dosage within the OT system, which is possibly correlated with the territorial behaviour of stickleback [49]. In summary, the above results suggest that both of WGD and lineage-specific expansion of particular genes generated genomic characteristics of each species.

It has been suggested that the 3R-WGD promoted speciation and thus led to an increase in species diversity among teleosts [12,50-54]; however, the association between WGD and species richness in teleosts remains controversial [55]. Even if the 3R-WGD is less directly associated with species richness, we propose that the lineage-specific evolution of many 3R-WGD-derived duplicates, that is, the differential sub-neofunctionalization or secondary loss between lineages, contributed to various aspects of genomic diversity that exists among teleosts. Such diversity may underlie the adaptive radiation of teleosts that inhabit a broad range of ecological zones – from marine to fresh water, from deep to shallow sea environments, and from equatorial to polar regions. To understand the diversity and evolution of teleosts, as well as vertebrates in general, it will be particularly valuable to characterize the differences between teleosts and the basal non-teleost fishes, which have not experienced the 3R-WGD [24-26], in relation to the genetic impact of the 3R-WGD.

Conclusion

More than 100 genes was systematically analyzed based on whole-genome sequences of multiple teleost fish species (zebrafish, medaka, stickleback, and Tetraodon) and state-of-art methods of phylogenetic inference. This revealed that, while many of the duplicate genes derived from the teleost-specific 3R-WGD has been rapidly lost after the 3R-WGD within about 75 MY, the remaining duplicates may have occupied more than half (56.5%) of the protein-coding loci in a common ancestor of the four teleosts. The following temporal pattern of loss/persistence of these duplicates is understood by sub-neofunctionalization model of duplicate gene evolution, and this view is supported by analysis of length of protein-coding sequences and numbers of interaction partner of proteins, which are likely associated with gene multifunctionality. Based on these results, derived from the first systematic assessment of evolutionary fate of more than 100 genes after the 3R-WGD, we propose that many of the 3R-WGD-derived duplicates, probably more than ten thousand of genes, have undergone lineage-specific evolution or secondary loss, and contributed to teleost diversity. Lineage-specific expansion of some portion of gene families also appears to have contributed to genomic characteristics of each lineages.

Methods

Identification of orthologous gene groups between tetrapods and teleosts

Network diagrams for LTP, TT, OT, and TCA, and the coding sequences of the human genes that comprise these networks were obtained from the KEGG pathway database [38]. The obtained human coding sequences were used as queries for a BLASTN search against the Ensembl genome database [27,41]. The following versions of the Ensembl genome database were used: human (Homo sapiens, NCBI 36, October 2005), chicken (Gallus gallus, WASHUC2, May 2006), clawed frog (Xenopus tropicalis, JGI 4.1, August 2005), zebrafish (Danio rerio, Zv7, April 2007), medaka (Oryzias latipes, HdrR, October 2005), stickleback (Gasterosteus aculeatus, BROAD S1, February 2006), green spotted puffer fish (Tetraodon nigroviridis, TETRAODON 7, April 2003), ascidian (Ciona intestinalis, JGI 2, March 2005), and fruit fly (Drosophila melanogaster, BDGP 4.3, July 2005). Since the two pufferfishes Fugu rubripes and T. nigroviridis are much closer to each other than to the other teleosts above, only the latter, for which the genome sequence was determined at relatively high coverage [16], was analyzed. The resulting BLAST hits were manually screened (E-value cut-off of < 10-3) and evaluated for their gene product length and Ensembl annotations to confirm their similarity to the human gene queries. When only a partial sequence was found in the Ensembl database, we predicted the full-length coding sequence from the genomic sequence using the program WISE2 [56].

The primary sequences of the proteins obtained by the above procedure were aligned using ClustalW [57]. All gap-containing sites were removed. For each alignment, a preliminary NJ analysis was performed based on Poisson-corrected genetic distances using MEGA 3.1 [58]. Based on the resultant NJ trees, those sequences comprising clades, which were apparently distinct from any human query gene, were excluded. The remaining sequences were again analyzed by NJ to identify ingroup and outgroup sequences. The selected primary sequences were then re-aligned and subjected to ML phylogenetic analysis using the program TreeFinder (version June 2007) [59,60] with 1,000 LR-ELW (the expected-likelihood weights applied to local rearrangements of tree topology) edge support tests [61]. In almost all cases in this analysis, the obtained tree topology was consistent between ML and second-ML tree with slight changes in the ML score and supporting values for the nodes of the tree. When the final ML tree was ambiguous, further ML analysis was performed based on nucleotide sequences. The phylogenetic information contained within nucleotide sequences is generally greater than that contained within amino acid sequences, and, in many cases, it increases the resolution of the resulting phylogeny. Such analysis require significant time and effort, and therefore could not be applied to all analyses undertaken in this study. The best-fitting amino acid and nucleotide substitution models were selected using programs ProtTest 1.4 [62] and ModelTest 3.06 [63], respectively.

Based on the final ML trees [see Additional file 2: Fig. S1–S63], the genes were classified into four categories: (i) genes with a 1:1 orthologous relationship [64] between human and teleosts (Figure 1A); (ii) genes that were duplicated in human or tetrapods but not in teleosts (Figure 1B); (iii) 3R-WGD-derived duplicates (Figure 1C); this identification was confirmed by data showing doubly conserved synteny in the medaka genome [20] [indicated in Additional file 2: Fig. S1–S63] and genome-to-genome conserved synteny analyses mentioned below [shown in Additional file 2: Fig. S64–S68]; and (iv) genes with no clear orthologous relationship between human and teleosts (Figure 1D). A full list of the genes used in our final ML analysis, including their Ensembl IDs, is given in [Additional file 3: appendix].

Additional file 3. Supplementary appendix. This PDF file includes an appendix.

Format: PDF Size: 384KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Synteny analysis

To confirm whether the putative 3R-WGD-derived duplicate genes were actually generated by the 3R-WGD, we analyzed conserved syntenies in the medaka genome. In Supplementary Table 15 of the report showing the draft genomic sequence of medaka [20], the authors presented information on conserved syntenies in the medaka genome. Specifically, they determined whether those medaka genes corresponding to 20,352 human protein-coding genes were located in a block of doubly conserved synteny (DCS) derived from the 3R-WGD. Using this information, we confirmed that the inferred teleost-specific duplicate genes were derived from the 3R-WGD. The results of this validation were shown in supplementary figures [see Additional file 2] (the data of gene locations are not shown).

For those teleost-specific duplicate genes that appeared to be derived from the 3R-WGD but were not located on DCSs in medaka, we analyzed the genomic regions around the loci in human, zebrafish, medaka, stickleback, and Tetraodon. Physical mapping data in the neighborhood of each locus were obtained via queries and BLAST searches (E-value threshold < 10-3) using the orthologous prediction section of the Ensembl genome database [27,41]. We picked up all identifiable genes described as putative orthologs of the queries. Their genomic locations were then used to rebuild the synteny maps [see Additional file 2: Fig. S64–S68].

Analysis of the temporal pattern of gene loss

The numbers of gene loss events after the 3R-WGD was estimated using parsimony based on the teleost phylogeny [28,30,31]. The inferred total number of 3R-WGD-derived duplicate pairs that persisted in each node of the phylogeny was matched with the corresponding divergence time estimate [33]. The derived two-dimensional data points, estimated occurrence time of the 3R-WGD event (average 353.0 MY ago) [15,17,19], and four data points representing the current state of the four teleost genomes were approximated to a neutral model of the loss-of-function of duplicated genes, αe-2μt [8,65,66], using least-squares fitting. The parameters μ, t, and α corresponds to the null mutation rate, number of generations since the time of duplication, and normalization factor, respectively. In this analysis, absolute time (years) was used to approximate the number of generations.

Estimation of the number of interacting partner of the investigated proteins

The number of interacting partners for the proteins in the LTP, TT, OT, and TCA was counted on the basis of human network diagrams obtained from the KEGG pathway database [38]. In the case of CAMK2B, CAMK2G, PRKACB, CAMK2A, CAMK2D, ADCY8, PLCB2, ITPR3, PRKACA, PRKACG, PRKX, PRKY, CaM, and PKA, those molecular interactions that served essentially the same function were counted as one. For example, adenylate cyclase 8 (ADCY8, EC: 4.6.1.1) is involved in LTP and TT, and it interacts with CaM and cAMP in the former, and with cAMP and G-protein in the latter; however, the total number of interacting partners for ADCY8 was counted as three (CaM, cAMP, and G-proteins). In the case of CaM, although CaM interacts with many of types of proteins, only two interacting partners were counted (Ca 2+ and another protein). This is because Ca 2+-activated CaM functions as a calcium sensor and signal transducer for a multitude of proteins. This view is supported by the fact that the amino acid sequences of CaM are extremely conserved among vertebrates [67] while the interacting partners of CaM are numerous. By similar reasoning, the number of interacting partners for PKA was counted as two (cAMP and another protein). All OR genes among the OT and T2R genes in the TT network were excluded from the analysis because the OR and T2R gene families have been enormously expanded by specific duplications in mammals, and the orthologous relationships between tetrapods and teleosts are less clear [46,47,68]. If large numbers of OR and T2R loci were incorporated into the current analysis, the results would be biased and beyond the scope of this study. The numbers loci of the taste receptor type 1 genes (T1R1, T1R2, and T1R3) are as previously reported [69,70].

List of abbreviations used

3R-WGD: third-round whole genome duplication; LTP: long-term potentiation; TT: taste transduction; OT: olfactory transduction; TCA: tricarboxylic acid cycle; OR: olfactory receptor; T2R: taste receptor type 2; NJ: neighbour joining; ML: maximum likelihood; PPP: protein phosphatase; CALM (gene name) and CaM (protein name): calmodulin; PRK: protein kinase; GUCA: guanylate cyclase activator; CLCA: calcium-dependent chloride channel; MY: million years; cAMP: cyclic adenosine monophosphate; DCS: doubly conserved synteny; CAMK: calcium/calmodulin-dependent protein kinase; ADCY: adenylate cyclase; PLCB: phospholipase C beta; ITPR: inositol 1,4,5-triphosphate receptor; T1R: taste receptor type 1; PKA: protein kinase A: LR-ELW, the expected-likelihood weights applied to local rearrangements of tree topology.

Authors' contributions

YS, YH, and MN designed the study. YS carried out the analyses and drafted the manuscript. MN participated in coordination and helped to draft the manuscript. All authors read and approved the final version of the manuscript.

Acknowledgements

We are grateful to the three anonymous reviewers for their constructive comments for improving earlier versions of the manuscript. The zebrafish, Tetraodon, stickleback, and medaka sequence data were produced by the Sanger Institute, Genoscope and the Broad Institute, the Broad Institute, and the National Institute of Genetics, Japan, respectively. The frog, chicken, human sequences were produced by the Joint Genome Institute (JGI), the Genome Sequencing Center at Washington University, St Louis, and the International Human Genome Sequencing Consortium, respectively. The fruitfly sequence data was provided by the Berkeley Drosophila Genome Project, FlyBase, and the Drosophila Heterochromatin Genome Project. The Ciona intestinalis sequence was produced by the JGI. This study was partially supported by Grants-in-Aid from the Japan Society for the Promotion of Science to M.N.

References

  1. Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution.

    Plant Cell 2004, 16:1679-1691. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Kellis M, Birren BW, Lander ES: Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae.

    Nature 2004, 428:617-624. PubMed Abstract | Publisher Full Text OpenURL

  3. Adams KL, Wendel JF: Polyploidy and genome evolution in plants.

    Curr Opin Plant Biol 2005, 8:135-141. PubMed Abstract | Publisher Full Text OpenURL

  4. Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate.

    PLoS Biol 2005, 3:e314. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Panopoulou G, Poustka AJ: Timing and mechanism of ancient vertebrate genome duplications – the adventure of a hypothesis.

    Trends Genet 2005, 21:559-567. PubMed Abstract | Publisher Full Text OpenURL

  6. Froschauer A, Braasch I, Volff JN: Fish genomes, comparative genomics and vertebrate evolution. [http://www.ingentaconnect.com/content/ben/cg/2006/00000007/00000001/art00005] webcite

    Curr Genomics 2006, 7:43-57. Publisher Full Text OpenURL

  7. Ohno S: Evolution by Gene Duplication. New York: Springer-Verlag; 1970. OpenURL

  8. Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations.

    Genetics 1999, 151:1531-1545. PubMed Abstract | PubMed Central Full Text OpenURL

  9. Lundin LG: Evolution of the vertebrate genome as reflected in paralogous chromosomal regions in man and the house mouse.

    Genomics 1993, 16:1-19. PubMed Abstract | Publisher Full Text OpenURL

  10. Holland PWH, Garcia-Fernandez J, Williams NA, Sidow A: Gene duplications and the origins of vertebrate development.

    Development 1994, (Suppl):125-133. OpenURL

  11. Amores A, Suzuki T, Yan YL, Pomeroy J, Singer A, Amemiya C, Postlethwait JH: Zebrafish hox clusters and vertebrate genome evolution.

    Science 1998, 282:1711-1714. PubMed Abstract | Publisher Full Text OpenURL

  12. Taylor JS, Braasch I, Frickey T, Meyer A, Peer Y: Genome duplication, a trait shared by 22000 species of ray-finned fish.

    Genome Res 2003, 13:382-390. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Abi-Rached L, Gilles A, Shiina T, Pontarotti P, Inoko H: Evidence of en bloc duplication in vertebrate genomes.

    Nat Genet 2002, 31:100-105. PubMed Abstract | Publisher Full Text OpenURL

  14. Panopoulou G, Hennig S, Groth D, Krause A, Poustka AJ, Herwig R, Vingron M, Lehrach H: New evidence for genome-wide duplications at the origin of vertebrates using an amphioxus gene set and completed animal genomes.

    Genome Res 2003, 13:1056-1066. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Christoffels A, Koh EG, Chia JM, Brenner S, Aparicio S, Venkatesh B: Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes.

    Mol Biol Evol 2004, 21:1146-1151. PubMed Abstract | Publisher Full Text OpenURL

  16. 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, Biémont 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, Guigó R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechavi M, Laudet V, Schachter V, Quétier 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-957. PubMed Abstract | Publisher Full Text OpenURL

  17. Vandepoele K, De Vos W, Taylor JS, Meyer A, Peer Y: Major events in the genome evolution of vertebrates: paranome age and size differ considerably between ray-finned fishes and land vertebrates.

    Proc Natl Acad Sci USA 2004, 101:1638-1643. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Hoegg S, Meyer A: Hox clusters as models for vertebrate genome evolution.

    Trends Genet 2005, 21:421-424. PubMed Abstract | Publisher Full Text OpenURL

  19. Christoffels A, Brenner S, Venkatesh B: Tetraodon genome analysis provides further evidence for whole-genome duplication in the ray-finned fish. Comp.

    Biochem Physiol Part D Genomics Proteomics 2006, 1:13-19. Publisher Full Text OpenURL

  20. Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, Jindo T, Kobayashi D, Shimada A, Toyoda A, Kuroki Y, Fujiyama A, Sasaki T, Shimizu A, Asakawa S, Shimizu N, Hashimoto S, Yang J, Lee Y, Matsushima K, Sugano S, Sakaizumi M, Narita T, Ohishi K, Haga S, Ohta F, Nomoto H, Nogata K, Morishita T, Endo T, Shin-I T, Takeda H, Morishita S, Kohara Y: The medaka draft genome and insights into vertebrate genome evolution.

    Nature 2007, 447:714-719. PubMed Abstract | Publisher Full Text OpenURL

  21. Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutiérrez EL, Dubchak I, Garcia-Fernàndez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, Toyoda A, Bronner-Fraser M, Fujiyama A, Holland LZ, Holland PW, Satoh N, Rokhsar DS: The amphioxus genome and the evolution of the chordate karyotype.

    Nature 2008, 453:1064-1071. PubMed Abstract | Publisher Full Text OpenURL

  22. Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, Gocayne JD, Amanatides P, Ballew RM, Huson DH, Wortman JR, Zhang Q, Kodira CD, Zheng XH, Chen L, Skupski M, Subramanian G, Thomas PD, Zhang J, Gabor Miklos GL, Nelson C, Broder S, Clark AG, Nadeau J, McKusick VA, Zinder N, Levine AJ, Roberts RJ, Simon M, Slayman C, Hunkapiller M, Bolanos R, Delcher A, Dew I, Fasulo D, Flanigan M, Florea L, Halpern A, Hannenhalli S, Kravitz S, Levy S, Mobarry C, Reinert K, Remington K, Abu-Threideh J, Beasley E, Biddick K, Bonazzi V, Brandon R, Cargill M, Chandramouliswaran I, Charlab R, Chaturvedi K, Deng Z, Di Francesco V, Dunn P, Eilbeck K, Evangelista C, Gabrielian AE, Gan W, Ge W, Gong F, Gu Z, Guan P, Heiman TJ, Higgins ME, Ji RR, Ke Z, Ketchum KA, Lai Z, Lei Y, Li Z, Li J, Liang Y, Lin X, Lu F, Merkulov GV, Milshina N, Moore HM, Naik AK, Narayan VA, Neelam B, Nusskern D, Rusch DB, Salzberg S, Shao W, Shue B, Sun J, Wang Z, Wang A, Wang X, Wang J, Wei M, Wides R, Xiao C, Yan C, Yao A, Ye J, Zhan M, Zhang W, Zhang H, Zhao Q, Zheng L, Zhong F, Zhong W, Zhu S, Zhao S, Gilbert D, Baumhueter S, Spier G, Carter C, Cravchik A, Woodage T, Ali F, An H, Awe A, Baldwin D, Baden H, Barnstead M, Barrow I, Beeson K, Busam D, Carver A, Center A, Cheng ML, Curry L, Danaher S, Davenport L, Desilets R, Dietz S, Dodson K, Doup L, Ferriera S, Garg N, Gluecksmann A, Hart B, Haynes J, Haynes C, Heiner C, Hladun S, Hostin D, Houck J, Howland T, Ibegwam C, Johnson J, Kalush F, Kline L, Koduru S, Love A, Mann F, May D, McCawley S, McIntosh T, McMullen I, Moy M, Moy L, Murphy B, Nelson K, Pfannkoch C, Pratts E, Puri V, Qureshi H, Reardon M, Rodriguez R, Rogers YH, Romblad D, Ruhfel B, Scott R, Sitter C, Smallwood M, Stewart E, Strong R, Suh E, Thomas R, Tint NN, Tse S, Vech C, Wang G, Wetter J, Williams S, Williams M, Windsor S, Winn-Deen E, Wolfe K, Zaveri J, Zaveri K, Abril JF, Guigó R, Campbell MJ, Sjolander KV, Karlak B, Kejariwal A, Mi H, Lazareva B, Hatton T, Narechania A, Diemer K, Muruganujan A, Guo N, Sato S, Bafna V, Istrail S, Lippert R, Schwartz R, Walenz B, Yooseph S, Allen D, Basu A, Baxendale J, Blick L, Caminha M, Carnes-Stine J, Caulk P, Chiang YH, Coyne M, Dahlke C, Mays A, Dombroski M, Donnelly M, Ely D, Esparham S, Fosler C, Gire H, Glanowski S, Glasser K, Glodek A, Gorokhov M, Graham K, Gropman B, Harris M, Heil J, Henderson S, Hoover J, Jennings D, Jordan C, Jordan J, Kasha J, Kagan L, Kraft C, Levitsky A, Lewis M, Liu X, Lopez J, Ma D, Majoros W, McDaniel J, Murphy S, Newman M, Nguyen T, Nguyen N, Nodell M, Pan S, Peck J, Peterson M, Rowe W, Sanders R, Scott J, Simpson M, Smith T, Sprague A, Stockwell T, Turner R, Venter E, Wang M, Wen M, Wu D, Wu M, Xia A, Zandieh A, Zhu X: The sequence of the human genome.

    Science 2001, 291:1304-1351. PubMed Abstract | Publisher Full Text OpenURL

  23. Wolfe KH: Yesterday's polyploids and the mystery of diploidization.

    Nat Rev Genet 2001, 2:333-341. PubMed Abstract | Publisher Full Text OpenURL

  24. Chiu CH, Dewar K, Wagner GP, Takahashi K, Ruddle F, Ledje C, Bartsch P, Scemama JL, Stellwag E, Fried C, Prohaska SJ, Stadler PF, Amemiya CT: Bichir hox A cluster sequence reveals surprising trends in ray-finned fish genomic evolution.

    Genome Res 2004, 14:11-17. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. Sato Y, Nishida M: Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites.

    BMC Evol Biol 2007, 7:204. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  27. Flicek P, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Eyre T, Fitzgerald S, Fernandez-Banet J, Gräf S, Haider S, Hammond M, Holland R, Howe KL, Howe K, Johnson N, Jenkinson A, Kähäri A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Megy K, Meidl P, Overduin B, Parker A, Pritchard B, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Slater G, Smedley D, Spudich G, Trevanion S, Vilella AJ, Vogel J, White S, Wood M, Birney E, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Herrero J, Hubbard TJ, Kasprzyk A, Proctor G, Smith J, Ureta-Vidal A, Searle S: Ensembl 2008.

    Nucleic Acids Res 2008, (36 Database):D707-714. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Inoue JG, Miya M, Tsukamoto K, Nishida M: Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish".

    Mol Phylogenet Evol 2003, 26:110-120. PubMed Abstract | Publisher Full Text OpenURL

  29. Inoue JG, Miya M, Venkatesh B, Nishida M: The mitochondrial genome of Indonesian coelacanth Latimeria menadoensis (Sarcopterygii: Coelacanthiformes) and divergence time estimation between the two coelacanths.

    Gene 2005, 349:227-235. PubMed Abstract | Publisher Full Text OpenURL

  30. Lavoué S, Miya M, Inoue JG, Saitoh K, Ishiguro NB, Nishida M: Molecular systematics of the gonorynchiform fishes (Teleostei) based on whole mitogenome sequences: Implications for higher-level relationships within the Otocephala.

    Mol Phylogenet Evol 2005, 37:165-177. PubMed Abstract | Publisher Full Text OpenURL

  31. Miya M, Satoh TP, Nishida M: The phylogenetic position of toadfishes (order Batrachoidiformes) in the higher ray-finned fish as inferred from partitioned Bayesian analysis of 102 whole mitochondrial genome sequences.

    Biol J Linn Soc Lond 2005, 85:289-306. Publisher Full Text OpenURL

  32. Yamanoue Y, Miya M, Inoue JG, Matsuura K, Nishida M: The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes.

    Genes Genet Syst 2006, 81:29-39. PubMed Abstract | Publisher Full Text OpenURL

  33. Azuma Y, Kumazawa Y, Miya M, K Mabuchi, Nishida M: Mitogenomic evaluation of the historical biogeography of cichlids toward reliable dating of teleostean divergences.

    BMC Evol Biol 2008, 8:215. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Hoegg S, Meyer A: Phylogenomic analyses of KCNA clusters in vertebrates: why do some clusters stay intact?

    BMC Evol Biol 2007, 7:139. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  35. Yu WP, Yew K, Rajasegaran V, Venkatesh B: Sequencing and comparative analysis of fugu protocadherin clusters reveal diversity of protocadherin genes among teleosts.

    BMC Evol Biol 2007, 7:49. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  36. 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-1816. PubMed Abstract | Publisher Full Text OpenURL

  37. Blomme T, Vandepoele K, De Bodt S, Simillion C, Maere S, Peer Y: The gain and loss of genes during 600 million years of vertebrate evolution.

    Genome Biol 2006, 7:R43. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  38. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome.

    Nucleic Acids Res 2004, (32 Database):D277-D280. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Aury JM, Jaillon O, Duret L, Noel B, Jubin C, Porcel BM, Ségurens B, Daubin V, Anthouard V, Aiach N, Arnaiz O, Billaut A, Beisson J, Blanc I, Bouhouche K, Câmara F, Duharcourt S, Guigo R, Gogendeau D, Katinka M, Keller AM, Kissmehl R, Klotz C, Koll F, Le Mouël A, Lepère G, Malinsky S, Nowacki M, Nowak JK, Plattner H, Poulain J, Ruiz F, Serrano V, Zagulski M, Dessen P, Bétermier M, Weissenbach J, Scarpelli C, Schächter 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. PubMed Abstract | Publisher Full Text OpenURL

  40. Nelson JS: Fishes of the World. 4th edition. New Jersey: John Wiley & Sons Inc; 2006. OpenURL

  41. Ensembl Genome Browser [http://www.ensembl.org/index.html] webcite

  42. Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization.

    Genetics 2000, 154:459-473. PubMed Abstract | PubMed Central Full Text OpenURL

  43. He X, Zhang : Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution.

    Genetics 2005, 169:1157-1164. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Rastogi S, Liberles DA: Subfunctionalization of duplicated genes as a transition state to neofunctionalization.

    BMC Evol Biol 2005, 5:28. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  45. Hughes T, Ekman D, Ardawatia H, Elofsson A, Liberles DA: Evaluating dosage compensation as a cause of duplicate gene retention in Paramecium tetraurelia.

    Genome Biol 2007, 8:213. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  46. Niimura Y, Nei M: Evolutionary dynamics of olfactory receptor genes in fishes and tetrapods.

    Proc Natl Acad Sci USA 2005, 102:6039-6044. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Gilad Y, Man O, Glusman G: A comparison of the human and chimpanzee olfactory receptor gene repertoires.

    Genome Res 2005, 15:224-230. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Jeon YH, Heo YS, Kim CM, Hyun YL, Lee TG, Ro S, Cho JM: Phosphodiesterase: overview of protein structures, potential therapeutic applications and recent progress in drug development.

    Cell Mol Life Sci 2005, 62:1198-1220. PubMed Abstract | Publisher Full Text OpenURL

  49. Sato Y, Hashiguchi Y, Nishida M: Evolution of multiple phosphodiesterase isoforms in stickleback involved in cAMP signal transduction pathway.

    BMC Syst Biol 2009, 3:23. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  50. Vogel G: Doubled genes may explain fish diversity.

    Science 1998, 281:1119-1121. PubMed Abstract | Publisher Full Text OpenURL

  51. Meyer A, Malaga-Trillo E: Vertebrate genomics: More fishy tales about Hox genes.

    Curr Biol 1999, 9:R210-R213. PubMed Abstract | Publisher Full Text OpenURL

  52. Meyer A, Schartl M: Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions.

    Curr Opin Cell Biol 1999, 11:699-704. PubMed Abstract | Publisher Full Text OpenURL

  53. Taylor JS, Peer Y, Meyer A: Genome duplication, divergent resolution and speciation.

    Trends Genet 2001, 17:299-301. PubMed Abstract | Publisher Full Text OpenURL

  54. Crow KD, Wagner GP: What is the role of genome duplication in the evolution of complexity and diversity?

    Mol Biol Evol 2006, 23:887-892. PubMed Abstract | Publisher Full Text OpenURL

  55. Donoghue PC, Purnell MA: Genome duplication, extinction and vertebrate evolution.

    Trends Ecol Evol 2005, 20:312-319. PubMed Abstract | Publisher Full Text OpenURL

  56. Birney E, Clamp M, Durbin R: GeneWise and Genomewise.

    Genome Res 2004, 14:988-995. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. 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-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Kumar S, Tamura K, Nei : MEGA3: integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment.

    Brief Bioinformatics 2004, 5:150-163. PubMed Abstract | Publisher Full Text OpenURL

  59. TREEFINDER [http://www.treefinder.de] webcite

  60. Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics.

    BMC Evol Biol 2004, 4:18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  61. Strimmer K, Rambaut A: Inferring confidence sets of possibly misspecified gene trees.

    Proc Roy Soc B 2002, 269:137-142. Publisher Full Text OpenURL

  62. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution.

    Bioinformatics 2005, 21:2104-2105. PubMed Abstract | Publisher Full Text OpenURL

  63. Posada D, Crandall KA: MODEL TEST: testing the model of DNA substitution.

    Bioinformatics Appl Note 1998, 14:817-818. Publisher Full Text OpenURL

  64. Mouse Genome Sequencing Consortium: Initial sequencing and comparative analysis of the mouse genome.

    Nature 2002, 420:520-562. PubMed Abstract | Publisher Full Text OpenURL

  65. Nei M, Roychoudhury KA: Probability of fixation of nonfunctional genes at duplicate loci.

    Am Nat 1973, 107:362-372. Publisher Full Text OpenURL

  66. Borenstein E, Shlomi T, Ruppin E, Sharan R: Gene loss rate: a probabilistic measure for the conservation of eukaryotic genes.

    Nucleic Acids Res 2007, 35:e7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Friedberg F, Taliaferro L: Calmodulin genes in zebrafish (revisited).

    Mol Biol Rep 2005, 32:55-60. PubMed Abstract | Publisher Full Text OpenURL

  68. Shi P, Zhang J: Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes.

    Mol Biol Evol 2006, 23:292-300. PubMed Abstract | Publisher Full Text OpenURL

  69. Ishimaru Y, Okada S, Naito H, Nagai T, Yasuoka A, Matsumoto I, Abe K: Two families of candidate taste receptors in fishes.

    Mech Dev 2005, 122:1310-1321. PubMed Abstract | Publisher Full Text OpenURL

  70. Hashiguchi Y, Furuta Y, Kawahara R, Nishida M: Diversification and adaptive evolution of putative sweet taste receptors in threespine stickleback.

    Gene 2007, 396:170-179. PubMed Abstract | Publisher Full Text OpenURL