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.     

Consequences of the historical demography on the global population structure of two highly migratory cosmopolitan marine fishes: the yellowfin tuna (<Emphasis Type="Italic">Thunnus albacares</Emphasis>) and the skipjack tuna (<Emphasis Type="Italic">Katsuwonus pelamis</Emphasis>)

Background  

Yellowfin and skipjack tuna are globally distributed in the world's tropical and sub-tropical oceans. Since little, if any, migration of these fishes occurs between the Atlantic and Indo-Pacific Oceans, one might expect to see genetic differences between sub-populations in these ocean basins. However, yellowfin and skipjack tuna have extremely large population sizes. Thus, the rate of genetic drift should be slower than that observed for other tunas.     

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.     

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     

Environmental preferences of bigeye tuna, <Emphasis Type="Italic">Thunnus obesus</Emphasis>, in the Indian Ocean: an application to a longline fishery

A survey of the fishing grounds for bigeye tuna, Thunnus obesus, in the Indian Ocean was carried out for a better understanding of the environmental preferences of bigeye tuna in a longline fishery. Catch rates of bigeye tuna were analyzed with respect to the ranges of depth, temperature, salinity, chlorophyll-a, and dissolved oxygen. The optimum capture depth, water temperature, and dissolved oxygen range of bigeye tuna were identified as 240.0 m to 279.9 m, 12.0°C to 13.9°C, and 2.00 mg·L⁻¹ to 2.99 mg·L⁻¹, respectively, in the study area of Indian Ocean. Neither salinity nor chlorophyll-a had a detectable effect on the vertical distribution of the adult bigeye tuna. The dissolved oxygen is the principal factor limiting the vertical distribution of bigeye tuna.     

Microsatellite DNA markers for population‐genetic studies of Atlantic bluefin tuna (Thunnus thynnus thynnus) and other species of genus Thunnus

Twenty‐five microsatellites from Atlantic bluefin tuna (Thunnus thynnus thynnus) were characterized. All 25 microsatellites were polymorphic; the number of alleles among up to 56 individuals surveyed ranged from two to 23. Atlantic bluefin tuna are highly exploited and major questions remain as to stock structure and abundance in the eastern and western North Atlantic. The microsatellites will be useful in testing stock‐structure hypotheses and in generating estimates of effective population size. The polymerase chain reaction primer sets developed also amplified identifiable alleles in three other species of genus Thunnus: T. albacares (yellowfin tuna), T. alalunga (albacore tuna) and T. obesus (bigeye tuna).     

Incorrect Likelihood Methods Were Used to Infer Scaling Laws of Marine Predator Search Behaviour

Background

Ecologists are collecting extensive data concerning movements of animals in marine ecosystems. Such data need to be analysed with valid statistical methods to yield meaningful conclusions.

Principal Findings

We demonstrate methodological issues in two recent studies that reached similar conclusions concerning movements of marine animals (Nature 451∶1098; Science 332∶1551). The first study analysed vertical movement data to conclude that diverse marine predators (Atlantic cod, basking sharks, bigeye tuna, leatherback turtles and Magellanic penguins) exhibited “Lévy-walk-like behaviour”, close to a hypothesised optimal foraging strategy. By reproducing the original results for the bigeye tuna data, we show that the likelihood of tested models was calculated from residuals of regression fits (an incorrect method), rather than from the likelihood equations of the actual probability distributions being tested. This resulted in erroneous Akaike Information Criteria, and the testing of models that do not correspond to valid probability distributions. We demonstrate how this led to overwhelming support for a model that has no biological justification and that is statistically spurious because its probability density function goes negative. Re-analysis of the bigeye tuna data, using standard likelihood methods, overturns the original result and conclusion for that data set. The second study observed Lévy walk movement patterns by mussels. We demonstrate several issues concerning the likelihood calculations (including the aforementioned residuals issue). Re-analysis of the data rejects the original Lévy walk conclusion.

Conclusions

