White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability

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 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 genomics illuminates karyotype and sex chromosome evolution of sharks

Chondrichthyes is an important lineage to reconstruct the evolutionary history of vertebrates. Here, we analyzed genome synteny for six chondrichthyan chromosome-level genomes. Our comparative analysis reveals a slow evolutionary rate of chromosomal changes, with infrequent but independent fusions observed in sharks, skates, and chimaeras. The chondrichthyan common ancestor had a proto-vertebrate-like karyotype, including the presence of 18 microchromosome pairs. The X chromosome is a conversed microchromosome shared by all sharks, suggesting a likely common origin of the sex chromosome at least 181 million years ago. We characterized the Y chromosomes of two sharks that are highly differentiated from the X except for a small young evolutionary stratum and a small pseudoautosomal region. We found that shark sex chromosomes lack global dosage compensation but that dosage-sensitive genes are locally compensated. Our study on shark chromosome evolution enhances our understanding of shark sex chromosomes and vertebrate chromosome evolution.

Genome analysis reveals three distinct lineages of the cosmopolitan white shark

The white shark (Carcharodon carcharias) (Linnaeus, 1758), an iconic apex predator occurring in all oceans,1,2 is classified as Vulnerable globally3-with global abundance having dropped to 63% of 1970s estimates,4-and as Critically Endangered in Europe.5 Identification of evolutionary significant units and their management are crucial for conservation,6 especially as the white shark is facing various but often region-specific anthropogenic threats.7,8,9,10,11 Assessing connectivity in a cosmopolitan marine species requires worldwide sampling and high-resolution genetic markers.12 Both are lacking for the white shark, with studies to date typified by numerous but geographically limited sampling, and analyses relying largely on relatively small numbers of nuclear microsatellites,13,14,15,16,17,18,19 which can be plagued by various genotyping artefacts and thus require cautious interpretation.20 Sequencing and computational advances are finally allowing genomes21,22,23 to be leveraged into population studies,24,25,26,27 with datasets comprising thousands of single-nucleotide polymorphisms (SNPs). Here, combining target gene capture (TGC)28 sequencing (89 individuals, 4,000 SNPs) and whole-genome re-sequencing (17 individuals, 391,000 SNPs) with worldwide sampling across most of the distributional range, we identify three genetically distinct allopatric lineages (North Atlantic, Indo-Pacific, and North Pacific). These diverged 100,000-200,000 years ago during the Penultimate Glaciation, when low sea levels, different ocean currents, and water temperatures produced significant biogeographic barriers. Our results show that without high-resolution genomic analyses of samples representative of a species' range,12 the true extent of diversity, presence of past and contemporary barriers to gene flow, subsequent speciation, and local evolutionary events will remain enigmatic.

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.

Fast, parallel, and cache-friendly suffix array construction

PURPOSE

String indexes such as the suffix array (SA) and the closely related longest common prefix (LCP) array are fundamental objects in bioinformatics and have a wide variety of applications. Despite their importance in practice, few scalable parallel algorithms for constructing these are known, and the existing algorithms can be highly non-trivial to implement and parallelize.

METHODS

In this paper we present CAPS-SA, a simple and scalable parallel algorithm for constructing these string indexes inspired by samplesort and utilizing an LCP-informed mergesort. Due to its design, CAPS-SA has excellent memory-locality and thus incurs fewer cache misses and achieves strong performance on modern multicore systems with deep cache hierarchies.

RESULTS

We show that despite its simple design, CAPS-SA outperforms existing state-of-the-art parallel SA and LCP-array construction algorithms on modern hardware. Finally, motivated by applications in modern aligners where the query strings have bounded lengths, we introduce the notion of a bounded-context SA and show that CAPS-SA can easily be extended to exploit this structure to obtain further speedups. We make our code publicly available at https://github.com/jamshed/CaPS-SA .

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

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

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

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

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

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

Whole-Genome Sequencing Analyses Reveal the Whip-like Tail Formation, Innate Immune Evolution, and DNA Repair Mechanisms of Eupleurogrammus muticus

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

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.

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.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information available, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprising 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system (immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes) exhibit expansion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum when compared to mammals. Our functional studies suggest that lack of the gc gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi, as well as the genome evolution in cartilaginous fishes.

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.

