Half a century of global decline in oceanic sharks and rays.

Overfishing is the primary cause of marine defaunation, yet declines in and increasing extinction risks of individual species are difficult to measure, particularly for the largest predators found in the high seas1-3. Here we calculate two well-established indicators to track progress towards Aichi Biodiversity Targets and Sustainable Development Goals4,5: the Living Planet Index (a measure of changes in abundance aggregated from 57 abundance time-series datasets for 18 oceanic shark and ray species) and the Red List Index (a measure of change in extinction risk calculated for all 31 oceanic species of sharks and rays). We find that, since 1970, the global abundance of oceanic sharks and rays has declined by 71% owing to an 18-fold increase in relative fishing pressure. This depletion has increased the global extinction risk to the point at which three-quarters of the species comprising this functionally important assemblage are threatened with extinction. Strict prohibitions and precautionary science-based catch limits are urgently needed to avert population collapse6,7, avoid the disruption of ecological functions and promote species recovery8,9.

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.

Global status and conservation potential of reef sharks.

Decades of overexploitation have devastated shark populations, leaving considerable doubt as to their ecological status1,2. Yet much of what is known about sharks has been inferred from catch records in industrial fisheries, whereas far less information is available about sharks that live in coastal habitats3. Here we address this knowledge gap using data from more than 15,000 standardized baited remote underwater video stations that were deployed on 371 reefs in 58 nations to estimate the conservation status of reef sharks globally. Our results reveal the profound impact that fishing has had on reef shark populations: we observed no sharks on almost 20% of the surveyed reefs. Reef sharks were almost completely absent from reefs in several nations, and shark depletion was strongly related to socio-economic conditions such as the size and proximity of the nearest market, poor governance and the density of the human population. However, opportunities for the conservation of reef sharks remain: shark sanctuaries, closed areas, catch limits and an absence of gillnets and longlines were associated with a substantially higher relative abundance of reef sharks. These results reveal several policy pathways for the restoration and management of reef shark populations, from direct top-down management of fishing to indirect improvement of governance conditions. Reef shark populations will only have a high chance of recovery by engaging key socio-economic aspects of tropical fisheries.

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.

Whale shark rhodopsin adapted to deep-sea lifestyle by a substitution associated with human disease.

Spectral tuning of visual pigments often facilitates adaptation to new environments, and it is intriguing to study the visual ecology of pelagic sharks with secondarily expanded habitats. The whale shark, which dives into the deep sea of nearly 2,000 meters besides near-surface filter feeding, was previously shown to possess the 'blue-shifted' rhodopsin (RHO), which is a signature of deep-sea adaptation. In this study, our spectroscopy of recombinant whale shark RHO mutants revealed that this blue shift is caused dominantly by an unprecedented spectral tuning site 94. In humans, the mutation at the site causes congenital stationary night blindness (CSNB) by reducing the thermal stability of RHO. Similarly, the RHO of deep-diving whale shark has reduced thermal stability, which was experimentally shown to be achieved by site 178 and 94. RHOs having the natural substitution at site 94 are also found in some Antarctic fishes, suggesting that the blue shift by the substitution at the CSNB site associated with the reduction in thermal stability might be allowed in cold-water deep-sea habitats.

Endothermic physiology of extinct megatooth sharks.

The evolution of the extinct megatooth shark, Otodus megalodon, and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark's gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

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.

Conservation successes and challenges for wide-ranging sharks and rays.

Overfishing is the most significant threat facing sharks and rays. Given the growth in consumption of seafood, combined with the compounding effects of habitat loss, climate change, and pollution, there is a need to identify recovery paths, particularly in poorly managed and poorly monitored fisheries. Here, we document conservation through fisheries management success for 11 coastal sharks in US waters by comparing population trends through a Bayesian state-space model before and after the implementation of the 1993 Fisheries Management Plan for Sharks. We took advantage of the spatial and temporal gradients in fishing exposure and fisheries management in the Western Atlantic to analyze the effect on the Red List status of all 26 wide-ranging coastal sharks and rays. We show that extinction risk was greater where fishing pressure was higher, but this was offset by the strength of management engagement (indicated by strength of National and Regional Plan of Action for sharks and rays). The regional Red List Index (which tracks changes in extinction risk through time) declined in all regions until the 1980s but then improved in the North and Central Atlantic such that the average extinction risk is currently half that in the Southwest. Many sharks and rays are wide ranging, and successful fisheries management in one country can be undone by poorly regulated or unregulated fishing elsewhere. Our study underscores that well-enforced, science-based management of carefully monitored fisheries can achieve conservation success, even for slow-growing species.

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.

