Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3′-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.     

Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes

With about 24,000 extant species, teleosts are the largest group of vertebrates. They constitute more than 99% of the ray-finned fishes (Actinopterygii) that diverged from the lobe-finned fish lineage (Sarcopterygii) about 450 MYA. Although the role of genome duplication in the evolution of vertebrates is now established, its role in structuring the teleost genomes has been controversial. At least two hypotheses have been proposed: a whole-genome duplication in an ancient ray-finned fish and independent gene duplications in different lineages. These hypotheses are, however, based on small data sets and lack adequate statistical and phylogenetic support. In this study, we have made a systematic comparison of the draft genome sequences of Fugu and humans to identify paralogous chromosomal regions ("paralogons") in the Fugu that arose in the ray-finned fish lineage ("fish-specific"). We identified duplicate genes in the Fugu by phylogenetic analyses of the Fugu, human, and invertebrate sequences. Our analyses provide evidence for 425 fish-specific duplicate genes in the Fugu and show that at least 6.6% of the genome is represented by fish-specific paralogons. We estimated the ages of Fugu duplicate genes and paralogons using the molecular clock. Remarkably, the ages of duplicate genes and paralogons are clustered, with a peak around 350 MYA. These data strongly suggest a whole-genome duplication event early during the evolution of ray-finned fishes, probably before the origin of teleosts.     

Reciprocal gene loss between Tetraodon and zebrafish after whole genome duplication in their ancestor

The whole genome duplication that occurred in ray-finned fish coincided with the radiation of teleost species; consequently, these two phenomena have often been linked. Using the Tetraodon and zebrafish complete genome sequences, we tested a molecular hypothesis that can relate whole genome duplication to speciation in teleosts. We estimate that thousands of genes that remained duplicated when Tetraodon and zebrafish diverged underwent reciprocal loss subsequently in these two species, probably contributing to reproductive isolation between them.     

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

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

A RAD-Tag Genetic Map for the Platyfish (Xiphophorus maculatus) Reveals Mechanisms of Karyotype Evolution Among Teleost Fish

Mammalian genomes can vary substantially in haploid chromosome number even within a small taxon (e.g., 3–40 among deer alone); in contrast, teleost fish genomes are stable (24–25 in 58% of teleosts), but we do not yet understand the mechanisms that account for differences in karyotype stability. Among perciform teleosts, platyfish (Xiphophorus maculatus) and medaka (Oryzias latipes) both have 24 chromosome pairs, but threespine stickleback (Gasterosteus aculeatus) and green pufferfish (Tetraodon nigroviridis) have just 21 pairs. To understand the evolution of teleost genomes, we made a platyfish meiotic map containing 16,114 mapped markers scored on 267 backcross fish. We tiled genomic contigs along the map to create chromosome-length genome assemblies. Genome-wide comparisons of conserved synteny showed that platyfish and medaka karyotypes remained remarkably similar with few interchromosomal translocations but with numerous intrachromosomal rearrangements (transpositions and inversions) since their lineages diverged ∼120 million years ago. Comparative genomics with platyfish shows how reduced chromosome numbers in stickleback and green pufferfish arose by fusion of pairs of ancestral chromosomes after their lineages diverged from platyfish ∼195 million years ago. Zebrafish and human genomes provide outgroups to root observed changes. These studies identify likely genome assembly errors, characterize chromosome fusion events, distinguish lineage-independent chromosome fusions, show that the teleost genome duplication does not appear to have accelerated the rate of translocations, and reveal the stability of syntenies and gene orders in teleost chromosomes over hundreds of millions of years.     

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

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

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

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

Molecular Cloning,Functional Characterization,and Evolutionary Analysis of Vitamin D Receptors Isolated from Basal Vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC₅₀ values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D₃ acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D₃ throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.     

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

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

From 2R to 3R: evidence for a fish-specific genome duplication (FSGD)

An important mechanism for the evolution of phenotypic complexity, diversity and innovation, and the origin of novel gene functions is the duplication of genes and entire genomes. Recent phylogenomic studies suggest that, during the evolution of vertebrates, the entire genome was duplicated in two rounds (2R) of duplication. Later, approximately 350 mya, in the stem lineage of ray-finned (actinopterygian) fishes, but not in that of the land vertebrates, a third genome duplication occurred-the fish-specific genome duplication (FSGD or 3R), leading, at least initially, to up to eight copies of the ancestral deuterostome genome. Therefore, the sarcopterygian (lobe-finned fishes and tetrapods) genome possessed originally only half as many genes compared to the derived fishes, just like the most-basal and species-poor lineages of extant fishes that diverged from the fish stem lineage before the 3R duplication. Most duplicated genes were secondarily lost, yet some evolved new functions. The genomic complexity of the teleosts might be the reason for their evolutionary success and astounding biological diversity.     

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

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.     

