Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

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.

Gene expression in the liver of the hagfish (Eptatretus burgeri) belonging to the Cyclostomata is ancestral to that of mammals

Although the liver of the hagfish, an earliest diverged lineage among vertebrates, has a histological architecture similar to that of mammals, its gene expression has not been explored yet. The present study was undertaken to comparatively characterize gene expression in the liver of the hagfish with that of the mouse, using in situ hybridization technique. Expression of alb (albumin) was detectable in all hepatocytes of the hagfish liver, but was negative in intrahepatic bile ducts. Their expression in abundant periportal ductules was weak. The expression pattern basically resembled that in mammalian livers, indicating that the differential expression of hepatocyte markers in hepatocytes and biliary cells may have been acquired in ancestral vertebrates. alb expression was almost homogeneous in the hagfish liver, whereas that in the mouse liver lobule was zonal. The glul (glutamate-ammonia ligase) expression was also homogeneously detectable in hepatocytes without zonation, and weakly so in biliary cells of the hagfish, which contrasted with its restricted pericentral expression in mouse livers. These findings indicated that the hagfish liver did not have mammalian-type zonation. Whereas tetrapods had Hnf (hepatocyte nuclear factor) 1a and Hnf1b genes encoding the transcription factors, the hagfish had a single gene of their orthologue hnf1. Although HNF1α and HNF1β were immunohistochemically detected in hepatocytes and biliary cells of the mouse, respectively, hnf1 was expressed in both hepatocytes and biliary cells of the hagfish. These data indicate that gene expression of hnf1 in the hagfish liver may be ancestral with that of alb and glul during vertebrate evolution.

Bioinformatic and fine-scale chromosomal mapping reveal the nature and evolution of eliminated chromosomes in the Japanese hagfish, Eptatretus burgeri, through analysis of repetitive DNA families

In the Japanese hagfish, Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n = 52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n = 36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome of E. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1-6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes of E. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7-10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1-10) by fluorescence in situ hybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2-10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

Phylogenetic and functional properties of hagfish neurohypophysial hormone receptors distinct from their jawed vertebrate counterparts

Vertebrate neurohypophysial hormones, i.e., vasopressin- and oxytocin-family peptides, exert versatile physiological actions via distinct G protein-coupled receptors. The neurohypophysial hormone receptor (NHR) family was classically categorized into four subtypes (V1aR, V1bR, V2R and OTR), while recent studies have identified seven subtypes (V1aR, V1bR, V2aR, V2bR, V2cR, V2dR and OTR; V2aR corresponds to the conventional V2R). The vertebrate NHR family were diversified via multiple gene duplication events at different scales. Despite intensive research effort in non-osteichthyes vertebrates such as cartilaginous fish and lamprey, the molecular phylogeny of the NHR family has not been fully understood. In the present study, we focused on the inshore hagfish (Eptatretus burgeri), another group of cyclostomes, and Arctic lamprey (Lethenteron camtschaticum) for comparison. Two putative NHR homologs, which were previously identified only in silico, were cloned from the hagfish and designated as ebV1R and ebV2R. In vitro, ebV1R, as well as two out of five Arctic lamprey NHRs, increased intracellular Ca2+ in response to exogenous neurohypophysial hormones. None of the examined cyclostome NHRs altered intracellular cAMP levels. Transcripts of ebV1R were detected in multiple tissues including the brain and gill, with intense hybridization signals in the hypothalamus and adenohypophysis, while ebV2R was predominantly expressed in the systemic heart. Similarly, Arctic lamprey NHRs showed distinct expression patterns, underscoring the multifunctionality of VT in the cyclostomes as in the gnathostomes. These results and exhaustive gene synteny comparisons provide new insights into the molecular and functional evolution of the neurohypophysial hormone system in vertebrates.

Absence of some cytochrome P450 (CYP) and hydroxysteroid dehydrogenase (HSD) enzymes in hagfishes

Corticosteroids are synthesized from cholesterol by steroidogenic enzyme catalysts belonging to two main families: the cytochrome p450s (CYPs) and hydroxysteroid dehydrogenases (HSDs). The action of these steroidogenic enzymes allows the genesis of the terminal active corticosteroids 11-deoxycortisol (S), 1ɑ-hydroxycorticosterone (1α-OH-B), or cortisol in different fish species. However, for Cyclostomes like hagfishes, the terminal corticosteroid is still undefined. In this study, we examined the presence or absence of CYPs and HSDs as traits in fishes to gain insight about the primary corticosteroid synthesis pathways of the hagfishes. We used published cytochrome c oxidase I (COXI) amino acid sequences to construct a phylogeny of fishes and then mapped the CYPs and HSDs as morphological traits onto the tree to predict the ancestral character states through ancestral character reconstruction (ACR). There is a clear phylogenetic signal for CYP (i.e., CYP11a1, 17, 21, and 11b) and HSD (i.e., 11-βHSD and 3β-HSD) derivatives of interest throughout the more derived fishes. Using trait-based ACR, we also found that hagfishes possess genes for 3β-HSD, CYP11a1, CYP17, and CYP21. Importantly, the presence of CYP21 implies that hagfish can synthesize 11-deoxycorticosterone (11-DOC) and S. Previous research demonstrated that despite hagfish having CYP21, neither 11-DOC nor S could be detected in hagfish. This discrepancy between the presence of steroidogenic enzymes and products brings into question the expression and/or function of CYP21 in hagfishes.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Novel developmental bases for the evolution of hypobranchial muscles in vertebrates

BACKGROUND

Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic.

RESULTS

Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle.

CONCLUSIONS

We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Molecular cloning and characterization of hagfish estrogen receptors

One or more distinct forms of the nuclear estrogen receptor (ER) have been isolated from many vertebrates to date. To better understand the molecular evolution of ERs, we cloned and characterized er cDNAs from the inshore hagfish, Eptatretus burgeri, a modern representative of the most primitive vertebrates, the agnathans. Two er cDNAs, er1 and er2, were isolated from the liver of a reproductive female hagfish. A phylogenetic analysis placed hagfish ER1 into a position prior to the divergence of vertebrate ERs. Conversely, hagfish ER2 was placed at the base of the vertebrate ERβ clade. The tissue distribution patterns of both ER subtype mRNAs appeared to be different, suggesting that each subtype has different physiological roles associated with estrogen actions. An estrogen responsive-luciferase reporter assay using mammalian HEK293 cells was used to functionally characterize these hagfish ERs. Both ER proteins displayed estrogen-dependent activation of transcription. These results clearly demonstrate that the hagfish has two functional ER subtypes.

