Transcriptomic analysis of the lesser spotted catshark (Scyliorhinus canicula) pancreas, liver and brain reveals molecular level conservation of vertebrate pancreas function

The chondrichthyan glucagon-like peptide 3 regulates hepatic ketone metabolism in the Pacific spiny dogfish Squalus suckleyi

Chondrichthyans have a novel proglucagon-derived peptide, glucagon-like peptide (GLP)-3, in addition to GLP-1 and GLP-2 that occur in other vertebrates. Given that the GLPs are important regulators of metabolic homeostasis across vertebrates, we sought to investigate whether GLP-3 displays functional actions on metabolism within a representative chondrichthyan, the Pacific spiny dogfish Squalus suckleyi. There were no observed effects of GLP-3 perfusion (10 nM for 15 min) on the rate of glucose or oleic acid acquisition at the level of the spiral valve nor were there any measured effects on intermediary metabolism within this tissue. Despite no effects on apparent glucose transport or glycolysis in the liver, a significant alteration to ketone metabolism occurred. Firstly, ketone flux through the perfused liver switched from a net endogenous production to consumption following hormone application. Accompanying this change, significant increases in mRNA transcript abundance of putative ketone transporters and in the activity of β-hydroxybutyrate dehydrogenase (a key enzyme regulating ketone flux in the liver) were observed. Overall, while these results show effects on hepatic metabolism, the physiological actions of GLP are distinct between this chondrichthyan and those of GLP-1 on teleost fishes. Whether this is the result of the particular metabolic dependency on ketone bodies in chondrichthyans or a differential function of a novel GLP remains to be fully elucidated.

Comparative Proteome Analysis of Four Stages of Spermatogenesis in the Small-Spotted Catshark (Scyliorhinus canicula), Using High-Resolution NanoLC-ESI-MS/MS

