Molecular phylogeny and divergence times of deuterostome animals

Evolution of gene regulatory network architectures: examples of subcircuit conservation and plasticity between classes of echinoderms

Developmental gene regulatory networks (GRNs) explain how regulatory states are established in particular cells during development and how these states then determine the final form of the embryo. Evolutionary changes to the sequence of the genome will direct reorganization of GRN architectures, which in turn will lead to the alteration of developmental programs. A comparison of GRN architectures must consequently reveal the molecular basis for the evolution of developmental programs among different organisms. This review highlights some of the important findings that have emerged from the most extensive direct comparison of GRN architectures to date. Comparison of the orthologous GRNs for endomesodermal specification in the sea urchin and sea star, provides examples of several discrete, functional GRN subcircuits and shows that they are subject to diverse selective pressures. This demonstrates that different regulatory linkages may be more or less amenable to evolutionary change. One of the more surprising findings from this comparison is that GRN-level functions may be maintained while the factors performing the functions have changed, suggesting that GRNs have a high capacity for compensatory changes involving transcription factor binding to cis regulatory modules.

Expression of multiple Src family kinases in sea urchin eggs and their function in Ca2+ release at fertilization

Egg activation at fertilization in deuterostomes requires a rise in intracellular Ca(2+), which is released from the egg's endoplasmic reticulum. In sea urchins, a Src Family Kinase (SpSFK1) is necessary for the PLCgamma-mediated signaling event that initiates this Ca(2+) release (Giusti, A.F., O'Neill, F.J., Yamasu, K., Foltz, K.R. and Jaffe, L.A., 2003. Function of a sea urchin egg Src family kinase in initiating Ca2+ release at fertilization. Dev. Biol. 256, 367-378.). Annotation of the Strongylocentrotus purpuratus genome sequence led to the identification of additional, predicted SFKs (Bradham, C.A., Foltz, D.R., Beane, W.S., Amone, M.I., Rizzo, F., Coffman, J.A., Mushegian, A., Goel, M., Morales, J., Geneviere, A.M., Lapraz, F., Robertson, A.J., Kelkar, H., Loza-Coll, M., Townley, I.K., Raisch, M., Roux, M.M., Lepage, T., Gache, C., McClay, D.R., Manning, G., 2006. The sea urchin kinome: a first look. Dev. Biol. 300, 180-193.; Roux, M.M., Townley, I.K., Raisch, M., Reade, A., Bradham, C., Humphreys, G., Gunaratne, H.J., Killian, C.E., Moy, G., Su, Y.H., Ettensohn, C.A., Wilt, F., Vacquier, V.D., Burke, R.D., Wessel, G. and Foltz, K.R., 2006. A functional genomic and proteomic perspective of sea urchin calcium signaling and egg activation. Dev. Biol. 300, 416-433.). Here, we describe the cloning and characterization of these 4 additional SFKs and test their function during the initial Ca(2+) release at fertilization using the dominant-interfering microinjection method coupled with Ca(2+) recording. While two of the new SFKs (SpFrk and SpSFK3) are necessary for Ca(2+) release, SpSFK5 appears dispensable for early egg to embryo transition events. Interestingly, SpSFK7 may be involved in preventing precocious release of Ca(2+). Binding studies indicate that only SpSFK1 is capable of direct interaction with PLCgamma. Immunolocalization studies suggest that one or more SpSFK and PLCgamma are localized to the egg cortex and at the site of sperm-egg interaction. Collectively, these data indicate that more than one SFK is involved in the Ca(2+) release pathway at fertilization.

Gnathostome phylogenomics utilizing lungfish EST sequences

The relationship between the Chondrichthyes (cartilaginous fishes), the Actinopterygii (ray-finned fishes), and the piscine Sarcopterygii (lobe-finned fishes) and how the Tetrapoda (four-limbed terrestrial vertebrates) are related to these has been a contentious issue for more than a century. A general consensus about the relationship of these vertebrate clades has gradually emerged among morphologists, but no molecular study has yet provided conclusive evidence for any specific hypothesis. In order to examine these relationships on the basis of more extensive sequence data, we have produced almost 1,000,000 bp of expressed sequence tags (ESTs) from the African marbled lungfish, Protopterus aethiopicus. This new data set yielded 771 transcribed nuclear sequences that had not been previously described. The lungfish EST sequences were combined with EST data from two cartilaginous fishes and whole genome data from an agnathan, four ray-finned fishes, and four tetrapods. Phylogenomic analysis of these data yielded, for the first time, significant maximum likelihood support for a traditional gnathostome tree with a split between the Chondrichthyes and remaining (bone) gnathostomes. Also, the sister group relationship between Dipnoi (lungfishes) and Tetrapoda received conclusive support. Previously proposed hypotheses, such as the monophyly of fishes, could be rejected significantly. The divergence time between lungfishes and tetrapods was estimated to 382-388 Ma by the current data set and six calibration points.

Origin of the genetic components of the vomeronasal system in the common ancestor of all extant vertebrates

Comparative genomics provides a valuable tool for inferring the evolutionary history of physiological systems, particularly when this information is difficult to ascertain by morphological traits. One such example is the vomeronasal system (VNS), a vertebrate nasal chemosensory system that is responsible for detecting intraspecific pheromonal cues as well as environmental odorants. The morphological components of the VNS are found only in tetrapods, but the genetic components of the system have been found in teleost fish, in addition to tetrapods. To determine when the genetic components of the VNS originated, we searched for the VNS-specific genes in the genomes of two early diverging vertebrate lineages: the sea lamprey from jawless fishes and the elephant shark from cartilaginous fishes. Genes encoding vomeronasal type 1 receptors (V1Rs) and Trpc2, two components of the vomeronasal signaling pathway, are present in the sea lamprey genome, and both are expressed in the olfactory organ, revealing that the genetic components of the present-day VNS existed in the common ancestor of all extant vertebrates. Additionally, all three VNS genes, Trpc2, V1Rs, and vomeronasal type 2 receptors (V2Rs), are found in the elephant shark genome. Because V1Rs and V2Rs are related to two families of taste receptors, we also searched the early diverging vertebrate genomes for taste system genes and found them in the shark genome but not in the lamprey. Coupled with known distributions of the genetic components of the vertebrate main olfactory system, our results suggest staggered origins of vertebrate sensory systems. These findings are important for understanding the evolution of vertebrate sensory systems and illustrate the utility of the genome sequences of early diverging vertebrates for uncovering the evolution of vertebrate-specific traits.

