Multilocus Bayesian Estimates of Intra-Oceanic Genetic Differentiation, Connectivity, and Admixture in Atlantic Swordfish (Xiphias gladius L.)

Previous genetic studies of Atlantic swordfish (Xiphias gladius L.) revealed significant differentiation among Mediterranean, North Atlantic and South Atlantic populations using both mitochondrial and nuclear DNA data. However, limitations in geographic sampling coverage, and the use of single loci, precluded an accurate placement of boundaries and of estimates of admixture. In this study, we present multilocus analyses of 26 single nucleotide polymorphisms (SNPs) within 10 nuclear genes to estimate population differentiation and admixture based on the characterization of 774 individuals representing North Atlantic, South Atlantic, and Mediterranean swordfish populations. Pairwise F ST values, AMOVA, PCoA, and Bayesian individual assignments support the differentiation of swordfish inhabiting these three basins, but not the current placement of the boundaries that separate them. Specifically, the range of the South Atlantic population extends beyond 5°N management boundary to 20°N-25°N from 45°W. Likewise the Mediterranean population extends beyond the current management boundary at the Strait of Gibraltar to approximately 10°W. Further, admixture zones, characterized by asymmetric contributions of adjacent populations within samples, are confined to the Northeast Atlantic. While South Atlantic and Mediterranean migrants were identified within these Northeast Atlantic admixture zones no North Atlantic migrants were identified respectively in these two neighboring basins. Owing to both, the characterization of larger number of loci and a more ample spatial sampling coverage, it was possible to provide a finer resolution of the boundaries separating Atlantic swordfish populations than previous studies. Finally, the patterns of population structure and admixture are discussed in the light of the reproductive biology, the known patterns of dispersal, and oceanographic features that may act as barriers to gene flow to Atlantic swordfish.


Introduction
The epipelagic realm of the world's oceans is a relatively homogeneous and contiguous environment characterized by the near absence of physical barriers. These features when coupled with the high dispersal behavior and demographic characteristics of many marine species result in high levels of gene flow [1]. Specifically, many epipelagic species of fishes are characterized by long distance potential dispersal of pelagic eggs and larvae facilitated by currents, high vagility of juveniles and adults, opportunistic feeding behaviors, eurythermal physiologies of adults, and large effective population sizes. Consequently, the levels of genetic differentiation in marine fishes are typically substantially lower compared to populations of anadromous and freshwater fishes [2]. It is thus not surprising that the distributions of many species of tunas, billfishes, dolphinfishes, lamnid sharks, are cosmopolitan [3,4], with some species, such as skipjack tuna (Katsuwonus pelamis) displaying no genetic differentiation among-oceans [5], and others, like wahoo (Acanthocybium solandri) showing only incipient levels inter-oceanic differentiation [6]. Significant inter-oceanic differentiation, however, has been documented in many highly migratory fishes, including bigeye tuna [7], Atlantic bluefin tuna [8,9], albacore tuna [10,11], and many species of billfishes [12] and swordfish [13]. Similarly, genetic investigations in several species of billfish have shown intra-oceanic genetic differentiation within the Pacific Ocean, including sailfish [14], white marlin, and blue marlin [12]. Within the Atlantic, no significant population substructure was detected in several epipelagic fishes; including shortfin mako (Isurus oxyrinchus) [15], blue marlin (Makaira nigricans) and sailfish (Istiophorus platypterus) [12], bigeye tuna (Thunnus obesus) [16], and in both yellowfin tuna (Thunnus albacares) and skipjack tuna [5]. By contrast, swordfish (Xiphias gladius L.) shows significant genetic differentiation within the Atlantic Ocean [13,17,18] and also within the Mediterranean Sea [19].
Swordfish is a large epipelagic monotypic cosmopolitan species capable of long distance movements, as documented using both conventional tags [20] and pop-up satellite tags (PSATs) [21][22][23][24]. Swordfish are subject to intensive commercial exploitation worldwide, and in the Atlantic Ocean they are managed by the International Commission for the Conservation of Tunas (ICCAT) as two separate stocks: the North Atlantic and the South Atlantic separated at 5°N. In turn, North Atlantic swordfish is managed separately from Mediterranean swordfish with a boundary at the Strait of Gibraltar. Previous genetic studies on Atlantic swordfish employing mitochondrial DNA (mtDNA) (see [25] for summary), single copy nuclear DNA (scnDNA) [18,26,27], and to lesser exent microsatellite loci [28][29][30], confirm the differentatiation among North Atlantic, South Atlantic and Mediterranean populations. Some of these studies also suggested that the range of these swordfish populations extend beyond current management boundaries [26,31]. However, limitations in geographic sampling coverage, and the use of single loci, precluded an accurate placement of boundaries and of estimates of admixture. Multilocus analyses of microsatellite data [29,30] attempted to refine these estimates by expanding sampling coverage and by utilizing Bayesian individual assignment to identify zones of admixture. While, microsatellite data confirmed the pronounced differences between Atlantic and Mediterranean swordfish [29], it failed to delineate the separation of North Atlantic and South Atlantic swordfish populations.
In recent years there has been an increase in the adoption of single nucleotide polymorphisms (SNPs) as an alternative to microsatellite data in Bayesian analyses of population admixture. In addition to a comparable differentiation power, SNPs have lower associated genotyping errors, simpler mutation models, higher amenability to automation and high throughput technologies, and increased amplification success with low-quality and ancient DNA as compared to microsatellites [32][33][34]. Further, the levels of genetic differentiation among Atlantic swordfish populations using SNPs was an entire order of magnitude greater compared to microsatellites [35,36]. Accordingly, in this study we evaluate the patterns of population differentiation and admixture of Atlantic and Mediterranean swordfish based on Bayesian analyses of multilocus nuclear SNPs using a more extensive and representative geographic sampling coverage compared to previous studies. We describe the patterns of differentiation and admixture among North Atlantic, South Atlantic and Mediterranean swordfish populations.

