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

Morphology and phylogenetic significance of the pelvic articular region in elasmobranchs (Chondrichthyes)

The morphology of paired fins is commonly overlooked in morphological studies, particularly the pelvic girdle and fins. Consequently, previous phylogenetic studies incorporating morphological data used few skeletal characters from this complex. In this paper, the phylogenetic significance of pelvic articular characters for elasmobranchs is discussed in light of the morphological variation observed in 130 species, the most comprehensive study exploring the morphology of the pelvic girdle done so far. The 10 morphological characters proposed herein for the pelvic articulation were incorporated into a molecular matrix of NADH2 sequences and submitted to an analysis of maximum parsimony employing extended implied weighting. The most stable tree was selected based on the distortion coefficients, SPR distances (subtree pruning and regrafting) and fit values. Some of the striking synapomorphies recovered within elasmobranchs include the presence of an articular surface for the first enlarged pelvic radial supporting Elasmobranchii and the pelvic articular region for the basipterygium extending from the posterolatral margin of the pelvic girdle over its lateral surface in Echinorhinus + Hexanchiformes. Additionally, the proposed characters and their distributions are discussed considering the relationships recovered and also compared with previous morphological and molecular phylogenetic hypotheses.

Worms and gills, plates and spines: the evolutionary origins and incredible disparity of deuterostomes revealed by fossils, genes, and development

Deuterostomes are the major division of animal life which includes sea stars, acorn worms, and humans, among a wide variety of ecologically and morphologically disparate taxa. However, their early evolution is poorly understood, due in part to their disparity, which makes identifying commonalities difficult, as well as their relatively poor early fossil record. Here, we review the available morphological, palaeontological, developmental, and molecular data to establish a framework for exploring the origins of this important and enigmatic group. Recent fossil discoveries strongly support a vermiform ancestor to the group Hemichordata, and a fusiform active swimmer as ancestor to Chordata. The diverse and anatomically bewildering variety of forms among the early echinoderms show evidence of both bilateral and radial symmetry. We consider four characteristics most critical for understanding the form and function of the last common ancestor to Deuterostomia: Hox gene expression patterns, larval morphology, the capacity for biomineralization, and the morphology of the pharyngeal region. We posit a deuterostome last common ancestor with a similar antero-posterior gene regulatory system to that found in modern acorn worms and cephalochordates, a simple planktonic larval form, which was later elaborated in the ambulacrarian lineage, the ability to secrete calcium minerals in a limited fashion, and a pharyngeal respiratory region composed of simple pores. This animal was likely to be motile in adult form, as opposed to the sessile origins that have been historically suggested. Recent debates regarding deuterostome monophyly as well as the wide array of deuterostome-affiliated problematica further suggest the possibility that those features were not only present in the last common ancestor of Deuterostomia, but potentially in the ur-bilaterian. The morphology and development of the early deuterostomes, therefore, underpin some of the most significant questions in the study of metazoan evolution.

Complete Mitochondrial DNA Genome of Nine Species of Sharks and Rays and Their Phylogenetic Placement among Modern Elasmobranchs

Chondrichthyes occupy a key position in the phylogeny of vertebrates. The complete sequence of the mitochondrial genome (mitogenome) of four species of sharks and five species of rays was obtained by whole genome sequencing (DNA-seq) in the Illumina HiSeq2500 platform. The arrangement and features of the genes in the assembled mitogenomes were identical to those found in vertebrates. Both Maximum Likelihood (ML) and Bayesian Inference (BI) analyses were used to reconstruct the phylogenetic relationships among 172 species (including 163 mitogenomes retrieved from GenBank) based on the concatenated dataset of 13 individual protein coding genes. Both ML and BI analyses did not support the “Hypnosqualea” hypothesis and confirmed the monophyly of sharks and rays. The broad notion in shark phylogeny, namely the division of sharks into Galeomorphii and Squalomorphii and the monophyly of the eight shark orders, was also supported. The phylogenetic placement of all nine species sequenced in this study produced high statistical support values. The present study expands our knowledge on the systematics, genetic differentiation, and conservation genetics of the species studied, and contributes to our understanding of the evolutionary history of Chondrichthyes.

The phylum Vertebrata: a case for zoological recognition

The group Vertebrata is currently placed as a subphylum in the phylum Chordata, together with two other subphyla, Cephalochordata (lancelets) and Urochordata (ascidians). The past three decades, have seen extraordinary advances in zoological taxonomy and the time is now ripe for reassessing whether the subphylum position is truly appropriate for vertebrates, particularly in light of recent advances in molecular phylogeny, comparative genomics, and evolutionary developmental biology. Four lines of current research are discussed here. First, molecular phylogeny has demonstrated that Deuterostomia comprises Ambulacraria (Echinodermata and Hemichordata) and Chordata (Cephalochordata, Urochordata, and Vertebrata), each clade being recognized as a mutually comparable phylum. Second, comparative genomic studies show that vertebrates alone have experienced two rounds of whole-genome duplication, which makes the composition of their gene family unique. Third, comparative gene-expression profiling of vertebrate embryos favors an hourglass pattern of development, the most conserved stage of which is recognized as a phylotypic period characterized by the establishment of a body plan definitively associated with a phylum. This mid-embryonic conservation is supported robustly in vertebrates, but only weakly in chordates. Fourth, certain complex patterns of body plan formation (especially of the head, pharynx, and somites) are recognized throughout the vertebrates, but not in any other animal groups. For these reasons, we suggest that it is more appropriate to recognize vertebrates as an independent phylum, not as a subphylum of the phylum Chordata.

A simulated heat wave shortens the telomere length and lifespan of a desert lizard

Understanding how organisms respond to warming contributes important information to the conservation of biodiversity that is threatened by climate warming. Here, we conducted experiments on a desert agama (Phrynocephalus przewalskii) to test the hypothesis that climate warming (an increase in both mean temperature and heat waves) would induce oxidative stress, shortening telomere length, and thereby decreasing survival. Our results demonstrated that one week of exposure to a simulated heat wave significantly shortened telomere length, and decreased the overwinter survival of lizards, but mean temperature increase did not affect the survival of lizards. However, the antioxidant capacity (anti-oxidative enzyme) was not affected by the warming treatments. Therefore, heat waves might have negative impacts on the desert agama, with shortened telomeres likely causing the lifespan of lizards to decrease under climate warming.

An evolutionary perspective on the role of mesencephalic astrocyte-derived neurotrophic factor (MANF): At the crossroads of poriferan innate immune and apoptotic pathways