RFamide peptides in agnathans and basal chordates

Since a peptide with a C-terminal Arg-Phe-NH2 (RFamide peptide) was first identified in the ganglia of the venus clam in 1977, RFamide peptides have been found in the nervous system of both invertebrates and vertebrates. In vertebrates, the RFamide peptide family includes gonadotropin-inhibitory hormone (GnIH), neuropeptide FF (NPFF), prolactin-releasing peptide (PrRP), pyroglutamylated RFamide peptide/26RFamide peptide (QRFP/26RFa), and kisspeptins (kiss1 and kiss2). They are involved in important functions such as the release of hormones, regulation of sexual or social behavior, pain transmission, reproduction, and feeding. In contrast to tetrapods and jawed fish, the information available on RFamide peptides in agnathans and basal chordates is limited, thus preventing further insights into the evolution of RFamide peptides in vertebrates. In this review, we focus on the previous research and recent advances in the studies on RFamide peptides in agnathans and basal chordates. In agnathans, the genes encoding GnIH, NPFF, and PrRP precursors and the mature peptides have been identified in lamprey (Petromyzon marinus) and hagfish (Paramyxine atami). Putative kiss1 and kiss2 genes have also been found in the genome database of lamprey. In basal chordates, namely, in amphioxus (Branchiostoma japonicum), a common ancestral form of GnIH and NPFF genes and their mature peptides, as well as the ortholog of the QRFP gene have been identified. The studies revealed that the number of orthologs of vertebrate RFamide peptides present in agnathans and basal chordates is greater than expected, suggesting that the vertebrate RFamide peptides might have emerged and expanded at an early stage of chordate evolution.

Bioinformatic Characterization of Genes and Proteins Involved in Blood Clotting in Lampreys

Lampreys and hagfish are the earliest diverging of extant vertebrates and are obvious targets for investigating the origins of complex biochemical systems found in mammals. Currently, the simplest approach for such inquiries is to search for the presence of relevant genes in whole genome sequence (WGS) assemblies. Unhappily, in the past a high-quality complete genome sequence has not been available for either lampreys or hagfish, precluding the possibility of proving gene absence. Recently, improved but still incomplete genome assemblies for two species of lamprey have been posted, and, taken together with an extensive collection of short sequences in the NCBI trace archive, they have made it possible to make reliable counts for specific gene families. Particularly, a multi-source tactic has been used to study the lamprey blood clotting system with regard to the presence and absence of genes known to occur in higher vertebrates. As was suggested in earlier studies, lampreys lack genes for coagulation factors VIII and IX, both of which are critical for the "intrinsic" clotting system and responsible for hemophilia in humans. On the other hand, they have three each of genes for factors VII and X, participants in the "extrinsic" clotting system. The strategy of using raw trace sequence "reads" together with partial WGS assemblies for lampreys can be used in studies on the early evolution of other biochemical systems in vertebrates.

Variable Lymphocyte Receptors: A Current Overview

Jawless vertebrates represented by lampreys and hagfish mount antigen-specific immune responses using variable lymphocyte receptors. These receptors generate diversity comparable to that of T-cell and B-cell receptors by assembling multiple leucine-rich repeat modules with highly variable sequences. Although it is true that jawed and jawless vertebrates have structurally unrelated antigen receptors, their adaptive immune systems have much in common. Most notable is the conservation of lymphocyte lineages. It appears that specialized lymphocyte lineages emerged in a common vertebrate ancestor and that jawed and jawless vertebrates co-opted different antigen receptors within the context of such lymphocyte lineages.

Distribution of mammalian-like melanopsin in cyclostome retinas exhibiting a different extent of visual functions

Mammals contain 1 melanopsin (Opn4) gene that is expressed in a subset of retinal ganglion cells to serve as a photopigment involved in non-image-forming vision such as photoentrainment of circadian rhythms. In contrast, most nonmammalian vertebrates possess multiple melanopsins that are distributed in various types of retinal cells; however, their functions remain unclear. We previously found that the lamprey has only 1 type of mammalian-like melanopsin gene, which is similar to that observed in mammals. Here we investigated the molecular properties and localization of melanopsin in the lamprey and other cyclostome hagfish retinas, which contribute to visual functions including image-forming vision and mainly to non-image-forming vision, respectively. We isolated 1 type of mammalian-like melanopsin cDNA from the eyes of each species. We showed that the recombinant lamprey melanopsin was a blue light-sensitive pigment and that both the lamprey and hagfish melanopsins caused light-dependent increases in calcium ion concentration in cultured cells in a manner that was similar to that observed for mammalian melanopsins. We observed that melanopsin was distributed in several types of retinal cells, including horizontal cells and ganglion cells, in the lamprey retina, despite the existence of only 1 melanopsin gene in the lamprey. In contrast, melanopsin was almost specifically distributed to retinal ganglion cells in the hagfish retina. Furthermore, we found that the melanopsin-expressing horizontal cells connected to the rhodopsin-containing short photoreceptor cells in the lamprey. Taken together, our findings suggest that in cyclostomes, the global distribution of melanopsin in retinal cells might not be related to the melanopsin gene number but to the extent of retinal contribution to visual function.

Characterization of novel RFamide peptides in the central nervous system of the brown hagfish: isolation, localization, and functional analysis

RFamide (RFa) peptides play various important roles in the central nervous system in both invertebrates and vertebrates. However, there is no evidence of the existence of any RFamide peptide in the brain of hagfish, one of the oldest lineages of vertebrates. In this study, we sought to identify novel RFamide peptides from the brains of hagfish (Paramyxine atami). We identified four novel RFamide peptides, which had the C-terminal Pro-Gln-Arg-Phe-NH2 structure. cDNA cloning revealed that the identified RFamide peptides are encoded in two types of cDNA. Molecular phylogenetic analysis of the two precursors indicated that the hagfish RFamide peptides belong to the PQRFamide peptide group that includes mammalian neuropeptide FF and AF. Based on immunohistochemistry and in situ hybridization, hagfish PQRFamide peptide precursor mRNA and its translated peptides were localized in the infundibular nucleus of the hypothalamus. Immunoreactive fibers were terminated on blood vessels in the infundibular nucleus. Dense immunoreactive fibers were also observed in other brain regions. We further showed that one of the hagfish PQRFamide peptides significantly stimulated the expression of gonadotropin-β mRNA in the cultured hagfish pituitary. These results indicate that the control mechanism of gonadotropin expression by a hypothalamic neuropeptide evolved in the agnathan brain. This is the first evidence describing the identification of RFamide peptides in the hagfish brain. This is also the first report showing the regulation of gonadotropin expression by a homolog of neuropeptide FF that belongs to the PQRFamide peptide group in any vertebrate.

