Skip to main content

Early vertebrate chromosome duplications and the evolution of the neuropeptide Y receptor gene regions

Abstract

Background

One of the many gene families that expanded in early vertebrate evolution is the neuropeptide (NPY) receptor family of G-protein coupled receptors. Earlier work by our lab suggested that several of the NPY receptor genes found in extant vertebrates resulted from two genome duplications before the origin of jawed vertebrates (gnathostomes) and one additional genome duplication in the actinopterygian lineage, based on their location on chromosomes sharing several gene families. In this study we have investigated, in five vertebrate genomes, 45 gene families with members close to the NPY receptor genes in the compact genomes of the teleost fishes Tetraodon nigroviridis and Takifugu rubripes. These correspond to Homo sapiens chromosomes 4, 5, 8 and 10.

Results

Chromosome regions with conserved synteny were identified and confirmed by phylogenetic analyses in H. sapiens, M. musculus, D. rerio, T. rubripes and T. nigroviridis. 26 gene families, including the NPY receptor genes, (plus 3 described recently by other labs) showed a tree topology consistent with duplications in early vertebrate evolution and in the actinopterygian lineage, thereby supporting expansion through block duplications. Eight gene families had complications that precluded analysis (such as short sequence length or variable number of repeated domains) and another eight families did not support block duplications (because the paralogs in these families seem to have originated in another time window than the proposed genome duplication events). RT-PCR carried out with several tissues in T. rubripes revealed that all five NPY receptors were expressed in the brain and subtypes Y2, Y4 and Y8 were also expressed in peripheral organs.

Conclusion

We conclude that the phylogenetic analyses and chromosomal locations of these gene families support duplications of large blocks of genes or even entire chromosomes. Thus, these results are consistent with two early vertebrate tetraploidizations forming a paralogon comprising human chromosomes 4, 5, 8 and 10 and one teleost tetraploidization. The combination of positional and phylogenetic data further strengthens the identification of orthologs and paralogs in the NPY receptor family.

Background

The evolutionary relationships of the NPY-receptor family receptors in vertebrates have been difficult to resolve due to uneven evolutionary rates and because some subtypes are missing in some classes of vertebrates. By using information on chromosomal location, initially in pig and human [1, 2], we suggested that chromosome duplications could account for the origin of several new family members. However, the relationships of the bony fish receptors called Y8a and Y8b, discovered in zebrafish and initially named Yc and Yb [3, 4], respectively, remained speculative [5] because they seemed to lack mammalian and bird orthologs.

Gene duplication by tetraploidization in the chordate lineage was proposed by Susumu Ohno in 1970 [6], based upon chromosome numbers and DNA content in different lineages. The first gene mapping data supporting a tetraploidization scenario emerged in 1987 when two human Hox clusters were mapped to human chromosomes Hsa7 and Hsa17 (Hsa for Homo sapiens) which resembled one another also with regard to other gene families [7]. Lundin described similarities in the other two Hox-bearing chromosomes, thereby identifying a quartet of related regions [8, 9]. The Hox chromosomes are now known to have involved duplication of more than 50 gene families [10–12].

In addition to the Hox-chromosome similarities, Lundin also reported resemblance within three other groups of human chromosomes. One group consisted of Hsa4 and Hsa5 [9], later found to contain NPY receptor genes [1] and extended to include Hsa8 and Hsa10 [13, 14]. Relationships between other chromosomes have been described by several authors, see for instance [11, 15–21]. Such groups of related, or paralogous, chromosome regions are called paralogons [22]. In tetrapod vertebrates, the paralogons are often comprised of quartets, consistent with a double tetraploidization scenario, called 2R for two rounds of genome doubling, before the origin of gnathostomes (jawed vertebrates) [23] although it is difficult to ascertain that the complete genome was quadrupled. Indeed, some regions do not have any paralogous counterparts [24]. More recently, a third tetraploidization has been identified in euteleost fish [25–28]. Several additional tetraploidizations have been described in specific lineages of for example fish and amphibians [29–32].

The sizes of the quadrupled paralogous gene regions have been difficult to determine because of numerous chromosomal rearrangements during the approximately 500 Myr since the tetraploidizations. Several vertebrate genome projects have recently been reported or are in progress, but due to incomplete assembly of the sequences into contigs or scaffolds, let alone chromosomes, these cannot always be used to analyze conserved synteny or paralogous gene regions. Another complicating factor has been the uneven divergence rates in some of the daughter genes after the duplications [10, 33–35] thereby aggravating the dating of the duplications. Indeed, inconsistent gene family phylogenies have been used as an argument against the tetraploidization hypothesis [36], although this can be seen as a natural consequence of uneven selection pressures or uneven re-diploidization rates after the two tetraploidizations, particularly as these may have taken place very close in time [10, 35, 37, 38].

Our laboratory has previously reported that the genes encoding NPY (neuropeptide Y)-family receptors, which belong to the superfamily of rhodopsin-like G-protein-coupled receptors (GPCRs), are located in the paralogon comprised of the human chromosomes Hsa4, 5 and 10 [1]. The fourth original chromosome member was shown to be partially represented by Hsa8 and Hsa2 [1, 14], although neither of the latter two chromosomes harbors NPY receptor genes. Our observation was based on a comparison of the human, mouse and pig chromosome regions [1] and has subsequently been supported by our analysis of the chicken NPY receptor genes [39]. However, neither the organization of NPY receptor genes in the recently reported euteleost fish genomes nor the extent of the chromosome regions comprising this paralogon or the phylogenetic relationships of the gene families involved, have been analyzed in detail.

We report here studies of 45 gene families whose members are located on Hsa 4, 5, 8/2/7 and 10. We have investigated conservation of synteny in human, mouse and three euteleost fishes, starting with the compact genomes of Tetraodon. nigroviridis and Takifugu. rubripes, and performed phylogenetic analyses of these gene families. This approach has been named "transitive homology" [40, 41] and allows for the identification of paralogous chromosomal segments despite the frequent loss of genes or rearrangement of gene order along chromosomes. The combined results of phylogenetic analyses and chromosomal locations reveal duplications of large chromosomal regions and are consistent with two basal vertebrate tetraploidizations, as well as the third tetraploidization in euteleosts. This analysis helps to clarify the evolutionary history of the chromosomal regions harboring the vertebrate NPY receptor genes and also further facilitates orthology/paralogy assignments of genes in the NPY receptor gene family.

Results

NPY receptors in Takifugu rubripes and Tetraodon nigroviridis

Both the T. nigroviridis genome and the T. rubripes genome was confirmed to harbor five NPY receptor genes previously found in D. rerio namely Y4 (Ya) [42], Y8a (Yc) [3] and Y8b (Yb) [4] belonging to the Y1 subfamily of receptors and Y2 and Y7 belonging to the Y2 subfamily. In D. rerio a Y1 receptor has been found recently that is not represented in the pufferfish genome databases [43]. So far no Y5 gene has been found in any of these three well-studied teleost fish species. The phylogenetic tree used to assign subfamily membership is shown in Fig. 1.

Figure 1
figure 1

Phylogenetic neighbor-joining tree of NPY receptors and closely related GPCRs. Numbers below branches show percent bootstrap support for each node. Nodes with values below 50 have been collapsed. Human bradykinin receptor B1 was used to root the tree.

RT-PCR and annotation of pufferfish genes

The RT-PCR carried out on eleven tissues in T. rubripes (Fig. 2) showed expression in brain and eye for all five receptors. Y8b was expressed in all eleven tissues while the other four receptors showed a more narrow expression pattern. Y8b also showed two distinct bands in some tissues because of alternative splicing. Interestingly, Y7 showed expression in several tissues in contrast to the very narrow expression observed for chicken Y7 [39].

Figure 2
figure 2

RT-PCR in Takifugu rubripes. Agarose gel showing the expression of the five Takifugu rubripes receptor genes in eleven different tissues. Actin was used as control to verify quality and content of samples.

The sequenced RT-PCR products revealed that the Y4 sequence contained an extended second extracellular loop in T. rubripes (Fig. 3A). This extension was also seen in the T. nigroviridis Y4 sequence obtained from the database. In addition, Y8a in T. rubripes has three novel introns. One intron comprising 102 bp is located in the end of the first extracellular loop. A second intron of 2.5 kb is present in the region encoding the middle of extracellular loop 2, and this intron contains an additional short (63 bp) exon that accounts for the extension of extracellular loop 2. A third 182 bp intron has been inserted in the region encoding intracellular loop 3 (Fig. 3B). Comparison with the T. nigroviridis Y8a sequence showed it to have the same overall organization (the zebrafish Y8a gene lacks introns in the coding region). The total length of Y8a is 450 aa in T. nigroviridis and 452 aa in T. rubripes. The T. rubripes Y8b gene lacks introns in the coding region apart from a cryptic intron spliced in some tissues (see Fig. 2 and 3C). The T. rubripes receptor sequences have been deposited to GenBank with the following accession numbers: [GenBank:EU104001 (Y2), GenBank:EU104002 (Y4), GenBank:EU104003 (Y7), GenBank:EU104004 (Y8a) and GenBank:EU104005 (Y8b)].

Figure 3
figure 3

Schematic picture of Takifugu rubripes NPY receptors. Picture depicting the Y4 (panel A), Y8a (panel B) and Y8b (panel C) receptors in Takifugu rubripes, indicating the extended second extracellular loop in Y4, the three extra introns present in the Y8a receptor gene and the alternative splicing of Y8b.

Conserved synteny and paralogous regions

