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

Phylogenetic Analysis of 590 Species Reveals Distinct Evolutionary Patterns of Intron-Exon Gene Structures Across Eukaryotic Lineages

Introns are highly prevalent in most eukaryotic genomes. Despite the accumulating evidence for benefits conferred by the possession of introns, their specific roles and functions, as well as the processes shaping their evolution, are still only partially understood. Here, we explore the evolution of the eukaryotic intron-exon gene structure by focusing on several key features such as the intron length, the number of introns, and the intron-to-exon length ratio in protein-coding genes. We utilize whole-genome data from 590 species covering the main eukaryotic taxonomic groups and analyze them within a statistical phylogenetic framework. We found that the basic gene structure differs markedly among the main eukaryotic groups, with animals, and particularly chordates, displaying intron-rich genes, compared with plants and fungi. Reconstruction of gene structure evolution suggests that these differences evolved prior to the divergence of the main phyla and have remained mostly conserved within groups. We revisit the previously reported association between the genome size and the mean intron length and report that this association differs considerably among phyla. Analyzing a large and diverse dataset of species with whole-genome information while applying advanced modeling techniques allowed us to obtain a global evolutionary perspective. Our findings may indicate that introns play different molecular and evolutionary roles in different organisms.

Structure and function of neurohypophysial hormones

Vasopressin (VP) and oxytocin (OXT) are neuropeptides that are synthesized in the hypothalamus and stored in/secreted from the neurohypophysis. Although VP and OXT were initially characterized as osmoregulatory and reproductive hormones, respectively, these peptides exert versatile actions not only in peripheral organs but also in the central nervous system via multiple G protein-coupled receptors. Orthologous peptides and receptors have been identified in various animal phyla, reflecting an ancient origin of this hormone family. The aim of this review is to provide basic information on this hormone family and to propose matters to be addressed in future studies. In the earlier sections of this review, we summarize the historical aspect of VP/OXT research as well as the basic features of hormonal peptides and corresponding receptors. The latter sections describe VP/OXT family peptides and their receptors in nonmammalian species, including invertebrates, to introduce the evolutionary aspect of this hormone family. By integrating knowledge from both general and comparative endocrinology perspectives, we highlight current and future research trends about the VP/OXT system.

Mitochondrial Genome and Phylogenetic Analysis of the Narrownose Smooth-Hound Shark Mustelus schmitti Springer, 1939

Southern Brazil is home to a large biodiversity of elasmobranchs from the Brazilian coast. Several genera and species of small sharks of the Triakidae family live in this marine environment. Studies on these shark species are scarce, with few genetic data and little information on animal population structures. The present study aimed to sequence the complete mitochondrial genome (mtDNA) of the endangered species Mustelus schmitti (narrownose smooth-hound shark) and to perform a phylogenetic analysis of the Triakidae family. The mtDNA sequenced here was 16,764 bp long and possessed the usual 13 mitochondrial protein coding genes (PCGs), 22 tRNAs, two rRNAs (12S and 16S) and a large D-loop DNA sequence, presenting an overall organization similar to other species from the genus Mustelus. Phylogenetic analyses were performed using a dataset containing this new mtDNA and 59 other mitochondrial genomes of the Carcharhiniformes species (including 14 from the Triakidae family), using the Maximum Likelihood (ML) method. All the species of the Triakidae family were clustered into a monophyletic topology group. In addition, polyphyly was observed in Galeorhinus galeus, Hemiatrakis japanica, Triakis megalopterus and Triakis semifasciata. In conclusion, this study contributes to a deeper understanding of the genetic diversity of sharks and represents an important step towards the conservation of these endangered animals.

The mechanoreceptors in hatchling and adult Elasmobranch skin

The skin is the most extensive organ in vertebrates, composed of two layers: the epidermis and the dermis. Sensory axons originating from the dorsal root ganglia innervate the skin mechanoreceptors in the dermis. Elasmobranchs, which appeared 380 million years ago, are characterized by rough skin composed of dermal denticles. While we know about the epidermis and dermis of elasmobranchs, we do not know much about the presence or abundance of mechanoreceptors in their skin. Using the classic histological hematoxylin and eosin method, we examined the skins of hatchling embryos and adults Batoidea (skates and rays) and Selachimorpha (modern sharks). Our histology findings provide substantial evidence to identify structures with similar morphology to traditional mammalian and reptilian mechanoreceptors like Pacinian and Meissner corpuscles. An interesting observation was the presence of Pacinian in the skin of Batoidea but not in the skin of a Selachimorpha Squalus shark.

Characterization of gonadotropins and their receptors in a chondrichthyan, Scyliorhinus canicula, fills a gap in the understanding of their coevolution

