Shark genome size evolution and its relationship with cellular, life-history, ecological, and diversity traits

Among vertebrates, sharks exhibit both large and heterogeneous genome sizes ranging from 2.86 to 17.05 pg. Aiming for a better understanding of the patterns and causalities of shark genome size evolution, we applied phylogenetic comparative methods to published genome-size estimates for 71 species representing the main phylogenetic lineages, life-histories and ecological traits. The sixfold range of genome size variation was strongly traceable throughout the phylogeny, with a major expansion preceding shark diversification during the late Paleozoic and an ancestral state (6.33 pg) close to the present-day average (6.72 pg). Subsequent deviations from this average occurred at higher rates in squalomorph than in galeomorph sharks and were unconnected to evolutionary changes in the karyotype architecture, which were dominated by descending disploidy events. Genome size was positively correlated with cell and nucleus sizes and negatively with metabolic rate. The metabolic constraints on increasing genome size also manifested at higher phenotypic scales, with large genomes associated with slow lifestyles and purely marine waters. Moreover, large genome sizes were also linked to non-placental reproductive modes, which may entail metabolically less demanding embryological developments. Contrary to ray-finned fishes, large genome size was associated neither with the taxonomic diversity of affected clades nor with low genetic diversity.

Genetic identification of three CITES-listed sharks using a paper-based Lab-on-a-Chip (LOC)

Threatened shark species are caught in large numbers by artisanal and commercial fisheries and traded globally. Monitoring both which shark species are caught and sold in fisheries, and the export of CITES-restricted products, are essential in reducing illegal fishing. Current methods for species identification rely on visual examination by experts or DNA barcoding techniques requiring specialist laboratory facilities and trained personnel. The need for specialist equipment and/or input from experts means many markets are currently not monitored. We have developed a paper-based Lab-on-a-Chip (LOC) to facilitate identification of three threatened and CITES-listed sharks, bigeye thresher (Alopias superciliosus), pelagic thresher (A. pelagicus) and shortfin mako shark (Isurus oxyrinchus) at market source. DNA was successfully extracted from shark meat and fin samples and combined with DNA amplification and visualisation using Loop Mediated Isothermal Amplification (LAMP) on the LOC. This resulted in the successful identification of the target species of sharks in under an hour, with a working positive and negative control. The LOC provided a simple "yes" or "no" result via a colour change from pink to yellow when one of the target species was present. The LOC serves as proof-of-concept (PoC) for field-based species identification as it does not require specialist facilities. It can be used by non-scientifically trained personnel, especially in areas where there are suspected high frequencies of mislabelling or for the identification of dried shark fins in seizures.

Molecular identification of sharks from the genus Sphyrna (Elasmobranchii: Chondrichthyes) in Maranhão Coast (Brazil)

Sharks of the genus Sphyrna are under intense exploitation globally. In Brazil's northern coast, this genus represents a high proportion of fisheries landings and comprises four species. However, due to difficulty of specific identification when specimens are landed, most of the records are limited to the genus level. Here we analyzed the effectiveness of ITS2 (Internal Transcribed Spacer 2 of rDNA) fragment length protocol (Abercrombie et al., 2005) for identifying hammerhead shark species, comparing with the analysis of COI (Cytochrome oxidase subunit I) and ITS2 sequences. We evaluated samples of muscle tissue acquired in the main fishing ports of Maranhão: Carutapera, Raposa e Tutóia. Sampling was conducted between March 2017 to March 2018 and complemented with material deposited in collection (2015). COI results indicated the occurrence of endangered species which are prohibited to be landed. These include Sphyrna mokarran (67%), S. lewini (15%), S. tudes (3%), and S. tiburo (15%). For the ITS2 marker, we investigated the optimization of the protocol developed by Abercrombie (2005) for to improve the use in this geographical area througout design of a new primers.

First molecular identification and phylogenetic illustration of Sarcocystis species infection in Red Sea shortfin mako shark (Isurus oxyrinchus Rafinesque, 1810)

BACKGROUND

members of the genus Sarcocystis are intracellular obligate protozoan parasites classified within the phylum Apicomplexa and have an obligate heteroxenous life cycle involving two hosts. A more comprehensive understanding of the prevalence and geographic range of different Sarcocystis species in marine ecosystems is needed globally and nationally. Hence, the objective of this study was to document the incidence of Sarcocystis infection in sharks within the aquarium ecosystem of Egypt and to identify the species through the characterization of the SSU rDNA gene.

METHODS

All organs of the mako shark specimen underwent macroscopic screening to detect the existence of a Sarcocystis cyst. Ten cysts were collected from the intestine and processed separately to extract the genomic DNA. The polymerase chain reaction (PCR) was accomplished by amplifying a specific 18S ribosomal RNA (rRNA) gene fragment. Subsequently, the resulting amplicons were subjected to purification and sequencing processes.

RESULTS

Macroscopic examination of the mako shark intestinal wall sample revealed the presence of Sarcocystis cysts of various sizes and shapes, and sequencing of the amplicons from Sarcocystis DNA revealed a 100% nucleotide identity with the sequence of Sarcocystis tenella recorded from sheep in Iran; The mako shark sequence has been deposited in the GeneBank with the accession number OQ721979. This study presents the first scientific evidence demonstrating the presence of the Sarcocystis parasite in sharks, thereby documenting this specific marine species as a novel intermediate host in the Sarcocystis life cycle.

CONCLUSIONS

This is the first identification of Sarcocystis infection in sharks, and we anticipate it will be an essential study for future screenings and establishing effective management measures for this disease in aquatic ecosystems.

De novo transcriptome for Chiloscyllium griseum, a long-tail carpet shark of the Indian waters

Sharks have thrived in the oceans for 400 million years, experienced five extinctions and evolved into today's apex predators. However, enormous genome size, poor karyotyping and limited tissue sampling options are the bottlenecks in shark research. Sharks of the family Orectolobiformes act as model species in transcriptome research with exceptionally high reproductive fecundity, catch prominence and oviparity. The present study illustrates a de novo transcriptome for an adult grey bamboo shark, Chiloscyllium griseum (Chondrichthyes; Hemiscyllidae) using paired-end RNA sequencing. Around 150 million short Illumina reads were obtained from five different tissues and assembled using the Trinity assembler. 70,647 hits on Uniprot by BLASTX was obtained after the transcriptome annotation. The data generated serve as a basis for transcriptome-based population genetic studies and open up new avenues in the field of comparative transcriptomics and conservation biology.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Diversity and genetic structure of freshwater shark Wallago attu: an emerging species of commercial interest

