Evolution of Hox gene clusters in gnathostomes: insights from a survey of a shark (Scyliorhinus canicula) transcriptome

The sensory shark: high-quality morphological, genomic and transcriptomic data for the small-spotted catshark Scyliorhinus canicula reveal the molecular bases of sensory organ evolution in jawed vertebrates

Cartilaginous fishes (chimaeras and elasmobranchs -sharks, skates and rays) hold a key phylogenetic position to explore the origin and diversifications of jawed vertebrates. Here, we report and integrate reference genomic, transcriptomic and morphological data in the small-spotted catshark Scyliorhinus canicula to shed light on the evolution of sensory organs. We first characterise general aspects of the catshark genome, confirming the high conservation of genome organisation across cartilaginous fishes, and investigate population genomic signatures. Taking advantage of a dense sampling of transcriptomic data, we also identify gene signatures for all major organs, including chondrichthyan specializations, and evaluate expression diversifications between paralogs within major gene families involved in sensory functions. Finally, we combine these data with 3D synchrotron imaging and in situ gene expression analyses to explore chondrichthyan-specific traits and more general evolutionary trends of sensory systems. This approach brings to light, among others, novel markers of the ampullae of Lorenzini electro-sensory cells, a duplication hotspot for crystallin genes conserved in jawed vertebrates, and a new metazoan clade of the Transient-receptor potential (TRP) family. These resources and results, obtained in an experimentally tractable chondrichthyan model, open new avenues to integrate multiomics analyses for the study of elasmobranchs and jawed vertebrates.

Interference with the retinoic acid signalling pathway inhibits the initiation of teeth and caudal primary scales in the small-spotted catshark Scyliorhinus canicula

The retinoic acid (RA) pathway was shown to be important for tooth development in mammals, and suspected to play a key role in tooth evolution in teleosts. The general modalities of development of tooth and "tooth-like" structures (collectively named odontodes) seem to be conserved among all jawed vertebrates, both with regard to histogenesis and genetic regulation. We investigated the putative function of RA signalling in tooth and scale initiation in a cartilaginous fish, the small-spotted catshark Scyliorhinus canicula. To address this issue, we identified the expression pattern of genes from the RA pathway during both tooth and scale development and performed functional experiments by exposing small-spotted catshark embryos to exogenous RA or an inhibitor of RA synthesis. Our results showed that inhibiting RA synthesis affects tooth but not caudal primary scale development while exposure to exogenous RA inhibited both. We also showed that the reduced number of teeth observed with RA exposure is probably due to a specific inhibition of tooth bud initiation while the observed effects of the RA synthesis inhibitor is related to a general delay in embryonic development that interacts with tooth development. This study provides data complementary to previous studies of bony vertebrates and support an involvement of the RA signalling pathway toolkit in odontode initiation in all jawed vertebrates. However, the modalities of RA signalling may vary depending on the target location along the body, and depending on the species lineage.

HOX-Gene Cluster Organization and Genome Duplications in Fishes and Mammals: Transcript Variant Distribution along the Anterior–Posterior Axis

Hox genes play a crucial role in morphogenesis, especially in anterior–posterior body axis patterning. The organization of Hox clusters in vertebrates is a result of several genome duplications: two rounds of duplication in the ancestors of all vertebrates and a third round that was specific for teleost fishes. Teleostei cluster structure has been significantly modified in the evolutionary processes by Hox gene losses and co-options, while mammals show no such tendency. In mammals, the Hox gene number in a single cluster is stable and generally large, and the numbers are similar to those in the Chondrichthyes. Hox gene alternative splicing activity slightly differs between fishes and mammals. Fishes and mammals have differences in their known alternative splicing activity for Hox gene distribution along the anterior–posterior body axis. The analyzed fish groups—the Coelacanthiformes, Chondrichthyes, and Teleostei—all have higher known alternative mRNA numbers from the anterior and posterior regions, whereas mammals have a more uniform Hox transcript distribution along this axis. In fishes, most Hox transcripts produce functioning proteins, whereas mammals have significantly more known transcripts that do not produce functioning proteins.