A total of 35 gene families were found with members in three or four of the regions harboring the T. nigroviridis NPY receptor genes. In addition to these, another 9 families were included due to linkage to the NPY receptors in T. rubripes. Phylogenetic analysis of these 45 families (including the NPY receptor family) confirmed 26 to be compatible with an expansion in vertebrate evolution before the origin of gnathostomes. Neighbor-joining trees are shown in Fig. 4A–F and 5A–F for twelve of these families. The full set of NJ and ML-trees are available as Additional files 1 and 2, respectively. The total number of genes located in the investigated regions linked to the NPY receptor genes in T. nigroviridis was 556 (Y4), 310 (Y8a), 375 (Y8b), 487 (Y2) and 370 (Y7) according to the 34.1d version of the Ensembl database. This number of genes (2098) represents approximately 7.5% of the total gene number in the T. nigroviridis genome (total gene number estimated to be 28005 in this release of the database). The corresponding human orthologs situated on chromosomes 4, 5, 8/2/7 and 10 are dispersed over a large portion of these chromosomes. The synteny group associated with T. nigroviridis Y8a/Y8b seems to have been broken up in the human genome because some families are located on chromosome 8 and a few are on chromosomes 2 and 7. In addition to the 35 gene families with at least three members on the T. nigroviridis chromosomes, 127 gene families were identified that are represented on two of the chromosomes [see Additional file 3]. The conservation of synteny for the four chromosomes in human and mouse and eight chromosomes in the three fish species is illustrated in Fig. 6, 7, 8, 9. A schematic view of the evolution of 16 of the investigated gene families and the NPY receptor family is shown in Fig. 10.

Figure 4
figure 4

A-F – Phylogenetic analysis of neighboring gene families. Examples of the phylogenetic relationship for six neighboring gene families A (ABLIM), B (ADAMTS), C (AP3), D (ZIMP), E (HNRNP) and F (LGI). The trees were constructed using the neighbor-joining method as implemented in MEGA 3.1 with pair-wise deletion of gaps and poisson-corrected distances [117] All datasets were bootstrapped 1000 times (percent bootstrap support are indicated below each node). Species abbreviations are: Dme (Drosophila melanogaster), Cel (Caenorhabditis elegans), Cin (Ciona intestinalis), Csa (Ciona savignyi), Tru (Takifugu rubripes), Tni (Tetraodon nigroviridis), Dre (Danio rerio), Mmu (Mus musculus), Hsa (Homo sapiens). Numbers refers to chromosome number or scaffold number (ScXX), letters after chromosomal/scaffold number was arbitrarily assigned to family members located on the same chromosome to tell them apart. The scale below every tree indicates numbers of substitutions/site.

Figure 5
figure 5

A-F – Phylogenetic analysis of neighboring gene families. A-F. Examples of the phylogenetic relationship for six additional neighboring gene families A (MAX), B (NKR), C (OGDH), D (PX19), E (SORB) and F (TSPAN). Trees were constructed and visualized in the same way as in Fig. 4.

Figure 6
figure 6

Conserved synteny among vertebrate species compared to human chromosome 4. Conservation of synteny for genes investigated in this study residing on human chromosome 4 compared to mouse, T. rubripes, T. nigroviridis and D. rerio. Numbers above boxes denotes chromosome numbers or scaffold names. Position of the genes (Mega base pairs) is given below each box. Genes are ordered according to their positions on the human chromosome. Genes linked in the other species are indicated with lines above boxes. Species name abbreviations are: Hsa (Homo sapiens), Mmu (Mus musculus), Tni (Tetraodon nigroviridis), Tru (Takifugu rubripes) and Dre (Danio rerio). The NPY receptors are indicated in darker color while loss of NPY receptor genes are indicated with striped boxes. White boxes indicate genes where the phylogenetic analysis is inconclusive but where position indicates orthology.

Figure 7
figure 7

Conserved synteny among vertebrate species compared to human chromosome 5. Conservation of synteny for genes investigated in this study residing on human chromosome 5 compared to mouse, T. rubripes, T. nigroviridis and D. rerio. Chromosomes, scaffolds, gene positions, gene order and species abbreviations are given in the same way as in Fig. 6.

Figure 8
figure 8

Conserved synteny among vertebrate species compared to human chromosomes 8, 2 and 7. Conservation of synteny for genes investigated in this study residing on human chromosome 8, 2 and 7 compared to mouse, T. rubripes, T. nigroviridis and D. rerio. Chromosomes, scaffolds, gene positions, gene order and species abbreviations are given in the same way as in Fig. 6.

Figure 9
figure 9

Conserved synteny among vertebrate species compared to human chromosome 10. Conservation of synteny for genes investigated in this study residing on human chromosome 10 compared to mouse, T. rubripes, T. nigroviridis and D. rerio. Chromosomes, scaffolds, gene positions, gene order and species abbreviations are given in the same way as in Fig. 6.

Figure 10
figure 10

Proposed evolutionary scenario for the NPY receptor gene chromosomes. Schematic picture showing the proposed evolutionary history for 16 gene families located on the same chromosomes as the NPY receptor genes in T. nigroviridis and human. Note that position of genes is shuffled to simplify the picture. Striped boxes indicate gene losses. Position for genes on Hsa2 and Hsa7 is indicated in the boxes.

Statistical testing of paralogous regions

The statistical test of the investigated paralogous regions based on a binomial test as used previously by Vienne et al. [14] shows that their positions are different from a random distribution (P << 0.05) in both the genomes of T. nigroviridis (98 paralogs in total with 26 outside the investigated regions) and human (74 paralogs with 5 paralogs outside of the investigated regions). These results are in agreement with previous results investigating other genes present on Hsa 4, 5, 8 and 10 [14].

Summary of gene families analyzed in this study

Descriptions of the 25 subfamilies linked to the NPY receptor genes investigated in detail in this study are given below, i.e., the families that have members on 3 or 4 of the paralogous chromosomal regions in T. nigroviridis and families known to be linked to the NPY receptor genes in the T. rubripes genome. Three families that also support block or chromosome duplications were recently analyzed by other groups and are therefore not described here, namely two subfamilies of adrenergic receptors and the ADAM family of metzincins [44, 45].

Eight families with members on these chromosomes showed a tree topology seemingly inconsistent with expansion in early vertebrate evolution i.e. with several outgroup sequences dispersed among the vertebrate sequences, suggesting earlier origin of these paralogs. In addition to these families, eight families had to be left out of the analysis because of short sequence length, high conservation and therefore uninformative alignments or difficulties to generate reasonable alignments due to varying numbers of repeated domains. For a complete list of all families in the regions investigated (including the 127 families with only 2 members in T. nigroviridis and families studied due to members being present on the NPY receptor scaffolds in T. rubripes) see Additional file 3. For several of the analyzed families the number of fish sequences predicted to be part of the family according to the Ensembl database is higher than the number in our final trees. This is because we had to exclude some family members that lack one or several domains in order to obtain alignments of sufficient quality for phylogenetic analysis. Thus the number of fish paralogs in some of the families analyzed is a conservative representation of the actual number. We believe that this will be improved with more refined versions of the fish genomes (especially the Danio rerio genome) as well as comparison with additional genomes not available at the onset of this study.

ABLIM

The actin-binding LIM family has four members in the human genome. They are characterized by the presence of 4 LIM-domains and one villin domain. These proteins have been implicated in modulation of cell shape and cell differentiation through interaction with the actin cytoskeleton [46, 47]. The orthologous protein found in C. elegans has also been shown to interact with actin and has been shown to mediate axon guidance [48]. Expression of the different members of this protein family has been observed in distinct areas of the nervous system for example retina (ABLIM-1) [46], caudate/putamen (ABLIM-2), olfactory bulb (ABLIM-3), hippocampus (ABLIM-2 and ABLIM-3) and cerebellum (ABLIM-2 and ABLIM-3) [47]. The topology of the tree is in agreement with expansion of this gene family in the vertebrate lineage.

ADAMTS

The ADAMTS protein family (a disintegrin-like and metalloprotease protein with trombospondin motifs) includes proteins with a metalloprotease domain, disintegrin and spacer domains and a number of trombospondin repeats. Members of this protein family have been implicated in several diseases [49]. The human genome has so far been shown to contain more than 20 members of this gene family. The initial phylogenetic analysis identified a subfamily with members on human chromosomes 4 (ADAMTS3), 5 (ADAMTS2) and 10 (ADAMTS14), also identified in other studies [49]. The final analysis of this subfamily is in agreement with an expansion of this gene family in the vertebrate lineage before the divergence of actinopterygians and sarcopterygians.

Ankyrin (ANK)

The ankyrin gene family is one of the families earlier found to be part of a paralogon [1, 14] with members on human chromosomes 4, 8 and 10. We recovered a similar topology to previous studies with high bootstrap support. In addition the tree showed a topology consistent with extra duplication events in the teleost lineage, i.e. two fish genes corresponding to one tetrapod gene. The 33 aa ankyrin repeat is one of the most common protein domains in the pfam database.

Annexin (ANX)

A core domain consisting of four repeated units, each about 70 amino acids long, characterizes the annexin family. Each repeat unit contains 5 alpha helices that usually contain a "type 2" motif for binding of Ca2+ ions [50]. The phylogenetic analysis recovered a topology expected by expansion in the early vertebrate lineage and also evidence for local duplications before the origin of vertebrates as well as additional fish specific and a few tetrapod specific duplications. In addition, the positions of the genes investigated in this study agree with block duplications of the chromosomal regions harboring these genes.

AP3

This family contains two members in the human genome (AP3M1 and AP3M2) on Hsa8 and Hsa10 consistent with early vertebrate block duplication. They are members of a larger family of proteins called adaptins that are important for vesicular transport common to all eukaryotes [51, 52]. The mammalian AP-3 complex has been demonstrated to interact with clathrin and has been implicated in lysosomal membrane protein trafficking, sorting of melanosomal proteins and neurotransmitter/synaptic vesicle formation [51].

CNNM

The CNNM protein family (also referred to as ACDP) contains four members in both the human and the mouse genomes [53, 54]. These genes share one highly conserved domain, the ACD domain, which is present in a large number of species. The functions of these proteins are largely unknown. Immunofluorescence studies showed all four members to be localized to the nucleus and it has been speculated that these genes are probably involved in cell cycle regulation due to their similarity to cyclins [53]. The topology of the tree indicates a local duplication before the split of actinopterygians and sarcopterygians followed by block duplications and also some additional local duplicates in the fish lineage.

DUSP

