Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

First tagging data on large Atlantic bluefin tuna returning to Nordic waters suggest repeated behaviour and skipped spawning.

Atlantic bluefin tuna (Thunnus thynnus; ABFT) is one of the most iconic fish species in the world. Recently, after being very rare for more than half a century, large bluefin tunas have returned to Nordic waters in late summer and autumn, marking the return of the largest predatory fish in Nordic waters. By tagging 18 bluefin tunas with electronic tags (pop-up satellite archival tags), we show that bluefin tuna observed in Nordic waters undertake different migration routes, with individuals migrating into the western Atlantic Ocean, while others stay exclusively in the eastern Atlantic and enter the Mediterranean Sea to spawn. We additionally present evidence of possible skipped spawning inferred from behavioural analyses. In Nordic waters, ABFT are primarily using the upper water column, likely reflecting feeding activity. The results support the hypothesis that ABFT migrating to Nordic waters return to the same general feeding area within the region on an annual basis. These observations may have important implications for management because (1) tunas that come into Nordic waters might represent only a few year classes (as evidenced by a narrow size range), and thus may be particularly vulnerable to area-specific exploitation, and (2) challenge the assumption of consecutive spawning in adult Atlantic bluefin tuna, as used in current stock assessment models. Without careful management and limited exploitation of this part of the ABFT population, the species' return to Nordic waters could be short-lived.

Genetic connectivity between Atlantic bluefin tuna larvae spawned in the Gulf of Mexico and in the Mediterranean Sea.

The highly migratory Atlantic bluefin tuna (ABFT) is currently managed as two distinct stocks, in accordance with natal homing behavior and population structuring despite the absence of barriers to gene flow. Larval fish are valuable biological material for tuna molecular ecology. However, they have hardly been used to decipher the ABFT population structure, although providing the genetic signal from successful breeders. For the first time, cooperative field collection of tuna larvae during 2014 in the main spawning area for each stock, the Gulf of Mexico (GOM) and the Mediterranean Sea (MED), enabled us to assess the ABFT genetic structure in a precise temporal and spatial frame exclusively through larvae. Partitioning of genetic diversity at nuclear microsatellite loci and in the mitochondrial control region in larvae spawned contemporarily resulted in low significant fixation indices supporting connectivity between spawners in the main reproduction area for each population. No structuring was detected within the GOM after segregating nuclear diversity in larvae spawned in two hydrographically distinct regions, the eastern GOM (eGOM) and the western GOM (wGOM), with the larvae from eGOM being more similar to those collected in the MED than the larvae from wGOM. We performed clustering of genetically characterized ABFT larvae through Bayesian analysis and by Discriminant Analysis of Principal Components (DAPC) supporting the existence of favorable areas for mixing of ABFT spawners from Western and Eastern stocks, leading to gene flow and apparent connectivity between weakly structured populations. Our findings suggest that the eastern GOM is more prone for the mixing of breeders from the two ABFT populations. Conservation of this valuable resource exploited for centuries calls for intensification of tuna ichthyoplankton research and standardization of genetic tools for monitoring population dynamics.

Projected effects of ocean warming on an iconic pelagic fish and its fishery.

Increasing sea temperature is a driver of change for many fish traits, particularly for fast-growing epipelagic species with short life spans. With warming, altered spawning phenology and faster growth may produce substantially larger body sizes of the new cohort, affecting fishery productivity. We present an individual-based model (IBM) that predicts the distribution of fish length at catch under observed and projected thermal scenarios, accounting for mortality, temperature-dependent spawning phenology, temperature- and photoperiod- dependent growth. This IBM was demonstrated with Coryphaena hippurus (common dolphinfish), a circumglobally-distributed and highly thermophilic species sustaining commercial and recreational fisheries where it is present. The model projected a 13.2% increase in the average length at catch under marine heatwave conditions compared to the current thermal regime (1995-2005 average). Projections under RCP scenarios 4.5 and 8.5 by the end of the century led to 5.1% and 12.8% increase in average length, respectively. Furthermore, these thermal scenarios affected spawning phenology differently, producing higher variance in body size under RCP 8.5 scenario with respect to marine heatwave conditions. This study highlights how the environmental effects of climate change can alter the distribution of species length at catch.

Accounting for ocean connectivity and hydroclimate in fish recruitment fluctuations within transboundary metapopulations.

