Molecular phylogeny and divergence times of deuterostome animals

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.

The phylogeny and evolution of host choice in the Hippoboscoidea (Diptera) as reconstructed using four molecular markers

Hippoboscoidea is a superfamily of Diptera that contains the Glossinidae or tsetse flies, the Hippoboscidae or louse flies, and two families of bat flies, the Streblidae and the Nycteribiidae. We reconstruct the phylogenetic relationships within Hippoboscoidea using maximum parsimony and Bayesian methods based on nucleotide sequences from fragments of four genes: nuclear 28S ribosomal DNA and the CPSase domain of CAD, and mitochondrial 16S rDNA and cytochrome oxidase I. We recover monophyly for most of the presently recognized groups within Hippoboscoidea including the superfamily as a whole, the Hippoboscidae, the Nycteribiidae, the bat flies, and the Pupipara (=Hippoboscidae+Nycteribiidae+Streblidae), as well as several subfamilies within the constituent families. Streblidae appear to be paraphyletic. Our phylogenetic hypothesis is well supported and decisive in that most competing topological hypotheses for the Hippoboscoidea require significantly longer trees. We confirm a single shift from a free-living fly to a blood-feeding ectoparasite of vertebrates and demonstrate that at least two host shifts from mammals to birds have occurred. Wings have been repeatedly lost, but never regained. The hippoboscoid ancestor also evolved adenotrophic viviparity and our cladogram is consistent with a gradual reduction in the motility of the deposited final instar larvae from active burrowing in the soil to true pupiparity where adult females glue the puparium within the confines of bat roosts.

A complement response may activate metamorphosis in the ascidian Boltenia villosa

Ascidian metamorphosis transforms a free-swimming larval chordate ascidian into a sessile adult through a distinct series of metamorphic events. Initially, larvae must become competent to respond to settlement cues. Settlement is then marked by dramatic body plan remodeling and may be accompanied by attachment to the substrate. Subtractive hybridization has revealed that many innate immunity transcripts are upregulated during metamorphosis in the ascidian Boltenia villosa. Several of these genes have well-known roles in the mannose-binding lectin (MBL)-complement pathway of innate immunity, including MBL and MBL-activated serine protease (MASP). MBL recognizes and binds to bacterial pathogens, activates MASP, and triggers the complement cascade. In B. villosa, larvae upregulate BvMASP at the time of competency to initiate settlement. We show that several bacterial strains can induce settlement and that the timing of BvMASP expression in the papillae-associated tissue (PAT) cells is tightly correlated with larval competency. We further demonstrate that serine protease inhibitors used to block the complement response also block metamorphosis, allowing tail resorption, but preventing further morphological changes. Based on these experiments, we propose that the MBL-complement pathway may be important for competency, bacterial substrate detection and body plan remodeling during metamorphosis.

Development of a chordate anterior–posterior axis without classical retinoic acid signaling

Developmental signaling by retinoic acid (RA) is thought to be an innovation essential for the origin of the chordate body plan. The larvacean urochordate Oikopleura dioica maintains a chordate body plan throughout life, and yet its genome appears to lack genes for RA synthesis, degradation, and reception. This suggests the hypothesis that the RA-machinery was lost during larvacean evolution, and predicts that Oikopleura development has become independent of RA-signaling. This prediction raises the problem that the anterior-posterior organization of a chordate body plan can be developed without the classical morphogenetic role of RA. To address this problem, we performed pharmacological treatments and analyses of developmental molecular markers to investigate whether RA acts in anterior-posterior axial patterning in Oikopleura embryos. Results revealed that RA does not cause homeotic posteriorization in Oikopleura as it does in vertebrates and cephalochordates, and showed that a chordate can develop the phylotypic body plan in the absence of the classical morphogenetic role of RA. A comparison of Oikopleura and ascidian evidence suggests that the lack of RA-induced homeotic posteriorization is a shared derived feature of urochordates. We discuss possible relationships of altered roles of RA in urochordate development to genomic events, such as rupture of the Hox-cluster, in the context of a new understanding of chordate phylogeny.

Divergent evolution of the vertebrate polysialyltransferase Stx and Pst genes revealed by fish-to-mammal comparison

Polysialic acid (PSA) is a developmentally regulated carbohydrate attached to the neural cell adhesion molecule (NCAM). PSA is involved in dynamic processes like cell migration, neurite outgrowth and neuronal plasticity. In mammals, polysialylation of NCAM is catalyzed independently by two polysialyltransferases, STX (ST8Sia II) and PST (ST8Sia IV), with STX mainly acting during early development and PST at later stages and into adulthood. Here, we functionally characterize zebrafish Stx and Pst homolog genes during fish development and evaluate their catalytic affinity for NCAM in vitro. Both genes have the typical gene architecture and share conserved synteny with their mammalian homologues. Expression analysis, gene-targeted knockdown experiments and in vitro catalytic assays indicate that zebrafish Stx is the principal--if not unique--polysialyltransferase performing NCAM-PSA modifications in both developing and adult fish. The knockdown of Stx exclusively affects PSA synthesis, producing defects in axonal growth and guidance. Zebrafish Pst is in principle capable of synthesizing PSA, however, our data argue against a fundamental function of the enzyme during development. Our findings reveal an important divergence of Stx and Pst enzymes in vertebrates, which is also characterized by a differential gene loss and rapid evolution of Pst genes within the bony-fish class.