The mesencephalic astrocyte-derived neurotrophic factor (MANF) belongs to a recently discovered family of neurotrophic factors. MANF can be secreted but is generally resident within the endoplasmic reticulum (ER) in neuronal and non-neuronal cells, where it is involved in the ER stress response with pro-survival effects. Here we report the discovery of the MANF homolog SDMANF in the sponge Suberites domuncula. The basal positioning of sponges (phylum Porifera) in the animal tree of life offers a unique vantage point on the early evolution of the metazoan-specific genetic toolkit and molecular pathways. Since sponges lack a conventional nervous system, SDMANF presents an enticing opportunity to investigate the evolutionary ancient role of these neurotrophic factors. SDMANF shares considerable sequence similarity with its metazoan homologs. It also comprises a putative protein binding domain with sequence similarities to the Bcl-2 family of apoptotic regulators. In Suberites, SDMANF is expressed in the vicinity of bacteriocytes, where it co-localizes with the toll-like receptor SDTLR. In transfected human cells, SDMANF was detected in both the organelle protein fraction and the cell culture medium. The intracellular SDMANF protein level was up-regulated in response to both a Golgi/ER transport inhibitor and bacterial lipopolysaccharides (LPS). Upon LPS challenge, transfected cells revealed a decreased caspase-3 activity and increased cell viability with no inducible Bax expression compared to the wild type. These results suggest a deep evolutionary original cytoprotective role of MANF, at the crossroads of innate immune and apoptotic pathways, of which a neurotrophic function might have arisen later in metazoan evolution.

Complete mitochondrial DNA genome of bonnethead shark, Sphyrna tiburo, and phylogenetic relationships among main superorders of modern elasmobranchs

Elasmobranchs are one of the most diverse groups in the marine realm represented by 18 orders, 55 families and about 1200 species reported, but also one of the most vulnerable to exploitation and to climate change. Phylogenetic relationships among main orders have been controversial since the emergence of the Hypnosqualean hypothesis by Shirai (1992) that considered batoids as a sister group of sharks. The use of the complete mitochondrial DNA (mtDNA) may shed light to further validate this hypothesis by increasing the number of informative characters. We report the mtDNA genome of the bonnethead shark Sphyrna tiburo, and compare it with mitogenomes of other 48 species to assess phylogenetic relationships. The mtDNA genome of S. tiburo, is quite similar in size to that of congeneric species but also similar to the reported mtDNA genome of other Carcharhinidae species. Like most vertebrate mitochondrial genomes, it contained 13 protein coding genes, two rRNA genes and 22 tRNA genes and the control region of 1086 bp (D-loop). The Bayesian analysis of the 49 mitogenomes supported the view that sharks and batoids are separate groups.

Cloning, expression and localization of DacaCSP2 and DacaCSP3 during different reproductive stages in Daphnia carinata

Daphnia carinata are unique freshwater crustaceans that undergo both sexual and asexual reproduction, depending on environmental factors. While the molecular mechanism behind the reproductive transformation has been unknown, chemosensory proteins (CSPs) may be involved. We have cloned the cDNA sequences of two CSP genes from D. carinata using primers based on homologous sequences and rapid amplification of cDNA ends (RACE). The full-length cDNA of DacaCSP2 (GenBank accession no: KM624608) was 632 bp, with an ORF (open reading frame) of 330 bp encoding a 12.02 kDa protein; and the full-length cDNA of DacaCSP3 (GenBank accession no: KM624609) was 935 bp, with an ORF of 342 bp encoding a 12.78kDa protein. Both CSPs encoded an N-terminal signal peptide, four conserved cysteines, an OS-D superfamily domain, a 2Fe-2S ferredoxin domain, an anaphylatoxin domain and an EGF-like domain. DaCaCSP2 and DaCaCSP3 proteins were most closely related to CSPs from Daphnia pulex and were more distantly related to CSPs from other insects. Using quantitative PCR, we found expression levels of DaCaCSP2 and DaCaCSP3 mRNA were highest in sexual females, followed by parthenogenetic females, and lowest in males. The expression levels of DaCaCSP2 and DaCaCSP3 mRNA also increased at lower temperatures, which suggested they could respond to environmental cues. Whole mount in situ hybridization (ISH) showed that DaCaCSP2 and DaCaCSP3 were expressed mainly in the ovaries, summer eggs, thoracic limbs, rectum and second antennae in sexual females; while they were expressed mainly in the ovaries, thoracic limbs, rectum and second antennae in parthenogenetic females. Together, these results suggest that DacaCSP2 and DacaCSP3 may respond to environmental cues and control the reproductive switch from sexual to asexual reproduction in D. carinata.

Molecular signatures that are distinctive characteristics of the vertebrates and chordates and supporting a grouping of vertebrates with the tunicates

Members of the phylum Chordata and the subphylum Vertebrata are presently distinguished solely on the basis of morphological characteristics. The relationship of the vertebrates to the two non-vertebrate chordate subphyla is also a subject of debate. Analyses of protein sequences have identified multiple conserved signature indels (CSIs) that are specific for Chordata or for Vertebrata. Five CSIs in 4 important proteins are specific for the Vertebrata, whereas two other CSIs are uniquely found in all sequenced chordate species including Ciona intestinalis and Oikapleura dioica (Tunicates) as well as Branchiostoma floridae (Cephalochordates). The shared presence of these molecular signatures by all vertebrates/chordate species, but in no other animal taxa, strongly indicates that the genetic changes represented by the identified CSIs diagnose monophyletic groups. Two other discovered CSIs are uniquely shared by different vertebrate species and by either one (Ciona intestinalis) or both tunicate (Ciona and Oikapleura) species, but they are not found in Branchiostoma or other animal species. Specific presence of these CSIs in different vertebrates and either one or both tunicate species provides strong independent evidence that the vertebrate species are more closely related to the urochordates (tunicates) than to the cephalochordates.

'Fish' (Actinopterygii and Elasmobranchii) diversification patterns through deep time

