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

Developmental constraints on fin diversity

The fish fin is a breathtaking repository full of evolutionary diversity, novelty, and convergence. Over 500 million years, the adaptation to novel habitats has provided landscapes of fin diversity. Although comparative anatomy of evolutionarily divergent patterns over centuries has highlighted the fundamental architectures and evolutionary trends of fins, including convergent evolution, the developmental constraints on fin evolution, which bias the evolutionary trajectories of fin morphology, largely remain elusive. Here, we review the evolutionary history, developmental mechanisms, and evolutionary underpinnings of paired fins, illuminating possible developmental constraints on fin evolution. Our compilation of anatomical and genetic knowledge of fin development sheds light on the canalized and the unpredictable aspects of fin shape in evolution. Leveraged by an arsenal of genomic and genetic tools within the working arena of spectacular fin diversity, evolutionary developmental biology embarks on the establishment of conceptual framework for developmental constraints, previously enigmatic properties of evolution.

Vision in sharks and rays: Opsin diversity and colour vision

The visual sense of elasmobranch fishes is poorly studied compared to their bony cousins, the teleosts. Nevertheless, the elasmobranch eye features numerous specialisations that have no doubt facilitated the diversification and evolutionary success of this fascinating taxon. In this review, I highlight recent discoveries on the nature and phylogenetic distribution of visual pigments in sharks and rays. Whereas most rays appear to be cone dichromats, all sharks studied to date are cone monochromats and, as a group, have likely abandoned colour vision on multiple occasions. This situation in sharks mirrors that seen in other large marine predators, the pinnipeds and cetaceans, which leads us to reassess the costs and benefits of multiple cone pigments and wavelength discrimination in the marine environment.

Molecular and cellular evolution of corticogenesis in amniotes

The cerebral cortex varies dramatically in size and complexity between amniotes due to differences in neuron number and composition. These differences emerge during embryonic development as a result of variations in neurogenesis, which are thought to recapitulate modifications occurred during evolution that culminated in the human neocortex. Here, we review work from the last few decades leading to our current understanding of the evolution of neurogenesis and size of the cerebral cortex. Focused on specific examples across vertebrate and amniote phylogeny, we discuss developmental mechanisms regulating the emergence, lineage, complexification and fate of cortical germinal layers and progenitor cell types. At the cellular level, we discuss the fundamental impact of basal progenitor cells and the advent of indirect neurogenesis on the increased number and diversity of cortical neurons and layers in mammals, and on cortex folding. Finally, we discuss recent work that unveils genetic and molecular mechanisms underlying this progressive expansion and increased complexity of the amniote cerebral cortex during evolution, with a particular focus on those leading to human-specific features. Whereas new genes important in human brain development emerged the recent hominid lineage, regulation of the patterns and levels of activity of highly conserved signaling pathways are beginning to emerge as mechanisms of central importance in the evolutionary increase in cortical size and complexity across amniotes.

Transcription factor AP2 controls cnidarian germ cell induction

Clonal animals do not sequester a germ line during embryogenesis. Instead, they have adult stem cells that contribute to somatic tissues or gametes. How germ fate is induced in these animals, and whether this process is related to bilaterian embryonic germline induction, is unknown. We show that transcription factor AP2 (Tfap2), a regulator of mammalian germ lines, acts to commit adult stem cells, known as i-cells, to the germ cell fate in the clonal cnidarian Hydractinia symbiolongicarpus Tfap2 mutants lacked germ cells and gonads. Transplanted wild-type cells rescued gonad development but not germ cell induction in Tfap2 mutants. Forced expression of Tfap2 in i-cells converted them to germ cells. Therefore, Tfap2 is a regulator of germ cell commitment across germ line-sequestering and germ line-nonsequestering animals.

Divergent genes encoding the putative receptors for growth hormone and prolactin in sea lamprey display distinct patterns of expression