The dual specificity phosphatases (or MAP kinase phosphatases) are evolutionarily conserved enzymes that are important in the regulation of apoptosis, cell proliferation and cell differentiation. They exert their effects by dephosphorylating and thereby inactivating MAP kinases [55]. The recovered topology in this study is in agreement with an expansion of this gene family in the vertebrate lineage. However the relationship of the different paralogs was not similar to that proposed earlier [14] possibly due to inclusion of other species in this analysis. Based on our analysis it is unlikely that the mouse gene named Dusp13 is the ortholog of human DUSP13.

Fibroblast growth factor receptor (FGFR)

The fibroblast growth factor receptors comprise one of the families already described as being a part of this paralogon [1, 14]. The topology of the phylogenetic tree is consistent with an expansion in the vertebrate lineage followed by further expansion in the teleost lineage. Interestingly, one of the families included in our initial list of gene families linked to the NPY receptor genes in T. nigroviridis is a subfamily of the fibroblast growth factors (FGFs) suggesting that both ligands and receptors were duplicated in the same time window. Hurst and Lercher observed that many ligand genes are linked to their receptor genes (including several of the FGFs and FGFRs) and that this could be due to block duplication that would keep the ratio of gene products at the same level as before duplication [56]. Our analysis of the FGF family was inconclusive in assigning orthology for the different members due to the rather short length of the sequences (~120aa). For detailed overviews of FGF and FGFR evolution with discussion about the role of gene and genome duplication see [57–59].

HNRNP

The HNRNP family is a family of proteins that contain two RBD or RRM domains. These proteins are known to interact with telomeric repeats d(TTAGGG) and to 3'splice sites r(UUAG/G) and have also been proposed to have a role in the regulation of mRNA stability. One of the human members, HNRPD has previously been mapped to Hsa 4q21 and been shown to have a highly conserved ortholog in the mouse genome [60]. Our present analysis recovered a subfamily with two members on Hsa4 and one member on Hsa5. The relationship of the fish sequences included was not resolved by this analysis, possibly due to loss of family members in human and mouse.

LGI

The gene LGI1 (leucine-rich gene – glioma inactivated) was discovered by positional cloning in 1998 and found to be mainly expressed in neural tissues, particularly in the brain. The gene is localized to human chromosome 10q24, a region found to be rearranged or deleted in several types of malignant brain tumors. Because of this LGI1 was proposed to be a tumor suppressor gene [61]. Three additional genes with similarity to LGI1 were found in the human genome and it was shown that these four genes constitute a subfamily of leucine rich repeat (LRR) genes [62]. The topology of the tree is consistent with duplications before the sarcopterygian-actinopterygian split, but one of the genes is located on Hsa19 rather than Hsa5.

MAPK

A subfamily of the large superfamily of mitogen-activated protein kinases (MAPK) was identified as being part of this paralogon. This subfamily contains the human members MAPK8, MAPK9 and MAPK10 (JNK1, JNK2 and JNK3). These proteins are involved in a wide variety of cellular processes like cell growth, proliferation, differentiation, immunity and development [63]. For a recent review of their role in vertebrate development see [64]. The gene duplications agree well with the block duplications prior to the sarcopterygian-actinopterygian divergence.

MXD

The MXD protein family has four members in the human genome. These proteins are characterized by a common basic-helix-loop-helix leucine zipper domain necessary for the formation of heterodimers with other proteins, such as Myc and Max [65]. The function of these proteins as regulators of Myc and Max activity with implications for tumour formation has recently been reviewed [66]. The topology of the phylogenetic tree and the positional information strongly supports the paralogy of the studied regions. It does not contain any extra fish specific duplicates.

NEF (Intermediate filaments/neurofilaments)

All intermediate filaments have a similar structural organization with a central alpha-helical rod domain that begins and ends with highly conserved aa-motifs necessary for correct assembly. Intermediate filament proteins have earlier been viewed as scaffolding structural proteins but more and more data suggest much more dynamic roles for these proteins [67]. The phylogeny of this protein family supports expansion in the vertebrate lineage. Interestingly, it seems like one subfamily in this gene family belongs to the HOX paralogon (with members on Hsa 2, 7, 12 and 17), which also contains the NPY peptide family. The tree also indicates a local duplication before the divergence of actinopterygians and sarcopterygians and additional fish specific duplications.

NKR

The tachykinin receptors are GPCRs that bind the amidated neuropeptides substance P, neurokinin A and neurokinin B as well as some other related peptides [68]. In vertebrates, three different tachykinin receptors have been described so far [68]. Our phylogenetic analysis of this receptor family shows a topology consistent with an expansion in the vertebrate lineage in agreement with the previous hypothesis [2], as well as an additional expansion consistent with block duplication in the teleost lineage.

OGDH

Oxoglutarate dehydrogenase (OGDH also designated E1k or lipoamide) makes up a part of the enzyme complex responsible for conversion of α-ketoglutarate to succinyl coenzyme A in the Krebs cycle [69]. The gene for OGDH and one related sequence (OGDHL) have previously been mapped to Hsa 7 and Hsa 10 respectively [70]. Both of the pufferfish species have an ortholog of OGDH and a presumed 3R co-ortholog (Fig. 8) and an OGDHL ortholog (Fig. 9). Oxoglutarate dehydrogenase has previously been shown to be involved in the production of reactive oxygen species in the brain of mice and thereby has been ascribed a role in neuronal cell death [71]. Apart from this, little is known about the function and evolution of this protein family.

PDLIM

The PDZ and LIM domain-containing family is a multi domain protein family characterized by the presence of one PDZ and one or several LIM domains. Both the ABLIM family described above and the PDLIM family described here belong to the same superfamily [72]. Our analysis uncovered a topology in accordance with expansion in the vertebrate lineage. This family contains a full quartet in the human genome as expected by 2R without loss of genes, although it does not display an (A,B)(C,D)-topology, presumably due to unequal evolutionary rates after duplication.

PX19

This family contains a conserved domain named MSF (after the protein MSF1 found in the yeast Saccaromyces cerevisiae). This family contains three closely related members of duplicated genes in the human lineage clearly duplicated after 2R but the two fish genes support the paralogon described here. The function of these proteins is unknown but they are thought to be involved in intra-mitochondrial protein sorting. One of the human genes is present on Hsa5 and earlier observed to be linked to one member of the MXD family (see above) [73].

SAMD8

The SAMD8 family (here named after one of its members) contains three different genes in the human genome: SAMD8 (referred to as sphingomyelin synthase-related protein 1 or sterile alpha motif domain-containing 8), TMEM23 or Mob (referred to as transmembrane protein 23, sphingomyelin synthase 1 or protein mob) and SMS2 (Sphingomyelin synthase 2). Proteins in this family all contain a SAM domain and 4–6 transmembrane domains and constitute a subset of a larger family of spingomyelin synthases [74]. Our analysis was indicative of an early local duplication event before the divergence between actinopterygians and sarcopterygians, with two genes still present on Hsa10 and one on Hsa4.

SFRP

The secreted frizzled related proteins (SFRP, also referred to as secreted apoptosis related proteins, SARPS) have five members in the human genome, three of which are located in the regions analyzed here. The SFRP family is characterized by a frizzled domain and a C-terminal netrin domain. Proteins of the SFRP family have been implicated in regulating Wnt-frizzled signalling either by interacting with Wnts or the frizzled receptors [75, 76]. Our phylogenetic analysis defines these three genes as a subgroup that has expanded during the evolution of vertebrates. Rattner et al. described this family in the mouse and showed that the members of this family was linked to some of the genes investigated in this study [77]. Of the fishes included here, only Danio rerio is indicative of any extra duplication.

SORB

The SORB or vinexin family has three members in the human genome namely vinexin, CAP (c-Cbl associated protein)/ponsin and ArgBP2 (Arg-binding protein 2). Common for all these proteins is that they contain a sorbin homology domain (pfamID PF02208) and three SH3 domains [78]. Members of this family have been implicated in the regulation of cell adhesion, and cytoskeletal organization and also growth factor signalling by functioning as adaptor proteins, connecting various other proteins [78]. This family supports an expansion in the proposed time window but does not show any evidence for a fish-specific duplication.

TIA

The two RNA binding proteins TIA1 (RNA-binding protein TIA-1) [79] and TIAR (TIA-1-related protein) [80] are the sole human members of a protein family that has been implicated in induction of apoptosis in certain cell types. Although the gene family only contains two members in the genomes of human and mouse, the phylogenetic and positional information support the conserved synteny of these chromosomes in this study. In the fishes this family supports the block duplication hypothesis for the chromosomes harboring the NPY8a and NPY8b receptors despite the failure to determine clear orthology relationships.

TSPAN

Tetraspanins are a large family of membrane proteins with four transmembrane domains (hence the name). These proteins are present in a wide variety of organisms and it has been proposed that the tetraspanin-like proteins present in plants share a common origin with tetraspanins in animals [81]. Tetraspanins are generally expressed in all cell types and usually several types are co-expressed. Their functions include various types of cell-cell and matrix-cell interactions and they have been implied in forming membrane microdomain structures and thereby working as "molecular facilitators" [82] or "molecular organizers" [81]. A total of 33 tetraspanins has previously been reported in human and 47 in bony fishes [81]. We identified a subfamily comprised of several vertebrate sequences as well as invertebrate outgroup sequences that support the block duplication hypothesis.

UNC5

The transmembrane protein UNC-5 was first characterized in C. elegans and implicated in regulation of netrin signalling since mutation of unc-5 causes neural migration defects while ectopic expression of unc-5 causes netrin-dependent redirection of axon growth in some neurons. Four different paralogs of UNC-5-like proteins have been found in vertebrates. These proteins contain two extracellular immunoglobulin-like domains and two extracellular trombospondin type 1 (TSP_1) domains [83] and three intracellular domains (ZU-5, DB and DD) [84]. The neuron guidance functions of this receptor family in response to netrins and possible interaction partners is reviewed in [84]. The topology of the tree is compatible with the 2R hypothesis and also shows some fish-specific duplicates. The fact that one rarely sees the perfect (A,B)(C,D)-topology for quartets of genes could be explained by unequal evolutionary rates following duplication (see also PDLIM above). However, when phylogenies and positional information are taken together the block duplication hypothesis is supported.