In Gnathostomes, reproduction is mainly controlled by the hypothalamic-pituitary-gonadal (HPG) axis, with the involvement of the pituitary gonadotropic hormones (GTH), follicle-stimulating hormone (FSH) and luteinizing hormone (LH), which activate their cognate receptors, FSHR and LHR, expressed in gonads. Each GTH consists of a common α subunit and of a specific FSHβ or LHβ subunit. Chondrichthyes (holocephalans and elasmobranchs) is a sister group of bony vertebrates. This position is highly favorable for the understanding of the evolution of endocrine regulations of reproduction among gnathostomes. Surprisingly, the characterization of gonadotropins and their receptors is still limited in chondrichthyes. In the present study, GTH and GTHR sequences have been identified from several chondrichthyan genomes, and their primary structures were analyzed relative to human orthologs. 3D models of GTH/GTHR interaction were built, highlighting the importance of the receptor hinge region for ligand recognition. Functional hormone-receptor interactions have been studied in HEK cells using the small-spotted catshark (Scyliorhinus canicula) recombinant proteins and showed that LHR was specifically activated by LH whereas FSHR was activated by both FSH and LH. Expression profiles of GTHs and their receptors were explored by real-time PCR, in situ hybridization and immunohistochemistry during spermatogenesis, along the male genital tract and other tissues, as well as in some female tissues for comparison. Tissue-expression analyses showed that the highest levels were observed for fshr transcripts in testis and ovary and for lhr in specific extragonadal tissues. The two receptors were expressed at all stages of spermatogenesis by both germ cells and somatic cells, including undifferentiated spermatogonia, spermatocytes, spermatids, somatic precursors and Sertoli cells; differentiated Leydig cells being absent in the testis of S. canicula. Receptors were also expressed by the lymphomyeloid epigonal tissue and the testicular tubules. These results, suggest a wide range of gonadotropin-regulated functions in Elasmobranchs, as well as functional redundancy during spermatogenesis. These extended functions are discussed in an evolutionary context in which the specificity of gonadotropin signaling must have contributed to the evolution of gonadal cells' morphology and function.

The Palaeozoic assembly of the holocephalan body plan far preceded post-Cretaceous radiations into the ocean depths

Among cartilaginous fishes, Holocephali represents the species-depauperate, morphologically conservative sister to sharks, rays and skates and the last survivor of a once far greater Palaeozoic and Mesozoic diversity. Currently, holocephalan diversity is concentrated in deep-sea species, suggesting that this lineage might contain relictual diversity that now persists in the ocean depths. However, the relationships of living holocephalans to their extinct relatives and the timescale of their diversification remain unclear. Here, we reconstruct the evolutionary history of holocephalans using comprehensive morphological and DNA sequence datasets. Our results suggest that crown holocephalans entered and diversified in deep (below 1000 m) ocean waters after the Cretaceous-Palaeogene mass extinction, contrasting with the hypothesis that this ecosystem has acted as a refugium of ancient cartilaginous fishes. These invasions were decoupled from the evolution of key features of the holocephalan body plan, including crushing dentition, a single frontal clasper, and holostylic jaw suspension, during the Palaeozoic Era. However, these invasions considerably postdated the appearance of extant holocephalan families 150 million years ago during a major period of biotic turnover in oceans termed the Mesozoic Marine Revolution. These results clarify the origins of living holocephalans as the recent diversification of a single surviving clade among numerous Palaeozoic lineages.

Understanding vertebrate immunity through comparative immunology

Evolutionary immunology has entered a new era. Classical studies, using just a handful of model animal species, combined with clinical observations, provided an outline of how innate and adaptive immunity work together to ensure tissue homeostasis and to coordinate the fight against infections. However, revolutionary advances in cellular and molecular biology, genomics and methods of genetic modification now offer unprecedented opportunities. They provide immunologists with the possibility to consider, at unprecedented scale, the impact of the astounding phenotypic diversity of vertebrates on immune system function. This Perspective is intended to highlight some of the many interesting, but largely unexplored, biological phenomena that are related to immune function among the roughly 60,000 existing vertebrate species. Importantly, hypotheses arising from such wide-ranging comparative studies can be tested in representative and genetically tractable species. The emerging general principles and the discovery of their evolutionarily selected variations may inspire the future development of novel therapeutic strategies for human immune disorders.

The evolution of tenascins

BACKGROUND

The evolution of extracellular matrix is tightly linked to the evolution of organogenesis in metazoans. Tenascins are extracellular matrix glycoproteins of chordates that participate in integrin-signaling and morphogenetic events. Single tenascins are encoded by invertebrate chordates, and multiple tenascin paralogs are found in vertebrates (designated tenascin-C, tenascin-R, tenascin-W and tenascin-X) yet, overall, the evolution of this family has remained unclear.

RESULTS