Growth hormone receptor (GHR) and prolactin receptor (PRLR) in jawed vertebrates were thought to arise after the divergence of gnathostomes from a basal vertebrate. In this study we have identified two genes encoding putative GHR and PRLR in sea lamprey (Petromyzon marinus) and Arctic lamprey (Lethenteron camtschaticum), extant members of one of the oldest vertebrate groups, agnathans. Phylogenetic analysis revealed that lamprey GHR and PRLR cluster at the base of gnathostome GHR and PRLR clades, respectively. This indicates that distinct GHR and PRLR arose prior to the emergence of the lamprey branch of agnathans. In the sea lamprey, GHR and PRLR displayed a differential but overlapping pattern of expression; GHR had high expression in liver and heart tissues, whereas PRLR was expressed highly in the brain and moderately in osmoregulatory tissues. Branchial PRLR mRNA levels were significantly elevated by stage 5 of metamorphosis and remained elevated through stage 7, whereas levels of GHR mRNA were only elevated in the final stage (7). Branchial expression of GHR increased following seawater (SW) exposure of juveniles, but expression of PRLR was not significantly altered. The results indicate that GHR and PRLR may both participate in metamorphosis and that GHR may mediate SW acclimation.

New Insights Into the Evolutionary History of Melatonin Receptors in Vertebrates, With Particular Focus on Teleosts

In order to improve our understanding of melatonin signaling, we have reviewed and revised the evolutionary history of melatonin receptor genes (mtnr) in vertebrates. All gnathostome mtnr genes have a conserved gene organization with two exons, except for mtnr1b paralogs of some teleosts that show intron gains. Phylogeny and synteny analyses demonstrate the presence of four mtnr subtypes, MTNR1A, MTNR1B, MTNR1C, MTNR1D that arose from duplication of an ancestral mtnr during the vertebrate tetraploidizations (1R and 2R). In tetrapods, mtnr1d was lost, independently, in mammals, in archosaurs and in caecilian amphibians. All four mtnr subtypes were found in two non-teleost actinopterygian species, the spotted gar and the reedfish. As a result of teleost tetraploidization (3R), up to seven functional mtnr genes could be identified in teleosts. Conservation of the mtnr 3R-duplicated paralogs differs among the teleost lineages. Synteny analysis showed that the mtnr1d was conserved as a singleton in all teleosts resulting from an early loss after tetraploidization of one of the teleost 3R and salmonid 4R paralogs. Several teleosts including the eels and the piranha have conserved both 3R-paralogs of mtnr1a, mtnr1b, and mtnr1c. Loss of one of the 3R-paralogs depends on the lineage: mtnr1ca was lost in euteleosts whereas mtnr1cb was lost in osteoglossomorphs and several ostariophysians including the zebrafish. We investigated the tissue distribution of mtnr expression in a large range of tissues in medaka. The medaka has conserved the four vertebrate paralogs, and these are expressed in brain and retina, and, differentially, in peripheral tissues. Photoperiod affects mtnr expression levels in a gene-specific and tissue-specific manner. This study provides new insights into the repertoire diversification and functional evolution of the mtnr gene family in vertebrates.

The retinal pigments of the whale shark (Rhincodon typus) and their role in visual foraging ecology

The spectral tuning properties of the whale shark (Rhincodon typus) rod (rhodopsin or Rh1) and long-wavelength-sensitive (LWS) cone visual pigments were examined to determine whether these retinal pigments have adapted to the broadband light spectrum available for surface foraging or to the narrowband blue-shifted light spectrum available at depth. Recently published whale shark genomes have identified orthologous genes for both the whale shark Rh1 and LWS cone opsins suggesting a duplex retina. Here, the whale shark Rh1 and LWS cone opsin sequences were examined to identify amino acid residues critical for spectral tuning. Surprisingly, the predicted absorbance maximum (λmax) for both the whale shark Rh1 and LWS visual pigments is near 500 nm. Although Rh1 λmax values near 500 nm are typical of terrestrial vertebrates, as well as surface foraging fish, it is uncommon for a vertebrate LWS cone pigment to be so greatly blue-shifted. We propose that the spectral tuning properties of both the whale shark Rh1 and LWS cone pigments are most likely adaptations to the broadband light spectrum available at the surface. Whale shark melanopsin (Opn4) deactivation kinetics was examined to better understand the underlying molecular mechanisms of the pupillary light reflex. Results show that the deactivation rate of whale shark Opn4 is similar to the Opn4 deactivation rate from vertebrates possessing duplex retinae and is significantly faster than the Opn4 deactivation rate from an aquatic rod monochromat lacking functional cone photoreceptors. The rapid deactivation rate of whale shark Opn4 is consistent with a functional cone class and would provide the animal with an exponential increase in the number of photons required for photoreceptor signaling when transitioning from photopic to scotopic light conditions, as is the case when diving.

