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

Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

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.

Loss of noggin1, a classic embryonic inducer gene, in elasmobranchs

Secreted proteins of the Noggin family serve as pivotal regulators of early development and cell differentiation in all multicellular animals, including vertebrates. Noggin1 was identified first among all Noggins. Moreover, it was described as the first known embryonic inducer specifically secreted by the Spemann organizer and capable of inducing a secondary body axis when expressed ectopically. In the classical default model of neural induction, Noggin1 is presented as an antagonist of BMP signalling, playing a role as a neural inducer. Additionally, Noggin1 is involved in the dorsalization of embryonic mesoderm and later controls the differentiation of various tissues, including muscles, bones, and neural crest derivatives. Hitherto, noggin1 was found in all studied vertebrates. Here, we report the loss of noggin1 in elasmobranchs (sharks, rays and skates), which is a unique case among vertebrates. noggin2 and noggin4 retained in this group and studied in the embryos of the grey bamboo shark Chiloscyllium griseum revealed similarities in expression patterns and functional properties with their orthologues described in other vertebrates. The loss of noggin1 in elasmobranchs may be associated with histological features of the formation of their unique internal cartilaginous skeleton, although additional research is required to establish functional connections between these events.

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.

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.

Ontogenetic morphometry of the brown smoothhound shark Mustelus henlei with implications for ecology and evolution

The central tenet of ecomorphology links ecological and morphological variation through the process of selection. Traditionally used to rationalise morphological differences between taxa, an ecomorphological approach is increasingly being utilised to study morphological differences expressed through ontogeny. Elasmobranchii (sharks, rays and skates) is one clade in which such ontogenetic shifts in body form have been reported. Such studies are limited to a relatively small proportion of total elasmobranch ecological and morphological diversity, and questions remain regarding the extent to which ecological selection are driving observed morphometric trends. In this study, we report ontogenetic growth trajectories obtained via traditional linear morphometrics from a large data set of the brown smoothhound shark (Mustelus henlei). We consider various morphological structures including the caudal, dorsal and pectoral fins, as well as several girth measurements. We use an ecomorphological approach to infer the broad ecological characteristics of this population and refine understanding of the selective forces underlying the evolution of specific morphological structures. We suggest that observed scaling trends in M. henlei are inconsistent with migratory behaviour, but do not contradict a putative trophic niche shift. We also highlight the role of predation pressure and sex-based ecological differences in driving observed trends in morphometry, a factor which has previously been neglected when considering the evolution of body form in sharks.

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.

Shark mandible evolution reveals patterns of trophic and habitat-mediated diversification

Environmental controls of species diversity represent a central research focus in evolutionary biology. In the marine realm, sharks are widely distributed, occupying mainly higher trophic levels and varied dietary preferences, mirrored by several morphological traits and behaviours. Recent comparative phylogenetic studies revealed that sharks present a fairly uneven diversification across habitats, from reefs to deep-water. We show preliminary evidence that morphological diversification (disparity) in the feeding system (mandibles) follows these patterns, and we tested hypotheses linking these patterns to morphological specialisation. We conducted a 3D geometric morphometric analysis and phylogenetic comparative methods on 145 specimens representing 90 extant shark species using computed tomography models. We explored how rates of morphological evolution in the jaw correlate with habitat, size, diet, trophic level, and taxonomic order. Our findings show a relationship between disparity and environment, with higher rates of morphological evolution in reef and deep-water habitats. Deep-water species display highly divergent morphologies compared to other sharks. Strikingly, evolutionary rates of jaw disparity are associated with diversification in deep water, but not in reefs. The environmental heterogeneity of the offshore water column exposes the importance of this parameter as a driver of diversification at least in the early part of clade history.

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.

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.

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.

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.