Egg Yolk Protein Homologs Identified in Live-Bearing Sharks: Co-Opted in the Lecithotrophy-to-Matrotrophy Shift?

Reproductive modes of vertebrates are classified into two major embryonic nutritional types: yolk deposits (i.e., lecithotrophy) and maternal investment (i.e., matrotrophy). Vitellogenin (VTG), a major egg yolk protein synthesized in the female liver, is one of the molecules relevant to the lecithotrophy-to-matrotrophy shift in bony vertebrates. In mammals, all VTG genes are lost following the lecithotrophy-to-matrotrophy shift, and it remains to be elucidated whether the lecithotrophy-to-matrotrophy shift in nonmammalians is also associated with VTG repertoire modification. In this study, we focused on chondrichthyans (cartilaginous fishes)-a vertebrate clade that underwent multiple lecithotrophy-to-matrotrophy shifts. For an exhaustive search of homologs, we performed tissue-by-tissue transcriptome sequencing for two viviparous chondrichthyans, the frilled shark Chlamydoselachus anguineus and the spotless smooth-hound Mustelus griseus, and inferred the molecular phylogeny of VTG and its receptor very low-density lipoprotein receptor (VLDLR), across diverse vertebrates. As a result, we identified either three or four VTG orthologs in chondrichthyans including viviparous species. We also showed that chondrichthyans had two additional VLDLR orthologs previously unrecognized in their unique lineage (designated as VLDLRc2 and VLDLRc3). Notably, VTG gene expression patterns differed in the species studied depending on their reproductive mode; VTGs are broadly expressed in multiple tissues, including the uterus, in the two viviparous sharks, and in addition to the liver. This finding suggests that the chondrichthyans VTGs do not only function as the yolk nutrient but also as the matrotrophic factor. Altogether, our study indicates that the lecithotrophy-to-matrotrophy shift in chondrichthyans was achieved through a distinct evolutionary process from mammals.

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.

Omics and Remote Homology Integration to Decipher Protein Functionality

In the recent years, several "omics" technologies based on specific biomolecules (from DNA, RNA, proteins, or metabolites) have won growing importance in the scientific field. Despite each omics possess their own laboratorial protocols, they share a background of bioinformatic tools for data integration and analysis. A recent subset of bioinformatic tools, based on available templates or remote homology protocols, allow computational fast and high-accuracy prediction of protein structures. The quickly predict of actually unsolved protein structures, together with late omics findings allow a boost of scientific advances in multiple fields such as cancer, longevity, immunity, mitochondrial function, toxicology, drug design, biosensors, and recombinant protein engineering. In this chapter, we assessed methodological approaches for the integration of omics and remote homology inferences to decipher protein functionality, opening the door to the next era of biological knowledge.

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.

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.

Age group DNA methylation differences in lemon sharks (Negaprion brevirostris): Implications for future age estimation tools

Age information is often non-existent for most shark populations due to a lack of measurable physiological and morphological traits that can be used to estimate age. Recently, epigenetic clocks have been found to accurately estimate age for mammals, birds, and fish. However, since these clocks rely, among other things, on the availability of reference genomes, their application is hampered in non-traditional model organisms lacking such molecular resources. The technique known as Methyl-Sensitive Amplified Polymorphism (MSAP) has emerged as a valid alternative for studying DNA methylation biomarkers when reference genome information is missing, and large numbers of samples need to be processed. Accordingly, the MSAP technique was used in the present study to characterize global DNA methylation patterns in lemon sharks from three different age groups (juveniles, subadults, and adults). The obtained results reveal that, while MSAP analyses lack enough resolution as a standalone approach to infer age in these organisms, the global DNA methylation patterns observed using this technique displayed significant differences between age groups. Overall, these results confer that DNA methylation does change with age in sharks like what has been seen for other vertebrates and that MSAP could be useful as part of an epigenetics pipeline to infer the broad range of ages found in large samples sizes.

The First Genome Survey and De Novo Assembly of the Short Mackerel (Rastrelliger brachysoma) and Indian Mackerel (Rastrelliger kanagurta)

Simple Summary