Shark tooth regeneration reveals common stem cell characters in both human rested lamina and ameloblastoma

The human dentition is a typical diphyodont mammalian system with tooth replacement of most positions. However, after dental replacement and sequential molar development, the dental lamina undergoes apoptosis and fragments, leaving scattered epithelial units (dental lamina rests; DLRs). DLRs in adult humans are considered inactive epithelia, thought to possess limited capacity for further regeneration. However, we show that these tissues contain a small proportion of proliferating cells (assessed by both Ki67 and PCNA) but also express a number of common dental stem cell markers (Sox2, Bmi1, β-catenin and PH3) similar to that observed in many vertebrates that actively, and continuously regenerate their dentition. We compared these human tissues with the dental lamina of sharks that regenerate their dentition throughout life, providing evidence that human tissues have the capacity for further and undocumented regeneration. We also assessed cases of human ameloblastoma to characterise further the proliferative signature of dental lamina rests. Ameloblastomas are assumed to derive from aberrant lamina rests that undergo changes, which are not well understood, to form a benign tumour. We suggest that dental lamina rests can offer a potential source of important dental stem cells for future dental regenerative therapy. The combined developmental genetic data from the shark dental lamina and ameloblastoma may lead to the development of novel methods to utilise these rested populations of adult lamina stem cells for controlled tooth replacement in humans.

The evolutionary origins and diversity of the neuromuscular system of paired appendages in batoids

Appendage patterning and evolution have been active areas of inquiry for the past two centuries. While most work has centred on the skeleton, particularly that of amniotes, the evolutionary origins and molecular underpinnings of the neuromuscular diversity of fish appendages have remained enigmatic. The fundamental pattern of segmentation in amniotes, for example, is that all muscle precursors and spinal nerves enter either the paired appendages or body wall at the same spinal level. The condition in finned vertebrates is not understood. To address this gap in knowledge, we investigated the development of muscles and nerves in unpaired and paired fins of skates and compared them to those of chain catsharks. During skate and shark embryogenesis, cell populations of muscle precursors and associated spinal nerves at the same axial level contribute to both appendages and body wall, perhaps representing an ancestral condition of gnathostome appendicular neuromuscular systems. Remarkably in skates, this neuromuscular bifurcation as well as colinear Hox expression extend posteriorly to pattern a broad paired fin domain. In addition, we identified migratory muscle precursors (MMPs), which are known to develop into paired appendage muscles with Pax3 and Lbx1 gene expression, in the dorsal fins of skates. Our results suggest that muscles of paired fins have evolved via redeployment of the genetic programme of MMPs that were already involved in dorsal fin development. Appendicular neuromuscular systems most likely have emerged as side branches of body wall neuromusculature and have been modified to adapt to distinct aquatic and terrestrial habitats.

Ancient Use of Ig Variable Domains Contributes Significantly to the TCRδ Repertoire

