Hox cluster organization in the jawless vertebrate Petromyzon marinus

Molecular Phylogenetic Analysis of the AIG Family in Vertebrates

Androgen-inducible genes (AIGs), which can be regulated by androgen level, constitute a group of genes characterized by the presence of the AIG/FAR-17a domain in its protein sequence. Previous studies on AIGs demonstrated that one member of the gene family, AIG1, is involved in many biological processes in cancer cell lines and that ADTRP is associated with cardiovascular diseases. It has been shown that the numbers of AIG paralogs in humans, mice, and zebrafish are 2, 2, and 3, respectively, indicating possible gene duplication events during vertebrate evolution. Therefore, classifying subgroups of AIGs and identifying the homologs of each AIG member are important to characterize this novel gene family further. In this study, vertebrate AIGs were phylogenetically grouped into three major clades, ADTRP, AIG1, and AIG-L, with AIG-L also evident in an outgroup consisting of invertebrsate species. In this case, AIG-L, as the ancestral AIG, gave rise to ADTRP and AIG1 after two rounds of whole-genome duplications during vertebrate evolution. Then, the AIG family, which was exposed to purifying forces during evolution, lost or gained some of its members in some species. For example, in eutherians, Neognathae, and Percomorphaceae, AIG-L was lost; in contrast, Salmonidae and Cyprinidae acquired additional AIG copies. In conclusion, this study provides a comprehensive molecular phylogenetic analysis of vertebrate AIGs, which can be employed for future functional characterization of AIGs.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERβ) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

Conservation of Hox gene clusters in the self-fertilizing fish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae)

In this study, whole Hox gene clusters in the self-fertilizing mangrove killifish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae), a unique hermaphroditic vertebrate in which both sex organs are functional at the same time, were identified from whole genome and transcriptome sequences. The aim was to increase the understanding of the evolutionary status of conservation of this Hox gene cluster across fish species.

Evolution and functional diversity of aquaporins

In this review, we provide a brief synopsis of the evolution and functional diversity of the aquaporin gene superfamily in prokaryotic and eukaryotic organisms. Based upon the latest data, we discuss the expanding list of molecules shown to permeate the central pore of aquaporins, and the unexpected diversity of water channel genes in Archaea and Bacteria. We further provide new insight into the origin by horizontal gene transfer of plant glycerol-transporting aquaporins (NIPs), and the functional co-option and gene replacement of insect glycerol transporters. Finally, we discuss the origins of four major grades of aquaporins in Eukaryota, together with the increasing repertoires of aquaporins in vertebrates.

Evolution of retinoic acid receptors in chordates: insights from three lamprey species, Lampetra fluviatilis, Petromyzon marinus, and Lethenteron japonicum

Background

Retinoic acid (RA) signaling controls many developmental processes in chordates, from early axis specification to late organogenesis. The functions of RA are chiefly mediated by a subfamily of nuclear hormone receptors, the retinoic acid receptors (RARs), that act as ligand-activated transcription factors. While RARs have been extensively studied in jawed vertebrates (that is, gnathostomes) and invertebrate chordates, very little is known about the repertoire and developmental roles of RARs in cyclostomes, which are extant jawless vertebrates. Here, we present the first extensive study of cyclostome RARs focusing on three different lamprey species: the European freshwater lamprey, Lampetra fluviatilis, the sea lamprey, Petromyzon marinus, and the Japanese lamprey, Lethenteron japonicum.

Results

We identified four rar paralogs (rar1, rar2, rar3, and rar4) in each of the three lamprey species, and phylogenetic analyses indicate a complex evolutionary history of lamprey rar genes including the origin of rar1 and rar4 by lineage-specific duplication after the lamprey-hagfish split. We further assessed their expression patterns during embryonic development by in situ hybridization. The results show that lamprey rar genes are generally characterized by dynamic and highly specific expression domains in different embryonic tissues. In particular, lamprey rar genes exhibit combinatorial expression domains in the anterior central nervous system (CNS) and the pharyngeal region.

Conclusions

Our results indicate that the genome of lampreys encodes at least four rar genes and suggest that the lamprey rar complement arose from vertebrate-specific whole genome duplications followed by a lamprey-specific duplication event. Moreover, we describe a combinatorial code of lamprey rar expression in both anterior CNS and pharynx resulting from dynamic and highly specific expression patterns during embryonic development. This ‘RAR code’ might function in regionalization and patterning of these two tissues by differentially modulating the expression of downstream effector genes during development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0016-4) contains supplementary material, which is available to authorized users.

The CXC chemokine receptors of fish: Insights into CXCR evolution in the vertebrates

This article will review current knowledge on CXCR in fish, that represent three distinct vertebrate groups: Agnatha (jawless fishes), Chondrichthyes (cartilaginous fishes) and Osteichthyes (bony fishes). With the sequencing of many fish genomes, information on CXCR in these species in particular has expanded considerably. In mammals, 6 CXCRs have been described, and their homologues will be initially reviewed before considering a number of atypical CXCRs and a discussion of CXCR evolution.

Evolution of endothelin receptors in vertebrates

