MetaHox gene clusters

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Evolution of vertebrate nicotinic acetylcholine receptors

Background

Many physiological processes are influenced by nicotinic acetylcholine receptors (nAChR), ranging from neuromuscular and parasympathetic signaling to modulation of the reward system and long-term memory. Due to the complexity of the nAChR family and variable evolutionary rates among its members, their evolution in vertebrates has been difficult to resolve. In order to understand how and when the nAChR genes arose, we have used a broad approach of analyses combining sequence-based phylogeny, chromosomal synteny and intron positions.

Results

Our analyses suggest that there were ten subunit genes present in the vertebrate predecessor. The two basal vertebrate tetraploidizations (1R and 2R) then expanded this set to 19 genes. Three of these have been lost in mammals, resulting in 16 members today. None of the ten ancestral genes have kept all four copies after 2R. Following 2R, two of the ancestral genes became triplicates, five of them became pairs, and three seem to have remained single genes. One triplet consists of CHRNA7, CHRNA8 and the previously undescribed CHRNA11, of which the two latter have been lost in mammals but are still present in lizards and ray-finned fishes. The other triplet consists of CHRNB2, CHRNB4 and CHRNB5, the latter of which has also been lost in mammals. In ray-finned fish the neuromuscular subunit gene CHRNB1 underwent a local gene duplication generating CHRNB1.2. The third tetraploidization in the predecessor of teleosts (3R) expanded the repertoire to a total of 31 genes, of which 27 remain in zebrafish. These evolutionary relationships are supported by the exon-intron organization of the genes.

Conclusion

The tetraploidizations explain all gene duplication events in vertebrates except two. This indicates that the genome doublings have had a substantial impact on the complexity of this gene family leading to a very large number of members that have existed for hundreds of millions of years.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1341-8) contains supplementary material, which is available to authorized users.

Evolution of the Muscarinic Acetylcholine Receptors in Vertebrates

The family of muscarinic acetylcholine receptors (mAChRs) consists of five members in mammals, encoded by the CHRM1-5 genes. The mAChRs are G-protein-coupled receptors, which can be divided into the following two subfamilies: M2 and M4 receptors coupling to G_i/o; and M1, M3, and M5 receptors coupling to G_q/11. However, despite the fundamental roles played by these receptors, their evolution in vertebrates has not yet been fully described. We have combined sequence-based phylogenetic analyses with comparisons of exon–intron organization and conserved synteny in order to deduce the evolution of the mAChR receptors. Our analyses verify the existence of two ancestral genes prior to the two vertebrate tetraploidizations (1R and 2R). After these events, one gene had duplicated, resulting in CHRM2 and CHRM4; and the other had triplicated, forming the CHRM1, CHRM3, and CHRM5 subfamily. All five genes are still present in all vertebrate groups investigated except the CHRM1 gene, which has not been identified in some of the teleosts or in chicken or any other birds. Interestingly, the third tetraploidization (3R) that took place in the teleost predecessor resulted in duplicates of all five mAChR genes of which all 10 are present in zebrafish. One of the copies of the CHRM2 and CHRM3 genes and both CHRM4 copies have gained introns in teleosts. Not a single separate (nontetraploidization) duplicate has been identified in any vertebrate species. These results clarify the evolution of the vertebrate mAChR family and reveal a doubled repertoire in zebrafish, inviting studies of gene neofunctionalization and subfunctionalization.

Evolution of neuropeptide signalling systems

Neuropeptides are a diverse class of neuronal signalling molecules that regulate physiological processes and behaviour in animals. However, determining the relationships and evolutionary origins of the heterogeneous assemblage of neuropeptides identified in a range of phyla has presented a huge challenge for comparative physiologists. Here, we review revolutionary insights into the evolution of neuropeptide signalling that have been obtained recently through comparative analysis of genome/transcriptome sequence data and by ‘deorphanisation’ of neuropeptide receptors. The evolutionary origins of at least 30 neuropeptide signalling systems have been traced to the common ancestor of protostomes and deuterostomes. Furthermore, two rounds of genome duplication gave rise to an expanded repertoire of neuropeptide signalling systems in the vertebrate lineage, enabling neofunctionalisation and/or subfunctionalisation, but with lineage-specific gene loss and/or additional gene or genome duplications generating complex patterns in the phylogenetic distribution of paralogous neuropeptide signalling systems. We are entering a new era in neuropeptide research where it has become feasible to compare the physiological roles of orthologous and paralogous neuropeptides in a wide range of phyla. Moreover, the ambitious mission to reconstruct the evolution of neuropeptide function in the animal kingdom now represents a tangible challenge for the future.

Summary: A review of the revolutionary advances in our knowledge of the evolution of neuropeptide signalling systems that have been enabled by comparative genomics and neuropeptide receptor deorphanisation.

Expansion of secretin-like G protein-coupled receptors and their peptide ligands via local duplications before and after two rounds of whole-genome duplication

In humans, the secretin-like G protein-coupled receptor (GPCR) family comprises 15 members with 18 corresponding peptide ligand genes. Although members have been identified in a large variety of vertebrate and nonvertebrate species, the origin and relationship of these proteins remain unresolved. To address this issue, we employed large-scale genome comparisons to identify genome fragments with conserved synteny and matched these fragments to linkage groups in reconstructed early gnathostome ancestral chromosomes (GAC). This genome comparison revealed that most receptor and peptide genes were clustered in three GAC linkage groups and suggested that the ancestral forms of five peptide subfamilies (corticotropin-releasing hormone-like, calcitonin-like, parathyroid hormone-like, glucagon-like, and growth hormone-releasing hormone-like) and their cognate receptor families emerged through tandem local gene duplications before two rounds (2R) of whole-genome duplication. These subfamily genes have, then, been amplified by 2R whole-genome duplication, followed by additional local duplications and gene loss prior to the divergence of land vertebrates and teleosts. This study delineates a possible evolutionary scenario for whole secretin-like peptide and receptor family members and may shed light on evolutionary mechanisms for expansion of a gene family with a large number of paralogs.

The cytochrome P450 genesis locus: the origin and evolution of animal cytochrome P450s