Materials and Methods Sampling
A total of 774 swordfish specimens from 18 localities in the Atlantic and western Mediterranean were collected from 1991-2006 (Table 1 and Fig 1). Muscle, liver, and heart tissue samples of were collected by observers on commercial longline vessels or obtained from collaborating researchers. Ethical approval was not required for the collection of adult swordfish, as these were collected during routine commercial operations. Swordfish adults were sacrificed by the fishermen by decapitation. In all instances the location of capture for each fish was recorded. Larval swordfish (n = 52), collected as part of icthyoplankton surveys in the northern Gulf of Mexico during the summers of 2005-2006 [37], were also characterized. Permits for fish larvae collections in the Gulf of Mexico were issued by the Highly Migratory Species Management Division of the National Oceanic and Atmospheric Administration.
recognition particle 54 (SRP54). Details on loci and corresponding SNPs selection, the minimization of ascertainment bias as well as the criteria to design primers and probes that target these, are given in Smith et al. [35]. Briefly, to minimize ascertainment bias potential loci and SNPs were selected after conducting a preliminary characterization with HRMA of 30 swordfish, 10 from each of three geographically disjunct swordfish populations (NW Atlantic, Mediterranean and Pacific). Those amplicons that generated distinct melting profiles, were sequenced, and electropherogram inspections were then used to validate the SNPs responsible for the polymorphic melting curves (see [35] for details). Sequence alignments were also used to evaluate the placement and design of new HRMA-genotyping primers and probes, when required. A total of 85 SNPs were initially identified but only those differing at frequencies >5% among reference samples were targeted for further HRMA primer and or probe design. PCR conditions of selected potential loci were further optimized and validated using additional samples (n = 40) from each of the three reference populations. This process ultimately yielded a total of 26 SNPs among 10 loci that we targeted for genotyping (see Table 2 in [35]). Procedures for the amplification and subsequent genotyping of SNPs for short amplicon high-resolution melting analysis (SA-HRMA) of ARP, ATPsβ, CaM, and SRP54 and unlabeled probe high-resolution melting analysis (UP-HRMA) of ANT, GpHR, ldhA, and Mlc2 are described in Smith et al. [35]. Two additional loci, AldB and Act2α, were scored as size polymorphisms as follows. A simple sequence repeat (SSR) identified in AldB [39] was amplified with the following primers; 5'-VIC-TGTGCCCAGTATAAGAAGGATGG-3' and 5'-CTGTGGAGAATCAGGGCTCC-3' (JX042447). Polymerase chain reactions (PCR) were conducted in 12.5 μL reactions containing 10 ng of genomic DNA, 1X Econotaq Green Master Mix (Lucigen), and 0.20 μM of each primer. Thermocycling was performed on an Eppendorf Mastercycler (Eppendorf) with an initial denaturation of 10 min at 95°C followed by 35 cycles denaturing for 1 min at 94°C, annealing for 1 min at 54°C, and extension for 1 min at 72°C. PCR products were diluted 1:10 and 1 μL of diluted PCR template and 0.2 μL of ROX size standard (Life Technologies) in 10 μL of Hi-Di formamide were run on an ABI 3130 Genetic Analyzer (Applied Biosystems). SSRs were scored by size and analyzed using GeneMapper v.4.0 (Applied Biosystems).
Locus Act2α contains an SNP that corresponds to a restriction site of endonuclease Hpy8I and was amplified using primers; 5'-GTCACCGGAGTCCAGGACG-3' and 5'-ATCTGGCAC CACACCTTCTACAA-3' (JX042448) using the same chemistry and thermocycling profile used in AldB amplification. Prior to digestion, the quality and quantity of PCR products were visualized in 1% Tris-acetate (TA) agarose gels via gel electrophoresis prior to digests. Restriction digests were conducted in 10 μL volumes containing 2-4 μL PCR product, 1X Buffer TANGO (Fermentas), and 0.4 U of Hpy8I (Fermentas), that were incubated for 16 hours at 37°C, followed by thermal inactivation at 80°C for 5 min. Restriction Fragment Length Polymorphisms (RFLPs) were separated in 2% TA agarose gels run at 100V for 30 min, and were visualized and captured on a Gel Doc XR (Bio-Rad) UV transilluminator using the Quality One v.4.6 (Bio-Rad) software.

