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

Contingency and chance erase necessity in the experimental evolution of ancestral proteins

The roles of chance, contingency, and necessity in evolution are unresolved because they have never been assessed in a single system or on timescales relevant to historical evolution. We combined ancestral protein reconstruction and a new continuous evolution technology to mutate and select proteins in the B-cell lymphoma-2 (BCL-2) family to acquire protein–protein interaction specificities that occurred during animal evolution. By replicating evolutionary trajectories from multiple ancestral proteins, we found that contingency generated over long historical timescales steadily erased necessity and overwhelmed chance as the primary cause of acquired sequence variation; trajectories launched from phylogenetically distant proteins yielded virtually no common mutations, even under strong and identical selection pressures. Chance arose because many sets of mutations could alter specificity at any timepoint; contingency arose because historical substitutions changed these sets. Our results suggest that patterns of variation in BCL-2 sequences – and likely other proteins, too – are idiosyncratic products of a particular and unpredictable course of historical events.

One of the most fundamental and unresolved questions in evolutionary biology is whether the outcomes of evolution are predictable. Is the diversity of life we see today the expected result of organisms adapting to their environment throughout history (also known as natural selection) or the product of random chance? Or did chance events early in history shape the paths that evolution could take next, determining the biological forms that emerged under natural selection much later?

These questions are hard to study because evolution happened only once, long ago. To overcome this barrier, Xie, Pu, Metzger et al. developed an experimental approach that can evolve reconstructed ancestral proteins that existed deep in the past. Using this method, it is possible to replay evolution multiple times, from various historical starting points, under conditions similar to those that existed long ago. The end products of the evolutionary trajectories can then be compared to determine how predictable evolution actually is.

Xie, Pu, Metzger et al. studied proteins belonging to the BCL-2 family, which originated some 800 million years ago. These proteins have diversified greatly over time in both their genetic sequences and their ability to bind to specific partner proteins called co-regulators. Xie, Pu, Metzger et al. synthesized BCL-2 proteins that existed at various times in the past. Each ancestral protein was then allowed to evolve repeatedly under natural selection to acquire the same co-regulator binding functions that evolved during history.

At the end of each evolutionary trajectory, the genetic sequence of the resulting BCL-2 proteins was recorded. This revealed that the outcomes of evolution were almost completely unpredictable: trajectories initiated from the same ancestral protein produced proteins with very different sequences, and proteins launched from different ancestral starting points were even more dissimilar.

Further experiments identified the mutations in each trajectory that caused changes in coregulator binding. When these mutations were introduced into other ancestral proteins, they did not yield the same change in function. This suggests that early chance events influenced each protein’s evolution in an unpredictable way by opening and closing the paths available to it in the future.

This research expands our understanding of evolution on a molecular level whilst providing a new experimental approach for studying evolutionary drivers in more detail. The results suggest that BCL-2 proteins, in all their various forms, are unique products of a particular, unpredictable course of history set in motion by ancient chance events.

Development and regeneration of the crushing dentition in skates (Rajidae)

Sharks and rays (elasmobranchs) have the remarkable capacity to continuously regenerate their teeth. The polyphyodont system is considered the ancestral condition of the gnathostome dentition. Despite this shared regenerative ability, sharks and rays exhibit dramatic interspecific variation in their tooth morphology. Ray (batoidea) teeth typically constitute crushing pads of flattened teeth, whereas shark teeth are pointed, multi-cuspid units. Although recent research has addressed the molecular development of the shark dentition, little is known about that of the ray. Furthermore, how dental diversity within the elasmobranch lineage is achieved remains unknown. Here, we examine dental development and regeneration in two Batoid species: the thornback skate (Raja clavata) and the little skate (Leucoraja erinacea). Using in situ hybridization and immunohistochemistry, we examine the expression of a core gnathostome dental gene set during early development of the skate dentition and compare it to development in the shark. Elasmobranch tooth development is highly conserved, with sox2 likely playing an important role in the initiation and regeneration of teeth. Alterations to conserved genes expressed in an enamel knot-like signalling centre may explain the morphological diversity of elasmobranch teeth, thereby enabling sharks and rays to occupy diverse dietary and ecological niches.