Comparison of hypoxia-inducible factor-1 alpha in hypoxia-sensitive and hypoxia-tolerant fish species

Levels of oxygen can vary dramatically in aquatic environments. Aquatic organisms, including fishes, have adapted accordingly to survive. As there are both phylogenetically closely related fish species with differing oxygen requirements and distantly related species with similar oxygen requirements, fishes are good candidates for examining oxygen-related functions in vertebrates. We set out to investigate if sequence variation in the hypoxia-inducible factor-1 alpha (HIF-1alpha) gene is associated with variations in oxygen requirements. Since the teleost HIF-1alpha sequences available in databases represent a very limited dataset both phylogenetically and with regard to oxygen requirements, we have sequenced the protein coding sequence for HIF-1alpha from an additional 9 fish species. Our results indicate that the deduced HIF-1alpha proteins of teleost fishes are somewhat shorter than those of tetrapods. Additionally, the results suggest that tetrapod sequences more closely resemble the ancestral form of the protein than do teleost sequences. No clear signatures which could be associated with the oxygen requirements of the species were found. This study suggests that if species-specific differences in HIF-1alpha function with regards to oxygen dependence have evolved, they do not occur in the protein coding sequence but at other levels of the HIF-1alpha pathway.

Halwaxiids and the early evolution of the lophotrochozoans

Halkieriids and wiwaxiids are cosmopolitan sclerite-bearing metazoans from the Lower and Middle Cambrian. Although they have similar scleritomes, their phylogenetic position is contested. A new scleritomous fossil from the Burgess Shale has the prominent anterior shell of the halkieriids but also bears wiwaxiid-like sclerites. This new fossil defines the monophyletic halwaxiids and indicates that they have a key place in early lophotrochozoan history.

Pax-Six-Eya-Dach network during amphioxus development: conservation in vitro but context specificity in vivo

The Drosophila retinal determination gene network occurs in animals generally as a Pax-Six-Eyes absent-Dachshund network (PSEDN). For amphioxus, we describe the complete network of nine PSEDN genes, four of which-AmphiSix1/2, AmphiSix4/5, AmphSix3/6, and AmphiEya-are characterized here for the first time. For amphioxus, in vitro interactions among the genes and proteins of the network resemble those of other animals, except for the absence of Dach-Eya binding. Amphioxus PSEDN genes are expressed in highly stage- and tissue-specific patterns (sometimes conspicuously correlated with the local intensity of cell proliferation) in the gastrular organizer, notochord, somites, anterior central nervous system, peripheral nervous system, pharyngeal endoderm, and the likely homolog of the vertebrate adenohypophysis. In this last tissue, the anterior region expresses all three amphioxus Six genes and is a zone of active cell proliferation, while the posterior region expresses only AmphiPax6 and is non-proliferative. In summary, the topologies of animal PSEDNs, although considerably more variable than originally proposed, are conserved enough to be recognizable among species and among developing tissues; this conservation may reflect indispensable involvement of PSEDNs during the critically important early phases of embryology (e.g. in the control of mitosis, apoptosis, and cell/tissue motility).

Hau-Pax3/7A is an early marker of leech mesoderm involved in segmental morphogenesis, nephridial development, and body cavity formation

Two genes of the Pax III subfamily, Hau-Pax3/7A and -Pax3/7B, were identified from the leech Helobdella, and the expression and function of Hau-Pax3/7A in development are described. Leech embryos undergo spiral cleavage, then produce a set of teloblastic stem cells that generate segmented mesoderm and ectoderm. Hau-Pax3/7A is present as a maternal transcript in both ectodermal and mesodermal progenitors, but this pool of early RNA disappears and is replaced by a pattern of zygotic transcription restricted to the blast cell progeny of the mesodermal M teloblasts. Each mesodermal blast cell clone goes through multiple phases of Hau-Pax3/7A expression, the last of which is associated with the organogenesis of the nephridia and other segment-specific structures. Morpholino-mediated knockdown of Hau-Pax3/7A expression causes the mesodermal blast cell clones to undergo irregular patterns of morphogenesis that disrupt the segmental organization of the germinal plate, and interferes with both the specification and morphological differentiation of the mesodermal nephridia. Knockdown of Hau-Pax3/7A in the mesoderm can also lead to abnormalities in the formation of the dorsal cavities, possibly through indirect effects of this germ layer on neighboring tissues. This is the first report of broad mesodermal Pax III expression outside of chordates, and raises the possibility that such expression may be a primitive trait inherited from the last common ancestor of the bilaterian superphyla.

