Phylogenetic relationships between tuna species of the genus Thunnus (Scombridae: Teleostei): Inconsistent implications from morphology,nuclear and mitochondrial genomes

In order to infer phylogenetic relationships between tuna species of the genus Thunnus, partial sequences of the mitochondrial cytochrome b and ATPase genes were determined in all eight species. Supplemental restriction analysis on the nuclear rRNA gene was also carried out. Pacific northern bluefin tuna (Thunnus thynnus orientalis) was found to have mtDNA distinct from that of the Atlantic subspecies (T. t. thynnus) but very similar to that from the species albacore (T. alaluga). In contrast, no differentiation in nuclear genome was observed between the Atlantic and Pacific northern bluefin tunas. The Atlantic northern bluefin and southern bluefin tunas possessed mtDNA sequences very similar to species of yellowfin tuna group and not so similar to albacore and bigeye tunas which were morphologically assigned to the bluefin tuna group. The molecular data indicate that (1) mtDNA from albacore has been incorporated into the Pacific population of northern bluefin tuna and has extensively displaced the original mtDNA, and (2) albacore is the earliest offshoot, followed by bigeye tuna in this genus, which is inconsistent with the phylogenetic relationships between these tuna species inferred from morphology. Correspondence to: S. Chow     

Population genetics of Atlantic bluefin tuna,Thunnus thynnus (Linnaeus, 1758), in the Mediterranean: implications for its conservation management

The movement of Atlantic bluefin tuna (ABFT) across international boundaries necessitates traceability strategies that would provide more accurate information needed for stock assessment. The Mediterranean Sea is one of the main contributors to ABFT reproduction and global population genetic diversity. In the present study this genetic variability was investigated using 193 samples of adult bluefin tuna from Spain, Turkey and Malta – a longitudinal distance of 3400 km. Analysed were 13 microsatellite loci (eight of which were newly‐tested) as genetic markers for the population study. Allele richness measured per locus and sampling location varied from 1.89 to 8.88, taking into account rarefaction. ABFT private alleles were detected in each of the three sampling sites. No significant spatial genetic divergence was found between pairs at the studied locations (F_ST values <0.0001; P‐values >0.05). Bayesian clustering analysis corroborated a single and panmictic ABFT population in the Mediterranean Sea. Statistical power analyses indicated a high probability of detecting genetic differentiation and population structure with the sample size and microsatellites used, even at an F_ST value of 0.005. From the results it may be postulated that migrating ABFT during the spawning season are allowing gene flow within the Mediterranean Sea. The complex interplay of movements, including plasticity in the selection of spawning sites with increasing age and environmental conditions, require multiple and new fisheries monitoring and management techniques in order to target the ABFT long‐term conservation effectively.     

Allozyme variation in Atlantic salmon from the British Isles: associations with geography and the environment

Data on geographical variation in allele frequencies at enzyme coding loci in Atlantic salmon from the British Isles were collated from published and unpublished sources. Statistically significant differences in allele frequencies were found among samples both within and among river systems, suggesting that the Atlantic salmon in the British Isles is not a panmictic population and that even within major river systems it cannot be treated as a single genetic stock for fisheries management purposes. Although there was some evidence of regional differences in the frequency of some rare alleles, most single‐locus variation did not show strong geographic patterns, with the exception of the AAT‐4 * locus at which allele frequencies had a significant latitudinal cline. There was some evidence for the existence of genetically‐distinct celtic and boreal races of Atlantic salmon in the British Isles as previously has been suggested. Multiple regression analyses revealed associations between genetic variation and local environmental conditions ( i.e . between variation at MEP‐2 * and both temperature and local river gradient), providing additional evidence for adaptive population divergence in the species.     

Comparative phylogeography of Atlantic bluefin tuna and swordfish: the combined effects of vicariance, secondary contact, introgression, and population expansion on the regional phylogenies of two highly migratory pelagic fishes