Different pressure resistance of lactate dehydrogenases from hagfish is dependent on habitat depth and caused by tetrameric structure dissociation

The effects of high hydrostatic pressure on lactate dehydrogenase (LDH) activities from two species of hagfish were examined. LDH from Eptatretus okinoseanus, a deep-sea species, retained 67% of the original activity even at 100 MPa. LDH activity from Eptatretus burgeri, a shallow-sea species, was completely lost at 50 MPa but recovered to the original value at 0.1 MPa. The tetrameric structure of LDH-A(4) from E. okinoseanus did not change at 50 MPa. In contrast, almost all LDH tetramers from E. burgeri dissociated to dimers and monomers at 50 MPa but reverted to tetramers at 0.1 MPa. These results show that the dissociation of tetramers caused the inactivation of E. burgeri LDH. The difference depends on the number 6 and 10 amino acids. The mechanism of the slight, gradual inactivation of E. okinoseanus LDH at high pressure differs and is probably due to the metamorphosis of its inner structures.

microRNAs reveal the interrelationships of hagfish, lampreys, and gnathostomes and the nature of the ancestral vertebrate

Hagfish and lampreys are the only living representatives of the jawless vertebrates (agnathans), and compared with jawed vertebrates (gnathostomes), they provide insight into the embryology, genomics, and body plan of the ancestral vertebrate. However, this insight has been obscured by controversy over their interrelationships. Morphological cladistic analyses have identified lampreys and gnathostomes as closest relatives, whereas molecular phylogenetic studies recover a monophyletic Cyclostomata (hagfish and lampreys as closest relatives). Here, we show through deep sequencing of small RNA libraries, coupled with genomic surveys, that Cyclostomata is monophyletic: hagfish and lampreys share 4 unique microRNA families, 15 unique paralogues of more primitive microRNA families, and 22 unique substitutions to the mature gene products. Reanalysis of morphological data reveals that support for cyclostome paraphyly was based largely on incorrect character coding, and a revised dataset is not decisive on the mono- vs. paraphyly of cyclostomes. Furthermore, we show fundamental conservation of microRNA expression patterns among lamprey, hagfish, and gnathostome organs, implying that the role of microRNAs within specific organs is coincident with their appearance within the genome and is conserved through time. Together, these data support the monophyly of cyclostomes and suggest that the last common ancestor of all living vertebrates was a more complex organism than conventionally accepted by comparative morphologists and developmental biologists.

Genome biology of the cyclostomes and insights into the evolutionary biology of vertebrate genomes

The jawless vertebrates (lamprey and hagfish) are the closest extant outgroups to all jawed vertebrates (gnathostomes) and can therefore provide critical insight into the evolution and basic biology of vertebrate genomes. As such, it is notable that the genomes of lamprey and hagfish possess a capacity for rearrangement that is beyond anything known from the gnathostomes. Like the jawed vertebrates, lamprey and hagfish undergo rearrangement of adaptive immune receptors. However, the receptors and the mechanisms for rearrangement that are utilized by jawless vertebrates clearly evolved independently of the gnathostome system. Unlike the jawed vertebrates, lamprey and hagfish also undergo extensive programmed rearrangements of the genome during embryonic development. By considering these fascinating genome biologies in the context of proposed (albeit contentious) phylogenetic relationships among lamprey, hagfish, and gnathostomes, we can begin to understand the evolutionary history of the vertebrate genome. Specifically, the deep shared ancestry and rapid divergence of lampreys, hagfish and gnathostomes is considered evidence that the two versions of programmed rearrangement present in lamprey and hagfish (embryonic and immune receptor) were present in an ancestral lineage that existed more than 400 million years ago and perhaps included the ancestor of the jawed vertebrates. Validating this premise will require better characterization of the genome sequence and mechanisms of rearrangement in lamprey and hagfish.

The evolution of early vertebrate photoreceptors

Meeting the challenge of sampling an ancient aquatic landscape by the early vertebrates was crucial to their survival and would establish a retinal bauplan to be used by all subsequent vertebrate descendents. Image-forming eyes were under tremendous selection pressure and the ability to identify suitable prey and detect potential predators was thought to be one of the major drivers of speciation in the Early Cambrian. Based on the fossil record, we know that hagfishes, lampreys, holocephalans, elasmobranchs and lungfishes occupy critical stages in vertebrate evolution, having remained relatively unchanged over hundreds of millions of years. Now using extant representatives of these 'living fossils', we are able to piece together the evolution of vertebrate photoreception. While photoreception in hagfishes appears to be based on light detection and controlling circadian rhythms, rather than image formation, the photoreceptors of lampreys fall into five distinct classes and represent a critical stage in the dichotomy of rods and cones. At least four types of retinal cones sample the visual environment in lampreys mediating photopic (and potentially colour) vision, a sampling strategy retained by lungfishes, some modern teleosts, reptiles and birds. Trichromacy is retained in cartilaginous fishes (at least in batoids and holocephalans), where it is predicted that true scotopic (dim light) vision evolved in the common ancestor of all living gnathostomes. The capacity to discriminate colour and balance the tradeoff between resolution and sensitivity in the early vertebrates was an important driver of eye evolution, where many of the ocular features evolved were retained as vertebrates progressed on to land.

Evolutionary origin and divergence of GnIH and its homologous peptides

Probing undiscovered hypothalamic neuropeptides that play important roles in the regulation of pituitary function in vertebrates is essential for the progress of neuroendocrinology. In 2000, we discovered a novel hypothalamic dodecapeptide inhibiting gonadotropin release in quail and termed it gonadotropin-inhibitory hormone (GnIH). GnIH acts on the pituitary and gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus via a novel G protein-coupled receptor for GnIH to inhibit gonadal development and maintenance by decreasing gonadotropin release and synthesis. Similar findings were observed in other avian species. Thus, GnIH is a key factor controlling avian reproduction. To give our findings a broader perspective, we also found GnIH homologous peptides in the hypothalamus of other vertebrates, such as mammals, reptiles, amphibians and teleosts. GnIH and its homologs share a common C-terminal LPXRFamide (X=L or Q) motif. A mammalian GnIH homolog also inhibited gonadotropin release in mammals like the GnIH action in birds. In contrast to higher vertebrates, hypophysiotropic activities of GnIH homologs were different in lower vertebrates. To clarify the evolutionary origin of GnIH and its homologs, we further sought to identify novel LPXRFamide peptides from the brain of sea lamprey and hagfish, two extant groups of the oldest lineage of vertebrates, Agnatha. In these agnathans, LPXRFamide peptide and its cDNA were identified only from the brain of hagfish. Based on these findings over the past decade, this paper summarizes the evolutionary origin and divergence of GnIH and its homologous peptides.

