The whale shark genome reveals patterns of vertebrate gene family evolution

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.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

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.

Genealogical Diversity of Endogenous Retrovirus in the Jawless Fish Genome

Retroviral integration into ancient vertebrate genomes left traces that can shed light on the early history of viruses. In this study, we explored the early evolution of retroviruses by isolating nine Spuma endogenous retroviruses (ERVs) and one Epsilon ERV from the genomes of Agnatha and Chondrichthyes. Phylogenetic analysis of protein sequences revealed a striking pattern of co-evolution between jawless fish ERV and their host, while shark ERV underwent ancient cross-class viral transmission with jawless fish, ray-finned fish, and amphibians. Nucleotide sequence analysis showed that jawless fish ERV emerged in the Palaeozoic period, relatively later than ray-finned fish ERV. Moreover, codon analysis suggested that the jawless fish ERV employed an infection strategy that mimics the host codon. The genealogical diversity of ERVs in the jawless fish genome highlights the importance of studying different viral species. Overall, our findings provide valuable insights into the evolution of retroviruses and their interactions with their hosts.

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.

Gene losses, parallel evolution and heightened expression confer adaptations to dedicated cleaning behaviour

BACKGROUND

Cleaning symbioses are captivating interspecific interactions in which a cleaner fish removes ectoparasites from its client, contributing to the health and diversity of natural fish communities and aquaculture systems. However, the genetic signatures underlying this specialized behaviour remain poorly explored. To shed light on this, we generated a high-quality chromosome-scale genome of the bluestreak cleaner wrasse Labroides dimidiatus, a dedicated cleaner with cleaning as primary feeding mechanism throughout its life.

RESULTS

Compared with facultative and non-cleaner wrasses, L. dimidiatus was found with notable contractions in olfactory receptors implying their limited importance in dedicated cleaning. Instead, given its distinct tactile pre-conflict strategies, L. dimidiatus may rely more heavily on touch sensory perception, with heightened gene expression in the brain in anticipation of cleaning. Additionally, a reduction in NLR family CARD domain-containing protein 3 might enhance innate immunity of L. dimidiatus, probably assisting to reduce the impacts from parasite infections. In addition, convergent substitutions for a taste receptor and bone development genes across cleaners (L. dimidiatus and facultative cleaners) may provide them with evolved food discrimination abilities and jaw morphology that differentiate them from non-cleaners. Moreover, L. dimidiatus may exhibit specialized neural signal transductions for cleaning, as evidenced by positive selection in genes related to the glutamatergic synapse pathway. Interestingly, numerous glutamate receptors also demonstrated significantly higher expression in L. dimidiatus not engaged in cleaning, as compared to those involved in cleaning. Besides, apparent contractions in L. dimidiatus for protocadherins, which are responsible for neuronal development, may further promote specialized neural signal transductions in this species.

CONCLUSIONS

This study reveals that L. dimidiatus harbours substantial losses in specific gene families, convergent evolutions across cleaners and a large-scale high gene expression in preparation for cleaning, allowing for adaptation to the dedicated cleaning behaviour.

Organized B cell sites in cartilaginous fishes reveal the evolutionary foundation of germinal centers

The absence of germinal centers (GCs) in cartilaginous fishes lies at odds with data showing that nurse sharks can produce robust antigen-specific responses and affinity mature their B cell repertoires. To investigate this apparent incongruity, we performed RNA sequencing on single nuclei, allowing us to characterize the cell types present in the nurse shark spleen, and RNAscope to provide in situ cellular resolution of key marker gene expression following immunization with R-phycoerythrin (PE). We tracked PE to the splenic follicles where it co-localizes with CXCR5high centrocyte-like B cells and a population of putative T follicular helper (Tfh) cells, surrounded by a peripheral ring of Ki67+ AID+ CXCR4+ centroblast-like B cells. Further, we reveal selection of mutations in B cell clones dissected from these follicles. We propose that the B cell sites identified here represent the evolutionary foundation of GCs, dating back to the jawed vertebrate ancestor.

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.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Why are whales big? Genes behind ocean giants

Gigantism is prevalent in animals, but it has never reached more extreme levels than in aquatic mammals such as whales, dolphins, and porpoises. A new study by Silva et al. has uncovered five genes underlying this gigantism, a phenotype with important connections to aging and cancer suppression in long-lived animals.

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.

Uncovering a 500 million year old history and evidence of pseudogenization for TLR15

Introduction

Toll like receptors (TLRs) are at the front line of pathogen recognition and host immune response. Many TLR genes have been described to date with some being found across metazoans while others are restricted to specific lineages. A cryptic member of the TLR gene family, TLR15, has a unique phylogenetic distribution. Initially described in extant species of birds and reptiles, an ortholog has been reported for cartilaginous fish.

Methods

Here, we significantly expanded the evolutionary analysis of TLR15 gene evolution, taking advantage of large genomic and transcriptomic resources available from different lineages of vertebrates. Additionally, we objectively search for TLR15 in lobe-finned and ray-finned fish, as well as in cartilaginous fish and jawless vertebrates.

Results and discussion