Expression of the AmphiTcf gene in amphioxus: insights into the evolution of the TCF/LEF gene family during vertebrate evolution

T-cell factor (TCF) and lymphoid enhancer factors (LEF) genes encode proteins that are transcription factors mediating beta-catenin/Wnt signaling. Whereas mammals have four such genes, the Florida amphioxus (Branchiostoma floridae) apparently has only one such gene (AmphiTcf). From cleavage through early gastrula, cytoplasmic maternal transcripts of this gene are localized toward the animal pole. In gastrulae, AmphiTcf expression begins in the mesendoderm. In neurulae, there is expression in the pharynx, hindgut, anterior notochord, somites, and at the anterior end of the neural plate. In early larvae, expression is detectable in the floor of the diencephalon, notochord, tail bud, forming somites, pharynx, and ciliated pit (a presumed homolog of the vertebrate adenohypophysis). Phylogenetic analysis of TCF/LEF proteins placed AmphiTcf as the sister group of a clade comprising vertebrate Tcf1, Lef1, Tcf3, and Tcf4. Comparison of developmental expression for amphioxus AmphiTcf and vertebrate TCF/LEF genes indicates that this gene family has undergone extensive subfunctionalization and neofunctionalization during vertebrate evolution.

The identification of lymphocyte-like cells and lymphoid-related genes in amphioxus indicates the twilight for the emergence of adaptive immune system

To seek evidence of a primitive adaptive immune system (AIS) before vertebrate, we examined whether lymphocytes or lymphocyte-like cells and the related molecules participating in the lymphocyte function existed in amphioxus. Anatomical analysis by electron microscopy revealed the presence of lymphocyte-like cells in gills, and these cells underwent morphological changes in response to microbial pathogens that are reminiscent of those of mammalian lymphocytes executing immune response to microbial challenge. In addition, a systematic comparative analysis of our cDNA database of amphioxus identified a large number of genes whose vertebrate counterparts are involved in lymphocyte function. Among these genes, several genes were found to be expressed in the vicinity of the lymphocyte-like cells by in situ hybridization and up-regulated after exposure to microbial pathogens. Our findings in the amphioxus indicate the twilight for the emergence of AIS before the invertebrate-vertebrate transition during evolution.

Axial patterning in cephalochordates and the evolution of the organizer

The organizer of the vertebrate gastrula is an important signalling centre that induces and patterns dorsal axial structures. Although a topic of long-standing interest, the evolutionary origin of the organizer remains unclear. Here we show that the gastrula of the cephalochordate amphioxus expresses dorsal/ventral (D/V) patterning genes (for example, bone morphogenetic proteins (BMPs), Nodal and their antagonists) in patterns reminiscent of those of their vertebrate orthlogues, and that amphioxus embryos, like those of vertebrates, are ventralized by exogenous BMP protein. In addition, Wnt-antagonists (for example, Dkks and sFRP2-like) are expressed anteriorly, whereas Wnt genes themselves are expressed posteriorly, consistent with a role for Wnt signalling in anterior/posterior (A/P) patterning. These results suggest evolutionary conservation of the mechanisms for both D/V and A/P patterning of the early gastrula. In light of recent phylogenetic analyses placing cephalochordates basally in the chordate lineage, we propose that separate signalling centres for patterning the D/V and A/P axes may be an ancestral chordate character.

Continuous culture of two lancelets and production of the second filial generations in the laboratory

The lancelet is considered to be a promising laboratory model animal. To establish laboratory colonies of lancelet, we collected parental lancelets of Branchiostoma belcheri and B. japonicum (previously named as B. belcheri tsingtauense) with fully developed gonads from Xiamen Rare Marine Creature Conservation Areas (Fujian, China) on dates just before their spawning in the field in 2005. Those parental lancelets spawned spontaneously in the laboratory and produced thousands upon thousands of fertilized eggs. After carefully hatching and maintaining for almost 1 year, we successfully obtained about 500 first generation (F1) adults of B. belcheri and 3,300 of B. japonicum. Part of those F1 lancelets ripened and spontaneously spawned in 2006, and several thousands of second generation (F2) individuals of both species were produced. The young F2 lancelets are growing in good condition and some of B. japonicum initiated gonad development in December, 2006. Our experience emphasizes that cleanness of settlement substratum and sufficient food supply are important factors for long-term culture of lancelets in the laboratory.

Profile of acute immune response in Chinese amphioxus upon Staphylococcus aureus and Vibrio parahaemolyticus infection