Whole-genome duplications in the ancestral vertebrate are detectable in the distribution of gene family sizes of tetrapod species

A clustering of all protein coding genes from the complete genomes of five tetrapod species into gene families shows a clear deviation from the expected power-law distribution of gene family size. We hypothesize that at least part of the deviation is the result of the two whole-genome duplications (WGDs) that are now known, with reasonable certainty, to have occurred prior to the fish-tetrapod split. We build a model of homologous gene family evolution and perform simulations to show that speciations alone cannot produce a distribution that resembles the empirical data. In order to replicate the features of the empirical distribution, the simulation must incorporate two WGD events. These WGDs must be such that a significant number of the gene duplicates generated in the WGDs have a higher retention rate than they do following small-scale duplication (SSD). This requirement is consistent with what is known about duplicate retention following a WGD, namely, that genes belonging to specific functional classes, such as genes regulating transcription, are much more likely to be retained following WGD than SSD. We conclude that the deviation from the power-law that we observe in the empirical data is the result of the two WGDs that occurred in the ancestral chordate. This implies that the two ancient WGDs continue to have a structural effect on gene families approximately 500 million years after the initial events. On the one hand, this is a surprising result, given the limited retention of duplicates generated by a WGD and the continual SSD, which further weakens the signal created by the fraction of duplicate pairs that are retained. On the other hand, WGD's capacity to fundamentally change the architecture of gene families in a profound and lasting way is consistent with the observed correlation between WGDs and important evolutionary transitions.

The quest for orthologs: finding the corresponding gene across genomes

Orthology is a key evolutionary concept in many areas of genomic research. It provides a framework for subjects as diverse as the evolution of genomes, gene functions, cellular networks and functional genome annotation. Although orthologous proteins usually perform equivalent functions in different species, establishing true orthologous relationships requires a phylogenetic approach, which combines both trees and graphs (networks) using reliable species phylogeny and available genomic data from more than two species, and an insight into the processes of molecular evolution. Here, we evaluate the available bioinformatics tools and provide a set of guidelines to aid researchers in choosing the most appropriate tool for any situation.

Evolutionary history of the alpha2,8-sialyltransferase (ST8Sia) gene family: tandem duplications in early deuterostomes explain most of the diversity found in the vertebrate ST8Sia genes

BACKGROUND

The animal sialyltransferases, which catalyze the transfer of sialic acid to the glycan moiety of glycoconjugates, are subdivided into four families: ST3Gal, ST6Gal, ST6GalNAc and ST8Sia, based on acceptor sugar specificity and glycosidic linkage formed. Despite low overall sequence identity between each sialyltransferase family, all sialyltransferases share four conserved peptide motifs (L, S, III and VS) that serve as hallmarks for the identification of the sialyltransferases. Currently, twenty subfamilies have been described in mammals and birds. Examples of the four sialyltransferase families have also been found in invertebrates. Focusing on the ST8Sia family, we investigated the origin of the three groups of alpha2,8-sialyltransferases demonstrated in vertebrates to carry out poly-, oligo- and mono-alpha2,8-sialylation.

RESULTS

We identified in the genome of invertebrate deuterostomes, orthologs to the common ancestor for each of the three vertebrate ST8Sia groups and a set of novel genes named ST8Sia EX, not found in vertebrates. All these ST8Sia sequences share a new conserved family-motif, named "C-term" that is involved in protein folding, via an intramolecular disulfide bridge. Interestingly, sequences from Branchiostoma floridae orthologous to the common ancestor of polysialyltransferases possess a polysialyltransferase domain (PSTD) and those orthologous to the common ancestor of oligosialyltransferases possess a new ST8Sia III-specific motif similar to the PSTD. In osteichthyans, we have identified two new subfamilies. In addition, we describe the expression profile of ST8Sia genes in Danio rerio.

CONCLUSION

Polysialylation appeared early in the deuterostome lineage. The recent release of several deuterostome genome databases and paralogons combined with synteny analysis allowed us to obtain insight into events at the gene level that led to the diversification of the ST8Sia genes, with their corresponding enzymatic activities, in both invertebrates and vertebrates. The initial expansion and subsequent divergence of the ST8Sia genes resulted as a consequence of a series of ancient duplications and translocations in the invertebrate genome long before the emergence of vertebrates. A second subset of ST8sia genes in the vertebrate genome arose from whole genome duplication (WGD) R1 and R2. Subsequent selective ST8Sia gene loss is responsible for the characteristic ST8Sia gene expression pattern observed today in individual species.

The Ediacaran emergence of bilaterians: congruence between the genetic and the geological fossil records

Unravelling the timing of the metazoan radiation is crucial for elucidating the macroevolutionary processes associated with the Cambrian explosion. Because estimates of metazoan divergence times derived from molecular clocks range from quite shallow (Ediacaran) to very deep (Mesoproterozoic), it has been difficult to ascertain whether there is concordance or quite dramatic discordance between the genetic and geological fossil records. Here, we show using a range of molecular clock methods that the major pulse of metazoan divergence times was during the Ediacaran, which is consistent with a synoptic reading of the Ediacaran macrobiota. These estimates are robust to changes in priors, and are returned with or without the inclusion of a palaeontologically derived maximal calibration point. Therefore, the two historical records of life both suggest that although the cradle of Metazoa lies in the Cryogenian, and despite the explosion of ecology that occurs in the Cambrian, it is the emergence of bilaterian taxa in the Ediacaran that sets the tempo and mode of macroevolution for the remainder of geological time.

Deciphering deuterostome phylogeny: molecular, morphological and palaeontological perspectives