Data analysis
Genetic diversity. The number of alleles, observed (H o ) and expected (H e ) heterozygosities under Hardy-Weinberg equilibrium (HWE), and inbreeding coefficients (F IS ) were calculated using GENALEX version 6.41 [40]. Departures from HWE were estimated using an exact probability test [41,42] and global linkage disequilibrium using a maximum likelihood-ratio test [43] in GENEPOP version 4.1 (Markov chain parameters: 10000 dememorization steps, 1000 batches, 10000 iterations per batch) [44,45]. Sequential Bonferroni corrections were applied to adjust the levels of statistical significance for multiple comparisons [46]. An F ST -outlier detection method, in which observed locus F ST values are compared to calculated global F ST values expected under neutrality using coalescent simulations [47] was performed in LOSI-TAN [48]. A total of 50,000 simulations were conducted to detect putative loci under selection. To remove potential selective loci when computing the initial mean F ST , the first simulation used the 'neutral' mean F ST option with a 0.995 confidence interval, a false discovery rate set at 0.1, and by pooling the sampling localities into the corresponding three populations (K = 3) identified by Bayesian clustering analyses (see below) as follows, North Atlantic (1-7, 10-11), South Atlantic (12)(13)(14)(15)(16)(17)(18), and Mediterranean (8)(9). The assignment power of each locus was ranked a posteriori to Bayesian analysis using the critical population method in WHICHLOCI v1.0 [49] using samples from localities identified as non-mixing areas.
Population differentiation. A hierarchical analysis of molecular variance (AMOVA) [50] was implemented in Arlequin version 3.5 [51] to estimate levels of population subdivision. Population differentiation was also evaluated using both global and pairwise F ST tests [52]. Initially, AMOVA groupings of two (MED and Atlantic) and three (MED, NA, and SA) adhering to the current management regions were evaluated. Alternative groupings in AMOVA were evaluated a posteriori to F ST and Bayesian analyses (see below). Slatkin's linearized F ST [53] was also calculated in Arlequin. Principal Coordinates Analysis (PCoA) [54], based on a standardized covariance matrix of genetic distances between each locality pair, was carried in Gen-AlEx version 6.4 [40].
Genetic clustering analysis. Genetic population structure and patterns of inter-population gene flow were assessed by Bayesian inference implemented in STRUCTURE version 2.3 [55]. Previous genetic studies using mtDNA and a single nuclear locus identified comparatively low levels of gene flow in between North and South Atlantic swordfish [17,26,56]. Accordingly, a no admixture ancestry with correlated allele frequencies model, as outlined in Falush et al. [57] was adopted. Compared to freshwater and anadromous fishes, the populations of marine fishes display weak levels of differentiation (F ST < 0.20) [58]. Accordingly, the LOCPRIOR option was implemented to detect weak signals of population structure using an a priori grouping as outlined in Hubisz et al. [59]. The number of clusters (K) was estimated using an ad hoc approach [55] by obtaining the mean posterior probability of the data (L(K)) and the ΔK approach of Evanno et al. [60] using STRUCTURE HARVESTER v0.6.92 [61]. Twenty independent runs for each K value (1-10) were performed using 100,000 Markov chain Monte Carlo (MCMC) iterations with a burn-in period of 100,000. Results from STRUCTURE were compared with the results from GENELAND [62], which incorporates individual spatial and genetic data to infer population structure and spatial boundaries between clusters. Individual spatial coordinates were set to the corresponding latitude and longitude of sampling. An uncertainty value on the spatial coordinates of each sample was set to 30 decimal degrees corresponding to the documented displacement of North Atlantic swordfish from tagging experiments [22,23]. Correlated allele and spatial models were implemented for twenty independent runs with 100,000 MCMC iterations and a burn-in of 10,000 iterations for each K value (1)(2)(3)(4)(5)(6)(7)(8). After estimating K, 25 independent runs with the optimal K value were run in GENELAND, 50 independent runs in STRUCTURE, and the optimal alignment of replicates were completed using CLUMPP version 1.1.12 [63].

Results
Hardy-Weinberg, linkage disequilibrium, and F ST outlier analysis A total of 774 swordfish were genotyped successfully at the ten loci (see supporting S5 Table). Five of these loci contained a single SNP and were bi-allelic and the rest contained both multiple SNPs and alleles. The number of alleles per locus ranged from 2 to 6, and in total 26 SNPs were characterized (see supporting S1 Table). Average (H o ) and expected (H e ) heterozygosities, and inbreeding coefficients (F IS ) are given in Table 2. None of the Fisher's exact probability tests for departure from HWE were significant (P < 0.05) after Bonferroni correction. Genotypic linkage was not significant among the 44 possible pairs of nuclear loci across all samples with (P < 0.01) and without (P < 0.05) Bonferroni correction. Likewise exact tests for genotypic linkage disequilibrium within each sample were also non-significant after Bonferroni correction (P < 0.05) (not shown). Initial F ST outlier analysis of the ten loci of all individuals from the 18 localities identified CaM as a potentially under selection (P < 0.05) (Fig 2A). However, since both STRUCTURE and GENELAND identified several localities (6)(7)(10)(11) in the North Atlantic as admixture zones there was the possibility that heterozygote deficits biased the outlier analysis. Consequently, this analysis was conducted excluding admixed samples resulting in no locus identified as an outlier (P < 0.05). This procedure also yielded higher F ST values for all loci as well as higher H e values for Act2α, AldB, ATPsβ, CaM, and Mlc2. (Fig 2B). The power assignment of each loci using WHICHLOCI v1.0 identified CaM as the locus with the highest power of discrimination for North and South Atlantic populations followed by Mlc2, SRP54, AldB, ARP, Act2α, ldhA, GpHR, ATPsβ, and ANT respectively (see supporting S2 Table). Interestingly, CaM displayed the lowest power of discrimination between North Atlantic and Mediterranean populations, with Mlc2 yielding the highest power of discrimination between them followed by ATPsβ, Act2α, ldhA, ANT, AldB, SRP54, ARP, and GpHR, respectively (see supporting S3 Table).