Actinopterygii (ray-finned fishes) and Elasmobranchii (sharks, skates and rays) represent more than half of today's vertebrate taxic diversity (approximately 33000 species) and form the largest component of vertebrate diversity in extant aquatic ecosystems. Yet, patterns of 'fish' evolutionary history remain insufficiently understood and previous studies generally treated each group independently mainly because of their contrasting fossil record composition and corresponding sampling strategies. Because direct reading of palaeodiversity curves is affected by several biases affecting the fossil record, analytical approaches are needed to correct for these biases. In this review, we propose a comprehensive analysis based on comparison of large data sets related to competing phylogenies (including all Recent and fossil taxa) and the fossil record for both groups during the Mesozoic-Cainozoic interval. This approach provides information on the 'fish' fossil record quality and on the corrected 'fish' deep-time phylogenetic palaeodiversity signals, with special emphasis on diversification events. Because taxonomic information is preserved after analytical treatment, identified palaeodiversity events are considered both quantitatively and qualitatively and put within corresponding palaeoenvironmental and biological settings. Results indicate a better fossil record quality for elasmobranchs due to their microfossil-like fossil distribution and their very low diversity in freshwater systems, whereas freshwater actinopterygians are diverse in this realm with lower preservation potential. Several important diversification events are identified at familial and generic levels for elasmobranchs, and marine and freshwater actinopterygians, namely in the Early-Middle Jurassic (elasmobranchs), Late Jurassic (actinopterygians), Early Cretaceous (elasmobranchs, freshwater actinopterygians), Cenomanian (all groups) and the Paleocene-Eocene interval (all groups), the latter two representing the two most exceptional radiations among vertebrates. For each of these events along with the Cretaceous-Paleogene extinction, we provide an in-depth review of the taxa involved and factors that may have influenced the diversity patterns observed. Among these, palaeotemperatures, sea-levels, ocean circulation and productivity as well as continent fragmentation and environment heterogeneity (reef environments) are parameters that largely impacted on 'fish' evolutionary history, along with other biotic constraints.

A higher level classification of all living organisms

We present a consensus classification of life to embrace the more than 1.6 million species already provided by more than 3,000 taxonomists' expert opinions in a unified and coherent, hierarchically ranked system known as the Catalogue of Life (CoL). The intent of this collaborative effort is to provide a hierarchical classification serving not only the needs of the CoL's database providers but also the diverse public-domain user community, most of whom are familiar with the Linnaean conceptual system of ordering taxon relationships. This classification is neither phylogenetic nor evolutionary but instead represents a consensus view that accommodates taxonomic choices and practical compromises among diverse expert opinions, public usages, and conflicting evidence about the boundaries between taxa and the ranks of major taxa, including kingdoms. Certain key issues, some not fully resolved, are addressed in particular. Beyond its immediate use as a management tool for the CoL and ITIS (Integrated Taxonomic Information System), it is immediately valuable as a reference for taxonomic and biodiversity research, as a tool for societal communication, and as a classificatory "backbone" for biodiversity databases, museum collections, libraries, and textbooks. Such a modern comprehensive hierarchy has not previously existed at this level of specificity.

Cloning, expression and localization of the Daphnia carinata transformer gene DcarTra during different reproductive stages

In this study, the full-length cDNA of the Transformer (Tra) gene from the common freshwater species Daphnia carinata (DcarTra; GenBank accession no. KJ735445) was cloned using primers based on homologous sequences and rapid amplification of cDNA ends (RACE). The relative expression and localization of DcarTra and the cellular abundance of the DcarTra protein during different sexual phases were subsequently investigated. The full-length DcarTra cDNA was 1620 bp with an ORF of 1143 bp encoding a 380 amino acid polypeptide. Phylogenetic analysis identified closely related genes in Daphnia magna and Daphnia pulex, and more distantly related genes in other insects. Quantitative PCR showed that DcarTra expression was highest in males, followed by sexual females, and lowest in parthenogenetic females. Whole-mount in situ hybridization showed that DcarTra was mainly expressed in the thoracic limbs, ovaries and rectum in parthenogenetic females, and in the joints of second antennae, ovaries, rectum and ventral processes in sexual females. Western blotting showed two differently phosphorylated forms of the Tra protein. When Tra is phosphorylated, DcarTra protein levels were much higher in males than in two females. Otherwise, when Tra is dephosphorylated, the highest Tra protein levels were in sexual females, which revealed that D. carinata can control the sexual transition via these two forms. Together these results suggest that DcarTra plays significant roles in the reproductive transformation of D. carinata and dephosphorylation of DcarTra may be the trigger for females to transform into males.

On the peculiar morphology and development of the hypoglossal, glossopharyngeal and vagus nerves and hypobranchial muscles in the hagfish

INTRODUCTION

The vertebrate body is characterized by its dual segmental organization: pharyngeal arches in the head and somites in the trunk. Muscular and nervous system morphologies are also organized following these metameric patterns, with distinct differences between head and trunk; branchiomeric nerves innervating pharyngeal arches are superficial to spinal nerves innervating somite derivatives. Hypobranchial muscles originate from rostral somites and occupy the "neck" at the head-trunk interface. Hypobranchial muscles, unlike ventral trunk muscles in the lateral body wall, develop from myocytes that migrate ventrally to occupy a space that is ventrolateral to the pharynx and unassociated with coelomic cavities. Occipitospinal nerves innervating these muscles also extend ventrally, thereby crossing the vagus nerve laterally.

RESULTS

In hagfishes, the basic morphological pattern of vertebrates is obliterated by the extreme caudal shift of the posterior part of the pharynx. The vagus nerve is found unusually medially, and occipitospinal nerves remain unfasciculated, appearing as metameric spinal nerves as in the posterior trunk region. Moreover, the hagfish exhibits an undifferentiated body plan, with the hypobranchial muscles not well dissociated from the abaxial muscles in the trunk. Comparative embryological observation showed that this hagfish-specific morphology is established by secondary modification of the common vertebrate embryonic pattern, and the hypobranchial muscle homologue can be found in the rostral part of the oblique muscle with pars decussata.

CONCLUSION

The morphological pattern of the hagfish represents an extreme case of heterotopy that led to the formation of the typical hypoglossal nerve, and can be regarded as an autapomorphic trait of the hagfish lineage.

Phylogenomic analysis of echinoderm class relationships supports Asterozoa