Evolutionary Transition in the Regulation of Vertebrate Pronephros Development: A New Role for Retinoic Acid

The anterior-posterior (AP) axis in chordates is regulated by a conserved set of genes and signaling pathways, including Hox genes and retinoic acid (RA), which play well-characterized roles in the organization of the chordate body plan. The intermediate mesoderm (IM), which gives rise to all vertebrate kidneys, is an example of a tissue that differentiates sequentially along this axis. Yet, the conservation of the spatiotemporal regulation of the IM across vertebrates remains poorly understood. In this study, we used a comparative developmental approach focusing on non-conventional model organisms, a chondrichthyan (catshark), a cyclostome (lamprey), and a cephalochordate (amphioxus), to assess the involvement of RA in the regulation of chordate and vertebrate pronephros formation. We report that the anterior expression boundary of early pronephric markers (Pax2 and Lim1), positioned at the level of somite 6 in amniotes, is conserved in the catshark and the lamprey. Furthermore, RA, driving the expression of Hox4 genes like in amniotes, regulates the anterior pronephros boundary in the catshark. We find no evidence for the involvement of this regulatory hierarchy in the AP positioning of the lamprey pronephros and the amphioxus pronephros homolog, Hatschek's nephridium. This suggests that despite the conservation of Pax2 and Lim1 expressions in chordate pronephros homologs, the responsiveness of the IM, and hence of pronephric genes, to RA- and Hox-dependent regulation is a gnathostome novelty.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

The Roles of lncRNA in Cutaneous Squamous Cell Carcinoma

Cutaneous squamous cell carcinoma derives from keratinocytes and is the second most common cause of non-melanoma skin cancer. Cutaneous squamous cell carcinoma (cSCC) develops rapidly and is also the leading cause of death in non-melanoma cancers. Lymph node metastasis occurs in 5% of cSCC patients, and some patients may even metastasize to the viscera. Patients with regional lymphatic metastasis or distant metastases have a <20% 10-year survival rate, indicating the substantial challenge in treating advanced and metastatic cSCC. Some lncRNAs have been found to be abnormally overexpressed in many tumor tissues, so that they can be considered as potential new biomarkers or targets that can be used in the diagnosis and treatment of cSCC in the future. In this review, we summarize the role of lncRNA in cutaneous squamous cell carcinoma to make a better understanding of mutations in cSCC and lay the foundation for effective target therapy of cSCC.

Skeletal Mineralization in Association with Type X Collagen Expression Is an Ancestral Feature for Jawed Vertebrates

In order to characterize the molecular bases of mineralizing cell evolution, we targeted type X collagen, a nonfibrillar network forming collagen encoded by the Col10a1 gene. It is involved in the process of endochondral ossification in ray-finned fishes and tetrapods (Osteichthyes), but until now unknown in cartilaginous fishes (Chondrichthyes). We show that holocephalans and elasmobranchs have respectively five and six tandemly duplicated Col10a1 gene copies that display conserved genomic synteny with osteichthyan Col10a1 genes. All Col10a1 genes in the catshark Scyliorhinus canicula are expressed in ameloblasts and/or odontoblasts of teeth and scales, during the stages of extracellular matrix protein secretion and mineralization. Only one duplicate is expressed in the endoskeletal (vertebral) mineralizing tissues. We also show that the expression of type X collagen is present in teeth of two osteichthyans, the zebrafish Danio rerio and the western clawed frog Xenopus tropicalis, indicating an ancestral jawed vertebrate involvement of type X collagen in odontode formation. Our findings push the origin of Col10a1 gene prior to the divergence of osteichthyans and chondrichthyans, and demonstrate its ancestral association with mineralization of both the odontode skeleton and the endoskeleton.

Asymmetric paralog evolution between the "cryptic" gene Bmp16 and its well-studied sister genes Bmp2 and Bmp4

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.

The role of HoxA11 and HoxA13 in the evolution of novel fin morphologies in a representative batoid (Leucoraja erinacea)

Background