Population differentiation
Pairwise F ST . Multilocus pairwise F ST values among the 18 populations ranged considerably (<0.001 to 0.121) ( Table 3). The comparison between Mediterranean and South Atlantic populations yielded the highest indices of differentiation (F ST >0.081 -0.121). Several sampling localities identified as admixture zones (6-7, 10-11) were excluded when estimating gene flow among North Atlantic (1-5), Mediterranean (8-9) and South Atlantic (12-18) populations. The lowest level of gene flow of about 5 migrants per generation (Nm) was estimated between Mediterranean and South Atlantic populations, compared to about 9 migrants per generation between North Atlantic and South Atlantic swordfish (Table 4). Slatkin's linearized F ST values (Fig 3) highlights the homogeneity between the Mediterranean sample collected to the east of the Alboran Sea (8) and the adjacent Northeast Atlantic area west of the Strait of Gibraltar (9). The magnitude of these fixation indices also underscore the high degree of differentiation of Mediterranean swordfish with respect to all the Atlantic localities (1-7, 10-18) including those inhabiting the adjacent NE Atlantic waters off both the Iberian Peninsula and the NW African coast. In addition, Fig 2 illustrates that the North Atlantic localities (1-5) are different from the South Atlantic localities (13-18) but also from locality 12 which corresponds geographically to the North Atlantic stock since it is located north of the 5°N management boundary. Linearized F ST values depict a geographic gradient of differentiation from the Northwest Atlantic to the Southeast Atlantic indicative of mixing between the South Atlantic and North Atlantic populations, primarily in the eastern localities of Iberian, Morocco, and Western Sahara (7,(10)(11) and to a lesser extent in the central North Atlantic (6) locality.
AMOVA. In all alternative hierarchical arrangements tested (two-and three-group; with all localities and excluding admixture areas) with AMOVA ( Table 5) the majority of the variation (>90%) in swordfish was contained within individuals. In all arrangements, amonggroups (F CT ) differentiation was highly significant (P<0.01). However, the proportions of variation explained by among-groups were slightly larger in two-group AMOVAs, with 6.63% with all localities and 7.05%, excluding admixture areas (localities 6, 10-11), than in threegroup arrangements that yielded 5.17% and 6.45% with and without admixture areas, respectively. In addition, the proportion of variance explained by differences among populations within groups (F SC ) was four to ten times larger in the two-groups arrangements than in the three-groups arrangements, with the lowest proportion of variation (0.19%) for the threegroups arrangement without admixture areas (F SC = 0.00207, P>0.05; North Atlantic: South Atlantic: Mediterranean). Principal coordinate analysis. The PCoA of the standardized covariance of pairwise population genetic distances provide concordant evidence of locality assignment to a minimum of three separate swordfish populations (Fig 4). The first three axes accounted for~93.86% of the variation among populations with Eigen values of 0.847 and 0.328 for the first and second axes, respectively. Based on the clustering patterns, three groups can be identified corresponding to North Atlantic, South Atlantic, and Mediterranean populations. Loadings in the first axis contrast the Mediterranean from the South Atlantic, whereas the second axis contrasts the Mediterranean with the North Atlantic. The North Atlantic cluster, is not as tightly grouped as other populations, and is comprised of the Gulf of Mexico (1), Southeast U.S. Coast (2),  (9), and together they form the Mediterranean group. Finally, the South Atlantic group included the samples from Cape Verde (12), collected north of the current boundary between northern and southern Atlantic swordfish stocks, which cluster tightly with the samples from Equatorial Brazil (13), Gulf of Guinea (14), Central South Atlantic (15), Brazil (16), Brazil-Uruguay (17), and Namibia  Table 3. Pairwise  Table 1). 1  (6) and Iberian (7) samples associating more closely with the North Atlantic group, whereas Morocco (10) and Western Sahara (11) samples more closely with the South Atlantic group.
Bayesian genetic clustering analyses STRUCTURE analysis. Evaluation of the mean posterior probabilities from multiple STRUCTURE analyses (K = 1-10) revealed that the mean ln P(D) increased from K = 1-3 and then sharply decreased when K 4 (see supporting S1 Fig) indicating that K = 3 reflects the majority of the genetic structure in the dataset [55]. Estimated ΔK's revealed two high ΔK peaks of similar magnitude corresponding to K = 2 (ΔK 2 = 35.887) and K = 3 (ΔK 3 = 33.6548), followed by very low ΔK (< 5.0) values (see supporting S1 Fig). Mediterranean localities (8)(9) were removed from the analysis to test whether the strong differentiation of the Mediterranean from the Atlantic was causing an underestimation of K. The evaluation including only North Atlantic and South Atlantic localities identified a single ΔK peak at K = 2 (not shown) congruent with the genetic heterogeneity between the North and South Atlantic populations of swordfish. While Evanno's ΔK method seeks to detect the uppermost hierarchical level of population structure the method is less reliable at lower levels of genetic differentiation and may underestimate K [64]. Therefore, the estimation of K using the ad hoc evaluation of mean posterior probabilities from multiple analyses of K [55] appears to be more appropriate in this instance than the ΔK approach [60].
Average population posterior probability memberships ( Q) of the 18 localities (see supporting S4 Table) and individual assignments at K = 3 were concordant with results from PCoA, AMOVA, and F ST analyses (Fig 5). Localities 1-5 had average values of ancestry probabilities ( Q ) belonging to the North Atlantic > 0.90 with no individuals of South Atlantic or Mediterranean origin. Swordfish collected in the central North Atlantic (6) clustered between the North Atlantic and the South Atlantic indicating a potential mixing area, however sample size was small (n = 20). The Iberian (7) (8) were identified by STRUCTURE as Mediterranean. Of the remaining individuals, three swordfish (3.6%) were identified as North Atlantic and another three (3.6%) as South Atlantic. Four Mediterranean swordfish were also identified in the Iberian sample (7) as follows: three at 37.5°N 12.5°W and one at 37.5°N 22.5°W, accounting for roughly 11% and 5% of the total number of individuals sampled respectively in these two locations.
The area off the Atlantic coast of Morocco (10) was also identified as an area of mixing with individuals of North Atlantic and South Atlantic ( Q= 0.52) origin. Although the sample of Western Sahara (11) contained a majority of individuals assigned to the South Atlantic ( Q= 0.57), it is a zone of admixture containing also individuals of North Atlantic origin. Swordfish GENELAND analysis. The K value with the greatest density on the MCMC chain after a burn-in of 10,000 was K = 4. However, inspection of map of posterior probabilities identified one cluster as a 'ghost' with no associated samples (not shown). Further, this ghost cluster was unstable, as it's placing or presence differed among replicate runs. The existence of 'ghost'  Genetic Differentiation in Atlantic Swordfish clusters has been reported as an artifact of GENELAND which tends to overestimate K in a comparison to other Bayesian genetic clustering algorithms [65]. Accordingly, GENELAND was run with the optimal value (K = 3) identified by STRUCTURE. Localities 1-5 and part of 7 comprised the North Atlantic ( Q> 0.90), localities 8 and 9 the Mediterranean ( Q> 0.90), and localities 12-18 the South Atlantic ( Q> 0.90). The individual posterior probabilities assigned in GENELAND (Fig 6) identified the region corresponding to the central North Atlantic (6) and Western Sahara (11) samples as areas of mixing between the North Atlantic and South Atlantic populations. The samples corresponding to the North Atlantic area west of the Strait of Gibraltar (8) have the genetic signature of the western Mediterranean (9), whereas swordfish from the Atlantic locality of Morocco (10) were identified as either Mediterranean or North Atlantic. The Iberian sample (7) is comprised of three subsamples located at 37.5°N 22.5°W (n = 19), 37.5°N 17.5°W (n = 5), and 37.55°N 12.5°W (n = 27), with a North Atlantic origin assigned by GENELAND of > 80%, >70%, and <50%, respectively. Since the precise geographic coordinates of capture for the individuals in these Atlantic subsamples were recorded only with a precision of 2.5°latitude by 2.5°longitude, and because the subsample sizes are small, these estimates are imprecise.