While some aspects of the phylogeny of the five living echinoderm classes are clear, the position of the ophiuroids (brittlestars) relative to asteroids (starfish), echinoids (sea urchins) and holothurians (sea cucumbers) is controversial. Ophiuroids have a pluteus-type larva in common with echinoids giving some support to an ophiuroid/echinoid/holothurian clade named Cryptosyringida. Most molecular phylogenetic studies, however, support an ophiuroid/asteroid clade (Asterozoa) implying either convergent evolution of the pluteus or reversals to an auricularia-type larva in asteroids and holothurians. A recent study of 10 genes from four of the five echinoderm classes used 'phylogenetic signal dissection' to separate alignment positions into subsets of (i) suboptimal, heterogeneously evolving sites (invariant plus rapidly changing) and (ii) the remaining optimal, homogeneously evolving sites. Along with most previous molecular phylogenetic studies, their set of heterogeneous sites, expected to be more prone to systematic error, support Asterozoa. The homogeneous sites, in contrast, support an ophiuroid/echinoid grouping, consistent with the cryptosyringid clade, leading them to posit homology of the ophiopluteus and echinopluteus. Our new dataset comprises 219 genes from all echinoderm classes; analyses using probabilistic Bayesian phylogenetic methods strongly support Asterozoa. The most reliable, slowly evolving quartile of genes also gives highest support for Asterozoa; this support diminishes in second and third quartiles and the fastest changing quartile places the ophiuroids close to the root. Using phylogenetic signal dissection, we find heterogenous sites support an unlikely grouping of Ophiuroidea + Holothuria while homogeneous sites again strongly support Asterozoa. Our large and taxonomically complete dataset finds no support for the cryptosyringid hypothesis; in showing strong support for the Asterozoa, our preferred topology leaves the question of homology of pluteus larvae open.

The biology of the germ line in echinoderms

The formation of the germ line in an embryo marks a fresh round of reproductive potential. The developmental stage and location within the embryo where the primordial germ cells (PGCs) form, however, differs markedly among species. In many animals, the germ line is formed by an inherited mechanism, in which molecules made and selectively partitioned within the oocyte drive the early development of cells that acquire this material to a germ-line fate. In contrast, the germ line of other animals is fated by an inductive mechanism that involves signaling between cells that directs this specialized fate. In this review, we explore the mechanisms of germ-line determination in echinoderms, an early-branching sister group to the chordates. One member of the phylum, sea urchins, appears to use an inherited mechanism of germ-line formation, whereas their relatives, the sea stars, appear to use an inductive mechanism. We first integrate the experimental results currently available for germ-line determination in the sea urchin, for which considerable new information is available, and then broaden the investigation to the lesser-known mechanisms in sea stars and other echinoderms. Even with this limited insight, it appears that sea stars, and perhaps the majority of the echinoderm taxon, rely on inductive mechanisms for germ-line fate determination. This enables a strongly contrasted picture for germ-line determination in this phylum, but one for which transitions between different modes of germ-line determination might now be experimentally addressed.

The development of lateral line placodes: taking a broader view

The lateral line system of anamniote vertebrates enables the detection of local water movement and weak bioelectric fields. Ancestrally, it comprises neuromasts - small sense organs containing mechanosensory hair cells - distributed in characteristic lines over the head and trunk, flanked on the head by fields of electroreceptive ampullary organs, innervated by afferent neurons projecting respectively to the medial and dorsal octavolateral nuclei in the hindbrain. Given the independent loss of the electrosensory system in multiple lineages, the development and evolution of the mechanosensory and electrosensory components of the lateral line must be dissociable. Nevertheless, the entire system arises from a series of cranial lateral line placodes, which exhibit two modes of sensory organ formation: elongation to form sensory ridges that fragment (with neuromasts differentiating in the center of the ridge, and ampullary organs on the flanks), or migration as collectives of cells, depositing sense organs in their wake. Intensive study of the migrating posterior lateral line placode in zebrafish has yielded a wealth of information concerning the molecular control of migration and neuromast formation in this migrating placode, in this cypriniform teleost species. However, our mechanistic understanding of neuromast and ampullary organ formation by elongating lateral line placodes, and even of other zebrafish lateral line placodes, is sparse or non-existent. Here, we attempt to highlight the diversity of lateral line development and the limits of the current research focus on the zebrafish posterior lateral line placode. We hope this will stimulate a broader approach to this fascinating sensory system.

The GC content of LSU rRNA evolves across topological and functional regions of the ribosome in all three domains of life

Large-subunit rRNA is the ribozyme that catalyzes protein synthesis by translation, and many of its features vary along a deep-to-superficial gradient. By measuring the G+C proportions in this rRNA in all three domains of life (60 bacteria, 379 eukaryote, and 23 archaean sequences), we tested whether the proportion of GC nucleotides varies along this in-out gradient. The rRNA regions used were several zones identified by Bokov and Steinberg (2009) as being arranged from deep to superficial within the LSU. To the Bokov-Steinberg zones, we added the most superficial zone of all, the divergent domains (expansion segments), which are greatly enlarged in eukaryotes. Regression lines constructed from the hundreds of species of organisms revealed the expected in-out gradient, showing that species with high %GC (or high %AT) in their rRNA distribute more of these abundant nucleotides into the peripheral zones. This could be explained by the evolutionary rates of replacement of all nucleotides (A, C, G, T), because these latter rates are fastest at the periphery and slowest near the conserved core. As an overall explanation, we propose that when extrinsic factors (whole-genome nucleotide composition, or environmental temperature) demand the percentage of GC in the rRNA of a species be high or low, then the deep-lying zones are buffered against GC variation because they are the slowest to evolve. The deep, conserved zones are also the most involved in translation, hinting that stabilizing selection there prevents a high GC variability that would diminish LSU rRNA's core functions. We found only a few domain-specific trends in rRNA-GC distribution, which relate to many Archaea living at high temperatures or to the highly complex genes and adaptations of Eukaryota. Use of rRNA sequences in molecular phylogenetic studies, for reconstructing the relationships of organisms across the tree of life, requires accurate models of how rRNA evolves. The demonstration that GC distributes in regular patterns across rRNA regions can improve these tree-reconstruction models in the future and should yield phylogenies of greater accuracy.

The evolution and development of vertebrate lateral line electroreceptors

Electroreception is an ancient vertebrate sense with a fascinating evolutionary history involving multiple losses as well as independent evolution at least twice within teleosts. We review the phylogenetic distribution of electroreception and the morphology and innervation of electroreceptors in different vertebrate groups. We summarise recent work from our laboratory that has confirmed the homology of ampullary electroreceptors in non-teleost jawed vertebrates by showing, in conjunction with previously published work, that these are derived embryonically from lateral line placodes. Finally, we review hypotheses to explain the distribution of electroreception within teleosts, including the hypothesis that teleost ampullary and tuberous electroreceptors evolved via the modification of mechanosensory hair cells in lateral line neuromasts. We conclude that further experimental work on teleost electroreceptor development is needed to test such hypotheses.

The evolutionary and genetic origins of consciousness in the Cambrian Period over 500 million years ago

