First molecular scombrid timetree (Percomorpha: Scombridae) shows recent radiation of tunas following invasion of pelagic habitat

Genomic resources for the Yellowfin tuna Thunnus albacares

BACKGROUND

The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares.

METHODS AND RESULTS

Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology.

CONCLUSIONS

The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.

A phylogenomic framework for pelagiarian fishes (Acanthomorpha: Percomorpha) highlights mosaic radiation in the open ocean

The fish clade Pelagiaria, which includes tunas as its most famous members, evolved remarkable morphological and ecological variety in a setting not generally considered conducive to diversification: the open ocean. Relationships within Pelagiaria have proven elusive due to short internodes subtending major lineages suggestive of rapid early divergences. Using a novel sequence dataset of over 1000 ultraconserved DNA elements (UCEs) for 94 of the 286 species of Pelagiaria (more than 70% of genera), we provide a time-calibrated phylogeny for this widely distributed clade. Some inferred relationships have clear precedents (e.g. the monophyly of 'core' Stromateoidei, and a clade comprising 'Gempylidae' and Trichiuridae), but others are unexpected despite strong support (e.g. Chiasmodontidae + Tetragonurus). Relaxed molecular clock analysis using node-based fossil calibrations estimates a latest Cretaceous origin for Pelagiaria, with crown-group families restricted to the Cenozoic. Estimated mean speciation rates decline from the origin of the group in the latest Cretaceous, although credible intervals for root and tip rates are broad and overlap in most cases, and there is higher-than-expected partitioning of body shape diversity (measured as fineness ratio) between clades concentrated during the Palaeocene-Eocene. By contrast, more direct measures of ecology show either no substantial deviation from a null model of diversification (diet) or patterns consistent with evolutionary constraint or high rates of recent change (depth habitat). Collectively, these results indicate a mosaic model of diversification. Pelagiarians show high morphological disparity and modest species richness compared to better-studied fish radiations in contrasting environments. However, this pattern is also apparent in other clades in open-ocean or deep-sea habitats, and suggests that comparative study of such groups might provide a more inclusive model of the evolution of diversity in fishes.

Phylotranscriptomic Insights into the Diversification of Endothermic Thunnus Tunas

Birds, mammals, and certain fishes, including tunas, opahs and lamnid sharks, are endothermic, conserving internally generated, metabolic heat to maintain body or tissue temperatures above that of the environment. Bluefin tunas are commercially important fishes worldwide, and some populations are threatened. They are renowned for their endothermy, maintaining elevated temperatures of the oxidative locomotor muscle, viscera, brain and eyes, and occupying cold, productive high-latitude waters. Less cold-tolerant tunas, such as yellowfin tuna, by contrast, remain in warm-temperate to tropical waters year-round, reproducing more rapidly than most temperate bluefin tuna populations, providing resiliency in the face of large-scale industrial fisheries. Despite the importance of these traits to not only fisheries but also habitat utilization and responses to climate change, little is known of the genetic processes underlying the diversification of tunas. In collecting and analyzing sequence data across 29,556 genes, we found that parallel selection on standing genetic variation is associated with the evolution of endothermy in bluefin tunas. This includes two shared substitutions in genes encoding glycerol-3 phosphate dehydrogenase, an enzyme that contributes to thermogenesis in bumblebees and mammals, as well as four genes involved in the Krebs cycle, oxidative phosphorylation, β-oxidation, and superoxide removal. Using phylogenetic techniques, we further illustrate that the eight Thunnus species are genetically distinct, but found evidence of mitochondrial genome introgression across two species. Phylogeny-based metrics highlight conservation needs for some of these species.

Scale diversity in bigeye tuna (Thunnus obesus): Fat-filled trabecular scales made of cellular bone

Tunas of the genus Thunnus possess many morphological and physiological adaptations for their high-performance epipelagic ecology. Although Thunnus anatomy has been studied, there are no quantitative studies on the structure of their scales. We investigated the scales of bigeye tuna (Thunnus obesus) from ten regions of the body using micro computed tomography (µCT)-scanning and histology to quantitatively and qualitatively compare regional scale morphology. We found a diversity of scale sizes and shapes across the body of bigeye tuna and discriminant function analysis on variables derived from µCT-data showed that scales across the body differ quantitatively in shape and size. We also report the discovery of a novel scale type in corselet, tail, and cheek regions. These modified scales are ossified shells supported by internal trabeculae, filled with fat, and possessing an internal blood supply. Histological analysis showed that the outer lamellar layers of these thickened scales are composed of cellular bone, unexpected for a perciform fish in which bone is typically acellular. In the fairing region of the anterior body, these fat-filled scales are stacked in layers up to five scales deep, forming a thickened bony casing. Cheek scales also possess a fat-filled internal trabecular structure, while most posterior body scales are more plate-like and similar to typical teleost scales. While the function of these novel fat-filled scales is unknown, we explore several possible hypotheses for their function.

Non-reef environments impact the diversification of extant jacks, remoras and allies (Carangoidei, Percomorpha)

