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

Functional characterization of follicle-stimulating hormone and luteinizing hormone receptors in cloudy catshark, Scyliorhinus torazame

The follicle-stimulating hormone receptor (FSHR) and luteinizing hormone receptor (LHR) in cloudy catshark were cloned, and recombinant FSHR and LHR were expressed for characterization. Ventral lobe extract (VLE) from the pituitary contains homologous FSH and LH, and it stimulated the cAMP signaling of FSHR and LHR dose-dependently. Two transcript variants of LHR (LHR-L with exon 10 and LHR-S without) were identified, and LHR-S was the dominant form with higher basal cAMP activity without VLE stimulation. Among various developmental stages of follicles, FSHR expression was mainly associated with the pre-vitellogenic and early white follicles. When follicles were recruited into vitellogenesis, the expression of FSHR decreased while of LHR was upregulated reciprocally, suggesting that LHR may also be responsible for the control of vitellogenesis in chondrichthyans. The expression of LHR-L was upregulated among maturing follicles before ovulation, indicating LHR-L could have a specific role in receiving the LH surge signal for final maturation. Plasma LH-like activity was transiently increased prior to the progesterone (P4)-surge and testosterone-drop at the beginning of P4-phase, supporting a pituitary control of follicle-maturation via LH signaling in chondrichthyans. The expression of follicular LHR was downregulated during the P4-phase when LH-like activity was high, indicating that the LH-dependent downregulation of LHR is conserved in chondrichthyans as it is in other vertebrate lineages. (213 words).

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.

Cost-effort analysis of Baited Remote Underwater Video (BRUV) and environmental DNA (eDNA) in monitoring marine ecological communities

Monitoring the diversity and distribution of species in an ecosystem is essential to assess the success of restoration strategies. Implementing biomonitoring methods, which provide a comprehensive assessment of species diversity and mitigate biases in data collection, holds significant importance in biodiversity research. Additionally, ensuring that these methods are cost-efficient and require minimal effort is crucial for effective environmental monitoring. In this study we compare the efficiency of species detection, the cost and the effort of two non-destructive sampling techniques: Baited Remote Underwater Video (BRUV) and environmental DNA (eDNA) metabarcoding to survey marine vertebrate species. Comparisons were conducted along the Sussex coast upon the introduction of the Nearshore Trawling Byelaw. This Byelaw aims to boost the recovery of the dense kelp beds and the associated biodiversity that existed in the 1980s. We show that overall BRUV surveys are more affordable than eDNA, however, eDNA detects almost three times as many species as BRUV. eDNA and BRUV surveys are comparable in terms of effort required for each method, unless eDNA analysis is carried out externally, in which case eDNA requires less effort for the lead researchers. Furthermore, we show that increased eDNA replication yields more informative results on community structure. We found that using both methods in conjunction provides a more complete view of biodiversity, with BRUV data supplementing eDNA monitoring by recording species missed by eDNA and by providing additional environmental and life history metrics. The results from this study will serve as a baseline of the marine vertebrate community in Sussex Bay allowing future biodiversity monitoring research projects to understand community structure as the ecosystem recovers following the removal of trawling fishing pressure. Although this study was regional, the findings presented herein have relevance to marine biodiversity and conservation monitoring programs around the globe.

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.

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.

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.

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.

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.