Deuterostomes are a monophyletic group of animals that include the vertebrates, invertebrate chordates, ambulacrarians and xenoturbellids. Fossil representatives from most major deuterostome groups, including some phylum-level crown groups, are found in the Lower Cambrian, suggesting that evolutionary divergence occurred in the Late Precambrian, in agreement with some molecular clock estimates. Molecular phylogenies, larval morphology and the adult heart/kidney complex all support echinoderms and hemichordates as a sister grouping (Ambulacraria). Xenoturbellids are a relatively newly discovered phylum of worm-like deuterostomes that lacks a fossil record, but molecular evidence suggests that these animals are a sister group to the Ambulacraria. Within the chordates, cephalochordates share large stretches of chromosomal synteny with the vertebrates, have a complete Hox complex and are sister group to the vertebrates based on ribosomal and mitochondrial gene evidence. In contrast, tunicates have a highly derived adult body plan and are sister group to the vertebrates based on the analyses of concatenated genomic sequences. Cephalochordates and hemichordates share gill slits and an acellular cartilage, suggesting that the ancestral deuterostome also shared these features. Gene network data suggest that the deuterostome ancestor had an anterior-posterior body axis specified by Hox and Wnt genes, a dorsoventral axis specified by a BMP/chordin gradient, and was bilaterally symmetrical with left-right asymmetry determined by expression of nodal.

Do vertebrate neural crest and cranial placodes have a common evolutionary origin?

Two embryonic tissues-the neural crest and the cranial placodes-give rise to most evolutionary novelties of the vertebrate head. These two tissues develop similarly in several respects: they originate from ectoderm at the neural plate border, give rise to migratory cells and develop into multiple cell fates including sensory neurons. These similarities, and the joint appearance of both tissues in the vertebrate lineage, may point to a common evolutionary origin of neural crest and placodes from a specialized population of neural plate border cells. However, a review of the developmental mechanisms underlying the induction, specification, migration and cytodifferentiation of neural crest and placodes reveals fundamental differences between the tissues. Taken together with insights from recent studies in tunicates and amphioxus, this suggests that neural crest and placodes have an independent evolutionary origin and that they evolved from the neural and non-neural side of the neural plate border, respectively.

Novel genes dramatically alter regulatory network topology in amphioxus

Domain rearrangements in the innate immune network of amphioxus suggests that domain shuffling has shaped the evolution of immune systems.

Background

Regulation in protein networks often utilizes specialized domains that 'join' (or 'connect') the network through specific protein-protein interactions. The innate immune system, which provides a first and, in many species, the only line of defense against microbial and viral pathogens, is regulated in this way. Amphioxus (Branchiostoma floridae), whose genome was recently sequenced, occupies a unique position in the evolution of innate immunity, having diverged within the chordate lineage prior to the emergence of the adaptive immune system in vertebrates.

Results

The repertoire of several families of innate immunity proteins is expanded in amphioxus compared to both vertebrates and protostome invertebrates. Part of this expansion consists of genes encoding proteins with unusual domain architectures, which often contain both upstream receptor and downstream activator domains, suggesting a potential role for direct connections (shortcuts) that bypass usual signal transduction pathways.

Conclusion

Domain rearrangements can potentially alter the topology of protein-protein interaction (and regulatory) networks. The extent of such arrangements in the innate immune network of amphioxus suggests that domain shuffling, which is an important mechanism in the evolution of multidomain proteins, has also shaped the development of immune systems.

Dual hindlimb control elements in the Tbx4 gene and region-specific control of bone size in vertebrate limbs

The Tbx4 transcription factor is crucial for normal hindlimb and vascular development, yet little is known about how its highly conserved expression patterns are generated. We have used comparative genomics and functional scanning in transgenic mice to identify a dispersed group of enhancers controlling Tbx4 expression in different tissues. Two independent enhancers control hindlimb expression, one located upstream and one downstream of the Tbx4 coding exons. These two enhancers, hindlimb enhancer A and hindlimb enhancer B (HLEA and HLEB), differ in their primary sequence, in their precise patterns of activity within the hindlimb, and in their degree of sequence conservation across animals. HLEB is highly conserved from fish to mammals. Although Tbx4 expression and hindlimb development occur at different axial levels in fish and mammals, HLEB cloned from either fish or mouse is capable of driving expression at the appropriate position of hindlimb development in mouse embryos. HLEA is highly conserved only in mammals. Deletion of HLEA from the endogenous mouse locus reduces expression of Tbx4 in the hindlimb during embryogenesis, bypasses the embryonic lethality of Tbx4-null mutations, and produces viable, fertile mice with characteristic changes in the size of bones in the hindlimb but not the forelimb. We speculate that dual hindlimb enhancers provide a flexible genomic mechanism for altering the strength and location of Tbx4 expression during normal development, making it possible to separately modify the size of forelimb and hindlimb bones during vertebrate evolution.

The amphioxus genome and the evolution of the chordate karyotype

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

An amphioxus orthologue of the estrogen receptor that does not bind estradiol: insights into estrogen receptor evolution

Background

The origin of nuclear receptors (NRs) and the question whether the ancestral NR was a liganded or an unliganded transcription factor has been recently debated. To obtain insight into the evolution of the ligand binding ability of estrogen receptors (ER), we comparatively characterized the ER from the protochordate amphioxus (Branchiostoma floridae), and the ER from lamprey (Petromyzon marinus), a basal vertebrate.

Results

Extensive phylogenetic studies as well as signature analysis allowed us to confirm that the amphioxus ER (amphiER) and the lamprey ER (lampER) belong to the ER group. LampER behaves as a "classical" vertebrate ER, as it binds to specific DNA Estrogen Responsive Elements (EREs), and is activated by estradiol (E₂), the classical ER natural ligand. In contrast, we found that although amphiER binds EREs, it is unable to bind E₂ and to activate transcription in response to E₂. Among the 7 natural and synthetic ER ligands tested as well as a large repertoire of 14 cholesterol derivatives, only Bisphenol A (an endocrine disruptor with estrogenic activity) bound to amphiER, suggesting that a ligand binding pocket exists within the receptor. Parsimony analysis considering all available ER sequences suggest that the ancestral ER was not able to bind E₂ and that this ability evolved specifically in the vertebrate lineage. This result does not support a previous analysis based on ancestral sequence reconstruction that proposed the ancestral steroid receptor to bind estradiol. We show that biased taxonomic sampling can alter the calculation of ancestral sequence and that the previous result might stem from a high proportion of vertebrate ERs in the dataset used to compute the ancestral sequence.