Spermatogenesis is a highly specialized process of cell proliferation and differentiation leading to the production of spermatozoa from spermatogonial stem cells. Due to its testicular anatomy, Scyliorhinus canicula is an interesting model to explore stage-based changes in proteins during spermatogenesis. The proteomes of four testicular zones corresponding to the germinative niche and to spermatocysts (cysts) with spermatogonia (zone A), cysts with spermatocytes (zone B), cysts with young spermatids (zone C), and cysts with late spermatids (zone D) have been analyzed by nanoLC-ESI-MS/MS. Gene ontology and KEGG annotations were also performed. A total of 3346 multiple protein groups were identified. Zone-specific protein analyses highlighted RNA-processing, chromosome-related processes, cilium organization, and cilium activity in zones A, D, C, and D, respectively. Analyses of proteins with zone-dependent abundance revealed processes related to cellular stress, ubiquitin-dependent degradation by the proteasome, post-transcriptional regulation, and regulation of cellular homeostasis. Our results also suggest that the roles of some proteins, such as ceruloplasmin, optineurin, the pregnancy zone protein, PA28β or the Culling-RING ligase 5 complex, as well as some uncharacterized proteins, during spermatogenesis could be further explored. Finally, the study of this shark species allows one to integrate these data in an evolutionary context of the regulation of spermatogenesis. Mass spectrometry data are freely accessible via iProX-integrated Proteome resources (https://www.iprox.cn/) for reuse purposes.

Sharks Provide Evidence for a Highly Complex TNFSF Repertoire in the Jawed Vertebrate Ancestor

Cytokines of the TNF superfamily (TNFSF) control many immunological processes and are implicated in the etiology of many immune disorders and diseases. Despite their obvious biological importance, the TNFSF repertoires of many species remain poorly characterized. In this study, we perform detailed bioinformatic, phylogenetic, and syntenic analyses of five cartilaginous fish genomes to identify their TNFSF repertoires. Strikingly, we find that shark genomes harbor ∼30 TNFSF genes, more than any other vertebrate examined to date and substantially more than humans. This is due to better retention of the ancestral jawed vertebrate TNFSF repertoire than any other jawed vertebrate lineage, combined with lineage-specific gene family expansions. All human TNFSFs appear in shark genomes, except for lymphotoxin-α (LTA; TNFSF1) and TNF (TNFSF2), and CD70 (TNFSF7) and 4-1BBL (TNFSF9), which diverged by tandem duplications early in tetrapod and mammalian evolution, respectively. Although lacking one-to-one LTA and TNF orthologs, sharks have evolved lineage-specific clusters of LTA/TNF co-orthologs. Other key findings include the presence of two BAFF (TNFSF13B) genes along with orthologs of APRIL (TNFSF13) and BALM (TNFSF13C) in sharks, and that all cartilaginous fish genomes harbor an ∼400-million-year-old cluster of multiple FASLG (TNFSF6) orthologs. Finally, sharks have retained seven ancestral jawed vertebrate TNFSF genes lost in humans. Taken together, our data indicate that the jawed vertebrate ancestor possessed a much larger and diverse TNFSF repertoire than previously hypothesized and oppose the idea that the cartilaginous fish immune system is “primitive” compared with that of mammals.

Evolution of calcitonin/calcitonin gene-related peptide family in chordates: Identification of CT/CGRP family peptides in cartilaginous fish genome

The calcitonin (CT)/CT gene-related peptide (CGRP) family is a peptide gene family that is widely found in bilaterians. CT, CGRP, adrenomedullin (AM), amylin (AMY), and CT receptor-stimulating peptide (CRSP) are members of the CT/CGRP family. In mammals, CT is involved in calcium homeostasis, while CGRP and AM primarily function in vasodilation. AMY and CRSP are associated with anorectic effects. Diversification of the molecular features and physiological functions of the CT/CGRP family in vertebrate lineages have been extensively reported. However, the origin and diversification mechanisms of the vertebrate CT/CGRP family of peptides remain unclear. In this review, the molecular characteristics of CT/CGRP family peptides and their receptors, along with their major physiological functions in mammals and teleosts, are introduced. Furthermore, novel candidates of the CT/CGRP family in cartilaginous fish are presented based on genomic information. The CT/CGRP family peptides and receptors in urochordates and cephalochordates, which are closely related to vertebrates, are also described. Finally, a putative evolutionary scenario of the CT/CGRP family peptides and receptors in chordates is discussed.

Sustained endocrine and exocrine function in the pancreas of the Pacific spiny dogfish post-feeding

Secretions of the exocrine pancreas contain digestive enzymes integral to the digestive process. The Pacific spiny dogfish (Squalus suckleyi) has a discrete pancreas, divided into two lobes termed the dorsal and ventral lobes. These lobes drain into the anterior intestine via a common duct to enable digestion. Previous studies have identified that the exocrine pancreas produces (co)lipases, chymotrypsin, carboxypeptidase, and low levels of chitinases; however, investigations into other digestive enzymes are limited. We detect the presence of lipase, trypsin, and carbohydrase and show that activities are equivalent between both lobes of the pancreas. Additionally, we sought to investigate the influence of a single feeding event (2% body weight ration of herring by gavage) on enzyme activities over an extended time course (0, 20, 48, 72, 168 h) post-feeding. The results indicate that there are no differences in pancreatic tissue digestive enzyme activities between fed or fasted states. Analysis of acinar cell circumference post-feeding demonstrates a significant increase at 20 and 48 h, that returns to fasting levels by 72 h. No significant changes were observed regarding whole-tissue insulin or glucagon mRNA abundance or with glucose transporter (glut) 1 or 3. Yet, a significant and transient decrease in glut4 and sodium glucose-linked transporter mRNA abundance was found at 48 h post-feeding. We propose that the constant enzyme activity across this relatively large organ, in combination with a relatively slow rate of digestion leads to an evenly distributed, sustained release of digestive enzymes regardless of digestive state.

Retrospective genomics highlights changes in genetic composition of tiger sharks (Galeocerdo cuvier) and potential loss of a south-eastern Australia population

Over the last century, many shark populations have declined, primarily due to overexploitation in commercial, artisanal and recreational fisheries. In addition, in some locations the use of shark control programs also has had an impact on shark numbers. Still, there is a general perception that populations of large ocean predators cover wide areas and therefore their diversity is less susceptible to local anthropogenic disturbance. Here we report on temporal genomic analyses of tiger shark (Galeocerdo cuvier) DNA samples that were collected from eastern Australia over the past century. Using Single Nucleotide Polymorphism (SNP) loci, we documented a significant change in genetic composition of tiger sharks born between ~1939 and 2015. The change was most likely due to a shift over time in the relative contribution of two well-differentiated, but hitherto cryptic populations. Our data strongly indicate a dramatic shift in the relative contribution of these two populations to the overall tiger shark abundance on the east coast of Australia, possibly associated with differences in direct or indirect exploitation rates.

The Evolution of Oxytocin and Vasotocin Receptor Genes in Jawed Vertebrates: A Clear Case for Gene Duplications Through Ancestral Whole-Genome Duplications

The neuronal and neuroendocrine peptides oxytocin (OT) and vasotocin (VT), including vasopressins, have six cognate receptors encoded by six receptor subtype genes in jawed vertebrates. The peptides elicit a broad range of responses that are specifically mediated by the receptor subtypes including neuronal functions regulating behavior and hormonal actions on reproduction and water/electrolyte balance. Previously, we have demonstrated that these six receptor subtype genes, which we designated VTR1A, VTR1B, OTR, VTR2A, VTR2B and VTR2C, arose from a syntenic ancestral gene pair, one VTR1/OTR ancestor and one VTR2 ancestor, through the early vertebrate whole-genome duplications (WGD) called 1R and 2R. This was supported by both phylogenetic and chromosomal conserved synteny data. More recently, other studies have focused on confounding factors, such as the OTR/VTR orthologs in cyclostomes, to question this scenario for the origin of the OTR/VTR gene family; proposing instead less parsimonious interpretations involving only one WGD followed by complex series of chromosomal or segmental duplications. Here, we have updated the phylogeny of the OTR/VTR gene family, including a larger number of vertebrate species, and revisited seven representative neighboring gene families from our previous conserved synteny analyses, adding chromosomal information from newer high-coverage genome assemblies from species that occupy key phylogenetic positions: the polypteriform fish reedfish (Erpetoichthys calabaricus), the cartilaginous fish thorny skate (Amblyraja radiata) and a more recent high-quality assembly of the Western clawed frog (Xenopus tropicalis) genome. Our analyses once again add strong support for four-fold symmetry, i.e., chromosome quadruplication in the same time window as the WGD events early in vertebrate evolution, prior to the jawed vertebrate radiation. Thus, the evolution of the OTR/VTR gene family can be most parsimoniously explained by two WGD events giving rise to the six ancestral genes, followed by differential gene losses of VTR2 genes in different lineages. We also argue for more coherence and clarity in the nomenclature of OT/VT receptors, based on the most parsimonious scenario.

Comparative eye and liver differentially expressed genes reveal monochromatic vision and cancer resistance in the shortfin mako shark (Isurus oxyrinchus)

The shortfin mako, Isurus oxyrinchus is an oceanic pelagic shark found worldwide in tropical and subtropical waters. However, the understanding of its biology at molecular level is still incipient. We sequenced the messenger RNA isolated from eye and liver tissues. De novo transcriptome yielded a total of 705,940 transcripts. A total of 3774 genes were differentially expressed (DEGs), with 1612 in the eye and 2162 in the liver. Most DEGs in the eye were related to structural and signaling functions, including nonocular and ocular opsin genes, whereas nine out of ten most overexpressed genes in the liver were related to tumor suppression, wound healing, and human diseases. Furthermore, DEGs findings provide insights on the monochromatic shark vision and a repertory of cancer-related genes, which may be insightful to elucidate shark resistance to cancer. Therefore, our results provide valuable sequence resources for future functional and population studies.

Lampreys, the jawless vertebrates, contain three Pax6 genes with distinct expression in eye, brain and pancreas

The transcription factor Pax6 is crucial for the development of the central nervous system, eye, olfactory system and pancreas, and is implicated in human disease. While a single Pax6 gene exists in human and chicken, Pax6 occurs as a gene family in other vertebrates, with two members in elephant shark, Xenopus tropicalis and Anolis lizard and three members in teleost fish such as stickleback and medaka. However, the complement of Pax6 genes in jawless vertebrates (cyclostomes), the sister group of jawed vertebrates (gnathostomes), is unknown. Using a combination of BAC sequencing and genome analysis, we discovered three Pax6 genes in lampreys. Unlike the paired-less Pax6 present in some gnathostomes, all three lamprey Pax6 have a highly conserved full-length paired domain. All three Pax6 genes are expressed in the eye and brain, with variable expression in other tissues. Notably, lamprey Pax6α transcripts are found in the pancreas, a vertebrate-specific organ, indicating the involvement of Pax6 in development of the pancreas in the vertebrate ancestor. Multi-species sequence comparisons revealed only a single conserved non-coding element, in the lamprey Pax6β locus, with similarity to the PAX6 neuroretina enhancer. Using a transgenic zebrafish enhancer assay we demonstrate functional conservation of this element over 500 million years of vertebrate evolution.

Transcriptome-Derived Microsatellites Demonstrate Strong Genetic Differentiation in Pacific White Sharks

Recent advances in genome-scale sequencing technology have allowed the development of high resolution genetic markers for the study of nonmodel taxa. In particular, transcriptome sequencing has proven to be highly useful in generating genomic markers for use in population genetic studies, allowing for insight into species connectivity, as well as local adaptive processes as many transcriptome-derived markers are found within or associated with functional genes. Herein, we developed a set of 30 microsatellite markers from a heart transcriptome for the white shark (Carcharodon carcharias), a widely distributed and globally vulnerable marine predator. Using these markers as well as 10 published anonymous genomic microsatellite loci, we provide 1) the first nuclear genetic assessment of the cross-Pacific connectivity of white sharks, and 2) a comparison of the levels of inferred differentiation across microsatellite marker sets (i.e., transcriptome vs. anonymous) to assess their respective utility to elucidate the population genetic dynamics of white sharks. Significant (FST = 0.083, P = 0.05; G″ST = 0.200; P = 0.001) genetic differentiation was found between Southwestern Pacific (n = 19) and Northeastern Pacific (n = 20) white sharks, indicating restricted, cross Pacific gene flow in this species. Transcriptome-derived microsatellite marker sets identified much higher (up to 2×) levels of genetic differentiation than anonymous genomic markers, underscoring potential utility of transcriptome markers in identifying subtle population genetic differences within highly vagile, globally distributed marine species.Subject areas: Population structure and phylogeography; Conservation genetics and biodiversity.

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.

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.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

