The challenge of implementing environmental DNA metabarcoding to detect elasmobranchs in a resource-limited marine protected area.

Elasmobranchs are threatened and eDNA metabarcoding is a powerful tool that can help efforts to better understand and conserve them. Nevertheless, the inter-calibration between optimal methodological practices and its implementation in resource-limited situations is still an issue. Based on promising results from recent studies, the authors applied a cost-effective protocol with parameters that could be easily replicated by any conservationist. Nonetheless, the results with fewer elasmobranchs detected than expected reveal that endorsed primers and sampling strategies still require further optimization, especially for applications in resource-limited conservation programmes.

Evidence of multiple paternity for the endemic Xingu River stingray.

Multiple paternity (MP) is a phenomenon observed for more than 30 elasmobranch species. The Batoidea is more specious than the Selachii, but only three studies of MP have been conducted on batoids. The occurrence of MP in freshwater stingrays was tested using microsatellite markers, which were developed for Potamotrygon leopoldi. Six mothers and their litters were genotyped, providing the first evidence of MP for Potamotrygonidae, with an MP frequency of 33%.

The Nursehound Scyliorhinus stellaris Mitochondrial Genome-Phylogeny, Relationships among Scyliorhinidae and Variability in Waters of the Balearic Islands.

The complete mitochondrial DNA sequence of the Nursehound Scyliorhinus stellaris has been determined for the first time and compared with congeneric species. The mitogenome sequence was 16,684 bp in length. The mitogenome is composed of 13 PCGs, 2 rRNAs, 22 transfer RNA genes and non-coding regions. The gene order of the newly sequenced mitogenome is analogous to the organization described in other vertebrate genomes. The typical conservative blocks in the control region were indicated. The phylogenetic analysis revealed a monophyletic origin of the Scyliorhininae subfamily, and within it, two subclades were identified. A significant divergence of Scyliorhinus spp. together with Poroderna patherinum in relation to the group of Cephaloscyllium spp. was observed, except for Scyliorhinus torazame, more related to this last cited clade. A hypothesis of a divergent evolution consequent to a selective pressure in different geographic areas, which lead to a global latitudinal diversity gradient, has been suggested to explain this phylogenetic reconstruction. However, convergent evolution on mitochondrial genes could also involve different species in some areas of the world.

Evolutionary developmental genetics of teeth and odontodes in jawed vertebrates: a perspective from the study of elasmobranchs.

Most extant vertebrates display a high variety of tooth and tooth-like organs (odontodes) that vary in shape, position over the body and nature of composing tissues. The development of these structures is known to involve similar genetic cascades and teeth and odontodes are believed to share a common evolutionary history. Gene expression patterns have previously been compared between mammalian and teleost tooth development but we highlight how the comparative framework was not always properly defined to deal with different tooth types or tooth developmental stages. Larger-scale comparative analyses also included cartilaginous fishes: sharks display oral teeth and dermal scales for which the gene expression during development started to be investigated in the small-spotted catshark Scyliorhinus canicula during the past decade. We report several descriptive approaches to analyse the embryonic tooth and caudal scale gene expressions in S. canicula. We compare these expressions wih the ones reported in mouse molars and teleost oral and pharyngeal teeth and highlight contributions and biases that arise from these interspecific comparisons. We finally discuss the evolutionary processes that can explain the observed intra and interspecific similarities and divergences in the genetic cascades involved in tooth and odontode development in jawed vertebrates.

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

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

Microsatellite primer development in elasmobranchs using next generation sequencing of enriched libraries.

Microsatellites are useful in studies of population genetics, sibship, and parentage. Here, we screened for microsatellites from multiple elasmobranch genomic libraries using an enrichment protocol followed by sequencing on an Illumina platform. We concurrently screened five and then nine genomes and describe the number of potential loci from each respective round of sequencing. To validate the efficacy of the protocol, we developed and tested primers for the pelagic thresher shark, Alopias pelagicus. The method described here is a cost-effective protocol to increase the pool of potential useful loci and allows the concurrent screening of multiple libraries.