Amphioxus is traditionally considered as the living invertebrate most closely related to vertebrate. However, no systematic study was performed about how the amphioxus defends against the microbial invasion. Here we reported a profile of gene transcription after Staphylococcus aureus (S.c) and Vibrio parahaemolyticus (V.p) challenged by suppression subtractive hybridization (SSH). When compared with mammals, amphioxus has the same acute immune defense genes (lectins, metalloproteinase, lysozymes and antimicrobial peptide, etc.) as well as a similar pattern and level of temporal gene expression. In contrast, amphioxus was demonstrated to have some novel acute immune response genes in response to the microbial challenge, such as apextrin and dermatopontin, which have a 3500-fold and 900-fold induction after the V.p infection, respectively, suggesting new functions in early immune system for these two genes. Our results reported for the first time a profile of primitive immune system defense against infection in protochordate.

Patterns of evolution of mitochondrial cytochrome c oxidase I and II DNA and implications for DNA barcoding

DNA barcoding has focused increasing attention on the use of specific regions of mitochondrial cytochrome c oxidase I and II genes (COI-COII) to diagnose and delimit species. However, our understanding of patterns of molecular evolution within these genes is limited. Here we examine patterns of nucleotide divergence in COI-COII within species and between species pairs of Lepidoptera and Diptera using a sliding window analysis. We found that: (1) locations of maximum divergence within COI-COII were highly variable among taxa surveyed in this study; (2) there was major overlap in divergence within versus between species, including within individual COI-COII profiles; (3) graphical DNA saturation analysis showed variation in percent nucleotide transitions throughout COI-COII and only limited association with levels of DNA divergence. Ultimately, no single optimally informative 600 bp location was found within the 2.3 kb of COI-COII, and the DNA barcoding region was no better than other regions downstream in COI. Consequently, we recommend that researchers should maximize sequence length to increase the probability of sampling regions of high phylogenetic informativeness, and to minimize stochastic variation in estimating total divergence.

Formation of the ascidian epidermal sensory neurons: insights into the origin of the chordate peripheral nervous system

The vertebrate peripheral nervous system (PNS) originates from neural crest and placodes. While its developmental origin is the object of intense studies, little is known concerning its evolutionary history. To address this question, we analyzed the formation of the larval tail PNS in the ascidian Ciona intestinalis. The tail PNS of Ciona is made of sensory neurons located within the epidermis midlines and extending processes in the overlying tunic median fin. We show that each midline corresponds to a single longitudinal row of epidermal cells and neurons sharing common progenitors. This simple organization is observed throughout the tail epidermis, which is made of only eight single-cell rows, each expressing a specific genetic program. We next demonstrate that the epidermal neurons are specified in two consecutive steps. During cleavage and gastrula stages, the dorsal and ventral midlines are independently induced by FGF9/16/20 and the BMP ligand ADMP, respectively. Subsequently, Delta/Notch-mediated lateral inhibition controls the number of neurons formed within these neurogenic regions. These results provide a comprehensive overview of PNS formation in ascidian and uncover surprising similarities between the fate maps and embryological mechanisms underlying formation of ascidian neurogenic epidermis midlines and the vertebrate median fin.

Dates from the molecular clock: how wrong can we be?

Large discrepancies have been found in dates of evolutionary events obtained using the molecular clock. Twofold differences have been reported between the dates estimated from molecular data and those from the fossil record; furthermore, different molecular methods can give dates that differ 20-fold. New software attempts to incorporate appropriate allowances for this uncertainty into the calculation of the accuracy of date estimates. Here, we propose that these innovations represent welcome progress towards obtaining reliable dates from the molecular clock, but warn that they are currently unproven, given that the causes and pattern of the discrepancies are the subject of ongoing research. This research implies that many previous studies, even some of those using recently developed methods, might have placed too much confidence in their date estimates, and their conclusions might need to be revised.

Analysis of cytoskeletal and motility proteins in the sea urchin genome assembly

The sea urchin embryo is a classical model system for studying the role of the cytoskeleton in such events as fertilization, mitosis, cleavage, cell migration and gastrulation. We have conducted an analysis of gene models derived from the Strongylocentrotus purpuratus genome assembly and have gathered strong evidence for the existence of multiple gene families encoding cytoskeletal proteins and their regulators in sea urchin. While many cytoskeletal genes have been cloned from sea urchin with sequences already existing in public databases, genome analysis reveals a significantly higher degree of diversity within certain gene families. Furthermore, genes are described corresponding to homologs of cytoskeletal proteins not previously documented in sea urchins. To illustrate the varying degree of sequence diversity that exists within cytoskeletal gene families, we conducted an analysis of genes encoding actins, specific actin-binding proteins, myosins, tubulins, kinesins, dyneins, specific microtubule-associated proteins, and intermediate filaments. We conducted ontological analysis of select genes to better understand the relatedness of urchin cytoskeletal genes to those of other deuterostomes. We analyzed developmental expression (EST) data to confirm the existence of select gene models and to understand their differential expression during various stages of early development.