This study examines the genomes of hemichordates, cephalochordates, tunicates, agnathans, cartilaginous fishes, lobe-finned fishes, ray-finned fishes and representative tetrapods to identify predicted tenascin proteins. We comprehensively assess their evolutionary relationships by sequence conservation, molecular phylogeny and examination of conservation of synteny of the encoding genes. The resulting new evolutionary model posits the origin of tenascin in an ancestral chordate, with tenascin-C-like and tenascin-R-like paralogs emerging after a whole genome duplication event in an ancestral vertebrate. Tenascin-X appeared following a second round of whole genome duplication in an ancestral gnathostome, most likely from duplication of the gene encoding the tenascin-R homolog. The fourth gene, encoding tenascin-W (also known as tenascin-N), apparently arose from a local duplication of tenascin-R.

CONCLUSIONS

The diversity of tenascin paralogs observed in agnathans and gnathostomes has evolved through selective retention of novel genes that arose from a combination of whole genome and local duplication events. The evolutionary appearance of specific tenascin paralogs coincides with the appearance of vertebrate-specific cell and tissue types where the paralogs are abundantly expressed, such as the endocranium and facial skeleton (tenascin-C), an expanded central nervous system (tenascin-R), and bone (tenascin-W).

Comparative Tissue Identification and Characterization of Long Non-Coding RNAs in the Globally Distributed Blue Shark Prionace glauca

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and serve crucial regulatory functions in both animals and plants. Nevertheless, there is limited understanding of lncRNAs and their patterns of expression and roles in sharks. In the current study, we systematically identified and characterized lncRNAs in the blue shark (Prionace glauca) from four tissues (liver, spleen, muscle, and kidney) using high-throughput sequencing and bioinformatics tools. A total of 21,932 high-confidence lncRNAs were identified, with 8984 and 3067 stably and tissue-specific expressed lncRNAs, respectively. In addition, a total of 45,007 differentially expressed (DE) lncRNAs were obtained among tissues, with kidney versus muscle having the largest numbers across tissues. DE lncRNAs trans target protein-coding genes were predicted, and functional gene ontology enrichment of these genes showed GO terms such as muscle system processes, cellular/metabolic processes, and stress and immune responses, all of which correspond with the specific biological functions of each tissue analyzed. These results advance our knowledge of lncRNAs in sharks and present novel data on tissue-specific lncRNAs, providing key information to support future functional shark investigations.

The sensory shark: high-quality morphological, genomic and transcriptomic data for the small-spotted catshark Scyliorhinus canicula reveal the molecular bases of sensory organ evolution in jawed vertebrates

Cartilaginous fishes (chimaeras and elasmobranchs -sharks, skates and rays) hold a key phylogenetic position to explore the origin and diversifications of jawed vertebrates. Here, we report and integrate reference genomic, transcriptomic and morphological data in the small-spotted catshark Scyliorhinus canicula to shed light on the evolution of sensory organs. We first characterise general aspects of the catshark genome, confirming the high conservation of genome organisation across cartilaginous fishes, and investigate population genomic signatures. Taking advantage of a dense sampling of transcriptomic data, we also identify gene signatures for all major organs, including chondrichthyan specializations, and evaluate expression diversifications between paralogs within major gene families involved in sensory functions. Finally, we combine these data with 3D synchrotron imaging and in situ gene expression analyses to explore chondrichthyan-specific traits and more general evolutionary trends of sensory systems. This approach brings to light, among others, novel markers of the ampullae of Lorenzini electro-sensory cells, a duplication hotspot for crystallin genes conserved in jawed vertebrates, and a new metazoan clade of the Transient-receptor potential (TRP) family. These resources and results, obtained in an experimentally tractable chondrichthyan model, open new avenues to integrate multiomics analyses for the study of elasmobranchs and jawed vertebrates.

Gluconeogenesis in the yolk syncytial layer-like tissue of cloudy catshark (Scyliorhinus torazame)

Glucose has important roles in the development of zebrafish, the vertebrate animal model; however, in most oviparous animals, the amount of maternally provided glucose in the yolk is scarce. For these reasons, developing animals need some ways to supplement glucose. Recently, it was found that developing zebrafish, a teleost fish, undergo gluconeogenesis in the yolk syncytial layer (YSL), an extraembryonic tissue that surrounds the yolk. However, teleost YSL is evolutionarily unique, and it is not clear how other vertebrates supplement glucose. In this study, we used cloudy catshark (or Torazame catshark), an elasmobranch species which possesses a YSL-like tissue during development, and sought for possible gluconeogenic activities in this tissue. In their yolk sac, glucose increased, and our isotope tracking analysis detected gluconeogenic activities with glycerol most preferred substrate. In addition, many of gluconeogenic genes were expressed at the YSL-like tissue, suggesting that cloudy catshark engages in gluconeogenesis in this tissue. The gluconeogenesis in teleost YSL and a similar tissue in elasmobranch species implies conserved mechanisms of yolk metabolism between these two lineages. Future studies on other vertebrate taxa will be helpful to understand the evolutionary changes in the modes of yolk metabolism that vertebrates have experienced.

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.

The origins and drivers of sexual size dimorphism in sharks