The fates of zebrafish Hox gene duplicates

In his 1970 book, Susumu Ohno stressed the importance of gene duplication in the evolution of the vertebrate genome and body plan. He elaborated the idea that duplication events provide novel genetic material on which evolution may act. Data are accumulating to show that extensive duplication events, perhaps incorporating the duplication of entire genomes, occurred in the lineage leading to teleost fishes. These duplications may have been pivotal in the explosive radiation of this highly successful vertebrate group. Thus, the teleosts provide us with an ideal opportunity to investigate the fates and functions of duplicated genes. A convenient system for these studies is the zebrafish, Danio rerio, which has become a popular genetic and embryological model.     

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

Background  

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

Birth and death of neuropeptide Y receptor genes in relation to the teleost fish tetraploidization

Extensive evidence exists for a genome duplication in the fish lineage leading to the species-rich clade of the teleosts, comprising > 99% of the known actinopterygian (ray-finned) fish species. Our previous studies of the neuropeptide Y receptor (NPYR) gene family suggested an ancestral gnathostome repertoire of 7 genes in 3 subfamilies. However, studies in the zebrafish have earlier identified only 5 NPYR genes, despite the expected increase in gene number due to the teleost tetraploidization. Notably, receptors Y(1), Y(5) and Y(6) were missing in the zebrafish genome database and only Y(8) had been duplicated. We report here an investigation of the evolutionary history of the Y(1) subfamily (Y(1), Y(4), Y(6) and Y(8)) and the Y(5) receptor. Seven basal actinopterygian species and a shark were investigated and a total of 22 gene fragments were cloned and analyzed. Our results show that subtypes Y(1), Y(5) and Y(6) still exist in species representing basal actinopterygian lineages (bichir, sturgeon, gar and bowfin) as well as in some basal teleost lineages. Surprisingly we identified a zebrafish Y(1) receptor, the first Y(1) receptor found in euteleosts. Thus, these findings confirm the ancestral gnathostome repertoire of 7 NPYR genes and show that many of these receptors are present in basal actinopterygians as well as some basal teleosts. NPYR losses seem to have occurred relatively recently in euteleosts because Y(1), Y(5) and Y(6) are absent in the genome databases of two pufferfishes as well as medaka and stickleback and Y(5) and Y(6) are absent in the zebrafish database. A duplicate of Y(8) seems to be the only remaining receptor gene resulting from the teleost tetraploidization. The unexpected absence of the two appetite-stimulating receptors Y(1) and Y(5) in some euteleosts, along with our discovery of duplicates of the peptide ligands NPY and PYY, has implications for the role of the NPY system in euteleost feeding behavior.     

鲤sPLA2-III基因结构、系统发育和表达特征

磷脂酶A₂在磷脂代谢、炎症反应、细胞增殖和动脉粥样硬化等生理病理过程中发挥作用。为了探究sPLA₂-III在鲤(Cyprinus carpio Linnaeus)中的生物学功能, 实验使用BLAST同源搜索和线性分析在鲤 (Cyprinus carpio)基因组中挖掘得到5个sPLA₂-III基因, 分别位于5个连锁群, 分别命名为Ccpla2g3a1、Ccpla2g3a2、Ccpla2g3b、Ccpla2g3c1和Ccpla2g3c2。基因结构和序列分析显示Ccpla2g3as含有7个外显子, Ccpla2g3b和Ccpla2g3cs则均含有4个外显子, 分别编码530、525、461、752和753个氨基酸, 均含有PLA₂_bee_venom_like region和sPLA₂-III特征序列。线性分析显示, Ccpla2g3a1和Ccpla2g3a2, Ccpla2g3c1和Ccpla2g3c2来自鲤特异的染色体加倍事件。从2R (Two rounds of genome duplication)鱼雀鳝 (Lepisosteus oculatus)到3R(Fish-specific genome duplication, FSGD)鱼斑马鱼(Danio rerio)等, 雀鳝的pla2g3a(b)加倍为3R鱼的pla2g3a和pla2g3b, pla2g3c因在3R鱼的一个连锁群上丢失, 故3R鱼中只保留一个基因。同源性分析表明已有鱼类Pla2g3a、Pla2g3b和Pla2g3c垂直同源蛋白之间相似性分别为48.8%—93.2%、37.6%—74.3%和49.6%—97.6%, 鲤和斑马鱼的垂直同源基因之间相似性最高。系统发育结果显示, 鱼类及其祖先雀鳝的Pla2g3c以100%的置信值归在一支, 雀鳝Pla2g3a(b)与除了鲤科鱼类 (鲤和斑马鱼) Pla2g3b外的Pla2g3a和Pla2g3b以96%置信值在一支, 鲤和斑马鱼Pla2g3b单独形成一支表明在进化过程中该基因变异较大。荧光定量PCR结果表明Ccpla2g3as在整个鲤早期发育阶段表达量均较低, 在成鱼肝脏中表达量最高 (P<0.01); Ccpla2g3b在鲤受精后0.5h的表达量显著高于之后各取样点 (P<0.05), 在成鱼卵巢中表达量最高 (P<0.01); Ccpla2g3cs在120h表达量达到最高, 在成鱼脑中表达量最高 (P<0.01)。实验比较全面地揭示了鲤sPLA₂-III亚家族基因结构、系统发育和表达特征。     