The loci encoding B and T cell Ag receptors are generally distinct in commonly studied mammals, with each receptor's gene segments limited to intralocus, cis chromosomal rearrangements. The nurse shark (Ginglymostoma cirratum) represents the oldest vertebrate class, the cartilaginous fish, with adaptive immunity provided via Ig and TCR lineages, and is one species among a growing number of taxa employing Ig-TCRδ rearrangements that blend these distinct lineages. Analysis of the nurse shark Ig-TCRδ repertoire found that these rearrangements possess CDR3 characteristics highly similar to canonical TCRδ rearrangements. Furthermore, the Ig-TCRδ rearrangements are expressed with TCRγ, canonically found in the TCRδ heterodimer. We also quantified BCR and TCR transcripts in the thymus for BCR (IgHV-IgHC), chimeric (IgHV-TCRδC), and canonical (TCRδV-TCRδC) transcripts, finding equivalent expression levels in both thymus and spleen. We also characterized the nurse shark TCRαδ locus with a targeted bacterial artifical chromosome sequencing approach and found that the TCRδ locus houses a complex of V segments from multiple lineages. An IgH-like V segment, nestled within the nurse shark TCRδ translocus, grouped with IgHV-like rearrangements we found expressed with TCRδ (but not IgH) rearrangements in our phylogenetic analysis. This distinct lineage of TCRδ-associated IgH-like V segments was termed "TAILVs." Our data illustrate a dynamic TCRδ repertoire employing TCRδVs, NARTCRVs, bona fide trans-rearrangements from shark IgH clusters, and a novel lineage in the TCRδ-associated Ig-like V segments.

Discovery of All Three Types in Cartilaginous Fishes Enables Phylogenetic Resolution of the Origins and Evolution of Interferons

Interferons orchestrate host antiviral responses in jawed vertebrates. They are categorized into three classes; IFN1 and IFN3 are the primary antiviral cytokine lineages, while IFN2 responds to a broader variety of pathogens. The evolutionary relationships within and between these three classes have proven difficult to resolve. Here, we reassess interferon evolution, considering key phylogenetic pitfalls including taxon sampling, alignment quality, model adequacy, and outgroup choice. We reveal that cartilaginous fishes, and hence the jawed vertebrate ancestor, possess(ed) orthologs of all three interferon classes. We show that IFN3 groups sister to IFN1, resolve the origins of the human IFN3 lineages, and find that intronless IFN3s emerged at least three times. IFN2 genes are highly conserved, except for IFN-γ-rel, which we confirm resulted from a teleost-specific duplication. Our analyses show that IFN1 phylogeny is highly sensitive to phylogenetic error. By accounting for this, we describe a new backbone IFN1 phylogeny that implies several IFN1 genes existed in the jawed vertebrate ancestor. One of these is represented by the intronless IFN1s of tetrapods, including mammalian-like repertoires of reptile IFN1s and a subset of amphibian IFN1s, in addition to newly-identified intron-containing shark IFN1 genes. IFN-f, previously only found in teleosts, likely represents another ancestral jawed vertebrate IFN1 family member, suggesting the current classification of fish IFN1s into two groups based on the number of cysteines may need revision. The providence of the remaining fish IFN1s and the coelacanth IFN1s proved difficult to resolve, but they may also be ancestral jawed vertebrate IFN1 lineages. Finally, a large group of amphibian-specific IFN1s falls sister to all other IFN1s and was likely also present in the jawed vertebrate ancestor. Our results verify that intronless IFN1s have evolved multiple times in amphibians and indicate that no one-to-one orthology exists between mammal and reptile IFN1s. Our data also imply that diversification of the multiple IFN1s present in the jawed vertebrate ancestor has occurred through a rapid birth-death process, consistent with functional maintenance over a 450-million-year host-pathogen arms race. In summary, this study reveals a new model of interferon evolution important to our understanding of jawed vertebrate antiviral immunity.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

Mitral cell development in the olfactory bulb of sharks: evidences of a conserved pattern of glutamatergic neurogenesis