Mackerel species are commercially important marine species in Southeast Asia, especially short mackerel and Indian mackerel. However, genomic information about them is still limited. Genome survey of these two mackerel species was reported in this study. Next-generation sequencing and comprehensive bioinformatics were performed to obtain the genetic information. The estimated genome size of both species is around 680 Mbp. The heterozygosity of these species was very similar, while the repeat content for Indian mackerel was slightly higher than for short mackerel. Functional annotation also was reported in this study. This is the first reported genome survey and assembly of species in the genus Rastrelliger and could be useful for future comparative genomic studies.

Abstract

Rastrelliger brachysoma (short mackerel) and Rastrelliger kanagurta (Indian mackerel) are commercially important marine species in Southeast Asia. In recent years, numbers of these two species have been decreasing in the wild, and genomic information about them is still limited. We conducted a genome survey of these two mackerel species to acquire essential genomic information using next-generation sequencing data. To obtain this genetic information, comprehensive bioinformatics analyses were performed, including de novo assembly, gene prediction, functional annotation, and phylogenetic analysis. The estimated genome sizes were around 680.14 Mbp (R. brachysoma) and 688.82 Mbp (R. kanagurta). The heterozygosity of these species was very similar (≈0.81), while the repeat content for R. kanagurta (9.30%) was slightly higher than for R. brachysoma (8.30%). Functional annotation indicated that most of the genes predicted in these two species shared very close average amino acid identities (94.06%). The phylogenetic analysis revealed close phylogenetic relationships between these two species and other scombrids. This is the first reported genome survey and assembly of species in the genus Rastrelliger and could be useful for future comparative genomic studies.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.

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.

Elasmobranchs as bioindicators of pollution in the marine environment

Bioindicator species are increasingly valuable in environmental pollution monitoring, and elasmobranch species include many suitable candidates for that role. By measuring contaminants and employing biomarkers of effect in relevant elasmobranch species, scientists may gain important insights about the impacts of pollution in marine ecosystems. This review compiles biomarkers applied in elasmobranchs to assess the effect of pollutants (e.g., metals, persistent organic pollutants, and plastics), and the environmental changes induced by anthropogenic activities (e.g., shifts in marine temperature, pH, and oxygenation). Over 30 biomarkers measured in more than 12 species were examined, including biotransformation biomarkers (e.g., cytochrome P450 1A), oxidative stress-related biomarkers (e.g., superoxide anion, lipid peroxidation, catalase, and vitamins), stress proteins (e.g., heat shock protein 70), reproductive and endocrine biomarkers (e.g., vitellogenin), osmoregulation biomarkers (e.g., trimethylamine N-oxide, Na⁺/K⁺-ATPase, and plasma ions), energetic and neurotoxic biomarkers (e.g., lactate dehydrogenase, lactate, and cholinesterases), and histopathological and morphologic biomarkers (e.g., tissue lesions and gross indices).

Increased risk of cancer in dogs and humans: a consequence of recent extension of lifespan beyond evolutionarily-determined limitations?

Cancer is among the most common causes of death for dogs (and cats) and humans in the developed world, even though it is uncommon in wildlife and other domestic animals. We provide a rationale for this observation based on recent advances in our understanding of the evolutionary basis of cancer. Over the course of evolutionary time, species have acquired and fine-tuned adaptive cancer protective mechanisms that are intrinsically related to their energy demands, reproductive strategies, and expected lifespan. These cancer protective mechanisms are general across species and/or specific to each species and their niche, and they do not seem to be limited in diversity. The evolutionarily acquired cancer-free longevity that defines a species' life history can explain why the relative cancer risk, rate, and incidence are largely similar across most species in the animal kingdom despite differences in body size and life expectancy. The molecular, cellular, and metabolic events that promote malignant transformation and cancerous growth can overcome these adaptive, species-specific protective mechanisms in a small proportion of individuals, while independently, some individuals in the population might achieve exceptional longevity. In dogs and humans, recent dramatic alterations in healthcare and social structures have allowed increasing numbers of individuals in both species to far exceed their species-adapted longevities (by 2-4 times) without allowing the time necessary for compensatory natural selection. In other words, the cancer protective mechanisms that restrain risk at comparable levels to other species for their adapted lifespan are incapable of providing cancer protection over this recent, drastic and widespread increase in longevity.

Evolution of brain-expressed biogenic amine receptors into olfactory trace amine-associated receptors