Pakistan has natural freshwater resources acting as a hotspot for diverse fish fauna. However, this aquatic fauna is declining at an alarming rate due to over-exploitation, habitat degradation, water pollution, climate change, and certain anthropogenic activities. The freshwater shark, Wallago attu, is a popular edible catfish inhabiting these freshwater ecosystems. Habitat degradation, overfishing, and human activities are heavily impacting the natural population of this species. So, sound knowledge about its population structure is necessary for its proper management in natural waters. The current study involves utilizing two mtDNA markers (COI, Cytb) to assess the genetic structure and differentiation among W. attu populations of Pakistani Rivers. Genetic variability analysis indicated a high haplotype (0.343 ± 0.046-0.870 ± 0.023) and low nucleotide diversity (0.0024 ± 0.012-0.0038 ± 0.018) among single and combined gene sequences, respectively. Overall, River Indus was populated with more diverse fauna of Wallago attu as compared to River Chenab and River Ravi. Population pairwise, Fst values (0.40-0.61) were found to be significantly different (p < 0.01) among three Riverine populations based upon combined gene sequences. The gene flow for the combined gene (COI + Cytb) dataset among three populations was less than 1.0. The transition/transversion bias value R (0.58) was calculated for testing of neutral evolution, and it declared low genetic polymorphism among natural riverine populations of Wallago attu. The current study's findings would be meaningful in planning the management and conservation of this economically important catfish in future.

DNA barcoding continues to identify endangered species of shark sold as food in a globally significant shark fin trade hub

Shark fins are a delicacy consumed throughout Southeast Asia. The life history characteristics of sharks and the challenges associated with regulating fisheries and the fin trade make sharks particularly susceptible to overfishing. Here, we used DNA barcoding techniques to investigate the composition of the shark fin trade in Singapore, a globally significant trade hub. We collected 505 shark fin samples from 25 different local seafood and Traditional Chinese Medicine shops. From this, we identified 27 species of shark, three species are listed as Critically Endangered, four as Endangered and ten as Vulnerable by the International Union for Conservation of Nature (IUCN). Six species are listed on CITES Appendix II, meaning that trade must be controlled in order to avoid utilization incompatible with their survival. All dried fins collected in this study were sold under the generic term "shark fin"; this vague labelling prevents accurate monitoring of the species involved in the trade, the effective implementation of policy and conservation strategy, and could unwittingly expose consumers to unsafe concentrations of toxic metals. The top five most frequently encountered species in this study are Rhizoprionodon acutus, Carcharhinus falciformis, Galeorhinus galeus, Sphyrna lewini and Sphyrna zygaena. Accurate labelling that indicates the species of shark that a fin came from, along with details of where it was caught, allows consumers to make an informed choice on the products they are consuming. Doing this could facilitate the avoidance of species that are endangered, and similarly the consumer can choose not to purchase species that are documented to contain elevated concentrations of toxic metals.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Northerly range expansion and first confirmed records of the smalltooth sand tiger shark, Odontaspis ferox, in the United Kingdom and Ireland

Three Odontaspis ferox (confirmed by mtDNA barcoding) were found in the English Channel and Celtic Sea in 2023 at Lepe, UK (50.7846, -1.3508), Kilmore Quay, Ireland (52.1714, -6.5937), and Lyme Bay, UK (50.6448, -2.9302). These are the first records of O. ferox in either country, and extend the species' range by over three degrees of latitude, to >52° N. They were ~275 (female), 433 (female), and 293 cm (male) total length, respectively. These continue a series of new records, possibly indicative of a climate change-induced shift in the species' range.

DNA metabarcoding of cloacal swabs provides insight into diets of highly migratory sharks in the Mid-Atlantic Bight

The abundances of migratory shark species observed throughout the Mid-Atlantic Bight (MAB) during productive summer months suggest that this region provides critical habitat and prey resources to these taxa. However, the principal prey assemblages sustaining migratory shark biomass in this region are poorly defined. We applied high-throughput DNA metabarcoding to shark feces derived from cloacal swabs across nine species of Carcharhinid and Lamnid sharks to (1) quantify the contribution of broad taxa (e.g., invertebrates, fishes) supporting shark biomass during seasonal residency in the MAB and (2) determine whether the species displayed distinct dietary preference indicative of resource partitioning. DNA metabarcoding resulted in high taxonomic (species-level) resolution of shark diets with actinopterygian and elasmobranch fishes as the dominant prey categories across the species. DNA metabarcoding identified several key prey groups consistent across shark taxa that are likely integral for sustaining their biomass in this region, including Atlantic menhaden (Brevoortia tyrannus), Atlantic mackerel (Scomber scombrus), and benthic elasmobranchs, including skates. Our results are consistent with previously published stomach content data for the shark species of similar size range in the Northwest Atlantic Ocean, supporting the efficacy of cloacal swab DNA metabarcoding as a minimally invasive diet reconstruction technique. The high reliance of several shark species on Atlantic menhaden could imply wasp-waist food-web conditions during the summer months, whereby high abundances of forage fishes sustain a diverse suite of migratory sharks within a complex, seasonal food web.

Endangered shark species traded as "cação" in São Paulo during the COVID-19 lockdown: DNA-barcoding a snapshot of products

BACKGROUND

Elasmobranch populations are declining, predominantly driven by overfishing, and over a third of global sharks, rays, and chimeras are estimated to be threatened with extinction. In terms of trade, Brazil is ranked the eleventh-largest shark producer and the top importer of shark meat in the world. Research has shown that elasmobranchs are sold in Brazil under the name "cação" (a generic designation for cartilaginous fish) to overcome consumer resistance.

METHODOLOGY AND RESULTS

This study used DNA barcoding to investigate the sale of sharks in the State of São Paulo during the COVID-19 lockdown. A total of 35 samples of "cação" were analysed, revealing six different shark species on sale, including Carcharhinus falciformis, Carcharhinus signatus, Carcharias taurus, Isurus oxyrinchus, and Isurus paucus, that are threatened with extinction according to the IUCN red list. This study demonstrates that vulnerable elasmobranchs are being commercialised under the label "cação" in the São Paulo State and Brazil.

CONCLUSIONS

Comparison of shark products traded before and during the COVID-19 pandemic showed no significant difference, suggesting lockdown did not affect patterns of species commercialisation. Effective fisheries and sale monitoring, correct product labelling legislation and increased consumer awareness that "cação" is shark are needed for appropriate conservation and management of shark populations in Brazil.

A singular shark bitter taste receptor provides insights into the evolution of bitter taste perception