Various factors may impact the processes of diversification of a clade. In the marine realm, it has been shown that coral reef environments have promoted diversification in various fish groups. With the exception of requiem sharks, all the groups showing a higher level of diversity in reefs than in non-reef habitats have diets based predominantly on plankton, algae or benthic invertebrates. Here we explore the pattern of diversification of carangoid fishes, a clade that includes numerous piscivorous species (e.g. trevallies, jacks and dolphinfishes), using time-calibrated phylogenies as well as ecological and morphological data from both extant and fossil species. The study of carangoid morphospace suggests that reef environments played a role in their early radiation during the Eocene. However, contrary to the hypothesis of a reef-association-promoting effect, we show that habitat shifts to non-reef environments have increased the rates of morphological diversification (i.e. size and body shape) in extant carangoids. Piscivory did not have a major impact on the tempo of diversification of this group. Through the ecological radiation of carangoid fishes, we demonstrate that non-reef environments may sustain and promote processes of diversification of different marine fish groups, at least those including a large proportion of piscivorous species.

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

Despite 400–450 million years of independent evolution, a strong phenotypic convergence has occurred between two groups of fish: tunas and lamnid sharks. This convergence is characterized by centralization of red muscle, a distinctive swimming style (stiffened body powered through tail movements) and elevated body temperature (endothermy). Furthermore, both groups demonstrate elevated white muscle metabolic capacities. All these traits are unusual in fish and more likely evolved to support their fast-swimming, pelagic, predatory behavior. Here, we tested the hypothesis that their convergent evolution was driven by selection on a set of metabolic genes. We sequenced white muscle transcriptomes of six tuna, one mackerel, and three shark species, and supplemented this data set with previously published RNA-seq data. Using 26 species in total (including 7,032 tuna genes plus 1,719 shark genes), we constructed phylogenetic trees and carried out maximum-likelihood analyses of gene selection. We inferred several genes relating to metabolism to be under selection. We also found that the same one gene, glycogenin-1, evolved under positive selection independently in tunas and lamnid sharks, providing evidence of convergent selective pressures at gene level possibly underlying shared physiology.

The historical biogeography of groupers: Clade diversification patterns and processes

Groupers (family Epinephelidae) are a clade of species-rich, biologically diverse reef fishes. Given their ecological variability and widespread distribution across ocean basins, it is important to scrutinize their evolutionary history that underlies present day distributions. This study investigated the patterns and processes by which grouper biodiversity has been generated and what factors have influenced their present day distributions. We reconstructed a robust, time-calibrated molecular phylogeny of Epinephelidae with comprehensive (∼87%) species sampling, whereby diversification rates were estimated and ancestral ranges were reconstructed. Our results indicate that groupers originated in what is now the East Atlantic during the mid-Eocene and diverged successively to form six strongly supported main clades. These clades differ in age (late Oligocene to mid-Miocene), geographic origin (West Atlantic to West Indo-Pacific) and temporal-spatial diversification pattern, ranging from constant rates of diversification to episodes of rapid radiation. Overall, divergence within certain biogeographic regions was most prevalent in groupers, while vicariant divergences were more common in Tropical Atlantic and East Pacific groupers. Our findings reveal that both biological and geographical factors have driven grouper diversification. They also underscore the importance of scrutinizing group-specific patterns to better understand reef fish evolution.

Mitochondrial DNA control region of three mackerels, genus Rastrelliger: structure, molecular diversity and phylogenetic relationship

The complete mitochondrial control regions (CR) of three mackerels (Rastrelliger spp.) were examined and analyzed. The CR contained three domains, in which three termination-associated sequences (TAS-I, TAS-II and TAS-III), two central conserved sequence blocks (CSB-E, CSB-D), three conserved sequence blocks (CSB-I, CSB-II, and CSB-III) and a putative promoter were detected. Molecular indices analyses of the aligned complete CR sequences showed high level of haplotype diversities and genetic divergences among the three species. The intraspecific divergence among species of this genus ranked from 0.25% to 1.62% and interspecific divergence from 1.90% to 4.30%. The phylogenetic tree shows monophyly with R. brachysoma as a basal species of Rastrelliger. Applying the average divergence rate for fish control regions, the results suggest that the time of separation among Rastrelligers could have occurred in the middle Pleistocene era.

Evolutionary origin of the Scombridae (tunas and mackerels): members of a paleogene adaptive radiation with 14 other pelagic fish families

Uncertainties surrounding the evolutionary origin of the epipelagic fish family Scombridae (tunas and mackerels) are symptomatic of the difficulties in resolving suprafamilial relationships within Percomorpha, a hyperdiverse teleost radiation that contains approximately 17,000 species placed in 13 ill-defined orders and 269 families. Here we find that scombrids share a common ancestry with 14 families based on (i) bioinformatic analyses using partial mitochondrial and nuclear gene sequences from all percomorphs deposited in GenBank (10,733 sequences) and (ii) subsequent mitogenomic analysis based on 57 species from those targeted 15 families and 67 outgroup taxa. Morphological heterogeneity among these 15 families is so extraordinary that they have been placed in six different perciform suborders. However, members of the 15 families are either coastal or oceanic pelagic in their ecology with diverse modes of life, suggesting that they represent a previously undetected adaptive radiation in the pelagic realm. Time-calibrated phylogenies imply that scombrids originated from a deep-ocean ancestor and began to radiate after the end-Cretaceous when large predatory epipelagic fishes were selective victims of the Cretaceous-Paleogene mass extinction. We name this clade of open-ocean fishes containing Scombridae “Pelagia” in reference to the common habitat preference that links the 15 families.