Global spatial risk assessment of sharks under the footprint of fisheries.

Effective ocean management and the conservation of highly migratory species depend on resolving the overlap between animal movements and distributions, and fishing effort. However, this information is lacking at a global scale. Here we show, using a big-data approach that combines satellite-tracked movements of pelagic sharks and global fishing fleets, that 24% of the mean monthly space used by sharks falls under the footprint of pelagic longline fisheries. Space-use hotspots of commercially valuable sharks and of internationally protected species had the highest overlap with longlines (up to 76% and 64%, respectively), and were also associated with significant increases in fishing effort. We conclude that pelagic sharks have limited spatial refuge from current levels of fishing effort in marine areas beyond national jurisdictions (the high seas). Our results demonstrate an urgent need for conservation and management measures at high-seas hotspots of shark space use, and highlight the potential of simultaneous satellite surveillance of megafauna and fishers as a tool for near-real-time, dynamic management.

Inner ear development in cyclostomes and evolution of the vertebrate semicircular canals.

Jawed vertebrates have inner ears with three semicircular canals, the presence of which has been used as a key to understanding evolutionary relationships. Ostracoderms, the jawless stem gnathostomes, had only two canals and lacked the lateral canal1-3. Lampreys, which are modern cyclostomes, are generally thought to possess two semicircular canals whereas the hagfishes-which are also cyclostomes-have only a single canal, which used to be regarded as a more primitive trait1,4. However, recent molecular and developmental analyses have strongly supported the monophyly of cyclostomes5-7, which has left the evolutionary trajectory of the vertebrate inner ear unclear8. Here we show the differentiation of the otic vesicle of the lamprey Lethenteron camtschaticum and inshore hagfish Eptatretus burgeri. This is the first time, to our knowledge, that the development of the hagfish inner ear is reported. We found that canal development in the lamprey starts with two depressions-which is reminiscent of the early developmental pattern of the inner ear in modern gnathostomes. These cyclostome otic vesicles show a pattern of expression of regulatory genes, including OTX genes, that is comparable to that of gnathosomes. Although two depressions appear in the lamprey vesicle, they subsequently fuse to form a single canal that is similar to that of hagfishes. Complete separation of the depressions results in anterior and posterior canals in gnathostomes. The single depression of the vesicle in hagfishes thus appears to be a secondarily derived trait. Furthermore, the lateral canal in crown gnathostomes was acquired secondarily-not by de novo acquisition of an OTX expression domain, but by the evolution of a developmental program downstream of the OTX genes.

Global collision-risk hotspots of marine traffic and the world's largest fish, the whale shark.

Marine traffic is increasing globally yet collisions with endangered megafauna such as whales, sea turtles, and planktivorous sharks go largely undetected or unreported. Collisions leading to mortality can have population-level consequences for endangered species. Hence, identifying simultaneous space use of megafauna and shipping throughout ranges may reveal as-yet-unknown spatial targets requiring conservation. However, global studies tracking megafauna and shipping occurrences are lacking. Here we combine satellite-tracked movements of the whale shark, Rhincodon typus, and vessel activity to show that 92% of sharks' horizontal space use and nearly 50% of vertical space use overlap with persistent large vessel (>300 gross tons) traffic. Collision-risk estimates correlated with reported whale shark mortality from ship strikes, indicating higher mortality in areas with greatest overlap. Hotspots of potential collision risk were evident in all major oceans, predominantly from overlap with cargo and tanker vessels, and were concentrated in gulf regions, where dense traffic co-occurred with seasonal shark movements. Nearly a third of whale shark hotspots overlapped with the highest collision-risk areas, with the last known locations of tracked sharks coinciding with busier shipping routes more often than expected. Depth-recording tags provided evidence for sinking, likely dead, whale sharks, suggesting substantial "cryptic" lethal ship strikes are possible, which could explain why whale shark population declines continue despite international protection and low fishing-induced mortality. Mitigation measures to reduce ship-strike risk should be considered to conserve this species and other ocean giants that are likely experiencing similar impacts from growing global vessel traffic.

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.

Embryonic development in the bonnethead (Sphyrna tiburo), a viviparous hammerhead shark.