Many animal and plant species synthesize toxic compounds as deterrent; thus, detection of these compounds is of vital importance to avoid their ingestion. Often, such compounds are recognized by taste 2 receptors that mediate bitter taste in humans. Until now, bitter taste receptors have only been found in bony vertebrates, where they occur as a large family already in coelacanth, a "living fossil" and the earliest-diverging extant lobe-finned fish. Here, we have revisited the evolutionary origin of taste 2 receptors (T2Rs) making use of a multitude of recently available cartilaginous fish genomes. We have identified a singular T2R in 12 cartilaginous fish species (9 sharks, 1 sawfish, and 2 skates), which represents a sister clade to all bony fish T2Rs. We have examined its ligands for two shark species, a catshark and a bamboo shark. The ligand repertoire of bamboo shark represents a subset of that of the catshark, with roughly similar thresholds. Amarogentin, one of the most bitter natural substances for humans, also elicited the highest signal amplitudes with both shark receptors. Other subsets of ligands are shared with basal bony fish T2Rs indicating an astonishing degree of functional conservation over nearly 500 mya of separate evolution. Both shark receptors respond to endogenous steroids as well as xenobiotic compounds, whereas separate receptors exist for xenobiotics both in early- and late-derived bony vertebrates (coelacanth, zebrafish, and human), consistent with the shark T2R reflecting the original ligand repertoire of the ancestral bitter taste receptor at the evolutionary origin of this family.

Evolution of reproductive modes in sharks and rays

The ecological and life history drivers of the diversification of reproductive modes in early vertebrates are not fully understood. Sharks, rays and chimaeras (group Chondrichthyes) have an unusually diverse variety of reproductive modes and are thus an ideal group to test the factors driving the evolution of reproductive complexity. Here, using 960 species representing all major Chondrichthyes taxa, we reconstruct the evolution of their reproduction modes and investigate the ecological and life history predictors of reproduction. We show that the ancestral Chondrichthyes state was egg-laying and find multiple independent transitions between egg-laying and live-bearing via an intermediate state of yolk-only live-bearing. Using phylogenetically informed analysis, we also show that live-bearing species have larger body size and larger offspring than egg-laying species. In addition, live-bearing species are distributed over shallow to intermediate depths, while egg-layers are typically found in deeper waters. This suggests that live-bearing is more closely associated with pelagic, rather than demersal habitats. Taken together, using a basal vertebrate group as a model, we demonstrat how reproductive mode co-evolves with environmental conditions and life-history traits.

Low mutation rate in epaulette sharks is consistent with a slow rate of evolution in sharks

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

Genetic stock structure of the silky shark Carcharhinus falciformis in the Indo-Pacific Ocean

The silky shark, Carcharhinus falciformis, is a cosmopolitan species commonly caught as a bycatch for longline fisheries. However, the genetic stock structure for the Indo-Pacific Ocean is not well-defined yet. Here, we used eight microsatellite loci to examine the genetic stock structure and effective population size of 307 silky sharks across 5 Indo-Pacific sampling locations. A major genetic break was found between Aceh and the remaining locations (FST = 0.0505-0.0828, p = 0.001). The Indian Ocean displayed a slightly lower effective population estimate (Ne) compared to the Pacific Ocean, potentially due to the higher fishing pressure in the Indian Ocean region. The lowest Ne was found in the Aceh population (Ne = 2.3), suggesting it might be a small and endemic population. These findings offer valuable information for the conservation and management of the silky shark. We suggest that the population around Aceh waters constitutes a distinct stock and should be managed independently. Further investigations into migratory and movement patterns are needed to define the boundaries of different stocks, ensuring effective management the silky shark across the Indo-Pacific region.

An Ancient MHC-Linked Gene Encodes a Nonrearranging Shark Antibody, UrIg, Convergent with IgG

Gnathostome adaptive immunity is defined by the Ag receptors, Igs and TCRs, and the MHC. Cartilaginous fish are the oldest vertebrates with these adaptive hallmarks. We and others have unearthed nonrearranging Ag receptor-like genes in several vertebrates, some of which are encoded in the MHC or in MHC paralogous regions. One of these genes, named UrIg, was detected in the class III region of the shark MHC that encodes a protein with typical V and C domains such as those found in conventional Igs and TCRs. As no transmembrane region was detected in gene models or cDNAs, the protein does not appear to act as a receptor. Unlike some other shark Ig genes, the UrIg V region shows no evidence of RAG-mediated rearrangement, and thus it is likely related to other V genes that predated the invasion of the RAG transposon. The UrIg gene is present in all elasmobranchs and evolves conservatively, unlike Igs and TCRs. Also, unlike Ig/TCR, the gene is not expressed in secondary lymphoid tissues, but mainly in the liver. Recombinant forms of the molecule form disulfide-linked homodimers, which is the form also detected in many shark tissues by Western blotting. mAbs specific for UrIg identify the protein in the extracellular matrix of several shark tissues by immunohistochemistry. We propose that UrIg is related to the V gene invaded by the RAG transposon, consistent with the speculation of emergence of Ig/TCR within the MHC or proto-MHC.

Genetic evidence for plastic reproductive philopatry and matrotrophy in blacktip reef sharks (Carcharhinus melanopterus) of the Moorea Island (French Polynesia)

The exploitation of sharks and the degradation of their habitats elevate the urgency to understand the factors that influence offspring survival and ultimately shark reproductive success. We monitored and sampled blacktip reef sharks (Carcharhinus melanopterus) in nursery habitats of Moorea Island (French Polynesia), to improve knowledge on shark reproductive behavior and biology. We sampled fin clips and morphometrics from 230 young-of-the-year sharks and used microsatellite DNA markers to process parentage analysis to study the reproductive philopatric behavior in female sharks and the matrotrophy within litters. These traits are driving the success of the local replenishment influencing selection through birth site and maternal reserves transmitted to pups. Parentage analysis revealed that some female sharks changed their parturition areas (inter-seasonally) while other female sharks came back to the same site for parturition, providing evidence for a plastic philopatric behavior. Morphometrics showed that there was no significant relationship between body condition indices and nursery locations. However, similarities and differences in body condition were observed between individuals sharing the same mother, indicating that resource allocation within some shark litters might be unbalanced. Our findings further our understanding of the reproductive biology and behavior that shape shark populations with the aim to introduce these parameters into future conservation strategies.

Philopatry influences the genetic population structure of the blacktip shark (Carcharhinus limbatus) at multiple spatial scales