Vertebrates evolved in the Cambrian Period before 520 million years ago, but we do not know when or how consciousness arose in the history of the vertebrate brain. Here we propose multiple levels of isomorphic or somatotopic neural representations as an objective marker for sensory consciousness. All extant vertebrates have these, so we deduce that consciousness extends back to the group's origin. The first conscious sense may have been vision. Then vision, coupled with additional sensory systems derived from ectodermal placodes and neural crest, transformed primitive reflexive systems into image forming brains that map and perceive the external world and the body's interior. We posit that the minimum requirement for sensory consciousness and qualia is a brain including a forebrain (but not necessarily a developed cerebral cortex/pallium), midbrain, and hindbrain. This brain must also have (1) hierarchical systems of intercommunicating, isomorphically organized, processing nuclei that extensively integrate the different senses into representations that emerge in upper levels of the neural hierarchy; and (2) a widespread reticular formation that integrates the sensory inputs and contributes to attention, awareness, and neural synchronization. We propose a two-step evolutionary history, in which the optic tectum was the original center of multi-sensory conscious perception (as in fish and amphibians: step 1), followed by a gradual shift of this center to the dorsal pallium or its cerebral cortex (in mammals, reptiles, birds: step 2). We address objections to the hypothesis and call for more studies of fish and amphibians. In our view, the lamprey has all the neural requisites and is likely the simplest extant vertebrate with sensory consciousness and qualia. Genes that pattern the proposed elements of consciousness (isomorphism, neural crest, placodes) have been identified in all vertebrates. Thus, consciousness is in the genes, some of which are already known.

Evolutionary relations of Hexanchiformes deep-sea sharks elucidated by whole mitochondrial genome sequences

Hexanchiformes is regarded as a monophyletic taxon, but the morphological and genetic relationships between the five extant species within the order are still uncertain. In this study, we determined the whole mitochondrial DNA (mtDNA) sequences of seven sharks including representatives of the five Hexanchiformes, one squaliform, and one carcharhiniform and inferred the phylogenetic relationships among those species and 12 other Chondrichthyes (cartilaginous fishes) species for which the complete mitogenome is available. The monophyly of Hexanchiformes and its close relation with all other Squaliformes sharks were strongly supported by likelihood and Bayesian phylogenetic analysis of 13,749 aligned nucleotides of 13 protein coding genes and two rRNA genes that were derived from the whole mDNA sequences of the 19 species. The phylogeny suggested that Hexanchiformes is in the superorder Squalomorphi, Chlamydoselachus anguineus (frilled shark) is the sister species to all other Hexanchiformes, and the relations within Hexanchiformes are well resolved as Chlamydoselachus, (Notorynchus, (Heptranchias, (Hexanchus griseus, H. nakamurai))). Based on our phylogeny, we discussed evolutionary scenarios of the jaw suspension mechanism and gill slit numbers that are significant features in the sharks.

Genome duplication in early vertebrates: insights from agnathan cytogenetics

Agnathans represent a remnant of a primitive offshoot of the vertebrates, and the long evolutionary separation between their 2 living groups, namely hagfishes and lampreys, could explain profound biological differences, also in karyotypes and genome sizes. Here, cytogenetic studies available on these vertebrates were summarized and data discussed with reference to the recently demonstrated monophyly of this group and to the 2 events of whole genome duplication (1R and 2R) characterizing the evolution of vertebrates. The comparison of cytogenetic data and phylogenetic relationships among agnathans and gnathostomes seems to support the hypothesis that 1R and 2R occurred before the evolutionary divergence between jawless and jawed vertebrates.

Revealing less derived nature of cartilaginous fish genomes with their evolutionary time scale inferred with nuclear genes

Cartilaginous fishes, divided into Holocephali (chimaeras) and Elasmoblanchii (sharks, rays and skates), occupy a key phylogenetic position among extant vertebrates in reconstructing their evolutionary processes. Their accurate evolutionary time scale is indispensable for better understanding of the relationship between phenotypic and molecular evolution of cartilaginous fishes. However, our current knowledge on the time scale of cartilaginous fish evolution largely relies on estimates using mitochondrial DNA sequences. In this study, making the best use of the still partial, but large-scale sequencing data of cartilaginous fish species, we estimate the divergence times between the major cartilaginous fish lineages employing nuclear genes. By rigorous orthology assessment based on available genomic and transcriptomic sequence resources for cartilaginous fishes, we selected 20 protein-coding genes in the nuclear genome, spanning 2973 amino acid residues. Our analysis based on the Bayesian inference resulted in the mean divergence time of 421 Ma, the late Silurian, for the Holocephali-Elasmobranchii split, and 306 Ma, the late Carboniferous, for the split between sharks and rays/skates. By applying these results and other documented divergence times, we measured the relative evolutionary rate of the Hox A cluster sequences in the cartilaginous fish lineages, which resulted in a lower substitution rate with a factor of at least 2.4 in comparison to tetrapod lineages. The obtained time scale enables mapping phenotypic and molecular changes in a quantitative framework. It is of great interest to corroborate the less derived nature of cartilaginous fish at the molecular level as a genome-wide phenomenon.

Molecular phylogeny of kinorhynchs

We reconstructed kinorhynch phylogeny using maximum-likelihood and Bayesian analyses of nuclear 18S and 28S rRNA gene sequences from 30 species in 13 genera (18S) and 23 species in 12 genera (28S), representing eight families and both orders (Cyclorhagida and Homalorhagida) currently recognized in the phylum. We analyzed the two genes individually (18S and 28S datasets) and in combination (18S+28S dataset). We detected four main clades (I-IV). Clade I consisted of family Echinoderidae. Clade II contained representatives of Zelinkaderidae, Antygomonidae, Semnoderidae, Centroderes, and Condyloderes, the latter two currently classified in Centroderidae; within Clade II, Zelinkaderidae, Antygomonidae, and Semnoderidae comprised a clade with strong nodal support. Clade III contained only two species in Campyloderes, also currently classified in the Centroderidae, indicating polyphyly for this family. Clades I-III, containing all representatives of Cyclorhagida included in the analysis except for Dracoderes abei, formed a clade with high nodal support in the 28S and 18S+28S trees. Clade IV, resolved in the 18S and 18S+28S trees with high nodal support, contained only species in order Homalorhagida, with the exception of the cyclorhagid Dracoderes abei. Order Cyclorhagida as it currently stands is thus polyphyletic, and order Homalorhagida paraphyletic. Our results indicate that Dracoderidae has been misplaced in Cyclorhagida based on homoplasious characters. Our analyses did not resolve the relationships among Clades I-III within Cyclorhagida. Neither gene alone nor the combined dataset resolved all nodes in trees, indicating that additional markers will be needed to reconstruct kinorhynch phylogeny.