Genetic and historical evidence of common sawsharks Pristiophorus cirratus in the waters of southern Queensland.

In 2011, a male pristiophorid was caught by a prawn trawler north east of Cape Moreton, Queensland, Australia. Molecular analyses confirmed the specimen to be the common sawshark Pristiophorus cirratus. Historical catch data indicate the occurrence of the species in the region but this is the first verified record of P. cirratus occurring in the waters of southern Queensland. Together, these records extend the recognised northern limit of P. cirratus by c. 500 km, which suggests that further investigation of its distribution is warranted.

Electroreception in marine fishes: chondrichthyans.

Electroreception in marine fishes occurs across a variety of taxa and is best understood in the chondrichthyans (sharks, skates, rays, and chimaeras). Here, we present an up-to-date review of what is known about the biology of passive electroreception and we consider how electroreceptive fishes might respond to electric and magnetic stimuli in a changing marine environment. We briefly describe the history and discovery of electroreception in marine Chondrichthyes, the current understanding of the passive mode, the morphological adaptations of receptors across phylogeny and habitat, the physiological function of the peripheral and central nervous system components, and the behaviours mediated by electroreception. Additionally, whole genome sequencing, genetic screening and molecular studies promise to yield new insights into the evolution, distribution, and function of electroreceptors across different environments. This review complements that of electroreception in freshwater fishes in this special issue, which provides a comprehensive state of knowledge regarding the evolution of electroreception. We conclude that despite our improved understanding of passive electroreception, several outstanding gaps remain which limits our full comprehension of this sensory modality. Of particular concern is how electroreceptive fishes will respond and adapt to a marine environment that is being increasingly altered by anthropogenic electric and magnetic fields.

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.

Role of elasmobranchs and holocephalans in understanding peptide evolution in the vertebrates: Lessons learned from gonadotropin releasing hormone (GnRH) and corticotropin releasing factor (CRF) phylogenies.

The cartilaginous fishes (Class Chondrichthyes) comprise two morphologically distinct subclasses; Elasmobranchii and Holocephali. Evidence indicates early divergence of these subclasses, suggesting monophyly of their lineage. However, such a phylogenetic understanding is not yet developed within two highly conserved peptide lineages, GnRH and CRF. Various GnRH forms exist across the Chondrichthyes. Although 4-7 immunoreactive forms have been described in Elasmobranchii, only one has been elucidated in Holocephali. In contrast, Chondrichthyan CRF phylogeny follows a pattern more consistent with vertebrate evolution. For example, three forms are expressed within the lamprey, with similar peptides present within the genome of the Callorhinchus milii, a holocephalan. Although these findings are consistent with recent evidence regarding the phylogenetic age of Chondrichthyan lineages, CRF evolution in vertebrates remains elusive. Assuming that the Elasmobranchii and Holocephali are part of a monocladistic clade within the Chondrichthyes, we interpret the findings of GnRH and CRF to be products of their respective lineages.

Synthetic capacity does not predict elasmobranchs' ability to maintain trimethylamine oxide without a dietary contribution.