Evolution of a core gene network for skeletogenesis in chordates

The skeleton is one of the most important features for the reconstruction of vertebrate phylogeny but few data are available to understand its molecular origin. In mammals the Runt genes are central regulators of skeletogenesis. Runx2 was shown to be essential for osteoblast differentiation, tooth development, and bone formation. Both Runx2 and Runx3 are essential for chondrocyte maturation. Furthermore, Runx2 directly regulates Indian hedgehog expression, a master coordinator of skeletal development. To clarify the correlation of Runt gene evolution and the emergence of cartilage and bone in vertebrates, we cloned the Runt genes from hagfish as representative of jawless fish (MgRunxA, MgRunxB) and from dogfish as representative of jawed cartilaginous fish (ScRunx1-3). According to our phylogenetic reconstruction the stem species of chordates harboured a single Runt gene and thereafter Runt locus duplications occurred during early vertebrate evolution. All newly isolated Runt genes were expressed in cartilage according to quantitative PCR. In situ hybridisation confirmed high MgRunxA expression in hard cartilage of hagfish. In dogfish ScRunx2 and ScRunx3 were expressed in embryonal cartilage whereas all three Runt genes were detected in teeth and placoid scales. In cephalochordates (lancelets) Runt, Hedgehog and SoxE were strongly expressed in the gill bars and expression of Runt and Hedgehog was found in endo- as well as ectodermal cells. Furthermore we demonstrate that the lancelet Runt protein binds to Runt binding sites in the lancelet Hedgehog promoter and regulates its activity. Together, these results suggest that Runt and Hedgehog were part of a core gene network for cartilage formation, which was already active in the gill bars of the common ancestor of cephalochordates and vertebrates and diversified after Runt duplications had occurred during vertebrate evolution. The similarities in expression patterns of Runt genes support the view that teeth and placoid scales evolved from a homologous developmental module.

Differential proteome analysis of hagfish dental and somatic skeletal muscles

Hagfish, the plesiomorphic sister group of all vertebrates, are deep-sea scavengers. The large musculus (m.) longitudinalis linguae (dental muscle) is a specialized element of the feeding apparatus that facilitates the efficient ingestion of food. In this article, we compare the protein expression in hagfish dental and somatic (the m. parietalis) skeletal muscles via two-dimensional gel electrophoresis and mass spectrometry in order to characterize the former muscle. Of the 500 proteins screened, 24 were identified with significant differential expression between these muscles. The proteins that were overexpressed in the dental muscle compared to the somatic muscle were troponin C (TnC), glycogen phosphorylase, beta-enolase, fructose-bisphosphate aldolase A (aldolase A), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). In contrast, myosin light chain 1 (MLC 1) and creatine kinase (CK) were over-expressed in the somatic muscle relative to the dental muscle. These results suggest that these two muscles have different energy sources and contractile properties and provide an initial representative map for comparative studies of muscle-protein expression in low craniates.

A degenerate ParaHox gene cluster in a degenerate vertebrate

The ParaHox genes consist of 3 homeobox gene families, Gsx, Xlox, and Cdx, all of which have fundamental roles in development. Xlox (known as IPF1 or PDX1 in vertebrates), for example, is crucial for development of the vertebrate pancreas and is also involved in regulation of insulin expression. The invertebrate amphioxus has a gene cluster containing one gene from each of the gene families, whereas in all vertebrates examined to date there are additional copies resultant from ParaHox gene cluster duplications at the base of the vertebrate lineage. Extant vertebrates basal to bony and cartilaginous fish are central to the question of when and how these multiple genes arose in the vertebrate genome. Here, we report the mapping of a ParaHox gene cluster in 2 species of hagfishes. Unexpectedly, these basal vertebrates have lost a functional Xlox gene from this cluster, unlike every other vertebrate examined to date. Furthermore, our phylogenetic analyses suggest that hagfishes may have diverged from the vertebrate lineage before the duplications, which created the multiple ParaHox clusters in jawed vertebrates.

Time scale for cyclostome evolution inferred with a phylogenetic diagnosis of hagfish and lamprey cDNA sequences

The Cyclostomata consists of the two orders Myxiniformes (hagfishes) and Petromyzoniformes (lampreys), and its monophyly has been unequivocally supported by recent molecular phylogenetic studies. Under this updated vertebrate phylogeny, we performed in silico evolutionary analyses using currently available cDNA sequences of cyclostomes. We first calculated the GC-content at four-fold degenerate sites (GC(4)), which revealed that an extremely high GC-content is shared by all the lamprey species we surveyed, whereas no striking pattern in GC-content was observed in any of the hagfish species surveyed. We then estimated the timing of diversification in cyclostome evolution using nucleotide and amino acid sequences. We obtained divergence times of 470-390 million years ago (Mya) in the Ordovician-Silurian-Devonian Periods for the interordinal split between Myxiniformes and Petromyzoniformes; 90-60 Mya in the Cretaceous-Tertiary Periods for the split between the two hagfish subfamilies, Myxininae and Eptatretinae; 280-220 Mya in the Permian-Triassic Periods for the split between the two lamprey subfamilies, Geotriinae and Petromyzoninae; and 30-10 Mya in the Tertiary Period for the split between the two lamprey genera, Petromyzon and Lethenteron. This evolutionary configuration indicates that Myxiniformes and Petromyzoniformes diverged shortly after the common ancestor of cyclostomes split from the future gnathostome lineage. Our results also suggest that intra-subfamilial diversification in hagfish and lamprey lineages (especially those distributed in the northern hemisphere) occurred in the Cretaceous or Tertiary Periods.

Hagfish embryology with reference to the evolution of the neural crest

Hagfish, which lack both jaws and vertebrae, have long been the subject of intense interest owing to their position at a crucial point in the evolutionary transition to a truly vertebrate body plan. However, unlike the comparatively well characterized vertebrate agnathan lamprey, little is known about hagfish development. The inability to analyse hagfish at early embryonic stages has frustrated attempts to resolve questions with important phylogenetic implications, including fundamental ones relating to the emergence of the neural crest. Here we report the obtainment of multiple pharyngula-stage embryos of the hagfish species Eptatretus burgeri and our preliminary analyses of their early development. We present histological evidence of putative neural crest cells, which appear as delaminated cells that migrate along pathways corresponding to neural crest cells in fish and amphibians. Molecular cloning studies further revealed the expression of several regulatory genes, including cognates of Pax6, Pax3/7, SoxEa and Sox9, suggesting that the hagfish neural crest is specified by molecular mechanisms that are general to vertebrates. We propose that the neural crest emerged as a population of de-epithelialized migratory cells in a common vertebrate ancestor, and suggest that the possibility of classical and molecular embryology in hagfish opens up new approaches to clarifying the evolutionary history of vertebrates.