Self/nonself discrimination at the basis of chordate evolution: limits on molecular conservation

All organisms rely on their capacity of self/nonself discrimination to rapidly detect approaching allogeneic cells as well as invading pathogenic microbes as foreign and to eliminate them. Failure to recognize nonself causes self-mating, germline parasitism and disease. Recent findings indicate that, in urochordates - the closest living relatives of vertebrates - different species use completely different molecules for allorecognition. Thanks to their phylogenetic position, these organisms might help us to understand the evolutionary origin of the vertebrate 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.

Ci-FoxA-a is the earliest zygotic determinant of the ascidian anterior ectoderm and directly activates Ci-sFRP1/5

This work focuses on the anteroposterior patterning of the ectoderm in the invertebrate chordate Ciona intestinalis. Previous work indicated that, by the eight-cell stage, the anterior and posterior animal blastomeres have acquired different properties, including a differential responsiveness to inducing signals from the underlying mesendoderm. Here, we investigated the molecular basis of this distinction. For this, we studied the regulation of the earliest marker specific for the anterior ectoderm, Ci-sFRP1/5, which is activated at the 64-cell stage. We first found that the activation of this marker in the anterior ectoderm does not involve communication with other lineages. We then identified, by phylogenetic footprinting and deletion analysis, a short conserved minimal enhancer driving the onset of expression of Ci-sFRP1/5. We showed that this enhancer was a direct target of the Ci-FoxA-a gene, a FoxA/HNF3 orthologue expressed in anterior ectodermal and mesendodermal lineages from the eight-cell stage. Gain- and loss-of-function experiments revealed that Ci-FoxA-a is necessary and sufficient within the ectoderm to impose an ectodermal anterior identity, and to repress the posterior programme. Thus, Ci-FoxA-a constitutes a major early zygotic anterior determinant for the ascidian ectoderm, acting autonomously in this territory, prior to the onset of vegetal inductions. Interestingly, while vertebrate FoxA2 are also involved in the regionalization of the ectoderm, they are thought to act during gastrulation to control, in the mesendoderm, the expression of organizer signals. We discuss the evolution of chordate ectodermal patterning in light of our findings.

Evolution of the Rho family of ras-like GTPases in eukaryotes

GTPases of the Rho family are molecular switches that play important roles in converting and amplifying external signals into cellular effects. Originally demonstrated to control the dynamics of the F-actin cytoskeleton, Rho GTPases have been implicated in many basic cellular processes that influence cell proliferation, differentiation, motility, adhesion, survival, or secretion. To elucidate the evolutionary history of the Rho family, we have analyzed over 20 species covering major eukaryotic clades from unicellular organisms to mammals, including platypus and opossum, and have reconstructed the ontogeny and the chronology of emergence of the different subfamilies. Our data establish that the 20 mammalian Rho members are structured into 8 subfamilies, among which Rac is the founder of the whole family. Rho, Cdc42, RhoUV, and RhoBTB subfamilies appeared before Coelomates and RhoJQ, Cdc42 isoforms, RhoDF, and Rnd emerged in chordates. In vertebrates, gene duplications and retrotranspositions increased the size of each chordate Rho subfamily, whereas RhoH, the last subfamily, arose probably by horizontal gene transfer. Rac1b, a Rac1 isoform generated by alternative splicing, emerged in amniotes, and RhoD, only in therians. Analysis of Rho mRNA expression patterns in mouse tissues shows that recent subfamilies have tissue-specific and low-level expression that supports their implication only in narrow time windows or in differentiated metabolic functions. These findings give a comprehensive view of the evolutionary canvas of the Rho family and provide guides for future structure and evolution studies of other components of Rho signaling pathways, in particular regulators of the RhoGEF family.

Meisetz and the birth of the KRAB motif

The largest family of transcription factors in mammals is of Cys(2)His(2) zinc finger-proteins, each with an NH(2)-terminal KRAB motif. Extensive expansions of this family have occurred in separate mammalian lineages, with approximately 400 such genes known in the human genome. Despite their widespread occurrence, the evolutionary provenance of the KRAB motif is unclear since previously it has not been found outside of the tetrapod vertebrates. Here, we show that homologues of the histone methyltransferase Meisetz are present within the sea urchin (Strongylocentrotus purpuratus) genome. Sea urchin and mammalian Meisetz sequences each contain an N-terminal KRAB motif, which thereby establishes an early origin of the KRAB motif prior to the divergence of echinoderm and chordate lineages. Finally, we present evidence that KRAB motifs derive from a novel family of KRI (KRAB Interior) motifs that were present in the last common ancestor of animals, plants and fungi.

Neofunctionalization in vertebrates: the example of retinoic acid receptors

