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

Maternal Effects and Trophodynamics Drive Interannual Larval Growth Variability of Atlantic Bluefin Tuna (Thunnus thynnus) from the Gulf of Mexico

Two cohorts of Atlantic bluefin tuna (Thunnus thynnus) larvae were sampled in 2017 and 2018 during the peak of spawning in the Gulf of Mexico (GOM). We examined environmental variables, daily growth, otolith biometry and stable isotopes and found that the GOM18 cohort grew at faster rates, with larger and wider otoliths. Inter and intra-population analyses (deficient vs. optimal growth groups) were carried out for pre- and post-flexion developmental stages to determine maternal and trophodynamic influences on larval growth variability based on larval isotopic signatures, trophic niche sizes and their overlaps. For the pre-flexion stages in both years, the optimal growth groups had significantly lower δ15N, implying a direct relationship between growth potential and maternal inheritance. Optimal growth groups and stages for both years showed lower C:N ratios, reflecting a greater energy investment in growth. The results of this study illustrate the interannual transgenerational trophic plasticity of a spawning stock and its linkages to growth potential of their offsprings in the GOM.

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.

Evidence of isotopic maternal transmission influence on bluefin tuna (Thunnus thynnus) larval growth

Pre-flexion stages of Atlantic bluefin tuna (Thunnus thynnus) larvae were collected in 2014 during the peak of spawning in the two main spawning areas: Gulf of Mexico (GOM) and Mediterranean Sea (MED). We examined daily growth, otolith biometry, and stable isotopes and found that the GOM grew at a faster rate, had larger otoliths, wider daily increments, and significantly lower values of δ15N when compared to the MED. In addition, an intra-population comparative analysis between slow- and fast-growing individuals (deficient vs. optimal growth groups, respectively) showed that optimal growth groups had significantly lower δ15N within each spawning area, implying a direct relationship between growth potential, development, and maternal transmission of isotopic signatures. A third pre-flexion larval group that was aquaculture-reared also exhibited the same pattern to the wild larval groups. In addition, for the first time, we estimated the maternal trophic niches using models developed with field-captured pre-flexion larvae. The estimated maternal trophic niches for the GOM were narrower than the MED, implying differences in the maternal trophodynamics from each nursery area. Overall, the inter-population (GOM vs. MED) and intra-population growth groups (deficient vs. optimal) grew faster and had narrower maternal niches. This study shows the advantages that larval SIA research can aid in the understanding of the trophodynamics of their breeders by examining the trophic relationship of a spawning stock jointly with the development of growth potential in offspring within the same breeding season.