The neighbourhoods of cytochrome P450 (CYP) genes in deuterostome genomes, as well as those of the cnidarians Nematostella vectensis and Acropora digitifera and the placozoan Trichoplax adhaerens were examined to find clues concerning the evolution of CYP genes in animals. CYP genes created by the 2R whole genome duplications in chordates have been identified. Both microsynteny and macrosynteny were used to identify genes that coexisted near CYP genes in the animal ancestor. We show that all 11 CYP clans began in a common gene environment. The evidence implies the existence of a single locus, which we term the 'cytochrome P450 genesis locus', where one progenitor CYP gene duplicated to create a tandem set of genes that were precursors of the 11 animal CYP clans: CYP Clans 2, 3, 4, 7, 19, 20, 26, 46, 51, 74 and mitochondrial. These early CYP genes existed side by side before the origin of cnidarians, possibly with a few additional genes interspersed. The Hox gene cluster, WNT genes, an NK gene cluster and at least one ARF gene were close neighbours to this original CYP locus. According to this evolutionary scenario, the CYP74 clan originated from animals and not from land plants nor from a common ancestor of plants and animals. The CYP7 and CYP19 families that are chordate-specific belong to CYP clans that seem to have originated in the CYP genesis locus as well, even though this requires many gene losses to explain their current distribution. The approach to uncovering the CYP genesis locus overcomes confounding effects because of gene conversion, sequence divergence, gene birth and death, and opens the way to understanding the biodiversity of CYP genes, families and subfamilies, which in animals has been obscured by more than 600 Myr of evolution.

Expansion of transducin subunit gene families in early vertebrate tetraploidizations

Hundreds of gene families expanded in the early vertebrate tetraploidizations including many gene families in the phototransduction cascade. We have investigated the evolution of the heterotrimeric G-proteins of photoreceptors, the transducins, in relation to these events using both phylogenetic analyses and synteny comparisons. Three alpha subunit genes were identified in amniotes and the coelacanth, GNAT1-3; two of these were identified in amphibians and teleost fish, GNAT1 and GNAT2. Most tetrapods have four beta genes, GNB1-4, and teleosts have additional duplicates. Finally, three gamma genes were identified in mammals, GNGT1, GNG11 and GNGT2. Of these, GNGT1 and GNGT2 were found in the other vertebrates. In frog and zebrafish additional duplicates of GNGT2 were identified. Our analyses show all three transducin families expanded during the early vertebrate tetraploidizations and the beta and gamma families gained additional copies in the teleost-specific genome duplication. This suggests that the tetraploidizations contributed to visual specialisations.

Whole-genome duplications spurred the functional diversification of the globin gene superfamily in vertebrates

It has been hypothesized that two successive rounds of whole-genome duplication (WGD) in the stem lineage of vertebrates provided genetic raw materials for the evolutionary innovation of many vertebrate-specific features. However, it has seldom been possible to trace such innovations to specific functional differences between paralogous gene products that derive from a WGD event. Here, we report genomic evidence for a direct link between WGD and key physiological innovations in the vertebrate oxygen transport system. Specifically, we demonstrate that key globin proteins that evolved specialized functions in different aspects of oxidative metabolism (hemoglobin, myoglobin, and cytoglobin) represent paralogous products of two WGD events in the vertebrate common ancestor. Analysis of conserved macrosynteny between the genomes of vertebrates and amphioxus (subphylum Cephalochordata) revealed that homologous chromosomal segments defined by myoglobin + globin-E, cytoglobin, and the α-globin gene cluster each descend from the same linkage group in the reconstructed proto-karyotype of the chordate common ancestor. The physiological division of labor between the oxygen transport function of hemoglobin and the oxygen storage function of myoglobin played a pivotal role in the evolution of aerobic energy metabolism, supporting the hypothesis that WGDs helped fuel key innovations in vertebrate evolution.

Extensive chordate and annelid macrosynteny reveals ancestral homeobox gene organization

Genes with the homeobox motif are crucial in developmental biology and widely implicated in the evolution of development. The Antennapedia (ANTP)-class is one of the two major classes of animal homeobox genes, and includes the Hox genes, renowned for their role in patterning the anterior-posterior axis of animals. The origin and evolution of the ANTP-class genes are a matter of some debate. A principal guiding hypothesis has been the existence of an ancient gene Mega-cluster deep in animal ancestry. This hypothesis was largely established from linkage data from chordates, and the Mega-cluster hypothesis remains to be seriously tested in protostomes. We have thus mapped ANTP-class homeobox genes to the chromosome level in a lophotrochozoan protostome. Our comparison of gene organization in Platynereis dumerilii and chordates indicates that the Mega-cluster, if it did exist, had already been broken up onto four chromosomes by the time of the protostome-deuterostome ancestor (PDA). These results not only elucidate an aspect of the genome organization of the PDA but also reveal high levels of macrosynteny between P. dumerilii and chordates. This implies a very low rate of interchromosomal genome rearrangement in the lineages leading to P. dumerilii and the chordate ancestor since the time of the PDA.

A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia

Background

Homeobox genes are a superclass of transcription factors with diverse developmental regulatory functions, which are found in plants, fungi and animals. In animals, several Antennapedia (ANTP)-class homeobox genes reside in extremely ancient gene clusters (for example, the Hox, ParaHox, and NKL clusters) and the evolution of these clusters has been implicated in the morphological diversification of animal bodyplans. By contrast, similarly ancient gene clusters have not been reported among the other classes of homeobox genes (that is, the LIM, POU, PRD and SIX classes).

Results

Using a combination of in silico queries and phylogenetic analyses, we found that a cluster of three PRD-class homeobox genes (Homeobrain (hbn), Rax (rx) and Orthopedia (otp)) is present in cnidarians, insects and mollusks (a partial cluster comprising hbn and rx is present in the placozoan Trichoplax adhaerens). We failed to identify this 'HRO' cluster in deuterostomes; in fact, the Homeobrain gene appears to be missing from the chordate genomes we examined, although it is present in hemichordates and echinoderms. To illuminate the ancestral organization and function of this ancient cluster, we mapped the constituent genes against the assembled genome of a model cnidarian, the sea anemone Nematostella vectensis, and characterized their spatiotemporal expression using in situ hybridization. In N. vectensis, these genes reside in a span of 33 kb with the same gene order as previously reported in insects. Comparisons of genomic sequences and expressed sequence tags revealed the presence of alternative transcripts of Nv-otp and two highly unusual protein-coding polymorphisms in the terminal helix of the Nv-rx homeodomain. A population genetic survey revealed the Rx polymorphisms to be widespread in natural populations. During larval development, all three genes are expressed in the ectoderm, in non-overlapping territories along the oral-aboral axis, with distinct temporal expression.

Conclusion

We report the first evidence for a PRD-class homeobox cluster that appears to have been conserved since the time of the cnidarian-bilaterian ancestor, and possibly even earlier, given the presence of a partial cluster in the placozoan Trichoplax. Very similar clusters comprising these three genes exist in Nematostella and diverse protostomes. Interestingly, in chordates, one member of the ancestral cluster (homeobrain) has apparently been lost, and there is no linkage between rx and orthopedia in any of the vertebrates. In Nematostella, the spatial expression of these three genes along the body column is not colinear with their physical order in the cluster but the temporal expression is, therefore, using the terminology that has been applied to the Hox cluster genes, the HRO cluster would appear to exhibit temporal but not spatial colinearity. It remains to be seen whether the mechanisms responsible for the evolutionary conservation of the HRO cluster are the same mechanisms responsible for cohesion of the Hox cluster and other ANTP-class homeobox clusters that have been widely conserved throughout animal evolution.