Batoids exhibit unique body plans with derived fin morphologies, such as the anteriorly expanded pectoral fins that fuse to the head, or distally extended anterior pelvic fin lobes used for a modified swimming technique utilized by skates (Rajidae). The little skate (Leucoraja erinacea), exhibits both of these unique fin morphologies. These fin modifications are not present in a typical shark body plan, and little is known regarding the mechanisms underlying their development. A recent study identified a novel apical ectodermal ridge (AER) associated with the development of the anterior pectoral fin in the little skate, but the role of the posterior HoxA genes was not featured during skate fin development.

Results

We present the first evidence for HoxA expression (HoxA11 and HoxA13) in novel AER domains associated with the development of three novel fin morphologies in a representative batoid, L. erinacea. We found HoxA13 expression associated with the recently described novel AER in the anterior pectoral fin, and HoxA11 expression in a novel AER domain in the anterior pelvic fin that we describe here. We find that both HoxA11 and HoxA13 are expressed in claspers, and while HoxA11 is expressed in pelvic fins and claspers, HoxA13 is expressed exclusively in developing claspers of males. Finally, HoxA11 expression is associated with the developing fin rays in paired fins.

Conclusion

Overall, these results indicate that the posterior HoxA genes play an important role in the morphological evolution of paired fins in a representative batoid. These data suggest that the batoids utilize a unique Hox code, where the posterior HoxA genes exhibit distinct expression patterns that are likely associated with specification of novel fin morphologies.

Electronic supplementary material

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

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy

BACKGROUND

Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy.

RESULTS

Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism.

CONCLUSIONS

This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Functional and Comparative Genomics of Hoxa2 Gene cis-Regulatory Elements: Evidence for Evolutionary Modification of Ancestral Core Element Activity

Hoxa2 is an evolutionarily conserved developmental regulatory gene that functions to specify rhombomere (r) and pharyngeal arch (PA) identities throughout the Osteichthyes. Japanese medaka (Oryzias latipes) hoxa2a, like orthologous Hoxa2 genes from other osteichthyans, is expressed during embryogenesis in r2–7 and PA2-7, whereas the paralogous medaka pseudogene, ψhoxa2b, is expressed in noncanonical Hoxa2 domains, including the pectoral fin buds. To understand the evolution of cis-regulatory element (CRE) control of gene expression, we conducted eGFP reporter gene expression studies with extensive functional mapping of several conserved CREs upstream of medaka hoxa2a and ψhoxa2b in transient and stable-line transgenic medaka embryos. The CREs tested were previously shown to contribute to directing mouse Hoxa2 gene expression in r3, r5, and PA2-4. Our results reveal the presence of sequence elements embedded in the medaka hoxa2a and ψhoxa2b upstream enhancer regions (UERs) that mediate expression in r4 and the PAs (hoxa2a r4/CNCC element) or in r3–7 and the PAs ψhoxa2b r3–7/CNCC element), respectively. Further, these elements were shown to be highly conserved among osteichthyans, which suggests that the r4 specifying element embedded in the UER of Hoxa2 is a deeply rooted rhombomere specifying element and the activity of this element has been modified by the evolution of flanking sequences that redirect its activity to alternative developmental compartments.

Tooth and scale morphogenesis in shark: an alternative process to the mammalian enamel knot system

BACKGROUND

The gene regulatory network involved in tooth morphogenesis has been extremely well described in mammals and its modeling has allowed predictions of variations in regulatory pathway that may have led to evolution of tooth shapes. However, very little is known outside of mammals to understand how this regulatory framework may also account for tooth shape evolution at the level of gnathostomes. In this work, we describe expression patterns and proliferation/apoptosis assays to uncover homologous regulatory pathways in the catshark Scyliorhinus canicula.

RESULTS

Because of their similar structural and developmental features, gene expression patterns were described over the four developmental stages of both tooth and scale buds in the catshark. These gene expression patterns differ from mouse tooth development, and discrepancies are also observed between tooth and scale development within the catshark. However, a similar nested expression of Shh and Fgf suggests similar signaling involved in morphogenesis of all structures, although apoptosis assays do not support a strictly equivalent enamel knot system in sharks. Similarities in the topology of gene expression pattern, including Bmp signaling pathway, suggest that mouse molar development is more similar to scale bud development in the catshark.