Endothelin receptors are G protein coupled receptors (GPCRs) of the β-group of rhodopsin receptors that bind to endothelin ligands, which are 21 amino acid long peptides derived from longer prepro-endothelin precursors. The most basal Ednr-like GPCR is found outside vertebrates in the cephalochordate amphioxus, but endothelin ligands are only present among vertebrates, including the lineages of jawless vertebrates (lampreys and hagfishes), cartilaginous vertebrates (sharks, rays, and chimaeras), and bony vertebrates (ray-finned fishes and lobe-finned vertebrates including tetrapods). A bona fide endothelin system is thus a vertebrate-specific innovation with important roles for regulating the cardiovascular system, renal and pulmonary processes, as well as for the development of the vertebrate-specific neural crest cell population and its derivatives. Expectedly, dysregulation of endothelin receptors and the endothelin system leads to a multitude of human diseases. Despite the importance of different types of endothelin receptors for vertebrate development and physiology, current knowledge on endothelin ligand-receptor interactions, on the expression of endothelin receptors and their ligands, and on the functional roles of the endothelin system for embryonic development and in adult vertebrates is very much biased towards amniote vertebrates. Recent analyses from a variety of vertebrate lineages, however, have shown that the endothelin system in lineages such as teleost fish and lampreys is more diverse and is divergent from the mammalian endothelin system. This diversity is mainly based on differential evolution of numerous endothelin system components among vertebrate lineages generated by two rounds of whole genome duplication (three in teleosts) during vertebrate evolution. Here we review current understanding of the evolutionary history of the endothelin receptor family in vertebrates supplemented with surveys on the endothelin receptor gene complement of newly available genome assemblies from phylogenetically informative taxa. Our assessment further highlights the diversity of the vertebrate endothelin system and calls for detailed functional and pharmacological analyses of the endothelin system beyond tetrapods.

Connectivity of vertebrate genomes: Paired-related homeobox (Prrx) genes in spotted gar, basal teleosts, and tetrapods

Teleost fish are important models for human biology, health, and disease. Because genome duplication in a teleost ancestor (TGD) impacts the evolution of teleost genome structure and gene repertoires, we must discriminate gene functions that are shared and ancestral from those that are lineage-specific in teleosts or tetrapods to accurately apply inferences from teleost disease models to human health. Generalizations must account both for the TGD and for divergent evolution between teleosts and tetrapods after the likely two rounds of genome duplication shared by all vertebrates. Progress in sequencing techniques provides new opportunities to generate genomic and transcriptomic information from a broad range of phylogenetically informative taxa that facilitate detailed understanding of gene family and gene function evolution. We illustrate here the use of new sequence resources from spotted gar (Lepisosteus oculatus), a rayfin fish that diverged from teleosts before the TGD, as well as RNA-Seq data from gar and multiple teleost lineages to reconstruct the evolution of the Paired-related homeobox (Prrx) transcription factor gene family, which is involved in the development of mesoderm and neural crest-derived mesenchyme. We show that for Prrx genes, the spotted gar genome and gene expression patterns mimic mammals better than teleosts do. Analyses force the seemingly paradoxical conclusion that regulatory mechanisms for the limb expression domains of Prrx genes existed before the evolution of paired appendages. Detailed evolutionary analyses like those reported here are required to identify fish species most similar to the human genome to optimally connect fish models to human gene functions in health and disease.

Essential roles of zebrafish rtn4/Nogo paralogues in embryonic development

Background

As a consequence of gene/genome duplication, the RTN4/Nogo gene has two counterparts in zebrafish: rtn4a and rtn4b. The shared presence of four specific amino acid motifs—M1 to M4—in the N-terminal region of mammalian RTN4, and zebrafish Rtn4b suggests that Rtn4b is the closest homologue of mammalian Nogo-A.

Results

To explore their combined roles in zebrafish development, we characterized the expression patterns of rtn4a and rtn4b in a comparative manner and performed morpholino-mediated knockdowns. Although both genes were coexpressed in the neural tube and developing brain at early stages, they progressively acquired distinct expression domains such as the spinal cord (rtn4b) and somites (rtn4a). Downregulation of rtn4a and rtn4b caused severe brain abnormalities, with rtn4b knockdown severely affecting the spinal cord and leading to immobility. In addition, the retinotectal projection was severely affected in both morphants, as the retina and optic tectum appeared smaller and only few retinal axons reached the abnormally reduced tectal neuropil. The neuronal defects were more persistent in rtn4b morphants. Moreover, the latter often lacked pectoral fins and lower jaws and had malformed branchial arches. Notably, these defects led to larval death in rtn4b, but not in rtn4a morphants.

Conclusions

In contrast to mammalian Nogo-A, its zebrafish homologues, rtn4a and particularly rtn4b, are essential for embryonic development and patterning of the nervous system.

Cloning, phylogeny, and regional expression of a Y5 receptor mRNA in the brain of the sea lamprey (Petromyzon marinus)

The NPY receptors known as Y receptors are classified into three subfamilies, Y1, Y2, and Y5, and are involved in different physiological functions. The Y5 receptor is the only member of the Y5 subfamily, and it is present in all vertebrate groups, except for teleosts. Both molecular and pharmacological studies show that Y5 receptor is highly conserved during vertebrate evolution. Furthermore, this receptor is widely expressed in the mammalian brain, including the hypothalamus, where it is thought to take part in feeding and homeostasis regulation. Lampreys belong to the agnathan lineage, and they are thought to have branched out between the two whole-genome duplications that occurred in vertebrates. Therefore, they are in a key position for studies on the evolution of gene families in vertebrates. Here we report the cloning, phylogeny, and brain expression pattern of the sea lamprey Y5 receptor. In phylogenetic studies, the lamprey Y5 receptor clusters in a basal position, together with Y5 receptors of other vertebrates. The mRNA of this receptor is broadly expressed in the lamprey brain, being especially abundant in hypothalamic areas. Its expression pattern is roughly similar to that reported for other vertebrates and parallels the expression pattern of the Y1 receptor subtype previously described by our group, as it occurs in mammals. Altogether, these results confirm that a Y5 receptor is present in lampreys, thus being highly conserved during the evolution of vertebrates, and suggest that it is involved in many brain functions, the only known exception being teleosts.

Evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum)

Cyclostomes, comprising jawless vertebrates such as lampreys and hagfishes, are the sister group of living jawed vertebrates (gnathostomes) and hence an important group for understanding the origin and diversity of vertebrates. In vertebrates and other metazoans, Hox genes determine cell fate along the anteroposterior axis of embryos and are implicated in driving morphological diversity. Invertebrates contain a single Hox cluster (either intact or fragmented), whereas elephant shark, coelacanth, and tetrapods contain four Hox clusters owing to two rounds of whole-genome duplication ("1R" and "2R") during early vertebrate evolution. By contrast, most teleost fishes contain up to eight Hox clusters because of an additional "teleost-specific" genome duplication event. By sequencing bacterial artificial chromosome (BAC) clones and the whole genome, here we provide evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum). This suggests that the lamprey lineage has experienced an additional genome duplication after 1R and 2R. The relative age of lamprey and human paralogs supports this hypothesis. Compared with gnathostome Hox clusters, lamprey Hox clusters are unusually large. Several conserved noncoding elements (CNEs) were predicted in the Hox clusters of lamprey, elephant shark, and human. Transgenic zebrafish assay indicated the potential of CNEs to function as enhancers. Interestingly, CNEs in individual lamprey Hox clusters are frequently conserved in multiple Hox clusters in elephant shark and human, implying a many-to-many orthology relationship between lamprey and gnathostome Hox clusters. Such a relationship suggests that the first two rounds of genome duplication may have occurred independently in the lamprey and gnathostome lineages.

Evolution of Hox gene clusters in deuterostomes

Hox genes, with their similar roles in animals as evolutionarily distant as humans and flies, have fascinated biologists since their discovery nearly 30 years ago. During the last two decades, reports on Hox genes from a still growing number of eumetazoan species have increased our knowledge on the Hox gene contents of a wide range of animal groups. In this review, we summarize the current Hox inventory among deuterostomes, not only in the well-known teleosts and tetrapods, but also in the earlier vertebrate and invertebrate groups. We draw an updated picture of the ancestral repertoires of the different lineages, a sort of “genome Hox bar-code” for most clades. This scenario allows us to infer differential gene or cluster losses and gains that occurred during deuterostome evolution, which might be causally linked to the morphological changes that led to these widely diverse animal taxa. Finally, we focus on the challenging family of posterior Hox genes, which probably originated through independent tandem duplication events at the origin of each of the ambulacrarian, cephalochordate and vertebrate/urochordate lineages.

Evolutionary origin of the structure and function of gonadotropin-inhibitory hormone: insights from lampreys

Gonadotropin (GTH)-inhibitory hormone (GnIH) is a novel hypothalamic neuropeptide that inhibits GTH secretion in mammals and birds by acting on gonadotropes and GnRH neurons within the hypothalamic-pituitary-gonadal axis. GnIH and its orthologs that have an LPXRFamide (X = L or Q) motif at the C terminus (LPXRFamide peptides) have been identified in representative species of gnathostomes. However, the identity of an LPXRFamide peptide had yet to be identified in agnathans, the most ancient lineage of vertebrates, leaving open the question of the evolutionary origin of GnIH and its ancestral function(s). In this study, we identified an LPXRFamide peptide gene encoding three peptides (LPXRFa-1a, LPXRFa-1b, and LPXRFa-2) from the brain of sea lamprey by synteny analysis and cDNA cloning, and the mature peptides by immunoaffinity purification and mass spectrometry. The expression of lamprey LPXRFamide peptide precursor mRNA was localized in the brain and gonad by RT-PCR and in the hypothalamus by in situ hybridization. Immunohistochemistry showed appositions of lamprey LPXRFamide peptide immunoreactive fibers in close proximity to GnRH-III neurons, suggesting that lamprey LPXRFamide peptides act on GnRH-III neurons. In addition, lamprey LPXRFa-2 stimulated the expression of lamprey GnRH-III protein in the hypothalamus and GTHβ mRNA expression in the pituitary. Synteny and phylogenetic analyses suggest that the LPXRFamide peptide gene diverged from a common ancestral gene likely through gene duplication in the basal vertebrates. These results suggest that one ancestral function of LPXRFamide peptides may be stimulatory compared with the inhibitory function seen in later-evolved vertebrates (birds and mammals).

Hox gene clusters of early vertebrates: do they serve as reliable markers for genome evolution?

Hox genes, responsible for regional specification along the anteroposterior axis in embryogenesis, are found as clusters in most eumetazoan genomes sequenced to date. Invertebrates possess a single Hox gene cluster with some exceptions of secondary cluster breakages, while osteichthyans (bony vertebrates) have multiple Hox clusters. In tetrapods, four Hox clusters, derived from the so-called two-round whole genome duplications (2R-WGDs), are observed. Overall, the number of Hox gene clusters has been regarded as a reliable marker of ploidy levels in animal genomes. In fact, this scheme also fits the situations in teleost fishes that experienced an additional WGD. In this review, I focus on cyclostomes and cartilaginous fishes as lineages that would fill the gap between invertebrates and osteichthyans. A recent study highlighted a possible loss of the HoxC cluster in the galeomorph shark lineage, while other aspects of cartilaginous fish Hox clusters usually mark their conserved nature. In contrast, existing resources suggest that the cyclostomes exhibit a different mode of Hox cluster organization. For this group of species, whose genomes could have differently responded to the 2R-WGDs from jawed vertebrates, therefore the number of Hox clusters may not serve as a good indicator of their ploidy level.

Analysis of lamprey clustered Fox genes: insight into Fox gene evolution and expression in vertebrates