Conclusion

Taken together, our results highlight the importance of comparative experimental approaches vs ancestral reconstructions for the evolutionary study of endocrine systems: comparative analysis of extant ERs suggests that the ancestral ER did not bind estradiol and that it gained the ability to be regulated by estradiol specifically in the vertebrate lineage, before lamprey split.

Early vertebrate whole genome duplications were predated by a period of intense genome rearrangement

Researchers, supported by data from polyploid plants, have suggested that whole genome duplication (WGD) may induce genomic instability and rearrangement, an idea which could have important implications for vertebrate evolution. Benefiting from the newly released amphioxus genome sequence (Branchiostoma floridae), an invertebrate that researchers have hoped is representative of the ancestral chordate genome, we have used gene proximity conservation to estimate rates of genome rearrangement throughout vertebrates and some of their invertebrate ancestors. We find that, while amphioxus remains the best single source of invertebrate information about the early chordate genome, its genome structure is not particularly well conserved and it cannot be considered a fossilization of the vertebrate preduplication genome. In agreement with previous reports, we identify two WGD events in early vertebrates and another in teleost fish. However, we find that the early vertebrate WGD events were not followed by increased rates of genome rearrangement. Indeed, we measure massive genome rearrangement prior to these WGD events. We propose that the vertebrate WGD events may have been symptoms of a preexisting predisposition toward genomic structural change.

The origin and evolution of the neural crest

Many of the features that distinguish the vertebrates from other chordates are derived from the neural crest, and it has long been argued that the emergence of this multipotent embryonic population was a key innovation underpinning vertebrate evolution. More recently, however, a number of studies have suggested that the evolution of the neural crest was less sudden than previously believed. This has exposed the fact that neural crest, as evidenced by its repertoire of derivative cell types, has evolved through vertebrate evolution. In this light, attempts to derive a typological definition of neural crest, in terms of molecular signatures or networks, are unfounded. We propose a less restrictive, embryological definition of this cell type that facilitates, rather than precludes, investigating the evolution of neural crest. While the evolutionary origin of neural crest has attracted much attention, its subsequent evolution has received almost no attention and yet it is more readily open to experimental investigation and has greater relevance to understanding vertebrate evolution. Finally, we provide a brief outline of how the evolutionary emergence of neural crest potentiality may have proceeded, and how it may be investigated.

NcDNAlign: plausible multiple alignments of non-protein-coding genomic sequences

Genome-wide multiple sequence alignments (MSAs) are a necessary prerequisite for an increasingly diverse collection of comparative genomic approaches. Here we present a versatile method that generates high-quality MSAs for non-protein-coding sequences. The NcDNAlign pipeline combines pairwise BLAST alignments to create initial MSAs, which are then locally improved and trimmed. The program is optimized for speed and hence is particulary well-suited to pilot studies. We demonstrate the practical use of NcDNAlign in three case studies: the search for ncRNAs in gammaproteobacteria and the analysis of conserved noncoding DNA in nematodes and teleost fish, in the latter case focusing on the fate of duplicated ultra-conserved regions. Compared to the currently widely used genome-wide alignment program TBA, our program results in a 20- to 30-fold reduction of CPU time necessary to generate gammaproteobacterial alignments. A showcase application of bacterial ncRNA prediction based on alignments of both algorithms results in similar sensitivity, false discovery rates, and up to 100 putatively novel ncRNA structures. Similar findings hold for our application of NcDNAlign to the identification of ultra-conserved regions in nematodes and teleosts. Both approaches yield conserved sequences of unknown function, result in novel evolutionary insights into conservation patterns among these genomes, and manifest the benefits of an efficient and reliable genome-wide alignment package. The software is available under the GNU Public License at http://www.bioinf.uni-leipzig.de/Software/NcDNAlign/.

Loss of egg yolk genes in mammals and the origin of lactation and placentation

Embryonic development in nonmammalian vertebrates depends entirely on nutritional reserves that are predominantly derived from vitellogenin proteins and stored in egg yolk. Mammals have evolved new resources, such as lactation and placentation, to nourish their developing and early offspring. However, the evolutionary timing and molecular events associated with this major phenotypic transition are not known. By means of sensitive comparative genomics analyses and evolutionary simulations, we here show that the three ancestral vitellogenin-encoding genes were progressively lost during mammalian evolution (until around 30–70 million years ago, Mya) in all but the egg-laying monotremes, which have retained a functional vitellogenin gene. Our analyses also provide evidence that the major milk resource genes, caseins, which have similar functional properties as vitellogenins, appeared in the common mammalian ancestor ∼200–310 Mya. Together, our data are compatible with the hypothesis that the emergence of lactation in the common mammalian ancestor and the development of placentation in eutherian and marsupial mammals allowed for the gradual loss of yolk-dependent nourishment during mammalian evolution.

Egg yolk contains the nutrients required for the development of the nonmammalian vertebrate embryo. These nutrients derive by and large from a single set of proteins, vitellogenins, which are produced in the liver and provide or transport amino acids, lipids, phosphorous, and calcium to the egg. Mammals have evolved new nutritional resources for their developing and early offspring, such as lactation and placentation. However, the evolutionary timing and molecular events associated with this major phenotypic transition are not well understood. In this study, we have investigated the evolutionary fate of the three ancestral vitellogenin-encoding genes in mammals. Using detailed evolutionary analyses of genomes from the three major mammalian lineages (eutherian “placental” mammals, marsupials, and monotremes), we found that these genes progressively lost their functions and became pseudogenes relatively recently during mammalian evolution (the most recent inactivation event occurred roughly 30–70 million years ago). Monotremes, which lactate yet lay small parchment-shelled eggs, even retained a functional vitellogenin gene, consistent with their intermediate reproductive state. Our analyses also provide evidence that the major milk resource genes, caseins, which have similar functional properties as vitellogenins, appeared in the common mammalian ancestor ∼200–310 million years ago. Based on our data, we suggest that the emergence of the alternative resources for the mammalian young—lactation and then placentation—only gradually reduced the need for egg yolk resources (and hence functional vitellogenin genes) in mammals.