In mammals, the development of the olfactory bulb (OB) relies in part on the expression of transcription factors involved in the specifications/differentiation of glutamatergic cells. In a previous study from our group, a high molecular similarity was reported between mammals and cartilaginous fishes regarding the neurogenic mechanisms underlying the development of glutamatergic cells in the telencephalon. However, information about the transcriptional program operating in the development of the glutamatergic system (mainly represented by mitral cells) in the OB is lacking in the catshark Scyliorhinus canicula, a cartilaginous fish. Using immunohistochemistry and in situ hybridization techniques, we have found that, previously to the appearance of the olfactory primordium (OP), proliferating cells expressing Pax6 with molecular hallmarks of progenitor radial glia were located in the ventrolateral pallial ventricular zone. Later in development, when the OP is recognizable, a stream of Pax6-positive cells were observed between the ventricular zone and the OP, where transcription factors involved in mitral cell development in mammals (ScTbr2, ScNeuroD, Tbr1) are expressed. Later in development, these transcription factors became expressed in a layered-like structure where ScVglut1, a marker of mitral cells, is also present. Our data suggest that the transcriptional program related with the specification/differentiation of glutamatergic cells in the telencephalon has been conserved throughout the evolution of vertebrates. These results, in combination with previous studies concerning GABAergic neurogenesis in sharks, have evidenced that the OB of mammals and sharks shares similarities in the timing and molecular programs of development.

Differential expression of myosin heavy chain isoforms in cardiac segments of gnathostome vertebrates and its evolutionary implications

Background

Immunohistochemical studies of hearts from the lesser spotted dogfish, Scyliorhinus canicula (Chondrichthyes) revealed that the pan-myosin heavy chain (pan-MyHC) antibody MF20 homogeneously labels all the myocardium, while the pan-MyHC antibody A4.1025 labels the myocardium of the inflow (sinus venosus and atrium) but not the outflow (ventricle and conus arteriosus) cardiac segments, as opposed to other vertebrates. We hypothesized that the conventional pattern of cardiac MyHC isoform distribution present in most vertebrates, i.e. MYH6 in the inflow and MYH7 in the outflow segments, has evolved from a primitive pattern that persists in Chondrichthyes. In order to test this hypothesis, we conducted protein detection techniques to identify the MyHC isoforms expressed in adult dogfish cardiac segments and to assess the pan-MyHC antibodies reactivity against the cardiac segments of representative species from different vertebrate groups.

Results

Western and slot blot results confirmed the specificity of MF20 and A4.1025 for MyHC in dogfish and their differential reactivity against distinct myocardial segments. HPLC-ESI-MS/MS and ESI-Quadrupole-Orbitrap revealed abundance of MYH6 and MYH2 in the inflow and of MYH7 and MYH7B in the outflow segments. Immunoprecipitation showed higher affinity of A4.1025 for MYH2 and MYH6 than for MYH7 and almost no affinity for MYH7B. Immunohistochemistry showed that A4.1025 signals are restricted to the inflow myocardial segments of elasmobranchs, homogeneous in all myocardial segments of teleosts and acipenseriforms, and low in the ventricle of polypteriforms.

Conclusions

The cardiac inflow and outflow segments of the dogfish show predominance of fast- and slow-twitch MyHC isoforms respectively, what can be considered a synapomorphy of gnathostomes. The myocardium of the dogfish contains two isomyosins (MYH2 and MYH7B) not expressed in the adult heart of other vertebrates. We propose that these isomyosins lost their function in cardiac contraction during the evolution of gnathostomes, the later acquiring a regulatory role in myogenesis through its intronic miRNA. Loss of MYH2 and MYH7B expression in the heart possibly occurred before the origin of Osteichthyes, being the latter reacquired in polypteriforms. We raise the hypothesis that the slow tonic MYH7B facilitates the peristaltic contraction of the conus arteriosus of fish with a primitive cardiac anatomical design and of the vertebrate embryo.

Evolutionary Plasticity in Detoxification Gene Modules: The Preservation and Loss of the Pregnane X Receptor in Chondrichthyes Lineages