In the human genome, members of the FoxC, FoxF, FoxL1, and FoxQ1 gene families are found in two paralagous clusters. One cluster contains the genes FOXQ1, FOXF2, FOXC1 and the second consists of FOXF1, FOXC2, and FOXL1. In jawed vertebrates these genes are known to be expressed in different pharyngeal tissues and all, except FoxQ1, are involved in patterning the early embryonic mesoderm. We have previously traced the evolution of this cluster in the bony vertebrates, and the gene content is identical in the dogfish, a member of the most basally branching lineage of the jawed vertebrates. Here we extend these analyses to jawless vertebrates. Using genomic searches and molecular approaches we have identified homologues of these genes from lampreys. We identify two FoxC genes, two FoxF genes, two FoxQ1 genes and single FoxL1 gene. We examine the embryonic expression of one predominantly mesodermally expressed gene family, FoxC, and the endodermally expressed member of the cluster, FoxQ1. We identified FoxQ1 transcripts in the pharyngeal endoderm, while the two FoxC genes are differentially expressed in the pharyngeal mesenchyme and ectoderm. Furthermore we identify conserved expression of lamprey FoxC genes in the paraxial and intermediate mesoderms. We interpret our results through a chordate-wide comparison of expression patterns and discuss gene content in the context of theories on the evolution of the vertebrate genome.

A family tree of vertebrate chemokine receptors for a unified nomenclature

Chemokines receptors are involved in the recruitment of various cell types in inflammatory and physiological conditions. There are 23 known chemokine receptor genes in the human genome. However, it is still unclear how many chemokine receptors exist in the genomes of various vertebrate species other than human and mouse. Moreover, the orthologous relationships are often obscure between the genes of higher and lower vertebrates. In order to provide a basis for a unified nomenclature system of the vertebrate chemokine receptor gene family, we have analysed the chemokine receptor genes from the genomes of 16 vertebrate species, and classify them into 29 orthologous groups using phylogenetic and comparative genomic analyses. The results reveal a continuous gene birth and death process during the vertebrate evolution and an interesting evolutionary history of the chemokine receptor genes after the emergence in agnathans.

Conserved synteny and the zebrafish genome

Zebrafish offers significant opportunities for the investigation of vertebrate development, evolution, physiology, and behavior and provides numerous models of human disease. Connecting zebrafish phenogenetic biology to that of humans and other vertebrates, however, requires the proper assignment of gene orthologies. Orthology assignments by phylogenetic analysis or by reciprocal best sequence similarity searches can lead to errors, especially in cases of gene duplication followed by gene loss or rapid lineage-specific gene evolution. Conserved synteny analysis provides a method that helps overcome such problems. Here we describe conserved synteny analysis for zebrafish genes and discuss the Synteny Database, a website specifically designed to identify conserved syntenies for zebrafish that takes into account the teleost genome duplication (TGD). We utilize the Synteny Database to demonstrate its power to resolve our understanding of the evolution of nerve growth factor receptor related genes, including Ngfr and the enigmatic Nradd. Finally, we compare conserved syntenies between zebrafish, stickleback, spotted gar, and human to understand the timing of chromosome rearrangements in teleost genome evolution. An improved understanding of gene histories that comes from the application of tools provided by the Synteny Database can facilitate the connectivity of zebrafish and human genomes.

Multiple neurofilament subunits are present in lamprey CNS

In mammals, there are three neurofilament (NF) subunits (NF-L, NF-M, and NF-H), but it was thought that only a single NF, NF180, exists in lamprey. However, NF180 lacked the ability to self-assemble, suggesting that like mammalian NFs, lamprey NFs are heteropolymers, and that additional NF subunits may exist. The present study provides evidence for the existence of a lamprey NF-L homolog (L-NFL). Genes encoding two new NF-M isoforms (NF132 and NF95) also have been isolated and characterized. With NF180, this makes three NF-M-like isoforms. In situ hybridization showed that all three newly cloned NFs are expressed in spinal cord neurons and in spinal-projecting neurons of the brainstem. Like NF180, there were no KSP multiphosphorylation repeat motifs in the tail regions of NF132 or NF95. NF95 was highly identical to homologous parts of NF180, sharing 2 common pieces of DNA with it. Northern blots suggested that NF95 may be expressed at very low levels in older larvae. The presence of L-NFL in lamprey CNS may support the hypothesis that as in mammals, NFs in lamprey are obligate heteropolymers, in which NF-L is a required subunit.

Speciation of polyploid Cyprinidae fish of common carp, crucian carp, and silver crucian carp derived from duplicated Hox genes

Recent studies on comparative genomics have suggested that a round of fish-specific whole genome duplication (3R) in ray-finned fishes might have occurred around 226-316 Mya. Additional genome duplication, specifically in cyprinids, may have occurred more recently after the divergence of the teleosts. The timing of this event, however, is unknown. To address this question, we sequenced four Hox genes from taxa representing the polyploid Cyprinidae fish, common carp (Cyprinus carpio, 2n=100), crucian carp (Carassius auratus auratus, 2n=100), and silver crucian carp (C. auratus gibelio, 2n=156), and then compared them with known sequences from the diploid Cyprinidae fish, blunt snout bream (Megalobrama amblycephala, 2n=48). Our results showed the presence of two distinct Hox duplicates in the genomes of common and crucian carp. Three distinct Hox sequences, one of them orthologous to a Hox gene in common carp and the other two orthologous to a Hox gene in crucian carp, were isolated in silver crucian carp, indicating a possible hybrid origin of silver crucian carp from crucian and common carp. The gene duplication resulting in the origin of the common ancestor of common and crucian carp likely occurred around 10.9-13.2 Mya. The speciations of common vs. crucian carp and silver crucian vs. crucian carp likely occurred around 8.1-11.4 and 2.3-3.0 Mya, respectively. Finally, nonfunctionalization resulting from point mutations in the coding region is a probable fate for some Hox duplicates. Taken together, these results suggested an evolutionary model for polyploidization in speciation and diversification of polyploid fish.