Analysis of pallial/cortical interneurons in key vertebrate models of Testudines, Anurans and Polypteriform fishes

The organization of the pallial derivatives across vertebrates follows a comparable elementary arrangement, although not all of them possess a layered cortical structure as sophisticated as the cerebral cortex of mammals. However, its expansion along evolution has only been possible by the development and coevolution of the cellular networks formed by excitatory neurons and inhibitory interneurons. Thus, the comparative analysis of interneuron types in vertebrate models of key evolutionary significance will provide important information, due to the extraordinary anatomical sophistication of their interneuron systems with simpler behavioral implications. Particularly in mammals, the main consensus for classifying interneuron types is based on non-overlapping markers, which do not form a single population, but consist of several distinct classes of inhibitory cells showing co-expression of other markers. In our study, we analyzed immunohistochemically the expression of the main markers like somatostatin (SOM), parvalbumin (PV), calretinin (CR), calbindin (CB), neuropeptide Y (NPY) and/or nitric oxide synthase (NOS) at the pallial regions of three different models of Osteichthyes. First, we selected two tetrapods, one amniote from the genus Pseudemys belonging to the order Testudine, at the base of the amniote diversification and with a three-layered simple cortex, and the Anuran Xenopus laevis, an anamniote tetrapod with a non-layered evaginated pallium, and finally the order Polypteriform, a small fish group at the base of the actinopterygian diversification with an everted telencephalon. SOM was the most conserved interneuron type in terms of its distribution and co-expression with other markers such as CR, in contrast to PV, which showed a different pattern between the models analyzed. In addition, the SOM expression supports a homological relationship between the medial pallial derivatives in all the models. CR and CB expressions in the tetrapods were observed, particularly, CR expressing cells were detected in the medial and the dorsal pallial derivatives, in contrast to CB, which appeared only in discrete scattered populations. However, the pallium of Polypteriforms fishes was almost devoid of CR cells, in contrast to the important number of CB cells observed in all the pallial regions. The NPY immunoreactivity was detected in all the pallial domains of all the models, as well as cells coexpressing CR. Finally, the pallial nitrergic expression was also conserved, which allows to postulate the homological relationships between the ventropallial and the amygdaloid derivatives. In summary, even in basal pallial models the neurochemically characterized interneurons indicate that their first appearance took place before the common ancestor of amniotes. Thus, our results suggest a shared pattern of interneuron types in the pallium of all Osteichthyes.

Reflections on the evolution of the vertebrate tissue inhibitors of metalloproteinases

Jawed vertebrates (Gnathostomes) have 4 tissue inhibitors of metalloproteinases (TIMPs), multifunctional proteins that all inhibit members of the large matrix metalloproteinase (MMP) family but differ in their other roles, including the regulation of pro-MMP activation, cell growth, apoptosis and angiogenesis, and the structure of extracellular matrices (ECMs). Molecular phylogeny analyses indicate that vertebrate TIMP genes arose from an invertebrate ancestor through 3 successive duplications, possibly including 2 whole genome duplications, during early vertebrate phylogeny. TIMPs from invertebrates also inhibit metalloproteinases, bind to pro-MMPs, and contribute to ECM structures but are not orthologs of any particular vertebrate TIMP. The most ancient vertebrate superclass, the Agnatha (jawless fish), seems to provide a snapshot of a stage in TIMP evolution preceding the third gene duplication. This review examines the structures of TIMPs from different vertebrate orders using information relating to the structural basis of their various functions. Provisional conclusions are that during their evolutionary divergence, various TIMPs lost inhibitory activity toward some metalloproteinases, specialized in effects on different pro-MMPs, and developed new interactions with discrete targets (including integrins and receptors), while recapitulating a role in ECM structure. The analysis is limited by the sparse information available regarding the functional properties of nonmammalian TIMPs.-Brew, K. Reflections on the evolution of the vertebrate tissue inhibitors of metalloproteinases.

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.