We consequently question the claimed existence of scaling laws of the search behaviour of marine predators and mussels, since such conclusions were reached using incorrect methods. We discourage the suggested potential use of “Lévy-like walks” when modelling consequences of fishing and climate change, and caution that any resulting advice to managers of marine ecosystems would be problematic. For reproducibility and future work we provide R source code for all calculations.     

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua)

Background  

Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution.     

Population genetics of the bigeye thresher shark <Emphasis Type="Italic">Alopias superciliosus</Emphasis> in the Atlantic and Indian Oceans: implications for conservation

Population structure and genetic connectivity are pivotal contributions to the establishment of conservation strategies for fisheries management, in particular for highly migratory species that are affected by commercial fisheries. This study used partial sequences of mitochondrial DNA control region to determine the genetic structure of the bigeye thresher shark Alopias superciliosus in the Atlantic and Indian Oceans. A total of 858 base pairs of mtDNA CR from 228 individuals were analyzed. The resulting nucleotide diversity (π) was 0.0011?±?0.0008 and the haplotype diversity (h) was 0.127?±?0.030. These are the lowest diversities registered in elasmobranchs with this genetic marker. Two genetically distinct lineages were identified, one of them represented by 3.9% of the analyzed individuals and none restricted to any particular area. Simulated scenarios of population structure, tested with AMOVA and pairwise Φ_ST did not result in significant values indicating high connectivity among all sampled groups. The absence of population structure, even between Atlantic and Indian Oceans, corroborates the high dispersal ability of this species. The low genetic diversity detected in this species and the identification of two historical lineages occurring in sympatry, one represented by a very small number of individuals, should be considered in the conservation efforts and management plans of A. superciliosus.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Tuna longline fishing around West and Central Pacific seamounts

Background

Seamounts have been identified as aggregating locations for pelagic biodiversity including tuna; however the topography and prevailing oceanography differ between seamounts and not all are important for tuna. Although a relatively common feature in oceanic ecosystems, little information is available that identifies those that are biologically important. Improved knowledge offers opportunities for unique management of these areas, which may advance the sustainable management of oceanic resources. In this study, we evaluate the existence of an association between seamounts and tuna longline fisheries at the ocean basin scale, identify significant seamounts for tuna in the western and central Pacific Ocean, and quantify the seamount contribution to the tuna longline catch.

Methodology/Principal Findings

We use data collected for the Western and Central Pacific Ocean for bigeye, yellowfin, and albacore tuna at the ocean basin scale. GLMs were applied to a coupled dataset of longline fisheries catch and effort, and seamount location information. The analyses show that seamounts may be associated with an annual longline combined catch of 35 thousand tonnes, with higher catch apparent for yellowfin, bigeye, and albacore tuna on 17%, 14%, and 14% of seamounts respectively. In contrast 14%, 18%, and 20% of seamounts had significantly lower catches for yellowfin, bigeye and albacore tuna respectively. Studying catch data in relation to seamount positions presents several challenges such as bias in location of seamounts, or lack of spatial resolution of fisheries data. Whilst we recognize these limitations the criteria used for detecting significant seamounts were conservative and the error in identification is likely to be low albeit unknown.

Conclusions/Significance

Seamounts throughout the study area were found to either enhance or reduce tuna catch. This indicates that management of seamounts is important Pacific-wide, but management approaches must take account of local conditions. Management of tuna and biodiversity resources in the region would benefit from considering such effects.     

High‐resolution genetic analysis reveals extensive gene flow within the jellyfish Pelagia noctiluca (Scyphozoa) in the North Atlantic and Mediterranean Sea