ZIMP

The ZIMP family of proteins contains a conserved SP-Ring/Miz domain, which they share with other PIAS proteins (protein inhibitors of activated STATS). In human, this family contains two members on Hsa7 and Hsa10 [85]. Both these proteins have been shown to regulate androgen receptor activity [86, 87]. Zimp10 has also been described to have a role in TGF-β/Smad signaling [88]. The phylogenetic analysis and position of the fish genes support the block duplication of the chromosomal segment harboring the Y8a and Y8b receptors.

ZINK

The ZINK family contains two proteins (zink finger protein 703 and zink finger protein 503) in the human genome. Not much is known about the functions of these proteins. They both contain a classical C2H2 zinc finger domain (pfam ID PF00096). The positional and phylogenetic information for these two genes are in agreement with block duplication before the actinopterygian-sarcopterygian split. The T. nigroviridis genome has two additional members resulting from the teleost fish tertaploidization on Tni2 and Tni Unr.

Discussion

The aim of the present study was to analyze in detail the phylogeny of gene families neighboring the NPY receptor genes to see if their evolutionary history was consistent with block/chromosome duplications [1, 2, 5]. We approached this by analyzing gene families on either side of the NPY-receptor family genes in the T. nigroviridis and T. rubripes genomes. We also used the information from the analyses of neighboring gene families to assign orthology and paralogy of the NPY receptor genes. The present study identified the receptors Y2, Y4, Y7, Y8a and Y8b in the two pufferfish genomes (Fig. 1). The chromosomal locations of the Y8a and Y8b genes show that they reside on two related fish chromosomes that most likely arose in the teleost tetraploidization (3R). Furthermore, the gene neighbors on these chromosomes strongly suggest that the corresponding ancestral teleost chromosome belongs to a quartet of ancestral gnathostome chromosomes that most likely arose in the proposed basal vertebrate tetraploidizations (Fig. 6, 7, 8, 9). The zebrafish Y8a gene is located on chromosome 17 according to the database version used in this study but has previously been mapped to chromosome 10 [42] (see Fig. 8). Thus, we can confirm and extend the proposed gene duplication scheme [5] so that it accounts for all of the NPY-family receptors in mammalian and teleost genomes (Fig. 10): an ancestral local triplication was followed by the basal vertebrate tetraploidizations whereupon several genes were probably lost, resulting in an ancestral gnathostome repertoire of seven NPY receptor genes.

In the tetrapod lineage, Y7 and Y8 are present in the frog Xenopus tropicalis (unpublished data) and Y7 is still present in chicken [39] while both the Y7 and Y8 genes seem to have been lost in the lineage leading to mammals. In the actinopterygian lineage leading to euteleosts, one additional copy arose in the teleost 3R tetraploidization, while two genes seem to have been lost in euteleosts, Y5 and Y6. We recently discovered the Y1 gene in the zebrafish genome but it has so far not been found in the genomes of pufferfishes, medaka or stickleback and may have been lost. For a description of the NPY receptor repertoire in basal teleosts see Salaneck et al. [43]. Finally, also the prolactin releasing hormone receptor gene shares the same degree of identity with the NPY-family receptors as the Y1–Y2–Y5 subfamilies display to one another [89], but the preferred ligand is not NPY and therefore we have not included it in the duplication scheme.

Our characterization of the five T. rubripes genes revealed that both Y4 and Y8a have received insertions in the coding region after divergence from the lineage leading to zebrafish, as seen for the melanocortin receptor genes MC2R and MC5R in pufferfish [90], despite the intron and general genome compaction of these genomes [91]. In the case of Y4 the insertion extends extracellular loop 2 with 74 amino acids (Fig. 3A). The Y8a gene has undergone no less than three insertions of introns whose positions are projected onto the protein structure shown in Fig. 3B. The presence of the insertion in Y4 mRNA and the removal of the three Y8a introns were confirmed by RACE and PCR in T. rubripes. The largest of the Y8a introns is 2.5 kb and carries a small exon that extends extracellular loop 2 with 21 amino acids. In the Y8b gene one cryptic intron is spliced and this shortened splice variant is present in a minor proportion of the mRNAs and presumably leads to a nonfunctional partial receptor protein. This cryptic intron splice site is also present in the medaka gene found in the genome database. We suggest that the Y4 insertion was also probably a reinserted intron that subsequently lost its splice signals, but thanks to its small size and the maintained reading frame it could be tolerated as a protein expansion in extracellular loop 2. Functional expression of the receptors will be necessary to see if the extensions of loop 2 in Y4 and Y8a affect the ligand-binding properties of the receptors. Studies of the anatomical distribution of the mRNAs for the five NPY receptors in T. rubripes, as detected by RT-PCR, show that all five are present in the brain and eye. Interestingly, Y8a and Y8b differ greatly in their distribution in that Y8b is expressed in all organs investigated whereas Y8a shows the narrowest distribution of all receptors, although Y8a and Y8b are most closely related to each other having originated in 3R. Functional information is still missing for many fish NPY-family receptors. The pufferfish NPY system is more complicated than in most other vertebrates because there are four peptide ligands due to 3R duplicates of both NPY and its relative PYY [92].

Several authors have discussed what criteria one should have for identifying paralogous regions and to safely infer that block duplications or chromosome duplications have occurred [93–96]. We argue that as many species as possible, with divergencies close in time to the proposed chromosome duplication events, should be included in the analyses to be able to date duplications and to reveal fluctuations of evolutionary rates among gene duplicates and across lineages. This also helps identify translocations and inversions as well as lineage-specific duplications and deletions. Due to high frequency of chromosomal rearrangements, inter- as well as intra-chromosomal, the gene content of chromosomal blocks can be considered sufficient to identify duplications on these time scales (several hundred million years). Given the high rate of deletion after duplication we want to emphasize the importance of combining map-based and phylogenetic approaches in order to understand the evolution of genomic regions. In this way, gene families with only two members (as opposed to the expected 4 as predicted by 2R) still give important information.

It has been suggested that the observed patterns of paralogy within genomes, often interpreted to be the remnants of large-scale duplications, could be the result of convergent evolution [97]. We believe this alternative explanation to be unlikely not only because polyploidization has been shown to be common in many lineages but mainly because reconstructions of ancestral chromosomal regions based on comparisons of vertebrate genomes shows that these paralogous regions span large genomic regions in many species [21]. If many small independent events produced the observed patterns of chromosomal similarity, one would need to infer that these events occurred in a relatively short period of time before the vertebrate radiation, otherwise one has to invoke several independent but identical duplication and translocation events giving the same chromosomal organization in different lineages. To our knowledge, no mechanism has been described in metazoans that could support such a scenario.

Among the gene neighbors of the NPY receptor genes, we analyzed 44 gene families with members on the chromosomes bearing NPY-receptor genes. In fact, chromosomes orthologous to three of the four human chromosomes were all found to be present in duplicate in the fish genomes, consistent with the 3R teleost tetraploidization (see Fig. 7, 8 and 9). This pattern has been referred to as "doubly conserved synteny" [26]. The phylogenetic relationships of 25 of these gene families are consistent with concomitant block or chromosome duplications. Another three gene families included in the initial table [see Additional file 3], have earlier been proposed to be part of this paralogon and recently been studied in detail by others (two subfamilies of adrenergic receptors [44, 98] and members of the large gene superfamily of metzincins other than the ADAMTS family included in our analysis [45]).

Eight gene families had to be excluded from analysis because they had too many family members to achieve a reasonably clear picture of their evolutionary relationships or problems in analyzing them due to too high sequence conservation to be informative or varying number of protein domains making it hard to obtain unequivocal alignments. More thorough analysis using diagnostic positions and comparisons of intron-exon structure could be used to clarify the relation of the different paralogs of these families. Thus, our study provides a minimal estimate of the number of gene families that expanded simultaneously with the NPY receptor genes.

Out of the 44 gene families a total of 26+3 support block/chromosome duplications whereas eight gene families are unsupportive of the block duplication scheme with invertebrate sequences located between the vertebrate paralogs. The reasons for this may be earlier expansions of these families, uneven evolutionary rates among the daughter genes, gene conversion between family members or simply translocation of distant family members to the chromosomes in this study.

We uncovered evidence for extra fish duplicates, in agreement with 3R, for 18 of these 26 gene families. However, because we were forced to exclude several gene family members as they lacked one or a few protein domains, this is certainly an underestimate of the true number (see above). It has been pointed out earlier that the number of genes retained after 3R is rather low as compared to after 1R and 2R [99]. One good example of this is the well studied Hox clusters that have lost many duplicates [100].

The fact that there are 127 gene families with members on two of the investigated T. nigroviridis chromosomes may further corroborate the conserved synteny between the chromosomal regions in this study, although they need to be analyzed phylogenetically as well. Among these 127 pairs of paralogs, all ten possible combinations of the five T. nigroviridis chromosomes are represented. The lack of sequence data from several important species is still a limiting factor in investigation of the early vertebrate chromosome duplications. In particular, sequences from the cephalochordate Branchiostoma floridae (first draft version of assembly recently released but with limited data on gene linkage) and representatives of the jawless and cartilaginous fishes will help to date the block/chromosome duplication events when more extensive gene contigs have been generated.

One interesting observation is that members of some of the gene families studied here have been reported to interact with each other in signaling networks or possess domains that commonly interact with each other. This opens for the possibility that blocks of genes that were duplicated simultaneously can co-evolve and therefore one might detect subfunctionalization of entire gene networks and not only of different paralogs within one family. Examples of such families according to the BioGRID database of protein interactions [101, 102] are the MAPK and the DUSP families, the ABLIM and NEF families, the ANX and NEF families, the CNNM and HNRNP families, the PDLIM and NEF families, the SAMD8 and ADAMTS families and the SAMD8 and the HNRNP families. However, these families contain many members and are known to have a diverse set of interaction partners and therefore functional experiments are needed to specifically test the interaction of the paralogs analyzed in this study. This observation also leads to speculation about differential retention of gene family members after duplication because of their different functions or involvement in signaling networks. A higher retention rate after duplication of genes and specialization of different paralogs that are tied up in signaling networks has been proposed earlier [103]. This also has been shown for genetic networks in yeast (Saccharomyces cerevisiae) duplicated by tetraploidization around 100 million years ago [104] as well as for duplicates in Arabidopsis thaliana produced by the most recent polyploidizations in this species that took place 20 to 60 million years ago [105, 106]. These processes are also linked to the predictions regarding gene fate after duplication deduced from models of subfunctionalization [107–110] stating that paralogs could be fixed in the genome by a partitioning of ancestral functions, something that would lead to a higher retention of gene duplicates that could subsequently evolve new functions [111, 112].