Understanding how interactions among microevolutionary forces generate genetic population structure of exploited species is vital to the implementation of management policies that facilitate persistence. Philopatry displayed by many coastal shark species can impact gene flow and facilitate selection, and has direct implications for the spatial scales of management. Here, genetic structure of the blacktip shark (Carcharhinus limbatus) was examined using a mixed-marker approach employing mitochondrial control region sequences and 4339 SNP-containing loci generated using ddRAD-Seq. Genetic variation was assessed among young-of-the-year sampled in 11 sites in waters of the United States in the western North Atlantic Ocean, including the Gulf of Mexico. Spatial and environmental analyses detected 68 nuclear loci putatively under selection, enabling separate assessments of neutral and adaptive genetic structure. Both mitochondrial and neutral SNP data indicated three genetically distinct units-the Atlantic, eastern Gulf, and western Gulf-that align with regional stocks and suggest regional philopatry by males and females. Heterogeneity at loci putatively under selection, associated with temperature and salinity, was observed among sites within Gulf units, suggesting local adaptation. Furthermore, five pairs of siblings were identified in the same site across timescales corresponding with female reproductive cycles. This indicates that females re-used a site for parturition, which has the potential to facilitate the sorting of adaptive variation among neighbouring sites. The results demonstrate differential impacts of microevolutionary forces at varying spatial scales and highlight the importance of conserving essential habitats to maintain sources of adaptive variation that may buffer species against environmental change.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Characteristics of the spiny dogfish (Squalus acanthias) nuclear genome

Sequenced shark nuclear genomes are underrepresented, with reference genomes available for only four out of nine orders so far. Here, we present the nuclear genome, with annotations, of the spiny dogfish (Squalus acanthias), a shark of interest to biomedical and conservation efforts, and the first representative of the second largest order of sharks (Squaliformes) with nuclear genome annotations available. Using Pacific Biosciences Continuous Long Read data in combination with Illumina paired-end and Hi-C sequencing, we assembled the genome de novo, followed by RNA-Seq-supported annotation. The final chromosome-level assembly is 3.7 Gb in size, has a BUSCO completeness score of 91.6%, and an error rate of less than 0.02%. Annotation predicted 33,283 gene models in the spiny dogfish's genome, of which 31,979 are functionally annotated.

A shark-inspired general model of tooth morphogenesis unveils developmental asymmetries in phenotype transitions

Developmental complexity stemming from the dynamic interplay between genetic and biomechanic factors canalizes the ways genotypes and phenotypes can change in evolution. As a paradigmatic system, we explore how changes in developmental factors generate typical tooth shape transitions. Since tooth development has mainly been researched in mammals, we contribute to a more general understanding by studying the development of tooth diversity in sharks. To this end, we build a general, but realistic, mathematical model of odontogenesis. We show that it reproduces key shark-specific features of tooth development as well as real tooth shape variation in small-spotted catsharks Scyliorhinus canicula. We validate our model by comparison with experiments in vivo. Strikingly, we observe that developmental transitions between tooth shapes tend to be highly degenerate, even for complex phenotypes. We also discover that the sets of developmental parameters involved in tooth shape transitions tend to depend asymmetrically on the direction of that transition. Together, our findings provide a valuable base for furthering our understanding of how developmental changes can lead to both adaptive phenotypic change and trait convergence in complex, phenotypically highly diverse, structures.

Ancient and nonuniform loss of olfactory receptor expression renders the shark nose a de facto vomeronasal organ

Cartilaginous fishes are renowned for a keen sense of smell, a reputation based on behavioral observations and supported by the presence of large and morphologically complex olfactory organs. At the molecular level, genes belonging to the four families coding for most olfactory chemosensory receptors in other vertebrates have been identified in a chimera and a shark, but it was unknown whether they actually code for olfactory receptors in these species. Here we describe the evolutionary dynamics of these gene families in cartilaginous fishes using genomes of a chimera, a skate, a sawfish and eight sharks. The number of putative OR, TAAR and V1R/ORA receptors is very low and stable whereas the number of putative V2R/OlfC receptors is higher and much more dynamic. In the catshark Scyliorhinus canicula, we show that many V2R/OlfC receptors are expressed in the olfactory epithelium in the sparsely distributed pattern characteristic for olfactory receptors. In contrast, the other three vertebrate olfactory receptor families are either not expressed (OR) or only represented with a single receptor (V1R/ORA and TAAR). The complete overlap of markers of microvillous olfactory sensory neurons with panneuronal marker HuC in the olfactory organ suggests the same cell type specificity of V2R/OlfC expression as for bony fishes, i.e. in microvillous neurons. The relatively low number of olfactory receptors in cartilaginous fishes compared to bony fishes could be the result of an ancient and constant selection in favor of a high olfactory sensitivity at the expense of a high discrimination capability.

DNA barcoding of fresh seafood in Australian markets reveals misleading labelling and sale of endangered species

Flake and shark samples were purchased from outlets in several coastal Australian regions and genetically barcoded using the cytochrome oxidase subunit 1 (CO1) gene to investigate labelling reliability and species-specific sources of ambiguously labelled fillets. Of the 41 shark fillet samples obtained, 23 yielded high-quality CO1 sequences, out of which 57% (n = 13) were labelled ambiguously (misleading) and 35% (n = 8) incorrectly. In contrast, barramundi fillets, which are widely available and sought after in Australian markets, were shown to be accurately labelled. Species identified from shark samples, including the shortfin mako (n = 3) and the scalloped hammerhead (n = 1), are assessed by the IUCN as endangered and critically endangered, respectively, with several others classified as vulnerable and near threatened.

IgNAR characterization and gene loci identification in whitespotted bamboo shark (Chiloscyllium plagiosum) genome

Single domain antibodies (sdAb) are promising candidates in cancer and anti-virus biotherapies for their unique structure characters. Though VHH and IgNAR have been discovered in camelidae and nurse shark (Ginlymostoma cirratum) respectively serval decades ago, expense of these large animals still limits the studies and applications of sdAb. Recently, IgNAR has been found in whitespotted bamboo shark (Chiloscyllium plagiosum), a small-sized sharks, while how to characterize and achieved the IgNAR of whitespotted bamboo shark is still unclear. In our research, we identified four IgNAR coding gene loci in whitespotted bamboo shark chromosome 44 (NC_057753.1), and primers were designed for single domain variable regions of IgNAR (VNAR) libraries preparation. Following sequencing results revealed that all plasmids constructed with our predicted VNAR libraries contained VNAR coding sequences, which confirmed the specificities of our primers in VNAR amplification. To our surprise, ≥90% VNAR sequences were encoded by IgNAR1, which suggests IgNAR1 is the most active IgNAR transcription locus in whitespotted bamboo shark. Interestingly, we found IgNAR(ΔC2-C3) were encoded by IgNAR3. Our findings gave a new sight of whitespotted bamboo shark IgNAR, which would broad the way of IgNAR studies and applications in biotherapies.

First evidence of multiple paternity in the blue shark (Prionace glauca)

Multiple paternity (MP) is defined as the behaviour in which females successfully mate with multiple males leading to offspring from different sires within the same litter. MP seems to be frequent and an evolutionary advantage in elasmobranchs. Here the authors report for the first time the occurrence of MP in the cosmopolitan blue shark Prionace glauca L. The evidence, gathered via microsatellite genotyping of pregnant females and their embryos, suggests that MP is very frequent in this species. Knowledge of MP in P. glauca should help describe more precisely its reproductive biology and contribute to the management of its populations.