While sexual size dimorphism (SSD) is abundant in nature, there is huge variation in both the intensity and direction of SSD. SSD results from a combination of sexual selection for large male size, fecundity selection for large female size and ecological selection for either. In most vertebrates, it is variation in the intensity of male-male competition that primarily underlies variation in SSD. In this study, we test four hypotheses regarding the adaptive value of SSD in sharks-considering the potential for each of fecundity, sexual, ecological selection and reproductive mode as the primary driver of variation in SSD between species. We also estimate past macroevolutionary shifts in SSD direction/intensity through shark phylogeny. We were unable to find evidence of significant SSD in early sharks and hypothesise that SSD is a derived state in this clade, that has evolved independently of SSD observed in other vertebrates. Moreover, there is no significant relationship between SSD and fecundity, testes mass or oceanic depth in sharks. However, there is evidence to support previous speculation that reproductive mode is an important determinant of interspecific variation in SSD in sharks. This is significant as in most vertebrates sexual selection is thought to be the primary driver of SSD trends, with evidence for the role of fecundity selection in other clades being inconsistent at best. While the phylogenetic distribution of SSD among sharks is superficially similar to that observed in other vertebrate clades, the relative importance of selective pressures underlying its evolution appears to differ.

Three foxg1 paralogues in lampreys and gnathostomes-brothers or cousins?

Foxg1 is a key regulator of the early development of the vertebrate forebrain and sensory organs. In this study, we describe for the first time three foxg1 paralogues in lamprey, representative of one of two basally diverged lineages of vertebrates-the agnathans. We also first describe three foxg1 genes in sterlet-representative of one of the evolutionarily ancient clades of gnathostomes. According to the analysis of local genomic synteny, three foxg1 genes of agnathans and gnathostomes have a common origin as a result of two rounds of genomic duplications in the early evolution of vertebrates. At the same time, it is difficult to reliably establish pairwise orthology between foxg1 genes of agnathans and gnathostomes based on the analysis of phylogeny and local genomic synteny, as well as our studies of the spatiotemporal expression of foxg1 genes in the river lamprey Lampetra fluviatilis and the sterlet Acipenser ruthenus. Thus, the appearance of three foxg1 paralogues in agnathans and gnathostomes could have occurred either as a result of two rounds of duplication of the vertebrate common ancestor genome (2R hypothesis) or as a result of the first common round followed by subsequent independent polyploidizations in two evolutionary lineages (1R hypothesis).

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.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Mass Spectrometry Analysis of Shark Skin Proteins

The mucus layer covering the skin of fish has several roles, including protection against pathogens and mechanical damage in which proteins play a key role. While proteins in the skin mucus layer of various common bony fish species have been explored, the proteins of shark skin mucus remain unexplored. In this pilot study, we examine the protein composition of the skin mucus in spiny dogfish sharks and chain catsharks through mass spectrometry (NanoLC-MS/MS). Overall, we identified 206 and 72 proteins in spiny dogfish (Squalus acanthias) and chain catsharks (Scyliorhinus retifer), respectively. Categorization showed that the proteins belonged to diverse biological processes and that most proteins were cellular albeit a significant minority were secreted, indicative of mucosal immune roles. The secreted proteins are reviewed in detail with emphasis on their immune potentials. Moreover, STRING protein-protein association network analysis showed that proteins of closely related shark species were more similar as compared to a more distantly related shark and a bony fish, although there were also significant overlaps. This study contributes to the growing field of molecular shark studies and provides a foundation for further research into the functional roles and potential human biomedical implications of shark skin mucus proteins.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage short reads 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 mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the 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 in line with the assembly length but slightly below the genome size of other representatives of 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 remarkable large coastal shark.

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.

Genomic Characteristics of Okamejei kenojei and the Implications to Its Evolutionary Biology Study

Okamejei kenojei is an important economic species widely distributed in shallow coastal waters of the western North Pacific. In this study, the whole-genome survey analysis of O. kenojei was conducted to reveal its genomic characteristics. The genome size was estimated to be 2027.44 Mb, the repeat sequence content was 44.90%, and the heterozygous ratio was 1.04%. The mitochondrial genome excavated from the sequencing data was 16,974 bp, and it can form the closed circular molecule. The phylogenetic tree based on 13 protein-coding gene sequences supported the validity of Okamejei and assisted the conclusion that Raja porosa was the junior synonym of O. kenojei. Plenty of potential microsatellite loci were identified, and the distribution frequency was estimated to be approximately 236.3 SSRs per Mb. Among all motif types of microsatellites, the dinucleotide repeats were dominant (82.59%), followed by the trinucleotide repeats (8.05%), tetranucleotide repeats (5.80%), pentanucleotide repeats (2.83%), and hexanucleotide repeats (0.72%). The results of the present study could not only provide useful information for understanding the genome structure and functional characteristics of O. kenojei, but also lay the foundation for the subsequent mapping of the whole genome.

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.

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

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

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

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Phylogenetic and functional properties of hagfish neurohypophysial hormone receptors distinct from their jawed vertebrate counterparts