Understanding the role of gene duplications in establishing vertebrate innovations is one of the main challenges of Evo-Devo (evolution of development) studies. Data on evolutionary changes in gene expression (i.e., evolution of transcription factor-cis-regulatory elements relationships) tell only part of the story; protein function, best studied by biochemical and functional assays, can also change. In this study, we have investigated how gene duplication has affected both the expression and the ligand-binding specificity of retinoic acid receptors (RARs), which play a major role in chordate embryonic development. Mammals have three paralogous RAR genes—RARα, β, and γ—which resulted from genome duplications at the origin of vertebrates. By using pharmacological ligands selective for specific paralogues, we have studied the ligand-binding capacities of RARs from diverse chordates species. We have found that RARβ-like binding selectivity is a synapomorphy of all chordate RARs, including a reconstructed synthetic RAR representing the receptor present in the ancestor of chordates. Moreover, comparison of expression patterns of the cephalochordate amphioxus and the vertebrates suggests that, of all the RARs, RARβ expression has remained most similar to that of the ancestral RAR. On the basis of these results together, we suggest that while RARβ kept the ancestral RAR role, RARα and RARγ diverged both in ligand-binding capacity and in expression patterns. We thus suggest that neofunctionalization occurred at both the expression and the functional levels to shape RAR roles during development in vertebrates.

In eukaryotic organisms, each gene is a stretch of DNA composed of control regions that bind transcription factors and coding regions that transcribe the mRNA that is later translated into proteins. At the molecular level, changes in control regions can affect the time and place at which a protein is synthesized, whereas changes in the coding region can alter the protein's function. Retinoic acid receptors (RARs) are chordate-specific transcription factors which, upon binding the natural morphogen retinoic acid, bind to and activate transcription from target genes. Here, the authors show how the ligand specificity of RARs has changed during vertebrate evolution in parallel with changes in expression. Through functional characterization of the RARs from several vertebrates, the chordate amphioxus, and the reconstructed ancestral RAR sequence, the authors show that of the three vertebrate RARs, RARβ has retained the ancestral characteristics in terms of both function and expression, while RARα and γ have evolved by acquiring new functions, both new binding specificity and new expression patterns. Thus both types of evolution have been important in the diversification of vertebrate RARs.

Evolution of the Neurotrophin Signaling System in Invertebrates

Nucleotide sequences encoding orthologs of neurotrophins and their receptors, p75(NTR) and Trk receptors, have been identified in the sea urchin Strongylocentrotus purpuratus, and the acorn worm, Saccoglossus kowalevskii, whereas the ascidian (sea squirt) species Ciona intestinalis and Ciona savignii appear to lack such orthologs. These results suggest that a functional neurotrophin system was already present at the beginning of deuterostome evolution, but was lost in ascidians. Remarkably, it appears that evolution of a p75(NTR) ortholog represented one of the earliest events in the expansion of the tumor necrosis factor receptor superfamily.

Formation and evolution of the chordate neurotrophin and Trk receptor genes

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

Genomic view of the evolution of the complement system

The recent accumulation of genomic information of many representative animals has made it possible to trace the evolution of the complement system based on the presence or absence of each complement gene in the analyzed genomes. Genome information from a few mammals, chicken, clawed frog, a few bony fish, sea squirt, fruit fly, nematoda and sea anemone indicate that bony fish and higher vertebrates share practically the same set of complement genes. This suggests that most of the gene duplications that played an essential role in establishing the mammalian complement system had occurred by the time of the teleost/mammalian divergence around 500 million years ago (MYA). Members of most complement gene families are also present in ascidians, although they do not show a one-to-one correspondence to their counterparts in higher vertebrates, indicating that the gene duplications of each gene family occurred independently in vertebrates and ascidians. The C3 and factor B genes, but probably not the other complement genes, are present in the genome of the cnidaria and some protostomes, indicating that the origin of the central part of the complement system was established more than 1,000 MYA.

A retinoic acid-Hox hierarchy controls both anterior/posterior patterning and neuronal specification in the developing central nervous system of the cephalochordate amphioxus

Retinoic acid (RA) mediates both anterior/posterior patterning and neuronal specification in the vertebrate central nervous system (CNS). However, the molecular mechanisms downstream of RA are not well understood. To investigate these mechanisms, we used the invertebrate chordate amphioxus, in which the CNS, although containing only about 20,000 neurons in adults, like the vertebrate CNS, has a forebrain, midbrain, hindbrain, and spinal cord and is regionalized by RA-signaling. Here we show, first, that domains of genes with expression normally limited to diencephalon and midbrain are generally not affected by altered RA-signaling, second, that contrary to previous reports, not only Hox1, 3, and 4, but also Hox2 and Hox6 are collinearly expressed in the amphioxus CNS, and third, that collinear expression of all these Hox genes is controlled by RA-signaling. Finally, we show that Hox1 is involved in mediating both the role of RA-signaling in regionalization of the hindbrain and in specification of hindbrain motor neurons. Thus, morpholino knock-down of the single amphioxus Hox1 mimics the effects of treatments with an RA-antagonist. This analysis establishes RA-dependent regulation of collinear Hox expression as a feature common to the chordate CNS and indicates that the RA-Hox hierarchy functions both in proper anterior/posterior patterning of the developing CNS and in specification of neuronal identity.