De novo assembly of the kidney and spleen transcriptomes of the cosmopolitan blue shark, Prionace glauca

Cartilaginous fishes (sharks, rays and chimaeras) comprise a highly diversified group of basal vertebrates occupying a plethora of ecological aquatic niches. They represent critical components of marine ecosystems and food webs, although numerous species are threatened and almost half are poorly known. Genomic resources emerging from this basal jawed vertebrate group have offered valuable insights into the evolution of vertebrate-specific traits. Yet, the taxon remains largely understudied. Here, we generated the first high-quality de novo assembly of kidney and spleen transcriptomes of the blue shark (Prionace glauca). A total of 32,917,412 and 52,666,542 reads were obtained for spleen and kidney, respectively, using RNA-Seq Illumina technology. De novo multi-tissue assembly resulted in 97,317 unigenes with an N50 of 1975bp, in which 87,571 were assigned to a particular tissue or combination of tissues based on the sequencing read mapping. Functional annotation generated 28,564 and 19,854 open reading frames in spleen and kidney, respectively. This dataset provides a significant resource for physiological and evolutionary studies, namely into the unique osmoregulatory system of Chondrichthyes and the evolution of the immune system in vertebrates.

The Calcitonin/Calcitonin Gene-Related Peptide Family in Invertebrate Deuterostomes