Two variable lymphocyte receptor genes of the inshore hagfish are located far apart on the same chromosome

Variable lymphocyte receptors (VLR) generate enormous diversity through assembling highly diverse leucine-rich repeat (LRR) modules and presumably function as antigen receptors in jawless vertebrates. The hagfish, which constitute major extant members of jawless vertebrates along with lampreys, have two VLR genes designated VLRA and VLRB, whereas only a single VLR gene has been identified in the lamprey. In the present study, we show by fluorescence in situ hybridization (FISH) that hagfish VLRA and VLRB are located on the same chromosome, but are far apart from each other. Analysis of available inshore hagfish complementary DNA sequences indicates that VLRA and VLRB do not share a LRR module with an identical nucleotide sequence. Physical separation of VLRA and VLRB is consistent with this observation and indicates that the two VLR genes function as separate units. The FISH protocol developed in this study should be useful for the analysis of the agnathan genome.

Eptatretus strickrotti n. sp. (Myxinidae): first hagfish captured from a hydrothermal vent

A single hagfish (Myxinidae, Eptatretus) specimen was recently captured at a hydrothermal vent site on the East Pacific Rise (38 degrees S). This is the first capture of a member of the jawless fishes (agnathans) from a hydrothermal vent site. The specimen differs from all congeners by the very slender body (depth 2.9% of total length), the paired and median ventral nasal sinus papillae, and the presence of 10 afferent branchial arteries on the medial ventral aorta. It is further unique because of a combination of the following features: slime pore counts; paired dorsal nasal sinus papillae; 12 gill pouches and gill apertures; posterior left side of body widely separated from pharyngocutaneous duct; 3/2 multicusp configuration; ventral aorta bifurcated anteriorly between 2nd and 3rd gill pouches (counted from the snout toward the heart); and pink coloration. The specimen is here described as a new species named Eptatretus strickrotti. Molecular 16S rRNA data places this new species as the basal-most species of Eptatretus, providing important new insight to the evolution of hagfishes as a whole.

Why did the acquired immune system of vertebrates evolve?

Rapidly expanding genomic information offers important insights into the mechanisms of acquired immunity and has fostered progress in comparative immunology. However, the question of why the acquired immune system, a trait restricted to jawed vertebrates, evolved has rarely been addressed. Here, I will examine three unique features of early vertebrate evolution, during the 50 million years when the acquired immune system evolved: co-evolution with specialised parasites, increased metabolic rates and genomic instability. These combine to shape a powerful multi-level and multi-causal selective scenario that I propose could have moulded the acquired immune system.

Hagfish and lancelet fibrillar collagens reveal that type II collagen-based cartilage evolved in stem vertebrates

The origin of vertebrates was defined by evolution of a skeleton; however, little is known about the developmental mechanisms responsible for this landmark evolutionary innovation. In jawed vertebrates, cartilage matrix consists predominantly of type II collagen (Col2alpha1), whereas that of jawless fishes has long been thought to be noncollagenous. We recently showed that Col2alpha1 is present in lamprey cartilage, indicating that type II collagen-based cartilage evolved earlier than previously recognized. Here, we investigate the origin of vertebrate cartilage, and we report that hagfishes, the sister group to lampreys, also have Col2alpha1-based cartilage, suggesting its presence in the common ancestor of crown-group vertebrates. We go on to show that lancelets, a sister group to vertebrates, possess an ancestral clade A fibrillar collagen (ColA) gene that is expressed in the notochord. Together, these results suggest that duplication and diversification of ColA genes at the chordate-vertebrate transition may underlie the evolutionary origin of vertebrate skeletal tissues.

Identification of a natriuretic peptide (NP) in cyclostomes (lamprey and hagfish): CNP-4 is the ancestral gene of the NP family

In bony fishes, natriuretic peptides (NPs) comprise a hormone family that is composed of seven subtypes; ANP, BNP, VNP that have an intramolecular ring and N- and C-terminal extensions, and four CNPs (CNP-1 to -4) that lack the C-terminal extension. To assess the ancestral molecule of the NP family, we determined the NP sequences in several species of two extant cyclosotome groups, lampreys and hagfishes. A cDNA encoding CNP was cloned from the heart and brain of three phylogenetically distant species of lampreys, Geotria australis, Lampetra japonica, and Petromyzon marinus. In the deduced prohormone sequence of each species, two potential processing signals, lysine-lysine (KK) that is commonly present in CNP precursors, and arginine-X-X-arginine (RXXR) for furin-like proprotein convertase (PC) that is typical for CNP-4 were present. The deduced mature peptides that are released at each signal were highly conserved among three species; 100% cleaved at KK and >92% processed at RXXR. In L. japonica, the CNP gene was expressed almost exclusively in the heart and brain. Meanwhile, a cDNA encoding NP with a C-terminal tail sequence was cloned from the heart and brain of three hagfish species in different genera, Myxine glutinosa, Eptatretus cirrhatus, and Paramyxine atami. The precursor sequences including the prosegment had >80% identity among the three hagfish species. A processing signal, RXXR, is also conserved in the prosegment of all hagfish NPs. The molecular phylogenetic analyses inferred that the lamprey CNP and hagfish NP belong to the CNP-4 group, even though the hagfish NP has a C-terminal sequence extended from the intramolecular ring. The presence of a processing signal, RXXR, in the prosegment of cyclostome NPs supports the above classification. Based on the current findings, we suggest that the ancestral gene of the NP family is CNP-4.

Variable domains in hagfish: NICIR is a polymorphic multigene family expressed preferentially in leukocytes and is related to lamprey TCR-like