Retracted: Evidence for duplicated Hox genes in polyploid Cyprinidae fish of common carp, crucian carp and silver crucian carp

Notice of Withdrawal: The following article from the Journal of Experimental Zoology Part B: Molecular and Developmental Evolution, "Evidence for duplicated Hox genes in polyploid Cyprinidae fish of common carp, crucian carp, and silver crucian carp" by Yuan J, He Z, Yuan X, Jiang X, Sun X, Zou S, published online on 29 Sept 2009 in Wiley InterScience (www.interscience.wiley.com), has been withdrawn from publication by agreement between the authors, the journal Editor-in-Chief, Gunter P. Wagner, and Wiley Periodicals, Inc.

Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes

We have sequenced and analyzed Hox gene clusters from elephant shark, a holocephalian cartilaginous fish. Elephant shark possesses 4 Hox clusters with 45 Hox genes that include orthologs for a higher number of ancient gnathostome Hox genes than the 4 clusters in tetrapods and the supernumerary clusters in teleost fishes. Phylogenetic analysis of elephant shark Hox genes from 7 paralogous groups that contain all of the 4 members indicated an ((AB)(CD)) topology for the order of Hox cluster duplication, providing support for the 2R hypothesis (i.e., 2 rounds of whole-genome duplication during the early evolution of vertebrates). Comparisons of noncoding sequences of the elephant shark and human Hox clusters have identified a large number of conserved noncoding elements (CNEs), which represent putative cis-regulatory elements that may be involved in the regulation of Hox genes. Interestingly, in fugu more than 50% of these ancient CNEs have diverged beyond recognition in the duplicated (HoxA, HoxB, and HoxD) as well as the singleton (HoxC) Hox clusters. Furthermore, the b-paralogs of the duplicated fugu Hox clusters are virtually devoid of unique ancient CNEs. In contrast to fugu Hox clusters, elephant shark and human Hox clusters have lost fewer ancient CNEs. If these ancient CNEs are indeed enhancers directing tissue-specific expression of Hox genes, divergence of their sequences in vertebrate lineages might have led to altered expression patterns and presumably the functions of their associated Hox genes.

Major genomic events and their consequences for vertebrate evolution and endocrinology

Comparative studies of proteins often face the problem of distinguishing a true orthologue (species homologue) from a paralogue (a gene duplicate). This identification task is particularly challenging for endocrine peptides and neuropeptides because they are short and usually have several invariant positions. For some peptide families, this has led to a terminology with peptide names relating to the first species where a specific peptide sequence was determined, such as chicken or salmon gonadotropin-releasing hormone, or names that highlight amino acid differences, e.g., Lys-vasopressin. With accumulating information from multiple species, such a terminology becomes almost impenetrable for nonexperts and difficult even for aficionados. The sequenced genomes offer a new way to distinguish orthologues and paralogues, namely by location of the genes relative to neighboring genes on the chromosomes. In addition, the genome databases can ideally provide a complete listing of the family members in each species. Many vertebrate gene families have expanded in the two basal tetraploidizations (2R) and the teleost fish third tetraploidization (3R), after which some vertebrate lineages have lost some of the duplicates. We review here some peptide families (neuropeptide Y, oxytocin-vasopressin, and somatostatin) where genomic information helps simplify nomenclature. This approach is useful also for other gene families, such as peptide receptors.

Insights into cyclostome phylogenomics: pre-2R or post-2R

Interest in understanding the transition from prevertebrates to vertebrates at the molecular level has resulted in accumulating genomic and transcriptomic sequence data for the earliest groups of extant vertebrates, namely, hagfishes (Myxiniformes) and lampreys (Petromyzontiformes). Molecular phylogenetic studies on species phylogeny have revealed the monophyly of cyclostomes and the deep divergence between hagfishes and lampreys (more than 400 million years). In parallel, recent molecular phylogenetic studies have shed light on the complex evolution of the cyclostome genome. This consists of whole genome duplications, shared at least partly with gnathostomes (jawed vertebrates), and cyclostome lineage-specific secondary modifications of the genome, such as gene gains and losses. Therefore, the analysis of cyclostome genomes requires caution in distinguishing between orthology and paralogy in gene molecular phylogeny at the gene family scale, as well as between apomorphic and plesiomorphic genomic traits in larger-scale analyses. In this review, we propose possible ways of improving the resolvability of these evolutionary events, and discuss probable scenarios for cyclostome genome evolution, with special emphasis on the hypothesis that two-round (2R) genome duplication events occurred before the divergence between cyclostomes and gnathostomes, and therefore that a post-2R state is a genomic synapomorphy for all extant vertebrates.

ParaHox cluster evolution--hagfish and beyond

Abstract The ParaHox genes comprise three Hox-related homeobox gene families, found throughout the animals. They were first discovered in the invertebrate chordate amphioxus, where they are tightly clustered. In this paper we carry out a comparative review of ParaHox gene cluster organization among the deuterostomes, and discuss how the recently published hagfish ParaHox clusters fit into current theories about the evolution of this group of genes.

Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after?