Parthenogenesis in an elasmobranch in the presence of conspecific males

Parthenogenesis has been observed in several elasmobranch species, primarily in public aquaria. The majority of cases of parthenogenesis have occurred either when females were held without males or once a male was removed from a female's habitat. Here we report a second instance of parthenogenesis in a zebra shark female that was housed with conspecific mature males. This study calls into question the conditions under which elasmobranch females undergo parthenogenesis.

A priority nursery area for the conservation of the scalloped hammerhead shark Sphyrna lewini in Mexico

This study describes the genetic diversity and level of genetic differentiation of the scalloped hammerhead shark Sphyrna lewini from eight putative nursery areas in the Mexican Pacific Ocean, using the mtDNA control region. Genetic analyses revealed a particular spatial divergence between La Reforma and all the remaining sites, with five exclusive haplotypes and the highest genetic diversity. This pattern may be interpreted as the signature of regional female philopatry, relatively to a particular female-mediated gene flow for La Reforma, which shows a strong subdivision in the Gulf of California.

Phylogeny explains capture mortality of sharks and rays in pelagic longline fisheries: a global meta-analytic synthesis

Apex and mesopredators such as elasmobranchs are important for maintaining ocean health and are the focus of conservation efforts to mitigate exposure to fishing and other anthropogenic hazards. Quantifying fishing mortality components such as at-vessel mortality (AVM) is necessary for effective bycatch management. We assembled a database for 61 elasmobranch species and conducted a global meta-synthesis to estimate pelagic longline AVM rates. Evolutionary history was a significant predictor of AVM, accounting for up to 13% of variance in Bayesian phylogenetic meta-regression models for Lamniformes and Carcharhiniformes clades. Phylogenetically related species may have a high degree of shared traits that explain AVM. Model-estimated posterior mean AVM rates ranged from 5% (95% HDI 0.1%-16%) for pelagic stingrays and 76% (95% HDI 49%-90%) for salmon sharks. Measures that reduce catch, and hence AVM levels, such as input controls, bycatch quotas and gear technology to increase selectivity are appropriate for species with higher AVM rates. In addition to reducing catchability, handling-and-release practices and interventions such as retention bans in shark sanctuaries and bans on shark finning and trade hold promise for species with lower AVM rates. Robust, and where applicable, phylogenetically-adjusted elasmobranch AVM rates are essential for evidence-informed bycatch policy.

[Identification and functional characterization of circular RNAs in the liver of whitespotted bamboo shark (Chiloscyllium plagiosum)]

This study aims to identify the circular RNAs (circRNAs) in the liver of whitespotted bamboo shark (Chiloscyllium plagiosum) and to explore the effect of the overexpression of circRNAs on the proliferation and migration of hepatocellular carcinoma HepG2 cells. We conducted high-throughput sequencing for prediction of the circRNAs and then designed forward and reverse primers to verify them. Further, we constructed overexpression vectors for transient transfection of circRNAs into HepG2 cells. Finally, we employed CCK-8 assay and scratch assay to measure the proliferation and migration of the treated HepG2 cells. A total of 4 558 circRNAs were predicted, among which 14 circRNAs were confirmed. The qRT-PCR showed that circRNA 13-566, circRNA 4-475, circRNA 5-402, circRNA 294-177, and circRNA 30-219 were transiently overexpressed in HepG2 cells. The overexpression of these five circRNAs inhibited the proliferation and migration of HepG2 cells to varying degrees, and circRNA 4-475 and circRNA 294-177 had especially notable effect. This study provided a basic database of circRNA genes that particularly active in whitespotted bamboo shark liver and demonstrated with functional studies of these circRNAs potentially involved in normal and malignant liver cells.

Exploration of the Nurse Shark (Ginglymostoma cirratum) Plasma Immunoproteome Using High-Resolution LC-MS/MS

Many animals of scientific importance lack species-specific reagents (e.g., monoclonal antibodies) for in-depth studies of immune proteins. Mass spectrometry (MS)-based proteomics has emerged as a useful method for monitoring changes in protein abundance and modifications in non-model species. It can be used to quantify hundreds of candidate immune molecules simultaneously without the generation of new reagents. Here, we used MS-based proteomics to identify and quantify candidate immune proteins in the plasma of the nurse shark (Ginglymostoma cirratum), a cartilaginous fish and representative of the most basal extant vertebrate lineage with an immunoglobulin-based immune system. Mass spectrometry-based LC-MS/MS was performed on the blood plasma of nurse sharks immunized with human serum albumin (n=4) or sham immunized (n=1), and sampled at days 0 (baseline control), 1, 2, 3, 5, 7, 14, 21, 28, 25, 42 and 49. An antigen-specific antibody response was experimentally confirmed post-immunization. To provide a high-quality reference to identify proteins, we assembled and annotated a multi-tissue de novo transcriptome integrating long- and short-read sequence data. This comprised 62,682 contigs containing open reading frames (ORFs) with a length >80 amino acids. Using this transcriptome, we reliably identified 626 plasma proteins which were broadly categorized into coagulation, immune, and metabolic functional groups. To assess the feasibility of performing LC-MS/MS proteomics in nurse shark in the absence of species-specific protein annotations, we compared the results to an alternative strategy, mapping peptides to proteins predicted in the genome assembly of a related species, the whale shark (Rhincodon typus). This approach reliably identified 297 proteins, indicating that useful data on the plasma proteome may be obtained in many instances despite the absence of a species-specific reference protein database. Among the plasma proteins defined against the nurse shark transcriptome, fifteen showed consistent changes in abundance across the immunized shark individuals, indicating a role in the immune response. These included alpha-2-macroglobulin (A2M) and a novel protein yet to be characterized in diverse vertebrate lineages. Overall, this study enhances genetic and protein-level resources for nurse shark research and vastly improves our understanding of the elasmobranch plasma proteome, including its remodelling following immune stimulation.

Regional variation in multiple paternity in the brown smooth-hound shark Mustelus henlei from the northeastern Pacific

Multiple paternity (MP) in the brown smooth-hound shark (Mustelus henlei) was assessed in 15 litters (15 mothers and 97 embryos) collected in the northern Gulf of California of which 86.7% were sired by more than one male (i.e., from 2 to 4 sires). When taken together with results from previous studies, this record indicates that there is regional variation in MP in M. henlei in the northeastern Pacific. This pattern is associated with variations in the reproductive traits of each population (e.g., female size and litter size). In the Gulf of California, the results of a generalized linear model (GLZ) indicated that the litters of larger females had a higher probability of MP compared to those of smaller females.