Calcitonin (CT)/CT gene-related peptide (CGRP) family peptides (CT/CGRP family peptides) including CT, CGRP, adrenomedullin, amylin, and CT receptor-stimulating peptide have been identified from various vertebrates and perform a variety of important physiological functions. These peptides bind to two types of receptors including CT receptor (CTR) and CTR-like receptor (CLR). Receptor recognition of CT/CGRP family peptides is determined by the heterodimer between CTR/CLR and receptor activity-modifying protein (RAMP). Comparative studies of the CT/CGRP family have been exclusively performed in vertebrates from teleost fishes to mammals and strongly manifest that the CGRP family system containing peptides, their receptors, and RAMPs was derived from a common ancestor. In addition, CT/CGRP family peptides and their receptors are also identified and inferred from various invertebrate species. However, the evolutionary process of the CT/CGRP family from invertebrates to vertebrates remains enigmatic. In this review, I principally summarize the CT/CGRP family peptides and their receptors in invertebrate deuterostomes, highlighting the study of invertebrate chordates including ascidians and amphioxi. The CT/CGRP family peptide that shows similar molecular structure and function with that of vertebrate CT has been identified from ascidian, Ciona intestinalis. Amphioxus, Branchiostoma floridae also possessed three CT/CGRP family peptides, one CTR/CLR receptor, and three RAMP-like proteins. The molecular function of the receptor complex formed by amphioxus CTR/CLR and a RAMP-like protein was clarified. Moreover, CT/CGRP family peptides have been identified in the superphylum Ambulacraria, which is close to Chordata. Finally, this review provides potential hypotheses of the evolution of CGRP family peptides and their receptors from invertebrates to vertebrates.