Trimethylamine oxide (TMAO) is an organic osmolyte and universal protein stabilizer. Its role as a cytoprotectant is particularly important in ureosmotic elasmobranchs that accumulate high levels of urea, a macromolecular perturbant. Feeding is a key component in the turnover and maintenance of these nitrogenous compounds. However, previous studies examining TMAO regulation have been largely completed using starved individuals, when nitrogen balance is altered. Here, under fed conditions, we test the importance of dietary TMAO on long-term maintenance in three elasmobranch species with differing endogenous synthetic capacities. Smoothhounds (Mustelus canis), spiny dogfish (Squalus acanthias), and little skates (Leucoraja erinacea) exhibited species- and tissue-specific differences in their ability to conserve TMAO when fed a low TMAO diet for 56days. Smoothhounds, a species with the capacity for endogenous production, exhibited a decrease in muscle TMAO. Spiny dogfish and little skates, species with no reported ability for synthesis, exhibited decreases in plasma and liver TMAO, respectively. Our findings are contrary to previous starvation studies demonstrating constant levels of TMAO for up to 56days in elasmobranchs. Further, the previously reported synthetic capacity of these species did not correlate with their ability to conserve TMAO and cannot be used to predict a species reliance on dietary contributions for prolonged maintenance. It is possible that all species rely to a degree on absorption of TMAO from the diet or that alternate synthetic or regulatory pathways play a larger role than previously thought.

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.

Distributional shift of urea production site from the extraembryonic yolk sac membrane to the embryonic liver during the development of cloudy catshark (Scyliorhinus torazame).

Urea is an essential osmolyte for marine cartilaginous fishes. Adult elasmobranchs and holocephalans are known to actively produce urea in the liver, muscle and other extrahepatic organs; however, osmoregulatory mechanisms in the developing cartilaginous fish embryo with an undeveloped urea-producing organ are poorly understood. We recently described the contribution of extraembryonic yolk sac membranes (YSM) to embryonic urea synthesis during the early developmental period of the oviparous holocephalan elephant fish (Callorhinchus milii). In the present study, to test whether urea production in the YSM is a general phenomenon among oviparous Chondrichthyes, we investigated gene expression and activities of ornithine urea cycle (OUC) enzymes together with urea concentrations in embryos of the elasmobranch cloudy catshark (Scyliorhinus torazame). The intracapsular fluid, in which the catshark embryo develops, had a similar osmolality to seawater, and embryos maintained a high concentration of urea at levels similar to that of adult plasma throughout development. Relative mRNA expressions and activities of catshark OUC enzymes were significantly higher in YSM than in embryos until stage 32. Concomitant with the development of the embryonic liver, the expression levels and activities of OUC enzymes were markedly increased in the embryo from stage 33, while those of the YSM decreased from stage 32. The present study provides further evidence that the YSM contributes to embryonic urea homeostasis until the liver and other extrahepatic organs become fully functional, and that urea-producing tissue shifts from the YSM to the embryonic liver in the late developmental period of oviparous marine cartilaginous fishes.

Molecular organization of the 5S rDNA gene type II in elasmobranchs.

The 5S rDNA gene is a non-coding RNA that can be found in 2 copies (type I and type II) in bony and cartilaginous fish. Previous studies have pointed out that type II gene is a paralog derived from type I. We analyzed the molecular organization of 5S rDNA type II in elasmobranchs. Although the structure of the 5S rDNA is supposed to be highly conserved, our results show that the secondary structure in this group possesses some variability and is different than the consensus secondary structure. One of these differences in Selachii is an internal loop at nucleotides 7 and 112. These mutations observed in the transcribed region suggest an independent origin of the gene among Batoids and Selachii. All promoters were highly conserved with the exception of BoxA, possibly due to its affinity to polymerase III. This latter enzyme recognizes a dT4 sequence as stop signal, however in Rajiformes this signal was doubled in length to dT8. This could be an adaptation toward a higher efficiency in the termination process. Our results suggest that there is no TATA box in elasmobranchs in the NTS region. We also provide some evidence suggesting that the complexity of the microsatellites present in the NTS region play an important role in the 5S rRNA gene since it is significantly correlated with the length of the NTS.

The diversity, distribution and status of deep-water elasmobranchs in the Rockall Trough, north-east Atlantic Ocean.