To appraise how evolutionary processes, such as gene duplication and loss, influence an organism's xenobiotic sensitivity is a critical question in toxicology. Of particular importance are gene families involved in the mediation of detoxification responses, such as members of the nuclear receptor subfamily 1 group I (NR1I), the pregnane X receptor (PXR), and the constitutive androstane receptor (CAR). While documented in multiple vertebrate genomes, PXR and CAR display an intriguing gene distribution. PXR is absent in birds and reptiles, while CAR shows a tetrapod-specific occurrence. More elusive is the presence of PXR and CAR gene orthologs in early branching and ecologically-important Chondrichthyes (chimaeras, sharks and rays). Therefore, we investigated various genome projects and use them to provide the first identification and functional characterization of a Chondrichthyan PXR from the chimaera elephant shark (Callorhinchus milii, Holocephali). Additionally, we substantiate the targeted PXR gene loss in Elasmobranchii (sharks and rays). Compared to other vertebrate groups, the chimaera PXR ortholog displays a diverse expression pattern (skin and gills) and a unique activation profile by classical xenobiotic ligands. Our findings provide insights into the molecular landscape of detoxification mechanisms and suggest lineage-specific adaptations in response to xenobiotics in gnathostome evolution.

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

'Genome skimming' with the MinION hand-held sequencer identifies CITES-listed shark species in India's exports market

Chondrichthyes - sharks, rays, skates, and chimeras, are among the most threatened and data deficient vertebrate species. Global demand for shark and ray derived products, drives unregulated and exploitative fishing practices, which are in turn facilitated by the lack of ecological data required for effective conservation of these species. Here, we describe a Next Generation Sequencing method (using the MinION, a hand-held portable sequencing device from Oxford Nanopore Technologies), and analyses pipeline for molecular ecological studies in Chondrichthyes. Using this method, the complete mitochondrial genome and nuclear intergenic and protein-coding sequences were obtained by direct sequencing of genomic DNA obtained from shark fin tissue. Recovered loci include mitochondrial barcode sequences- Cytochrome oxidase I, NADH2, 16S rRNA and 12S rRNA- and nuclear genetic loci such as 5.8S rRNA, Internal Transcribed Spacer 2, and 28S rRNA regions, which are commonly used for taxonomic identification. Other loci recovered were the nuclear protein-coding genes for antithrombin or SerpinC, Immunoglobulin lambda light chain, Preprogehrelin, selenium binding protein 1(SBP1), Interleukin-1 beta (IL-1β) and Recombination-Activating Gene 1 (RAG1). The median coverage across all genetic loci was 20x and sequence accuracy was ≥99.8% compared to reference sequences. Analyses of the nuclear ITS2 region and the mitochondrial protein-encoding loci allowed accurate taxonomic identification of the shark specimen as Carcharhinus falciformis, a CITES Appendix II species. MinION sequencing provided 1,152,211 bp of new shark genome, increasing the number of sequenced shark genomes to five. Phylogenetic analyses using both mitochondrial and nuclear loci provided evidence that Prionace glauca is nested within Carcharhinus, suggesting the need for taxonomic reassignment of P. glauca. We increased genomic information about a shark species for ecological and population genetic studies, enabled accurate identification of the shark tissue for biodiversity indexing and resolved phylogenetic relationships among multiple taxa. The method was independent of amplification bias, and adaptable for field assessments of other Chondrichthyes and wildlife species in the future.

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.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Regulation of Gastrointestinal Motility by Motilin and Ghrelin in Vertebrates

The energy balance of vertebrates is regulated by the difference in energy input and energy expenditure. Generally, most vertebrates obtain their energy from nutrients of foods through the gastrointestinal (GI) tract. Therefore, food intake and following food digestion, including motility of the GI tract, secretion and absorption, are crucial physiological events for energy homeostasis. GI motility changes depending on feeding, and GI motility is divided into fasting (interdigestive) and postprandial (digestive) contraction patterns. GI motility is controlled by contractility of smooth muscles of the GI tract, extrinsic and intrinsic neurons (motor and sensory) and some hormones. In mammals, ghrelin (GHRL) and motilin (MLN) stimulate appetite and GI motility and contribute to the regulation of energy homeostasis. GHRL and MLN are produced in the mucosal layer of the stomach and upper small intestine, respectively. GHRL is a multifunctional peptide and is involved in glucose metabolism, endocrine/exocrine functions and cardiovascular and reproductive functions, in addition to feeding and GI motility in mammals. On the other hand, the action of MLN is restricted and species such as rodentia, including mice and rats, lack MLN peptide and its receptor. From a phylogenetic point of view, GHRL and its receptor GHS-R1a have been identified in various vertebrates, and their structural features and various physiological functions have been revealed. On the other hand, MLN or MLN-like peptide (MLN-LP) and its receptors have been found only in some fish, birds and mammals. Here, we review the actions of GHRL and MLN with a focus on contractility of the GI tract of species from fish to mammals.