Evolution of vertebrate rod and cone phototransduction genes

Vertebrate cones and rods in several cases use separate but related components for their signal transduction (opsins, G-proteins, ion channels, etc.). Some of these proteins are also used differentially in other cell types in the retina. Because cones, rods and other retinal cell types originated in early vertebrate evolution, it is of interest to see if their specific genes arose in the extensive gene duplications that took place in the ancestor of the jawed vertebrates (gnathostomes) by two tetraploidizations (genome doublings). The ancestor of teleost fishes subsequently underwent a third tetraploidization. Our previously reported analyses showed that several gene families in the vertebrate visual phototransduction cascade received new members in the basal tetraploidizations. We here expand these data with studies of additional gene families and vertebrate species. We conclude that no less than 10 of the 13 studied phototransduction gene families received additional members in the two basal vertebrate tetraploidizations. Also the remaining three families seem to have undergone duplications during the same time period but it is unclear if this happened as a result of the tetraploidizations. The implications of the many early vertebrate gene duplications for functional specialization of specific retinal cell types, particularly cones and rods, are discussed.

TFCat: the curated catalog of mouse and human transcription factors

TFCat is a catalog of mouse and human transcription factors based on a reliable core collection of annotations obtained by expert review of the scientific literature

Unravelling regulatory programs governed by transcription factors (TFs) is fundamental to understanding biological systems. TFCat is a catalog of mouse and human TFs based on a reliable core collection of annotations obtained by expert review of the scientific literature. The collection, including proven and homology-based candidate TFs, is annotated within a function-based taxonomy and DNA-binding proteins are organized within a classification system. All data and user-feedback mechanisms are available at the TFCat portal .

Domain duplication, divergence, and loss events in vertebrate Msx paralogs reveal phylogenomically informed disease markers

Background

Msx originated early in animal evolution and is implicated in human genetic disorders. To reconstruct the functional evolution of Msx and inform the study of human mutations, we analyzed the phylogeny and synteny of 46 metazoan Msx proteins and tracked the duplication, diversification and loss of conserved motifs.

Results

Vertebrate Msx sequences sort into distinct Msx1, Msx2 and Msx3 clades. The sister-group relationship between MSX1 and MSX2 reflects their derivation from the 4p/5q chromosomal paralogon, a derivative of the original "MetaHox" cluster. We demonstrate physical linkage between Msx and other MetaHox genes (Hmx, NK1, Emx) in a cnidarian. Seven conserved domains, including two Groucho repression domains (N- and C-terminal), were present in the ancestral Msx. In cnidarians, the Groucho domains are highly similar. In vertebrate Msx1, the N-terminal Groucho domain is conserved, while the C-terminal domain diverged substantially, implying a novel function. In vertebrate Msx2 and Msx3, the C-terminal domain was lost. MSX1 mutations associated with ectodermal dysplasia or orofacial clefting disorders map to conserved domains in a non-random fashion.

Conclusion

Msx originated from a MetaHox ancestor that also gave rise to Tlx, Demox, NK, and possibly EHGbox, Hox and ParaHox genes. Duplication, divergence or loss of domains played a central role in the functional evolution of Msx. Duplicated domains allow pleiotropically expressed proteins to evolve new functions without disrupting existing interaction networks. Human missense sequence variants reside within evolutionarily conserved domains, likely disrupting protein function. This phylogenomic evaluation of candidate disease markers will inform clinical and functional studies.

Evolution of vertebrate opioid receptors

The opioid peptides and receptors have prominent roles in pain transmission and reward mechanisms in mammals. The evolution of the opioid receptors has so far been little studied, with only a few reports on species other than tetrapods. We have investigated species representing a broader range of vertebrates and found that the four opioid receptor types (delta, kappa, mu, and NOP) are present in most of the species. The gene relationships were deduced by using both phylogenetic analyses and chromosomal location relative to 20 neighboring gene families in databases of assembled genomes. The combined results show that the vertebrate opioid receptor gene family arose by quadruplication of a large chromosomal block containing at least 14 other gene families. The quadruplication seems to coincide with, and, therefore, probably resulted from, the two proposed genome duplications in early vertebrate evolution. We conclude that the quartet of opioid receptors was already present at the origin of jawed vertebrates approximately 450 million years ago. A few additional opioid receptor gene duplications have occurred in bony fishes. Interestingly, the ancestral receptor gene duplications coincide with the origin of the four opioid peptide precursor genes. Thus, the complete vertebrate opioid system was already established in the first jawed vertebrates.

Phylogenetic and chromosomal analyses of multiple gene families syntenic with vertebrate Hox clusters

Background

Ever since the theory about two rounds of genome duplication (2R) in the vertebrate lineage was proposed, the Hox gene clusters have served as the prime example of quadruplicate paralogy in mammalian genomes. In teleost fishes, the observation of additional Hox clusters absent in other vertebrate lineages suggested a third tetraploidization (3R). Because the Hox clusters occupy a quite limited part of each chromosome, and are special in having position-dependent regulation within the multi-gene cluster, studies of syntenic gene families are needed to determine the extent of the duplicated chromosome segments. We have analyzed in detail 14 gene families that are syntenic with the Hox clusters to see if their phylogenies are compatible with the Hox duplications and the 2R/3R scenario. Our starting point was the gene family for the NPY family of peptides located near the Hox clusters in the pufferfish Takifugu rubripes, the zebrafish Danio rerio, and human.

Results

Seven of the gene families have members on at least three of the human Hox chromosomes and two families are present on all four. Using both neighbor-joining and quartet-puzzling maximum likelihood methods we found that 13 families have a phylogeny that supports duplications coinciding with the Hox cluster duplications. One additional family also has a topology consistent with 2R but due to lack of urochordate or cephalocordate sequences the time window when these duplications could have occurred is wider. All but two gene families also show teleost-specific duplicates.

Conclusion

Based on this analysis we conclude that the Hox cluster duplications involved a large number of adjacent gene families, supporting expansion of these families in the 2R, as well as in the teleost 3R tetraploidization. The gene duplicates presumably provided raw material in early vertebrate evolution for neofunctionalization and subfunctionalization.

The early ANTP gene repertoire: insights from the placozoan genome