Evolutionary history of the T cell receptor complex as revealed by small-spotted catshark (Scyliorhinus canicula)

In every jawed vertebrate species studied so far, the T cell receptor (TCR) complex is composed of two different TCR chains (α/β or γ/δ) and a number of CD3 subunits responsible for transmitting signals into the T cell. In this study, we characterised all of the TCR and CD3 genes of small-spotted catshark (Scyliorhinus canicula) and analysed their expression in a broad range of tissues. While the TCR complex is highly conserved across jawed vertebrates, we identified a number of differences in catshark, most notably the presence of two copies of both TCRβ and CD3γδ, and the absence of a functionally-important proline rich region from CD3ε. We also demonstrate that TCRβ has duplicated independently multiple times in jawed vertebrate evolution, bringing additional diversity to the TCR complex. This study reveals new insights about the evolutionary history of the TCR complex and raises new avenues for future exploration.

Multi-tissue RNA-seq and transcriptome characterisation of the spiny dogfish shark (Squalus acanthias) provides a molecular tool for biological research and reveals new genes involved in osmoregulation

The spiny dogfish shark (Squalus acanthias) is one of the most commonly used cartilaginous fishes in biological research, especially in the fields of nitrogen metabolism, ion transporters and osmoregulation. Nonetheless, transcriptomic data for this organism is scarce. In the present study, a multi-tissue RNA-seq experiment and de novo transcriptome assembly was performed in four different spiny dogfish tissues (brain, liver, kidney and ovary), providing an annotated sequence resource. The characterization of the transcriptome greatly increases the scarce sequence information for shark species. Reads were assembled with the Trinity de novo assembler both within each tissue and across all tissues combined resulting in 362,690 transcripts in the combined assembly which represent 289,515 Trinity genes. BUSCO analysis determined a level of 87% completeness for the combined transcriptome. In total, 123,110 proteins were predicted of which 78,679 and 83,164 had significant hits against the SwissProt and Uniref90 protein databases, respectively. Additionally, 61,215 proteins aligned to known protein domains, 7,208 carried a signal peptide and 15,971 possessed at least one transmembrane region. Based on the annotation, 81,582 transcripts were assigned to gene ontology terms and 42,078 belong to known clusters of orthologous groups (eggNOG). To demonstrate the value of our molecular resource, we show that the improved transcriptome data enhances the current possibilities of osmoregulation research in spiny dogfish by utilizing the novel gene and protein annotations to investigate a set of genes involved in urea synthesis and urea, ammonia and water transport, all of them crucial in osmoregulation. We describe the presence of different gene copies and isoforms of key enzymes involved in this process, including arginases and transporters of urea and ammonia, for which sequence information is currently absent in the databases for this model species. The transcriptome assemblies and the derived annotations generated in this study will support the ongoing research for this particular animal model and provides a new molecular tool to assist biological research in cartilaginous fishes.