Despite the importance of gelatinous zooplankton as components of marine ecosystems, both ecologically and socio‐economically, relatively little information is known about population persistence or connectivity in jellyfish. In the present study, we employed a combination of nuclear microsatellite markers and sequence data from the mitochondrial cytochrome oxidase I (COI) gene to determine levels and patterns of population genetic structuring in the holoplanktonic jellyfish Pelagia noctiluca across the northeast Atlantic Ocean and Mediterranean Sea. Our results indicate a high degree of connectivity in P. noctiluca, with little evidence of geographical structuring of genetic variation. A small but significant differentiation of Atlantic Ocean and Mediterranean stocks was detected based on the microsatellite data, but no evidence of differentiation was observed with the mtDNA, probably due to the higher power of the microsatellites to detect low levels of genetic structuring. Two clearly distinct groups of genotypes were observed within the mtDNA COI, which probably diverged in the early Pleistocene, but with no evidence of geographical structuring. Palaeodistribution modelling of P. noctiluca at the Last Glacial Maximum (LGM; c. 21 Kya) indicated large areas of suitable habitat south of the species’ current‐day distribution, with little reduction in area. The congruent evidence for minimal genetic differentiation from the nuclear microsatellites and the mtDNA, coupled with the results of the palaeodistribution modelling, supports the idea of long‐term population stability and connectivity, thus providing key insights into the population dynamics and demography of this important species.     

Discrimination of juvenile yellowfin (Thunnus albacares) and bigeye (T. obesus) Tunas using mitochondrial DNA control region and liver morphology

Yellowfin tuna, Thunnus albacares (Bonnaterre, 1788) and bigeye tuna, Thunnus obesus (Lowe, 1839) are two of the most economically important tuna species in the world. However, identification of their juveniles, especially at sizes less than 40 cm, is very difficult, often leading to misidentification and miscalculation of their catch estimates. Here, we applied the mitochondrial DNA control region D-loop, a recently validated genetic marker used for identifying tuna species (Genus Thunnus), to discriminate juvenile tunas caught by purse seine and ringnet sets around fish aggregating devices (FADs) off the Southern Iloilo Peninsula in Central Philippines. We checked individual identifications using the Neighbor-Joining Method and compared results with morphometric analyses and the liver phenotype. We tested 48 specimens ranging from 13 to 31 cm fork length. Morpho-meristic analyses suggested that 12 specimens (25%) were bigeye tuna and 36 specimens (75%) were yellowfin tuna. In contrast, the genetic and liver analyses both showed that 5 specimens (10%) were bigeye tuna and 43 (90%) yellowfin tuna. This suggests that misidentification can occur even with highly stringent morpho-meristic characters and that the mtDNA control region and liver phenotype are excellent markers to discriminate juveniles of yellowfin and bigeye tunas.     

Population genetic structure in Atlantic and Pacific Ocean common murres (Uria aalge): natural replicate tests of post‐Pleistocene evolution

Understanding the factors that influence population differentiation in temperate taxa can be difficult because the signatures of both historic and contemporary demographics are often reflected in population genetic patterns. Fortunately, analyses based on coalescent theory can help untangle the relative influence of these historic and contemporary factors. Common murres (Uria aalge) are vagile seabirds that breed in the boreal and low arctic waters of the Northern Hemisphere. Previous analyses revealed that Atlantic and Pacific populations are genetically distinct; however, less is known about population genetic structure within ocean basins. We employed the mitochondrial control region, four microsatellite loci and four intron loci to investigate population genetic structure throughout the range of common murres. As in previous studies, we found that Atlantic and Pacific populations diverged during the Pleistocene and do not currently exchange migrants. Therefore, Atlantic and Pacific murre populations can be used as natural replicates to test mechanisms of population differentiation. While we found little population genetic structure within the Pacific, we detected significant east–west structuring among Atlantic colonies. The degree that population genetic structure reflected contemporary population demographics also differed between ocean basins. Specifically, while the low levels of population differentiation in the Pacific are at least partially due to high levels of contemporary gene flow, the east–west structuring of populations within the Atlantic appears to be the result of historic fragmentation of populations rather than restricted contemporary gene flow. The contrasting results in the Atlantic and Pacific Oceans highlight the necessity of carefully considering multilocus nonequilibrium population genetic approaches when reconstructing the demographic history of temperate Northern Hemisphere taxa.     