The jawless vertebrates, represented by hagfish and lampreys, are the most advanced animals that apparently lack T cell and B cell receptors. As such, they offer unique opportunities for understanding the evolution of antigen receptors and variable (V)-type immunoglobulin (Ig)-like domains. In the present study, we describe four hagfish Ig superfamily (IgSF) members carrying V-type domains. None of them appeared to have direct counterparts in jawed vertebrates, indicating that many IgSF molecules have either evolved independently in jawed and jawless vertebrates or diverged to the extent that clear homology is no longer recognizable. One of the members encoded a molecule closely related to the previously described membrane protein designated novel ITAM (immunoreceptor tyrosine-based activation motif)-containing IgSF receptors (NICIR). We show here that NICIR is a polymorphic multigene family with at least three members and is expressed predominantly in peripheral blood leukocytes. Phylogenetic analysis indicates that among known proteins, NICIR is most closely related to the lamprey molecule recently proposed to be a potential ancestor of T cell receptors.

Genes coding for intermediate filament proteins closely related to the hagfish "thread keratins (TK)" alpha and gamma also exist in lamprey, teleosts and amphibians

The "thread keratins (TK)" alpha and gamma so far have been considered highly specialized intermediate filament (IF) proteins restricted to hagfish. From lamprey, we now have sequenced five novel IF proteins closely related to TKalpha and TKgamma, respectively. Moreover, we have detected corresponding sequences in EST and genomic databases of teleosts and amphibians. The structure of the TKalpha genes and the positions of their deduced amino acid sequences in a phylogenetic tree clearly support their classification as type II keratins. The genes encoding TKgamma show a structure typical for type III IF proteins, whereas their positions in phylogenetic trees favor a close relationship to the type I keratins. Considering that most keratin-like sequences detected in the lancelet also exhibit a gene structure typical for type III IF proteins, it seems likely that the keratin gene(s) originated from an ancient type III IF protein gene. According to EST analyses, the expression of the thread keratins in teleost fish and amphibians may be particularly restricted to larval stages, which, in conjunction with the observed absence of TKalpha and TKgamma genes in any of the available Amniota databases, indicates a thread keratin function closely related to larval development in an aquatic environment.

cDNA cloning of a mannose-binding lectin-associated serine protease (MASP) gene from hagfish (Eptatretus burgeri)

Hagfish, agnathan cyclostome, is the most primitive extant vertebrate and its complement (C) system seems to be a primordial system in comparison with a well-developed C system in gnathostome vertebrates. From a phylogenic perspective of defense mechanisms, we have isolated complement C3 from the serum of hagfish (Eptatretus burgeri). In this study, we first attempted to identify a hagfish Bf or C2 as a C3 convertase by RT-PCR using degenerative primers designed on the basis of the conserved amino acid stretches among the several kinds of serine proteases. Contrary to our expectation, homology search of cloned RT-PCR product suggested that there was a partial cDNA encoding the homologue of neither Bf nor C2 but a mannose-binding lectin-associated serine protease (MASP). Analyses of a full-length cDNA clone isolated from a hagfish liver cDNA library by using the partial cDNA as a probe indicated that this cDNA encoded hagfish MASP 1. This evidence strongly suggests that the hagfish defends itself against pathogens at least by the complement system composed of lectin pathway.

[Chromatin diminution at the border of the XX and XXI centuries]

The size of genomes in eukaryotic organisms is one of the greatest mysteries of biology. As known from the middle of the XX century, the level of organization of a particular organism, does not depend on its genome size, i. e. on DNA amount in the nucleus. We believe that an actual function of non-coding DNA stands behind the phenomenon of chromatin diminution, known already for 100 years. Diminution of chromatin normally takes place in cells involved in body building and never occurs in developmental precursors of germ cells. Apparently, the former are cells, in which non-coding DNA is functionally significant. We cloned a fraction of DNA eliminated during chromatin diminution of Cyclops kolensis (Cyclopoida, Crustascea) and sequenced 90 clones totally making 32 kb. Taken together, the provided evidence has demonstrated a high organization ordering of DNA sequences restricted to the germ line. Chromatin diminution never takes place in human cells and in cells of the majority of animals. These cells may isolate non-coding DNA in other ways, making it unreactable for most enzymes and thus functionally cut off. Thus, a certain part of genome with a particular size and structure may serve for genetic isolation of species as shellfish or junk DNA are vital components rather than pieces of garbage.

Long-term inheritance of the 28S rDNA-specific retrotransposon R2

R2 is a non-long-terminal-repeat (LTR) retrotransposon that inserts specifically into 28S rDNA. R2 has been identified in many species of arthropods and three species of chordates. R2 may be even more widely distributed in animals, and its origin may be traceable to early animal evolution. In this study, we identified R2 elements in medaka fish, White Cloud Mountain minnow, Reeves' turtle, hagfish, sea lilies, and some arthropod species, using degenerate polymerase chain reaction methods. We also identified two R2 elements from the public genomic sequence database of the bloodfluke Schistosoma mansoni. One of the two bloodfluke R2 elements has two zinc-finger motifs at the N-terminus; this differs from other known R2 elements, which have one or three zinc-finger motifs. Phylogenetic analysis revealed that the whole phylogeny of R2 can be divided into 11 parts (subclades), in which the local R2 phylogeny and the corresponding host phylogeny are consistent. Divergence-versus-age analysis revealed that there is no reliable evidence for the horizontal transfer of R2 but supports the proposition that R2 has been vertically transferred since before the divergence of the deuterostomes and protostomes. The seeming inconsistency between the R2 phylogeny and the phylogeny of their hosts is due to the existence of paralogous lineages. The number of N-terminal zinc-finger motifs is consistent with the deep phylogeny of R2 and indicates that the common ancestor of R2 had three zinc-finger motifs at the N-terminus. This study revealed the long-term vertical inheritance and the ancient origin of sequence specificity of R2, both of which seem applicable to some other non-LTR retrotransposons.

Hagfish leukocytes express a paired receptor family with a variable domain resembling those of antigen receptors

Jawed vertebrates are equipped with TCR and BCR with the capacity to rearrange their V domains. By contrast, jawless vertebrates, represented by hagfish and lampreys, apparently lack such receptors. We describe in this study a family of hagfish genes carrying a single V-type domain resembling those of TCR/BCR. This multigene family, which we call agnathan paired receptors resembling Ag receptors (APAR), is expressed in leukocytes and predicted to encode a group of membrane glycoproteins with organizations characteristic of paired Ig-like receptors, consisting of activating and inhibitory forms. APAR has a J region in its V-type domain, and its V and J regions are encoded in a single exon. Thus, APAR is a member of the emerging families of diversified, innate immune-type receptors with TCR/BCR-like V-type domains and has many of the features expected for a primordial TCR/BCR-like receptor. The extracellular domain of APAR may be descended from a V-type domain postulated to have acquired recombination signal sequences in a jawed vertebrate lineage.

Construction of a bacterial artificial chromosome library from the inshore hagfish, Eptatretus burgeri: A resource for the analysis of the agnathan genome

