De novo assembly of the kidney and spleen transcriptomes of the cosmopolitan blue shark, Prionace glauca

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

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

A multi-tissue de novo transcriptome assembly and relative gene expression of the vulnerable freshwater salmonid Thymallus ligericus

Freshwater ecosystems are among the most endangered ecosystems worldwide. While numerous taxa are on the verge of extinction as a result of global changes and direct or indirect anthropogenic activity, genomic and transcriptomic resources represent a key tool for comprehending species' adaptability and serve as the foundation for conservation initiatives. The Loire grayling, Thymallus ligericus, is a freshwater European salmonid endemic to the upper Loire River basin. The species is comprised of fragmented populations that are dispersed over a small area and it has been identified as a vulnerable species. Here, we provide a multi-tissue de novo transcriptome assembly of T. ligericus. The completeness and integrity of the transcriptome were assessed before and after redundancy removal with lineage-specific libraries from Eukaryota, Metazoa, Vertebrata, and Actinopterygii. Relative gene expression was assessed for each of the analyzed tissues, using the de novo assembled transcriptome and a genome-based analysis using the available T. thymallus genome as a reference. The final assembly, with a contig N50 of 1221 and Benchmarking Universal Single-Copy Orthologs (BUSCO) scores above 94%, is made accessible along with structural and functional annotations and relative gene expression of the five tissues (NCBI SRA and FigShare databases). This is the first transcriptomic resource for this species, which provides a foundation for future research on this and other salmonid species that are increasingly exposed to environmental stressors.

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.

De novo assembly and annotation of the whole transcriptome of Muraenesox cinereus

In this study, we generated the transcriptome of Muraenesox cinereus from four combined tissues (muscle, sexual, liver and heart) using high-seq sequencing technology. De novo assembly was performed using Trinity software and a total of 62,125,296 high-quality clean reads were assembled and clustered into 75,862 unigene with an N50 of 2034 nt. After annotation, 43,157 unigenes had significant hit in Nr database. And then, 24,510 unigenes were annotated into three GO categories: biological processes, cellular components, and molecular functions. Moreover, 33,032 unigenes were mapped 25 different clusters of eukaryotic proteins. Furthermore, we predicted the structure of all unigenes using HMMER and MISA software, respectively. The result showed that a total of 33,183 nucleotide sequences of coding regions (direction of the sequences is 5'- > 3') were confirmed to the protein database and a total of 29,487 simple sequence repeats (SSR) were identified. The whole transcriptome is an important foundation for future genomic research on the M. cinereus, it could provide comprehensively understanding and further characterizations of transcriptomes of non-model organisms.

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.

Comparative eye and liver differentially expressed genes reveal monochromatic vision and cancer resistance in the shortfin mako shark (Isurus oxyrinchus)

The shortfin mako, Isurus oxyrinchus is an oceanic pelagic shark found worldwide in tropical and subtropical waters. However, the understanding of its biology at molecular level is still incipient. We sequenced the messenger RNA isolated from eye and liver tissues. De novo transcriptome yielded a total of 705,940 transcripts. A total of 3774 genes were differentially expressed (DEGs), with 1612 in the eye and 2162 in the liver. Most DEGs in the eye were related to structural and signaling functions, including nonocular and ocular opsin genes, whereas nine out of ten most overexpressed genes in the liver were related to tumor suppression, wound healing, and human diseases. Furthermore, DEGs findings provide insights on the monochromatic shark vision and a repertory of cancer-related genes, which may be insightful to elucidate shark resistance to cancer. Therefore, our results provide valuable sequence resources for future functional and population studies.

Conservation genetics of elasmobranchs of the Mexican Pacific Coast, trends and perspectives

One of the most critical threats to biodiversity is the high extinction rate driven by human activities. Reducing extinction rates requires the implementation of conservation programmes based on robust scientific data. Elasmobranchs are important ecological components of the ocean, and several species sustain substantial economic activities. Unfortunately, elasmobranchs are one of the most threatened and understudied animal taxa. The Mexican Pacific Coast (MPC) is a region with high elasmobranch diversity and is the seat of major elasmobranch fisheries. But it is also a developing region with several conservation and management challenges which require national and international attention. Here, we review the conservation genetics literature of elasmobranchs from the MPC. We present a synthesis of the works using samples from the region and emphasize the main gaps and biases in these data. In addition, we discuss the benefits and challenges of generating genomic information to improve the management and conservation of an elasmobranch biodiversity hotspot in a developing country. We found 47 elasmobranch genetic articles that cover <30% of the elasmobranch diversity in the region. These studies mainly used mitochondrial DNA sequences to analyse the genetic structure of commercially important and abundant species of the order Carcharhiniformes. Some of these papers also assessed mating systems, demographic parameters, and taxonomic uncertainties, all of which are important topics for efficient management decisions. In terms of conservation genetics, elasmobranchs from the MPC remain understudied. However, high-throughput sequencing technologies have increased the power and accessibility of genomic tools, even in developing countries such as Mexico. The tools described here provide information relevant for biodiversity conservation. Therefore, we strongly suggest that investment in genomic research will assist implementation of efficient management strategies. In time, this will reduce the extinction risk of the unique elasmobranch biodiversity from the MPC.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.