Comparative phylogeography has revealed remarkable patterns of concordance in the maternal phylogenies of many species. The phylogeography and historical demography of the mitochondrial control region I for 607 Atlantic bluefin tuna (Thunnus thynnus) and 275 swordfish (Xiphias gladius) were analyzed to clarify the complex phylogenetic signals in the North Atlantic-Mediterranean region where they are sympatric. Atlantic bluefin tuna mtDNA is polyphyletic, and includes rare sequences sister to Pacific bluefin tuna (Thunnus orientalis) and introgressed albacore (Thunnus alalunga) sequences. There is no geographic partitioning between Atlantic and Mediterranean samples of Atlantic bluefin tuna (Phi(ST)=0.002). In contrast, Atlantic and Mediterranean swordfish are differentiated (Phi(ST)=0.091) due to the combined effects of vicariance, secondary contact, and dissimilar regional demographic histories. Mediterranean swordfish has substantially less variation, and a more recent history (tau=2.42) than that of Atlantic swordfish (tau=7.02). In spite of the discordant phylogenetic and phylogeographic signals, the demographic history of Atlantic swordfish and Atlantic bluefin tuna (tau=7.51) suggests concordance in the timeline of population expansion. Possible scenarios of cladogenesis, expansion, and contraction, influenced by glacial cycles during the Pleistocene, are formulated.     

Fast versus slow growing tuna species: age,growth, and implications for population dynamics and fisheries management

Growth models describe the change in length or weight as a function of age. Growth curves in tunas can take different forms from relatively simple von Bertalanffy growth curves (Atlantic bluefin, albacore tunas) to more complex two- or three-stanza growth curves (yellowfin, bigeye, skipjack, southern bluefin tunas). We reviewed the growth of the principal market tunas (albacore, bigeye, skipjack, yellowfin and the three bluefin tuna species) in all oceans to ascertain the different growth rates among tuna species and their implications for population productivity and resilience. Tunas are among the fastest-growing of all fishes. Compared to other species, tunas exhibit rapid growth (i.e., relatively high K) and achieve large body sizes (i.e., high L _∞ ). A comparison of their growth functions reveals that tunas have evolved different growth strategies. Tunas attain asymptotic sizes (L _∞ ), ranging from 75 cm FL (skipjack tuna) to 400 cm FL (Atlantic bluefin tuna), and reach L _∞ at different rates (K), varying from 0.95 year^?1 (skipjack tuna) to 0.05 year^?1 (Atlantic bluefin tuna). Skipjack tuna (followed by yellowfin tuna) is considered the “fastest growing” species of all tunas. Growth characteristics have important implications for population dynamics and fisheries management outcomes since tunas, and other fish species, with faster growth rates generally support higher estimates of Maximum Sustainable Yield (MSY) than species with slower growth rates.     

Genetic identity of YOY bluefin tuna from the eastern and western Atlantic spawning areas

We used 320 young-of-the-year (YOY) specimens of the highly migratory and overfished Atlantic bluefin tuna, Thunnus thynnus, Linnaeus 1758, to evaluate the hypothesis that Atlantic bluefin tuna comprises 2 stocks with spawning grounds in the Gulf of Mexico and in the Mediterranean Sea. Significant genetic differentiation at 8 nuclear microsatellite loci (F(ST) = 0.0059, P = 0.0005) and at the mitochondrial control region (Phi(ST) = 0.0129, P = 0.0139) was detected among YOY Atlantic bluefin tuna captured on spawning grounds in the Gulf of Mexico (n = 40) versus the western (n = 255) and eastern (n = 25) basins of the Mediterranean Sea. The genetic divergence among spawning populations, combined with the extensive trans-Atlantic movements reported for juvenile and adult Atlantic bluefin tuna, indicates a high degree of spawning site fidelity. Recognition of genetically distinct populations necessitates independent management of Atlantic bluefin tuna on spawning grounds and warrants evaluation of the level of mixing of populations on feeding grounds. The genetic pattern might not have been detected unless juvenile specimens or actively spawning adults had been sampled.     

Genetic structuring and migration patterns of Atlantic bigeye tuna, <Emphasis Type="Italic">Thunnus obesus</Emphasis>(Lowe, 1839)

Background  

Large pelagic fishes are generally thought to have little population genetic structuring based on their cosmopolitan distribution, large population sizes and high dispersal capacities. However, gene flow can be influenced by ecological (e.g. homing behaviour) and physical (e.g. present-day ocean currents, past changes in sea temperature and levels) factors. In this regard, Atlantic bigeye tuna shows an interesting genetic structuring pattern with two highly divergent mitochondrial clades (Clades I and II), which are assumed to have been originated during the last Pleistocene glacial maxima. We assess genetic structure patterns of Atlantic bigeye tuna at the nuclear level, and compare them with mitochondrial evidence.     

Isolation and characterization of seven tetranucleotide microsatellite loci from Atlantic northern bluefin tuna Thunnus thynnus thynnus