Outgroup, alignment and modelling improvements indicate that two TNFSF13-like genes existed in the vertebrate ancestor

The molecular machinery required for lymphocyte development and differentiation appears to have emerged concomitantly with distinct B- and T-like lymphocyte subsets in the ancestor of all vertebrates. The TNFSF superfamily (TNFSF) members BAFF (TNFSF13/Blys) and APRIL (TNFSF13) are key regulators of B cell development survival, and activation in mammals, but the temporal emergence of these molecules, and their precise relationship to the newly identified TNFSF gene BALM (BAFF and APRIL-like molecule), have not yet been elucidated. Here, to resolve the early evolutionary history of this family, we improved outgroup sampling and alignment quality, and applied better fitting substitution models compared to past studies. Our analyses reveal that BALM is a definitive TNFSF13 family member, which split from BAFF in the gnathostome (jawed vertebrate) ancestor. Most importantly, however, we show that both the APRIL and BAFF lineages existed in the ancestors of all extant vertebrates. This implies that APRIL has been lost, or is yet to be found, in cyclostomes (jawless vertebrates). Our results suggest that lineage-specific gene duplication and loss events have caused lymphocyte regulation, despite shared origins, to become secondarily distinct between gnathostomes and cyclostomes. Finally, the structure of lamprey BAFF-like, and its phylogenetic placement as sister to BAFF and BALM, but not the more slowly evolving APRIL, indicates that the primordial lymphocyte regulator was more APRIL-like than BAFF-like.

Evidence from cyclostomes for complex regionalization of the ancestral vertebrate brain

The vertebrate brain is highly complex, but its evolutionary origin remains elusive. Because of the absence of certain developmental domains generally marked by the expression of regulatory genes, the embryonic brain of the lamprey, a jawless vertebrate, had been regarded as representing a less complex, ancestral state of the vertebrate brain. Specifically, the absence of a Hedgehog- and Nkx2.1-positive domain in the lamprey subpallium was thought to be similar to mouse mutants in which the suppression of Nkx2-1 leads to a loss of the medial ganglionic eminence. Here we show that the brain of the inshore hagfish (Eptatretus burgeri), another cyclostome group, develops domains equivalent to the medial ganglionic eminence and rhombic lip, resembling the gnathostome brain. Moreover, further investigation of lamprey larvae revealed that these domains are also present, ruling out the possibility of convergent evolution between hagfish and gnathostomes. Thus, brain regionalization as seen in crown gnathostomes is not an evolutionary innovation of this group, but dates back to the latest vertebrate ancestor before the divergence of cyclostomes and gnathostomes more than 500 million years ago.

An ancient dental gene set governs development and continuous regeneration of teeth in sharks

The evolution of oral teeth is considered a major contributor to the overall success of jawed vertebrates. This is especially apparent in cartilaginous fishes including sharks and rays, which develop elaborate arrays of highly specialized teeth, organized in rows and retain the capacity for life-long regeneration. Perpetual regeneration of oral teeth has been either lost or highly reduced in many other lineages including important developmental model species, so cartilaginous fishes are uniquely suited for deep comparative analyses of tooth development and regeneration. Additionally, sharks and rays can offer crucial insights into the characters of the dentition in the ancestor of all jawed vertebrates. Despite this, tooth development and regeneration in chondrichthyans is poorly understood and remains virtually uncharacterized from a developmental genetic standpoint. Using the emerging chondrichthyan model, the catshark (Scyliorhinus spp.), we characterized the expression of genes homologous to those known to be expressed during stages of early dental competence, tooth initiation, morphogenesis, and regeneration in bony vertebrates. We have found that expression patterns of several genes from Hh, Wnt/β-catenin, Bmp and Fgf signalling pathways indicate deep conservation over ~450 million years of tooth development and regeneration. We describe how these genes participate in the initial emergence of the shark dentition and how they are redeployed during regeneration of successive tooth generations. We suggest that at the dawn of the vertebrate lineage, teeth (i) were most likely continuously regenerative structures, and (ii) utilised a core set of genes from members of key developmental signalling pathways that were instrumental in creating a dental legacy redeployed throughout vertebrate evolution. These data lay the foundation for further experimental investigations utilizing the unique regenerative capacity of chondrichthyan models to answer evolutionary, developmental, and regenerative biological questions that are impossible to explore in classical models.