Lampreys, the jawless vertebrates, contain only two ParaHox gene clusters

ParaHox genes (Gsx, Pdx, and Cdx) are an ancient family of developmental genes closely related to the Hox genes. They play critical roles in the patterning of brain and gut. The basal chordate, amphioxus, contains a single ParaHox cluster comprising one member of each family, whereas nonteleost jawed vertebrates contain four ParaHox genomic loci with six or seven ParaHox genes. Teleosts, which have experienced an additional whole-genome duplication, contain six ParaHox genomic loci with six ParaHox genes. Jawless vertebrates, represented by lampreys and hagfish, are the most ancient group of vertebrates and are crucial for understanding the origin and evolution of vertebrate gene families. We have previously shown that lampreys contain six Hox gene loci. Here we report that lampreys contain only two ParaHox gene clusters (designated as α- and β-clusters) bearing five ParaHox genes (Gsxα, Pdxα, Cdxα, Gsxβ, and Cdxβ). The order and orientation of the three genes in the α-cluster are identical to that of the single cluster in amphioxus. However, the orientation of Gsxβ in the β-cluster is inverted. Interestingly, Gsxβ is expressed in the eye, unlike its homologs in jawed vertebrates, which are expressed mainly in the brain. The lamprey Pdxα is expressed in the pancreas similar to jawed vertebrate Pdx genes, indicating that the pancreatic expression of Pdx was acquired before the divergence of jawless and jawed vertebrate lineages. It is likely that the lamprey Pdxα plays a crucial role in pancreas specification and insulin production similar to the Pdx of jawed vertebrates.

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.

The Neuroendocrine Regulation of Food Intake in Fish: A Review of Current Knowledge

Fish are the most diversified group of vertebrates and, although progress has been made in the past years, only relatively few fish species have been examined to date, with regards to the endocrine regulation of feeding in fish. In fish, as in mammals, feeding behavior is ultimately regulated by central effectors within feeding centers of the brain, which receive and process information from endocrine signals from both brain and peripheral tissues. Although basic endocrine mechanisms regulating feeding appear to be conserved among vertebrates, major physiological differences between fish and mammals and the diversity of fish, in particular in regard to feeding habits, digestive tract anatomy and physiology, suggest the existence of fish- and species-specific regulating mechanisms. This review provides an overview of hormones known to regulate food intake in fish, emphasizing on major hormones and the main fish groups studied to date.

Extracting DNA from 'jaws': high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so-called bio-swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high-yield sources of DNA for genomic-scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.

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.

Alternative adaptive immunity strategies: coelacanth, cod and shark immunity

The advent of high throughput sequencing has permitted to investigate the genome and the transcriptome of novel non-model species with unprecedented depth. This technological advance provided a better understanding of the evolution of adaptive immune genes in gnathostomes, revealing several unexpected features in different fish species which are of particular interest. In the present paper, we review the current understanding of the adaptive immune system of the coelacanth, the elephant shark and the Atlantic cod. The study of coelacanth, the only living extant of the long thought to be extinct Sarcopterygian lineage, is fundamental to bring new insights on the evolution of the immune system in higher vertebrates. Surprisingly, coelacanths are the only known jawed vertebrates to lack IgM, whereas two IgD/W loci are present. Cartilaginous fish are of great interest due to their basal position in the vertebrate tree of life; the genome of the elephant shark revealed the lack of several important immune genes related to T cell functions, which suggest the existence of a primordial set of TH1-like cells. Finally, the Atlantic cod lacks a functional major histocompatibility II complex, but balances this evolutionary loss with the expansion of specific gene families, including MHC I, Toll-like receptors and antimicrobial peptides. Overall, these data point out that several fish species present an unconventional adaptive immune system, but the loss of important immune genes is balanced by adaptive evolutionary strategies which still guarantee the establishment of an efficient immune response against the pathogens they have to fight during their life.