Vertebrate neurohypophysial hormones, i.e., vasopressin- and oxytocin-family peptides, exert versatile physiological actions via distinct G protein-coupled receptors. The neurohypophysial hormone receptor (NHR) family was classically categorized into four subtypes (V1aR, V1bR, V2R and OTR), while recent studies have identified seven subtypes (V1aR, V1bR, V2aR, V2bR, V2cR, V2dR and OTR; V2aR corresponds to the conventional V2R). The vertebrate NHR family were diversified via multiple gene duplication events at different scales. Despite intensive research effort in non-osteichthyes vertebrates such as cartilaginous fish and lamprey, the molecular phylogeny of the NHR family has not been fully understood. In the present study, we focused on the inshore hagfish (Eptatretus burgeri), another group of cyclostomes, and Arctic lamprey (Lethenteron camtschaticum) for comparison. Two putative NHR homologs, which were previously identified only in silico, were cloned from the hagfish and designated as ebV1R and ebV2R. In vitro, ebV1R, as well as two out of five Arctic lamprey NHRs, increased intracellular Ca2+ in response to exogenous neurohypophysial hormones. None of the examined cyclostome NHRs altered intracellular cAMP levels. Transcripts of ebV1R were detected in multiple tissues including the brain and gill, with intense hybridization signals in the hypothalamus and adenohypophysis, while ebV2R was predominantly expressed in the systemic heart. Similarly, Arctic lamprey NHRs showed distinct expression patterns, underscoring the multifunctionality of VT in the cyclostomes as in the gnathostomes. These results and exhaustive gene synteny comparisons provide new insights into the molecular and functional evolution of the neurohypophysial hormone system in vertebrates.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

N-terminal domain influences steroid activation of the Atlantic sea lamprey corticoid receptor

Lampreys are jawless fish that evolved about 550 million years ago at the base of the vertebrate line. Modern lampreys contain a corticoid receptor (CR), the common ancestor of the glucocorticoid receptor (GR) and mineralocorticoid receptor (MR), which first appear in cartilaginous fish, such as sharks. Until recently, 344 amino acids at the amino terminus of adult lamprey CR were not present in the lamprey CR sequence in GenBank. A search of the recently sequenced lamprey germline genome identified two CR sequences, CR1 and CR2, containing the 344 previously un-identified amino acids. CR1 also contains a novel four amino acid insertion in the DNA-binding domain (DBD). We studied corticosteroid and progesterone activation of CR1 and CR2 and found their strongest response was to 11-deoxycorticosterone and 11-deoxycortisol, the two circulating corticosteroids in lamprey. Based on steroid specificity, both CRs are close to elephant shark MR and distant from elephant shark GR. HEK293 cells that were transfected with full-length CR1 or CR2 and the MMTV promoter have about 3-fold higher steroid-mediated activation compared to HEK293 cells transfected with these CRs and the TAT3 promoter. Deletion of the amino-terminal domain (NTD) of lamprey CR1 and CR2 to form truncated CRs decreased transcriptional activation by about 70% in HEK293 cells that were transfected with MMTV, but increased transcription by about 6-fold in cells transfected with TAT3. This indicated that the promoter has an important effect on NTD regulation of transcriptional activation of the CR by steroids. Our results also indicate that the entire lamprey CR sequence is needed for an accurate determination of steroid-mediated transcription.

Nitrogen transporters along the intestinal spiral valve of cloudy catshark (Scyliorhinus torazame): Rhp2, Rhbg, UT

As part of their osmoregulatory strategy, marine elasmobranchs retain large quantities of urea to balance the osmotic pressure of the marine environment. The main source of nitrogen used to synthesize urea comes from the digestion and absorption of food across the gastrointestinal tract. In this study we investigated possible mechanisms of nitrogen movement across the spiral valve of the cloudy catshark (Scyliorhinus torazame) through the molecular identification of two Rhesus glycoprotein ammonia transporters (Rhp2 and Rhbg) and a urea transporter (UT). We used immunohistochemistry to determine the cellular localizations of Rhp2 and UT. Within the spiral valve, Rhp2 was expressed along the apical brush-border membrane, and UT was expressed along the basolateral membrane and the blood vessels. The mRNA abundance of Rhp2 was significantly higher in all regions of the spiral valve of fasted catsharks compared to fed catsharks. The mRNA abundance of UT was significantly higher in the anterior spiral valve of fasted catsharks compared to fed. The mRNA transcript of four ornithine urea cycle (OUC) enzymes were detected along the length of the spiral valve and in the renal tissue, indicating the synthesis of urea via the OUC occurs in these tissues. The presence of Rhp2, Rhbg, and UT along the length of the spiral valve highlights the importance of ammonia and urea movement across the intestinal tissues, and increases our understanding of the mechanisms involved in maintaining whole-body nitrogen homeostasis in the cloudy catshark.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Evolution of Vertebrate Hormones and Their Receptors: Insights from Non-Osteichthyan Genomes