The family of trace amine-associated receptors (TAARs) is distantly related to G protein-coupled biogenic aminergic receptors. TAARs are found in the brain as well as in the olfactory epithelium where they detect biogenic amines. However, the functional relationship of receptors from distinct TAAR subfamilies and in different species is still uncertain. Here, we perform a thorough phylogenetic analysis of 702 TAAR-like (TARL) and TAAR sequences from 48 species. We show that a clade of Tarl genes has greatly expanded in lampreys, whereas the other Tarl clade consists of only one or two orthologs in jawed vertebrates and is lost in amniotes. We also identify two small clades of Taar genes in sharks related to the remaining Taar genes in bony vertebrates, which are divided into four major clades. We further identify ligands for 61 orphan TARLs and TAARs from sea lamprey, shark, ray-finned fishes, and mammals, as well as novel ligands for two 5-hydroxytryptamine receptor 4 orthologs, a serotonin receptor subtype closely related to TAARs. Our results reveal a pattern of functional convergence and segregation: TARLs from sea lamprey and bony vertebrate olfactory TAARs underwent independent expansions to function as chemosensory receptors, whereas TARLs from jawed vertebrates retain ancestral response profiles and may have similar functions to TAAR1 in the brain. Overall, our data provide a comprehensive understanding of the evolution and ligand recognition profiles of TAARs and TARLs.

Cas9-Based Local Enrichment and Genomics Sequence Revision of Megabase-Sized Shark IgNAR Loci

The 0.8-Mb Ig new Ag receptor (IgNAR) region of the whitespotted bamboo shark (Chiloscyllium plagiosum) is incompletely assembled in Chr_44 of the reference genome. Here we used Cas9-assisted targeting of chromosome segments (CATCH) to enrich the 2 Mb region of the Chr_44 IgNAR loci and sequenced it by PacBio and next-generation sequencing. A fragment >3.13 Mb was isolated intact from the RBCs of sharks. The target was enriched 245.531-fold, and sequences had up to 94% coverage with a 255× mean depth. Compared with the previously published sequences, 20 holes were filled, with a total length of 3508 bp. In addition, we report five potential germline V alleles of IgNAR1 from six sharks that may belong to two clusters of the IgNAR. Our results provide a new method to research the germline of large Ig gene segments, as well as provide the enhanced bamboo shark IgNAR gene loci with fewer gaps.

Bamboo Shark as a Small Animal Model for Single Domain Antibody Production

The development of shark single domain antibodies (sdAbs) is hindered by the high cost and tediousness of large-sized shark farming. Here, we demonstrated white-spotted bamboo sharks (Chiloscyllium plagiosum) being cultivated commercially as a promising small animal model to produce sdAbs. We found that immunoglobulin new antigen receptor (IgNAR) presented in bamboo shark genome, transcriptome, and plasma. Four complete IgNAR clusters including variable domains (vNARs) were discovered in the germline, and the Variable–Joining pair from IgNAR1 cluster was dominant from immune repertoires in blood. Bamboo sharks developed effective immune responses upon green fluorescent protein (GFP), near-infrared fluorescent protein iRFP713, and Freund’s adjuvant immunization revealed by elevated lymphocyte counts and antigen specific IgNAR. Before and after immunization, the complementarity determining region 3 (CDR3) of IgNAR were the major determinant of IgNAR diversity revealed by 400-bp deep sequencing. To prove that bamboo sharks could produce high-affinity IgNAR, we isolated anti-GFP and anti-iRFP713 vNARs with up to 0.3 and 3.8 nM affinities, respectively, from immunized sharks. Moreover, we constructed biparatopic vNARs with the highest known affinities (20.7 pM) to GFP and validated the functions of anti-GFP vNARs as intrabodies in mammalian cells. Taken together, our study will accelerate the discovery and development of bamboo shark sdAbs for biomedical industry at low cost and easy operation.

Origin and evolutionary landscape of Nr2f transcription factors across Metazoa

BACKGROUND

Nuclear Receptor Subfamily 2 Group F (Nr2f) orphan nuclear hormone transcription factors (TFs) are fundamental regulators of many developmental processes in invertebrates and vertebrates. Despite the importance of these TFs throughout metazoan development, previous work has not clearly outlined their evolutionary history.

RESULTS