We confirm the presence of TLR15 in early branching jawed vertebrates - the cartilaginous fish, as well as in basal Sarcopterygii - in lungfish. However, within cartilaginous fish, the gene is present in Holocephalans (all three families) but not in Elasmobranchs (its sister-lineage). Holocephalans have long TLR15 protein sequences that disrupt the typical TLR structure, and some species display a pseudogene sequence due to the presence of frameshift mutations and early stop codons. Additionally, TLR15 has low expression levels in holocephalans when compared with other TLR genes. In turn, lungfish also have long TLR15 protein sequences but the protein structure is not compromised. Finally, TLR15 presents several sites under negative selection. Overall, these results suggest that TLR15 is an ancient TLR gene and is experiencing ongoing pseudogenization in early-branching vertebrates.

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.

Exploration of the Nurse Shark (Ginglymostoma cirratum) Plasma Immunoproteome Using High-Resolution LC-MS/MS

Many animals of scientific importance lack species-specific reagents (e.g., monoclonal antibodies) for in-depth studies of immune proteins. Mass spectrometry (MS)-based proteomics has emerged as a useful method for monitoring changes in protein abundance and modifications in non-model species. It can be used to quantify hundreds of candidate immune molecules simultaneously without the generation of new reagents. Here, we used MS-based proteomics to identify and quantify candidate immune proteins in the plasma of the nurse shark (Ginglymostoma cirratum), a cartilaginous fish and representative of the most basal extant vertebrate lineage with an immunoglobulin-based immune system. Mass spectrometry-based LC-MS/MS was performed on the blood plasma of nurse sharks immunized with human serum albumin (n=4) or sham immunized (n=1), and sampled at days 0 (baseline control), 1, 2, 3, 5, 7, 14, 21, 28, 25, 42 and 49. An antigen-specific antibody response was experimentally confirmed post-immunization. To provide a high-quality reference to identify proteins, we assembled and annotated a multi-tissue de novo transcriptome integrating long- and short-read sequence data. This comprised 62,682 contigs containing open reading frames (ORFs) with a length >80 amino acids. Using this transcriptome, we reliably identified 626 plasma proteins which were broadly categorized into coagulation, immune, and metabolic functional groups. To assess the feasibility of performing LC-MS/MS proteomics in nurse shark in the absence of species-specific protein annotations, we compared the results to an alternative strategy, mapping peptides to proteins predicted in the genome assembly of a related species, the whale shark (Rhincodon typus). This approach reliably identified 297 proteins, indicating that useful data on the plasma proteome may be obtained in many instances despite the absence of a species-specific reference protein database. Among the plasma proteins defined against the nurse shark transcriptome, fifteen showed consistent changes in abundance across the immunized shark individuals, indicating a role in the immune response. These included alpha-2-macroglobulin (A2M) and a novel protein yet to be characterized in diverse vertebrate lineages. Overall, this study enhances genetic and protein-level resources for nurse shark research and vastly improves our understanding of the elasmobranch plasma proteome, including its remodelling following immune stimulation.

Molecular evolution of hatching enzymes and their paralogous genes in vertebrates

BACKGROUND

Hatching is identified as one of the most important events in the reproduction of oviparous vertebrates. The genes for hatching enzymes, which are vital in the hatching process, are conserved among vertebrates. However, especially in teleost, it is difficult to trace their molecular evolution in detail due to the presence of other C6astacins, which are the subfamily to which the genes for hatching enzymes belong and are highly diverged. In particular, the hatching enzyme genes are diversified with frequent genome translocations due to retrocopy.

RESULTS

In this study, we took advantage of the rapid expansion of whole-genome data in recent years to examine the molecular evolutionary process of these genes in vertebrates. The phylogenetic analysis and the genomic synteny analysis revealed C6astacin genes other than the hatching enzyme genes, which was previously considered to be retained only in teleosts, was also retained in the genomes of basal ray-finned fishes, coelacanths, and cartilaginous fishes. These results suggest that the common ancestor of these genes can be traced back to at least the common ancestor of the Gnathostomata. Moreover, we also found that many of the C6astacin genes underwent multiple gene duplications during vertebrate evolution, and the results of gene expression analysis in frogs implied that genes derived from hatching enzyme genes underwent neo-functionalization.

CONCLUSIONS

In this study, we describe in detail the molecular evolution of the C6astacin gene in vertebrates, which has not been summarized previously. The results revealed the presence of the previously unknown C6astacin gene in the basal-lineage of jawed vertebrates and large-scale gene duplication of hatching enzyme genes in amphibians. The comprehensive investigation reported in this study will be an important basis for studying the molecular evolution of the vertebrate C6astacin genes, hatching enzyme, and its paralogous genes and for identifying these genes without the need for gene expression and functional analysis.

Immunological Diversity Is a Cornerstone of Organismal Defense and Allorecognition across Metazoa

The ongoing arms race between hosts and microbes has fueled the evolution of novel strategies for diversifying the molecules involved in immune responses. Characterization of immune systems from an ever-broadening phylogenetic range of organisms reveals that there are many mechanisms by which this diversity can be generated and maintained. Diversification strategies operate at the level of populations, genomes, genes, and even individual transcripts. Lineage-specific innovations have been cataloged within the immune systems of both invertebrates and vertebrates. Furthermore, somatic diversification of immune receptor genes has now been described in jawless vertebrates and some invertebrate species. In addition to pathogen detection, immunological diversity plays important roles in several distinct allorecognition systems. In this Brief Review, we highlight some of the evolutionary innovations employed by a variety of metazoan species to generate the molecular diversity required to detect a vast array of molecules in the context of both immune response and self/nonself-recognition.