Building divergent body plans with similar genetic pathways

Deuterostome animals exhibit widely divergent body plans. Echinoderms have either radial or bilateral symmetry, hemichordates include bilateral enteropneust worms and colonial pterobranchs, and chordates possess a defined dorsal-ventral axis imposed on their anterior-posterior axis. Tunicates are chordates only as larvae, following metamorphosis the adults acquire a body plan unique for the deuterostomes. This paper examines larval and adult body plans in the deuterostomes and discusses two distinct ways of evolving divergent body plans. First, echinoderms and hemichordates have similar feeding larvae, but build a new adult body within or around their larvae. In hemichordates and many direct-developing echinoderms, the adult is built onto the larva, with the larval axes becoming the adult axes and the larval mouth becoming the adult mouth. In contrast, indirect-developing echinoderms undergo radical metamorphosis where adult axes are not the same as larval axes. A second way of evolving a divergent body plan is to become colonial, as seen in hemichordates and tunicates. Early embryonic development and gastrulation are similar in all deuterostomes, but, in chordates, the anterior-posterior axis is established at right angles to the animal-vegetal axis, in contrast to hemichordates and indirect-developing echinoderms. Hox gene sequences and anterior-posterior expression patterns illuminate deuterostome phylogenetic relationships and the evolution of unique adult body plans within monophyletic groups. Many genes that are considered vertebrate 'mesodermal' genes, such as nodal and brachyury T, are likely to ancestrally have been involved in the formation of the mouth and anus, and later were evolutionarily co-opted into mesoderm during vertebrate development.

A Mitogenomic Timescale for Birds Detects Variable Phylogenetic Rates of Molecular Evolution and Refutes the Standard Molecular Clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.

A Gbx homeobox gene in amphioxus: insights into ancestry of the ANTP class and evolution of the midbrain/hindbrain boundary

In the vertebrate central nervous system (CNS), mutual antagonism between posteriorly expressed Gbx2 and anteriorly expressed Otx2 positions the midbrain/hindbrain boundary (MHB), but does not induce MHB organizer genes such as En, Pax2/5/8 and Wnt1. In the CNS of the cephalochordate amphioxus, Otx is also expressed anteriorly, but En, Pax2/5/8 and Wnt1 are not expressed near the caudal limit of Otx, raising questions about the existence of an MHB organizer in amphioxus. To investigate the evolutionary origins of the MHB, we cloned the single amphioxus Gbx gene. Fluorescence in situ hybridization showed that, as in vertebrates, amphioxus Gbx and the Hox cluster are on the same chromosome. From analysis of linked genes, we argue that during evolution a single ancestral Gbx gene duplicated fourfold in vertebrates, with subsequent loss of two duplicates. Amphioxus Gbx is expressed in all germ layers in the posterior 75% of the embryo, and in the CNS, the Gbx and Otx domains abut at the boundary between the cerebral vesicle (forebrain/midbrain) and the hindbrain. Thus, the genetic machinery to position the MHB was present in the protochordate ancestors of the vertebrates, but is insufficient for induction of organizer genes. Comparison with hemichordates suggests that anterior Otx and posterior Gbx domains were probably overlapping in the ancestral deuterostome and came to abut at the MHB early in the chordate lineage before MHB organizer properties evolved.

Hox genes are not always Colinear

The deuterostomes are the clade of animals for which we have the most detailed understanding of Hox cluster organisation. With the Hox cluster of amphioxus (Branchiostoma floridae) we have the best prototypical, least derived Hox cluster for the group, whilst the urochordates present us with some of the most highly derived and disintegrated clusters. Combined with the detailed mechanistic understanding of vertebrate Hox regulation, the deuterostomes provide much of the most useful data for understanding Hox cluster evolution. Considering both the prototypical and derived deuterostome Hox clusters leads us to hypothesize that Temporal Colinearity is the main constraining force on Hox cluster organisation, but until we have a much deeper understanding of the mechanistic basis for this phenomenon, and know how widespread across the Bilateria the mechanism(s) is/are, then we cannot know how the Hox cluster of the last common bilaterian operated and what have been the major evolutionary forces operating upon the Hox gene cluster.

Amphioxus and tunicates as evolutionary model systems

One important question in evolutionary biology concerns the origin of vertebrates from invertebrates. The current consensus is that the proximate ancestor of vertebrates was an invertebrate chordate. Today, the invertebrate chordates comprise cephalochordates (amphioxus) and tunicates (each a subphylum in the phylum Chordata, which also includes the vertebrate subphylum). It was widely accepted that, within the chordates, tunicates represent the sister group of a clade of cephalochordates plus vertebrates. However, recent studies suggest that the evolutionary positions of tunicates and cephalochordates should be reversed, the implications of which are considered here. We also review the two major groups of invertebrate chordates and compare relative advantages (and disadvantages) of each as model systems for elucidating the origin of the vertebrates.