Y chromosome phylogeny for cutthroat trout (Oncorhynchus clarkii) subspecies is generally concordant with those of other markers

Sequence divergence was evaluated in the non-recombining, male-specific OmyY1 region of the Y chromosome among the subspecies of cutthroat trout (Oncorhynchus clarkii) in the western United States. This evaluation identified subspecies-discriminating OmyY1-haplotypes within a ∼1200bp region of the OmyY1 locus and localized the region to the end of the Y chromosome by FISH analysis. OmyY1 sequences were aligned and used to reconstruct a phylogeny of the cutthroat trout subspecies and related species via maximum-parsimony and Bayesian analyses. In the Y-haplotype phylogeny, clade distributions generally corresponded to the geographic distributions of the recognized subspecies. This phylogeny generally corresponded to a mitochondrial tree obtained for these subspecies in a previous study. Both support a clade of trout vs. Pacific salmon, of rainbow trout, and of a Yellowstone cutthroat group within the cutthroat trout. In our OmyY1 tree, however, the cutthroat "clade", although present topologically, was not statistically significant. Some key differences were found between trees obtained from the paternally-inherited OmyY1 vs. maternally-inherited mitochondrial haplotypes in cutthroat trout compared to rainbow trout. Other findings are: The trout OmyY1 region evolves between 3 and 13 times slower than the trout mitochondrial regions that have been studied. The Lahontan cutthroat trout had a fixed OmyY1 sequence throughout ten separate populations, suggesting this subspecies underwent a severe population bottleneck prior to its current dispersal throughout the Great Basin during the pluvial phase of the last ice age. The Yellowstone group is the most derived among the cutthroat trout and consists of the Yellowstone, Bonneville, Colorado, Rio Grande and greenback subspecies. Identification of subspecies and sex with this Y-chromosome marker may prove useful in conservation efforts.

An analytical approach for estimating fossil record and diversification events in sharks, skates and rays

Background

Modern selachians and their supposed sister group (hybodont sharks) have a long and successful evolutionary history. Yet, although selachian remains are considered relatively common in the fossil record in comparison with other marine vertebrates, little is known about the quality of their fossil record. Similarly, only a few works based on specific time intervals have attempted to identify major events that marked the evolutionary history of this group.

Methodology/Principal Findings

Phylogenetic hypotheses concerning modern selachians’ interrelationships are numerous but differ significantly and no consensus has been found. The aim of the present study is to take advantage of the range of recent phylogenetic hypotheses in order to assess the fit of the selachian fossil record to phylogenies, according to two different branching methods. Compilation of these data allowed the inference of an estimated range of diversity through time and evolutionary events that marked this group over the past 300 Ma are identified. Results indicate that with the exception of high taxonomic ranks (orders), the selachian fossil record is by far imperfect, particularly for generic and post-Triassic data. Timing and amplitude of the various identified events that marked the selachian evolutionary history are discussed.

Conclusion/Significance

Some identified diversity events were mentioned in previous works using alternative methods (Early Jurassic, mid-Cretaceous, K/T boundary and late Paleogene diversity drops), thus reinforcing the efficiency of the methodology presented here in inferring evolutionary events. Other events (Permian/Triassic, Early and Late Cretaceous diversifications; Triassic/Jurassic extinction) are newly identified. Relationships between these events and paleoenvironmental characteristics and other groups’ evolutionary history are proposed.

Rates of molecular evolution vary in vertebrates for insulin-like growth factor-1 (IGF-1), a pleiotropic locus that regulates life history traits

Insulin-like growth factor-1 (IGF-1) is a member of the vertebrate insulin/insulin-like growth factor/relaxin gene family necessary for growth, reproduction, and survival at both the cellular and organismal level. Its sequence, protein structure, and function have been characterized in mammals, birds, and fish; however, a notable gap in our current knowledge of the function of IGF-1 and its molecular evolution is information in ectothermic reptiles. To address this disparity, we sequenced the coding region of IGF-1 in 11 reptile species-one crocodilian, three turtles, three lizards, and four snakes. Complete sequencing of the full mRNA transcript of a snake revealed the Ea-isoform, the predominant isoform of IGF-1 also reported in other vertebrate groups. A gene tree of the IGF-1 protein-coding region that incorporated sequences from diverse vertebrate groups showed similarity to the species phylogeny, with the exception of the placement of Testudines as sister group to Aves, due to their high nucleotide sequence similarity. In contrast, long-branch lengths indicate more rapid divergence in IGF-1 among lizards and snakes. Additionally, lepidosaurs (i.e., lizards and snakes) had higher rates of non-synonymous:synonymous substitutions (dN/dS) relative to archosaurs (i.e., birds and crocodilians) and turtles. Tests for positive selection on specific codons within branches and evaluation of the changes in the amino acid properties, suggested positive selection in lepidosaurs on the C domain of IGF-1, which is involved in binding affinity to the IGF-1 receptor. Predicted structural changes suggest that major alterations in protein structure and function may have occurred in reptiles. These data propose new insights into the molecular co-evolution of IGF-1 and its receptors, and ultimately the evolution of IGF-1's role in regulating life-history traits across vertebrates.

Plated Cambrian bilaterians reveal the earliest stages of echinoderm evolution

Echinoderms are unique in being pentaradiate, having diverged from the ancestral bilaterian body plan more radically than any other animal phylum. This transformation arises during ontogeny, as echinoderm larvae are initially bilateral, then pass through an asymmetric phase, before giving rise to the pentaradiate adult. Many fossil echinoderms are radial and a few are asymmetric, but until now none have been described that show the original bilaterian stage in echinoderm evolution. Here we report new fossils from the early middle Cambrian of southern Europe that are the first echinoderms with a fully bilaterian body plan as adults. Morphologically they are intermediate between two of the most basal classes, the Ctenocystoidea and Cincta. This provides a root for all echinoderms and confirms that the earliest members were deposit feeders not suspension feeders.

The rhizome of life: what about metazoa?