In earlier analyses of these paralogous regions using several other gene families, the order of duplication events has been inferred based on the phylogenetic trees giving highest support for the (Hsa8, Hsa10), (Hsa4, Hsa5) topology [14]. In our dataset this is also the most frequent relationship observed for the human chromosomes, possibly reflecting the order of duplication events.

Conclusion

In summary, we have characterized the NPY receptor repertoire in the two pufferfishes T. rubripes and T. nigroviridis and compared the chromosomal regions where the receptor genes are located in one additional fish and two mammalian species. The conserved synteny shows that many of the gene families were located together in the same chromosome regions of the common ancestor of gnathostomes more than 400 Myr ago [113] (for a summary see Fig. 10). Our results are in line with the tetraploidizations in early vertebrate evolution as well as an additional tetraploidization in teleosts. Although gene losses are frequent after duplication it is possible to infer paralogy and orthology in this way by analyzing both phylogentic and positional information simultaneously. This "transitive homology" approach [40, 41] in combination with dating of duplication events in relation to speciation events is in our opinion, as shown by the present study, a more reliable way to unravel the evolutionary history of gene families in cases where phylogenetic analyses alone are not fully informative.

Methods

Identification and analysis of NPY receptor genes in Tetraodon nigroviridis and Takifugu rubripes

BLAST searches were carried out on the Ensembl database version 35.1d [114] using human and zebrafish NPY receptor sequences in order to identify all NPY receptor sequences in the genomes of the two pufferfishes T. rubripes and T. nigroviridis. The sequences found were aligned with previously known NPY receptor sequences and closely related peptide binding receptors using the Windows version of ClustalX 1.81 [115, 116]. The alignment was manually edited to remove poorly aligned residues. Thereafter an initial phylogenetic analysis was performed using the neighbor-joining method in MEGA3.1 [117] with standard settings in order to assign which sub-families of receptors that were represented in the pufferfish genomes. In addition to the NJ-tree a quartet-puzzling tree was constructed with Treepuzzle 5.2 [118] using the same alignment. This tree was made with the following settings: 9 categories of sites (8 gamma + 1 invariant, parameters estimated from the dataset using the "exact-slow" option) with 10000 puzzling steps using the JTT substitution matrix (tree not shown).

The NJ-tree was bootstrapped 1000 times. Several closely related human sequences were included in the analyses (see Fig. 1). Human bradykinin B1 receptor was used to root the tree and nodes with bootstrap support values below 50% were collapsed (see Fig. 1).

RT-PCR in Takifugu rubripes

Total RNA was isolated from eleven T. rubripes tissues using TRIzol reagent (Invitrogen, USA) according to manufacturer's protocol. Purified total RNA was reverse transcribed and single-strand 5'RACE-ready cDNA was prepared using SMART RACE cDNA Amplification Kit (Clontech, USA). Primers used for the receptors as well as an internal actin control are listed in Table 1. The PCR was carried out according to the following protocol: a denaturation step at 95°C for 2 min, 35 cycles of 95°C 30 sec, 55°C for 1 min, 72°C for 1 min followed by a final elongation step at 72°C for 5 min. Identity of representative RT-PCR products were confirmed by sequencing on an Applied Biosystems 3700 DNA Analyzer using dye-terminator chemistry.

Table 1 Primers used for RT-PCR on eleven Takifugu rubripes tissues.

Phylogenetic analysis of neighboring genes

Starting from the confirmed receptor genes in the T. nigroviridis genome, all Ensembl Gene IDs for genes positioned four megabases on each side of the receptor genes were downloaded and saved in an Excel file. In addition to Ensembl Gene IDs the file contained information on Ensembl Family ID and family description as well as information about chromosomal position and human orthologs. This list was sorted based on Ensembl Family ID and all multiple entries for the same gene due to multiple transcripts were removed. The gene families containing members close to three or more of the T. nigroviridis NPY receptor genes were used for phylogenetic analysis. In addition to gene families found in this way, several families with members on the NPY receptor gene harboring scaffolds in T. rubripes were included [see Additional file 3]. All amino acid sequences of members included in the Ensembl families were downloaded from T. nigroviridis, T. rubripes, Danio rerio, Homo sapiens and Mus musculus. Invertebrate sequences representing at least one outgroup species were included in order to root the trees. For the main part of the trees both Drosophila melanogaster and Ciona intestinalis was used as outgroups. In cases where no clear orthologs could be found in these two species we used sequences from other invertebrate genomes available in the Ensembl database [114] (see figure legends and additional files for complete description of outgroups). Sequences for each family were aligned using MEGA 3.1 with default settings. All alignments were manually inspected and short and poorly aligned sequences were removed. Sequence alignments were further adjusted with the aid of the pfam database [119] for prediction of protein domains. Relevant literature describing the families was also used to find description of the domains. Thereafter the cut sequences were realigned and neighbor-joining trees was constructed using MEGA 3.1 for each family with pairwise deletion of gaps, 1000 bootstrap replicates and poisson-corrected distances. Alignments for the 26 families analyzed in detail are available [see Additional file 4].

Initial phylogenetic trees were inspected for topologies consistent with an expansion of the family in vertebrate evolution i.e. the outgroup sequence/sequences rooting several of the sequences residing on the particular chromosomes under study (for example Hsa 4, 5, 8 and 10). In cases where a large number of sequences were included in the Ensembl family, the initial phylogenetic tree was used to find relevant sub families represented by multiple vertebrate sequences and at least one outgroup sequence. The identified sub families were realigned and subjected to neighbor-joining analysis as described above. In addition to gene families represented on three or more of the T. nigroviridis chromosomes all families with members on two of the chromosomes were saved in a table [see Additional file 3]. Phylogenetic trees with a topology consistent with expansion in vertebrate evolution before the origin of gnathostomes using the NJ-method were examined using the Quartet-puzzling method in the Windows version of Treepuzzle 5.2 [118] to further support the result. The settings used for each analysis was the same as mentioned above for the NPY receptor family but with varying number of puzzling steps (1000–25000) depending on the size of the family.

Statistical testing of paralogous regions

The position of paralogs were statistically tested for random distribution in both the human and T. nigroviridis genomes using the binomial test as described earlier by Vienne et al. [14]. Genes belonging to the same family residing close to each other on the same chromosome that grouped together in the phylogenetic trees were counted as one member because they are most probably the result of recent lineage specific local duplications and thereby not of interest in our analysis. The chromosomes tested for random distribution were the ones shown by phylogenetic analysis to contain several paralogs (i.e human chromosomes 4, 5, 8/2/7 and 10). For details on statistical analysis: [see Additional file 5].