Comparative genetic stock structure in three species of commercially exploited Indo-Malay Carangidae (Teleosteii,Perciformes)

We examine genetic structuring in three commercially important species of the teleost family Carangidae from Malaysian waters: yellowtail scad Atule mate, bigeye scad Selar crumenophthalmus and yellowstripe scad Selaroides leptolepis, from the Indo-Malay Archipelago. In view of their distribution across contrasting habitats, we tested the hypothesis that pelagic species display less genetic divergence compared with demersal species, due to their potential to undertake long-distance migrations in oceanic waters. To evaluate population genetic structure, we sequenced two mitochondrial (mt)DNA [650 bp of cytochrome oxidase I (coI), 450 bp of control region (CR)] and one nuclear gene (910 bp of rag1) in each species. One hundred and eighty samples from four geographical regions within the Indo-Malay Archipelago including a population of yellowtail from Kuwait were examined. Findings revealed that the extent of genetic structuring among populations in the semi-pelagic and pelagic, yellowtail and bigeye were lower than demersal yellowstripe, consistent with the hypothesis that pelagic species display less genetic divergence compared with demersal species. The yellowtail phylogeny identified three distinct clades with bootstrap values of 86%–99% in mtDNA and 63%–67% in rag1. However, in bigeye, three clades were also observed from mtDNA data while only one clade was identified in rag1 dataset. In yellowstripe, the mtDNA tree was split into three closely related clades and two clades in rag1 tree with bootstraps value of 73%–99% and 56% respectively. However, no geographic structure appears in both mtDNA and rag1 datasets. Hierarchical molecular variance analysis (AMOVA), pair wise F_ST comparisons and the nearest-neighbour statistic (S_nn) showed significant genetic differences among Kuwait and Indo-Malay yellowtail. Within the Indo-Malay Archipelago itself, two distinct mitochondrial lineages were detected in yellowtail suggesting potential cryptic species. Findings suggests varying degrees of genetic structuring, key information relevant to management of exploited stocks, though more rapidly evolving genetic markers should be used in future to better delimit the nature and dynamics of putative stock boundaries.     

Global isolation by distance despite strong regional phylogeography in a small metazoan

Background  

Small vagile eukaryotic organisms, which comprise a large proportion of the Earth's biodiversity, have traditionally been thought to lack the extent of population structuring and geographic speciation observed in larger taxa. Here we investigate the patterns of genetic diversity, amongst populations of the salt lake microscopic metazoan Brachionus plicatilis s. s. (sensu stricto) (Rotifera: Monogononta) on a global scale. We examine the phylogenetic relationships of geographic isolates from four continents using a 603 bp fragment of the mitochondrial COI gene to investigate patterns of phylogeographic subdivision in this species. In addition we investigate the relationship between genetic and geographic distances on a global scale to try and reconcile the paradox between the high vagility of this species and the previously reported patterns of restricted gene flow, even over local spatial scales.     

Evidence for panmixia despite barriers to gene flow in the southern African endemic, <Emphasis Type="Italic">Caffrogobius caffer</Emphasis> (Teleostei: Gobiidae)

Background  

Oceanography and life-history characteristics are known to influence the genetic structure of marine species, however the relative role that these factors play in shaping phylogeographic patterns remains unresolved. The population genetic structure of the endemic, rocky shore dwelling Caffrogobius caffer was investigated across a known major oceanographic barrier, Cape Agulhas, which has previously been shown to strongly influence genetic structuring of South African rocky shore and intertidal marine organisms. Given the variable and dynamic oceanographical features of the region, we further sought to test how the pattern of gene flow between C. caffer populations is affected by the dominant Agulhas and Benguela current systems of the southern oceans.     

A Validated Methodology for Genetic Identification of Tuna Species (Genus Thunnus)

Background