Microsatellite‐enriched genomic libraries were obtained from the Atlantic northern bluefin tuna, Thunnus thynnus thynnus and seven tetranucleotide markers were successfully isolated and characterized from this library. These markers were found to have between 1 and 17 alleles in Atlantic northern bluefin tuna and heterozygosity ranged from 0 to 0.85. No deviations from the expectation of Hardy‐Weinburg equilibrium were found for any marker. Several of these markers amplify reliably in other tuna species.     

The genetic population structure of Thunnus thynnus (Linnaeus, 1758) in the Mediterranean Sea,a controversial issue

The Atlantic bluefin tuna (ABFT), Thunnus thynnus (Linnaeus, 1758), is an important commercial species managed as two different stocks, western and eastern Atlantic, with their spawning grounds in the Gulf of Mexico and in the Mediterranean Sea, respectively. The eastern Atlantic stock has been overexploited in the last decades, leading to the application of specific management measures introduced by the International Commission for the Atlantic Tuna (ICCAT). A clear understanding of the genetic structure of ABFT Mediterranean population should be pursued in order to support management decisions. To date the genetic studies on the Mediterranean ABFT, carried out with different molecular markers and sampling procedures, have produced unclear results. Here, we analysed ABFT samples from central and western Mediterranean Sea with mitochondrial sequences and 11 microsatellite loci to investigate, among the others, the area of the Strait of Messina, where environmental conditions seem to support a resident population of ABFT. Furthermore, genetic analyses of mitochondrial sequences were carried out including nucleotide sequences of Adriatic ABFT wild larvae retrieved from GenBank. Among the investigated areas a genetic differentiation was detected between the Strait of Messina and the Tyrrhenian Sea with microsatellite loci according to the exact G test, but not to the Bayesian analyses carried out with STRUCTURE. The analyses with mitochondrial sequences do not reveal any differentiation among sampled areas, however, a highly significant genetic divergence was observed between the Adriatic mitochondrial sequences retrieved from GenBank and the central‐western Mediterranean sequences obtained in the present work. Our results provide some evidence of population structure of Mediterranean ABFT adding pieces to a still unclear picture.     

Microsatellite and mitochondrial DNA analyses of Atlantic bluefin tuna (Thunnus thynnus thynnus) population structure in the Mediterranean Sea

Genetic variation was surveyed at nine microsatellite loci and the mitochondrial control region (868 bp) to test for the presence of genetic stock structure in young-of-the-year Atlantic bluefin tuna (Thunnus thynnus thynnus) from the Mediterranean Sea. Bluefin tuna were sampled over a period of 5 years from the Balearic and Tyrrhenian seas in the western basin of the Mediterranean Sea, and from the southern Ionian Sea in the eastern basin of the Mediterranean Sea. Analyses of multilocus microsatellite genotypes and mitochondrial control region sequences revealed no significant heterogeneity among collections taken from the same location in different years; however, significant spatial genetic heterogeneity was observed across all samples for both microsatellite markers and mitochondrial control region sequences (FST=0.0023, P=0.038 and PhiST=0.0233, P=0.000, respectively). Significant genetic differentiation between the Tyrrhenian and Ionian collections was found for both microsatellite and mitochondrial markers (FST=0.0087, P=0.015 and PhiST=0.0367, P=0.030, respectively). These results suggest the possibility of a genetically discrete population in the eastern basin of the Mediterranean Sea.     

Seasonal Movements,Aggregations and Diving Behavior of Atlantic Bluefin Tuna (Thunnus thynnus) Revealed with Archival Tags

Electronic tags were used to examine the seasonal movements, aggregations and diving behaviors of Atlantic bluefin tuna (Thunnus thynnus) to better understand their migration ecology and oceanic habitat utilization. Implantable archival tags (n = 561) were deployed in bluefin tuna from 1996 to 2005 and 106 tags were recovered. Movement paths of the fish were reconstructed using light level and sea-surface-temperature-based geolocation estimates. To quantify habitat utilization we employed a weighted kernel estimation technique that removed the biases of deployment location and track length. Throughout the North Atlantic, high residence times (167±33 days) were identified in four spatially confined regions on a seasonal scale. Within each region, bluefin tuna experienced distinct temperature regimes and displayed different diving behaviors. The mean diving depths within the high-use areas were significantly shallower and the dive frequency and the variance in internal temperature significantly higher than during transit movements between the high-use areas. Residence time in the more northern latitude high-use areas was significantly correlated with levels of primary productivity. The regions of aggregation are associated with areas of abundant prey and potentially represent critical foraging habitats that have seasonally abundant prey. Throughout the North Atlantic mean diving depth was significantly correlated with the depth of the thermocline, and dive behavior changed in relation to the stratification of the water column. In this study, with numerous multi-year tracks, there appear to be repeatable patterns of clear aggregation areas that potentially are changing with environmental conditions. The high concentrations of bluefin tuna in predictable locations indicate that Atlantic bluefin tuna are vulnerable to concentrated fishing efforts in the regions of foraging aggregations.     