Novel duplication remnant in the first complete mitogenome of Hemitriakis japanica and the unique phylogenetic position of family Triakidae

In this study, we firstly determined the complete mitogenome of the Japanese topeshark (Hemitriakis japonica), which belong to the family Triakidae and was assessed as Endangered A2d on the IUCN Red List in 2021. The mitogenome is 17,301 bp long, has a high AT content (60.0%), and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, a control region and specially a 594 bp-long non-coding region between Cytb gene and tRNA-Thr gene. The novel non-coding region share high sequence similarity with segments of the former and latter genes, so it was recognized as a duplication remnant. In addition, the Cytb gene and tRNA-Thr gene tandemly duplicated twice while accompanied by being deleted once at least. This is the first report of mitogenomic gene-arrangement in Triakidae. The phylogenetic trees were constructed using Bayesian inference (BI) and maximum likelihood (ML) methods based on the mitogenomic data of 51 shark species and two outgroups. In summary, basing on a novel type of gene rearrangements in houndshark mitogenome, the possibly rearranged process was analyzed and contributed further insight of shark mitogenomes evolution and phylogeny.

Site fidelity and shallow genetic structure in the common smooth-hound shark Mustelus mustelus confirmed by tag-recapture and genetic data

The common smooth-hound shark, Mustelus mustelus, is a widely distributed demersal shark under heavy exploitation from various fisheries throughout its distribution range. To assist in the development of appropriate management strategies, the authors evaluate stock structure, site fidelity and movement patterns along the species' distribution in southern Africa based on a combination of molecular and long-term tag-recapture data. Eight species-specific microsatellite markers (N = 73) and two mitochondrial genes, nicotinamide adenine dehydrogenase subunit 4 and control region (N = 45), did not reveal any significant genetic structure among neighbouring sites. Nonetheless, tagging data demonstrate a remarkable degree of site fidelity with 76% of sharks recaptured within 50 km of the original tagging location. On a larger geographic scale, dispersal is governed by oceanographic features as demonstrated by the lack of movements across the Benguela-Agulhas transition zone separating the South-East Atlantic Ocean (SEAO) and South-West Indian Ocean (SWIO) populations. Microsatellite data supported very shallow ocean-based structure (SEAO and SWIO) and historical southward gene flow following the Agulhas Current, corroborating the influence of this dynamic oceanographic system on gene flow. Moreover, no movements between Namibia and South Africa were observed, indicating that the Lüderitz upwelling formation off the Namibian coast acts as another barrier to dispersal and gene flow. Overall, these results show that dispersal and stock structure of M. mustelus are governed by a combination of behavioural traits and oceanographic features such as steep temperature gradients, currents and upwelling systems.

Preparation of Immune and Synthetic VNAR Libraries as Sources of High-Affinity Binders

The shark-derived autonomous variable antibody domains known as VNARs are attractive tools for therapeutic and diagnostic applications due to their favorable properties like small size (approximately 12 kDa), high thermal and chemical stability, and good tissue penetration. Currently, different techniques have been reported to generate VNAR domains against targets of therapeutic interest. Here, we describe methods for the preparation of an immune VNAR library based on bacteriophage display, and for the preparation of a synthetic library of VNAR domains using a modified protocol based on Kunkel mutagenesis. Finally, we describe procedures for in silico maturation of a VNAR using a bioinformatic approach to obtain higher affinity binders.

Phylogeography and population genetics of the cryptic bonnethead shark Sphyrna aff. tiburo in Brazil and the Caribbean inferred from mtDNA markers

Resolving the identity, phylogeny and distribution of cryptic species within species complexes is an essential precursor to management. The bonnethead shark, Sphyrna tiburo, is a small coastal shark distributed in the Western Atlantic from North Carolina (U.S.A.) to southern Brazil. Genetic analyses based on mitochondrial markers revealed that bonnethead sharks comprise a species complex with at least two lineages in the Northwestern Atlantic and the Caribbean (S. tiburo and Sphyrna aff. tiburo, respectively). The phylogeographic and phylogenetic analysis of two mitochondrial markers [control region (mtCR) and cytochrome oxidase I (COI)] showed that bonnethead sharks from southeastern Brazil correspond to S. aff. tiburo, extending the distribution of this cryptic species >5000 km. Bonnethead shark populations are only managed in the U.S.A. and in the 2000s were considered to be regionally extinct or collapsed in southeast Brazil. The results indicate that there is significant genetic differentiation between S. aff. tiburo from Brazil and other populations from the Caribbean (Φ_ST  = 0.9053, P < 0.000), which means that collapsed populations in the former are unlikely to be replenished from Caribbean immigration. The species identity of bonnethead sharks in the Southwest Atlantic and their relationship to North Atlantic and Caribbean populations still remains unresolved. Taxonomic revision and further sampling are required to reevaluate the status of the bonnethead shark complex through its distribution range.

Seeing through sedimented waters: environmental DNA reduces the phantom diversity of sharks and rays in turbid marine habitats

BACKGROUND

Sharks and rays are some of the most threatened marine taxa due to the high levels of bycatch and significant demand for meat and fin-related products in many Asian communities. At least 25% of shark and ray species are considered to be threatened with extinction. In particular, the density of reef sharks in the Pacific has declined to 3-10% of pre-human levels. Elasmobranchs are thought to be sparse in highly urbanised and turbid environments. Low visibility coupled with the highly elusive behaviour of sharks and rays pose a challenge to diversity estimation and biomonitoring efforts as sightings are limited to chance encounters or from carcasses ensnared in nets. Here we utilised an eDNA metabarcoding approach to enhance the precision of elasmobranch diversity estimates in urbanised marine environments.

RESULTS

We applied eDNA metabarcoding on seawater samples to detect elasmobranch species in the hyper-urbanised waters off Singapore. Two genes-vertebrate 12S and elasmobranch COI-were targeted and amplicons subjected to Illumina high-throughput sequencing. With a total of 84 water samples collected from nine localities, we found 47 shark and ray molecular operational taxonomic units, of which 16 had species-level identities. When data were compared against historical collections and contemporary sightings, eDNA detected 14 locally known species as well as two potential new records.

CONCLUSIONS

Local elasmobranch richness uncovered by eDNA is greater than the seven species sighted over the last two decades, thereby reducing phantom diversity. Our findings demonstrate that eDNA metabarcoding is effective in detecting shark and ray species despite the challenges posed by the physical environment, granting a more consistent approach to monitor these highly elusive and threatened species.

The whale shark genome reveals patterns of vertebrate gene family evolution