Evaluating Genome Assemblies and Gene Models Using gVolante

In daily practice of de novo genome assembly and gene prediction, it would be a natural urge to evaluate their products. Different programs and parameter settings give rise to variable outputs, which leaves a decision of which output to adopt for downstream analysis for addressing biological questions. Instead of superficial assessment of length-based statistics of output sequences (e.g., N50 scaffold length), completeness assessment by means of scoring the coverage of reference orthologs has been increasingly utilized.We previously launched a web service, gVolante ( https://gvolante.riken.jp /), to provide a user-friendly interface and a uniform environment for completeness assessment with the pipelines CEGMA and BUSCO. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths, allowing quality control in multiple aspects. This chapter focuses on the procedure for such assessment and provides technical tips for higher accuracy.

De novo transcriptome analyses provide insights into opsin-based photoreception in the lanternshark Etmopterus spinax

The velvet belly lanternshark (Etmopterus spinax) is a small deep-sea shark commonly found in the Eastern Atlantic and the Mediterranean Sea. This bioluminescent species is able to emit a blue-green ventral glow used in counter-illumination camouflage, mainly. In this study, paired-end Illumina HiSeqTM technology has been employed to generate transcriptome data from eye and ventral skin tissues of the lanternshark. About 64 and 49 million Illumina reads were generated from skin and eye tissues respectively. The assembly allowed us to predict 119,749 total unigenes including 94,569 for the skin transcriptome and 94,365 for the eye transcriptome while 74,753 were commonly found in both transcriptomes. A taxonomy filtering was applied to extract a reference transcriptome containing 104,390 unigenes among which 38,836 showed significant similarities to known sequences in NCBI non-redundant protein sequences database. Around 58% of the annotated unigenes match with predicted genes from the Elephant shark (Callorhinchus milii) genome. The transcriptome completeness has been evaluated by successfully capturing around 98% of orthologous genes of the « Core eukaryotic gene dataset » within the E. spinax reference transcriptome. We identified potential "light-interacting toolkit" genes including multiple genes related to ocular and extraocular light perception processes such as opsins, phototransduction actors or crystallins. Comparative gene expression analysis reveals eye-specific expression of opsins, ciliary phototransduction actors, crystallins and vertebrate retinoid pathway actors. In particular, mRNAs from a single rhodopsin gene and its potentially associated peropsin were detected in the eye transcriptome, only, confirming a monochromatic vision of the lanternshark. Encephalopsin mRNAs were mainly detected in the ventral skin transcriptome. In parallel, immunolocalization of the encephalopsin within the ventral skin of the shark suggests a functional relation with the photophores, i.e. epidermal light-producing organs. We hypothesize that extraocular photoreception might be involved in the bioluminescence control possibly acting on the shutter opening and/or the photocyte activity itself. The newly generated reference transcriptome provides a valuable resource for further understanding of the shark biology.

Potential Human Health Applications from Marine Biomedical Research with Elasmobranch Fishes

Members of the subclass of fishes collectively known as elasmobranchs (Class Chondrichthyes, Subclass Elasmobranchii) include sharks, skates, rays, guitarfish, and sawfish. Having diverged from the main line of vertebrate evolution some 400 million years ago, these fishes have continued to be successful in our ever-changing oceans. Much of their success must be attributed to their uncanny ability to remain healthy. Based on decades of basic research, some of their secrets may be very close to benefitting man. In this short review, some of the molecular and cellular biological areas that show promise for potential human applications are presented. With a brief background and current status of relevant research, these topics include development of new antibiotics and novel treatments for cancer, macular degeneration, viral pathogens, and Parkinson’s disease; potentially useful genomic information from shark transcriptomes; shark antibody-derived drug delivery systems; and immune cell-derived compounds as potential cancer therapeutic agents.