Homeostatic control and reproductive functions of humans are regulated at the molecular levels largely by peptide hormones secreted from endocrine and/or neuroendocrine cells in the central nervous system and peripheral organs. Homologs of those hormones and their receptors function similarly in many vertebrate species distantly related to humans, but the evolutionary history of the endocrine system involving those factors has been obscured by the scarcity of genome DNA sequence information of some taxa that potentially contain their orthologs. Focusing on non-osteichthyan vertebrates, namely jawless and cartilaginous fishes, this article illustrates how investigating genome sequence information assists our understanding of the diversification of vertebrate gene repertoires in four broad themes: (a) the presence or absence of genes, (b) multiplication and maintenance of paralogs, (c) differential fates of duplicated paralogs, and (d) the evolutionary timing of gene origins.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Segment-Dependent Gene Expression Profiling of the Cartilaginous Fish Nephron Using Laser Microdissection for Functional Characterization of Nephron at Segment Levels

For adaptation to a high salinity marine environment, cartilaginous fishes have evolved a ureosmotic strategy. They have a highly elaborate "four-loop nephron" in the kidney, which is considered to be important for reabsorption of urea from the glomerular filtrate to maintain a high concentration of urea in the body. However, the function and regulation, generally, of the "four-loop nephron" are still largely unknown due to the complicated configuration of the nephron and its many subdivided segments. Laser microdissection (LMD) followed by RNA-sequencing (RNA-seq) analysis is a powerful technique to obtain segment-dependent gene expression profiles. In the present study, using the kidney of cloudy catshark, Scyliorhinus torazame, we tested several formaldehyde-free and formaldehyde-based fixatives to optimize the fixation methods. Fixation by 1% neutral buffered formalin for 15 min resulted in sufficient RNA and structural integrities, which allowed LMD clipping of specific nephron segments and subsequent RNA-seq analysis. RNA-seq from the LMD samples of the second-loop, the fourth-loop, and the five tubular segments in the bundle zone revealed a number of specific membrane transporter genes that can characterize each segment. Among them, we examined expressions of the Na ⁺ -coupled cotransporters abundantly expressed in the second loop samples. Although the proximal II segment of the second loop is known for the elimination of excess solutes, the present results imply that the PII segment is also crucial for reabsorption of valuable solutes. Looking ahead to future studies, the segment-dependent gene expression profiling will be a powerful technique for unraveling the renal mechanisms and regulation in euryhaline elasmobranchs.

The first crystal structure of CD8αα from a cartilaginous fish

Introduction

Cartilaginous fishes are the most evolutionary-distant vertebrates from mammals and possess an immunoglobulin (Ig)- and T cell-mediated adaptive immunity. CD8 is the hallmark receptor of cytotoxic T cells and is required for the formation of T cell receptor-major histocompatibility complex (TCR-MHC) class I complexes.

Methods

RACE PCR was used to obtain gene sequences. Direct dilution was applied for the refolding of denatured recombinant CD8 protein. Hanging-drop vapor diffusion method was performed for protein crystallization.

Results

In this study, CD8α and CD8β orthologues (termed ScCD8α and ScCD8β) were identified in small-spotted catshark (Scyliorhinus canicula). Both ScCD8α and ScCD8β possess an extracellular immunoglobulin superfamily (IgSF) V domain as in previously identified CD8 proteins. The genes encoding CD8α and CD8β are tandemly linked in the genomes of all jawed vertebrates studied, suggesting that they were duplicated from a common ancestral gene before the divergence of cartilaginous fishes and other vertebrates. We determined the crystal structure of the ScCD8α ectodomain homodimer at a resolution of 1.35 Å and show that it exhibits the typical topological structure of CD8α from endotherms. As in mammals, the homodimer formation of ScCD8αα relies upon interactions within a hydrophobic core although this differs in position and amino acid composition. Importantly, ScCD8αα shares the canonical cavity required for interaction with peptide-loaded MHC I in mammals. Furthermore, it was found that ScCD8α can co-immunoprecipitate with ScCD8β, indicating that it can form both homodimeric and heterodimeric complexes.

Conclusion

Our results expand the current knowledge of vertebrate CD8 dimerization and the interaction between CD8α with p/MHC I from an evolutionary perspective.

Characterization of mitochondrial genome of Indian Ocean blue-spotted maskray, Neotrygon indica and its phylogenetic relationship within Dasyatidae Family