The major egg yolk genes, those that express vitellogenins, appear to have progressively lost their functionality during mammalian evolution, probably due to the emergence of the mammalian-specific developmental nourishment resources, lactation, and placentation.

Gene structure evolution of the Na+-Ca2+ exchanger (NCX) family

Background

The Na⁺-Ca²⁺ exchanger (NCX) is an important regulator of cytosolic Ca²⁺ levels. Many of its structural features are highly conserved across a wide range of species. Invertebrates have a single NCX gene, whereas vertebrate species have multiple NCX genes as a result of at least two duplication events. To examine the molecular evolution of NCX genes and understand the role of duplicated genes in the evolution of the vertebrate NCX gene family, we carried out phylogenetic analyses of NCX genes and compared NCX gene structures from sequenced genomes and individual clones.

Results

A single NCX in invertebrates and the protochordate Ciona, and the presence of at least four NCX genes in the genomes of teleosts, an amphibian, and a reptile suggest that a four member gene family arose in a basal vertebrate. Extensive examination of mammalian and avian genomes and synteny analysis argue that NCX4 may be lost in these lineages. Duplicates for NCX1, NCX2, and NCX4 were found in all sequenced teleost genomes. The presence of seven genes encoding NCX homologs may provide teleosts with the functional specialization analogous to the alternate splicing strategy seen with the three NCX mammalian homologs.

Conclusion

We have demonstrated that NCX4 is present in teleost, amphibian and reptilian species but has been secondarily and independently lost in mammals and birds. Comparative studies on conserved vertebrate homologs have provided a possible evolutionary route taken by gene duplicates subfunctionalization by minimizing homolog number.

The evolution of core proteins involved in microRNA biogenesis

BACKGROUND

MicroRNAs (miRNAs) are a recently discovered class of non-coding RNAs (ncRNAs) which play important roles in eukaryotic gene regulation. miRNA biogenesis and activation is a complex process involving multiple protein catalysts and involves the large macromolecular RNAi Silencing Complex or RISC. While phylogenetic analyses of miRNA genes have been previously published, the evolution of miRNA biogenesis itself has been little studied. In order to better understand the origin of miRNA processing in animals and plants, we determined the phyletic occurrences and evolutionary relationships of four major miRNA pathway protein components; Dicer, Argonaute, RISC RNA-binding proteins, and Exportin-5.

RESULTS

Phylogenetic analyses show that all four miRNA pathway proteins were derived from large multiple protein families. As an example, vertebrate and invertebrate Argonaute (Ago) proteins diverged from a larger family of PIWI/Argonaute proteins found throughout eukaryotes. Further gene duplications among vertebrates after the evolution of chordates from urochordates but prior to the emergence of fishes lead to the evolution of four Ago paralogues. Invertebrate RISC RNA-binding proteins R2D2 and Loquacious are related to other RNA-binding protein families such as Staufens as well as vertebrate-specific TAR (HIV trans-activator RNA) RNA-binding protein (TRBP) and protein kinase R-activating protein (PACT). Export of small RNAs from the nucleus, including miRNA, is facilitated by three closely related karyopherin-related nuclear transporters, Exportin-5, Exportin-1 and Exportin-T. While all three exportins have direct orthologues in deutrostomes, missing exportins in arthropods (Exportin-T) and nematodes (Exportin-5) are likely compensated by dual specificities of one of the other exportin paralogues.

CONCLUSION

Co-opting particular isoforms from large, diverse protein families seems to be a common theme in the evolution of miRNA biogenesis. Human miRNA biogenesis proteins have direct, orthologues in cold-blooded fishes and, in some cases, urochordates and deutrostomes. However, lineage specific expansions of Dicer in plants and invertebrates as well as Argonaute and RNA-binding proteins in vertebrates suggests that novel ncRNA regulatory mechanisms can evolve in relatively short evolutionary timeframes. The occurrence of multiple homologues to RNA-binding and Argonaute/PIWI proteins also suggests the possible existence of further pathways for additional types of ncRNAs.

Multiple waves of recent DNA transposon activity in the bat, Myotis lucifugus

DNA transposons, or class 2 transposable elements, have successfully propagated in a wide variety of genomes. However, it is widely believed that DNA transposon activity has ceased in mammalian genomes for at least the last 40 million years. We recently reported evidence for the relatively recent activity of hAT and Helitron elements, two distinct groups of DNA transposons, in the lineage of the vespertilionid bat Myotis lucifugus. Here, we describe seven additional families that have also been recently active in the bat lineage. Early vespertilionid genome evolution was dominated by the activity of Helitrons, mariner-like and Tc2-like elements. This was followed by the colonization of Tc1-like elements, and by a more recent explosion of hAT-like elements. Finally, and most recently, piggyBac-like elements have amplified within the Myotis genome and our results indicate that one of these families is probably still expanding in natural populations. Together, these data suggest that there has been tremendous recent activity of various DNA transposons in the bat lineage that far exceeds those previously reported for any mammalian lineage. The diverse and recent populations of DNA transposons in genus Myotis will provide an unprecedented opportunity to study the impact of this class of elements on mammalian genome evolution and to better understand what makes some species more susceptible to invasion by genomic parasites than others.

The search for non-chordate retinoic acid signaling: lessons from chordates

Signaling by retinoic acid (RA) is an important pathway in the development and homeostasis of vertebrate and invertebrate chordates, with a critical role in mesoderm patterning. Classical studies on the distribution of nuclear receptors of animals suggested that the family of RA receptors (RARs/NR1B) was restricted to chordates, while the family of RA X receptors (RXR/NR2B) was distributed from cnidarians to chordates. However, the accumulation of data from genome projects and studies in non-model species is questioning this traditional view. Here we discuss the evidence for non-chordate RA signaling systems in the light of recent advances in our understanding of carotene (pro-Vitamin A) metabolism and of the identification of potential RARs and members of the NR1 family in echinoderms and lophotrochozoan trematodes, respectively. We conclude, as have others before (Bertrand et al., 2004. Mol Biol Evol 21(10):1923-1937), that signaling by RA is more likely an ancestral feature of bilaterians than a chordate innovation.

The human phylome