Tuna species of the genus Thunnus, such as the bluefin tunas, are some of the most important and yet most endangered trade fish in the world. Identification of these species in traded forms, however, may be difficult depending on the presentation of the products, which may hamper conservation efforts on trade control. In this paper, we validated a genetic methodology that can fully distinguish between the eight Thunnus species from any kind of processed tissue.

Methodology

After testing several genetic markers, a complete discrimination of the eight tuna species was achieved using Forensically Informative Nucleotide Sequencing based primarily on the sequence variability of the hypervariable genetic marker mitochondrial DNA control region (mtDNA CR), followed, in some specific cases, by a second validation by a nuclear marker rDNA first internal transcribed spacer (ITS1). This methodology was able to distinguish all tuna species, including those belonging to the subgenus Neothunnus that are very closely related, and in consequence can not be differentiated with other genetic markers of lower variability. This methodology also took into consideration the presence of introgression that has been reported in past studies between T. thynnus, T. orientalis and T. alalunga. Finally, we applied the methodology to cross-check the species identity of 26 processed tuna samples.

Conclusions

Using the combination of two genetic markers, one mitochondrial and another nuclear, allows a full discrimination between all eight tuna species. Unexpectedly, the genetic marker traditionally used for DNA barcoding, cytochrome oxidase 1, could not differentiate all species, thus its use as a genetic marker for tuna species identification is questioned.     

Large‐scale distribution of tuna species in a warming ocean

Tuna are globally distributed species of major commercial importance and some tuna species are a major source of protein in many countries. Tuna are characterized by dynamic distribution patterns that respond to climate variability and long‐term change. Here, we investigated the effect of environmental conditions on the worldwide distribution and relative abundance of six tuna species between 1958 and 2004 and estimated the expected end‐of‐the‐century changes based on a high‐greenhouse gas concentration scenario (RCP8.5). We created species distribution models using a long‐term Japanese longline fishery dataset and two‐step generalized additive models. Over the historical period, suitable habitats shifted poleward for 20 out of 22 tuna stocks, based on their gravity centre (GC) and/or one of their distribution limits. On average, tuna habitat distribution limits have shifted poleward 6.5 km per decade in the northern hemisphere and 5.5 km per decade in the southern hemisphere. Larger tuna distribution shifts and changes in abundance are expected in the future, especially by the end‐of‐the‐century (2080–2099). Temperate tunas (albacore, Atlantic bluefin, and southern bluefin) and the tropical bigeye tuna are expected to decline in the tropics and shift poleward. In contrast, skipjack and yellowfin tunas are projected to become more abundant in tropical areas as well as in most coastal countries' exclusive economic zones (EEZ). These results provide global information on the potential effects of climate change in tuna populations and can assist countries seeking to minimize these effects via adaptive management.     

Genetic Structure of Bluefin Tuna in the Mediterranean Sea Correlates with Environmental Variables

Background

Atlantic Bluefin Tuna (ABFT) shows complex demography and ecological variation in the Mediterranean Sea. Genetic surveys have detected significant, although weak, signals of population structuring; catch series analyses and tagging programs identified complex ABFT spatial dynamics and migration patterns. Here, we tested the hypothesis that the genetic structure of the ABFT in the Mediterranean is correlated with mean surface temperature and salinity.

Methodology

We used six samples collected from Western and Central Mediterranean integrated with a new sample collected from the recently identified easternmost reproductive area of Levantine Sea. To assess population structure in the Mediterranean we used a multidisciplinary framework combining classical population genetics, spatial and Bayesian clustering methods and a multivariate approach based on factor analysis.

Conclusions

F_ST analysis and Bayesian clustering methods detected several subpopulations in the Mediterranean, a result also supported by multivariate analyses. In addition, we identified significant correlations of genetic diversity with mean salinity and surface temperature values revealing that ABFT is genetically structured along two environmental gradients. These results suggest that a preference for some spawning habitat conditions could contribute to shape ABFT genetic structuring in the Mediterranean. However, further studies should be performed to assess to what extent ABFT spawning behaviour in the Mediterranean Sea can be affected by environmental variation.