The present study characterized complete mitochondrial genome of Blue-spotted maskray, Neotrygon indica and studied the evolutionary relationship of the species within the Dasyatidae family. The total length of the mitogenome was 17,974 bp including 37 genes and a non-coding control region. The average frequency of nucleotides in protein-coding genes was A: 29.1 %, T: 30.2 %, G: 13.0 % and C: 27.7 % with AT content of 59.3 %. The values of AT and GC skewness were -0.018 and -0.338, respectively. Comparative analyses showed a large number of average synonymous substitutions per synonymous site (Ks) in gene NADH4 (5.07) followed by NADH5 (4.72). High values of average number of non-synonymous substitutions per non-synonymous site (Ka) were observed in genes ATPase8 (0.54) and NADH2 (0.44). Genes NADH4L and NADH2 showed high interspecific genetic distance values of 0.224 ± 0.001 and 0.213 ± 0.002, respectively. Heat map analysis showed variation in codon usage among different species of the Dasyatidae family. The phylogenetic tree showed a sister relationship between the Dasyatinae and the Neotrygoninae subfamilies. Neotrygon indica formed as a sister species to the clade consisting of N. varidens and N. orientalis. Based on the present results, Neotrygon indica could have diverged from the common ancestor of the two latter in the Plio-Pleistocene. The present study showed distinct characteristics of N. indica from its congeners through comparative mitogenomics.

Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

In vitro and in vivo gene transfer in the cloudy catshark Scyliorhinus torazame

Cartilaginous fishes have various unique physiological features such as a cartilaginous skeleton and a urea-based osmoregulation strategy for adaptation to their marine environment. Also, because they are a sister group of bony vertebrates, understanding their unique features is important from an evolutionary perspective. However, genetic engineering based on gene functions as well as cellular behavior has not been effectively utilized in cartilaginous fishes. This is partly because their reproductive strategy involves internal fertilization, which results in difficulty in microinjection into fertilized eggs at the early developmental stage. Here, to identify efficient gene transfer methods in cartilaginous fishes, we examined the effects of various methods both in vitro and in vivo using the cloudy catshark, a candidate model cartilaginous fish species. In all methods, green fluorescent protein (GFP) expression was used to evaluate exogenous gene transfer. First, we examined gene transfer into primary cultured cells from cloudy catshark embryos by lipofection, polyethylenimine (PEI) transfection, adenovirus infection, baculovirus infection, and electroporation. Among the methods tested, lipofection, electroporation, and baculovirus infection enabled the successful transfer of exogenous genes into primary cultured cells. We then attempted in vivo transfection into cloudy catshark embryos by electroporation and baculovirus infection. Although baculovirus-injected groups did not show GFP fluorescence, electroporation successfully introduced GFP into muscle cells. Furthermore, we succeeded in GFP transfer into adult tissues by electroporation. The in vitro and in vivo gene transfer methods that worked in this study may open ways for genetic manipulation including knockout experiments and cellular lineage analysis in cartilaginous fishes.

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

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

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

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

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.

Little skate genome provides insights into genetic programs essential for limb-based locomotion

The little skate Leucoraja erinacea, a cartilaginous fish, displays pelvic fin driven walking-like behavior using genetic programs and neuronal subtypes similar to those of land vertebrates. However, mechanistic studies on little skate motor circuit development have been limited, due to a lack of high-quality reference genome. Here, we generated an assembly of the little skate genome, with precise gene annotation and structures, which allowed post-genome analysis of spinal motor neurons (MNs) essential for locomotion. Through interspecies comparison of mouse, skate and chicken MN transcriptomes, shared and divergent gene expression profiles were identified. Comparison of accessible chromatin regions between mouse and skate MNs predicted shared transcription factor (TF) motifs with divergent ones, which could be used for achieving differential regulation of MN-expressed genes. A greater number of TF motif predictions were observed in MN-expressed genes in mouse than in little skate. These findings suggest conserved and divergent molecular mechanisms controlling MN development of vertebrates during evolution, which might contribute to intricate gene regulatory networks in the emergence of a more sophisticated motor system in tetrapods.

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

BACKGROUND

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is listed as Vulnerable on the International Union for Conservation of Nature Red List (v2.3) due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts for both wild and ex situpopulations. A high-quality genome would also open avenues to investigate the genetic basis of the species' exceptionally long life span.

FINDINGS

We produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture. We produced a 2.37-Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNA sequencing-assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. To assess the utility of the high-quality assembly for species conservation, we performed a low-coverage resequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. We further identified putatively deleterious mutations to be monitored.

CONCLUSIONS

We establish a high-quality chromosome-level reference genome for A. gigantea and one of the most complete turtle genomes available. We show that low-coverage whole-genome resequencing, for which alignment to the reference genome is a necessity, is a powerful tool to assess the population structure of the wild population and reveal the geographic origins of ex situ individuals relevant for genetic diversity management and rewilding efforts.

Long-term monitoring of egg-laying cycle using ultrasonography reveals the reproductive dynamics of circulating sex steroids in an oviparous catshark, Scyliorhinus torazame