BACKGROUND: The hammerhead sharks (family Sphyrnidae) are an immediately recognizable group of sharks due to their unique head shape. Though there has long been an interest in hammerhead development, there are currently no explicit staging tables published for any members of the group. The bonnethead Sphyrna tiburo is the smallest member of Sphyrnidae and is abundant in estuarine and nearshore waters in the Gulf of Mexico and Western North Atlantic Ocean. Due to their relative abundance, close proximity to shore, and brief gestation period, it has been possible to collect and document multiple embryonic specimens at progressive stages of development. RESULTS: We present the first comprehensive embryonic staging series for the Bonnethead, a viviparous hammerhead shark. Our stage series covers a period of development from stages that match the vertebrate phylotypic period, from Stage 23, through stages of morphological divergence to complete development at birth-Stage 35). Notably, we use a variety of techniques to document crucial stages that lead to their extreme craniofacial diversity, resulting in the formation of one of the most distinctive characters of any shark species, the cephalofoil or hammer-like head. CONCLUSION: Documenting the development of hard-to-access vertebrates, like this viviparous shark species, offers important information about how new and diverse morphologies arise that otherwise may remain poorly studied. This work will serve as a platform for future comparative developmental research both within sharks and across the phylogeny of vertebrates, underpinning the extreme potential of craniofacial development and morphological diversity in vertebrate animals.

Collagen IV of basement membranes: II. Emergence of collagen IVα345 enabled the assembly of a compact GBM as an ultrafilter in mammalian kidneys.

The collagen IVα345 (Col-IVα345) scaffold, the major constituent of the glomerular basement membrane (GBM), is a critical component of the kidney glomerular filtration barrier. In Alport syndrome, affecting millions of people worldwide, over two thousand genetic variants occur in the COL4A3, COL4A4, and COL4A5 genes that encode the Col-IVα345 scaffold. Variants cause loss of scaffold, a suprastructure that tethers macromolecules, from the GBM or assembly of a defective scaffold, causing hematuria in nearly all cases, proteinuria, and often progressive kidney failure. How these variants cause proteinuria remains an enigma. In a companion paper, we found that the evolutionary emergence of the COL4A3, COL4A4, COL4A5, and COL4A6 genes coincided with kidney emergence in hagfish and shark and that the COL4A3 and COL4A4 were lost in amphibians. These findings opened an experimental window to gain insights into functionality of the Col-IVα345 scaffold. Here, using tissue staining, biochemical analysis and TEM, we characterized the scaffold chain arrangements and the morphology of the GBM of hagfish, shark, frog, and salamander. We found that α4 and α5 chains in shark GBM and α1 and α5 chains in amphibian GBM are spatially separated. Scaffolds are distinct from one another and from the mammalian Col-IVα345 scaffold, and the GBM morphologies are distinct. Our findings revealed that the evolutionary emergence of the Col-IVα345 scaffold enabled the genesis of a compact GBM that functions as an ultrafilter. Findings shed light on the conundrum, defined decades ago, whether the GBM or slit diaphragm is the primary filter.

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 pan-neuronal marker HuC in the olfactory organ suggests the same cell-type specificity of V2R/OlfC expression as for bony fishes, that is, in microvillous neurons. The relatively low number of olfactory receptors in cartilaginous fishes compared with 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.

Development of Shark Single Domain Antibodies Specific for Human α-Fetoprotein and the Multimerization Strategy in Serum Detection.

Sensitive detection of cancer biomarkers can contribute to the timely diagnosis and treatment of diseases. In this study, the whitespotted bamboo sharks were immunized with human α-fetoprotein (AFP), and a phage-displayed variable new antigen receptor (VNAR) single domain antibody library was constructed. Then four unique VNARs (VNAR1, VNAR11, VNAR21, and VNAR25) against AFP were isolated from the library by biopanning for the first time. All of the sequences belong to type II of VNAR, and the VNAR11 was much different from the rest of the three sequences. Then VNAR1 and VNAR11 were selected to fuse with the C4-binding protein α chain (C4bpα) sequence and efficiently expressed in the Escherichia coli system. Furthermore, a VNAR-C4bpα-mediated sandwich chemiluminescence immunoassay (VSCLIA) was developed for the detection of AFP in human serum samples. After optimization, the VSCLIA showed a limit of detection of 0.74 ng/mL with good selectivity and accuracy. Moreover, the results of clinical serum samples detected by the VSCLIA were confirmed by an automatic immunoanalyzer in the hospital, indicating its practical application in actual samples. In conclusion, the novel antibody element VNAR exhibits great potential for immunodiagnosis, and this study also provides a new direction and experimental basis for AFP detection.