Chondrichthyes (cartilaginous fishes) are fundamental for understanding vertebrate evolution, yet their genomes are understudied. We report long-read sequencing of the whale shark genome to generate the best gapless chondrichthyan genome assembly yet with higher contig contiguity than all other cartilaginous fish genomes, and studied vertebrate genomic evolution of ancestral gene families, immunity, and gigantism. We found a major increase in gene families at the origin of gnathostomes (jawed vertebrates) independent of their genome duplication. We studied vertebrate pathogen recognition receptors (PRRs), which are key in initiating innate immune defense, and found diverse patterns of gene family evolution, demonstrating that adaptive immunity in gnathostomes did not fully displace germline-encoded PRR innovation. We also discovered a new toll-like receptor (TLR29) and three NOD1 copies in the whale shark. We found chondrichthyan and giant vertebrate genomes had decreased substitution rates compared to other vertebrates, but gene family expansion rates varied among vertebrate giants, suggesting substitution and expansion rates of gene families are decoupled in vertebrate genomes. Finally, we found gene families that shifted in expansion rate in vertebrate giants were enriched for human cancer-related genes, consistent with gigantism requiring adaptations to suppress cancer.

Mitochondrial DNA genome evidence for the existence of a third divergent lineage in the western Atlantic Ocean for the bull shark (Carcharhinus leucas)

We report for the first time a highly divergent lineage in the Caribbean Sea for the bull shark (Carcharhinus leucas) based on the analysis of 51 mitochondrial DNA genomes of individuals collected in the western North Atlantic. When comparing the mtDNA control region obtained from the mitogenomes to sequences reported previously for Brazil, the Caribbean lineage remained highly divergent. These results support the existence of a discrete population in Central America due to a phylogeographic break separating the Caribbean Sea from the western North Atlantic, Gulf of Mexico and South America.

Identifying shark species responsible for fisheries depredation off Southeast Queensland, Australia

Anecdotal reports from fishers in Southeast Queensland, Australia suggest that shark depredation is a significant issue, however little is known about which species are responsible for depredating catches. This research aimed to identify depredating species in Southeast Queensland line based fisheries, by undertaking a genetic analysis of depredated samples collected by commercial, charter and recreational fishers. The genetic analysis successfully identified ten depredating sharks, all from the genus Carcharhinus (19.2% success). The species identified using mitochondrial DNA included five C. leucas (bull sharks), two C. plumbeus (sandbar sharks), one C. amboinensis (pigeye shark), one C. brevipinna (spinner shark) and one unconfirmed C. plumbeus/C. altimus (bignose shark). While many species of Carcharhinus have been found to depredate catches in Australia, C. leucas has not been highlighted until this research as a potential problematic species. The optimised protocol allowed for the confident identification of shark species responsible for depredation in fisheries using frozen fish samples donated by fishers.

Lost structural and functional inter-relationships between Ig and TCR loci in mammals revealed in sharks

Immunoglobulins and T cell receptors (TCR) have obvious structural similarities as well as similar immunogenetic diversification and selection mechanisms. Nevertheless, the two receptor systems and the loci that encode them are distinct in humans and classical murine models, and the gene segments comprising each repertoire are mutually exclusive. Additionally, while both B and T cells employ recombination-activating genes (RAG) for primary diversification, immunoglobulins are afforded a supplementary set of activation-induced cytidine deaminase (AID)-mediated diversification tools. As the oldest-emerging vertebrates sharing the same adaptive B and T cell receptor systems as humans, extant cartilaginous fishes allow a potential view of the ancestral immune system. In this review, we discuss breakthroughs we have made in studies of nurse shark (Ginglymostoma cirratum) T cell receptors demonstrating substantial integration of loci and diversification mechanisms in primordial B and T cell repertoires. We survey these findings in this shark model where they were first described, while noting corroborating examples in other vertebrate groups. We also consider other examples where the gnathostome common ancestry of the B and T cell receptor systems have allowed dovetailing of genomic elements and AID-based diversification approaches for the TCR. The cartilaginous fish seem to have retained this T/B cell plasticity to a greater extent than more derived vertebrate groups, but representatives in all vertebrate taxa except bony fish and placental mammals show such plasticity.

A method for making alignments of related protein sequences that share very little similarity; shark interleukin 2 as an example

An optimized alignment of related protein sequences helps to see their important shared features and to deduce their phylogenetic relationships. At low levels of sequence similarity, there are no suitable computer programs for making the best possible alignment. This review summarizes some guidelines for how in such instances, nevertheless, insightful alignments can be made. The method involves, basically, the understanding of molecular family features at both the protein and intron-exon level, and the collection of many related sequences so that gradual differences may be observed. The method is exemplified by identifying and aligning interleukin 2 (IL-2) and related sequences in Elasmobranchii (sharks/rays) and coelacanth, as other authors have expressed difficulty with their identification. From the point of general immunology, it is interesting that the unusual long "leader" sequence of IL-15, already known in other species, is even more impressively conserved in cartilaginous fish. Furthermore, sequence comparisons suggest that IL-2 in cartilaginous fish has lost its ability to bind an IL-2Rα/15Rα receptor chain, which would prohibit the existence of a mechanism for regulatory T cell regulation identical to mammals.

Evolution of vertebrate gill covers via shifts in an ancient Pou3f3 enhancer

Whereas the gill chambers of jawless vertebrates open directly into the environment, jawed vertebrates evolved skeletal appendages that drive oxygenated water unidirectionally over the gills. A major anatomical difference between the two jawed vertebrate lineages is the presence of a single large gill cover in bony fishes versus separate covers for each gill chamber in cartilaginous fishes. Here, we find that these divergent patterns correlate with the pharyngeal arch expression of Pou3f3 orthologs. We identify a deeply conserved Pou3f3 arch enhancer present in humans through sharks but undetectable in jawless fish. Minor differences between the bony and cartilaginous fish enhancers account for their restricted versus pan-arch expression patterns. In zebrafish, mutation of Pou3f3 or the conserved enhancer disrupts gill cover formation, whereas ectopic pan-arch Pou3f3b expression generates ectopic skeletal elements resembling the multimeric covers of cartilaginous fishes. Emergence of this Pou3f3 arch enhancer >430 Mya and subsequent modifications may thus have contributed to the acquisition and diversification of gill covers and respiratory strategies during gnathostome evolution.

Genetic connectivity of the scalloped hammerhead shark Sphyrna lewini across Indonesia and the Western Indian Ocean