The jawless fish occupy an important phylogenetic position for understanding the evolution of body plans, the origin of adaptive immunity and genome evolution in chordates. We describe here the construction of a large-insert bacterial artificial chromosome (BAC) library from the inshore hagfish, Eptatretus burgeri. The BAC library contains 93,978 clones with an average insert size of 100 kb and is estimated to represent threefold genome-equivalent coverage. The library was organized in three-dimensional pools to facilitate screening by PCR. We have screened this library by PCR and isolated several BAC clones; the average number of positive clones was compatible with the estimated genome coverage of the library. This BAC library, constructed for the first time from the jawless fish, should serve as a useful resource for the scientific community.

Evidence for independent Hox gene duplications in the hagfish lineage: a PCR-based gene inventory of Eptatretus stoutii

Hox genes code for transcription factors that play a major role in the development of all animal phyla. In invertebrates these genes usually occur as tightly linked cluster, with a few exceptions where the clusters have been dissolved. Only in vertebrates multiple clusters have been demonstrated which arose by duplication from a single ancestral cluster. This history of Hox cluster duplications, in particular during the early elaboration of the vertebrate body plan, is still poorly understood. In this paper we report the results of a PCR survey on genomic DNA of the pacific hagfish Eptatretus stoutii. Hagfishes are one of two clades of recent jawless fishes that are an offshoot of the early radiation of jawless vertebrates. Our data provide evidence for at least 33 distinct Hox genes in the hagfish genome, which is most compatible with the hypothesis of multiple Hox clusters. The largest number, seven, of distinct homeobox fragments could be assigned to paralog group 9, which could imply that the hagfish has more than four clusters. Quartet mapping reveals that within each paralog group the hagfish sequences are statistically more closely related to gnathostome Hox genes than with either amphioxus or lamprey genes. These results support two assumptions about the history of Hox genes: (1) The association of hagfish homeobox sequences with gnathostome sequences suggests that at least one Hox cluster duplication event happened in the stem of vertebrates, i.e., prior to the most recent common ancestor of jawed and jawless vertebrates. (2) The high number of paralog group 9 sequences in hagfish and the phylogenetic position of hagfish suggests that the hagfish lineage underwent additional independent Hox cluster/-gene duplication events.

Transcriptome analysis of hagfish leukocytes: a framework for understanding the immune system of jawless fishes

Jawless fishes occupy a critical phylogenetic position in understanding the origin of the adaptive immune system. Here, we performed large-scale expressed sequence tag analysis of leukocytes isolated from the inshore hagfish Eptatretus burgeri. Although we found many immunity-related genes such as those involved in lymphocyte or hematopoietic cell signaling and development as well as cytokine and cytokine receptor genes, MHC molecules or antigen receptors were not identified. We characterized two hagfish cDNAs that closely resembled mammalian proteins with essential roles in adaptive immunity, one encoding a GATA3-like molecule and another encoding a Bruton's tyrosine kinase (Btk)-like molecule. The GATA3-like gene of hagfish was equidistant from GATA3 and GATA2 in jawed vertebrates. Similarly, the hagfish Btk-like molecule was not Btk itself, but qualified as a pre-duplicated form of Btk and Bmx in jawed vertebrates. In total, our work provides circumstantial evidence that adaptive immunity is unique to jawed vertebrates.

A fast algorithm for determining the best combination of local alignments to a query sequence

Background

Existing sequence alignment algorithms assume that similarities between DNA or amino acid sequences are linearly ordered. That is, stretches of similar nucleotides or amino acids are in the same order in both sequences. Recombination perturbs this order. An algorithm that can reconstruct sequence similarity despite rearrangement would be helpful for reconstructing the evolutionary history of recombined sequences.

Results

We propose a graph-based algorithm for combining multiple local alignments to a query sequence into the single combination of alignments that either covers the maximal portion of the query or results in the single highest alignment score to the query. This algorithm can help study the process of genome rearrangement, improve functional gene annotation, and reconstruct the evolutionary history of recombined proteins. The algorithm takes O(n²) time, where n is the number of local alignments considered.

Conclusions

We discuss two example applications of the algorithm. The algorithm is able to provide useful reconstructions of the metazoan mitochondrial genome. It is also able to increase the percentage of a query sequence's amino acid residues for which similar stretches of amino acids can be found in sequence databases.

Elements of the major myofibrillar binding peptide motif are present in the earliest of true muscle type creatine kinases

Most vertebrates possess two genes for cytoplasmic creatine kinase (CK) coding for muscle (M-CK) and brain (B-CK) isoforms which assemble into homo-dimeric (MM, BB) and hetero-dimeric (MB) active enzymes. In mammals and birds, a significant fraction of MM-CK is bound to the myofibrillar M-line where it is thought to facilitate energy buffering and transport. Myofibrillar binding is mediated by major and minor lysine charge clamp motifs (K104/K115 [major] and K8/K24 [minor] in chicken M-CK) located in the N-terminal region [J. Cell Biol. 149 (2000) 1225]. We have obtained the cDNA and deduced amino acid sequences for cytoplasmic CKs from two hagfish, Myxine glutinosa and Eptatretus stoutii, non-vertebrate craniates, and the sequences for two cytoplasmic CKs from the lamprey Lampetra japonica, a jawless true vertebrate. All four cDNAs code for CKs consisting of approximately 380 residues. Phylogenetic analyses showed that the hagfish and lamprey CKs are coded for by genes which are clearly muscle type (M) creatine kinases. Two of these four M-CKs have the K104/K115-equivalent residues of the major myofibrillar binding region while the other two have the K115 equivalent but lack the corresponding K104 residue. All four M-CKs lack the K8/K24 equivalent elements of the minor myofibrillar binding region. Comparison of these sequences to corresponding sequences of cytoplasmic CKs from two protochordates (tunicate, amphioxus) and M- and B-CKs from true fish and above reveal a pattern of acquisition (and loss) of key lysine residues consistent with the physiological context in which these enzymes operate.

Ikaros family members from the agnathan Myxine glutinosa and the urochordate Oikopleura dioica: emergence of an essential transcription factor for adaptive immunity