The increase in huge number of genomic sequences in recent years has contributed to various genetic events such as horizontal gene transfer (HGT), gene duplication and hybridization of species. Among them HGT has played an important role in the genome evolution and was believed to occur only in Bacterial and Archaeal genomes. As a result, genomes were found to be chimeric and the evolution of life was represented in different forms such as forests, networks and species evolution was described more like a rhizome, rather than a tree. However, in the last few years, HGT has also been evidenced in other group such as metazoa (for example in root-knot nematodes, bdelloid rotifers and mammals). In addition to HGT, other genetic events such as transfer by retrotransposons and hybridization between more closely related lineages are also well established. Therefore, in the light of such genetic events, whether the evolution of metazoa exists in the form of a tree, network or rhizome is highly questionable and needs to be determined. In the current review, we will focus on the role of HGT, retrotransposons and hybridization in the metazoan evolution.

Adding resolution to ordinal level relationships of tapeworms (Platyhelminthes: Cestoda) with large fragments of mtDNA

The construction of a stable phylogeny for the Cestoda, indicating the interrelationships of recognised orders and other major lineages, has proceeded iteratively since the group first received attention from phylogenetic systematists. Molecular analyses using nuclear ribosomal RNA gene fragments from the small (ssrDNA) and large (lsrDNA) subunits have been used to test competing evolutionary scenarios based on morphological data but could not arbitrate between some key conflicting hypotheses. To the ribosomal data, we have added a contiguous fragment of mitochondrial (mt) genome data (mtDNA) of partial nad1-trnN-trnP-trnI-trnK-nad3-trnS-trnW-cox1-trnT-rrnL-trnC-partial rrnS, spanning 4034-4447 bp, where new data for this region were generated for 18 species. Bayesian analysis of mtDNA and rDNA as nucleotides, and where appropriate as amino acids, demonstrated that these two classes of genes provide complementary signal across the phylogeny. In all analyses, except when using mt amino acids only, the Gyrocotylidea is sister group to all other Cestoda (Nephroposticophora), and Amphilinidea forms the sister group to the Eucestoda. However, an earliest-diverging position of Amphilinidea is strongly supported in the mt amino acid analysis. Amphilinidea exhibit a unique tRNA arrangement (nad1-trnI-trnL2-trnP-trnK-trnV-trnA-trnN-nad3), whereas Gyrocotylidea shares that of the derived lineages, providing additional evidence of the uniqueness of amphilinid genes and genomes. The addition of mtDNA to the rDNA genes supported the Caryophyllidea as the sister group to (Spathebothriidea+remaining Eucestoda), a hypothesis consistently supported by morphology. This relationship suggests a history of step-wise evolutionary transitions from simple monozoic, unsegmented tapeworms to the more familiar polyzoic, externally segmented (strobilate) forms. All our data partitions recovered Haplobothriidea as the sister group to Diphyllobothriidae. The sister-group relationship between Diphyllidea and Trypanorhyncha, as previously established using rDNA, is not supported by the mt data, although it is supported by the combined mt and rDNA analysis. With regards to the more derived taxa, in all except the mt amino acid analysis, the following topology is supported: (Bothriocephalidea (Litobothriidea (Lecanicephalidea (Rhinebothriidea (Tetraphyllidea, (Acanthobothrium, Proteocephalidea), (Nippotaeniidea, Mesocestoididae, Tetrabothriidea, Cyclophyllidea)))))), where the Tetraphyllidea are paraphyletic. Evidence from the mt data provides strong (nucleotides) to moderate (amino acids) support for Tetraphyllidea forming a group to the inclusion of Proteocephalidea, with the latter consistently forming the sister group to Acanthobothrium. The interrelationships among Nippotaeniidea, Mesocestoididae, Tetrabothriidea and Cyclophyllidea remain ambiguous and require further systematic attention. Mitochondrial and nuclear rDNA data provide conflicting signal for certain parts of the cestode tree. In some cases mt data offer results in line with morphological evidence, such as the interrelationships of the early divergent lineages. Also, Tetraphyllidea, although remaining paraphyletic with the inclusion of the Proteocephalidea, does not include the most derived cestodes; a result which has consistently been obtained with rDNA.

Phylogenetics of Chondrichthyes and the problem of rooting phylogenies with distant outgroups

Erroneous estimates of ingroup relationships can be caused by attributes in the outgroup chosen to root the tree. Phylogenetic analyses of DNA sequences frequently yield incorrect estimates of ingroup relationships when the outgroup used to "root" the tree is highly divergent from the ingroup. This is especially the case when the outgroup has a different base composition than the ingroup. Unfortunately, in many instances, alternative less divergent outgroups are not available. In such cases, investigators must either target genes with attributes that minimize the problem (slowly evolving genes with stationary base compositions--which are often not ideal for estimating relationships among the more closely related ingroup taxa) or use inference models that are explicitly tailored to deal with an attenuated historical signal with a superimposed non-stationary base composition. In this paper we explore the problem both empirically and through simulation. For the empirical component we looked at the phylogenetic relationships among elasmobranch fishes (sharks and rays), a group whose closest living outgroup, the holocephalan Ghost fishes, are separated from the elasmobranchs by more than 100 million years of evolution. We compiled a data set for analysis comprising 10 single-copy nuclear protein-coding genes (12,096 bp) for representatives of the major lineages within elasmobranchs and holocephalans. For the simulation, we used an evolutionary model on a fixed tree topology to generate DNA sequence data sets which varied both in their distance to the outgroup, and in their base compositional difference between ingroup and outgroup. Results from both the empirical data set and the simulation, support the idea that deviation from base compositional stationarity, in conjunction with distance from the root can act in concert to compromise accuracy of estimated relationships within the ingroup. We tested several approaches to mitigate such problems. We found, that excluding genes with overall faster rates and heterogeneous base compositions, while the least sophisticated of the methods evaluated, seemed to be the most effective.

Resolving phylogenetic signal from noise when divergence is rapid: a new look at the old problem of echinoderm class relationships

Resolving evolutionary relationships in groups that underwent fast radiation in deep time is a problem for molecular phylogeny, as the scant phylogenetic signal that characterises short internal branches is generally swamped by more recent substitutions. We implement an approach, that maps how the support for rival phylogenies changes when analysing subsets of sites with either faster and more heterogeneous rates or slower and more homogeneous rates, to address a long-standing problem in deuterostome phylogeny - the interrelationships of the eleutherozoan echinoderm classes. We show that miRNA genes are phylogenetically uninformative as to the relationships of asteroids, echinoids and ophiuroids, consistent with a rapid radiation of these groups as suggested by their fossil record. Using three nuclear rRNAs and seven nuclear housekeeping genes, we map the support for the three possible phylogenetic arrangements of asteroids, ophiuroids and echinoids when moving between subsets of the data with very similar or very different rates of evolution. Only one of the three possible topologies (asteroids (ophiuroids+echinoids)) strengthens when the most rate-homogeneous subset of data are analysed. The other two possible pairings become stronger in a less reliable data subset, which includes the fastest and thus homoplasy-rich data in our alignment. Thus, while superficial analysis of our concatenated alignment identifies asteroids and ophiuroids as sister taxa, more thorough analyses suggest that ophiuroids may be more closely related to echinoids. Divergence of these echinoderm groups, using a relaxed molecular clock, is estimated to have occurred within ≈ 5 million years. Our results illustrate that the analytic approach of phylogenetic signal dissection can be a powerful tool to investigate rapid radiations in deep geologic time.