Nuclear loci evaluation
Genetic population structure studies using nuclear DNA have traditionally excluded adaptive markers to avoid introducing bias to estimates of migration rates and effective population sizes  [66]. In this study we attempted to identify variation in non-coding nuclear introns to maximize the potential of characterizing neutral loci. Of the 26 SNPs assayed 25 were contained in the introns of 10 nuclear genes, with the exception of a SNP corresponding to a third codon synonymous substitution in one exon of the ATPsβ gene. Although introns are generally  assumed to be neutral, some may contain transcription regulation or splicing control elements that alter gene expression [67]. Introns may also experience purifying and/or positive selection [68] or the effects of hitch-hiking (i.e. linkage between a non-coding variant and a gene under selection) as reported in Atlantic cod [69]. Therefore deviations from HWE, linkage disequilibrium tests, and outlier analysis are commonly applied to detect potentially adaptive loci. In our study exact HWE [42] and linkage equilibrium [43] tests were not significant after Bonferroni corrections. Similarly the the F ST outlier analysis failed to identify any loci characterized as subject to selection. Given that the number of loci analyzed in the current study is small (n = 10), the associated error in the calculation of global neutral F ST could be large. Accordingly, future studies utilizing genomic technologies will certainly yield greater numbers of loci to conduct more conclusive outlier tests. Alternatively, outlier and non-neutral SNPs can provide insight into adaptive genetic diversity within and amoung populations [70,71] which may be useful for resolving fine-scale population connectivity, admixture, mixed stock analysis, and identifying local adaptation for conservation or resource management [72][73][74]. Studies utilizing adaptive loci in Atlantic cod [75], Atlantic herring [74,76], Atlantic salmon [77], and Pacific lamprey [78] reported that fewer loci and smaller sample sizes were required to detect population differentiation in marine fishes when using outlier loci as compared to neutral loci. Muths, Le Couls (79) hypothesized that the inclusion of adaptive loci in future genetic analses of large pelagic fishes, specifically swordfish, may identify additional evolutionary significant population units than analyses relying on only neutral loci. However, since the inclusion of outlier loci may violate model assumptions used for individual assignment, it is important that outlier loci included in such analyses are carefully evaluated and presented with comparative results using neutral loci [74]. While the ranking analysis identified CaM as the most informative locus for distinguishing North and South Atlantic swordfish populations, this locus alone fails to differentiate between North Atlantic and Mediterranean swordfish populations. Further, STRUCTURE results without CaM (see supporting S2 Fig) show that the remaining loci were capable to discriminate between North and South Atlantic swordfish. Accordingly, the high resolution of population differentiation reported in this study among North Atlantic, South Atlantic and Mediterranean swordfish cannot be attributed to CaM alone. In fact, we demonstrate that the inclusion of the additional loci increases the power of population discrimination.