The human phylome, which includes evolutionary relationships of all human proteins and their homologs among thirty-nine fully sequenced eukaryotes, is reconstructed.

Background:

Phylogenomics analyses serve to establish evolutionary relationships among organisms and their genes. A phylome, the complete collection of all gene phylogenies in a genome, constitutes a valuable source of information, but its use in large genomes still constitutes a technical challenge. The use of phylomes also requires the development of new methods that help us to interpret them.

Results:

We reconstruct here the human phylome, which includes the evolutionary relationships of all human proteins and their homologs among 39 fully sequenced eukaryotes. Phylogenetic techniques used include alignment trimming, branch length optimization, evolutionary model testing and maximum likelihood and Bayesian methods. Although differences with alternative topologies are minor, most of the trees support the Coelomata and Unikont hypotheses as well as the grouping of primates with laurasatheria to the exclusion of rodents. We assess the extent of gene duplication events and their relationship with the functional roles of the protein families involved. We find support for at least one, and probably two, rounds of whole genome duplications before vertebrate radiation. Using a novel algorithm that is independent from a species phylogeny, we derive orthology and paralogy relationships of human proteins among eukaryotic genomes.

Conclusion:

Topological variations among phylogenies for different genes are to be expected, highlighting the danger of gene-sampling effects in phylogenomic analyses. Several links can be established between the functions of gene families duplicated at certain phylogenetic splits and major evolutionary transitions in those lineages. The pipeline implemented here can be easily adapted for use in other organisms.

Development and evolution of chordate cartilage

Deuterostomes are a monophyletic group of animals containing vertebrates, lancelets, tunicates, hemichordates, echinoderms, and xenoturbellids. Four out of these six extant groups-vertebrates, lancelets, tunicates, and hemichordates-have pharyngeal gill slits. All groups of deuterostome animals that have pharyngeal gill slits also have a pharyngeal skeleton supporting the pharyngeal openings, except tunicates. We previously found that pharyngeal cartilage in hemichordates and cephalochordates contains a fibrillar collagen protein similar to vertebrate type II collagen, but unlike vertebrate cartilage, the invertebrate deuterostome cartilages are acellular. We found SoxE and fibrillar collagen expression in the pharyngeal endodermal cells adjacent to where the cartilages form. These same endodermal epithelial cells also express Pax1/9, a marker of pharyngeal endoderm in vertebrates, lancelets, tunicates, and hemichordates. In situ experiments with a cephalochordate fibrillar collagen also showed expression in pharyngeal endoderm, as well as the ectoderm and the mesodermal coelomic pouches lining the gill bars. These results indicate that the pharyngeal endodermal cells are responsible for secretion of the cartilage in hemichordates, whereas in lancelets, all the pharyngeal cells surrounding the gill bars, ectodermal, endodermal, and mesodermal may be responsible for cartilage formation. We propose that endoderm secretion was primarily the ancestral mode of making pharyngeal cartilages in deuterostomes. Later the evolutionary origin of neural crest allowed co-option of the gene network for the secretion of pharyngeal cartilage matrix in the new migratory neural crest cell populations found in vertebrates.

Molecular evolution of fibrillar collagen in chordates, with implications for the evolution of vertebrate skeletons and chordate phylogeny

Vertebrates have seven types of fibrillar collagens that are encoded by 11 genes. Types I, V, and XXIV collagens are components of mineralized bone, whereas types II, XI, and XXVII collagens are components of cartilage. In this study, we traced the molecular evolutionary history of chordate collagen genes and examined how gene duplications gave rise to the collagen genes used for skeletons. Our analyses of deuterostome collagen genes, including one amphioxus gene that we identified in this study, suggest that the common ancestors of deuterostomes possessed three fibrillar collagen genes. Expression analyses of chordate fibrillar collagen genes suggest that in the ancestors of chordates, fibrillar collagen was co-opted to the formation of the notochord sheath independently in three clades. Our results also imply that co-option of collagen genes to cartilage occurred in clade A (col2A1), clade B (col11A1, 11A2), and clade C (COL27A1). Similarly, some fibrillar collagen genes have been co-opted for mineralized bone independently from clade A genes (col1A1, 1A2, 5A2), clade B genes (col5A1), and clade C genes (COL24A1). These frequent co-options for notochord, cartilage, and mineralized bone must have been accompanied by the rapid evolution of cis-regulatory elements for transcription. In addition, we found that one of the ascidian fibrillar collagen genes possesses an amino acid insertion at the identical site of the C-terminal noncollagenous domain in vertebrate fibrillar collagen genes. This observation raises a suspicion about the relatively well-accepted phylogeny of the close relationship between amphioxus and vertebrates.

Evolution of the vertebrate eye: opsins, photoreceptors, retina and eye cup

Charles Darwin appreciated the conceptual difficulty in accepting that an organ as wonderful as the vertebrate eye could have evolved through natural selection. He reasoned that if appropriate gradations could be found that were useful to the animal and were inherited, then the apparent difficulty would be overcome. Here, we review a wide range of findings that capture glimpses of the gradations that appear to have occurred during eye evolution, and provide a scenario for the unseen steps that have led to the emergence of the vertebrate eye.

A conserved role for FGF signaling in chordate otic/atrial placode formation

The widely held view that neurogenic placodes are vertebrate novelties has been challenged by morphological and molecular data from tunicates suggesting that placodes predate the vertebrate divergence. Here, we examine requirements for the development of the tunicate atrial siphon primordium, thought to share homology with the vertebrate otic placode. In vertebrates, FGF signaling is required for otic placode induction and for later events following placode invagination, including elaboration and patterning of the inner ear. We show that results from perturbation of the FGF pathway in the ascidian Ciona support a similar role for this pathway: inhibition with MEK or Fgfr inhibitor at tailbud stages in Ciona results in a larva which fails to form atrial placodes; inhibition during metamorphosis disrupts development of the atrial siphon and gill slits, structures which form where invaginated atrial siphon ectoderm apposes pharyngeal endoderm. We show that laser ablation of atrial primordium ectoderm also results in a failure to form gill slits in the underlying endoderm. Our data suggest interactions required for formation of the atrial siphon and highlight the role of atrial ectoderm during gill slit morphogenesis.