The many diverse reproductive strategies of elasmobranchs (sharks, skates and rays) from lecithotrophic oviparity to matrotrophic viviparity have attracted significant research attention. However, the endocrine control of elasmobranch reproduction is less well-documented largely due to their reproductive characteristics, such as a long reproductive cycle, and/or repeated internal fertilization using stored sperm in oviparous species. In the present study, for the first time, we succeeded in non-invasive monitoring of the continuing egg-laying cycle of the cloudy catshark Scyliorhinus torazame using portable ultrasound devices. Furthermore, long-term simultaneous monitoring of the egg-laying cycle and measurement of plasma sex steroids revealed cycling patterns of estradiol-17β (E2), testosterone (T) and progesterone (P4). In particular, a decline in T followed by a reciprocal surge in plasma P4 were consistently observed prior to the appearance of the capsulated eggs, implying that P4 is likely associated with the ovulation and/or egg-case formation. While the cycling pattern of E2 was not as apparent as those of T and P4, threshold levels of E2 (>5 ng/mL) and T (>1 ng/mL) appeared to be crucial in the continuation of egg-laying cycle. The possibility to trace the dynamics of plasma sex steroids in a single individual throughout the reproductive cycles makes the catshark a useful model for regulatory and mechanistic studies of elasmobranch reproduction.

First observation of the mating, egg-laying and hatching behaviour of a captive female Kong skate, Okamejei kenojei (Müller & Henle, 1841)

This study provides a novel record of the reproductive behaviour of the Kong skate (Okamejei kenojei) in captivity. These skates were found to mate and deposit eggs at a temperature of 16.5 ± 0.5°C. The results showed that 76.13% of the eggs possessed one yolk, 0.77% of the eggs possessed two yolks and 23.11% of the eggs had no yolk (N = 1043). The deposition of non-yolk and double-yolk eggs was random. A total of 100 eggs were collected. After nearly 92 ± 5 days of incubation, 28 eggs failed to hatch, and 72 skates were successfully hatched with a female-to-male ratio of 1:1 (P > 0.05). The results enrich our knowledge of the reproduction in cartilaginous fishes and can inform management and conservation strategies for this species.

Squalomix: shark and ray genome analysis consortium and its data sharing platform

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Intron size minimisation in teleosts

Background

Spliceosomal introns are parts of primary transcripts that are removed by RNA splicing. Although introns apparently do not contribute to the function of the mature transcript, in vertebrates they comprise the majority of the transcribed region increasing the metabolic cost of transcription. The persistence of long introns across evolutionary time suggests functional roles that can offset this metabolic cost. The teleosts comprise one of the largest vertebrate clades. They have unusually compact and variable genome sizes and provide a suitable system for analysing intron evolution.

Results

We have analysed intron lengths in 172 vertebrate genomes and show that teleost intron lengths are relatively short, highly variable and bimodally distributed. Introns that were long in teleosts were also found to be long in mammals and were more likely to be found in regulatory genes and to contain conserved sequences. Our results argue that intron length has decreased in parallel in a non-random manner throughout teleost evolution and represent a deviation from the ancestral state.

Conclusion

Our observations indicate an accelerated rate of intron size evolution in the teleosts and that teleost introns can be divided into two classes by their length. Teleost intron sizes have evolved primarily as a side-effect of genome size evolution and small genomes are dominated by short introns (<256 base pairs). However, a non-random subset of introns has resisted this process across the teleosts and these are more likely have functional roles in all vertebrate clades.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08760-w).

Parallel evolution of ameloblastic scpp genes in bony and cartilaginous vertebrates

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5′ exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates.

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.

Remarkable diversity of vomeronasal type 2 receptor (OlfC) genes of basal ray-finned fish and its evolutionary trajectory in jawed vertebrates

The vomeronasal type 2 receptor (V2R, also called OlfC) multigene family is found in a broad range of jawed vertebrates from cartilaginous fish to tetrapods. V2Rs encode receptors for food-related amino acids in teleost fish, whereas for peptide pheromones in mammals. In addition, V2Rs of teleost fish are phylogenetically distinct from those of tetrapods, implying a drastic change in the V2R repertoire during terrestrial adaptation. To understand the process of diversification of V2Rs in vertebrates from "fish-type" to "tetrapod-type", we conducted an exhaustive search for V2Rs in cartilaginous fish (chimeras, sharks, and skates) and basal ray-finned fish (reedfish, sterlet, and spotted gar), and compared them with those of teleost, coelacanth, and tetrapods. Phylogenetic and synteny analyses on 1897 V2Rs revealed that basal ray-finned fish possess unexpectedly higher number of V2Rs compared with cartilaginous fish, implying that V2R gene repertoires expanded in the common ancestor of Osteichthyes. Furthermore, reedfish and sterlet possessed various V2Rs that belonged to both "fish-type" and "tetrapod-type", suggesting that the common ancestor of Osteichthyes possess "tetrapod-type" V2Rs although they inhabited underwater environments. Thus, the unexpected diversity of V2Rs in basal ray-finned fish may provide insight into how the olfaction of osteichthyan ancestors adapt from water to land.