Scalloped Hammerhead shark (Sphyrna lewini) is an endangered species which its populations have been declining globally including in Indonesia, the world’s top shark fishing country. However, there is a lack of information on the recent population structure of this species to promote proper management and its conservation status. This study aimed to investigate the genetic diversity, population structure, and connectivity of the S. lewini population, in three major shark landing sites: Aceh (n = 41), Balikpapan (n = 30), and Lombok (n = 29). Meanwhile, additional sequences were retrieved from West Papua (n = 14) and the Western Indian Ocean (n = 65) populations. From the analyses of the mitochondrial CO1 gene, a total of 179 sequences of S. lewini, with an average size of 594 bp, and 40 polymorphic loci in four and eight haplotypes for the Indonesian population and the Western Indian Ocean population were identified. The overall values of genetic diversity were high (h = 0.717; π = 0.013), with the highest values recorded in Aceh (h = 0.668; π = 0.002) and the lowest in Papua (h = 0.143; π = 0.000). On the contrary, the overall value was fairly low in the Western Indian Ocean (h = 0.232; π = 0.001). Furthermore, AMOVA and F_ST showed three significant subdivisions in Indonesia (F_ST = 0.442; P < 0.001), with separated populations for Aceh and West Papua, and mixed between Balikpapan and Lombok (F_ST = 0.044; P = 0.091). In contrast, genetic homogeneity was observed within the population of the Western Indian Ocean (F_ST = –0.013; P = 0.612). The establishment of a haplotype network provided evidence of a significantly different population and a limited genetic distribution between the Indonesian and the Western Indian Ocean populations (F_ST = 0.740; P < 0.001). This study showed the presence of a complex population of S. lewini with limited connectivity only in Indonesia separated from the Western Indian Ocean and requiring specific management measures based on the population structure at the regional level.

Novel developmental bases for the evolution of hypobranchial muscles in vertebrates

BACKGROUND

Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic.

RESULTS

Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle.

CONCLUSIONS

We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Growth trajectories of prenatal embryos of the deep-sea shark Chlamydoselachus anguineus (Chondrichthyes)

Chlamydoselachus anguineus, Garman 1884, commonly called the frilled shark, is a deep-sea shark species occurring up to depths of 1300 m. It is assumed to represent an ancient morphotype of sharks (e.g., terminal mouth opening, more than five gill slits) and thus is often considered to represent plesiomorphic traits for sharks. Therefore, its early ontogenetic developmental traits are important for understanding the evolution of its particular phenotype. Here, we established six stages for prenatal embryos and used linear measurements and geometric morphometrics to analyse changes in shape and size as well as their timing during different embryonic stages. Our results show a change in head shape and a relocation of the mouth opening at a late stage of development. We also detected a negative allometric growth of the head and especially the eye compared to the rest of the body and a sexual dimorphism in total body length, which differs from the known data for adults. A multivariate analysis of covariance shows a significant interaction of shape related to the logarithm of centroid size and developmental stage. Geometric morphometrics results indicate that the head shape changes as a covariate of body size while not accounting for differences between sexes. The growth pattern of stages 32 and 33 indicates a shift in head shape, thus highlighting the moment in development when the jaws start to elongate anteriorly to finally achieve the adult condition of terminal mouth opening rather than retaining the early embryonic subterminal position as is typical for sharks. Thus, the antero-terminal mouth opening of the frilled shark has to be considered a derived feature.

Sequence structure character of IgNAR Sec in whitespotted bamboo shark (Chiloscyllium plagiosum)

Whitespotted bamboo shark (Chiloscyllium plagiosum) is a demersal cartilaginous fish with an adaptive immune system founded upon immunoglobulins. In this manuscript, we characterize the IgNAR of the whitespotted bamboo shark. A newly discovered alternative splicing form of IgNAR Sec (IgNAR_short (ΔC2-C3) Sec) was identified, in which the C1 domain was spliced directly to the C4 domain, the process resulted in a molecule containing three constant domains. However, a single unpaired cysteine remains in the highly flexible hinge region, contributing in the formation of an interchain disulfide bond. Two types of C1 domain were found, and the one lacking a short α-helix showed lower proportion. This finding suggests that short α-helices might be important to the stability of IgNAR. High-throughput sequencing revealed that the percentage of V_NAR types significantly vary between the diverse species of sharks. The variable region of IgNAR (the V_NAR) with small size and stabilization is a potential candidate for immunotherapeutic agents. The structure and stability analysis in this manuscript may be useful in future biomedical applications.

An Ancient, MHC-Linked, Nonclassical Class I Lineage in Cartilaginous Fish

Cartilaginous fishes, or chondrichthyans, are the oldest jawed vertebrates that have an adaptive immune system based on the MHC and Ig superfamily-based AgR. In this basal group of jawed vertebrates, we identified a third nonclassical MHC class I lineage (UDA), which is present in all species analyzed within the two major cartilaginous subclasses, Holocephali (chimaeras) and Elasmobranchii (sharks, skates, and rays). The deduced amino acid sequences of UDA have eight out of nine typically invariant residues that bind to the N and C termini of bound peptide found in most vertebrae classical class I (UAA); additionally, the other predicted 28 peptide-binding residues are perfectly conserved in all elasmobranch UDA sequences. UDA is distinct from UAA in its differential tissue distribution and its lower expression levels and is mono- or oligomorphic unlike the highly polymorphic UAA UDA has a low copy number in elasmobranchs but is multicopy in the holocephalan spotted ratfish (Hydrolagus colliei). Using a nurse shark (Ginglymostoma cirratum) family, we found that UDA is MHC linked but separable by recombination from the tightly linked cluster of UAA, TAP, and LMP genes, the so-called class I region found in most nonmammalian vertebrates. UDA has predicted structural features that are similar to certain nonclassical class I genes in other vertebrates, and, unlike polymorpic classical class I, we anticipate that it may bind to a conserved set of specialized peptides.

Spatio-temporal genetic tagging of a cosmopolitan planktivorous shark provides insight to gene flow, temporal variation and site-specific re-encounters

Migratory movements in response to seasonal resources often influence population structure and dynamics. Yet in mobile marine predators, population genetic consequences of such repetitious behaviour remain inaccessible without comprehensive sampling strategies. Temporal genetic sampling of seasonally recurring aggregations of planktivorous basking sharks, Cetorhinus maximus, in the Northeast Atlantic (NEA) affords an opportunity to resolve individual re-encounters at key sites with population connectivity and patterns of relatedness. Genetic tagging (19 microsatellites) revealed 18% of re-sampled individuals in the NEA demonstrated inter/multi-annual site-specific re-encounters. High genetic connectivity and migration between aggregation sites indicate the Irish Sea as an important movement corridor, with a contemporary effective population estimate (Ne) of 382 (CI = 241-830). We contrast the prevailing view of high gene flow across oceanic regions with evidence of population structure within the NEA, with early-season sharks off southwest Ireland possibly representing genetically distinct migrants. Finally, we found basking sharks surfacing together in the NEA are on average more related than expected by chance, suggesting a genetic consequence of, or a potential mechanism maintaining, site-specific re-encounters. Long-term temporal genetic monitoring is paramount in determining future viability of cosmopolitan marine species, identifying genetic units for conservation management, and for understanding aggregation structure and dynamics.