The evolution of ANTP genes in the Metazoa has been the subject of conflicting hypotheses derived from full or partial gene sequences and genomic organization in higher animals. Whole genome sequences have recently filled in some crucial gaps for the basal metazoan phyla Cnidaria and Porifera. Here we analyze the complete genome of Trichoplax adhaerens, representing the basal metazoan phylum Placozoa, for its set of ANTP class genes. The Trichoplax genome encodes representatives of Hox/ParaHox-like, NKL, and extended Hox genes. This repertoire possibly mirrors the condition of a hypothetical cnidarian-bilaterian ancestor. The evolution of the cnidarian and bilaterian ANTP gene repertoires can be deduced by a limited number of cis-duplications of NKL and “extended Hox” genes and the presence of a single ancestral “ProtoHox” gene.

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

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.

Inferring ancestral gene order

To explain the evolutionary mechanisms by which populations of organisms change over time, it is necessary to first understand the pathways by which genomes have changed over time. Understanding genome evolution requires comparing modern genomes with ancestral genomes, which thus necessitates the reconstruction of those ancestral genomes. This chapter describes automated approaches to infer the nature of ancestral genomes from modern sequenced genomes. Because several rounds of whole genome duplication have punctuated the evolution of animals with backbones, and current methods for ortholog calling do not adequately account for such events, we developed ways to infer the nature of ancestral chromosomes after genome duplication. We apply this method here to reconstruct the ancestors of a specific chromosome in the zebrafish Danio rerio.

The zebrafish genome in context: ohnologs gone missing

Some zebrafish genes appear to lack an ortholog in the human genome and researchers often call them "novel" genes. The origin of many so-called "novel" genes becomes apparent when considered in the context of genome duplication events that occurred during evolution of the phylum Chordata, including two rounds at about the origin of the subphylum Vertebrata (R1 and R2) and one round before the teleost radiation (R3). Ohnologs are paralogs stemming from such genome duplication events, and some zebrafish genes said to be "novel" are more appropriately interpreted as "ohnologs gone missing", cases in which ohnologs are preserved differentially in different evolutionary lineages. Here we consider ohnologs present in the zebrafish genome but absent from the human genome. Reasonable hypotheses are that lineage-specific loss of ohnologs can play a role in establishing lineage divergence and in the origin of developmental innovations. How does the evolution of ohnologs differ from the evolution of gene duplicates arising from other mechanisms, such as tandem duplication or retrotransposition? To what extent do different major vertebrate lineages or different teleost lineages differ in ohnolog content? What roles do differences in ohnolog content play in the origin of developmental mechanisms that differ among lineages? This review explores these questions.

The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone, Nematostella vectensis

BACKGROUND

Homeodomain transcription factors are key components in the developmental toolkits of animals. While this gene superclass predates the evolutionary split between animals, plants, and fungi, many homeobox genes appear unique to animals. The origin of particular homeobox genes may, therefore, be associated with the evolution of particular animal traits. Here we report the first near-complete set of homeodomains from a basal (diploblastic) animal.

RESULTS

Phylogenetic analyses were performed on 130 homeodomains from the sequenced genome of the sea anemone Nematostella vectensis along with 228 homeodomains from human and 97 homeodomains from Drosophila. The Nematostella homeodomains appear to be distributed among established homeodomain classes in the following fashion: 72 ANTP class; one HNF class; four LIM class; five POU class; 33 PRD class; five SINE class; and six TALE class. For four of the Nematostella homeodomains, there is disagreement between neighbor-joining and Bayesian trees regarding their class membership. A putative Nematostella CUT class gene is also identified.

CONCLUSION

The homeodomain superclass underwent extensive radiations prior to the evolutionary split between Cnidaria and Bilateria. Fifty-six homeodomain families found in human and/or fruit fly are also found in Nematostella, though seventeen families shared by human and fly appear absent in Nematostella. Homeodomain loss is also apparent in the bilaterian taxa: eight homeodomain families shared by Drosophila and Nematostella appear absent from human (CG13424, EMXLX, HOMEOBRAIN, MSXLX, NK7, REPO, ROUGH, and UNC4), and six homeodomain families shared by human and Nematostella appear absent from fruit fly (ALX, DMBX, DUX, HNF, POU1, and VAX).

The Mitochondrial Subgenomes of the Nematode Globodera pallida Are Mosaics: Evidence of Recombination in an Animal Mitochondrial Genome

We sequenced four mitochondrial subgenomes from the potato cyst nematode Globodera pallida, previously characterized as one of the few animals to have a multipartite mitochondrial genome. The sequence data indicate that three of these subgenomic mitochondrial circles are mosaics, comprising long, multigenic fragments derived from fragments of the other circles. This pattern is consistent with the operation of intermitochondrial recombination, a process generally considered absent in animal mitochondria. We also report that many of the duplicated genes contain deleterious mutations, ones likely to render the gene nonfunctional; gene conversion does not appear to be homogenizing the different gene copies. The proposed nonfunctional copies are clustered on particular circles, whereas copies that are likely to code functional gene products are clustered on others.

Evidence for Evolving Toll-IL-1 Receptor-Containing Adaptor Molecule Function in Vertebrates

In mammals, Toll-IL-1R-containing adaptor molecule 1 (TICAM1)-dependent TLR pathways induce NF-kappaB and IFN-beta responses. TICAM1 activates NF-kappaB through two different pathways involving its interactions with TNFR-associated factor 6 and receptor-interacting protein 1. It also activates IFN regulatory factor 3/7 through its interaction with TANK-binding kinase-1, leading to the robust up-regulation of IFN-beta. In this study, we describe the role of zebrafish (Danio rerio) TICAM1 in activating NF-kappaB and zebrafish type I IFN. Zebrafish IFN is unique in that it cannot be categorized as being alpha- or beta-like. Through comprehensive sequence, phylogenetic, and syntenic analyses, we fully describe the identification of a zebrafish TICAM1 ortholog. Zebrafish TICAM1 exhibits sequence divergence from its mammalian orthologs and our data demonstrate that these sequence differences have functional consequences. Zebrafish TICAM1 activates zebrafish IFN; however, it does so in an apparently IFN regulatory factor 3/7-independent manner. Furthermore, zebrafish TICAM1 does not interact with zebrafish TNFR-associated factor 6, thus NF-kappaB activation is dependent upon its interaction with receptor-interacting protein 1. Comparative genome analysis suggests that TICAM1 and TICAM2 evolved from a common vertebrate TICAM ancestor following a gene duplication event and that TICAM2 was lost in teleosts following the divergence of the rayfin and lobefin fishes 450 million years ago. These studies provide evidence, for the first time, of the evolving function of a vertebrate TLR pathway.

Ancient complexity of the non-Hox ANTP gene complement in the anthozoan Nematostella vectensis: implications for the evolution of the ANTP superclass