Organization of the nitrergic neuronal system in the primitive bony fishes Polypterus senegalus and Erpetoichthys calabaricus (Actinopterygii: Cladistia)

Cladistians are a group of basal actinopterygian fishes that constitute a good model for studying primitive brain features, most likely present in the ancestral bony fishes. The analysis of the nitrergic neurons (with the enzyme nitric oxide synthase; NOS) has helped in understanding important aspects of brain organization in all vertebrates studied. We investigated the nitrergic system of two cladistian species by means of specific antibodies against NOS and NADPH-diaphorase (NADPH-d) histochemistry, which, with the exception of the primary olfactory and terminal nerve fibers, labeled only for NADPH-d, yielded identical results. Double immunohistochemistry was conducted for simultaneous detection of NOS with tyrosine hydroxylase, choline acetyltransferase, calbindin, calretinin, and serotonin, to establish accurately the localization of the nitrergic neurons and fibers and to assess possible interactions between these neuroactive substances. The pattern of distribution in both species showed only subtle differences in the density of labeled cells. Distinct groups of NOS-immunoreactive cells were observed in pallial and subpallial areas, paraventricular region, tuberal and retromammillary hypothalamic areas, posterior tubercle, prethalamic and thalamic areas, optic tectum, torus semicircularis, mesencephalic tegmentum, interpeduncular nucleus, superior and middle reticular nuclei, magnocellular vestibular nucleus, solitary tract nucleus, nucleus medianus magnocellularis, the spinal cord and amacrine cells in the retina. Large neurons in cranial nerve sensory ganglia were also labeled. The comparison of these results with those from other vertebrates, using a neuromeric analysis, reveals a conserved pattern of organization of the nitrergic system from this primitive fish group to amniotes, including mammals.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0022-6) contains supplementary material, which is available to authorized users.

The evolution of tenascins and fibronectin

Tenascins are extracellular matrix glycoproteins that act both as integrin ligands and as modifiers of fibronectin-integrin interactions to regulate cell adhesion, migration, proliferation and differentiation. In tetrapods, both tenascins and fibronectin bind to integrins via RGD and LDV-type tripeptide motifs found in exposed loops in their fibronectin-type III domains. We previously showed that tenascins appeared early in the chordate lineage and are represented by single genes in extant cephalochordates and tunicates. Here we have examined the genomes of the coelacanth Latimeria chalumnae, the elephant shark Callorhinchus milii as well as the lampreys Petromyzon marinus and Lethenteron japonicum to learn more about the evolution of the tenascin gene family as well as the timing of the appearance of fibronectin during chordate evolution. The coelacanth has 4 tenascins that are more similar to tetrapod tenascins than are tenascins from ray-finned fishes. In contrast, only 2 tenascins were identified in the elephant shark and the Japanese lamprey L. japonicum. An RGD motif exposed to integrin binding is observed in tenascins from many, but not all, classes of chordates. Tetrapods that lack this RGD motif in tenascin-C have a similar motif in the paralog tenascin-W, suggesting the potential for some overlapping function. A predicted fibronectin with the same domain organization as the fibronectin from tetrapods is found in the sea lamprey P. marinus but not in tunicates, leading us to infer that fibronectin first appeared in vertebrates. The motifs that recognize LDV-type integrin receptors are conserved in fibronectins from a broad spectrum of vertebrates, but the RGD integrin-binding motif may have evolved in gnathostomes.

Evolution of endothelin receptors in vertebrates