Characterization, developmental expression and evolutionary features of the huntingtin gene in the amphioxus Branchiostoma floridae

BACKGROUND

Huntington's disease is an inherited neurodegenerative disorder that is caused by the expansion of an N-terminal polyQ stretch in the huntingtin protein. In order to investigate the hypothesis that huntingtin was already involved in development of the nervous system in the last common ancestor of chordates, we isolated and characterised the huntingtin homologue from the amphioxus Branchiostoma floridae. In the present paper the amphioxus general term must be referred to Branchiostoma floridae.

RESULTS

In this report, we show that the exon-intron organization of the amphioxus huntingtin gene is highly conserved with that of other vertebrates species. The AmphiHtt protein has two glutamine residues in the position of the typical vertebrate polyQ tract. Sequence conservation is greater along the entire length of the protein than in a previously identified Ciona huntingtin. The first three N-terminal HEAT repeats are highly conserved in vertebrates and amphioxus, although exon rearrangement has occurred in this region. AmphiHtt expression is detectable by in situ hybridization starting from the early neurula stage, where it is found in cells of the neural plate. At later stages, it is retained in the neural compartment but also it appears in limited and well-defined groups of non-neural cells. At subsequent larval stages, AmphiHtt expression is detected in the neural tube, with the strongest signal being present in the most anterior part.

CONCLUSION

The cloning of amphioxus huntingtin allows to infer that the polyQ in huntingtin was already present 540 million years ago and provides a further element for the study of huntingtin function and its evolution along the deuterostome branch.

Unifying themes in host defence effector polypeptides

It is said that nature is the greatest innovator, yet molecular conservation can be equally powerful. One key requirement for the survival of any host is its ability to defend against infection, predation and competition. Recent discoveries, including the presence of a multidimensional structural signature, have revealed a previously unforeseen structural and functional congruence among host defence effector molecules spanning all kingdoms of life. Antimicrobial peptides, kinocidins, polypeptide venoms and other molecules that were once thought to be distinct in form and function now appear to be members of an ancient family of host defence effectors. These molecules probably descended from archetype predecessors that emerged during the beginning of life on earth. Understanding how nature has sustained these host defence molecules with a potent efficacy in the face of dynamic microbial evolution should provide new opportunities to prevent or treat life-threatening infections.

Fossils Impact as Hard as Living Taxa in Parsimony Analyses of Morphology

Systematists disagree whether data from fossils should be included in parsimony analyses. In a handful of well-documented cases, the addition of fossil data radically overturns a hypothesis of relationships based on extant taxa alone. Fossils can break up long branches and preserve character combinations closer in time to deep splitting events. However, fossils usually require more interpretation than extant taxa, introducing greater potential for spurious codings. Moreover, because fossils often have more "missing" codings, they are frequently accused of increasing numbers of MPTs, frustrating resolution and reducing support. Despite the controversy, remarkably little is known about the effects of fossils more generally. Here we provide the first systematic study, investigating empirically the behavior of fossil and extant taxa in 45 published morphological data sets. First-order jackknifing is used to determine the effects that each terminal has on inferred relationships, on the number of MPTs, and on CI' and RI as measures of homoplasy. Bootstrap leaf stabilities provide a proxy for the contribution of individual taxa to the branch support in the rest of the tree. There is no significant difference in the impact of fossil versus extant taxa on relationships, numbers of MPTs, and CI' or RI. However, adding individual fossil taxa is more likely to reduce the total branch support of the tree than adding extant taxa. This must be weighed against the superior taxon sampling afforded by including judiciously coded fossils, providing data from otherwise unsampled regions of the tree. We therefore recommend that investigators should include fossils, in the absence of compelling and case specific reasons for their exclusion.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications

Background

The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes.

Results

We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons.

Conclusion

There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular – but possibly clusters of genes more generally – might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

The amphioxus SoxB family: implications for the evolution of vertebrate placodes

Cranial placodes are regions of thickened ectoderm that give rise to sense organs and ganglia in the vertebrate head. Homologous structures are proposed to exist in urochordates, but have not been found in cephalochordates, suggesting the first chordates lacked placodes. SoxB genes are expressed in discrete subsets of vertebrate placodes. To investigate how placodes arose and diversified in the vertebrate lineage we isolated the complete set of SoxB genes from amphioxus and analyzed their expression in embryos and larvae. We find that while amphioxus possesses a single SoxB2 gene, it has three SoxB1 paralogs. Like vertebrate SoxB1 genes, one of these paralogs is expressed in non-neural ectoderm destined to give rise to sensory cells. When considered in the context of other amphioxus placode marker orthologs, amphioxus SoxB1 expression suggests a diversity of sensory cell types utilizing distinct placode-type gene programs was present in the first chordates. Our data supports a model for placode evolution and diversification whereby the full complement of vertebrate placodes evolved by serial recruitment of distinct sensory cell specification programs to anterior pre-placodal ectoderm.

Chordate ancestry of the neural crest: New insights from ascidians

This article reviews new insights from ascidians on the ancestry of vertebrate neural crest (NC) cells. Ascidians have neural crest-like cells (NCLC), which migrate from the dorsal midline, express some of the typical NC markers, and develop into body pigment cells. These characters suggest that primordial NC cells were already present in the common ancestor of the vertebrates and urochordates, which have been recently inferred as sister groups. The primitive role of NCLC may have been in pigment cell dispersal and development. Later, additional functions may have appeared in the vertebrate lineage, resulting in the evolution of definitive NC cells.

Cis-regulation of the amphioxus engrailed gene: Insights into evolution of a muscle-specific enhancer