Evolution of invertebrate deuterostomes and Hox/ParaHox genes

Transcription factors encoded by Antennapedia-class homeobox genes play crucial roles in controlling development of animals, and are often found clustered in animal genomes. The Hox and ParaHox gene clusters have been regarded as evolutionary sisters and evolved from a putative common ancestral gene complex, the ProtoHox cluster, prior to the divergence of the Cnidaria and Bilateria (bilaterally symmetrical animals). The Deuterostomia is a monophyletic group of animals that belongs to the Bilateria, and a sister group to the Protostomia. The deuterostomes include the vertebrates (to which we belong), invertebrate chordates, hemichordates, echinoderms and possibly xenoturbellids, as well as acoelomorphs. The studies of Hox and ParaHox genes provide insights into the origin and subsequent evolution of the bilaterian animals. Recently, it becomes apparent that among the Hox and ParaHox genes, there are significant variations in organization on the chromosome, expression pattern, and function. In this review, focusing on invertebrate deuterostomes, I first summarize recent findings about Hox and ParaHox genes. Next, citing unsolved issues, I try to provide clues that might allow us to reconstruct the common ancestor of deuterostomes, as well as understand the roles of Hox and ParaHox genes in the development and evolution of deuterostomes.

Phylogenetic aspects of carbamoyl phosphate synthetase in lungfish: a transitional enzyme in transitional fishes

Carbamoyl phosphate synthetase (CPS) catalyses the formation of carbamoyl phosphate from glutamine or ammonia, bicarbonate and ATP. There are three different isoforms of CPS that play vital roles in two metabolic pathways, pyrimidine biosynthesis (CPS II) and arginine/urea biosynthesis (CPS I and CPS III). Gene duplication has been proposed as the evolutionary mechanism creating this gene family with CPS II likely giving rise to the CPS I/III clade. In the evolutionary history of this gene family it is still undetermined when CPS I diverged from CPS III on the path to terrestriality in the vertebrates. Transitional organisms such as lungfishes are of particular interest because they are capable of respiring via gills and with lungs and therefore can be found in both aquatic and terrestrial environments. Notably, enzymatic characterization of the mitochondrial CPS isoforms in this transitional group has not led to clear conclusions. In order to determine which CPS isoform is present in transitional animals, we examined partial sequences for liver CPS amplified from five species of lungfish, and a larger fragment of CPS from one lungfish species (Protopterus annectens) and compared them to CPS isoforms from other fish and mammals. Enzyme activities for P. annectens liver were also examined. While enzyme activities did not yield a clear distinction between isoforms (virtually equal activities were obtained for either CPS I or III), CPS sequences from the lungfishes formed a monophyletic clade within the CPS I clade and separate from the CPS III clade of other vertebrates. This finding implies that the mitochondrial isoform of CPS in lungfish is derived from CPS I and is likely to have a physiological function similar to CPS I. This finding is important because it supports the hypothesis that lungfish employ a urea cycle similar to terrestrial air-breathing vertebrates.

Environmental sex determination in the branchiopod crustacean Daphnia magna: deep conservation of a Doublesex gene in the sex-determining pathway

Sex-determining mechanisms are diverse among animal lineages and can be broadly divided into two major categories: genetic and environmental. In contrast to genetic sex determination (GSD), little is known about the molecular mechanisms underlying environmental sex determination (ESD). The Doublesex (Dsx) genes play an important role in controlling sexual dimorphism in genetic sex-determining organisms such as nematodes, insects, and vertebrates. Here we report the identification of two Dsx genes from Daphnia magna, a freshwater branchiopod crustacean that parthenogenetically produces males in response to environmental cues. One of these genes, designated DapmaDsx1, is responsible for the male trait development when expressed during environmental sex determination. The domain organization of DapmaDsx1 was similar to that of Dsx from insects, which are thought to be the sister group of branchiopod crustaceans. Intriguingly, the molecular basis for sexually dimorphic expression of DapmaDsx1 is different from that of insects. Rather than being regulated sex-specifically at the level of pre–mRNA splicing in the coding region, DapmaDsx1 exhibits sexually dimorphic differences in the abundance of its transcripts. During embryogenesis, expression of DapmaDsx1 was increased only in males and its transcripts were primarily detected in male-specific structures. Knock-down of DapmaDsx1 in male embryos resulted in the production of female traits including ovarian maturation, whereas ectopic expression of DapmaDsx1 in female embryos resulted in the development of male-like phenotypes. Expression patterns of another D. magna Dsx gene, DapmaDsx2, were similar to those of DapmaDsx1, but silencing and overexpression of this gene did not induce any clear phenotypic changes. These results establish DapmaDsx1 as a key regulator of the male phenotype. Our findings reveal how ESD is implemented by selective expression of a fundamental genetic component that is functionally conserved in animals using GSD. We infer that there is an ancient, previously unidentified link between genetic and environmental sex determination.

Sex determination is a fundamental biological process that can be broadly divided into two major categories. In genetic sex determination (GSD), sex-specific differentiation results from intrinsic genetic differences between males and females, whereas environmental sex determination (ESD) relies on environmental signals to induce male or female sex determination. In contrast to model organisms that utilize GSD system, environmental sex-determining organisms are poor genetic models. Therefore, although candidate genes involved in ESD have been found in vertebrates, their functions have remained largely unknown, impairing our understanding of ESD and making the comparison of sex-determining genes between both systems difficult. Here, we report the identification of a gene responsible for the production of males during environmental sex determination in the crustacean Daphnia. This gene is homologous to the Doublesex gene that is functionally conserved in animals that use GSD. Expression of Doublesex was increased primarily in male-specific structures. Gain- and loss-of-function analyses established that Daphnia Doublesex gene is a major effector that regulates the male phenotype in Daphnia. We infer that there is an ancient, previously unidentified link between genetic and environmental sex determination.