We integrated molecular phylogeny with comparisons of intron/exon structure, domain architecture, and syntenic conservation to define critical evolutionary events that distinguish the Nr2f gene family in Metazoa. Our data indicate that a single ancestral eumetazoan Nr2f gene predated six main Bilateria subfamilies, which include single Nr2f homologs, here referred to as Nr2f1/2/5/6, that are present in invertebrate protostomes and deuterostomes, Nr2f1/2 homologs in agnathans, and Nr2f1, Nr2f2, Nr2f5, and Nr2f6 orthologs that are found in gnathostomes. Four cnidarian Nr2f1/2/5/6 and three agnathan Nr2f1/2 members are each due to independent expansions, while the vertebrate Nr2f1/Nr2f2 and Nr2f5/Nr2f6 members each form paralogous groups that arose from the established series of whole-genome duplications (WGDs). Nr2f6 members are the most divergent Nr2f subfamily in gnathostomes. Interestingly, in contrast to the other gnathostome Nr2f subfamilies, Nr2f5 has been independently lost in numerous vertebrate lineages. Furthermore, our analysis shows there are differential expansions and losses of Nr2f genes in teleosts following their additional rounds of WGDs.

CONCLUSION

Overall, our analysis of Nr2f gene evolution helps to reveal the origins and previously unrecognized relationships of this ancient TF family, which may allow for greater insights into the conservation of Nr2f functions that shape Metazoan body plans.

On the inference of complex phylogenetic networks by Markov Chain Monte-Carlo

For various species, high quality sequences and complete genomes are nowadays available for many individuals. This makes data analysis challenging, as methods need not only to be accurate, but also time efficient given the tremendous amount of data to process. In this article, we introduce an efficient method to infer the evolutionary history of individuals under the multispecies coalescent model in networks (MSNC). Phylogenetic networks are an extension of phylogenetic trees that can contain reticulate nodes, which allow to model complex biological events such as horizontal gene transfer, hybridization and introgression. We present a novel way to compute the likelihood of biallelic markers sampled along genomes whose evolution involved such events. This likelihood computation is at the heart of a Bayesian network inference method called SnappNet, as it extends the Snapp method inferring evolutionary trees under the multispecies coalescent model, to networks. SnappNet is available as a package of the well-known beast 2 software. Recently, the MCMC_BiMarkers method, implemented in PhyloNet, also extended Snapp to networks. Both methods take biallelic markers as input, rely on the same model of evolution and sample networks in a Bayesian framework, though using different methods for computing priors. However, SnappNet relies on algorithms that are exponentially more time-efficient on non-trivial networks. Using simulations, we compare performances of SnappNet and MCMC_BiMarkers. We show that both methods enjoy similar abilities to recover simple networks, but SnappNet is more accurate than MCMC_BiMarkers on more complex network scenarios. Also, on complex networks, SnappNet is found to be extremely faster than MCMC_BiMarkers in terms of time required for the likelihood computation. We finally illustrate SnappNet performances on a rice data set. SnappNet infers a scenario that is consistent with previous results and provides additional understanding of rice evolution.

The whale shark genome reveals patterns of vertebrate gene family evolution

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

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Using an Ensemble to Identify and Classify Macroalgae Antimicrobial Peptides

The rapid spread of multi-drug resistant microbes has lead researchers to discover natural alternative remedies such as antimicrobial peptides (AMPs). In the first line of defense, AMPs display a broad spectrum of potent activity against multi-resistant pathogenic bacteria, viruses, fungi, and even cancer. AMPs can be further characterised into families according to amino acid composition, secondary structure, and function. However, despite recent advancements in rapid computational methods for AMP prediction from various mammalian, aquatic, and terrestrial species, there is limited information regarding their presence, functional roles, and family type from marine macroalgae. In this paper, we present a promising two-tier ensemble of heterogeneous machine learning models that integrates seven well-known machine learning classifiers to predict AMPs from macroalgae. The first tier of the ensemble consists of a suite of binary classifiers that identify AMPs from protein sequence data which are then forwarded to a second-tier multi-class ensemble to characterise their functional family type. The two-tier ensemble was successfully used to identify 39 putative AMP sequences in 12 macroalgae species from three different phyla groups. The approach we describe is not limited to AMPs and can also be applied to search sequence data for other types of proteins.

Ecomorphological, space, and mineral relations of dermal denticles in angular angel shark (Squatina guggenheim)