The origin and evolution of ANTP superclass genes has raised controversial discussions. While recent evidence suggests that a true Hox cluster emerged after the cnidarian bilaterian split, the origin of the ANTP superclass as a whole remains unclear. Based on analyses of bilaterian genomes, it seems very likely that clustering has once been a characteristic of all ANTP homeobox genes and that their ancestors have emerged through several series of cis-duplications from the same genomic region. Since the diploblastic Cnidaria possess orthologs of some non-Hox ANTP genes, at least some steps of the expansion of this hypothetical homeobox gene array must have occurred in the last common ancestor of both lineages--but it is unknown to what extent. By screening the unassembled Nematostella genome, we have identified unambiguous orthologs to almost all non-Hox ANTP genes which are present in Bilateria--with the exception of En, Tlx and (possibly) Vax. Furthermore, Nematostella possesses ANTP genes that are missing in some bilaterian lineages, like the rough gene or NK7. In addition, several ANTP homeobox gene families have been independently duplicated in Nematostella. We conclude that the last cnidarian/bilaterian ancestor already harboured the almost full complement of non-Hox ANTP genes before the Hox system evolved.

Ancient linkage of a POU class 6 and an anterior hox-like gene in cnidaria: implications for the evolution of homeobox genes

Linkage analyses in metazoan genomes suggest two ancestral arrays for the majority of homeobox genes. The related homeobox genes and chromosomal regions that are dispersed in extant species derived possibly from only two single common ancestor regions. One proposed ancestral array, designated as ANTP mega-array, contains most of the ANTP class homeobox genes; the second, named the contraHox super-paralogon, would consist of the classes PRD, POU, LIM, CUT, prospero, TALE and SIX. Here, we report the tight linkage of a POU class 6 gene to an anterior Hox-like gene in the hydrozoan Eleutheria dichotoma and discuss its possible significance for the evolution of homeobox genes. POU class 6 genes also seem to be ancestrally linked to the HoxC and A clusters in vertebrates, despite POU homeobox genes belonging to the contraHox paralogon. Hence, the much tighter linkage of a POU class 6 gene to an anterior Hox-like gene in a cnidarian is possibly the evolutionary echo of an ancestral genomic region from which most metazoan homeobox classes emerged.

Formation and evolution of the chordate neurotrophin and Trk receptor genes

Neurotrophins are structurally related neurotrophic polypeptide factors that regulate neuronal differentiation and are essential for neuronal survival, neurite growth and plasticity. It has until very recently been thought that the neurotrophin system appeared with the vertebrate species, but identification of a cephalochordate neurotrophin receptor (Trk), and more recently neurotrophin sequences in several genomes of deuterostome invertebrates, show that the system already existed at the stem of the deuterostome group. Comparative genomics supports the hypothesis that two whole genome duplications produced many of the vertebrate gene families, among those the neurotrophin and Trk families. It remains to be proven to what extent the whole genome duplications have driven macroevolutionary change, but it appears certain that the formation of the multi-gene copy neurotrophin and Trk receptor families at the stem of vertebrates has provided a foundation from which the various functions and pleiotropic effects produced by each of the four extant neurotrophins have evolved.

Genetic organization and embryonic expression of the ParaHox genes in the sea urchin S. purpuratus: insights into the relationship between clustering and colinearity

The ANTP family of homeodomain transcription factors consists of three major groups, the NKL, the extended Hox, and the Hox/ParaHox family. Hox genes and ParaHox genes are often linked in the genome forming two clusters of genes, the Hox cluster and the ParaHox cluster, and are expressed along the major body axis in a nested fashion, following the relative positions of the genes within these clusters, a property called colinearity. While the presences of a Hox cluster and a ParaHox cluster appear to be primitive for bilaterians, few taxa have actually been examined for spatial and temporal colinearity, and, aside from chordates, even fewer still manifest it. Here we show that the ParaHox genes of the sea urchin Strongylocentrotus purpuratus show both spatial and temporal colinearity, but with peculiarities. Specifically, two of the three ParaHox genes-discovered through the S. purpuratus genome project-Sp-lox and Sp-Cdx, are expressed in the developing gut with nested domains in a spatially colinear manner. However, transcripts of Sp-Gsx, although anterior of Sp-lox, are detected in the ectoderm and not in the gut. Strikingly, the expression of the three ParaHox genes would follow temporal colinearity if they were clustered in the same order as in chordates, but each ParaHox gene is actually found on a different genomic scaffold (>300 kb each), which suggests that they are not linked into a single coherent cluster. Therefore, ParaHox genes are dispersed in the genome and are used during embryogenesis in a temporally and spatially coherent manner, whereas the Hox genes, now fully sequenced and annotated, are still linked and are employed as a complex only during the emergence of the adult body plan in the larva.

Neuropeptide Y-family receptors Y6 and Y7 in chicken. Cloning, pharmacological characterization, tissue distribution and conserved synteny with human chromosome region

The peptides of the neuropeptide Y (NPY) family exert their functions, including regulation of appetite and circadian rhythm, by binding to G-protein coupled receptors. Mammals have five subtypes, named Y1, Y2, Y4, Y5 and Y6, and recently Y7 has been discovered in fish and amphibians. In chicken we have previously characterized the first four subtypes and here we describe Y6 and Y7. The genes for Y6 and Y7 are located 1 megabase apart on chromosome 13, which displays conserved synteny with human chromosome 5 that harbours the Y6 gene. The porcine PYY radioligand bound the chicken Y6 receptor with a K(d) of 0.80 +/- 0.36 nm. No functional coupling was demonstrated. The Y6 mRNA is expressed in hypothalamus, gastrointestinal tract and adipose tissue. Porcine PYY bound chicken Y7 with a K(d) of 0.14 +/- 0.01 nm (mean +/- SEM), whereas chicken PYY surprisingly had a much lower affinity, with a Ki of 41 nm, perhaps as a result of its additional amino acid at the N terminus. Truncated peptide fragments had greatly reduced affinity for Y7, in agreement with its closest relative, Y2, in chicken and fish, but in contrast to Y2 in mammals. This suggests that in mammals Y2 has only recently acquired the ability to bind truncated PYY. Chicken Y7 has a much more restricted tissue distribution than other subtypes and was only detected in adrenal gland. Y7 seems to have been lost in mammals. The physiological roles of Y6 and Y7 remain to be identified, but our phylogenetic and chromosomal analyses support the ancient origin of these Y receptor genes by chromosome duplications in an early (pregnathostome) vertebrate ancestor.

A low diversity of ANTP class homeobox genes in Placozoa

Homeobox genes of the ANTP and PRD classes play important roles in body patterning of metazoans, and a large diversity of these genes have been described in bilaterian animals and cnidarians. Trichoplax adhaerens (Phylum Placozoa) is a small multicellular marine animal with one of the simplest body organizations of all metazoans, showing no symmetry and a small number of distinct cell types. Only two ANTP class genes have been described from Trichoplax: the Hox/ParaHox gene Trox-2 and a gene related to the Not family. Here we report an extensive screen for ANTP class genes in Trichoplax, leading to isolation of three additional ANTP class genes. These can be assigned to the Dlx, Mnx and Hmx gene families. Sequencing approximately 12-20 kb around each gene indicates that none are part of tight gene clusters, and in situ hybridization reveals that at least two have spatially restricted expression around the periphery of the animal. The low diversity of ANTP class genes isolated in Trichoplax can be reconciled with the low anatomical complexity of this animal, although the finding that these genes are assignable to recognized gene families is intriguing.