Atlantic swordfish population differentiation
Analyses of multilocus nuclear SNP data indicate that Atlantic swordfish can be subdivided into at least three populations, namely Mediterranean, North Atlantic, and South Atlantic. The presence of these populations is supported by linearized pairwise F ST [52,53], AMOVA [50], PCoA [54], and Bayesian genetic clustering using both in STRUCTURE [55,57,59] and GENELAND [62]. Though no geographic barriers prevent gene flow among North Atlantic, South Atlantic and Mediterranean swordfish populations, the number of migrants among populations is small (Table 4). Additionally, the observed temporal stability in allele frequencies within the Northwest Atlantic, South Atlantic, and Mediterranean populations also suggests limited gene flow among these swordfish populations. This temporal stability is perhaps best evident in the parity of genetic signature between the larval Gulf of Mexico samples and the samples of adult Northwest Atlantic swordfish collected over a decade prior. Based on these results the North Atlantic population extends from the tropical spawning areas in northwestern Atlantic waters [80,81] to 50°W. In the feeding areas north of 40°N the range of the North Atlantic population extends to the Azores and may extend at that latitude as far east as 15°W. The South Atlantic population extends from 50°S through the South Atlantic tropical and equatorial spawning areas [17] to 20°N at 40°W and 25°N along the African coast.
The Mediterranean population extends past the Strait of Gibraltar, possibly as far as 8°W, with little evidence of Atlantic swordfish entering the Mediterranean Sea. These results are concordant with the phylogeographical association of mtDNA lineages previously reported [7,13,17,56] and also with evidence derived from fisheries data [82]. The Mediterranean population is identified here as the most divergent of the three populations, and these findings agree with other studies using mtDNA, scnDNA (see [25] for summary), and microsatellite loci [29,30].
The Bayesian clustering analyses based on multilocus SNPs presented here are not consistent with analyses of global swordfish population structure using microsatellites [28,30]. Those studies reported extremely low levels of differentiation between North and South Atlantic populations (F ST 0.0015), precluding individual cluster assignment. By contrast, Kasapidis et al. [29] argued that the considerably higher estimates of differentiation reported with nuclear SNPs, and specifically with CaM, is strong adaptive selection. Such explanation, however, fails to explain the concordance between mtDNA and nuclear SNPs data (reviewed in [25]), including the current study. Instead, such agreement suggests that the differentiation between North Atlantic and South Atlantic swordfish is the consequence of population subdivision and low levels of gene flow.
The shallow levels of differentiation between North and South Atlantic using microsatellites could be due to allele size homoplasy in combination with large effective sizes [83] of swordfish populations and the high mutation rates characteristic of these type of loci [84]. Alternatively, the failure to find differences using microsatellite analyses within the Atlantic [28][29][30] could be due to the small number of loci (4-10) characterized in those studies, or because the samples employed to represent the North Atlantic came from areas identified here as admixture zones within the Northeast Atlantic. Accordingly, a stronger signal of differentiation between South Atlantic and North Atlantic swordfish may be obtained by increasing the number of microsatellite loci aimed to characterize representative samples from reproductive areas in the Northwest Atlantic. Yet, in other highly migratory marine fishes the strong differentiation revealed with mtDNA has not been obtained with microsatellites, and such discordance was interpreted as evidence of male mediated gene flow in bigeye tuna (Thunnus obesus) [16], bluefin tuna (T. thynnus) and yellowfin tuna (T. albacares) [85], sandbar shark (Carcharhinus plumbeus) [86], and shortfin mako (Isurus oxyrinchus) [15], without considering alternative hypotheses. In blue marlin (Makaira nigricans) mitochondrial cytochrome b gene sequences [87] and RFLP data of the mtDNA molecule [88] revealed a pronounced signal of inter-oceanic differentiation that was considerably higher compared to allozyme and scnDNA data [89]. Such disparity was attributed to the enhanced effect of genetic drift under migration-drift equilibrium resulting from the fourfold lower effective population size of mitochondrial than nuclear molecules, and while the effect of male biased inter-ocean exchange could not be rejected, it was not required to account for the observed disparity. Computer simulations and power analyses conducted later on the blue marlin allozyme and scnDNA data together with newly generated microsatellite and whole molecule mtDNA RFLP data [90], indicated that the large range of F ST values among marker classes and among loci within a marker class were consistent with the same population level N e m, given differences in both mode of inheritance and genomic sampling variance among estimates. Further, one microsatellite locus (Mn8) yielded a pattern concordant with mtDNA regarding both the strong inter-oceanic differentiation, and inferred unidirectional gene flow of Indo-Pacific blue marlin into the South Atlantic. In swordfish, the similarity in levels of differentiation between North and South Atlantic populations using mtDNA (F ST = 0.0244; [17]) and multilocus nuclear SNPs (F ST = 0.0393; this study), is not consistent with gender-biased gene flow. Considering that both swordfish and blue marlin females attain larger sizes than males and are more abundant at higher latitudes [81,91,92] where inter-oceanic gene flow would be facilitated (e.g., around southern Africa), it would be highly unlikely that male biased dispersal explains the observed inter-oceanic allele frequency discrepancies among nuclear and mitochondrial markers in both species. Further, in the Atlantic the differences between the sex ratios at-size seem to support the hypothesis that mature swordfish females are more mobile, tending to cover (in terms of probability) broader geographic areas to sustain the energetic requirements required to support the high production of eggs, whereas the migration of males would be generally more restricted and conditioned by the lower average body biomass and the lower energy requirements in comparison to the females to produce spawning products [92]. Thus in spite of the higher mobility and gene flow potential of females no gender-biased gene flow was detected in Atlantic swordfish.

Atlantic swordfish population subdivision and admixture
Previous studies of Atlantic swordfish population structure only broadly identified potential areas of admixture between the North Atlantic and Mediterranean populations [7,31] and between the North Atlantic and South Atlantic populations [17,18,26] and offered no estimates of the contribution of each population. The Strait of Gibraltar is the current management boundary utilized by ICCAT to separate the Mediterranean and the North Atlantic populations. Multilocus analyses of SNPs with STRUCTURE ( Fig 5) and GENELAND (Fig 6C) suggests that the Mediterranean stock extends past the management boundary to around 10°W into North Atlantic waters. However, given that the mapping algorithm of GENELAND extrapolates posterior probabilities into areas of limited or no sample coverage some of the posterior probability curves should be viewed with caution. Specifically, since no samples off the northern coasts of Portugal or Spain were characterized, the posterior probability curves in this region should not be interpreted as evidence of the Mediterranean population extending around the Iberian Peninsula into the Bay of Biscay and the Celtic Sea (Fig 6C). Taken all together, Mediterranean swordfish exiting through the Strait of Gibraltar appear to remain primarily in the shallower Atlantic shelf waters off the northwest African and Iberian coasts, with occasional migrants venturing farther into North Atlantic waters, perhaps as far as the Azores [30, this study]. Movements of mature swordfish inhabiting the shelf area west of the Strait of Gibraltar into the Mediterranean, putatively associated with spawning behavior around the Balearic Islands, occur from May to the beginning of July [82]. Later in the year, between August-November a migration associated with feeding behavior in the opposite direction (east to west) across the Strait takes place with no reproductively active swordfish and with 15% of the catch corresponding to juveniles [82]. Given that our samples are limited to summer and fall when migrations are believed to take place, the expected reduction in frequency of Mediterranean swordfish in this Atlantic region during the winter and spring cannot be confirmed at this time. In a previous study using mtDNA CR-I sequence data, Alvarado Bremer et al. [93] documented no Atlantic swordfish east of the Strait of Gibraltar. Similarly, in this study, no Atlantic swordfish were identified within the Mediterranean Sea (n = 59). It should be noted that in this study Mediterranean swordfish were only sampled in the western Mediterranean, specifically in the area east of the Alboran Sea and south of the Balearic Islands, where the presence of swordfish juveniles has been documented [94] and reproductive behavior is expected to occur [82]. If North Atlantic swordfish were to enter the Mediterranean, their presence would be best documented in this region. Further, the low frequency of movements of Atlantic swordfish into the Mediterranean is also supported by the Bayesian analysis of microsatellite data from 602 western and eastern Mediterranean swordfish in which only three individuals of putative Atlantic origin were identified [30]. Altogether it appears that oceanographic features and or behavioral cues limit the extent of admixture of Atlantic and Mediterranean swordfish populations, with the majority of the contact confined to the shallower waters of both the northwest African shelf and the Gulf of Cadiz.
While the Strait of Gibraltar and associated oceanographic and geophysical features may constrain the admixture between North Atlantic and Mediterranean populations, no such physical barrier separates the populations of North Atlantic and South Atlantic swordfish with admixture occurring over a considerably broader geographic area. A proposal to move the current management boundary separating the North Atlantic from the South Atlantic populations from 5°N to about 15°N was suggested by Chow et al. [26] on the basis of significant differences in allele frequency of CaM. Our analyses with PCoA (Fig 4), STRUCTURE (Fig 5), and GENE-LAND (Fig 6) are in agreement with the extension of the South Atlantic swordfish stock past the current 5°N boundary but suggest an extension even farther north, past Cape Verde (12) to around 20°N. This boundary does not conforms to a horizontal line extending from the Caribbean Sea to the African coast. Instead, the 0.90 isoline of GENELAND is more akin to a step extending from about 0°N 45°W north to 25°N 45°W, and from there almost horizontally towards the African coast (Fig 6). To the north and west of this boundary, an admixture zone begins around Western Sahara (11) and continues north to Morocco (10) and the Iberian sea (7) extending west towards the central North Atlantic (6) and then south towards the northern coast of Brazil. An interesting feature of the pattern of admixture between northern and southern stocks is that the relative proportion contribution to neighboring populations is asymmetrical, with South Atlantic swordfish moving past this boundary towards the north, but with no swordfish of North Atlantic origin moving into the range of the South Atlantic population, reminiscent of the asymmetric contact between the North Atlantic and Mediterranean swordfish populations.
Posterior probability membership contours generated by GENELAND (Fig 6) could be used for mixed stock allocation in the identified areas of stock admixture, with the caveat that the sensitivity of spatial Bayesian clustering, particularly when identifying linear boundaries, is contingent upon sampling coverage [95] as well as initial model assumptions. However, altering initial model assumptions did not alter the general pattern of genetic differentiation reported (not shown). Seasonal displacements of North and South Atlantic swordfish, as well as sampling in regions not properly represented here (e.g., 20-25°N and 30-40°W), may change the extent and geographical placement of admixture, but are not expected to alter the general pattern of differentiation separating these two populations (see notes below regarding reduced abundance in these areas derived from catches).