Endothelin receptors are G protein coupled receptors (GPCRs) of the β-group of rhodopsin receptors that bind to endothelin ligands, which are 21 amino acid long peptides derived from longer prepro-endothelin precursors. The most basal Ednr-like GPCR is found outside vertebrates in the cephalochordate amphioxus, but endothelin ligands are only present among vertebrates, including the lineages of jawless vertebrates (lampreys and hagfishes), cartilaginous vertebrates (sharks, rays, and chimaeras), and bony vertebrates (ray-finned fishes and lobe-finned vertebrates including tetrapods). A bona fide endothelin system is thus a vertebrate-specific innovation with important roles for regulating the cardiovascular system, renal and pulmonary processes, as well as for the development of the vertebrate-specific neural crest cell population and its derivatives. Expectedly, dysregulation of endothelin receptors and the endothelin system leads to a multitude of human diseases. Despite the importance of different types of endothelin receptors for vertebrate development and physiology, current knowledge on endothelin ligand-receptor interactions, on the expression of endothelin receptors and their ligands, and on the functional roles of the endothelin system for embryonic development and in adult vertebrates is very much biased towards amniote vertebrates. Recent analyses from a variety of vertebrate lineages, however, have shown that the endothelin system in lineages such as teleost fish and lampreys is more diverse and is divergent from the mammalian endothelin system. This diversity is mainly based on differential evolution of numerous endothelin system components among vertebrate lineages generated by two rounds of whole genome duplication (three in teleosts) during vertebrate evolution. Here we review current understanding of the evolutionary history of the endothelin receptor family in vertebrates supplemented with surveys on the endothelin receptor gene complement of newly available genome assemblies from phylogenetically informative taxa. Our assessment further highlights the diversity of the vertebrate endothelin system and calls for detailed functional and pharmacological analyses of the endothelin system beyond tetrapods.

Organization of the orexin/hypocretin system in the brain of two basal actinopterygian fishes, the cladistians Polypterus senegalus and Erpetoichthys calabaricus

Cladistians are primitive actinopterygian fishes mostly neglected in neuroanatomical studies. In the present study, the detailed neuroanatomical distribution of orexin (hypocretin)-like immunoreactive (OX-ir) cell bodies and fibers was analyzed in the brain of two species representative of the two extant genera of cladistians. Antibodies against mammalian orexin-A and orexin-B peptides were used. Simultaneous detection of orexins with neuropeptide Y (NPY), tyrosine hydroxylase (TH), and serotonin (5-HT) was used to establish accurately the topography of the orexin system and to evaluate the possible interactions with NPY and monoaminergic systems. A largely common pattern of OX-ir distribution in the two cladistian species was observed. Most OX-ir cells were located in the suprachiasmatic nucleus and tuberal hypothalamus, whereas scarce cells were observed in the posterior tubercle. In addition, a population of OX-ir cells was found in the preoptic area only in Polypterus and some cells also contained TH. The observed widespread distribution of OX-ir fibers was especially abundant in the retrobulbar area, subpallial areas, preoptic area, suprachiasmatic nucleus, tuberal hypothalamic area, prethalamus, thalamus, pretectum, optic tectum, and tegmentum. Low innervation was found in relation to monoaminergic cell groups, whereas a high NPY innervation was observed in all OX-ir cell groups. These relationships would represent the anatomical substrate for the functional interdependence between these systems. The organization of the orexin system in cladistians revealed a pattern largely consistent with those reported for all studied groups of vertebrates, suggesting that the primitive organization of this peptidergic system occurred in the common ancestor of gnathostome vertebrates.

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.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

Transcriptome profiling reveals higher vertebrate orthologous of intra-cytoplasmic pattern recognition receptors in grey bamboo shark