CONCLUSIONS

These results support the fact that no enamel knot, as described in mammalian teeth, can be described in the morphogenesis of shark teeth or scales. However, homologous signaling pathways are involved in growth and morphogenesis with variations in their respective expression patterns. We speculate that variations in this topology of expression are also a substrate for tooth shape evolution, notably in regulating the growth axis and symmetry of the developing structure.

Molecular footprinting of skeletal tissues in the catshark Scyliorhinus canicula and the clawed frog Xenopus tropicalis identifies conserved and derived features of vertebrate calcification

Understanding the evolutionary emergence and subsequent diversification of the vertebrate skeleton requires a comprehensive view of the diverse skeletal cell types found in distinct developmental contexts, tissues, and species. To date, our knowledge of the molecular nature of the shark calcified extracellular matrix, and its relationships with osteichthyan skeletal tissues, remain scarce. Here, based on specific combinations of expression patterns of the Col1a1, Col1a2, and Col2a1 fibrillar collagen genes, we compare the molecular footprint of endoskeletal elements from the chondrichthyan Scyliorhinus canicula and the tetrapod Xenopus tropicalis. We find that, depending on the anatomical location, Scyliorhinus skeletal calcification is associated to cell types expressing different subsets of fibrillar collagen genes, such as high levels of Col1a1 and Col1a2 in the neural arches, high levels of Col2a1 in the tesserae, or associated to a drastic Col2a1 downregulation in the centrum. We detect low Col2a1 levels in Xenopus osteoblasts, thereby revealing that the osteoblastic expression of this gene was significantly reduced in the tetrapod lineage. Finally, we uncover a striking parallel, from a molecular and histological perspective, between the vertebral cartilage calcification of both species and discuss the evolutionary origin of endochondral ossification.

[ENCODE apophenia or a panglossian analysis of the human genome]

In September 2012, a batch of more than 30 articles presenting the results of the ENCODE (Encyclopaedia of DNA Elements) project was released. Many of these articles appeared in Nature and Science, the two most prestigious interdisciplinary scientific journals. Since that time, hundreds of other articles dedicated to the further analyses of the Encode data have been published. The time of hundreds of scientists and hundreds of millions of dollars were not invested in vain since this project had led to an apparent paradigm shift: contrary to the classical view, 80% of the human genome is not junk DNA, but is functional. This hypothesis has been criticized by evolutionary biologists, sometimes eagerly, and detailed refutations have been published in specialized journals with impact factors far below those that published the main contribution of the Encode project to our understanding of genome architecture. In 2014, the Encode consortium released a new batch of articles that neither suggested that 80% of the genome is functional nor commented on the disappearance of their 2012 scientific breakthrough. Unfortunately, by that time many biologists had accepted the idea that 80% of the genome is functional, or at least, that this idea is a valid alternative to the long held evolutionary genetic view that it is not. In order to understand the dynamics of the genome, it is necessary to re-examine the basics of evolutionary genetics because, not only are they well established, they also will allow us to avoid the pitfall of a panglossian interpretation of Encode. Actually, the architecture of the genome and its dynamics are the product of trade-offs between various evolutionary forces, and many structural features are not related to functional properties. In other words, evolution does not produce the best of all worlds, not even the best of all possible worlds, but only one possible world.

Genome-wide analysis of homeobox genes from Mesobuthus martensii reveals Hox gene duplication in scorpions

