Substitutions in the Glycogenin-1 Gene Are Associated with the Evolution of Endothermy in Sharks and Tunas

Genomic basis of evolutionary adaptation in a warm-blooded fish

Few fishes have evolved elevated body temperatures compared with ambient temperatures, and only in opah (Lampris spp) is the entire body affected. To understand the molecular basis of endothermy, we analyzed the opah genome and identified 23 genes with convergent amino acid substitutions across fish, birds, and mammals, including slc8b1, which encodes the mitochondrial Na⁺/Ca²⁺ exchanger and is essential for heart function and metabolic heat production. Among endothermic fishes, 44 convergent genes with suggestive metabolic functions were identified, such as glrx3, encoding a crucial protein for hemoglobin maturation. Numerous genes involved in the production and retention of metabolic heat were also found to be under positive selection. Analyses of opah's unique inner-heat-producing pectoral muscle layer (PMI), an evolutionary key innovation, revealed that many proteins were co-opted from dorsal swimming muscles for thermogenesis and oxidative phosphorylation. Thus, the opah genome provides valuable resources and opportunities to uncover the genetic basis of thermal adaptations in fish.

De novo transcriptome assembly, functional annotation and characterization of the Atlantic bluefin tuna (Thunnus thynnus) larval stage

In the present work, we assembled and characterized a de novo larval transcriptome of the Atlantic bluefin tuna Thunnus thynnus by taking advantage of publicly available databases with the goal of better understanding its larval development. The assembled transcriptome comprised 37,117 protein-coding transcripts, of which 13,633 full-length (>80% coverage), with an Ex90N50 of 3061 bp and 76% of complete and single-copy core vertebrate genes orthologues. Of these transcripts, 34,980 had a hit against the EggNOG database and 14,983 with the KEGG database. Codon usage bias was identified in processes such as translation and muscle development. By comparing our data with a set of representative fish species, 87.1% of tuna transcripts were included in orthogroups with other species and 5.1% in assembly-specific orthogroups, which were enriched in terms related to muscle and bone development, visual system and ion transport. Following this comparative approach, protein families related to myosin, extracellular matrix and immune system resulted significantly expanded in the Atlantic bluefin tuna. Altogether, these results provide a glimpse of how the Atlantic bluefin tuna might have achieved early physical advantages over competing species in the pelagic environment. The information generated lays the foundation for future research on the more detailed exploration of physiological responses at the molecular level in different larval stages and paves the way to evolutionary studies on the Atlantic bluefin tuna.

Skeletal muscle and cardiac transcriptomics of a regionally endothermic fish, the Pacific bluefin tuna, Thunnus orientalis

Background

The Pacific bluefin tuna (Thunnus orientalis) is a regionally endothermic fish that maintains temperatures in their swimming musculature, eyes, brain and viscera above that of the ambient water. Within their skeletal muscle, a thermal gradient exists, with deep muscles, close to the backbone, operating at elevated temperatures compared to superficial muscles near the skin. Their heart, by contrast, operates at ambient temperature, which in bluefin tunas can range widely. Cardiac function in tunas reduces in cold waters, yet the heart must continue to supply blood for metabolically demanding endothermic tissues. Physiological studies indicate Pacific bluefin tuna have an elevated cardiac capacity and increased cold-tolerance compared to warm-water tuna species, primarily enabled by increased capacity for sarcoplasmic reticulum calcium cycling within the cardiac muscles.

Results

Here, we compare tissue-specific gene-expression profiles of different cardiac and skeletal muscle tissues in Pacific bluefin tuna. There was little difference in the overall expression of calcium-cycling and cardiac contraction pathways between atrium and ventricle. However, expression of a key sarcoplasmic reticulum calcium-cycling gene, SERCA2b, which plays a key role maintaining intracellular calcium stores, was higher in atrium than ventricle. Expression of genes involved in aerobic metabolism and cardiac contraction were higher in the ventricle than atrium. The two morphologically distinct tissues that derive the ventricle, spongy and compact myocardium, had near-identical levels of gene expression. More genes had higher expression in the cool, superficial muscle than in the warm, deep muscle in both the aerobic red muscle (slow-twitch) and anaerobic white muscle (fast-twitch), suggesting thermal compensation.

Conclusions

We find evidence of widespread transcriptomic differences between the Pacific tuna ventricle and atrium, with potentially higher rates of calcium cycling in the atrium associated with the higher expression of SERCA2b compared to the ventricle. We find no evidence that genes associated with thermogenesis are upregulated in the deep, warm muscle compared to superficial, cool muscle. Heat generation may be enabled by by the high aerobic capacity of bluefin tuna red muscle.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.

Transcriptome-Derived Microsatellites Demonstrate Strong Genetic Differentiation in Pacific White Sharks

Recent advances in genome-scale sequencing technology have allowed the development of high resolution genetic markers for the study of nonmodel taxa. In particular, transcriptome sequencing has proven to be highly useful in generating genomic markers for use in population genetic studies, allowing for insight into species connectivity, as well as local adaptive processes as many transcriptome-derived markers are found within or associated with functional genes. Herein, we developed a set of 30 microsatellite markers from a heart transcriptome for the white shark (Carcharodon carcharias), a widely distributed and globally vulnerable marine predator. Using these markers as well as 10 published anonymous genomic microsatellite loci, we provide 1) the first nuclear genetic assessment of the cross-Pacific connectivity of white sharks, and 2) a comparison of the levels of inferred differentiation across microsatellite marker sets (i.e., transcriptome vs. anonymous) to assess their respective utility to elucidate the population genetic dynamics of white sharks. Significant (FST = 0.083, P = 0.05; G″ST = 0.200; P = 0.001) genetic differentiation was found between Southwestern Pacific (n = 19) and Northeastern Pacific (n = 20) white sharks, indicating restricted, cross Pacific gene flow in this species. Transcriptome-derived microsatellite marker sets identified much higher (up to 2×) levels of genetic differentiation than anonymous genomic markers, underscoring potential utility of transcriptome markers in identifying subtle population genetic differences within highly vagile, globally distributed marine species.Subject areas: Population structure and phylogeography; Conservation genetics and biodiversity.

Review of Current Conservation Genetic Analyses of Northeast Pacific Sharks

Conservation genetics is an applied science that utilizes molecular tools to help solve problems in species conservation and management. It is an interdisciplinary specialty in which scientists apply the study of genetics in conjunction with traditional ecological fieldwork and other techniques to explore molecular variation, population boundaries, and evolutionary relationships with the goal of enabling resource managers to better protect biodiversity and identify unique populations. Several shark species in the northeast Pacific (NEP) have been studied using conservation genetics techniques, which are discussed here. The primary methods employed to study population genetics of sharks have historically been nuclear microsatellites and mitochondrial (mt) DNA. These markers have been used to assess genetic diversity, mating systems, parentage, relatedness, and genetically distinct populations to inform management decisions. Novel approaches in conservation genetics, including next-generation DNA and RNA sequencing, environmental DNA (eDNA), and epigenetics are just beginning to be applied to elasmobranch evolution, physiology, and ecology. Here, we review the methods and results of past studies, explore future directions for shark conservation genetics, and discuss the implications of molecular research and techniques for the long-term management of shark populations in the NEP.