Potential barriers to gene flow
In the absence of a geographic barrier to explain the asymmetric pattern of admixture separating north and south populations of Atlantic swordfish, biological, behavioral, and oceanographic impediments to gene flow need to be examined. Conventional tag-recapture data [20,96], though biased in its coverage to primarily the Northwest Atlantic fishery, indicates that swordfish tagged west of the Azores generally follow similar patterns of directional migration, west by southwest, back to the primary breeding grounds of the North Atlantic indicative of spawning and feeding ground fidelity, whereas swordfish tagged east of the Azores (of which there are far fewer records) tend to move east by southeast. Interestingly, of the few swordfish tagged between 15°N and 30°N, movement generally leads east by northeast. While conventional tagging data may mirror some of the same molecular distribution trends inferred in this study, these data are heavily biased to areas of high fishing effort [20] and the reporting tag-recapture rates among fleets. Although limited by their relative short duration, Pop-up satellite archival tag (PSAT) data has yielded similar evidence of seasonal spawning and feeding ground migrations in the Northwest Atlantic [22,23].
Recently, Abascal et al. [24] reported the PSAT tracks for 10 Atlantic swordfish tagged in the North Atlantic as part of a tagging study aimed to characterize diel vertical diel behaviors and horizontal seasonal migrations and habitat preferences of this species. Seven tracks corresponded to swordfish tagged around 35°N-45°N in latitude and 40°-50°W in longitude. The other three tracks were of swordfish tagged in a southeastern area of the North Atlantic to the north of Cape Verde archipelago, around 20°N-23°N in latitude and 20°W-25°W in longitude. The movements of these 10 swordfish are concordant with the results obtained with STRUC-TURE, with four tracks indicating movements within areas (two tracks each) identified here as North Atlantic (localities 3 and 5) or South Atlantic (locality 12) populations, whereas six swordfish moved in and out of sampling areas where the posterior probabilities suggest admixture zones (i.e., locality 6, for the fish tagged in the northern area, and 6, 7 and 11, for fish tagged in the southern area). However, one PSAT deployed in the northern area after spending several months in the central North Atlantic admixture zone, spent the last six months of a year at liberty within locality 12 (Cape Verde), identified here as an area belonging to the South Atlantic population. Altogether, only one of 10 swordfish tracked were in conflict with the probability assignments (>0.95) of STRUCTURE. By contrast, concordance of the PSAT tracks with GENELAND's North-South Atlantic population boundaries (0.10-0.90 confidence limit) occurred in seven out of ten swordfish. Two swordfish tagged in the northern area, crossed the 0.90 isoline of assignment to the South Atlantic population. Conversely, one of three tracks of swordfish tagged in the southern area crossed the 0.90 confidence isoline delimiting the North Atlantic population. However, none of the swordfish tagged either in the northern or southern areas of the North Atlantic crossed into the respective 1.0 probability isolines delimiting the North Atlantic and South Atlantic swordfish populations, respectively.