An ancient Wnt-Dickkopf antagonism in Hydra

The dickkopf (dkk) gene family encodes secreted antagonists of Wnt signalling proteins, which have important functions in the control of cell fate, proliferation, and cell polarity during development. Here, we report the isolation, from a regeneration-specific signal peptide screen, of a novel dickkopf gene from the fresh water cnidarian Hydra. Comparative sequence analysis demonstrates that the Wnt antagonistic subfamily Dkk1/Dkk2/Dkk4 and the non-modulating subfamily Dkk3 separated prior to the divergence of cnidarians and bilaterians. In steady-state Hydra, hydkk1/2/4-expression is inversely related to that of hywnt3a. hydkk1/2/4 is an early injury and regeneration responsive gene, and hydkk1/2/4-expressing gland cells are essential for head regeneration in Hydra, although once the head has regenerated they are excluded from it. Activation of Wnt/beta-Catenin signalling leads to the complete downregulation of hydkk1/2/4 transcripts. When overexpressed in Xenopus, HyDkk1/2/4 has similar Wnt-antagonizing activity to the Xenopus gene Dkk1. Based on the corresponding expression patterns of hydkk1/2/4 and neuronal genes, we suggest that the body column of Hydra is a neurogenic environment suppressing Wnt signalling and facilitating neurogenesis.

Hox genes and their candidate downstream targets in the developing central nervous system

1. Homeobox (Hox) genes were originally discovered in the fruit fly Drosophila, where they function through a conserved homeodomain as transcriptional regulators to control embryonic morphogenesis. Since then over 1000 homeodomain proteins have been identified in several species. In vertebrates, 39 Hox genes have been identified as homologs of the original Drosophila complex, and like their Drosophila counterparts they are organized within chromosomal clusters. Vertebrate Hox genes have also been shown to play a critical role in embryonic development as transcriptional regulators. 2. Both the Drosophila and vertebrate Hox genes have been shown to interact with various cofactors, such as the TALE homeodomain proteins, in recognition of consensus sequences within regulatory elements of their target genes. These protein-protein interactions are believed to contribute to enhancing the specificity of target gene recognition in a cell-type or tissue- dependent manner. The regulatory activity of a particular Hox protein on a specific regulatory element is highly variable and dependent on its interacting partners within the transcriptional complex. 3. In vertebrates, Hox genes display spatially restricted patterns of expression within the developing CNS, both along the anterioposterior and dorsoventral axis of the embryo. Their restricted gene expression is suggestive of a regulatory role in patterning of the CNS, as well as in cell specification. Determining the precise function of individual Hox genes in CNS morphogenesis through classical mutational analyses is complicated due to functional redundancy between Hox genes. 4. Understanding the precise mechanisms through which Hox genes mediate embryonic morphogenesis requires the identification of their downstream target genes. Although Hox genes have been implicated in the regulation of several pathways, few target genes have been shown to be under their direct regulatory control. Development of methodologies used for the isolation of target genes and for the analysis of putative targets will be beneficial in establishing the genetic pathways controlled by Hox factors. 5. Within the developing CNS various cell adhesion molecules and signaling molecules have been identified as candidate downstream target genes of Hox proteins. These targets play a role in processes such as cell migration and differentiation, and are implicated in contributing to neuronal processes such as plasticity and/or specification. Hence, Hox genes not only play a role in patterning of the CNS during early development, but may also contribute to cell specification and identity.

An evolutionary history of the FGF superfamily

Fibroblast growth factors (FGF) are associated with multiple developmental and metabolic processes in triploblasts, and perhaps also in diploblasts. The evolution of the FGF superfamily has accompanied the major morphological and functional innovations of metazoan species. The study of FGFs throughout species shows that the FGF superfamily can be subdivided in eight families in present-day organisms and has evolved through phases of gene duplications and gene losses. At least two major expansions of the superfamily can be recognized: a first expansion increased the number of FGFs from one or few archeo-FGFs to eight proto-FGFs, prototypic of the eight families. A second expansion, which took place during euchordate evolution, is associated with genome duplications. It increased the number of members in the families. Subsequent losses reduced that number to the present-day figures.

Conserved synteny between the Ciona genome and human paralogons identifies large duplication events in the molecular evolution of the insulin-relaxin gene family

The aims of the study were to outline the sequence of events that gave rise to the vertebrate insulin-relaxin gene family and the chromosomal regions in which they reside. We analyzed the gene content surrounding the human insulin/relaxin genes with respect to what family they belonged to and if the duplication history of investigated families parallels the evolution of the insulin-relaxin family members. Markov Clustering and phylogenetic analysis were used to determine family identity. More than 15% of the genes belonged to families that have paralogs in the regions, defining two sets of quadruplicate paralogy regions. Thereby, the localization of insulin/relaxin genes in humans is in accordance with those regions on human chromosomes 1, 11, 12, 19q (insulin/insulin-like growth factors) and 1, 6p/15q, 9/5, 19p (insulin-like factors/relaxins) were formed during two genome duplications. We compared the human genome with that of Ciona intestinalis, a species that split from the vertebrate lineage before the two suggested genome duplications. Two insulin-like orthologs were discovered in addition to the already described Ci-insulin gene. Conserved synteny between the Ciona regions hosting the insulin-like genes and the two sets of human paralogons implies their common origin. Linkage of the two human paralogons, as seen in human chromosome 1, as well as the two regions hosting the Ciona insulin-like genes suggests that a segmental duplication gave rise to the region prior to the genome doublings. Thus, preserved gene content provides support that genome duplication(s) in addition to segmental and single-gene duplications shaped the genomes of extant vertebrates.

A pair as a minimum: the two fibroblast growth factors of the nematode Caenorhabditis elegans

Fibroblast growth factors (FGFs) regulate many important developmental and homeostatic physiological events. The FGF superfamily contains several families. In this review, we present recent findings on the two FGFs of the nematode Caenorhabditis elegans from both functional and phylogenic points of view. C. elegans has a single FGFR (EGL-15) with two functionally exclusive isoforms, and two FGFs (LET-756 and EGL-17), which play distinct roles: an essential function for the former, and guidance of the migrating sex myoblasts for the latter. Regulation of homeostasis by control of the fluid balance could be the basis for the essential function of LET-756. Phylogenetic and functional studies suggest that LET-756, like vertebrate FGF9, -16, and -20, belongs to the FGF9 family, whereas EGL-17, like vertebrate FGF8, -17, and -18, could be included in the FGF8 family.