Dispersal routes and habitat utilization of juvenile Atlantic bluefin tuna, Thunnus thynnus, tracked with mini PSAT and archival tags

Between 2005 and 2009, we deployed 58 miniature pop-up satellite archival tags (PSAT) and 132 implanted archival tags on juvenile Atlantic bluefin tuna (age 2-5) in the northwest Atlantic Ocean. Data returned from these efforts (n?=?26 PSATs, 1 archival tag) revealed their dispersal routes, horizontal and vertical movements and habitat utilization. All of the tagged bluefin tuna remained in the northwest Atlantic for the duration observed, and in summer months exhibited core-use of coastal seas extending from Maryland to Cape Cod, MA, (USA) out to the shelf break. Their winter distributions were more spatially disaggregated, ranging south to the South Atlantic Bight, northern Bahamas and Gulf Stream. Vertical habitat patterns showed that juvenile bluefin tuna mainly occupied shallow depths (mean=?5-12 m, sd?=?15-23.7 m) and relatively warm water masses in summer (mean=?17.9-20.9°C, sd=?4.2-2.6°C) and had deeper and more variable depth patterns in winter (mean=?41-58 m, sd=?48.9-62.2 m). Our tagging results reveal annual dispersal patterns, behavior and oceanographic associations of juvenile Atlantic bluefin tuna that were only surmised in earlier studies. Fishery independent profiling from electronic tagging also provide spatially and temporally explicit information for evaluating dispersals rates, population structure and fisheries catch patterns.     

Contrasting levels of genetic differentiation among populations of wolverines (<Emphasis Type="Italic">Gulo gulo</Emphasis>) from northern Canada revealed by nuclear and mitochondrial loci

Habitat loss, fragmentation, overharvest, and other anthropogenic factors have resulted in population and distribution declines for North American wolverines (Gulo gulo). Currently, wolverines east of the Hudson Bay are endangered and possibly extinct, whereas the status of wolverines throughout the remaining Holarctic is vulnerable. Three previous studies using nuclear loci have detected little to no significant structuring among wolverines sampled across northern Canada. Based on these results it has been suggested that wolverines in northern Canada represent a single, panmictic population. However, as has been shown in numerous studies, in cases of female site fidelity, it is possible to have demographically autonomous populations even with male-biased gene flow. To better assess the genetic structure of wolverines in northern Canada, we examined nine microsatellite loci and DNA sequence variation from a 200 bp fragment of the mitochondrial (mtDNA) control region for 270 wolverines from nine collecting areas representing three regions of northern Canada. In agreement with previous studies, microsatellite analyses revealed a lack of significant population substructure (F _ST=0.0004). However, analysis of molecular variance, comparisons of pairwise F _ST values and nested-clade analysis of the mtDNA data revealed considerable genetic structuring among samples of wolverines from these three regions of northern Canada. These mitochondrial data provide evidence that wolverines in Canada are genetically structured due to female philopatry. The contrasting patterns of genetic differentiation based on nuclear and mitochondrial data highlight the importance of examining both nuclear and mitochondrial loci when attempting to elucidate patterns of genetic structure.     

Genetic differentiation among North Atlantic killer whale populations

Population genetic structure of North Atlantic killer whale samples was resolved from differences in allele frequencies of 17 microsatellite loci, mtDNA control region haplotype frequencies and for a subset of samples, using complete mitogenome sequences. Three significantly differentiated populations were identified. Differentiation based on microsatellite allele frequencies was greater between the two allopatric populations than between the two pairs of partially sympatric populations. Spatial clustering of individuals within each of these populations overlaps with the distribution of particular prey resources: herring, mackerel and tuna, which each population has been seen predating. Phylogenetic analyses using complete mitogenomes suggested two populations could have resulted from single founding events and subsequent matrilineal expansion. The third population, which was sampled at lower latitudes and lower density, consisted of maternal lineages from three highly divergent clades. Pairwise population differentiation was greater for estimates based on mtDNA control region haplotype frequencies than for estimates based on microsatellite allele frequencies, and there were no mitogenome haplotypes shared among populations. This suggests low or no female migration and that gene flow was primarily male mediated when populations spatially and temporally overlap. These results demonstrate that genetic differentiation can arise through resource specialization in the absence of physical barriers to gene flow.     