Marine resources stewardships are progressively becoming more receptive to an effective incorporation of both ecosystem and environmental complexities into the analytical frameworks of fisheries assessment. Understanding and predicting marine fish production for spatially and demographically complex populations in changing environmental conditions is however still a difficult task. Indeed, fisheries assessment is mostly based on deterministic models that lack realistic parameterizations of the intricate biological and physical processes shaping recruitment, a cornerstone in population dynamics. We use here a large metapopulation of a harvested fish, the European hake (Merluccius merluccius), managed across transnational boundaries in the northwestern Mediterranean, to model fish recruitment dynamics in terms of physics-dependent drivers related to dispersal and survival. The connectivity among nearby subpopulations is evaluated by simulating multi-annual Lagrangian indices of larval retention, imports, and self-recruitment. Along with a proxy of the regional hydroclimate influencing early life stages survival, we then statistically determine the relative contribution of dispersal and hydroclimate for recruitment across contiguous management units. We show that inter-annual variability of recruitment is well reproduced by hydroclimatic influences and synthetic connectivity estimates. Self-recruitment (i.e., the ratio of retained locally produced larvae to the total number of incoming larvae) is the most powerful metric as it integrates the roles of retained local recruits and immigrants from surrounding subpopulations and is able to capture circulation patterns affecting recruitment at the scale of management units. We also reveal that the climatic impact on recruitment is spatially structured at regional scale due to contrasting biophysical processes not related to dispersal. Self-recruitment calculated for each management unit explains between 19% and 32.9% of the variance of recruitment variability, that is much larger than the one explained by spawning stock biomass alone, supporting an increase of consideration of connectivity processes into stocks assessment. By acknowledging the structural and ecological complexity of marine populations, this study provides the scientific basis to link spatial management and temporal assessment within large marine metapopulations. Our results suggest that fisheries management could be improved by combining information of physical oceanography (from observing systems and operational models), opening new opportunities such as the development of short-term projections and dynamic spatial management.

Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next-generation sequencing.

The Atlantic bluefin tuna is a highly migratory species emblematic of the challenges associated with shared fisheries management. In an effort to resolve the species' stock dynamics, a genomewide search for spatially informative single nucleotide polymorphisms (SNPs) was undertaken, by way of sequencing reduced representation libraries. An allele frequency approach to SNP discovery was used, combining the data of 555 larvae and young-of-the-year (LYOY) into pools representing major geographical areas and mapping against a newly assembled genomic reference. From a set of 184,895 candidate loci, 384 were selected for validation using 167 LYOY. A highly discriminatory genotyping panel of 95 SNPs was ultimately developed by selecting loci with the most pronounced differences between western Atlantic and Mediterranean Sea LYOY. The panel was evaluated by genotyping a different set of LYOY (n = 326), and from these, 77.8% and 82.1% were correctly assigned to western Atlantic and Mediterranean Sea origins, respectively. The panel revealed temporally persistent differentiation among LYOY from the western Atlantic and Mediterranean Sea (FST  = 0.008, p = .034). The composition of six mixed feeding aggregations in the Atlantic Ocean and Mediterranean Sea was characterized using genotypes from medium (n = 184) and large (n = 48) adults, applying population assignment and mixture analyses. The results provide evidence of persistent population structuring across broad geographic areas and extensive mixing in the Atlantic Ocean, particularly in the mid-Atlantic Bight and Gulf of St. Lawrence. The genomic reference and genotyping tools presented here constitute novel resources useful for future research and conservation efforts.

Ecology of the Atlantic black skipjack Euthynnus alletteratus (Osteichthyes: Scombridae) in the western Mediterranean Sea inferred by parasitological analysis.

Between 2008 and 2011, the head of 150 Euthynnus alletteratus (Osteichthyes: Scombridae) caught inshore off the southeastern Iberian coast (western Mediterranean Sea) were examined for parasites. Two monogeneans, four didymozoid trematodes and four copepods were found. Parasite abundance showed a positive relationship with the annual sea surface temperature, except for Pseudocycnus appendiculatus, but negative with the sea depth (Capsala manteri, Neonematobothrium cf. kawakawa and Caligus bonito). Prevalences and mean abundances differed significantly among sampling areas, except for C. manteri, Oesophagocystis sp. 2 and Ceratocolax euthynni, and sampling years (Melanocystis cf. kawakawa, N.cf. kawakawa, P. appendiculatus and Unicolax collateralis). Results indicate that the parasite abundances of E. alletteratus in the western Mediterranean Sea depend mainly on regional environmental variables, which can show interannual variations. The presence of pelagic parasites, i.e. didymozoids and P. appendiculatus, could indicate that E. alletteratus migrates between inshore and offshore pelagic domains. The different parasite faunas reported in E. alletteratus populations from the western Atlantic Ocean and the Mediterranean Sea appear to point out the geographical host isolation. These results suggest that E. alletteratus inhabiting the western Mediterranean Sea performs inshore-offshore small-scale migrations, and not transoceanic migrations between the western Atlantic Ocean and Mediterranean Sea.

Molecular Identification of Atlantic Bluefin Tuna (Thunnus thynnus, Scombridae) Larvae and Development of a DNA Character-Based Identification Key for Mediterranean Scombrids.