Data from a scientific deep-water trawl fisheries survey in the north-east Atlantic were analysed to determine the spatial and bathymetric distribution of elasmobranch species and assess the change in relative abundance over the period 1998-2013. During this period, commercial fisheries for deep-water sharks went from being entirely unregulated, to being briefly managed, to being completely prohibited. A total of 22 species of shark and 10 species of skate were recorded between depths of 300 and 2030 m. All showed strong species-specific depth-related trends in abundance. Out of the 11 more common species, five showed no change in relative abundance over time, two (Centrophorus squamosus and Centroselachus crepidater) declined significantly and four increased in relative abundance (Apristurus aphyodes, Apristurus microps, Galeus melastomus and Deania calcea). Assuming these populations were depleted by fisheries in the past, the current data do not suggest there has been an overall recovery. Positive signs for some species in the most recent years suggest movement or recruitment back into the area; however, it is of concern that two species continued to decline. There is a continued need to have precautionary management of these elasmobranch species, and the current ban on landing these species in European waters remains appropriate.

Molecular phylogeny of elasmobranchs inferred from mitochondrial and nuclear markers.

The elasmobranchs (sharks, rays and skates) being the extant survivors of one of the earliest offshoots of the vertebrate evolutionary tree are good model organisms to study the primitive vertebrate conditions. They play a significant role in maintaining the ecological balance and have high economic value. Due to over-exploitation and illegal fishing worldwide, the elasmobranch stocks are being decimated at an alarming rate. Appropriate management measures are necessary for restoring depleted elasmobranch stocks. One approach for restoring stocks is implementation of conservation measures and these measures can be formulated effectively by knowing the evolutionary relationship among the elasmobranchs. In this study, a total of 30 species were chosen for molecular phylogeny studies using mitochondrial cytochrome c oxidase subunit I, 12S ribosomal RNA gene and nuclear Internal Transcribed Spacer 2. Among different genes, the combined dataset of COI and 12S rRNA resulted in a well resolved tree topology with significant bootstrap/posterior probabilities values. The results supported the reciprocal monophyly of sharks and batoids. Within Galeomorphii, Heterodontiformes (bullhead sharks) formed as a sister group to Lamniformes (mackerel sharks): Orectolobiformes (carpet sharks) and to Carcharhiniformes (ground sharks). Within batoids, the Myliobatiformes formed a monophyly group while Pristiformes (sawfishes) and Rhinobatiformes (guitar fishes) formed a sister group to all other batoids.

Holocephalan embryos provide evidence for gill arch appendage reduction and opercular evolution in cartilaginous fishes.

Chondrichthyans possess endoskeletal appendages called branchial rays that extend laterally from their hyoid and gill-bearing (branchial) arches. Branchial ray outgrowth, like tetrapod limb outgrowth, is maintained by Sonic hedgehog (Shh) signaling. In limbs, distal endoskeletal elements fail to form in the absence of normal Shh signaling, whereas shortened duration of Shh expression correlates with distal endoskeletal reduction in naturally variable populations. Chondrichthyans also exhibit natural variation with respect to branchial ray distribution--elasmobranchs (sharks and batoids) possess a series of ray-supported septa on their hyoid and gill arches, whereas holocephalans (chimaeras) possess a single hyoid arch ray-supported operculum. Here we show that the elongate hyoid rays of the holocephalan Callorhinchus milii grow in association with sustained Shh expression within an opercular epithelial fold, whereas Shh is only transiently expressed in the gill arches. Coincident with this transient Shh expression, branchial ray outgrowth is initiated in C. milii but is not maintained, yielding previously unrecognized vestigial gill arch branchial rays. This is in contrast to the condition seen in sharks, where sustained Shh expression corresponds to the presence of fully formed branchial rays on the hyoid and gill arches. Considered in light of current hypotheses of chondrichthyan phylogeny, our data suggest that the holocephalan operculum evolved in concert with gill arch appendage reduction by attenuation of Shh-mediated branchial ray outgrowth, and that chondrichthyan branchial rays and tetrapod limbs exhibit parallel developmental mechanisms of evolutionary reduction.

Proteomic characterization of the acute-phase response of yellow stingrays Urobatis jamaicensis after injection with a Vibrio anguillarum-ordalii bacterin.