The Ikaros multigene family encodes a number of zinc finger transcription factors that play key roles in vertebrate hemopoietic stem cell differentiation and the generation of B, T, and NK cell lineages. In this study, we describe the identification and characterization of an Ikaros family-like (IFL) protein from the agnathan hagfish Myxine glutinosa and the marine urochordate Oikopleura dioica, both of which lie on the evolutionary boundary between the vertebrates and invertebrates. The IFL molecules identified in these animals displayed high conservation in the zinc finger motifs critical for DNA binding and dimerization in comparison with those of jawed vertebrates. Expression of the IFL gene in hagfish was strongest in blood, intestine, and gills. In O. dioica, transcription from the IFL gene was initiated at or around the time of hatching and maintained throughout the life span of the animal. In situ hybridization localized O. dioica IFL expression to the Fol cells, which are responsible for generating the food filter of the house. Biochemical analysis of the DNA binding and dimerization domains from M. glutinosa and O. dioici IFLs showed that M. glutinosa behaves as a true Ikaros family member. Taken together, these results indicate that the properties associated with the Ikaros family preceded the emergence of the jawed vertebrates and thus adaptive immunity.

Hagfish intestinal antimicrobial peptides are ancient cathelicidins

Three potent broad-spectrum antimicrobial peptides (HFIAP-1, -2, and -3) isolated from intestinal tissues of Myxine glutinosa (Atlantic hagfish) are identified as ancient members of the cathelicidin family of antimicrobial peptides, hitherto known only from mammals. In situ hybridization reveals that HFIAPs are produced in nests of myeloid cells within the loose connective tissue of the gut wall, a tissue reminiscent of both gut-associated lymphoid tissue (GALT) and vertebrate spleen. We suggest that this tissue organization provides local defense of the hagfish gastrointestinal tract via innate immunity and possibly served as the architectural plan upon which the adaptive immune system evolved.

Phylogeny of hagfish based on the mitochondrial 16S rRNA gene

The phylogenetic relationships among the species belonging to the family Myxinidae are still debatable. The mitochondrial DNA sequences from the large ribosomal RNA gene may be of great value for systematic and phylogenetic studies within families. Partial sequences of the 16S rRNA gene were obtained for comparisons among the following hagfish species, Paramyxine nelsoni, Paramyxine sheni, Paramyxine taiwanae, Paramyxine yangi, Paramyxine cheni, Eptatretus burgeri, Eptatretus stouii, Eptatretus cirrhatus, Myxine glutinosa, Myxine formosana, Myxine circifrons, Myxine sp1, and Myxine sp2. The boundary of four Paramyxine species (P. sheni, P. taiwanae, P. nelsoni, and P. yangi) from 16S rRNA sequences is ambiguous, however, they are valid based on our unpublished isozyme data as well as the gill aperture arrangement pattern. Both NJ and MP trees constructed from the present molecular data indicate that the genus Paramyxine is diphyletic and Eptatretus paraphyletic. The complexity of Eptatretus and Paramyxine in the clade would not be solved until the farther departed P. cheni is included to form a new clade under the genus Eptatretus. The other clade of Myxininae contains but single genus Myxine.

A single and novel natriuretic peptide is expressed in the heart and brain of the most primitive vertebrate, the hagfish (Eptatretus burgeri)

In teleost fish and tetrapods, the natriuretic peptide (NP) family consists of ANP (atrial natriuretic peptide), BNP (brain natriuretic peptide) and VNP (ventricular natriuretic peptide) that are secreted from the heart, and C-type natriuretic peptide (CNP) that is found in the brain. However, CNP is the only NP identified in the heart and brain of elasmobranchs, suggesting that it is the ancestral type of the NP family and that ANP, BNP and VNP appeared later in the vertebrate phylogeny. To delineate more clearly the molecular evolution of this hormone family, we determined the sequence of NP molecule(s) in evolutionarily the oldest vertebrate group, the cyclostomes. We have cloned a novel NP cDNA from the heart and brain of hagfish, Eptatretus burgeri, using the RACE method and degenerate primers that amplify all known types of NP cDNAs. The novel NP, named EbuNP after the scientific name of this hagfish, appears to be the only NP in the heart and brain, as no other NP cDNAs were amplified even after specific removal of the cloned EbuNP mRNA from the mRNA pool, except for a minor alternatively spliced EbuNP cDNA with a truncated 3'-untranslated sequence. The EbuNP was equally similar to known NPs but was not considered to be a CNP because of the presence of a C-terminal tail sequence. The EbuNP gene was abundantly expressed in the cardiac atrium, ventricle, portal heart and brain but scarcely in the intestine; no expression was observed in the gill and kidney. Mass spectrometry of affinity-purified EbuNP in plasma, heart and brain revealed a 68 amino acid peptide circulating in the blood and stored in the heart, which is cleaved at the typical cleavage signal of a processing enzyme, furin, as observed in mammalian BNP. The C-terminal Gly residue was used for amidation as is the case in eel ANP. The immunoreactive EbuNP was not detected in the brain, suggesting the presence of a different processing form in the brain. These results show that the molecular evolution of the NP family in vertebrates is more complex than previously thought.

Macrophage migration inhibitory factor (MIF) of jawed and jawless fishes: implications for its evolutionary origin

The macrophage migration inhibitory factor (MIF) is a cytokine produced by T lymphocytes and macrophages in response to inflammatory stimuli. We sequenced MIF cDNA clones of two jawless fishes, the sea lamprey (Petromyzon marinus) and the North Atlantic hagfish (Myxine glutinosa), as well as of the jawed (cichlid) fish Paralabidochromis chilotes. The fish MIF-encoding genes have the same exon-intron organization as the mammalian MIF genes and are present in one copy per haploid genome. Secondary and tertiary structure predictions suggest that the fish MIF proteins have a topology characteristic of the entire MIF-family of proteins. Phylogenetic analysis separates the known nematode members of the family into two groups, one having a sister group relationship with the mammalian D-dopachrome tautomerase (DDT) proteins and the other being related to vertebrate MIFs. It also reveals a high degree of convergent evolution among the members of the family. Finally, it suggests that the divergence of MIF and DDT occurred before the emergence of nematodes in metazoan evolution.

Molecular phylogeny of early vertebrates: monophyly of the agnathans as revealed by sequences of 35 genes

Extant vertebrates are divided into three major groups: hagfishes (Hyperotreti, myxinoids), lampreys (Hyperoartia, petromyzontids), and jawed vertebrates (Gnathostomata). The phylogenetic relationships among the groups and within the jawed vertebrates are controversial, for both morphological and molecular studies have rendered themselves to conflicting interpretations. Here, we use the sequences of 35 nuclear protein-encoding genes to provide definitive evidence for the monophyly of the Agnatha (jawless vertebrates, a group encompassing the hagfishes and lampreys). Our analyses also give a strong support for the separation of Chondrichthyes (cartilaginous fishes) before the divergence of Osteichthyes (bony fishes) from the other gnathostomes.

Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.