From an immunologist perspective, sharks are an important group of jawed cartilaginous fishes and survey of the public database revealed a great gap in availability of large-scale sequence data for the group of Chondrichthyans the elasmobranchs. In an attempt to bridge this deficit we generated the transcriptome from the spleen and kidney tissues (a total of 1,606,172 transcripts) of the shark, Chiloscyllium griseum using the Illumina HiSeq2000 platform. With a cut off of > = 300 bp and an expression value of >1RPKM we used 43,385 transcripts for BLASTX analysis which revealed 17,548 transcripts matching to the NCBI nr database with an E-value of < = 10(-5) and similarity score of 40%. The longest transcript was 16,974 bases with matched to HECT domain containing E3 ubiqutin protein ligase. MEGAN4 annotation pipeline revealed immune and signalling pathways including cell adhesion molecules, cytokine-cytokine receptor interaction, T-cell receptor signalling pathway and chemokine signaling pathway to be highly expressed in spleen, while different metabolism pathways such as amino acid metabolism, carbohydrate metabolism, lipid metabolism and xenobiotic biodegradation were highly expressed in kidney. Few of the candidate genes were selected to analyze their expression levels in various tissues by real-time PCR and also localization of a receptor by in-situ PCR to validate the prediction. We also predicted the domains structures of some of the identified pattern recognition receptors, their phylogenetic relationship with lower and higher vertebrates and the complete downstream signaling mediators of classical dsRNA signaling pathway. The generated transcriptome will be a valuable resource to further genetic and genomic research in elasmobranchs.

Definition of a third VLR gene in hagfish

Jawless vertebrates (cyclostomes) have an alternative adaptive immune system in which lymphocytes somatically diversify their variable lymphocyte receptors (VLR) through recombinatorial use of leucine-rich repeat cassettes during VLR gene assembly. Three types of these anticipatory receptors in lampreys (VLRA, VLRB, and VLRC) are expressed by separate lymphocyte lineages. However, only two VLR genes (VLRA and VLRB) have been found in hagfish. Here we have identified a third hagfish VLR, which undergoes somatic assembly to generate sufficient diversity to encode a large repertoire of anticipatory receptors. Sequence analysis, structural comparison, and phylogenetic analysis indicate that the unique hagfish VLR is the counterpart of lamprey VLRA and the previously identified hagfish "VLRA" is the lamprey VLRC counterpart. The demonstration of three orthologous VLR genes in both lampreys and hagfish suggests that this anticipatory receptor system evolved in a common ancestor of the two cyclostome lineages around 480 Mya.

Development of the head and trunk mesoderm in the dogfish, Scyliorhinus torazame: II. Comparison of gene expression between the head mesoderm and somites with reference to the origin of the vertebrate head

The vertebrate mesoderm differs distinctly between the head and trunk, and the evolutionary origin of the head mesoderm remains enigmatic. Although the presence of somite-like segmentation in the head mesoderm of model animals is generally denied at molecular developmental levels, the appearance of head cavities in elasmobranch embryos has not been explained, and the possibility that they may represent vestigial head somites once present in an amphioxus-like ancestor has not been ruled out entirely. To examine whether the head cavities in the shark embryo exhibit any molecular signatures reminiscent of trunk somites, we isolated several developmentally key genes, including Pax1, Pax3, Pax7, Pax9, Myf5, Sonic hedgehog, and Patched2, which are involved in myogenic and chondrogenic differentiation in somites, and Pitx2, Tbx1, and Engrailed2, which are related to the patterning of the head mesoderm, from an elasmobranch species, Scyliorhinus torazame. Observation of the expression patterns of these genes revealed that most were expressed in patterns that resembled those found in amniote embryos. In addition, the head cavities did not exhibit an overt similarity to somites; that is, the similarity was no greater than that of the unsegmented head mesoderm in other vertebrates. Moreover, the shark head mesoderm showed an amniote-like somatic/visceral distinction according to the expression of Pitx2, Tbx1, and Engrailed2. We conclude that the head cavities do not represent a manifestation of ancestral head somites; rather, they are more likely to represent a derived trait obtained in the lineage of gnathostomes.

Late development of hagfish vertebral elements