Shark skin is predominantly specialized for swimming and protection, with the dermal denticle being the main structure associated with these abilities. The dermal denticle is a mineral structure with a unique morphology for each species, which allows its use as a taxonomic tool. Few studies have investigated the microscopy aspects of skin and dermal denticles, considering the high diversity of sharks. Here, we investigated the three-dimensional morphoquantitative aspects and mineral composition of dermal denticles in different regions of the angular angel shark, Squatina guggenheim, using scanning electron microscopy and dispersive energy system. With the microscopy, we were able to observe that the dermal denticle morphology changes according to the area it is located. It was possible to describe the dermal denticles individually, from root to the crown, highlighting all of their individualities. Through the dispersive energy system, we showed the proportions of each mineral found in the denticle, by area, demonstrating the composition and the particularities of crown, body, and root, where whitlockite was described for the first time in elasmobranchs. In this way, the present study presented the specificities of the dermal denticles of S. guggenheim, as well sought to understand the different structure functions for the animal, thus assisting future research in animal morphology.

An inconvenient tooth: Evaluating female choice in multiple paternity using an evolutionarily and ecologically important vertebrate clade

Understanding mating systems is a pillar of behavioural ecology, placing the complex interactions between females and males into a reproductive context. The field of multiple paternity, the phenomenon whereby many sires contribute to an individual litter, has traditionally viewed females as passive players in a male–male competitive framework. With the emergence of feminist perspectives in ecological fields, novel alternative mechanisms and evolutionary theories across invertebrate and vertebrate taxa recognize females are active stakeholders in the reproductive process. Despite their evolutionary significance, ecological diversity and myriad reproductive modes, elasmobranch (sharks, skates and rays) research lags behind other fields regarding complex biological processes, such as multiple paternity which is often ascribed to convenience polyandry. Here, we layout hypotheses and resynthesize multiple paternity literature from a female and life history perspective to highlight how alternative mechanisms influence the predominance of multiple paternity across elasmobranchs. We draw upon parallels in other invertebrate and vertebrate taxa to demonstrate how female elasmobranchs can influence multiple paternity outcomes that benefit their reproductive success. Our article challenges dogma that has resulted from years of dismissing the female perspective as important and provides a framework for future advancement using more holistic approaches to studying mating systems.

New Species Can Broaden Myelin Research: Suitability of Little Skate, Leucoraja erinacea

Although myelinated nervous systems are shared among 60,000 jawed vertebrates, studies aimed at understanding myelination have focused more and more on mice and zebrafish. To obtain a broader understanding of the myelination process, we examined the little skate, Leucoraja erinacea. The reasons behind initiating studies at this time include: the desire to study a species belonging to an out group of other jawed vertebrates; using a species with embryos accessible throughout development; the availability of genome sequences; and the likelihood that mammalian antibodies recognize homologs in the chosen species. We report that the morphological features of myelination in a skate hatchling, a stage that supports complex behavioral repertoires needed for survival, are highly similar in terms of: appearances of myelinating oligodendrocytes (CNS) and Schwann cells (PNS); the way their levels of myelination conform to axon caliber; and their identity in terms of nodal and paranodal specializations. These features provide a core for further studies to determine: axon-myelinating cell communication; the structures of the proteins and lipids upon which myelinated fibers are formed; the pathways used to transport these molecules to sites of myelin assembly and maintenance; and the gene regulatory networks that control their expressions.

Microbiome Analyses Demonstrate Specific Communities Within Five Shark Species

Profiles of symbiotic microbial communities (“microbiomes”) can provide insight into the natural history and ecology of their hosts. Using high throughput DNA sequencing of the 16S rRNA V4 region, microbiomes of five shark species in South Florida (nurse, lemon, sandbar, Caribbean reef, and tiger) have been characterized for the first time. The microbiomes show species specific microbiome composition, distinct from surrounding seawater. Shark anatomical location (gills, teeth, skin, cloaca) affected the diversity of microbiomes. An in-depth analysis of teeth communities revealed species specific microbial communities. For example, the genus Haemophilus, explained 7.0% of the differences of the teeth microbiomes of lemon and Caribbean reef sharks. Lemon shark teeth communities (n = 11) contained a high abundance of both Vibrio (10.8 ± 26.0%) and Corynebacterium (1.6 ± 5.1%), genera that can include human pathogenic taxa. The Vibrio (2.8 ± 6.34%) and Kordia (3.1 ± 6.0%) genera and Salmonella enterica (2.6 ± 6.4%) were the most abundant members of nurse shark teeth microbial communities. The Vibrio genus was highly represented in the sandbar shark (54.0 ± 46.0%) and tiger shark (5.8 ± 12.3%) teeth microbiomes. The prevalence of genera containing potential human pathogens could be informative in shark bite treatment protocols and future research to confirm or deny human pathogenicity. We conclude that South Florida sharks host species specific microbiomes that are distinct from their surrounding environment and vary due to differences in microbial community composition among shark species and diversity and composition among anatomical locations. Additionally, when considering the confounding effects of both species and location, microbial community diversity and composition varies.

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

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

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

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

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

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

