Gene Loss and Evolutionary Rates Following Whole-Genome Duplication in Teleost Fishes Frédéric G. Brunet,* Hugues Roest Crollius, Mathilde Paris,* Jean-Marc Aury,à Patricia Gibert,§ Olivier Jaillon,à Vincent Laudet,* and Marc Robinson-Rechavik Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals. Introduction Whole-genome duplication (WGD) is expected to have a large impact on the evolution of lineages in which it has occurred. Understanding this impact necessitates unraveling a complex network of causes and consequences, including which genes are retained in duplicate, and how the duplication modifies their evolutionary constraints. The most progress on these and related issues has been made in the yeast Saccharomyces cerevisiae (Kellis et al. 2004; Davis and Petrov 2005; Scannell et al. 2006) as well as in the plant Arabidopsis thaliana (Blanc and Wolfe 2004; Maere et al. 2005). Studies in animals are less advanced, may be because the species that are the easiest to study did not experience such genome duplication (fruit flies and nematodes). Three WGD events have been proposed in ancient vertebrate history: 2 at the origin of the group and a third specific to fishes (Meyer and Van de Peer 2005). Although the 2 rounds of duplication at the origin of vertebrates remain controversial (but see Dehal and Boore 2005), fishes have provided the first clear demonstration of ancient genome duplication in vertebrate evolution (Jaillon et al. 2004; Woods et al. 2005). More recent tetraploids are known in many vertebrate lineages (Otto and Whitton 2000), but they do not give insight into the longterm consequences of duplication and cannot explain the eventual success of tetraploid lineages, in the way that more ancient duplications do. Thus, fishes provide a unique platform to study genome duplication in vertebrates. A genome duplication was originally suggested in ray-finned fishes based on the finding that zebrafish and medaka possess 7 Hox clusters (Amores et al. 1998; Wittbrodt Key words: genome duplication, gene loss, neofunctionalization, Gene Ontology, evolutionary rates, selection. E-mail: marc.robinson-rechavi@unil.ch. Mol. Biol. Evol. 23(9):1808–1816. 2006 doi:10.1093/molbev/msl049 Advance Access publication June 29, 2006 et al. 1998; Naruse et al. 2000), against 4 in mammals and 1 in most invertebrates, as well as by comparative mapping (Postlethwait et al. 2000), but it was only confirmed with the release of 2 pufferfishes’ genome sequences. First, 2 studies identified many short duplicated groups of linked genes (paralogons) in the Takifugu rubripes genome and dated duplications with a molecular clock to a window between divergence of ray-finned fishes from tetrapods and diversification of teleost fishes (Christoffels et al. 2004; Vandepoele et al. 2004). Second, all chromosomes of Tetraodon nigroviridis were assigned to paralogons, demonstrating the genome scale of the duplication, and each pair of paralogons was shown to be homologous to 1 (nonduplicated) human chromosomal region (Jaillon et al. 2004). It has been estimated that 76% to 80% of paralogs were secondarily lost after the WGD (Jaillon et al. 2004; Woods et al. 2005). When did this duplication occur? Direct dating of fish gene duplications has relied on molecular clock methods with large error bars, using 2 closely related Tetraodontiformes (Tetraodon and Takifugu). Molecular clock dating can be biased by saturation of synonymous changes at the timescales considered as well as by changes in evolutionary rates between mammals and fishes and in duplicate genes (Robinson-Rechavi and Laudet 2001; Van de Peer et al. 2001; Jaillon et al. 2004; Jordan et al. 2004; Venkatesh et al. 2005). Less sensitive to these problems, phylogenies of a few tens of gene families have shown a high frequency of gene duplications to be a common feature among sampled teleosts, but not among other fishes (RobinsonRechavi et al. 2001; Hoegg et al. 2004), confirmed by the distribution of Hox complexes (Crow et al. 2006). Finally, comparative mapping shows that paralogons are homologous between Tetraodon and zebrafish (a Cypriniforme), confirming that the WGD occurred before the divergence of these 2 lineages of teleosts (Woods et al. 2005). The aim of this study is to provide a framework for understanding the consequences of WGD in teleost Ó 2006 The Authors This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/ by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 *Laboratoire de Biologie Moléculaire de la Cellule, INRA LA 1237, CNRS UMR5161, IFR 128 BioSciences Lyon-Gerland, Ecole Normale Supérieure de Lyon, Lyon, France; Dyogen Group, CNRS UMR8541, Ecole Normale Supérieure, Paris, France; àStructure et évolution des génomes, CNRS UMR8030, Genoscope, Evry, France; §Laboratoire de Biométrie et Biologie Evolutive, CNRS UMR5558, Université Claude Bernard Lyon 1, Villeurbanne, France; and kDepartment of Ecology and Evolution, Biophore, University of Lausanne, Lausanne, Switzerland Consequences of Whole Genome Duplication in Fishes 1809 Table 1 Distribution of Tetraodon Gene Duplications by Phylogenetic Relative Dating Ancestral Fish-Specific Duplication Chromosome Distribution Strong Evidencea 14 247 261 12 213 Unrootedb Weak Supportc 11 120 378 2 20 Recent Duplication Total 10 26 36 49 626 675 NOTE.—Classification of all genes that have at least 2 paralogs in the Tetraodon nigroviridis genome, which share the same human ortholog, and for which both paralogs are mapped to T. nigroviridis chromosomes. Duplications are considered ‘‘ancestral to vertebrates’’ if they occurred in vertebrates before the actinopterygian/tetrapod divergence (fig. 1A). Duplications are considered ‘‘ancestral to fishes’’ if at least 1 order of fishes outside percomorphs shares the duplication (fig. 1B). There are 447 gene phylogenies that could not be classified because of low phylogenetic resolution or because only Tetraodon and Takifugu sequences were available among fishes, making phylogenetic dating impossible. a Phylogenies with a clear outgroup, Neighbor-Joining bootstrap support over 75%, and consistent results with different methods. b Phylogenies with good support, as in the footnote ‘‘a,’’ but no good outgroup sequences available for rooting. c Phylogenies with low bootstrap support (under 75%) in Neighbor-Joining or different results depending on methods. fishes. Duplicate genes were characterized by phylogeny, evolutionary rates, and functional classification. We used the 2 published fish genomes as well as other available fish data. Materials and Methods We combined predicted genes from genome sequences and gene sequences from databases. Genes from the freshwater pufferfish T. nigroviridis were from the published version of the genome (Jaillon et al. 2004) (8.2X coverage). Genes from the pufferfish T. rubripes (Aparicio et al. 2002) and the unpublished genome of the zebrafish Danio rerio were from Ensembl (Hubbard et al. 2002), versions 22.2c and 22.3b, respectively. These genomic data were complemented with all actinopterygian fish proteins from the Hovergen ‘‘clean prot’’ database version 44 (Duret et al. 1994). This database classifies vertebrate genes in homologous families, excluding partial sequences with less than 80% coverage of the alignment of complete proteins, and provided us directly with protein alignments to compare with the genes from genomes. In addition, human and mouse proteins were recovered from Swissprot and its complement TrEMBL (Boeckmann et al. 2003), as outgroups to actinopterygian fishes. Amino acid sequences were aligned using ClustalW (Thomson et al. 1994), followed by manual adjustments if needed. We defined groups of paralogs as 2 or 3 Tetraodon genes that have the same human best hit from Swissprot, using Blast with the settings used for ‘‘exofish’’ gene detection (Roest Crollius et al. 2000; Jaillon et al. 2004); ‘‘best hit’’ was defined according to Blast score (Altschul et al. 1997). To reinforce phylogenetic with mapping analysis, we used only groups with at least 2 genes mapped to Tetraodon chromosomes. This defines a set of 704 pairs with both paralogs assigned to chromosomes, plus 419 triplets with at least 2 paralogs assigned to chromosomes, making a total of 1,123 gene families to investigate. For each corresponding amino acid alignment, a preliminary phylogenetic tree was estimated by Neighbor-Joining (Saitou and Nei 1987) with rate heterogeneity between sites corrected for by a gamma law of parameter alpha set to 0.8 and 500 bootstrap replicates. Phylogenies were then refined with regard to species sampling and partial sequences in Phylo_win (Galtier et al. 1996) and interpreted manually. All phylogenies are available in the Supplementary Material online. The rate of gene loss after genome duplication was computed using the subset of genes that are mapped to Tetraodon chromosomes and have a clear orthology relation between human and Tetraodon. Orthology was defined either from the phylogeny or by reciprocal best Blast hit, for cases of one-to-one homologs. There are 2,371 such human genes, of which 364 have 2 Tetraodon co-orthologs, which were shown in our analysis to be duplicated in the fish-specific genome duplication. This gives a rate of gene retention after the duplication of 15% (364/2,371). The figure of 378 in table 1, instead of 364, includes fish genes for which human orthology is not clear. For evolutionary rate computations, amino acid alignments were extended to the nucleotide coding sequences. For each pair, the number of synonymous substitutions per synonymous site, Ks, and the number of nonsynonymous substitutions per nonsynonymous site, Ka (Li 1993; Pamilo and Bianchi 1993), were calculated using the ‘‘diverge’’ option of the GCG software. We removed pairs with null Ka or Ks values and pairs for which Ka or Ks were not computable because sites were too few. We also manually removed predicted proteins containing clear frameshifts. Thus, evolutionary rates were computed for 254 WGD and 1,666 singleton genes (categories defined in fig. 2). For evolutionary rate computations, we updated the predicted Takifugu genes to version 31.2 from Ensembl, and we added human and mouse genes from Ensembl version 36. Human orthologs of mapped duplicated Tetraodon genes, as determined by phylogenetic analysis, were compared for over- or underrepresentation of Gene Ontology (GO) terms (Ashburner et al. 2000), using GOToolBox (Martin et al. 2004), with the hypergeometric test, and the Benjamini and Hochberg (1995) correction for multiple testing. The reference data set was composed of all human orthologs of mapped Tetraodon genes. Results Phylogenetic and Chromosomal Classification of Duplicated Genes We identified 1,123 pairs or triplets of Tetraodon paralogs with the same human best hit and for which at Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 Same Different All Ancestral Vertebrate Duplication 1810 Brunet et al. consistent with the findings of Woods et al. (2005), based on linkage mapping. Genome duplication in fishes was followed by massive gene loss because most genes in the Tetraodon or Takifugu genomes are not found in characteristic gene pairs. Out of 2,371 human genes with a clear orthology relation to mapped Tetraodon genes, 364 have a phylogenetic profile typical of the WGD in Tetraodon (see Materials and Methods). This gives an estimate of 15% of gene pairs retained after the WGD, that is, loss of one of the paralogs in 85% of pairs. This is slightly higher than the estimates of 76% to 80% reached by more approximate methods (Jaillon et al. 2004; Woods et al. 2005). It is possible that 85% is still an underestimate because we neglect the fast-evolving genes for which a human ortholog is not detected and which may be more likely to lose a paralog after duplication. The exact figure of gene loss remains to be determined, ideally by comparison to a closer nonduplicated genome, such as that of nonteleost fish. To improve our understanding of the process of gene loss and retention, we compared the characteristics of paralogs from the WGD to those of nonduplicated genes. Evolutionary Rates of Duplicate Genes We first compared the substitution rates between Tetraodon and Takifugu orthologs. As they all diverged at the time of the Tetraodon/Takifugu speciation, we can directly convert differences between evolutionary distances (Ka and Ks) into differences between evolutionary rates. We see 2 clear trends (fig. 2A): WGD paralogs differ from singleton genes and there is asymmetry between the 2 paralogs from WGD. The most significant differences are in selective pressure, as measured by the Ka/Ks ratio. But these differences can be due either to changes in selective patterns after duplication, or to bias in the genes kept after duplication, or to a combination of both. To distinguish these possible causes, we tested rates that are phylogenetically independent of the duplication. Comparing the evolutionary rates between mouse and human orthologs of the fish genes, we find results similar to the slower WGD paralogs in pufferfishes, with notably lower Ka in the orthologs of pufferfish genes that were kept in double (fig. 2A and B). Mammalian orthologs of fish WGD pairs evolve 23% slower than orthologs of singletons. These differences in evolutionary rates between orthologs of singletons and those of duplicates cannot be due to the fish WGD. Instead, they show that the fish genes that were already under the strongest selective pressure were retained as duplicates after the genome duplication. Another feature to take into account is that teleost fish genes seem to evolve faster, under weaker constraint, than mammalian genes (Robinson-Rechavi and Laudet 2001; Jaillon et al. 2004; Steinke et al. 2006). Indeed, in our data, the Ka/Ks ratio is systematically higher in pufferfishes than in mammalian orthologs (fig. 2A). The difference is highly significant for singleton genes and fastevolving paralogs (paired t-test: P , 2.2 3 10ÿ16) but not for slow-evolving paralogs. The latter are thus under relatively high selective pressure, typical of mammalian rather than pufferfish genes. Of note, the divergence time Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 least 2 genes were mapped to Tetraodon chromosomes. We conducted a phylogenetic analysis of each of the corresponding gene families, of which 675 could be classified with respect to gene duplication in fishes. Phylogenetic reconstruction was used to date duplication events, relatively to speciation events (fig. 1), as opposed to molecular clock dating. Phylogenies of duplicate genes can be biased by rate differences when few species and simple methods are used (Van de Peer et al. 2002; Fares et al. 2006). Although we can never guarantee the accuracy of all reconstructions, care was taken to maximize taxonomic sampling, to take into account the rate heterogeneity among sites, and to check automatic alignment results. Altogether, we analyzed sequences from 112 species of ray-finned fishes, sampling 13 lineages (Percomorpha 1 12 actinopterygian orders). The most notable contribution was from zebrafish, which is included in 90% of the phylogenies; 79% of phylogenies include only zebrafish and the 2 Tetraodontiformes. We also checked for long-branch attraction a la Fares et al. (2006) by comparing evolutionary rates: there is no significant difference between the ‘‘ancestral vertebrate’’ paralogs and the ‘‘ancestral fish-specific’’ paralogs (not shown). Among paralogs that diverged in the teleost fishes (table 1), most duplicated before the divergence of main teleost fish lineages (e.g., fig. 1B). We take these pairs to be the best markers of the teleost fish–specific WGD in the following analyses. A previous classification of Tetraodon paralogs, based on a simple synonymous molecular clock (Jaillon et al. 2004), provided evidence for paralogons, that is, large duplicated regions that often comprise whole chromosomes. We repeated this analysis, but using exclusively the paralogs whose period of duplication has been phylogenetically verified to be fish specific (Supplementary Figure I, Supplementary Material online). This confirms the coverage of the Tetraodon genome by paralogons in a robust manner. Twothirds of all ancient fish-specific paralogs are part of a duplicated block of at least 5 genes. Conversely, recent duplicates are the only class of paralogs that are frequently neighbors on the same chromosome (table 1). Consistency between chromosomal order and phylogeny supports the validity of our phylogenetic analysis in identifying the paralogs from WGD (Fares et al. 2006). This also provides a verification of the ‘‘ancestral vertebrate’’ paralogs (table 1): 80% are not included in fish-specific paralogons, consistent with a separate origin. Interestingly, these ‘‘ancestral vertebrate’’ paralogs define 6 duplicated blocks (of at least 5 genes), which are distinct from the fish-specific paralogons; these may be remnants of ancestral vertebrate genome duplication (Dehal and Boore 2005). We verified orthology relations between zebrafish and Tetraodon chromosomes (Supplementary Figure II, Supplementary Material online): most zebrafish duplicated linkage groups are indeed orthologous to Tetraodon paralogons. The duplicated blocks from 1 species for which we did not identify an homologous duplicated block in the other species are those defined by the smallest number of gene pairs. Thus, the combined phylogenetic and mapping analysis confirms that the paralogons covering Tetraodon and zebrafish chromosomes find their origin in a WGD that took place before the divergence of these species. This is Consequences of Whole Genome Duplication in Fishes 1811 A TR_155330 89 100 93 TN_05197_1 Q7T2N7 Fugu B Tetraodon 100 53 TR_159362 Fugu TN_19087_7 Tetraodon goldfish DR_03264_u mammals AAQ91267 FGFI_HUMAN zebrafish 100 Q7SX66 zebrafish 100 100 TR_148211 100 mouse Fugu TN_31158_1 Tetraodon Fugu REQU_MOUSE TN_24466_20 Tetraodon BAC36859 TR_142024 100 mammals 100 REQU_HUMAN 100 Q80ZL6 99 100 O57341 mouse zebrafish Q8CAD8 mouse O42278 0.05 C zebrafish 0.02 AAH52348 DR_15947_u zebrafish 100 TR_150018 83 Fugu TN_32094_10 Tetraodon 99 TN_17687_5 COG4_HUMAN mammals 0.05 COG4_MOUSE FIG. 1.—Examples of relative dating of gene duplications by phylogeny. Representative examples of the main phylogenetic scenarios observed. Neighbor-Joining phylogenetic trees, with substitutions corrected for heterogeneity between sites by a gamma law (a 5 0.8); Maximum Likelihood gives identical topologies (not shown). Branch length is proportional to the number of substitutions per site (see measure bar for each tree); numbers at nodes are support in percentage of 500 bootstrap replicates. The stars represent the duplication that led to 2 Tetraodon paralogs in each case. (A) Fibroblast growth factor-18 precursor, duplicated before the tetrapod–actinopterygian split (ancestral vertebrate duplication), with secondary loss of 1 copy in mammals, leading to 2 copies in fishes but only 1 in human; clade BAC36859–O42278 is a more distant paralog used as an outgroup. (B) Zinc-finger protein ubi-d4 (Requiem), duplicated anciently in fishes; clade Q8CAD8–AAH52348 is a more distant paralog used as an outgroup. (C) Conserved oligomeric Golgi complex component 4, duplicated specifically in Tetraodon, after the divergence with other fishes. between Tetraodon and Takifugu was recently reassessed at 85 MYA (Yamanoue et al. 2006), similar to the divergence between human and mouse. This affects neither relative rate tests (Robinson-Rechavi and Laudet 2001; Steinke et al. 2006) nor comparisons of Ka/Ks ratios (this study), which are both independent of divergence time. Once we have established contributions to evolutionary rate, which are not a direct result of duplication (biased retention and species-specific differences), can we also measure the direct effect of duplication? Lower purifying selection is expected on redundant genes after duplication. In principle, both paralogs could experience released Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 100 zebrafish 77 FGFI_MOUSE 1812 Brunet et al. selection and higher rates, or the relaxation could be restricted to only one paralog, whereas the other would conserve ancestral constraints and evolutionary rates. We first tested whether the observed difference between WGD paralogs is meaningful. Indeed 2 homologous gene pairs will never have exactly the same rate because of stochastic as well as mechanistic factors, such as variations in mutation rates along the genome. For this, we contrast variation in Ka to variation in Ks, used as a proxy for neutral variation in evolutionary rate. The relative variation in Ka is significantly larger than the variation in Ks (figs. 2A and 3), supporting a selective difference between the paralogs. A striking 36% of paralogous pairs have higher Ka variation than predicted by the neutral distribution (DKa higher than the 97.5th percentile of DKs distribution). And this is a conservative estimate, considering that for 81% of paralogs, DKa is higher than DKs, indicating potential selective differences. The results are similar if we exclude pairs with mean Ka , 0.05, to account for possible stochastic effects on small numbers of substitutions (not shown). Thus, the rate differences between the 2 paralogs can be largely interpreted in terms of selective differences. If we compare the rates of these paralogs to those of mammalian orthologs on the one hand and to singleton pufferfish genes on the other, it seems that the slower copy has retained the strong constraint of its ancestor (Ka2/Ks2 Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 FIG. 2.—Evolutionary rates of singleton and duplicate genes in Tetraodon nigroviridis. (A) Genes are classified according to duplication history. Tn 5 T. nigroviridis; Tr 5 Takifugu rubripes; Hs 5 Homo sapiens; Mm 5 Mus musculus. Ka 5 mean number of nonsynonymous substitutions per nonsynonymous site, measuring protein evolutionary rate; Ks 5 mean number of synonymous substitutions per synonymous site, measuring neutral evolutionary rate; Ka/Ks 5 mean ratio of nonsynonymous to synonymous substitutions, measuring selective pressure. All means are 6 standard deviation. The arrows indicate which sequences were compared. In red and brown, the comparison between pufferfish orthologs; for WGD paralogs, we report separately the values for ortholog pairs with highest Ka (in brown) and with lowest Ka (in red); stars indicate a significant difference to the pufferfish singleton values (unpaired t-test: * P , 0.05; ** P , 10ÿ4; *** P , 10ÿ6). In blue, the comparison between mouse and human orthologs; stars indicate a significant difference to the mammalian singleton values. (B) Frequency distributions of Ka values between human and mouse orthologous genes, according to their duplication history in fishes. Consequences of Whole Genome Duplication in Fishes 1813 not significantly different from mammalian ortholog Ka/ Ks), whereas the other copy has experienced a large relaxation of selective constraint (Ka1/Ks1 highest of any subset of data, P , 10ÿ6 in all comparisons). We have only identified 36 unambiguous recently duplicated genes in Tetraodon, and Ka and Ks can be computed for both fishes and mammals for only 23 of them. This low number limits statistical power, but we observe that on average they have intermediary evolutionary rates and selective pressures, between those of paralogs from the WGD and those of singletons (i.e., Ka of mammalian orthologs 5 0.052). Interestingly, there appear to be selective differences between these recent paralogs, as between WGD paralogs, although sample size limits significance (DKa 5 39%; DKs 5 22%; P 5 0.056). Functional Classification of Duplicate Genes We examined the possibility that the genome duplication may have retained duplicate genes in specific functional classes, by using the GO classification (Ashburner et al. 2000) of human orthologs of Tetraodon genes. Orthologs of WGD paralogs are significantly enriched in terms that relate to development, signaling, behavior, and regulation (table 2). This enrichment is confirmed if we consider the less detailed level 2 GO terms (not shown), with notably an excess of genes associated to the biological process ‘‘development’’ (P 5 0.0065). Conversely, they are depleted in terms describing fundamental processes of the cell, such as metabolism or catalytic activity. For the most part, the same terms are enriched or depleted in the recently duplicated genes, but low sample size again limits significance (not shown). Discussion We interpret our results in light of the now classical ‘‘duplication–degeneration–complementation’’(DDC)model for the functional consequences of gene duplication (Force et al. 1999). The DDC model proposes 3 possible fates for the new paralogs: nonfunctionalization, in which 1 copy is lost and the situation reverts to its preduplication state; neofunctionalization, in which 1 copy acquires a new function while the other keeps the ancestral function; and subfunctionalization, in which each copy loses part of the ancestral function, both copies being required then to keep the full function active. It is clear from genomic studies that nonfunctionalization is the main fate of duplicate genes: gene loss appears to have been 88% in about 80 Myr since genome duplication in yeasts (Kellis et al. 2004), 70% in 86 Myr in Arabidopsis (Bowers et al. 2003), and 79% in about 61–67 Myr in cereals (Paterson et al. 2004). The figure of 85% we obtain for gene loss in Tetraodon is similar to the observations in other lineages, Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 FIG. 3.—Variation in evolutionary rates after WGD in Tetraodontiformes. Distributions of DKs (red) and DKa (blue) values between WGD paralogs. Ka1, Ka2, Ks1, and Ks2 as defined in figure 2A. The arrow shows the range of variation in DKs for 95% of gene pairs, defining a neutral expectation of rate variation between paralogs. 1814 Brunet et al. Table 2 Functional Classification of Tetraodon Duplicated Genes Ontology Biological process Cellular component Observed Gene Number Enrichment/ Depletion P Value Cell communication Organismal physiological process Death Locomotory behavior System development Mesoderm development Tissue development Regulation of cellular process Response to stimulus Regulation of physiological process Metabolism Cellular physiological process Pattern binding Neurotransmitter binding Structural constituent of eye lens Ion binding Carbohydrate binding Channel or pore class transporter activity Catalytic activity Nucleic acid binding Oxidoreductase activity Ligase activity Hydrolase activity Nucleotide binding Connexon complex Extracellular matrix (sensu Metazoa) Membrane Membrane-bound organelle Intracellular organelle Intracellular 43 24 10 6 9 3 5 29 19 28 54 77 6 4 2 30 6 4 41 13 4 1 13 13 3 9 35 39 45 58 31.79 32.04 32.86 33.21 32.14 35.15 32.67 31.24 31.35 31.22 O1.31 O1.06 33.16 33.74 38.42 31.33 32.53 32.59 O1.28 O1.54 O2.23 O4.87 O1.45 O1.43 38.22 32.54 31.28 O1.30 O1.23 O1.13 3.7 3 10ÿ06 1.2 3 10ÿ04 9.5 3 10ÿ04 0.0054 0.012 0.012 0.023 0.037 0.039 0.042 4.4 3 10ÿ05 0.050 0.0059 0.013 0.014 0.017 0.018 0.045 0.0045 0.018 0.027 0.029 0.031 0.034 0.0018 0.0038 0.017 0.0029 0.0051 0.015 NOTE.—Statistics on human orthologs of Tetraodon genes, which were classified according to duplication history as in figure 2. The observed gene number is the number of human genes annotated for each term. Only GO terms of level 3 or higher with a significant enrichment or depletion in paralogs from WGD are reported (P  0.05 after correction for multiple tests). Enrichment (3 n) and depletion (O n) are reported as the ratio between expected and observed frequency of terms, based on expectations from human genes that were defined as orthologs of Tetraodon genes in the combined phylogenetic and mapping analysis. despite the greater age of the event; even less gene loss has been reported based on similarity hits without phylogenetic analysis (Jaillon et al. 2004; Woods et al. 2005). This is probably because most gene loss occurs rapidly after duplication (Lynch and Conery 2000; Jaillon et al. 2004; Scannell et al. 2006), so that subsequent evolution does not change the figure significantly. Slowly evolving genes are more likely to be found duplicated in the yeast S. cerevisiae and the nematode worm Caenorhabditis elegans (Davis and Petrov 2004). The bias is similar in both yeast and worm and is maintained over evolutionary time, indicating that gene retention was also biased after the WGD in yeast. Jordan et al. (2004) also showed that genes that have paralogs evolve slower than singletons, in several eukaryotic and prokaryotic genomes, but did not use a phylogenetic framework with independent estimates of duplication and evolutionary rate. Using phylogenetically independent estimates, we show that such a bias exists in duplicated genes in a vertebrate genome (fig. 2): nonduplicated mammalian orthologs of gene pairs retained from a WGD in fishes evolve 23% slower in mammals than orthologs of singleton genes. This is a comparable figure to observations for nematode and yeast genes (Davis and Petrov [2004] computed differences relative to the slower rate; this gives 30% in fishes, compared to 25% in nematode and 50% in yeast). Because these gene pairs result from 1 WGD event, where all genes were duplicated, the bias must come from a differential retention of the duplicates, rather than a biased generation. Unraveling the contributions of biased retention, different rates in fishes and mammals, and duplication itself, we find evidence for asymmetric acceleration of protein evolution after WGD. A recent study found evidence for different rates of evolution between fish-specific paralogs in 24 genes as well as for increased rates of radical amino acid substitution (Steinke et al. 2006), although it is not clear to what extent this constitutes evidence for positive selection. In yeast, differences in evolutionary rate between WGD paralogs have been found to be pervasive and associated to some evidence for positive selection on amino acid sites (Fares et al. 2006). In that study, it was estimated that 19% of duplicate genes have fixed amino acids by positive selection. We find considerably higher differences in Ka than in Ks between paralogs, which also indicates a change in selective pressure, for at least 36% of gene pairs. The difference appears due to acceleration of nonsynonymous rate in 1 paralog of the pair, whereas the other keeps the ancestral constraint. Therefore, as in yeast, fly, and nematode (Conant and Wagner 2003; Kellis et al. 2004), results are suggestive of neofunctionalization of duplicated genes in vertebrates, insofar as this can be indicated by evolutionary rates. This would Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 Molecular function GO Term Consequences of Whole Genome Duplication in Fishes 1815 Supplementary Material Supplementary Table I and Figures I and II are available at Molecular Biology and Evolution online (http:// www.mbe.oxfordjournals.org/). Acknowledgments We are grateful to the Sanger Institute for making zebrafish genome data publicly available. Laurent Duret and Simon Penel helped with fish data from Hovergen. F.G.B. and V.L. acknowledge support from Région Rhone-Alpes, Ministére de l’éducation nationale, de la recherche et de la technologie, European Union, and Institut national de la recherche agronomique. Funding to pay the Open Access publication charges for this article was provided by Etat de Vaud. Literature Cited Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389–402. Amores A, Force A, Yan YL, et al. (13 co-authors). 1998. Zebrafish hox clusters and vertebrate genome evolution. Science 282:1711–4. Aparicio S, Chapman J, Stupka E, et al. (41 co-authors). 2002. Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science 297:1301–10. Ashburner M, Ball CA, Blake JA, et al. (20 co-authors). 2000. Gene ontology: tool for the unification of biology. Nat Genet 25:25–9. Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc 57:289–300. Blanc G, Wolfe KH. 2004. Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell 16:1679–91. Boeckmann B, Bairoch A, Apweiler R, et al. (12 co-authors). 2003. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res 31:365–70. Bowers JE, Chapman BA, Rong J, Paterson AH. 2003. Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422:433–8. Christoffels A, Koh EG, Chia JM, Brenner S, Aparicio S, Venkatesh B. 2004. Fugu genome analysis provides evidence for a wholegenome duplication early during the evolution of ray-finned fishes. Mol Biol Evol 21:1146–51. Conant GC, Wagner A. 2003. Asymmetric sequence divergence of duplicate genes. Genome Res 13:2052–8. Crow KD, Stadler PF, Lynch VJ, Amemiya C, Wagner GP. 2006. The ‘‘fish-specific’’ Hox cluster duplication is coincident with the origin of teleosts. Mol Biol Evol 23:121–36. Davis JC, Petrov DA. 2004. Preferential duplication of conserved proteins in eukaryotic genomes. PLoS Biol 2:e55. Davis JC, Petrov DA. 2005. Do disparate mechanisms of duplication add similar genes to the genome? Trends Genet 21:548–51. Dehal P, Boore JL. 2005. Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol 3:e314. Duret L, Mouchiroud D, Gouy M. 1994. HOVERGEN: a database of homologous vertebrate genes. Nucleic Acids Res 22:2360–5. Fares MA, Byrne KP, Wolfe KH. 2006. Rate asymmetry after genome duplication causes substantial long-branch attraction artifacts in the phylogeny of Saccharomyces species. Mol Biol Evol 23:245–53. Force A, Lynch M, Pickett FB, Amores A, Yan Y-l, Postlethwait J. 1999. Preservation of duplicate genes by complementary, degenerative mutations. Genetics 151:1531–1545. Galtier N, Gouy M, Gautier C. 1996. SEAVIEW and PHYLO_ WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci 12:543–8. He X, Zhang J. 2005. Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics 169:1157–64. Hoegg S, Brinkmann H, Taylor JS, Meyer A. 2004. Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish. J Mol Evol 59: 190–203. Hubbard T, Barker D, Birney E, et al. (35 co-authors). 2002. The Ensembl genome database project. Nucleic Acids Res 30: 38–41. Jaillon O, Aury J-M, Brunet F, et al. (61 co-authors). 2004. Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature 431: 946–57. Jordan IK, Wolf YI, Koonin EV. 2004. Duplicated genes evolve slower than singletons despite the initial rate increase. BMC Evol Biol 4:22. Kellis M, Birren BW, Lander ES. 2004. Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature 428:617–24. Li WH. 1993. Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J Mol Evol 36:96–9. Lynch M, Conery JS. 2000. The evolutionary fate and consequences of duplicate genes. Science 290:1151–5. Maere S, De Bodt S, Raes J, Casneuf T, Van Montagu M, Kuiper M, Van de Peer Y. 2005. Modeling gene and genome Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 be in agreement with the original evolutionary theory of Ohno (1970) as well as with a more recent model of subfunctionalization followed by neofunctionalization (He and Zhang 2005). Retention of duplicates is also biased with regard to function (table 2). The excess of development and signaling functions supports the much discussed idea that genome duplication may have played an important role in establishing developmental diversity of fishes (Vogel 1998; Wittbrodt et al. 1998; Meyer and Schartl 1999; Venkatesh 2003; Postlethwait et al. 2004; Meyer and Van de Peer 2005; Volff 2005). The duplication may also have been important in establishing behavioral diversity. Interestingly, by a quite different approach, Steinke et al. (2006) found an excess of transcription factors among genes with accelerated evolution or specific duplication in fishes. The enrichment in communication and developmental genes is consistent with observations in insects and yeasts (Jordan et al. 2004) and Arabidopsis (Maere et al. 2005), although another study found an enrichment of metabolic genes in Arabidopsis (Blanc and Wolfe 2004), and contrasting results were found in mammals (Nembaware et al. 2002; Jordan et al. 2004). In conclusion, the consequences of WGD in fishes have been mediated both by biased gene loss, resulting in enrichment of development, signaling, or behavioral genes compared with mammals, and by changes in selective pressure, asymmetric between duplicates. The functional consequences of these changes remain to be explored. 1816 Brunet et al. Scannell DR, Byrne KP, Gordon JL, Wong S, Wolfe KH. 2006. Multiple rounds of speciation associated with reciprocal gene loss in polyploid yeasts. Nature 440:341–5. Steinke D, Salzburger W, Braasch I, Meyer A. 2006. Many genes in fish have species-specific asymmetric rates of molecular evolution. BMC Genomics 7:20. Thomson JD, Higgins DG, Gibson TJ. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673–80. Van de Peer Y, Frickey T, Taylor JS, Meyer A. 2002. Dealing with saturation at the amino acid level: a case study based on anciently duplicated zebrafish genes. Gene 295:205–11. Van de Peer Y, Taylor JS, Braasch I, Meyer A. 2001. The ghosts of selection past: rates of evolution and functional divergence in anciently duplicated genes. J Mol Evol 53:434–44. Vandepoele K, De Vos W, Taylor JS, Meyer A, Van de Peer Y. 2004. 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 101:1638–43. Venkatesh B. 2003. Evolution and diversity of fish genomes. Curr Opin Genet Dev 13:588–92. Venkatesh B, Tay A, Dandona N, Patil JG, Brenner S. 2005. A compact cartilaginous fish model genome. Curr Biol 15:R82–3. Vogel G. 1998. Doubled genes may explain fish diversity. Science 281:1119–21. Volff JN. 2005. Genome evolution and biodiversity in teleost fish. Heredity 94:280–94. Wittbrodt J, Meyer A, Schartl M. 1998. More genes in fish? Bioessays 20:511–5. Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS. 2005. The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res 15:1307–14. Yamanoue Y, Miya M, Inoue JG, Matsuura K, Nishida M. 2006. The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes. Genes Genet Syst 81:29–39. Kenneth Wolfe, Associate Editor Accepted June 23, 2006 Downloaded from https://academic.oup.com/mbe/article/23/9/1808/1014301 by EVES-Escola Valenciana dÉstudis de la Salut user on 08 June 2022 duplications in eukaryotes. Proc Natl Acad Sci USA 102: 5454–9. Martin D, Brun C, Remy E, Mouren P, Thieffry D, Jacq B. 2004. GOToolBox: functional analysis of gene datasets based on Gene Ontology. Genome Biol 5:R101. Meyer A, Schartl M. 1999. 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 1:699–704. Meyer A, Van de Peer Y. 2005. From 2R to 3R: evidence for a fishspecific genome duplication (FSGD). Bioessays 27:937–45. Naruse K, Fukamachi S, Mitani H, et al. (20 co-authors). 2000. A detailed linkage map of medaka, Oryzias latipes: comparative genomics and genome evolution. Genetics 154:1773–84. Nembaware V, Crum K, Kelso J, Seoighe C. 2002. Impact of the presence of paralogs on sequence divergence in a set of mousehuman orthologs. Genome Res 12:1370–6. Ohno S. 1970. Evolution by gene duplication. Heidelberg: Springer-Verlag. Otto SP, Whitton J. 2000. Polyploid incidence and evolution. Annu Rev Genet 34:401–37. Pamilo P, Bianchi NO. 1993. Evolution of the Zfx and Zfy genes: rates and interdependence between the genes. Mol Biol Evol 10:271–81. Paterson AH, Bowers JE, Chapman BA. 2004. Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics. Proc Natl Acad Sci USA 101:9903–8. Postlethwait J, Amores A, Cresko W, Singer A, Yan Y-L. 2004. Subfunction partitioning, the teleost radiation and the annotation of the human genome. Trends Genet 20:481–90. Postlethwait JH, Woods IG, Ngo-Hazelett P, Yan YL, Kelly PD, Chu F, Huang H, Hill-Force A, Talbot WS. 2000. Zebrafish comparative genomics and the origins of vertebrate chromosomes. Genome Res 10:1890–902. Robinson-Rechavi M, Laudet V. 2001. Evolutionary rates of duplicate genes in fish and mammals. Mol Biol Evol 18:681–3. Robinson-Rechavi M, Marchand O, Escriva H, Bardet P-L, Zelus D, Hughes S, Laudet V. 2001. Euteleost fish genomes are characterized by expansion of gene families. Genome Res 11:781–8. Roest Crollius H, Jaillon O, Bernot A, et al. (12 co-authors). 2000. Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence. Nat Genet 25:235–8. Saitou N, Nei M. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–25.