The genomic environment around the Aromatase gene: evolutionary insights

Background

The cytochrome P450 aromatase (CYP19), catalyses the aromatisation of androgens to estrogens, a key mechanism in vertebrate reproductive physiology. A current evolutionary hypothesis suggests that CYP19 gene arose at the origin of vertebrates, given that it has not been found outside this clade. The human CYP19 gene is located in one of the proposed MHC-paralogon regions (HSA15q). At present it is unclear whether this genomic location is ancestral (which would suggest an invertebrate origin for CYP19) or derived (genomic location with no evolutionary meaning). The distinction between these possibilities should help to clarify the timing of the CYP19 emergence and which taxa should be investigated.

Results

Here we determine the "genomic environment" around CYP19 in three vertebrate species Homo sapiens, Tetraodon nigroviridis and Xenopus tropicalis. Paralogy studies and phylogenetic analysis of six gene families suggests that the CYP19 gene region was structured through "en bloc" genomic duplication (as part of the MHC-paralogon formation). Four gene families have specifically duplicated in the vertebrate lineage. Moreover, the mapping location of the different paralogues is consistent with a model of "en bloc" duplication. Furthermore, we also determine that this region has retained the same gene content since the divergence of Actinopterygii and Tetrapods. A single inversion in gene order has taken place, probably in the mammalian lineage. Finally, we describe the first invertebrate CYP19 sequence, from Branchiostoma floridae.

Conclusion

Contrary to previous suggestions, our data indicates an invertebrate origin for the aromatase gene, given the striking conservation pattern in both gene order and gene content, and the presence of aromatase in amphioxus. We propose that CYP19 duplicated in the vertebrate lineage to yield four paralogues, followed by the subsequent loss of all but one gene in vertebrate evolution. Finally, we suggest that agnathans and lophotrocozoan protostomes should be investigated for the presence of aromatase.

Extensive duplications of phototransduction genes in early vertebrate evolution correlate with block (chromosome) duplications

Many gene families in mammals have members that are expressed more or less uniquely in the retina or differentially in specific retinal cell types. We describe here analyses of nine such gene families with regard to phylogenetic relationships and chromosomal location. The families are opsins, G proteins (alpha, beta, and gamma subunits), phosphodiesterases type 6, cyclic nucleotide-gated channels, G-protein-coupled receptor kinases, arrestins, and recoverins. The results suggest that multiple new gene copies arose in all of these families very early in vertebrate evolution during a period with extensive gene duplications. Many of the new genes arose through duplications of large chromosome regions (blocks of genes) or even entire chromosomes, as shown by linkage with other gene families. Some of the phototransduction families belong to the same duplicated regions and were thus duplicated simultaneously. We conclude that gene duplications in early vertebrate evolution probably helped facilitate the specialization of the retina and the subspecialization of different retinal cell types.

Phylogenetic analysis of Ciona intestinalis gene superfamilies supports the hypothesis of successive gene expansions

Understanding the formation of metazoan multigene families is a good approach to reconstitute the evolution of the chordate genome. In this attempt, the analysis of the genome of selected species provides valuable information. Ciona intestinalis belongs to the urochordates, whose lineage separated from the chordate lineage that later gave birth to vertebrates. We have searched available sequences from the small marine ascidian C. intestinalis for orthologs of members of five vertebrate superfamilies, including tyrosine kinase receptors, ETS, FOX and SOX transcription factors, and WNT secreted regulatory factors, and conducted phylogenetic analyses. We have found that most vertebrate subfamilies have a single C. intestinalis ortholog. Our results support the hypothesis of a gene expansion prior the base of chordate ancestry followed by another gene expansion during vertebrate evolution. They also indicate that Ciona intestinalis genome will be a very valuable tool for evolutionary analyses.

Molecular evolution of NPY receptor subtypes

The neuropeptide Y (NPY) system consists in mammals of three peptides and 4-5 G-protein-coupled receptors called Y receptors that are involved in a variety of physiological functions such as appetite regulation, circadian rhythm and anxiety. Both the receptor family and the peptide family display unexpected evolutionary complexity and flexibility as shown by information from different classes of vertebrates. The vertebrate ancestor most likely had a single peptide gene and three Y receptor genes, the progenitors of the Y1, Y2 and Y5 subfamilies. The receptor genes were probably located in the same chromosomal segment. Additional gene copies arose through the chromosome quadruplication that took place before the emergence of jawed vertebrates (gnathostomes) whereupon differential losses of the gene copies ensued. The inferred ancestral gnathostome gene repertoire most likely consisted of two peptide genes, NPY and PYY, and no less than seven Y receptor genes: four Y1-like (Y1, Y4/a, Y6, and Yb), two Y2-like (Y2 and Y7), and a single Y5 gene. Whereas additional peptide genes have arisen in various lineages, the most common trend among the Y receptor genes has been further losses. Mammals have lost Yb and Y7 (the latter still exists in frogs) and Y6 is a pseudogene in several mammalian species but appears to be still functional in some. One challenge is to find out if mammals have been deprived of any functions through these gene losses. Teleost fishes like zebrafish and pufferfish, on the other hand, have lost the two major appetite-stimulating receptors Y1 and Y5. Nevertheless, teleost fishes seem to respond to NPY with increased feeding why some other subtype probably mediates this effect. Another challenge is to deduce how Y2 and Y4 came to evolve an inhibitory effect on appetite. Changes in anatomical distribution of receptor expression may have played an important part in such functional switching along with changes in receptor structures and ligand preferences.

Paired-type homeodomain transcription factors are imported into the nucleus by karyopherin 13

We report that the paired homeodomain transcription factor Pax6 is imported into the nucleus by the Karyopherin beta family member Karyopherin 13 (Kap13). Pax6 was identified as a potential cargo for Kap13 by a yeast two-hybrid screen. Direct binding of Pax6 to Kap13 was subsequently confirmed by in vitro assays with recombinant proteins, and binding in vivo was shown by coimmunoprecipitation. Ran-dependent import of Pax6 by Kap13 was shown to occur by using a digitonin-permeabilized cells assay. Kap13 binds to Pax6 via a nuclear localization sequence (NLS), which is located within a segment of 80 amino acid residues that includes the homeodomain. Kap13 showed reduced binding to Pax6 when either region located at each end of the homeodomain (208 to 214 and 261 to 267) was deleted. The paired-type homeodomain transcription factor family includes more than 20 members. All members contain a region similar to the NLS found in Pax6 and are therefore likely to be imported by Kap13. We confirmed this hypothesis for Pax3 and Crx, which bind to and are imported by Kap13.

Recent developments in computational approaches for uncovering genomic homology