It has been demonstrated recently that hagfishes, one of two groups of extant jawless vertebrates, have cartilaginous vertebral elements. Embryological and gene expression analyses have also shown that this group of animals develops a sclerotome, the potential primordium of the axial skeleton. However, it has not been shown unequivocally that the hagfish sclerotome truly differentiates into cartilage, because access to late-stage embryos and information about the cartilaginous extracellular matrix (ECM) are lacking for these animals. Here we investigated the expression patterns of the biglycan/decorin (BGN/DCN) gene in the inshore hagfish, Eptatretus burgeri. The homologue of this gene encodes the major noncollagenous component of the cartilaginous ECM among gnathostomes. We clearly identified the expression of this gene in adult vertebral tissues and in embryonic mesenchymal cells on the ventral aspect of the notochord. Taking into account that the sclerotome in the gnathostomes expresses BGN/DCN gene during the chondrogenesis, it is highly expected the hagfish BGN/DCN-positive mesenchymal cells are derived from the sclerotomes. We propose that hagfishes and gnathostomes share conserved developmental mechanisms not only in their somite differentiation, but also in chondrogenesis of their vertebral elements.

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.

Saltatory evolution of the ectodermal neural cortex gene family at the vertebrate origin

The ectodermal neural cortex (ENC) gene family, whose members are implicated in neurogenesis, is part of the kelch repeat superfamily. To date, ENC genes have been identified only in osteichthyans, although other kelch repeat-containing genes are prevalent throughout bilaterians. The lack of elaborate molecular phylogenetic analysis with exhaustive taxon sampling has obscured the possible link of the establishment of this gene family with vertebrate novelties. In this study, we identified ENC homologs in diverse vertebrates by means of database mining and polymerase chain reaction screens. Our analysis revealed that the ENC3 ortholog was lost in the basal eutherian lineage through single-gene deletion and that the triplication between ENC1, -2, and -3 occurred early in vertebrate evolution. Including our original data on the catshark and the zebrafish, our comparison revealed high conservation of the pleiotropic expression pattern of ENC1 and shuffling of expression domains between ENC1, -2, and -3. Compared with many other gene families including developmental key regulators, the ENC gene family is unique in that conventional molecular phylogenetic inference could identify no obvious invertebrate ortholog. This suggests a composite nature of the vertebrate-specific gene repertoire, consisting not only of de novo genes introduced at the vertebrate origin but also of long-standing genes with no apparent invertebrate orthologs. Some of the latter, including the ENC gene family, may be too rapidly evolving to provide sufficient phylogenetic signals marking orthology to their invertebrate counterparts. Such gene families that experienced saltatory evolution likely remain to be explored and might also have contributed to phenotypic evolution 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.

Comparative transcriptome analysis between planarian Dugesia japonica and other platyhelminth species

Background

Planarians are considered to be among the extant animals close to one of the earliest groups of organisms that acquired a central nervous system (CNS) during evolution. Planarians have a bilobed brain with nine lateral branches from which a variety of external signals are projected into different portions of the main lobes. Various interneurons process different signals to regulate behavior and learning/memory. Furthermore, planarians have robust regenerative ability and are attracting attention as a new model organism for the study of regeneration. Here we conducted large-scale EST analysis of the head region of the planarian Dugesia japonica to construct a database of the head-region transcriptome, and then performed comparative analyses among related species.

Results

A total of 54,752 high-quality EST reads were obtained from a head library of the planarian Dugesia japonica, and 13,167 unigene sequences were produced by de novo assembly. A new method devised here revealed that proteins related to metabolism and defense mechanisms have high flexibility of amino-acid substitutions within the planarian family. Eight-two CNS-development genes were found in the planarian (cf. C. elegans 3; chicken 129). Comparative analysis revealed that 91% of the planarian CNS-development genes could be mapped onto the schistosome genome, but one-third of these shared genes were not expressed in the schistosome.

Conclusions

We constructed a database that is a useful resource for comparative planarian transcriptome studies. Analysis comparing homologous genes between two planarian species showed that the potential of genes is important for accumulation of amino-acid substitutions. The presence of many CNS-development genes in our database supports the notion that the planarian has a fundamental brain with regard to evolution and development at not only the morphological/functional, but also the genomic, level. In addition, our results indicate that the planarian CNS-development genes already existed before the divergence of planarians and schistosomes from their common ancestor.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.