Exceptionally preserved shark fossils from Mexico elucidate the long-standing enigma of the Cretaceous elasmobranch Ptychodus.

The fossil fish Ptychodus Agassiz, 1834, characterized by a highly distinctive grinding dentition and an estimated gigantic body size (up to around 10 m), has remained one of the most enigmatic extinct elasmobranchs (i.e. sharks, skates and rays) for nearly two centuries. This widespread Cretaceous taxon is common in Albian to Campanian deposits from almost all continents. However, specimens mostly consist of isolated teeth or more or less complete dentitions, whereas cranial and post-cranial skeletal elements are very rare. Here we describe newly discovered material from the early Late Cretaceous of Mexico, including complete articulated specimens with preserved body outline, which reveals crucial information on the anatomy and systematic position of Ptychodus. Our phylogenetic and ecomorphological analyses indicate that ptychodontids were high-speed (tachypelagic) durophagous lamniforms (mackerel sharks), which occupied a specialized predatory niche previously unknown in fossil and extant elasmobranchs. Our results support the view that lamniforms were ecomorphologically highly diverse and represented the dominant group of sharks in Cretaceous marine ecosystems. Ptychodus may have fed predominantly on nektonic hard-shelled prey items such as ammonites and sea turtles rather than on benthic invertebrates, and its extinction during the Campanian, well before the end-Cretaceous crisis, might have been related to competition with emerging blunt-toothed globidensine and prognathodontine mosasaurs.

Development of branchial ionocytes in embryonic and larval stages of cloudy catshark, Scyliorhinus torazame.

In teleost fish, branchial ionocytes are important sites for osmoregulation and acid-base regulation by maintaining ionic balance in the body fluid. During the early developmental stages before the formation of the gills, teleost ionocytes are localized in the yolk-sac membrane and body skin. By comparing with teleost fish, much less is known about ionocytes in developing embryos of elasmobranch fish. The present study investigated the development of ionocytes in the embryo and larva of cloudy catshark, Scyliorhinus torazame. We first observed ionocyte distribution by immunohistochemical staining with anti-Na+/K+-ATPase (NKA) and anti-vacuolar-type H+-ATPase (V-ATPase) antibodies. The NKA- and V-ATPase-rich ionocytes appeared as single cells in the gill filaments from stage 31, the stage of pre-hatching, while the ionocytes on the body skin and yolk-sac membrane were also observed. From stage 32, in addition to single ionocytes on the gill filaments, some outstanding follicular structures of NKA-immunoreactive cells were developed to fill the inter-filament region of the gill septa. The follicular ionocytes possess NKA in the basolateral membrane and Na+/H+ exchanger 3 in the apical membrane, indicating that they are involved in acid-base regulation like single NKA-rich ionocytes. Three-dimensional analysis and whole-mount immunohistochemistry revealed that the distribution of follicular ionocytes was limited to the rostral side of gill septum. The rostral sides of gill septum might be exposed to faster water flow than caudal side because the gills of sharks gently curved backward. This dissymmetric distribution of follicular ionocytes is considered to facilitate efficient body-fluid homeostasis of catshark embryo.

Northeast Atlantic elasmobranch community on the move: Functional reorganization in response to climate change.

While spatial distribution shifts have been documented in many marine fishes under global change, the responses of elasmobranchs have rarely been studied, which may have led to an underestimation of their potential additional threats. Given their irreplaceable role in ecosystems and their high extinction risk, we used a 24-year time series (1997-2020) of scientific bottom trawl surveys to examine the effects of climate change on the spatial distribution of nine elasmobranch species within Northeast Atlantic waters. Using a hierarchical modeling of species communities, belonging to the joint species distribution models, we found that suitable habitats for four species increased on average by a factor of 1.6 and, for six species, shifted north-eastwards and/or to deeper waters over the past two decades. By integrating species traits, we showed changes in habitat suitability led to changes in the elasmobranchs trait composition. Moreover, communities shifted to deeper waters and their mean trophic level decreased. We also note an increase in the mean community size at maturity concurrent with a decrease in fecundity. Because skates and sharks are functionally unique and dangerously vulnerable to both climate change and fishing, we advocate for urgent considerations of species traits in management measures. Their use would make it better to identify species whose loss could have irreversible impacts in face of the myriad of anthropogenic threats.