Phylogenetic relationships among Thunnus species inferred from rDNA ITS1 sequence

Intra and interspecific nucleotide sequence variation of rDNA first internal transcribed spacer (ITS1) was analysed using all eight species of the genus Thunnus plus two out‐group species within the same family, skipjack tuna Katsuwonus pelamis and striped bonito Sarda orientalis . Intraspecific nucleotide sequence variation in ITS1, including intra‐genomic variation, was low, ranging from 0·003 to 0·014 [Kimura's two parameter distance (K2P)], whereas variation between species within the genus Thunnus ranged from 0·009 to 0·05. The Atlantic and Pacific northern bluefin tunas Thunnus thynnus thynnus and Thunnus thynnus orientalis , recently proposed to be distinct species, were found to share nearly identical ITS1 sequences (mean K2P = 0·006) well within the range of intraspecific variation. The northern bluefin tuna appeared to be a sister group to albacore Thunnus alalunga , with all other Thunnus species in a distinct clade. The ITS1 phylogeny was consistent with mtDNA phylogeny in clustering the three tropical Thunnus species ( T. albacares , T. atlanticus and T. tonggol ). Southern bluefin Thunnus maccoyii and bigeye Thunnus obesus tunas showed a closer affinity to this tropical tuna group than to the northern bluefin tuna and albacore. The molecular data supported mitochondrial introgression between species and contradicted morphological subdivision of the genus into two subgenera Neothunnus and Thunnus .     

Spawning Behaviour and Post-Spawning Migration Patterns of Atlantic Bluefin Tuna (Thunnus thynnus) Ascertained from Satellite Archival Tags

Spawning behaviour of Atlantic bluefin tuna (Thunnus thynnus) was investigated using electronic satellite tags deployed in the western Mediterranean spawning ground, around the Balearic Islands (years 2009-2011). All the fish were tagged underwater and released within schools. In general, the fish tagged in the same year/school displayed common migratory trends. Following extended residency around the Balearic Islands, most tagged tuna crossed the Strait of Gibraltar heading for the North Atlantic. Discrepancies between the migratory tracks reconstructed from this and previous electronic tagging studies suggest that the bluefin tuna Mediterranean population may comprise distinct units exhibiting differing migratory behaviours. The diving behaviour varied between oceanic regions throughout the migratory pathways, the shallowest distribution taking place in the spawning ground and the deepest at the Strait of Gibraltar. A unique diving pattern was found on the majority of nights while the fish stayed at the spawning ground; it consisted of frequent and brief oscillatory movements up and down through the mixed layer, resulting in thermal profiles characterized by oscillations about the thermocline. Such a pattern is believed to reflect recent courtship and spawning activity. Reproductive parameters inferred from the analysis of vertical profiles are consistent with those estimated in previous studies based on biological samples.     

Mitochondrial DNA analysis reveals three stocks of yellowfin tuna Thunnus albacares (Bonnaterre, 1788) in Indian waters

Yellowfin tuna (Thunnus albacares) is an epipelagic, oceanic species of family Scombridae found in tropical and subtropical region of Pacific, Indian and Atlantic Ocean. It is commercially important fish and accounts for 19 % of total tuna catches in Indian waters. In present study, population structure of yellowfin tuna was examined using sequence analysis of mitochondrial DNA from seven geographically distinct locations along the Indian coast. A 500 bp segment of D-loop region was sequenced and analysed for 321 yellowfin samples. Hierarchical analysis of molecular variance showed significant genetic differentiation among three groups (VE); (AG); (KO, TU, PO, VI, PB) analyzed (Φ _ST  = 0.03844, P ≤ 0.001). In addition, spatial analysis of molecular variance identified three genetically heterogeneous groups of yellowfin tuna in Indian waters. Results were further corroborated by significant value of nearest neighbour statistic (S _nn = 0.261, P ≤ 0.001). Thus finding of this study rejects the null hypothesis of single panmictic population of yellowfin tuna in Indian waters.     

Genetic diversity and historical demography of Atlantic bigeye tuna (Thunnus obesus)

Bigeye (Thunnus obesus) is a large, pelagic, and migratory species of tuna that inhabits tropical and temperate marine waters worldwide. Previous studies based on mitochondrial RFLP data have shown that bigeye tunas from the Atlantic Ocean are the most interesting from a genetic point of view. Two highly divergent mitochondrial haplotype clades (I and II) coexist in the Atlantic Ocean. One is almost exclusive of the Atlantic Ocean whereas the other is also found in the Indo-Pacific Ocean. Bigeye tuna from the Atlantic Ocean is currently managed as a single stock, although this assumption remains untested at the genetic level. Therefore, genetic diversity was determined at the mitochondrial control region to test the null hypothesis of no population structure in bigeye tuna from the Atlantic Ocean. A total of 331 specimens were sampled from four locations in the Atlantic Ocean (Canada, Azores, Canary Islands, and Gulf of Guinea), and one in the Indian and Pacific Oceans, respectively. The reconstructed neighbor-joining phylogeny confirmed the presence of Clades I and II throughout the Atlantic Ocean. No apparent latitudinal gradient of the proportions of both clades in the different collection sites was observed. Hierarchical AMOVA tests and pairwise phi(ST) comparisons involving Atlantic Ocean Clades I and II were consistent with a single stock of bigeye tuna in the Atlantic Ocean. Population genetic analyses considering phylogroups independently supported gene flow within Clade II throughout the Atlantic Ocean, and within Clade I between Atlantic and Indo-Pacific Oceans. The latter result suggests present uni-directional gene flow from the Indo-Pacific into the Atlantic Ocean. Moreover, mismatch analyses dated divergence of Clades I and II during the Pleistocene, as previously proposed. In addition, migration rates were estimated using coalescent methods, and showed a net migration from Atlantic Ocean feeding grounds towards the Gulf of Guinea, the best-known spawning ground of Atlantic bigeye tuna.     

Implications of isolation and low genetic diversity in peripheral populations of an amphi-Atlantic coral

Limited dispersal and connectivity in marine organisms can have negative fitness effects in populations that are small and isolated, but reduced genetic exchange may also promote the potential for local adaptation. Here, we compare the levels of genetic diversity and connectivity in the coral Montastraea cavernosa among both central and peripheral populations throughout its range in the Atlantic. Genetic data from one mitochondrial and two nuclear loci in 191 individuals show that M. cavernosa is subdivided into three genetically distinct regions in the Atlantic: Caribbean-North Atlantic, Western South Atlantic (Brazil) and Eastern Tropical Atlantic (West Africa). Within each region, populations have similar allele frequencies and levels of genetic diversity; indeed, no significant differentiation was found between populations separated by as much as 3000 km, suggesting that this coral species has the ability to disperse over large distances. Gene flow within regions does not, however, translate into connectivity across the entire Atlantic. Instead, substantial differences in allele frequencies across regions suggest that genetic exchange is infrequent between the Caribbean, Brazil and West Africa. Furthermore, markedly lower levels of genetic diversity are observed in the Brazilian and West African populations. Genetic diversity and connectivity may contribute to the resilience of a coral population to disturbance. Isolated peripheral populations may be more vulnerable to human impacts, disease or climate change relative to those in the genetically diverse Caribbean-North Atlantic region.     

Electronic Tagging of Atlantic Bluefin Tuna (Thunnus thynnus,L.) Reveals Habitat Use and Behaviors in the Mediterranean Sea

We analyzed the movements of Atlantic tuna (Thunnus thynnus L.) in the Mediterranean Sea using data from 2 archival tags and 37 pop-up satellite archival tags (PAT). Bluefin tuna ranging in size from 12 to 248 kg were tagged on board recreational boats in the western Mediterranean and the Adriatic Sea between May and September during two different periods (2000 to 2001 and 2008 to 2012). Although tuna migrations between the Mediterranean Sea and the Atlantic Ocean have been well reported, our results indicate that part of the bluefin tuna population remains in the Mediterranean basin for much of the year, revealing a more complex population structure. In this study we demonstrate links between the western Mediterranean, the Adriatic and the Gulf of Sidra (Libya) using over 4336 recorded days of location and behavior data from tagged bluefin tuna with a maximum track length of 394 days. We described the oceanographic preferences and horizontal behaviors during the spawning season for 4 adult bluefin tuna. We also analyzed the time series data that reveals the vertical behavior of one pop-up satellite tag recovered, which was attached to a 43.9 kg tuna. This fish displayed a unique diving pattern within 16 days of the spawning season, suggesting a use of the thermocline as a thermoregulatory mechanism compatible with spawning. The results obtained hereby confirm that the Mediterranean is clearly an important habitat for this species, not only as spawning ground, but also as an overwintering foraging ground.