References

  1. Wraith A, Tornsten A, Chardon P, Harbitz I, Chowdhary BP, Andersson L, Lundin LG, Larhammar D: Evolution of the neuropeptide Y receptor family: gene and chromosome duplications deduced from the cloning and mapping of the five receptor subtype genes in pig. Genome research. 2000, 10: 302-310. 10.1101/gr.10.3.302.

    PubMed Central  CAS  PubMed  Google Scholar 

  2. Larhammar D, Fredriksson R, Larson ET, Salaneck E: Phylogeny of NPY-family peptides and their receptors. "Neuropeptide Y and Related Peptides", Handbook of Experimental Pharmacology. Edited by: Michel MC. 2004, Berlin-Heidelberg: Springer-Verlag, 75-100.

    Google Scholar 

  3. Ringvall M, Berglund MM, Larhammar D: Multiplicity of neuropeptide Y receptors: cloning of a third distinct subtype in the zebrafish. Biochem Biophys Res Commun. 1997, 241: 749-755. 10.1006/bbrc.1997.7886.

    CAS  PubMed  Google Scholar 

  4. Lundell I, Berglund MM, Starback P, Salaneck E, Gehlert DR, Larhammar D: Cloning and characterization of a novel neuropeptide Y receptor subtype in the zebrafish. DNA Cell Biol. 1997, 16: 1357-1363.

    CAS  PubMed  Google Scholar 

  5. Larhammar D, Salaneck E: Molecular evolution of NPY receptor subtypes. Neuropeptides. 2004, 38: 141-151. 10.1016/j.npep.2004.06.002.

    CAS  PubMed  Google Scholar 

  6. Ohno S: Evolution by gene duplication. 1970, Berlin: Springer Verlag

    Google Scholar 

  7. Hart CP, Fainsod A, Ruddle FH: Sequence analysis of the murine Hox-2.2, -2.3, and -2.4 homeo boxes: evolutionary and structural comparisons. Genomics. 1987, 1: 182-195. 10.1016/0888-7543(87)90011-5.

    CAS  PubMed  Google Scholar 

  8. Lundin L-G: Gene homologies and the nerve system. Genetics of Neuropsychiatric Diseases. Wenner-Gren Internatl. Symp. Ser. Edited by: Wetterberg L. 1989, London: MacMillan Press, 51: 43-58.

    Google Scholar 

  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. 10.1006/geno.1993.1133.

    CAS  PubMed  Google Scholar 

  10. Larhammar D, Lundin LG, Hallbook F: The human Hox-bearing chromosome regions did arise by block or chromosome (or even genome) duplications. Genome research. 2002, 12: 1910-1920. 10.1101/gr.445702.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Popovici C, Leveugle M, Birnbaum D, Coulier F: Homeobox gene clusters and the human paralogy map. FEBS letters. 2001, 491: 237-242. 10.1016/S0014-5793(01)02187-1.

    CAS  PubMed  Google Scholar 

  12. Abbasi AA, Grzeschik KH: An insight into the phylogenetic history of HOX linked gene families in vertebrates. BMC Evol Biol. 2007, 7: 239-10.1186/1471-2148-7-239.

    PubMed Central  PubMed  Google Scholar 

  13. Pebusque MJ, Coulier F, Birnbaum D, Pontarotti P: Ancient large-scale genome duplications: phylogenetic and linkage analyses shed light on chordate genome evolution. Mol Biol Evol. 1998, 15 (9): 1145-1159.

    CAS  PubMed  Google Scholar 

  14. Vienne A, Rasmussen J, Abi-Rached L, Pontarotti P, Gilles A: Systematic phylogenomic evidence of en bloc duplication of the ancestral 8p11.21–8p21.3-like region. Molecular biology and evolution. 2003, 20: 1290-1298. 10.1093/molbev/msg127.

    CAS  PubMed  Google Scholar 

  15. Popovici C, Leveugle M, Birnbaum D, Coulier F: Coparalogy: physical and functional clusterings in the human genome. Biochem Biophys Res Commun. 2001, 288: 362-370. 10.1006/bbrc.2001.5794.

    CAS  PubMed  Google Scholar 

  16. Hokamp K, McLysaght A, Wolfe KH: The 2R hypothesis and the human genome sequence. J Struct Funct Genomics. 2003, 3: 95-110. 10.1023/A:1022661917301.

    CAS  PubMed  Google Scholar 

  17. Lundin LG, Larhammar D, Hallbook F: Numerous groups of chromosomal regional paralogies strongly indicate two genome doublings at the root of the vertebrates. J Struct Funct Genomics. 2003, 3: 53-63. 10.1023/A:1022600813840.

    CAS  PubMed  Google Scholar 

  18. Abi-Rached L, Gilles A, Shiina T, Pontarotti P, Inoko H: Evidence of en bloc duplication in vertebrate genomes. Nature genetics. 2002, 31: 100-105. 10.1038/ng855.

    CAS  PubMed  Google Scholar 

  19. Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate. PLoS biology. 2005, 3: e314-10.1371/journal.pbio.0030314.

    PubMed Central  PubMed  Google Scholar 

  20. Olinski RP, Lundin LG, Hallbook F: Conserved synteny between the Ciona genome and human paralogons identifies large duplication events in the molecular evolution of the insulin-relaxin gene family. Molecular biology and evolution. 2006, 23: 10-22. 10.1093/molbev/msj002.

    CAS  PubMed  Google Scholar 

  21. Nakatani Y, Takeda H, Kohara Y, Morishita S: Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates. Genome research. 2007, 17: 1254-1265. 10.1101/gr.6316407.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Coulier F, Popovici C, Villet R, Birnbaum D: MetaHox gene clusters. J Exp Zool. 2000, 288: 345-351. 10.1002/1097-010X(20001215)288:4<345::AID-JEZ7>3.0.CO;2-Y.

    CAS  PubMed  Google Scholar 

  23. Furlong RF, Holland PW: Were vertebrates octoploid?. Philos Trans R Soc Lond B Biol Sci. 2002, 357: 531-544. 10.1098/rstb.2001.1035.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Bromee T, Venkatesh B, Brenner S, Postlethwait JH, Yan YL, Larhammar D: Uneven evolutionary rates of bradykinin B1 and B2 receptors in vertebrate lineages. Gene. 2006, 373: 100-108. 10.1016/j.gene.2006.01.017.

    CAS  PubMed  Google Scholar 

  25. 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. Molecular biology and evolution. 2004, 21: 1146-1151. 10.1093/molbev/msh114.

    CAS  PubMed  Google Scholar 

  26. 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, Biemont C, Skalli Z, Cattolico L, Poulain J, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.

    PubMed  Google Scholar 

  27. 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, et al: The medaka draft genome and insights into vertebrate genome evolution. Nature. 2007, 447: 714-719. 10.1038/nature05846.

    CAS  PubMed  Google Scholar 

  28. Postlethwait JH, Woods IG, Ngo-Hazelett P, Yan YL, Kelly PD, Chu F, Huang H, Hill-Force A, Talbot WS: Zebrafish comparative genomics and the origins of vertebrate chromosomes. Genome research. 2000, 10: 1890-1902. 10.1101/gr.164800.

    CAS  PubMed  Google Scholar 

  29. Uyeno T, Smith GR: Tetraploid origin of the karyotype of catostomid fishes. Science. 1972, 175: 644-646. 10.1126/science.175.4022.644.

    CAS  PubMed  Google Scholar 

  30. Allendorf FW, Thorgaard GH: Tetraploidy and the Evolution of Salmonid Fishes. The Evolutionary Genetics of Fishes. Edited by: Turner BJ. 1984, Plenum Press, 1-53.

    Google Scholar 

  31. Otto SP, Whitton J: Polyploid incidence and evolution. Annu Rev Genet. 2000, 34: 401-437. 10.1146/annurev.genet.34.1.401.

    CAS  PubMed  Google Scholar 

  32. Le Comber SC, Smith C: Polyploidy in fishes: patterns and processes. Biological Journal of the Linnean Society. 2004, 82: 431-442. 10.1111/j.1095-8312.2004.00330.x.

    Google Scholar 

  33. Wagner A: Asymmetric functional divergence of duplicate genes in yeast. Mol Biol Evol. 2002, 19 (10): 1760-1768.

    CAS  PubMed  Google Scholar 

  34. Conant GC, Wagner A: Asymmetric sequence divergence of duplicate genes. Genome research. 2003, 13: 2052-2058. 10.1101/gr.1252603.

    PubMed Central  CAS  PubMed  Google Scholar 

  35. Aburomia R, Khaner O, Sidow A: Functional evolution in the ancestral lineage of vertebrates or when genomic complexity was wagging its morphological tail. J Struct Funct Genomics. 2003, 3: 45-52. 10.1023/A:1022648729770.

    CAS  PubMed  Google Scholar 

  36. Hughes AL, da Silva J, Friedman R: Ancient genome duplications did not structure the human Hox-bearing chromosomes. Genome research. 2001, 11: 771-780. 10.1101/gr.GR-1600R.

    PubMed Central  CAS  PubMed  Google Scholar 

  37. Nadeau JH, Sankoff D: Comparable rates of gene loss and functional divergence after genome duplications early in vertebrate evolution. Genetics. 1997, 147: 1259-1266.

    PubMed Central  CAS  PubMed  Google Scholar 

  38. Gibson TJ, Spring J: Evidence in favour of ancient octaploidy in the vertebrate genome. Biochemical Society transactions. 2000, 28: 259-264.

    CAS  PubMed  Google Scholar 

  39. Bromee T, Sjodin P, Fredriksson R, Boswell T, Larsson TA, Salaneck E, Zoorob R, Mohell N, Larhammar D: Neuropeptide Y-family receptors Y6 and Y7 in chicken. Cloning, pharmacological characterization, tissue distribution and conserved synteny with human chromosome region. Febs J. 2006, 273: 2048-2063. 10.1111/j.1742-4658.2006.05221.x.

    CAS  PubMed  Google Scholar 

  40. Simillion C, Vandepoele K, Van Montagu MC, Zabeau M, Peer Van de Y: The hidden duplication past of Arabidopsis thaliana. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99: 13627-13632. 10.1073/pnas.212522399.

    PubMed Central  CAS  PubMed  Google Scholar 

  41. Vandepoele K, Simillion C, Peer Van de Y: Detecting the undetectable: uncovering duplicated segments in Arabidopsis by comparison with rice. Trends Genet. 2002, 18: 606-608. 10.1016/S0168-9525(02)02796-8.

    CAS  PubMed  Google Scholar 

  42. Starback P, Lundell I, Fredriksson R, Berglund MM, Yan YL, Wraith A, Soderberg C, Postlethwait JH, Larhammar D: Neuropeptide Y receptor subtype with unique properties cloned in the zebrafish: the zYa receptor. Brain Res Mol Brain Res. 1999, 70: 242-252. 10.1016/S0169-328X(99)00152-7.

    CAS  PubMed  Google Scholar 

  43. Salaneck E, Larsson TA, Larson ET, Larhammar D: Birth and death of neuropeptide Y receptor genes in relation to the teleost fish tetraploidization. Gene. 2008, 409: 61-71. 10.1016/j.gene.2007.11.011.

    CAS  PubMed  Google Scholar 

  44. Postlethwait JH: The zebrafish genome in context: ohnologs gone missing. Journal of experimental zoology Part B. 2007, 308: 563-577. 10.1002/jez.b.21137.

    Google Scholar 

  45. Huxley-Jones J, Clarke TK, Beck C, Toubaris G, Robertson DL, Boot-Handford RP: The evolution of the vertebrate metzincins; insights from Ciona intestinalis and Danio rerio. BMC Evol Biol. 2007, 7: 63-10.1186/1471-2148-7-63.

    PubMed Central  PubMed  Google Scholar 

  46. Roof DJ, Hayes A, Adamian M, Chishti AH, Li T: Molecular characterization of abLIM, a novel actin-binding and double zinc finger protein. The Journal of cell biology. 1997, 138: 575-588. 10.1083/jcb.138.3.575.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Barrientos T, Frank D, Kuwahara K, Bezprozvannaya S, Pipes GC, Bassel-Duby R, Richardson JA, Katus HA, Olson EN, Frey N: Two novel members of the ABLIM protein family, ABLIM-2 and -3, associate with STARS and directly bind F-actin. J Biol Chem. 2007, 282: 8393-8403. 10.1074/jbc.M607549200.

    CAS  PubMed  Google Scholar 

  48. Lundquist EA, Herman RK, Shaw JE, Bargmann CI: UNC-115, a conserved protein with predicted LIM and actin-binding domains, mediates axon guidance in C. elegans. Neuron. 1998, 21: 385-392. 10.1016/S0896-6273(00)80547-4.

    CAS  PubMed  Google Scholar 

  49. Nicholson AC, Malik SB, Logsdon JM, Van Meir EG: Functional evolution of ADAMTS genes: evidence from analyses of phylogeny and gene organization. BMC Evol Biol. 2005, 5: 11-10.1186/1471-2148-5-11.

    PubMed Central  PubMed  Google Scholar 

  50. Moss SE, Morgan RO: The annexins. Genome biology. 2004, 5: 219-10.1186/gb-2004-5-4-219.

    PubMed Central  PubMed  Google Scholar 

  51. Boehm M, Bonifacino JS: Genetic analyses of adaptin function from yeast to mammals. Gene. 2002, 286: 175-186. 10.1016/S0378-1119(02)00422-5.

    CAS  PubMed  Google Scholar 

  52. Odorizzi G, Cowles CR, Emr SD: The AP-3 complex: a coat of many colours. Trends Cell Biol. 1998, 8: 282-288. 10.1016/S0962-8924(98)01295-1.

    CAS  PubMed  Google Scholar 

  53. Wang CY, Shi JD, Yang P, Kumar PG, Li QZ, Run QG, Su YC, Scott HS, Kao KJ, She JX: Molecular cloning and characterization of a novel gene family of four ancient conserved domain proteins (ACDP). Gene. 2003, 306: 37-44. 10.1016/S0378-1119(02)01210-6.

    CAS  PubMed  Google Scholar 

  54. Wang CY, Yang P, Shi JD, Purohit S, Guo D, An H, Gu JG, Ling J, Dong Z, She JX: Molecular cloning and characterization of the mouse Acdp gene family. BMC genomics. 2004, 5: 7-10.1186/1471-2164-5-7.

    PubMed Central  PubMed  Google Scholar 

  55. Theodosiou A, Ashworth A: MAP kinase phosphatases. Genome biology. 2002, 3: REVIEWS3009-10.1186/gb-2002-3-7-reviews3009.

    PubMed Central  PubMed  Google Scholar 

  56. Hurst LD, Lercher MJ: Unusual linkage patterns of ligands and their cognate receptors indicate a novel reason for non-random gene order in the human genome. BMC Evol Biol. 2005, 5: 62-10.1186/1471-2148-5-62.

    PubMed Central  PubMed  Google Scholar 

  57. Coulier F, Pontarotti P, Roubin R, Hartung H, Goldfarb M, Birnbaum D: Of worms and men: an evolutionary perspective on the fibroblast growth factor (FGF) and FGF receptor families. Journal of molecular evolution. 1997, 44: 43-56. 10.1007/PL00006120.

    CAS  PubMed  Google Scholar 

  58. Popovici C, Roubin R, Coulier F, Birnbaum D: An evolutionary history of the FGF superfamily. Bioessays. 2005, 27: 849-857. 10.1002/bies.20261.

    CAS  PubMed  Google Scholar 

  59. Itoh N, Ornitz DM: Evolution of the Fgf and Fgfr gene families. Trends Genet. 2004, 20: 563-569. 10.1016/j.tig.2004.08.007.

    CAS  PubMed  Google Scholar 

  60. Dempsey LA, Li MJ, DePace A, Bray-Ward P, Maizels N: The human HNRPD locus maps to 4q21 and encodes a highly conserved protein. Genomics. 1998, 49: 378-384. 10.1006/geno.1998.5237.

    CAS  PubMed  Google Scholar 

  61. Chernova OB, Somerville RP, Cowell JK: A novel gene, LGI1, from 10q24 is rearranged and downregulated in malignant brain tumors. Oncogene. 1998, 17: 2873-2881. 10.1038/sj.onc.1202481.

    CAS  PubMed  Google Scholar 

  62. Gu W, Wevers A, Schroder H, Grzeschik KH, Derst C, Brodtkorb E, de Vos R, Steinlein OK: The LGI1 gene involved in lateral temporal lobe epilepsy belongs to a new subfamily of leucine-rich repeat proteins. FEBS letters. 2002, 519: 71-76. 10.1016/S0014-5793(02)02713-8.

    CAS  PubMed  Google Scholar 

  63. Roux PP, Blenis J: ERK and p38 MAPK-activated protein kinases: a family of protein kinases with diverse biological functions. Microbiol Mol Biol Rev. 2004, 68: 320-344. 10.1128/MMBR.68.2.320-344.2004.

    PubMed Central  CAS  PubMed  Google Scholar 

  64. Krens SF, Spaink HP, Snaar-Jagalska BE: Functions of the MAPK family in vertebrate-development. FEBS letters. 2006, 580: 4984-4990. 10.1016/j.febslet.2006.08.025.

    CAS  PubMed  Google Scholar 

  65. Hurlin PJ, Queva C, Koskinen PJ, Steingrimsson E, Ayer DE, Copeland NG, Jenkins NA, Eisenman RN: Mad3 and Mad4: novel Max-interacting transcriptional repressors that suppress c-myc dependent transformation and are expressed during neural and epidermal differentiation. Embo J. 1995, 14: 5646-5659.

    PubMed Central  CAS  PubMed  Google Scholar 

  66. Hurlin PJ, Huang J: The MAX-interacting transcription factor network. Semin Cancer Biol. 2006, 16: 265-274. 10.1016/j.semcancer.2006.07.009.

    CAS  PubMed  Google Scholar 

  67. Toivola DM, Tao GZ, Habtezion A, Liao J, Omary MB: Cellular integrity plus: organelle-related and protein-targeting functions of intermediate filaments. Trends Cell Biol. 2005, 15: 608-617. 10.1016/j.tcb.2005.09.004.

    CAS  PubMed  Google Scholar 

  68. Pennefather JN, Lecci A, Candenas ML, Patak E, Pinto FM, Maggi CA: Tachykinins and tachykinin receptors: a growing family. Life Sci. 2004, 74: 1445-1463. 10.1016/j.lfs.2003.09.039.

    CAS  PubMed  Google Scholar 

  69. Koike K, Urata Y, Goto S: Cloning and nucleotide sequence of the cDNA encoding human 2-oxoglutarate dehydrogenase (lipoamide). Proceedings of the National Academy of Sciences of the United States of America. 1992, 89: 1963-1967. 10.1073/pnas.89.5.1963.

    PubMed Central  CAS  PubMed  Google Scholar 

  70. Szabo P, Cai X, Ali G, Blass JP: Localization of the gene (OGDH) coding for the E1k component of the alpha-ketoglutarate dehydrogenase complex to chromosome 7p13-p11.2. Genomics. 1994, 20: 324-326. 10.1006/geno.1994.1178.

    CAS  PubMed  Google Scholar 

  71. Sadakata T, Furuichi T: Identification and mRNA expression of Ogdh, QP-C, and two predicted genes in the postnatal mouse brain. Neuroscience letters. 2006, 405: 217-222. 10.1016/j.neulet.2006.07.008.

    CAS  PubMed  Google Scholar 

  72. Te Velthuis AJ, Isogai T, Gerrits L, Bagowski CP: Insights into the Molecular Evolution of the PDZ/LIM Family and Identification of a Novel Conserved Protein Motif. PLoS ONE. 2007, 2: e189-10.1371/journal.pone.0000189.

    PubMed Central  PubMed  Google Scholar 

  73. Fox EJ, Stubbs SA, Kyaw Tun J, Leek JP, Markham AF, Wright SC: PRELI (protein of relevant evolutionary and lymphoid interest) is located within an evolutionarily conserved gene cluster on chromosome 5q34-q35 and encodes a novel mitochondrial protein. The Biochemical journal. 2004, 378: 817-825. 10.1042/BJ20031504.

    PubMed Central  CAS  PubMed  Google Scholar 

  74. Huitema K, Dikkenberg van den J, Brouwers JF, Holthuis JC: Identification of a family of animal sphingomyelin synthases. Embo J. 2004, 23: 33-44. 10.1038/sj.emboj.7600034.

    PubMed Central  CAS  PubMed  Google Scholar 

  75. Melkonyan HS, Chang WC, Shapiro JP, Mahadevappa M, Fitzpatrick PA, Kiefer MC, Tomei LD, Umansky SR: SARPs: a family of secreted apoptosis-related proteins. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94: 13636-13641. 10.1073/pnas.94.25.13636.

    PubMed Central  CAS  PubMed  Google Scholar 

  76. Jones SE, Jomary C: Secreted Frizzled-related proteins: searching for relationships and patterns. Bioessays. 2002, 24: 811-820. 10.1002/bies.10136.

    CAS  PubMed  Google Scholar 

  77. Rattner A, Hsieh JC, Smallwood PM, Gilbert DJ, Copeland NG, Jenkins NA, Nathans J: A family of secreted proteins contains homology to the cysteine-rich ligand-binding domain of frizzled receptors. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94: 2859-2863. 10.1073/pnas.94.7.2859.

    PubMed Central  CAS  PubMed  Google Scholar 

  78. Kioka N, Ueda K, Amachi T: Vinexin, CAP/ponsin, ArgBP2: a novel adaptor protein family regulating cytoskeletal organization and signal transduction. Cell Struct Funct. 2002, 27: 1-7. 10.1247/csf.27.1.

    CAS  PubMed  Google Scholar 

  79. Kawakami A, Tian Q, Streuli M, Poe M, Edelhoff S, Disteche CM, Anderson P: Intron-exon organization and chromosomal localization of the human TIA-1 gene. J Immunol. 1994, 152: 4937-4945.

    CAS  PubMed  Google Scholar 

  80. Kawakami A, Tian Q, Duan X, Streuli M, Schlossman SF, Anderson P: Identification and functional characterization of a TIA-1-related nucleolysin. Proceedings of the National Academy of Sciences of the United States of America. 1992, 89: 8681-8685. 10.1073/pnas.89.18.8681.

    PubMed Central  CAS  PubMed  Google Scholar 

  81. Huang S, Yuan S, Dong M, Su J, Yu C, Shen Y, Xie X, Yu Y, Yu X, Chen S, et al: The phylogenetic analysis of tetraspanins projects the evolution of cell-cell interactions from unicellular to multicellular organisms. Genomics. 2005, 86: 674-684. 10.1016/j.ygeno.2005.08.004.

    CAS  PubMed  Google Scholar 

  82. Maecker HT, Todd SC, Levy S: The tetraspanin superfamily: molecular facilitators. Faseb J. 1997, 11: 428-442.

    CAS  PubMed  Google Scholar 

  83. Leonardo ED, Hinck L, Masu M, Keino-Masu K, Ackerman SL, Tessier-Lavigne M: Vertebrate homologues of C. elegans UNC-5 are candidate netrin receptors. Nature. 1997, 386: 833-838. 10.1038/386833a0.

    CAS  PubMed  Google Scholar 

  84. Barallobre MJ, Pascual M, Del Rio JA, Soriano E: The Netrin family of guidance factors: emphasis on Netrin-1 signalling. Brain Res Brain Res Rev. 2005, 49: 22-47. 10.1016/j.brainresrev.2004.11.003.

    CAS  PubMed  Google Scholar 

  85. Beliakoff J, Sun Z: Zimp7 and Zimp10, two novel PIAS-like proteins, function as androgen receptor coregulators. Nucl Recept Signal. 2006, 4: e017-

    PubMed Central  PubMed  Google Scholar 

  86. Huang CY, Beliakoff J, Li X, Lee J, Li X, Sharma M, Lim B, Sun Z: hZimp7, a novel PIAS-like protein, enhances androgen receptor-mediated transcription and interacts with SWI/SNF-like BAF complexes. Molecular endocrinology. 2005, 19: 2915-2929. 10.1210/me.2005-0097.

    CAS  PubMed  Google Scholar 

  87. Sharma M, Li X, Wang Y, Zarnegar M, Huang CY, Palvimo JJ, Lim B, Sun Z: hZimp10 is an androgen receptor co-activator and forms a complex with SUMO-1 at replication foci. Embo J. 2003, 22: 6101-6114. 10.1093/emboj/cdg585.

    PubMed Central  CAS  PubMed  Google Scholar 

  88. Li X, Thyssen G, Beliakoff J, Sun Z: The novel PIAS-like protein hZimp10 enhances Smad transcriptional activity. J Biol Chem. 2006, 281: 23748-23756. 10.1074/jbc.M508365200.

    CAS  PubMed  Google Scholar 

  89. Lagerstrom MC, Fredriksson R, Bjarnadottir TK, Fridmanis D, Holmquist T, Andersson J, Yan YL, Raudsepp T, Zoorob R, Kukkonen JP, Lundin LG, Klovins J, Chowdhary BP, Postlethwait JH, Schioth HB: Origin of the prolactin-releasing hormone (PRLH) receptors: evidence of coevolution between PRLH and a redundant neuropeptide Y receptor during vertebrate evolution. Genomics. 2005, 85: 688-703. 10.1016/j.ygeno.2005.02.007.

    PubMed  Google Scholar 

  90. Klovins J, Haitina T, Fridmanis D, Kilianova Z, Kapa I, Fredriksson R, Gallo-Payet N, Schioth HB: The melanocortin system in Fugu: determination of POMC/AGRP/MCR gene repertoire and synteny, as well as pharmacology and anatomical distribution of the MCRs. Molecular biology and evolution. 2004, 21: 563-579. 10.1093/molbev/msh050.

    CAS  PubMed  Google Scholar 

  91. Aparicio S, Chapman J, Stupka E, Putnam N, Chia JM, Dehal P, Christoffels A, Rash S, Hoon S, Smit A, Gelpke MD, Roach J, Oh T, Ho IY, Wong M, Detter C, Verhoef F, Predki P, Tay A, Lucas S, Richardson P, Smith SF, Clark MS, Edwards YJ, Doggett N, Zharkikh A, Tavtigian SV, Pruss D, Barnstead M, Evans C, et al: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science. 2002, 297: 1301-1310. 10.1126/science.1072104.

    CAS  PubMed  Google Scholar 

  92. Sundström G, Larsson TA, Brenner S, Venkatesh B, Larhammar D: Evolution of the neuropeptide Y family: New genes by chromosome duplications in early vertebrates and in teleost fishes. Gen Comp Endocrinol. 2008, 155: 705-716. 10.1016/j.ygcen.2007.08.016.

    PubMed  Google Scholar 

  93. Skrabanek L, Wolfe KH: Eukaryote genome duplication – where's the evidence?. Curr Opin Genet Dev. 1998, 8: 694-700. 10.1016/S0959-437X(98)80039-7.

    CAS  PubMed  Google Scholar 

  94. Peer Van de Y: Computational approaches to unveiling ancient genome duplications. Nat Rev Genet. 2004, 5: 752-763. 10.1038/nrg1449.

    PubMed  Google Scholar 

  95. Durand D, Hoberman R: Diagnosing duplications – can it be done?. Trends Genet. 2006, 22: 156-164. 10.1016/j.tig.2006.01.002.

    CAS  PubMed  Google Scholar 

  96. Simillion C, Vandepoele K, Peer Van de Y: Recent developments in computational approaches for uncovering genomic homology. Bioessays. 2004, 26: 1225-1235. 10.1002/bies.20127.

    CAS  PubMed  Google Scholar 

  97. Martin N, Ruedi EA, Leduc R, Sun FJ, Caetano-Anolles G: Gene-interleaving patterns of synteny in the Saccharomyces cerevisiae genome: are they proof of an ancient genome duplication event?. Biology direct. 2007, 2: 23-10.1186/1745-6150-2-23.

    PubMed Central  PubMed  Google Scholar 

  98. Ruuskanen JO, Xhaard H, Marjamaki A, Salaneck E, Salminen T, Yan YL, Postlethwait JH, Johnson MS, Larhammar D, Scheinin M: Identification of duplicated fourth alpha2-adrenergic receptor subtype by cloning and mapping of five receptor genes in zebrafish. Molecular biology and evolution. 2004, 21: 14-28. 10.1093/molbev/msg224.

    CAS  PubMed  Google Scholar 

  99. Blomme T, Vandepoele K, De Bodt S, Simillion C, Maere S, Peer Van de Y: The gain and loss of genes during 600 million years of vertebrate evolution. Genome biology. 2006, 7: R43-10.1186/gb-2006-7-5-r43.

    PubMed Central  PubMed  Google Scholar 

  100. Hoegg S, Meyer A: Hox clusters as models for vertebrate genome evolution. Trends Genet. 2005, 21: 421-424. 10.1016/j.tig.2005.06.004.

    CAS  PubMed  Google Scholar 

  101. The BioGRID. [http://www.thebiogrid.org]

  102. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic acids research. 2006, 34: D535-539. 10.1093/nar/gkj109.

    PubMed Central  CAS  PubMed  Google Scholar 

  103. Shimeld SM: Gene function, gene networks and the fate of duplicated genes. Semin Cell Dev Biol. 1999, 10: 549-553. 10.1006/scdb.1999.0336.

    CAS  PubMed  Google Scholar 

  104. Conant GC, Wolfe KH: Functional partitioning of yeast co-expression networks after genome duplication. PLoS biology. 2006, 4: e109-10.1371/journal.pbio.0040109.

    PubMed Central  PubMed  Google Scholar 

  105. Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004, 16: 1679-1691. 10.1105/tpc.021410.

    PubMed Central  CAS  PubMed  Google Scholar 

  106. Gachon CM, Langlois-Meurinne M, Henry Y, Saindrenan P: Transcriptional co-regulation of secondary metabolism enzymes in Arabidopsis: functional and evolutionary implications. Plant Mol Biol. 2005, 58: 229-245. 10.1007/s11103-005-5346-5.

    CAS  PubMed  Google Scholar 

  107. 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 Central  CAS  PubMed  Google Scholar 

  108. Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290 (5494): 1151-1155. 10.1126/science.290.5494.1151.

    CAS  PubMed  Google Scholar 

  109. Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization. Genetics. 2000, 154: 459-473.

    PubMed Central  CAS  PubMed  Google Scholar 

  110. Lynch M, O'Hely M, Walsh B, Force A: The probability of preservation of a newly arisen gene duplicate. Genetics. 2001, 159: 1789-1804.

    PubMed Central  CAS  PubMed  Google Scholar 

  111. He X, Zhang J: Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics. 2005, 169: 1157-1164. 10.1534/genetics.104.037051.

    PubMed Central  PubMed  Google Scholar 

  112. Rastogi S, Liberles DA: Subfunctionalization of duplicated genes as a transition state to neofunctionalization. BMC Evol Biol. 2005, 5: 28-10.1186/1471-2148-5-28.

    PubMed Central  PubMed  Google Scholar 

  113. Benton MJ, Donoghue PC: Paleontological evidence to date the tree of life. Molecular biology and evolution. 2007, 24: 26-53. 10.1093/molbev/msl150.

    CAS  PubMed  Google Scholar 

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

  115. Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ: Multiple sequence alignment with Clustal X. Trends Biochem Sci. 1998, 23: 403-405. 10.1016/S0968-0004(98)01285-7.

    CAS  PubMed  Google Scholar 

  116. 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 research. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.

    PubMed Central  CAS  PubMed  Google Scholar 

  117. Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.

    CAS  PubMed  Google Scholar 

  118. Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.

    CAS  PubMed  Google Scholar 

  119. The Pfam database. [http://pfam.sanger.ac.uk]

Download references

Acknowledgements

The authors would like to thank Christina Bergqvist for excellent technical assistance and Susanne Dreborg for help with preparing illustrations. This work was supported by grants from the Swedish Research Council and Carl Trygger's Foundation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Dan Larhammar.

Additional information

Authors' contributions

TAL performed the phylogenetic and chromosome analyses and drafted and coordinated the manuscript, FO performed many of the initial analyses, GS participated in the phylogenetic and chromosome analyses and contributed to the manuscript. L–GL spawned the concept of duplication of these chromosomal regions and contributed conceptually. SB participated in the design of experiments and contributed conceptually. BV performed the experimental work in T. rubripes and wrote part of the manuscript. DL conceived and initiated the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Larsson, T.A., Olsson, F., Sundstrom, G. et al. Early vertebrate chromosome duplications and the evolution of the neuropeptide Y receptor gene regions. BMC Evol Biol 8, 184 (2008). https://doi.org/10.1186/1471-2148-8-184

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-8-184

Keywords