To gain insights into the relation between evolution of cis-regulatory DNA and evolution of gene function, we identified tissue-specific enhancers of the engrailed gene of the basal chordate amphioxus (Branchiostoma floridae) and compared their ability to direct expression in both amphioxus and its nearest chordate relative, the tunicate Ciona intestinalis. In amphioxus embryos, the native engrailed gene is expressed in three domains - the eight most anterior somites, a few cells in the central nervous system (CNS) and a few ectodermal cells. In contrast, in C. intestinalis, in which muscle development is highly divergent, engrailed expression is limited to the CNS. To characterize the tissue-specific enhancers of amphioxus engrailed, we first showed that 7.8kb of upstream DNA of amphioxus engrailed directs expression to all three domains in amphioxus that express the native gene. We then identified the amphioxus engrailed muscle-specific enhancer as the 1.2kb region of upstream DNA with the highest sequence identity to the mouse en-2 jaw muscle enhancer. This amphioxus enhancer directed expression to both the somites in amphioxus and to the larval muscles in C. intestinalis. These results show that even though expression of the native engrailed has apparently been lost in developing C. intestinalis muscles, they express the transcription factors necessary to activate transcription from the amphioxus engrailed enhancer, suggesting that gene networks may not be completely disrupted if an individual component is lost.

Phylogenetic position of a whale-fall lancelet (Cephalochordata) inferred from whole mitochondrial genome sequences

Background

The lancelet Asymmetron inferum (subphylum Cephalochordata) was recently discovered on the ocean floor off the southwest coast of Japan at a depth of 229 m, in an anaerobic and sulfide-rich environment caused by decomposing bodies of the sperm whale Physeter macrocephalus. This deep sulfide-rich habitat of A. inferum is unique among the lancelets. The distinguishing adaptation of this species to such an extraordinary habitat can be considered in a phylogenetic framework. As the first step of reconstruction of the evolutionary processes in this species, we investigated its phylogenetic position based on 11 whole mitochondrial genome sequences including the newly determined ones of the whale-fall lancelet A. inferum and two coral-reef congeners.

Results

Our phylogenetic analyses showed that extant lancelets are clustered into two major clades, the Asymmetron clade and the Epigonichthys + Branchiostoma clade. A. inferum was in the former and placed in the sister group to A. lucayanum complex. The divergence time between A. inferum and A. lucayanum complex was estimated to be 115 Mya using the penalized likelihood (PL) method or 97 Mya using the nonparametric rate smoothing (NPRS) method (the middle Cretaceous). These are far older than the first appearance of large whales (the middle Eocene, 40 Mya). We also discovered that A. inferum mitogenome (mitochondrial genome) has been subjected to large-scale gene rearrangements, one feature of rearrangements being unique among the lancelets and two features shared with A. lucayanum complex.

Conclusion

Our study supports the monophyly of genus Asymmetron assumed on the basis of the morphological characters. Furthermore, the features of the A. inferum mitogenome expand our knowledge of variation within cephalochordate mitogenomes, adding a new case of transposition and inversion of the trnQ gene. Our divergence time estimation suggests that A. inferum remained a member of the Mesozoic and the early Cenozoic large vertebrate-fall communities before shifting to become a whale-fall specialist.

Evolution of Class I cytokine receptors

BACKGROUND

The Class I cytokine receptors have a wide range of actions, including a major role in the development and function of immune and blood cells. However, the evolution of the genes encoding them remains poorly understood. To address this we have used bioinformatics to analyze the Class I receptor repertoire in sea squirt (Ciona intestinalis) and zebrafish (Danio rerio).

RESULTS

Only two Class I receptors were identified in sea squirt, one with homology to the archetypal GP130 receptor, and the other with high conservation with the divergent orphan receptor CLF-3. In contrast, 36 Class I cytokine receptors were present in zebrafish, including representative members for each of the five structural groups found in mammals. This allowed the identification of 27 core receptors belonging to the last common ancestor of teleosts and mammals.

CONCLUSION

This study suggests that the majority of diversification of this receptor family occurred after the divergence of urochordates and vertebrates approximately 794 million years ago (MYA), but before the divergence of ray-finned from lobe-finned fishes around 476 MYA. Since then, only relatively limited lineage-specific diversification within the different Class I receptor structural groups has occurred.

Ribosomal RNA genes and deuterostome phylogeny revisited: More cyclostomes, elasmobranchs, reptiles, and a brittle star

This is an expanded study of the relationships among the deuterostome animals based on combined, nearly complete 28S and 18S rRNA genes (>3925 nt.). It adds sequences from 20 more taxa to the approximately 45 sequences used in past studies. Seven of the new taxa were sequenced here (brittle star Ophiomyxa, lizard Anolis, turtle Chrysemys, sixgill shark Hexanchus, electric ray Narcine, Southern Hemisphere lamprey Geotria, and Atlantic hagfish Myxine for 28S), and the other 13 were from GenBank and the literature (from a chicken, dog, rat, human, three lungfishes, and several ray-finned fishes, or Actinopterygii). As before, our alignments were based on secondary structure but did not account for base pairing in the stems of rRNA. The new findings, derived from likelihood-based tree-reconstruction methods and by testing hypotheses with parametric bootstrapping, include: (1) brittle star joins with sea star in the echinoderm clade, Asterozoa; (2) with two hagfishes and two lampreys now available, the cyclostome (jawless) fishes remain monophyletic; (3) Hexanchiform sharks are monophyletic, as Hexanchus groups with the frilled shark, Chlamydoselachus; (4) turtle is the sister taxon of all other amniotes; (5) bird is closer to the lizard than to the mammals; (6) the bichir Polypterus is in a monophyletic Actinopterygii; (7) Zebrafish Danio is the sister taxon of the other two teleosts we examined (trout and perch); (8) the South American and African lungfishes group together to the exclusion of the Australian lungfish. Other findings either upheld those of the previous rRNA-based studies (e.g., echinoderms and hemichordates group as Ambulacraria; orbitostylic sharks; batoids are not derived from any living lineage of sharks) or were obvious (monophyly of mammals, gnathostomes, vertebrates, echinoderms, etc.). Despite all these findings, the rRNA data still fail to resolve the relations among the major groups of deuterostomes (tunicates, Ambulacraria, cephalochordates and vertebrates) and of gnathostomes (chondrichthyans, lungfishes, coelacanth, actinopterygians, amphibians, and amniotes), partly because tunicates and lungfishes are rogue taxa that disrupt the tree. Nonetheless, parametric bootstrapping showed our RNA-gene data are only consistent with these dominant hypotheses: (1) deuterostomes consist of Ambulacraria plus Chordata, with Chordata consisting of tunicates and 'vertebrates plus cephalochordates'; and (2) lungfishes are the closest living relatives of tetrapods.