The Atlantic bluefin tuna, Thunnus thynnus, is a commercially important species that has been severely over-exploited in the recent past. Although the eastern Atlantic and Mediterranean stock is now showing signs of recovery, its current status remains very uncertain and as a consequence their recovery is dependent upon severe management informed by rigorous scientific research. Monitoring of early life history stages can inform decision makers about the health of the species based upon recruitment and survival rates. Misidentification of fish larvae and eggs can lead to inaccurate estimates of stock biomass and productivity which can trigger demands for increased quotas and unsound management conclusions. Herein we used a molecular approach employing mitochondrial and nuclear genes (CO1 and ITS1, respectively) to identify larvae (n = 188) collected from three spawning areas in the Mediterranean Sea by different institutions working with a regional fisheries management organization. Several techniques were used to analyze the genetic sequences (sequence alignments using search algorithms, neighbour joining trees, and a genetic character-based identification key) and an extensive comparison of the results is presented. During this process various inaccuracies in related publications and online databases were uncovered. Our results reveal important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. While less than half of larvae provided were bluefin tuna, other dominant taxa were bullet tuna (Auxis rochei), albacore (Thunnus alalunga) and little tunny (Euthynnus alletteratus). We advocate an expansion of expertise for a new generation of morphology-based taxonomists, increased dialogue between morphology-based and molecular taxonomists and increased scrutiny of public sequence databases.

Trophic ecology of Atlantic Bluefin Tuna (Thunnus thynnus) [corrected] larvae from the Gulf of Mexico and NW Mediterranean spawning grounds: A Comparative Stable Isotope Study.

The present study uses stable isotopes of nitrogen and carbon (δ15Nandδ13C) as trophic indicators for Atlantic bluefin tuna larvae (BFT) (6-10 mm standard length) in the highly contrasting environmental conditions of the Gulf of Mexico (GOM) and the Balearic Sea (MED). These regions are differentiated by their temperature regime and relative productivity, with the GOM being significantly warmer and more productive. MED BFT larvae showed the highest δ15N signatures, implying an elevated trophic position above the underlying microzooplankton baseline. Ontogenetic dietary shifts were observed in the BFT larvae from the GOM and MED which indicates early life trophodynamics differences between these spawning habitats. Significant trophic differences between the GOM and MED larvae were observed in relation to δ15N signatures in favour of the MED larvae, which may have important implications in their growth during their early life stages.These low δ15N levels in the zooplankton from the GOM may be an indication of a shifting isotopic baseline in pelagic food webs due to diatrophic inputs by cyanobacteria. Lack of enrichment for δ15N in BFT larvae compared to zooplankton implies an alternative grazing pathway from the traditional food chain of phytoplankton-zooplankton-larval fish. Results provide insight for a comparative characterization of the trophic pathways variability of the two main spawning grounds for BFT larvae.

Spatial scale, means and gradients of hydrographic variables define pelagic seascapes of bluefin and bullet tuna spawning distribution.

Seascape ecology is an emerging discipline focused on understanding how features of the marine habitat influence the spatial distribution of marine species. However, there is still a gap in the development of concepts and techniques for its application in the marine pelagic realm, where there are no clear boundaries delimitating habitats. Here we demonstrate that pelagic seascape metrics defined as a combination of hydrographic variables and their spatial gradients calculated at an appropriate spatial scale, improve our ability to model pelagic fish distribution. We apply the analysis to study the spawning locations of two tuna species: Atlantic bluefin and bullet tuna. These two species represent a gradient in life history strategies. Bluefin tuna has a large body size and is a long-distant migrant, while bullet tuna has a small body size and lives year-round in coastal waters within the Mediterranean Sea. The results show that the models performance incorporating the proposed seascape metrics increases significantly when compared with models that do not consider these metrics. This improvement is more important for Atlantic bluefin, whose spawning ecology is dependent on the local oceanographic scenario, than it is for bullet tuna, which is less influenced by the hydrographic conditions. Our study advances our understanding of how species perceive their habitat and confirms that the spatial scale at which the seascape metrics provide information is related to the spawning ecology and life history strategy of each species.

Metazoan parasites on the gills of the skipjack tuna Katsuwonus pelamis (Osteichthyes: Scombridae) from the Alboran Sea (western Mediterranean Sea).

The gills of 31 skipjack tuna Katsuwonus pelamis (L.) caught in the Alboran Sea (western Mediterranean Sea) were examined for metazoan parasites, and the gills of 4 specimens from the Balearic Sea (also western Mediterranean Sea) were analysed for comparative purposes. Nine -species of parasites were found, including 8 didymozoid trematodes (Atalostrophion cf. bio-varium, Didymocylindrus filiformis, Didymocylindrus simplex, Didymocystis reniformis, Didymoproblema fusiforme, Didymozoon longicolle, Koellikeria sp. and Lobatozoum multisacculatum) and 1 caligid copepod (Caligus bonito). Koellikeria sp. and L. multisacculatum were not recorded in the Balearic Sea. Most of the parasites (79.2% of all specimens) were didymozoids. Didymozoon longicolle was the dominant species; A. cf. biovarium, D. simplex, D. fusiforme and L. multisacculatum are reported from the Mediterranean Sea for the first time. No correlation was found between the intensity of infection of any parasite species and host size or sex. Most of the parasites, particularly didymozoids, showed a high site-specificity. Significant differences were found between the parasite assemblages of K. pelamis from the Alboran Sea and from the Atlantic Ocean. D. fusiforme, D. longicolle and L. multisacculatum are suggested as potential tags to follow skipjack tuna migrations between the Atlantic Ocean and Mediterranean Sea.