A SINE in the genome of the cephalochordate amphioxus is an Alu element

Transposable elements of about 300 bp, termed "short interspersed nucleotide elements or SINEs are common in eukaryotes. However, Alu elements, SINEs containing restriction sites for the AluI enzyme, have been known only from primates. Here I report the first SINE found in the genome of the cephalochordate, amphioxus. It is an Alu element of 375 bp that does not share substantial identity with any genomic sequences in vertebrates. It was identified because it was located in the FoxD regulatory region in a cosmid derived from one individual, but absent from the two FoxD alleles of BACs from a second individual. However, searches of sequences of BACs and genomic traces from this second individual gave an estimate of 50-100 copies in the amphioxus genome. The finding of an Alu element in amphioxus raises the question of whether Alu elements in amphioxus and primates arose by convergent evolution or by inheritance from a common ancestor. Genome-wide analyses of transposable elements in amphioxus and other chordates such as tunicates, agnathans and cartilaginous fishes could well provide the answer.

Tunicates and not cephalochordates are the closest living relatives of vertebrates

Tunicates or urochordates (appendicularians, salps and sea squirts), cephalochordates (lancelets) and vertebrates (including lamprey and hagfish) constitute the three extant groups of chordate animals. Traditionally, cephalochordates are considered as the closest living relatives of vertebrates, with tunicates representing the earliest chordate lineage. This view is mainly justified by overall morphological similarities and an apparently increased complexity in cephalochordates and vertebrates relative to tunicates. Despite their critical importance for understanding the origins of vertebrates, phylogenetic studies of chordate relationships have provided equivocal results. Taking advantage of the genome sequencing of the appendicularian Oikopleura dioica, we assembled a phylogenomic data set of 146 nuclear genes (33,800 unambiguously aligned amino acids) from 14 deuterostomes and 24 other slowly evolving species as an outgroup. Here we show that phylogenetic analyses of this data set provide compelling evidence that tunicates, and not cephalochordates, represent the closest living relatives of vertebrates. Chordate monophyly remains uncertain because cephalochordates, albeit with a non-significant statistical support, surprisingly grouped with echinoderms, a hypothesis that needs to be tested with additional data. This new phylogenetic scheme prompts a reappraisal of both morphological and palaeontological data and has important implications for the interpretation of developmental and genomic studies in which tunicates and cephalochordates are used as model animals.

When Brachyury meets Smad1: the evolution of bilateral symmetry during gastrulation

Understanding the events that led to the emergence of the bilaterians is a daunting task, impaired by the huge evolutionary gap separating us from the pre-Cambrian. During gastrulation, the expression of the transcription factor Brachyury is remarkably well conserved around the blastopore of bilaterians and cnidarians. Only the bilaterian Brachyury proteins, however, share a distinctive N-terminal sequence not found in outgroups such as cnidarians, sponges or placozoans. We now know that, in vertebrates, this N-terminal domain confers specific transcriptional activity, by recruiting Smad1, the first identified co-factor for Brachyury. Smad1 is an effector of the BMP pathway, and has been isolated in bilaterians and cnidarians. Here, I propose that the protein-protein interaction between Brachyury and Smad1 represents an evolutionary novelty of the Urbilateria. The gain of the N-terminal domain might have been selected to spatially modulate the activity of Brachyury, thereby facilitating the establishment of bilateral symmetry during gastrulation movements.

Evolution and development of the chordates: collagen and pharyngeal cartilage

Chordates evolved a unique body plan within deuterostomes and are considered to share five morphological characters, a muscular postanal tail, a notochord, a dorsal neural tube, an endostyle, and pharyngeal gill slits. The phylum Chordata typically includes three subphyla, Cephalochordata, Vertebrata, and Tunicata, the last showing a chordate body plan only as a larva. Hemichordates, in contrast, have pharyngeal gill slits, an endostyle, and a postanal tail but appear to lack a notochord and dorsal neural tube. Because hemichordates are the sister group of echinoderms, the morphological features shared with the chordates must have been present in the deuterostome ancestor. No extant echinoderms share any of the chordate features, so presumably they have lost these structures evolutionarily. We review the development of chordate characters in hemichordates and present new data characterizing the pharyngeal gill slits and their cartilaginous gill bars. We show that hemichordate gill bars contain collagen and proteoglycans but are acellular. Hemichordates and cephalochordates, or lancelets, show strong similarities in their gill bars, suggesting that an acellular cartilage may have preceded cellular cartilage in deuterostomes. Our evidence suggests that the deuterostome ancestor was a benthic worm with gill slits and acellular gill cartilages.