Two rounds of whole-genome duplications are thought to have played an important role in the establishment of gene repertoires in vertebrates. These events occurred during chordate evolution after the split of the urochordate and cephalochordate lineages but before the radiation of extant gnathostomes (jawed vertebrates). During this interval, diverse agnathans (jawless fishes), including cyclostomes (hagfishes and lampreys), diverged. However, there is no solid evidence for the timing of these genome duplications in relation to the divergence of cyclostomes from the gnathostome lineage. We conducted cDNA sequencing in diverse early vertebrates for members of homeobox-containing (Dlx and ParaHox) and other gene families that would serve as landmarks for genome duplications. Including these new sequences, we performed a molecular phylogenetic census using the maximum likelihood method for 55 gene families. In most of these gene families, we detected many more gene duplications before the cyclostome-gnathostome split, than after. Many of these gene families (e.g., visual opsins, RAR, Notch) have multiple paralogs in conserved, syntenic genomic regions that must have been generated by large-scale duplication events. Taken together, this indicates that the genome duplications occurred before the cyclostome-gnathostome split. We propose that the redundancy in gene repertoires possessed by all vertebrates, including hagfishes and lampreys, was introduced primarily by genome duplications. Apart from subsequent lineage-specific modifications, these ancient genome duplication events might serve generally to distinguish vertebrates from invertebrates at the genomic level.

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.

Genome duplication and the origin of the vertebrate skeleton

During vertebrate embryonic development, tissue patterning and differentiation are regulated by members of multigene families. Evolutionary expansion of these families is thought to have played a role in the evolution of anatomical complexity, including the origins of new cell and tissue types. A defining feature of vertebrates is an endoskeleton, the primary components of which are cartilage and bone. The molecular control of skeletal patterning has been the subject of intensive investigation for over two decades. More recently, comparative studies of organisms at key phylogenetic positions have highlighted the importance of gene duplication in the evolutionary diversification of connective tissues. Understanding the natural histories of gene families involved in skeletogenesis is therefore central to the issue of vertebrate skeletal evolution.

Retracted: Gene duplication and functional evolution of Hox genes in fishes

With their power to shape animal morphology, few genes have captured the imagination of biologists as much as the evolutionarily conserved members of the Hox clusters. Hox genes encode transcription factors that play a key role in specifying the body plan in metazoans and are therefore essential in explaining patterns of evolutionary diversity. While each Hox cluster contains the same genes among the different mammalian species, this does not happen in ray-finned fish, in which both the number and organization of Hox genes and even Hox clusters are variable. Teleost fishes provide the first unambiguous support for ancient whole-genome duplication (third round) in an animal lineage. The number of genes differs in each cluster as a result of increased freedom to mutate after duplication. This has also allowed them to diverge and to adopt novel developmental roles. In this review, the authors have firstly focused on broadly outlining the duplication of Hoxgenes in fishes and discussing how comparative genomics is elucidating the molecular changes associated with the evolution of Hox genes expression and developmental function in the teleost fishes.Additional related research aspects, such as imaging of roles of microRNAs, chromatin regulation and evolutionary findings are also discussed.

Noncanonical role of Hox14 revealed by its expression patterns in lamprey and shark

Hox genes are arranged in uninterrupted clusters in vertebrate genomes, and the nested patterns of their expression define spatial identities in multiple embryonic tissues. The ancestral Hox cluster of vertebrates has long been thought to consist of, maximally, 13 Hox genes. However, recently, Hox14 genes were discovered in three chordate lineages, the coelacanth, cartilaginous fishes, and amphioxus, but their expression patterns have not yet been analyzed. We isolated Hox14 cDNAs from the Japanese lamprey and cloudy catshark. These genes were not expressed in the central nervous systems, somites, or fin buds/folds but were expressed in a restricted cell population surrounding the hindgut. The lack of Hox14 expression in most of the embryonic axial elements, where nested Hox expressions define spatial identities, suggests a decoupling of Hox14 genes' regulation from the ancestral regulatory mechanism. The relaxation of preexisting constraint for collinear expression may have permitted the secondary losses of this Hox member in the tetrapod and teleost lineages.

Ventralized zebrafish embryo rescue by overexpression of Zic2a

The neuroectoderm arises during gastrulation as a population of undifferentiated proliferating neuroepithelial cells. As development continues, neuroepithelial cells leave the cell cycle and differentiate into neurons and glia of the functioning central nervous system. What processes establish the spatial distribution of proliferating neuroepithelial cells? To investigate this question, zic2a was isolated from zebrafish, a homolog of the Drosophila pair-rule gene odd-paired, which is involved in nervous system patterning. At shield stage, zic2a was expressed in the zebrafish organizer and the blastoderm margin, and became restricted to the axial mesoderm in mid-gastrula. Expression of zic2a appeared in the prospective neuroectoderm during gastrulation, and later demarcated the presumptive forebrain. This expression pattern suggests that zic2a may function early in the organizer and later in the neural plate to demarcate the population of proliferating neuroectoderm. Consistent with a function for zic2a in transducing signals from the organizer, overexpression of zic2a resulted in an expansion of proliferating neuroectoderm. Furthermore, zic2a overexpression rescued the ventralized phenotype of chordino mutant embryos, which lack a functional chordin gene. Early expression of zic2 in the zebrafish organizer, and the phenotype resulting from overexpression, show a role for zic2a downstream of chordin or other secreted organizer proteins in establishing the initial size of the population of neuroectoderm cells.

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.

Hox gene clusters in blunt snout bream, Megalobrama amblycephala and comparison with those of zebrafish, fugu and medaka genomes