Tetraodon genome analysis provides further evidence for whole-genome duplication in the ray-finned fish lineage

A whole-genome duplication in the ray-finned fish lineage has been supported by the analyses of the genome sequence of the Japanese pufferfish, Fugu rubripes. Recently, genome sequence of a second teleost fish, the freshwater pufferfish, Tetraodon nigroviridis, was completed. Comparisons of long-range synteny between the Tetraodon and human genomes provided additional evidence for the whole-genome duplication in the ray-finned fish lineage. In the present study, we conducted phylogenetic analysis of the Tetraodon and human proteins to identify ray-finned fish lineage-specific (‘fish-specific’) duplicate genes in the Tetraodon genome. Our analyses provide evidence for 1087 well defined fish-specific duplicate genes in Tetraodon. We also analyzed the Fugu proteome that was predicted in the recent Fugu genome assembly, and identified 346 duplicate genes in addition to the 425 duplicates previously identified. We estimated the ages of duplicate genes using the molecular clock. The ages of duplicate genes in the two pufferfishes independently support a large-scale gene duplication around 380–400 Myr ago. In addition, a burst of recent gene duplications was evident in the Tetraodon lineage. These findings provide further evidence for a whole-genome duplication early in the evolution of ray-finned fishes, and suggest that independent gene duplications have occurred recently in the Tetraodon lineage.     

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

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

Multiple Kisspeptin Receptors in Early Osteichthyans Provide New Insights into the Evolution of This Receptor Family

Deorphanization of GPR54 receptor a decade ago led to the characterization of the kisspeptin receptor (Kissr) in mammals and the discovery of its major role in the brain control of reproduction. While a single gene encodes for Kissr in eutherian mammals including human, other vertebrates present a variable number of Kissr genes, from none in birds, one or two in teleosts, to three in an amphibian, xenopus. In order to get more insight into the evolution of Kissr gene family, we investigated the presence of Kissr in osteichthyans of key-phylogenetical positions: the coelacanth, a representative of early sarcopterygians, the spotted gar, a non-teleost actinopterygian, and the European eel, a member of an early group of teleosts (elopomorphs). We report the occurrence of three Kissr for the first time in a teleost, the eel. As measured by quantitative RT-PCR, the three eel Kissr were differentially expressed in the brain-pituitary-gonadal axis, and differentially regulated in experimentally matured eels, as compared to prepubertal controls. Subfunctionalisation, as shown by these differences in tissue distribution and regulation, may have represented significant evolutionary constraints for the conservation of multiple Kissr paralogs in this species. Furthermore, we identified four Kissr in both coelacanth and spotted gar genomes, providing the first evidence for the presence of four Kissr in vertebrates. Phylogenetic and syntenic analyses supported the existence of four Kissr paralogs in osteichthyans and allowed to propose a clarified nomenclature of Kissr (Kissr-1 to -4) based on these paralogs. Syntenic analysis suggested that the four Kissr paralogs arose through the two rounds of whole genome duplication (1R and 2R) in early vertebrates, followed by multiple gene loss events in the actinopterygian and sarcopterygian lineages. Due to gene loss there was no impact of the teleost-specific whole genome duplication (3R) on the number of Kissr paralogs in current teleosts.     

Duplication of phospholipase C-delta gene family in fish genomes

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