Systemic inflammatory responses of mammals and bony fish are primarily driven by coordinated up-regulation and down-regulation of plasma acute-phase proteins. Although this general principle is believed to be universal among vertebrates, it remains relatively unexplored in elasmobranchs. The objective of this study was to characterize acute changes in the plasma proteome of three yellow stingrays Urobatis jamaicensis following intraperitoneal injection with a commercial Vibrio bacterin. Changes in plasma protein levels were analyzed immediately prior to vaccination (time 0) and at 24 and 72 h post-injection by isobaric tags for relative and absolute quantitation (iTRAQ 4-plex) using shotgun-based nano liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and de novo peptide sequencing. Pooled 2D-LC-MS/MS and de novo sequencing data revealed differential expression of 156 distinct plasma proteins between time 0 and at least one post-vaccination time point. Using 1.5-fold change in expression as physiologically significant, 14/156 (9.0%) proteins were upregulated in at least one stingray through at least one experimental timepoint. Upregulated proteins included complement factors, Mx-protein, hemopexin, factor X and prothrombin. Seventy-six of 156 (48.7%) proteins were downregulated in the acute-phase response, including transferrin, apolipoprotein B, heparin cofactor 2, alpha2-macroglobulin, and various growth factors. Other differentially upregulated or downregulated proteins included intracellular, cell binding and structural proteins, proteins involved in physiologic processes, and unknown/hypothetical proteins. Selected bioactive factors are discussed for their putative roles in the elasmobranchs acute-phase response. These findings contribute to our understanding of disease processes in elasmobranchs, immunologic phylogeny in vertebrates, and begin the search for potential biomarkers of disease in these ecologically important fish.

Sonic hedgehog function in chondrichthyan fins and the evolution of appendage patterning.

The genetic mechanisms regulating tetrapod limb development are well characterized, but how they were assembled during evolution and their function in basal vertebrates is poorly understood. Initial studies report that chondrichthyans, the most primitive extant vertebrates with paired appendages, differ from ray-finned fish and tetrapods in having Sonic hedgehog (Shh)-independent patterning of the appendage skeleton. Here we demonstrate that chondrichthyans share patterns of appendage Shh expression, Shh appendage-specific regulatory DNA, and Shh function with ray-finned fish and tetrapods. These studies demonstrate that some aspects of Shh function are deeply conserved in vertebrate phylogeny, but also highlight how the evolution of Shh regulation may underlie major morphological changes during appendage evolution.

Mitochondrial haplotypes indicate parapatric-like phylogeographic structure in blue-spotted maskray (Neotrygon kuhlii) from the Coral Triangle region.

Phylogeographic structure was investigated in the blue-spotted maskray, Neotrygon kuhlii, focusing on the Coral Triangle region. We used as genetic marker a 519-bp fragment of the cytochrome c-oxidase subunit I (COI) gene, sequenced in a total of 147 individuals from 26 sampling locations. The parsimony network of COI haplotypes was split into seven distinct clades within the Coral Triangle region. Different clades had exclusive but contiguous geographic distributions, indicating parapatric-like phylogeographic structure. Strong genetic differences were also inferred between local populations within a clade, where reciprocal monophyly between geographically adjacent samples was observed on several instances. Nearly 25% of the total molecular variance could be ascribed to differences between geographic samples within a clade, whereas interclade variation accounted for >65% of the total variance. The strong phylogeographic structure observed within a clade can be explained by either sedentarity or female philopatry. We interpret the parapatric distribution of clades as the joint result of 1) expansion from refuge populations at times of low sea level, and 2) possible enhanced competition between individuals from different clades, or assortative mating, or hybrid zones, along lines of secondary contact. The parapatric-like structure uncovered in the present study parallels regional differences at nuclear marker loci, thus pointing to incipient speciation within Coral Triangle N. kuhlii.