Hox genes encode transcription factors that play a key role in specifying the body plan in metazoans and are therefore essential in explaining patterns of evolutionary diversity. While each Hox cluster contains the same genes among the different mammalian species, this does not happen in ray-finned fish, in which both the number and organization of Hox genes and even Hox clusters are variables. Here we reveal the organization of Hox genes loci in blunt snout bream. Forty-nine Hox genes including a pseudogene A9b in total have been found in seven clusters as follows: 8 Hox genes in the Aa cluster; 5 in Ab; 10 in Ba; 4 in Bb; 11 in Ca; 4 in Cb; and 7 in Da. In terms of gene content, clusters organization and sequence similarities of putative amino acids, blunt snout bream is more closely related to zebrafish than to fugu and medaka. In contrast to the situation in fugu and medaka, both blunt snout bream and zebrafish have duplicated HoxC cluster but only a single copy of the HoxD cluster. The result implies that the loss of the second HoxD cluster might be a shared feature of the Ostariophysi, to which zebrafish and blunt snout bream both belong. Phylogenetic analysis bases on the paralogous genes from twin clusters supports the duplication-first model, i.e., four original clusters may have duplicated in an event before the divergence of the blunt snout bream-plus-zebrafish lineage and the fugu-plus-medaka lineage. Additionally, the relationship between the decrease of GC level and the loss of conservation and function of one of the paralogous genes from twin clusters is discussed.

Hox gene expression patterns in Lethenteron japonicum embryos--insights into the evolution of the vertebrate Hox code

The Hox code of jawed vertebrates is characterized by the colinear and rostrocaudally nested expression of Hox genes in pharyngeal arches, hindbrain, somites, and limb/fin buds. To gain insights into the evolutionary path leading to the gnathostome Hox code, we have systematically analyzed the expression pattern of the Hox gene complement in an agnathan species, Lethenteron japonicum (Lj). We have isolated 15 LjHox genes and assigned them to paralogue groups (PG) 1-11, based on their deduced amino acid sequences. LjHox expression during development displayed gnathostome-like spatial patterns with respect to the PG numbers. Specifically, lamprey PG1-3 showed homologous expression patterns in the rostral hindbrain and pharyngeal arches to their gnathostome counterparts. Moreover, PG9-11 genes were expressed specifically in the tailbud, implying its posteriorizing activity as those in gnathostomes. We conclude that these gnathostome-like colinear spatial patterns of LjHox gene expression can be regarded as one of the features already established in the common ancestor of living vertebrates. In contrast, we did not find evidence for temporal colinearity in the onset of LjHox expression. The genomic and developmental characteristics of Hox genes from different chordate species are also compared, focusing on evolution of the complex body plan of vertebrates.

Evolutionary biology of CD1

The recognition more than a decade ago that lipids presented by CD1 could function as T cell antigens revealed a startling and previously unappreciated complexity to the adaptive immune system. The initial novelty of lipid antigen presentation by CD1 has since given way to a broader perspective of the immune system's capacity to sense and respond to a diverse array of macromolecules. Some immune recognition systems such as Toll-like receptors can trace their origins back into the deep history of sea urchins and arthropods. Others such as the major histocompatibility complex (MHC) appear relatively recently and interestingly, only in animals that also possess a jaw. The natural history of CD1 is thus part of the wider story of immune system evolution and should be considered in this context. Most evidence indicates that CD1 probably evolved from a classical MHC class I (MHC I) gene at some point during vertebrate evolution. This chapter reviews the evidence for this phylogenetic relationship and attempts to connect CD1 to existing models of MHC evolution. This endeavor is facilitated today by the recent availability of whole genome sequence data from a variety of species. Investigators have used these data to trace the ultimate origin of the MHC to a series of whole genome duplications that occurred roughly 500 million years ago. Sequence data have also revealed homologs of the mammalian MHC I and MHC II gene families in virtually all jawed vertebrates including sharks, bony fishes, reptiles, and birds. In contrast, CD1 genes have thus far been found only in a subset of these animal groups. This pattern of CD1 occurrence in the genomes of living species suggests the emergence of CD 1 in an early terrestrial vertebrate.

Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4x coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element-like and long interspersed element-like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.

Vertebrate vitellogenin gene duplication in relation to the "3R hypothesis": correlation to the pelagic egg and the oceanic radiation of teleosts

The spiny ray-finned teleost fishes (Acanthomorpha) are the most successful group of vertebrates in terms of species diversity. Their meteoric radiation and speciation in the oceans during the late Cretaceous and Eocene epoch is unprecedented in vertebrate history, occurring in one third of the time for similar diversity to appear in the birds and mammals. The success of marine teleosts is even more remarkable considering their long freshwater ancestry, since it implies solving major physiological challenges when freely broadcasting their eggs in the hyper-osmotic conditions of seawater. Most extant marine teleosts spawn highly hydrated pelagic eggs, due to differential proteolysis of vitellogenin (Vtg)-derived yolk proteins. The maturational degradation of Vtg involves depolymerization of mainly the lipovitellin heavy chain (LvH) of one form of Vtg to generate a large pool of free amino acids (FAA 150-200 mM). This organic osmolyte pool drives hydration of the ooctye while still protected within the maternal ovary. In the present contribution, we have used Bayesian analysis to examine the evolution of vertebrate Vtg genes in relation to the "3R hypothesis" of whole genome duplication (WGD) and the functional end points of LvH degradation during oocyte maturation. We find that teleost Vtgs have experienced a post-R3 lineage-specific gene duplication to form paralogous clusters that correlate to the pelagic and benthic character of the eggs. Neo-functionalization allowed one paralogue to be proteolyzed to FAA driving hydration of the maturing oocytes, which pre-adapts them to the marine environment and causes them to float. The timing of these events matches the appearance of the Acanthomorpha in the fossil record. We discuss the significance of these adaptations in relation to ancestral physiological features, and propose that the neo-functionalization of duplicated Vtg genes was a key event in the evolution and success of the teleosts in the oceanic environment.

Evidence that mechanisms of fin development evolved in the midline of early vertebrates