Multimodal Imaging and Analysis of the Neuroanatomical Organization of the Primary Olfactory Inputs in the Brownbanded Bamboo Shark, Chiloscyllium punctatum

There is currently a limited understanding of the morphological and functional organization of the olfactory system in cartilaginous fishes, particularly when compared to bony fishes and terrestrial vertebrates. In this fish group, there is a clear paucity of information on the characterization, density, and distribution of olfactory receptor neurons (ORNs) within the sensory olfactory epithelium lining the paired olfactory rosettes, and their functional implications with respect to the hydrodynamics of incurrent water flow into the nares. This imaging study examines the brownbanded bamboo shark Chiloscyllium punctatum (Elasmobranchii) and combines immunohistochemical labeling using antisera raised against five G-protein α-subunits (Gα_s/olf, Gα_q/11/14, Gα_i–1/2/3, Gα_i–3, Gα_o) with light and electron microscopy, to characterize the morphological ORN types present. Three main ORNs (“long”, “microvillous” and “crypt-like”) are confirmed and up to three additional microvilli-bearing types are also described; “Kappe-like” (potential or homologous “Kappe” as in teleosts), “pear-shaped” and “teardrop-shaped” cells. These morphotypes will need to be confirmed molecularly in the future. Using X-ray diffusible iodine-based contrast-enhanced computed tomography (diceCT), high-resolution scans of the olfactory rosettes, olfactory bulbs (OBs), peduncles, and telencephalon reveal a lateral segregation of primary olfactory inputs within the OBs, with distinct medial and lateral clusters of glomeruli, suggesting a potential somatotopic organization. However, most ORN morphotypes are found to be ubiquitously distributed within the medial and lateral regions of the olfactory rosette, with at least three microvilli-bearing ORNs labeled with anti-Gα_o found in significantly higher densities in lateral lamellae [in lateral lamellae] and on the anterior portion of lamellae (facing the olfactory cavity). These microvilli-bearing ORN morphotypes (microvillous, “Kappe-like,” “pear-shaped,” and “teardrop-shaped”) are the most abundant across the olfactory rosette of this species, while ciliated ORNs are less common and crypt cells are rare. Spatial simulations of the fluid dynamics of the incurrent water flow into the nares and within the olfactory cavities indicate that the high densities of microvilli-bearing ORNs located within the lateral region of the rosette are important for sampling incoming odorants during swimming and may determine subsequent tracking behavior.

Histologic Observations of Dermal Wound Healing in a Free-Ranging Blacktip Shark from the Southeastern U.S. Atlantic Coast: A Case Report

Despite indications that sharks have an exceptional capacity to heal from traumatic injuries, no detailed microscopic observations of integumental wounds have been reported for sharks. This study details the histopathological features of such wounds in a free-ranging shark. An adult male Blacktip Shark Carcharhinus limbatus was collected in 2017 during fisheries-independent sampling efforts in the coastal southeastern U.S. Atlantic. The shark had numerous lesions on his head, torso, and left pectoral fin that were compatible with shark bites. Representative samples from two wounds on the head were examined by light microscopy. The epidermal changes included hyperplasia and hypertrophy; intracellular edema; and the absence of goblet cells, denticles, and chromatophores. In the dermis, fibrinocellular exudation, granulation tissue, and marked skeletal muscle necrosis and regeneration were observed. The above features were comparable to wound healing in bony fish, albeit minor differences were found. Although this case documents exceptionally good regeneration of skeletal muscle in the shark, we found no evidence of unique morphological healing patterns. Further studies on wound healing are needed because recent molecular and genetic findings do suggest evolutionary adaptations enhancing healing in sharks.