Homeobox genes belong to a large gene group, which encodes the famous DNA-binding homeodomain that plays a key role in development and cellular differentiation during embryogenesis in animals. Here, one hundred forty-nine homeobox genes were identified from the Asian scorpion, Mesobuthus martensii (Chelicerata: Arachnida: Scorpiones: Buthidae) based on our newly assembled genome sequence with approximately 248 × coverage. The identified homeobox genes were categorized into eight classes including 82 families: 67 ANTP class genes, 33 PRD genes, 11 LIM genes, five POU genes, six SINE genes, 14 TALE genes, five CUT genes, two ZF genes and six unclassified genes. Transcriptome data confirmed that more than half of the genes were expressed in adults. The homeobox gene diversity of the eight classes is similar to the previously analyzed Mandibulata arthropods. Interestingly, it is hypothesized that the scorpion M. martensii may have two Hox clusters. The first complete genome-wide analysis of homeobox genes in Chelicerata not only reveals the repertoire of scorpion, arachnid and chelicerate homeobox genes, but also shows some insights into the evolution of arthropod homeobox genes.

The coelacanth: Can a "living fossil" have active transposable elements in its genome?

The coelacanth has long been regarded as a "living fossil," with extant specimens looking very similar to fossils dating back to the Cretaceous period. The hypothesis of a slowly or even not evolving genome has been proposed to account for this apparent morphological stasis. While this assumption seems to be sustained by different evolutionary analyses on protein-coding genes, recent studies on transposable elements have provided more conflicting results. Indeed, the coelacanth genome contains many transposable elements and has been shaped by several major bursts of transposition during evolution. In addition, comparison of orthologous genomic regions from the genomes of the 2 extant coelacanth species L. chalumnae and L. menadoensis revealed multiple species-specific insertions, indicating transposable element recent activity and contribution to post-speciation genome divergence. These observations, which do not support the genome stasis hypothesis, challenge either the impact of transposable elements on organismal evolution or the status of the coelacanth as a "living fossil." Closer inspection of fossil and molecular data indicate that, even if coelacanths might evolve more slowly than some other lineages due to demographic and/or ecological factors, this variation is still in the range of a "non-fossil" vertebrate species.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.

SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Chondrichthyan fishes are a diverse class of gnathostomes that provide a valuable perspective on fundamental characteristics shared by all jawed and limbed vertebrates. Studies of phylogeny, species diversity, population structure, conservation, and physiology are accelerated by genomic, transcriptomic and protein sequence data. These data are widely available for many sarcopterygii (coelacanth, lungfish and tetrapods) and actinoptergii (ray-finned fish including teleosts) taxa, but limited for chondrichthyan fishes.  In this study, we summarize available data for chondrichthyes and describe resources for one of the largest projects to characterize one of these fish, Leucoraja erinacea, the little skate.  SkateBase ( http://skatebase.org) serves as the skate genome project portal linking data, research tools, and teaching resources.

Evolution of Hox gene clusters in deuterostomes

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

Heterogeneous conservation of Dlx paralog co-expression in jawed vertebrates

Background

The Dlx gene family encodes transcription factors involved in the development of a wide variety of morphological innovations that first evolved at the origins of vertebrates or of the jawed vertebrates. This gene family expanded with the two rounds of genome duplications that occurred before jawed vertebrates diversified. It includes at least three bigene pairs sharing conserved regulatory sequences in tetrapods and teleost fish, but has been only partially characterized in chondrichthyans, the third major group of jawed vertebrates. Here we take advantage of developmental and molecular tools applied to the shark Scyliorhinus canicula to fill in the gap and provide an overview of the evolution of the Dlx family in the jawed vertebrates. These results are analyzed in the theoretical framework of the DDC (Duplication-Degeneration-Complementation) model.

Results

The genomic organisation of the catshark Dlx genes is similar to that previously described for tetrapods. Conserved non-coding elements identified in bony fish were also identified in catshark Dlx clusters and showed regulatory activity in transgenic zebrafish. Gene expression patterns in the catshark showed that there are some expression sites with high conservation of the expressed paralog(s) and other expression sites with events of paralog sub-functionalization during jawed vertebrate diversification, resulting in a wide variety of evolutionary scenarios within this gene family.

Conclusion

Dlx gene expression patterns in the catshark show that there has been little neo-functionalization in Dlx genes over gnathostome evolution. In most cases, one tandem duplication and two rounds of vertebrate genome duplication have led to at least six Dlx coding sequences with redundant expression patterns followed by some instances of paralog sub-functionalization. Regulatory constraints such as shared enhancers, and functional constraints including gene pleiotropy, may have contributed to the evolutionary inertia leading to high redundancy between gene expression patterns.

Development of head and trunk mesoderm in the dogfish, Scyliorhinus torazame: I. Embryology and morphology of the head cavities and related structures

Vertebrate head segmentation has attracted the attention of comparative and evolutionary morphologists for centuries, given its importance for understanding the developmental body plan of vertebrates and its evolutionary origin. In particular, the segmentation of the mesoderm is central to the problem. The shark embryo has provided a canonical morphological scheme of the head, with its epithelialized coelomic cavities (head cavities), which have often been regarded as head somites. To understand the evolutionary significance of the head cavities, the embryonic development of the mesoderm was investigated at the morphological and histological levels in the shark, Scyliorhinus torazame. Unlike somites and some enterocoelic mesodermal components in other vertebrates, the head cavities in S. torazame appeared as irregular cyst(s) in the originally unsegmented mesenchymal head mesoderm, and not via segmentation of an undivided coelom. The mandibular cavity appeared first in the paraxial part of the mandibular mesoderm, followed by the hyoid cavity, and the premandibular cavity was the last to form. The prechordal plate was recognized as a rhomboid roof of the preoral gut, continuous with the rostral notochord, and was divided anteroposteriorly into two parts by the growth of the hypothalamic primordium. Of those, the posterior part was likely to differentiate into the premandibular cavity, and the anterior part disappeared later. The head cavities and somites in the trunk exhibited significant differences, in terms of histological appearance and timing of differentiation. The mandibular cavity developed a rostral process secondarily; its homology to the anterior cavity reported in some elasmobranch embryos is discussed.

[Tree thinking: the evolution of vertebrate as a case study]

In the last decades, the reconstruction of phylogenetic relationships and their representation in the form of a phylogenetic tree have become a powerful tool for biology. These methods involve the use of sophisticated computer programs that are unfamiliar to most of biologists. However, our experience as teacher and researcher in evolutionary biology prompted us to realize that the main and most common source of confusion in depicting the evolution of various traits comes from the misunderstanding of the basis of tree reading. We have identified, not only in the work of our students, but also in scientific literature, some common mistakes that reveal the persistency of the concept of the scale of beings that unfortunately maintained as frame of reference to analyse phylogenies.

The Dlx genes as clues to vertebrate genomics and craniofacial evolution

The group of Dlx genes belongs to the homeobox-containing superfamily, and its members are involved in various morphogenetic processes. In vertebrate genomes, Dlx genes exist as multiple paralogues generated by tandem duplication followed by whole genome duplications. In this review, we provide an overview of the Dlx gene phylogeny with an emphasis on the chordate lineage. Referring to the Dlx gene repertoire, we discuss the establishment and conservation of the nested expression patterns of the Dlx genes in craniofacial development. Despite the accumulating genomic sequence resources in diverse vertebrates, embryological analyses of Dlx gene expression and function remain limited in terms of species diversity. By supplementing our original analysis of shark embryos with previous data from other osteichthyans, such as mice and zebrafish, we support the previous speculation that the nested Dlx expression in the pharyngeal arch is likely a shared feature among all the extant jawed vertebrates. Here, we highlight several hitherto unaddressed issues regarding the evolution and function of Dlx genes, with special reference to the craniofacial development of vertebrates.

Sequencing and analysis of full-length cDNAs, 5'-ESTs and 3'-ESTs from a cartilaginous fish, the elephant shark (Callorhinchus milii)

Cartilaginous fishes are the most ancient group of living jawed vertebrates (gnathostomes) and are, therefore, an important reference group for understanding the evolution of vertebrates. The elephant shark (Callorhinchus milii), a holocephalan cartilaginous fish, has been identified as a model cartilaginous fish genome because of its compact genome (∼910 Mb) and a genome project has been initiated to obtain its whole genome sequence. In this study, we have generated and sequenced full-length enriched cDNA libraries of the elephant shark using the 'oligo-capping' method and Sanger sequencing. A total of 6,778 full-length protein-coding cDNA and 10,701 full-length noncoding cDNA were sequenced from six tissues (gills, intestine, kidney, liver, spleen, and testis) of the elephant shark. Analysis of their polyadenylation signals showed that polyadenylation usage in elephant shark is similar to that in mammals. Furthermore, both coding and noncoding transcripts of the elephant shark use the same proportion of canonical polyadenylation sites. Besides BLASTX searches, protein-coding transcripts were annotated by Gene Ontology, InterPro domain, and KEGG pathway analyses. By comparing elephant shark genes to bony vertebrate genes, we identified several ancient genes present in elephant shark but differentially lost in tetrapods or teleosts. Only ∼6% of elephant shark noncoding cDNA showed similarity to known noncoding RNAs (ncRNAs). The rest are either highly divergent ncRNAs or novel ncRNAs. In addition to full-length transcripts, 30,375 5'-ESTs and 41,317 3'-ESTs were sequenced and annotated. The clones and transcripts generated in this study are valuable resources for annotating transcription start sites, exon-intron boundaries, and UTRs of genes in the elephant shark genome, and for the functional characterization of protein sequences. These resources will also be useful for annotating genes in other cartilaginous fishes whose genomes have been targeted for whole genome sequencing.

Transcriptional enhancers in protein-coding exons of vertebrate developmental genes

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA.

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

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

A natural deletion of the HoxC cluster in elasmobranch fishes

Hox proteins are a metazoan-specific family of transcription factors that are required for developmental patterning. The genomic arrangement of Hox genes into four paralogous clusters is a primitive feature of jawed vertebrates. By using high-throughput sequencing, we demonstrate the absence of all HoxC transcripts from embryos of the shark Scyliorhinus canicula and the skate Leucoraja erinacea and the absence of all HoxC genes and two HoxC-associated microRNAs from the genome of L. erinacea. These data suggest a loss of the entire HoxC cluster in elasmobranch fishes and represent evidence for the natural deletion of an entire Hox cluster in vertebrates.

First analysis of multiple paternity in an oviparous shark, the small-spotted catshark (Scyliorhinus canicula L.)

Multiple paternity (MP) has been demonstrated in a variety of sharks, although its prevalence and the number of sires per litter vary considerably among species. To date, such analyses have focused on viviparous species that possess only part of the wide spectrum of reproductive strategies developed in elasmobranchs. We analyzed MP in an oviparous species, the small-spotted catshark (Scyliorhinus canicula). In total, 150 neonates originating from 13 different mothers were genotyped using 12 microsatellite loci. MP was commonplace, with progeny from 92% of females sired by multiple males. This result is consistent with the reproductive biology of the species, particularly its protracted breeding season and potential for long-term sperm storage. The significance of these findings is discussed in light of small-spotted catshark behavior, which suggests that the cost of avoiding mating attempts initiated by males may be high and is therefore supportive of convenience polyandry as an explanation for MP. Eggs were followed from the time they were laid to when they hatched, offering a rare opportunity to investigate juvenile development in more detail.

Revisiting the origin of the vertebrate Hox14 by including its relict sarcopterygian members

Bilaterian Hox genes play pivotal roles in the specification of positional identities along the anteroposterior axis. Particularly in vertebrates, their regulation is tightly coordinated by tandem arrays of genes [paralogy groups (PGs)] in four gene clusters (HoxA-D). Traditionally, the uninterrupted Hox cluster (Hox1-14) of the invertebrate chordate amphioxus was regarded as an archetype of the vertebrate Hox clusters. In contrast to Hox1-13 that are globally regulated by the "Hox code" and are often phylogenetically conserved, vertebrate Hox14 members were only recently revealed to be present in an African lungfish, a coelacanth, chondrichthyans and a lamprey, and decoupled from the Hox code. In this study we performed a PCR-based search of Hox14 members from diverse vertebrates, and identified one in the Australian lungfish, Neoceratodus forsteri. Based on a molecular phylogenetic analysis, this gene was designated NfHoxA14. Our real-time RT-PCR suggested its hindgut-associated expression, previously observed also in cloudy catshark HoxD14 and lamprey Hox14α. It is likely that this altered expression scheme was established before the Hox cluster quadruplication, probably at the base of extant vertebrates. To investigate the origin of vertebrate Hox14, by including this sarcopterygian Hox14 member, we performed focused phylogenetic analyses on its relationship with other vertebrate posterior Hox PGs (Hox9-13) as well as amphioxus posterior Hox genes. Our results confirmed the hypotheses previously proposed by other studies that vertebrate Hox14 does not have any amphioxus ortholog, and that none of 1-to-1 pairs of vertebrate and amphioxus posterior Hox genes, based on their relative location in the clusters, is orthologous.

Overview of the transcriptome profiles identified in hagfish, shark, and bichir: current issues arising from some nonmodel vertebrate taxa

Because of their crucial phylogenetic positions, hagfishes, sharks, and bichirs are recognized as key taxa in our understanding of vertebrate evolution. The expression patterns of the regulatory genes involved in developmental patterning have been analyzed in the context of evolutionary developmental studies. However, in a survey of public sequence databases, we found that the large-scale sequence data for these taxa are still limited. To address this deficit, we used conventional Sanger DNA sequencing and a next-generation sequencing technology based on 454 GS FLX sequencing to obtain expressed sequence tags (ESTs) of the Japanese inshore hagfish (Eptatretus burgeri; 161,482 ESTs), cloudy catshark (Scyliorhinus torazame; 165,819 ESTs), and gray bichir (Polypterus senegalus; 34,336 ESTs). We deposited the ESTs in a newly constructed database, designated the "Vertebrate TimeCapsule." The ESTs include sequences from genes that can be effectively used in evolutionary developmental studies; for instance, several encode cartilaginous extracellular matrix proteins, which are central to an understanding of the ways in which evolutionary processes affected the skeletal elements, whereas others encode regulatory genes involved in craniofacial development and early embryogenesis. Here, we discuss how hagfishes, sharks, and bichirs contribute to our understanding of vertebrate evolution, we review the current status of the publicly available sequence data for these three taxa, and we introduce our EST projects and newly developed database.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

Regionalization of the shark hindbrain: a survey of an ancestral organization

Cartilaginous fishes (chondrichthyans) represent an ancient radiation of vertebrates currently considered the sister group of the group of gnathostomes with a bony skeleton that gave rise to land vertebrates. This out-group position makes chondrichthyans essential in assessing the ancestral organization of the brain of jawed vertebrates. To gain knowledge about hindbrain evolution we have studied its development in a shark, the lesser spotted dogfish Scyliorhinus canicula by analyzing the expression of some developmental genes and the origin and distribution of specific neuronal populations, which may help to identify hindbrain subdivisions and boundaries and the topology of specific cell groups. We have characterized three developmental periods that will serve as a framework to compare the development of different neuronal systems and may represent a suitable tool for comparing the absolute chronology of development among vertebrates. The expression patterns of Pax6, Wnt8, and HoxA2 genes in early embryos of S. canicula showed close correspondence to what has been described in other vertebrates and helped to identify the anterior rhombomeres. Also in these early embryos, the combination of Pax6 with protein markers of migrating neuroblasts (DCX) and early differentiating neurons (general: HuC/D; neuron type specific: GAD, the GABA synthesizing enzyme) revealed the organization of S. canicula hindbrain in both transverse segmental units corresponding to visible rhombomeres and longitudinal columns. Later in development, when the interrhombomeric boundaries fade away, accurate information about S. canicula hindbrain subdivisions was achieved by comparing the expression patterns of Pax6 and GAD, serotonin (serotoninergic neurons), tyrosine hydroxylase (catecholaminergic neurons), choline acetyltransferase (cholinergic neurons), and calretinin (a calcium-binding protein). The patterns observed revealed many topological correspondences with other vertebrates and led to reconsideration of the current view of the elasmobranch hindbrain segmentation as peculiar among vertebrates.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.