Identifying genomic homology within and between genomes is essential when studying genome evolution. In the past years, different computational techniques have been developed to detect homology even when the actual similarity between homologous segments is low. Depending on the strategy used, these methods search for pairs of chromosomal segments between which either both gene content and order are conserved or gene content only. However, due to fact that, after their divergence, homologous segments can lose a different set of genes, these methods still often fail to detect genomic homology. Recently, more advanced approaches have been developed that can combine gene order and content information of multiple genomic segments.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and AmphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.

A genomewide survey of developmentally relevant genes in Ciona intestinalis. II. Genes for homeobox transcription factors

Homeobox-containing genes play crucial roles in various developmental processes, including body-plan specification, pattern formation and cell-type specification. The present study searched the draft genome sequence and cDNA/EST database of the basal chordate Ciona intestinalis to identify 83 homeobox-containing genes in this animal. This number of homeobox genes in the Ciona genome is smaller than that in the Caenorhabditis elegans, Drosophila melanogaster, human and mouse genomes. Of the 83 genes, 76 have possible human orthologues and 7 may be unique to Ciona. The ascidian homeobox genes were classified into 11 classes, including Hox class, NK class, Paired class, POU class, LIM class, TALE class, SIX class, Prox class, Cut class, ZFH class and HNF1 class, according to the classification scheme devised for known homeobox genes. As to the Hox cluster, the Ciona genome contains single copies of each of the paralogous groups, suggesting that there is a single Hox cluster, if any, but genes orthologous to Hox7, 8, 9 and 11 were not found in the genome. In addition, loss of genes had occurred independently in the Ciona lineage and was noticed in Gbx of the EHGbox subclass, Sax, NK3, Vax and vent of the NK class, Cart, Og9, Anf and Mix of the Paired class, POU-I, III, V and VI of the POU class, Lhx6/7 of the LIM class, TGIF of the TALE class, Cux and SATB of the Cut class, and ZFH1 of the ZFH class, which might have reduced the number of Ciona homeobox genes. Interestingly, one of the newly identified Ciona intestinalis genes and its vertebrate counterparts constitute a novel subclass of HNF1 class homeobox genes. Furthermore, evidence for the gene structures and expression of 54 of the 83 homeobox genes was provided by analysis of ESTs, suggesting that cDNAs for these 54 genes are available. The present data thus reveal the repertoire of homeodomain-containing transcription factors in the Ciona genome, which will be useful for future research on the development and evolution of chordates.

ParaDB: a tool for paralogy mapping in vertebrate genomes

We present ParaDB (http://abi.marseille.inserm.fr/paradb/), a new database for large-scale paralogy studies in vertebrate genomes. We intended to collect all information (sequence, mapping and phylogenetic data) needed to map and detect new paralogous regions, previously defined as Paralogons. The AceDB database software was used to generate graphical objects and to organize data. General data were automatically collated from public sources (Ensembl, GadFly and RefSeq). ParaDB provides access to data derived from whole genome sequences (Homo sapiens, Mus musculus and Drosophila melanogaster): cDNA and protein sequences, positional information, bibliographical links. In addition, we provide BLAST results for each protein sequence, InParanoid orthologs and 'In-Paralogs' data, previously established paralogy data, and, to compare vertebrates and Drosophila, orthology data.

The human Hox-bearing chromosome regions did arise by block or chromosome (or even genome) duplications

Many chromosome regions in the human genome exist in four similar copies, suggesting that the entire genome was duplicated twice in early vertebrate evolution, a concept called the 2R hypothesis. Forty-two gene families on the four Hox-bearing chromosomes were recently analyzed by others, and 32 of these were reported to have evolutionary histories incompatible with duplications concomitant with the Hox clusters, thereby contradicting the 2R hypothesis. However, we show here that nine of the families have probably been translocated to the Hox-bearing chromosomes more recently, and that three of these belong to other chromosome quartets where they actually support the 2R hypothesis. We consider 13 families too complex to shed light on the chromosome duplication hypothesis. Among the remaining 20 families, 14 display phylogenies that support or are at least consistent with the Hox-cluster duplications. Only six families seem to have other phylogenies, but these trees are highly uncertain due to shortage of sequence information. We conclude that all relevant and analyzable families support or are consistent with block/chromosome duplications and that none clearly contradicts the 2R hypothesis.

Were vertebrates octoploid?

It has long been suggested that gene and genome duplication play important roles in the evolution of organismal complexity. For example, work by Ohno proposed that two rounds of whole genome doubling (tetraploidy) occurred during the evolution of vertebrates: the extra genes permitting an increase in physiological and anatomical complexity. Several modifications of this 'two tetraploidies' hypothesis have been proposed, taking into account accumulating data, and there is wide acceptance of the basic scheme. In the past few years, however, several authors have raised doubts, citing lack of direct support or even evidence to the contrary. Here, we review the evidence for and against the occurrence of tetraploidies in early vertebrate evolution, and present a new compilation of molecular phylogenetic data for amphioxus. We argue that evidence in favour of tetraploidy, based primarily on genome and gene family analyses, is strong. Furthermore, we show that two observations used as evidence against genome duplication are in fact compatible with the hypothesis: but only if the genome doubling occurred by two closely spaced sequential rounds of autotetraploidy. We propose that early vertebrates passed through an autoautooctoploid phase in the evolution of their genomes.

Missense mutations of human homeoboxes: A review

The homeodomain (encoded by the homeobox) is the DNA-binding domain of a large variety of transcriptional regulators involved in controlling cell fate decisions and development. Mutations of homeobox-containing genes cause several diseases in humans. A variety of missense mutations giving rise to human diseases have been described. These mutations are an excellent model to better understand homeodomain molecular functions. To this end, homeobox missense mutations giving rise to human diseases are reviewed. Seventy-four independent homeobox mutations have been observed in 17 different genes. In the same genes, 30 missense mutations outside the homeobox have been observed, indicating that the homeodomain is more easily affected by single amino acids changes than the rest of the protein. Most missense mutations have dominant effects. Several data indicate that dominance is mostly due to haploinsufficiency. Among proteins having the homeodomain as the only DNA-binding domain, three "hot spot" regions can be delineated: 1) at codon encoding for Arg5; 2) at codon encoding for Arg31; and 3) at codons encoding for amino acids of recognition helix. In the latter, mutations at codons encoding for Arg residues at positions 52 and 53 are prevalent. In the recognition helix, Arg residues at positions 52 and 53 establish contacts with phosphates in the DNA backbone. Missense mutations of amino acids that contribute to sequence discrimination (such as those at positions 50 and 54) are present only in a minority of cases. Similar data have been obtained when missense mutations of proteins possessing an additional DNA-binding domain have been analyzed. The only exception is observed in the POU1F1 (PIT1) homeodomain, in which Arg58 is a "hot spot" for mutations, but is not involved in DNA recognition.