Within-ocean differences in patterns of genetic differentiation
Strong natal fidelity, with admixture primarily confined to feeding grounds, has been proposed as an explanation for the unparalleled intra-oceanic differentiation separating North Atlantic and South Atlantic swordfish populations [17]. By contrast, within the Pacific Ocean [97] and along the Mediterranean Sea [31] differentiation of mtDNA appears to conform to isolation by distance (IBD), although such a pattern in the Pacific Ocean is not supported by the multilocus analysis of nDNA SNP data [98]. Conversely, in the Indian Ocean the presence of at least two potential genetic populations was suggested [18,99,100], but this interpretation was rejected recently [79] in favor of panmixia since the frequency distribution of mitochondrial and microsatellite markers using a more ample sampling coverage of that basin including multiple seasons was homogeneous. Further Muths et al. [79] propose that the presence of distinct swordfish populations respectively within the Atlantic and Pacific Oceans could be partly explained by the existence in those basins of two separate hemispheric oceanographic systems with corresponding discrete feeding areas for the respective populations with limited or no admixture. The association of genetic differentiation with the existence of separate feeding areas as proposed by Muths et al. [79] contrasts dramatically with the interpretation of spawning site fidelity as the driving force behind the observed genetic heterogeneity of swordfish populations. Accordingly, maximum levels of genetic differentiation using both mtDNA and nuclear SNP data were obtained when separate breeding areas were compared [17]. In the current study restricting the analysis to spawning areas yielded higher F ST values for all loci while at the same time increased the H e values for several loci (e.g. Act2α, CaM, Mlc2). This result supports the hypothesis that samples from feeding grounds in both the NW Atlantic and the South Atlantic Ocean (i.e., Gulf of Guinea) may be more prone to population admixture as suggested by Alvarado Bremer et al. [17]. Similarly, spawning site fidelity is also supported by the evidence of the admixture of North Atlantic and Mediterranean swordfish populations in the feeding area to the west of the Strait of Gibraltar by both mtDNA [25] and nDNA data (this study). The concordance of mtDNA and nDNA data was also reported for the separation of South Atlantic from Indian Ocean swordfish [79]. Such concordance of nuclear and mtDNA data in the separation of swordfish populations inhabiting the Mediterranean, North Atlantic, South Atlantic and Indian Ocean basins, is particularly relevant given reported differences in rates of dispersal of swordfish males and females. While fidelity to spawning and feeding grounds may explain low levels of gene flow between North and South Atlantic swordfish, it doesn't fully explain why admixture only seems to occur in the Northeast Atlantic waters east of about 50°W and south of 30°N. Chow et al. [26] suggested that changes in oceanographic parameters (e.g. salinity, dissolved oxygen, temperature, etc.) in deep waters might explain the strong genetic break that their CaM data revealed since swordfish are deep divers. Nutrient upwelling, caused by prevailing winds off the western coast of Africa, creates a large hypoxic zone in the eastern tropical Atlantic [101] that apparently correlates with the inferred area of population differentiation identified in this study (see supporting S3 Fig). Sensitivity to oceanographic conditions in deep waters has been documented in other billfish [102] and swordfish are known to dive substantially deeper than other billfish [21], and thus may be more sensitive to changes in concentrations in DO at depth, which might be used as a cue by North Atlantic swordfish not to move south to the areas occupied by South Atlantic swordfish.
Other oceanographic features are known to acts as barriers to dispersal in proximity to coastal areas. For instance while recognizing that the Amazon discharge would not be expected to affect the movements of adult highly migratory pelagic fish, Chow et al. [26] hypothesized it may act as a freshwater barrier to passive larval dispersal between the Caribbean and Brazilian provinces. This feature geographically coincides with the posterior probability break between North Atlantic and South Atlantic populations identified by STRUCTURE and GENELAND, as well as an inferred reported break in potential spawning areas [17] (Fig 7). Further, the same region coincides with areas where historical catches of adult swordfish are low, and if catches properly reflect abundance (Fig 6), then the low density of adults in this region together with the Amazon discharge preventing the admixture of larvae from northern and southern breeding areas, as hypothesized by Chow et al. [26] could together promote the observed patterns of genetic differentiation between northern and southern populations of swordfish. Finally, the population structure of Atlantic swordfish may be influenced by sea surface currents as reported in other pelagic species (e.g. blue marlin [103]). The North Atlantic gyre, the Canary current, and the Atlantic equatorial counter current all associate with the population differentiation of North Atlantic and South Atlantic swordfish, and may explain why population admixture is primarily confined to Northeast Atlantic waters.

Conclusions and Future Studies
The findings of this study can be summarized as follows: first, significant differentiation among Mediterranean, North Atlantic, and South Atlantic swordfish populations is confirmed with very low rates of gene flow, particularly when compared to other highly migratory fishes. Second, population admixture between the North Atlantic and the South Atlantic population, and to a much lesser extent with the Mediterranean, is primarily confined to the Northeast Atlantic and in both cases appears to be unidirectional, from the Mediterranean to the NE Atlantic and from the South Atlantic to the NE Atlantic. Third, the genetic partitioning of these populations differs from currently recognized management boundaries. Fourth, and lastly, while the oceanographic and/or behavioral barriers to gene flow are not yet identified, future tagging, genetic, oceanographic, or physiological studies should focus efforts to areas of admixture in the Northeast Atlantic and areas where individuals were assigned with high posterior probabilities, i.e. North and South Atlantic spawning and feeding grounds. Additionally, future genomic studies utilizing large numbers of SNPs may provide insight into whether adaptive genetic variation is responsible for the fine-scale population structure here reported.  [17]. Shaded in dark blue are the Northwest Atlantic and South Atlantic areas where eggs and or early life history stages have been collected. The areas shaded in lighter blue correspond to areas where females with highgonadal indices and/or hydrated oocytes have been collected that coincide to areas of biased male to female sex-ratios [104]. Supporting Information S1 Fig. Bayesian cluster (K) estimation. The estimation of the number of clusters (K) in the STRUCTURE v2.3 [55] analysis of Atlantic and Mediterranean swordfish using (A) the ad hoc approach of Prichard et al. [55] and the mean posterior probability of the data (L(K)) and (B) ΔK approach of Evanno et al. [60]. Twenty independent runs for each K value (1-10) were performed using 100,000 MCMC iterations with a burn in period of 100,000.  Table. SNP allelic data. The alleles for the ten loci and all 774 swordfish analyzed in this study. Allelic numbering corresponds to reported GenBank alleles for all loci except ALDB. Allelic numbering for ALDB corresponds to STR amplicon size scoring. (XLSX) possible. Additional special thanks to L. Janson, Dennis Lee, Jaime Mejuto, Jay Rooker, Bernard Stequert, Marc Griffths, Favio Hazin, Carla Marques, Carles Pla, NMFS-SEFSC, and Insituto Español de Oceanografia for providing swordfish tissue samples for analysis.