The origin of paired appendages was a major evolutionary innovation for vertebrates, marking the first step towards fin- (and later limb-) driven locomotion. The earliest vertebrate fossils lack paired fins but have well-developed median fins, suggesting that the mechanisms of fin development were assembled first in the midline. Here we show that shark median fin development involves the same genetic programs that operate in paired appendages. Using molecular markers for different cell types, we show that median fins arise predominantly from somitic (paraxial) mesoderm, whereas paired appendages develop from lateral plate mesoderm. Expression of Hoxd and Tbx18 genes, which specify paired limb positions, also delineates the positions of median fins. Proximodistal development of median fins occurs beneath an apical ectodermal ridge, the structure that controls outgrowth of paired appendages. Each median fin bud then acquires an anteroposteriorly-nested pattern of Hoxd expression similar to that which establishes skeletal polarity in limbs. Thus, despite their different embryonic origins, paired and median fins utilize a common suite of developmental mechanisms. We extended our analysis to lampreys, which diverged from the lineage leading to gnathostomes before the origin of paired appendages, and show that their median fins also develop from somites and express orthologous Hox and Tbx genes. Together these results suggest that the molecular mechanisms for fin development originated in somitic mesoderm of early vertebrates, and that the origin of paired appendages was associated with re-deployment of these mechanisms to lateral plate mesoderm.

Lamprey type II collagen and Sox9 reveal an ancient origin of the vertebrate collagenous skeleton

Type II collagen is the major cartilage matrix protein in the jawed vertebrate skeleton. Lampreys and hagfishes, by contrast, are thought to have noncollagenous cartilage. This difference in skeletal structure has led to the hypothesis that the vertebrate common ancestor had a noncollagenous skeleton, with type II collagen becoming the predominant cartilage matrix protein after the divergence of jawless fish from the jawed vertebrates approximately 500 million years ago. Here we report that lampreys have two type II collagen (Col2alpha1) genes that are expressed during development of the cartilaginous skeleton. We also demonstrate that the adult lamprey skeleton is rich in Col2alpha1 protein. Furthermore, we have isolated a lamprey orthologue of Sox9, a direct transcriptional regulator of Col2alpha1 in jawed vertebrates, and show that it is coexpressed with both Col2alpha1 genes during skeletal development. These results reveal that the genetic pathway for chondrogenesis in lampreys and gnathostomes is conserved through the activation of cartilage matrix molecules and suggest that a collagenous skeleton evolved surprisingly early in vertebrate evolution.

Defining the origins and evolution of the chemokine/chemokine receptor system

The chemokine system has a critical role in mammalian immunity, but the evolutionary history of chemokines and chemokine receptors are ill-defined. We used comparative whole genome analysis of fruit fly, sea urchin, sea squirt, pufferfish, zebrafish, frog, and chicken to identify chemokines and chemokine receptors in each species. We report 127 chemokine and 70 chemokine receptor genes in the 7 species, with zebrafish having the most chemokines, 63, and chemokine receptors, 24. Fruit fly, sea urchin, and sea squirt have no identifiable chemokines or chemokine receptors. This study represents the most comprehensive analysis of the chemokine system to date and the only complete characterization of chemokine systems outside of mouse and human. We establish a clear evolutionary model of the chemokine system and trace the origin of the chemokine system to approximately 650 million years ago, identifying critical steps in their evolution and demonstrating a more extensive chemokine system in fish than previously thought.

A PCR survey for posterior Hox genes in amphibians

Hox genes encode transcription factors that play a key role in specifying the body plan in metazoans and are therefore essential in explaining patterns of evolutionary diversity. As an ancient tetrapod group with diverse limb types, amphibians are important for understanding the origin and diversification of limbs in land vertebrates. We conducted a PCR survey in two species of each amphibian order to identify Hox-9 to Hox-13, known to function in limb development. Fifteen distinct posterior Hox genes and one retro-pseudogene were identified, and the former confirm the existence of four Hox clusters in each amphibian order. Some genes expected to occur in all tetrapods, based on the posterior Hox complement of mammals, fishes and coelacanth, were not recovered from our survey, and may have been lost. Hoxd-12 is absent in frogs and possibly other amphibians. Considering its function in autopodial development, the loss of this gene may be related to the absence of the fifth finger in frogs and salamanders.

Organisation of the lamprey (Lampetra fluviatilis) embryonic brain: insights from LIM-homeodomain, Pax and hedgehog genes

To investigate the embryonic development of the central nervous system of the lamprey Lampetra fluviatilis, we have isolated and analysed the expression patterns of members of the LIM-homeodomain, Pax, Hedgehog and Nkx2.1 families. Using degenerate RT-PCR, single representatives of Lhx1/Lhx5, Lhx2/Lhx9, Pax3/Pax7 and Hedgehog families could be isolated in L. fluviatilis. Expression analysis revealed that the lamprey forebrain presents a clear neuromeric pattern. We describe the existence of 4 embryonic diencephalic prosomeres whose boundaries can be identified by the combined and relative expressions of LfPax37, LfLhx15 and LfLhx29. This suggests that the embryonic lamprey and gnathostome forebrain are patterned in a highly similar manner. Moreover, analysis of the LfHh gene, which is expressed in the hypothalamus, zona limitans intrathalamica and floor plate, reveals the possible molecular origin of this neuromeric brain pattern. By contrast, LfHh and LfNkx2.1 expressions suggest major differences in patterning mechanisms of the ventral telencephalon when compared to gnathostomes. In summary, our findings highlight a neuromeric organisation of the embryonic agnathan forebrain and point to the possible origin of this organisation, which is thus a truly vertebrate character. They also suggest that Hh/Shh midline signalling might act as a driving force for forebrain evolution.