# Gene flow at major transitional areas in sea bass (Dicentrarchus labrax) and the possible emergence of a hybrid swarm

Nolwenn Quéré, Erick Desmarais, Costas Tsigenopoulos, Khalid Belkhir, François Bonhomme, Bruno Guinand

<figure><img src="image_1.png" /></figure>

## To cite this version:

Nolwenn Quéré, Erick Desmarais, Costas Tsigenopoulos, Khalid Belkhir, François Bonhomme, et al.. Gene flow at major transitional areas in sea bass (Dicentrarchus labrax) and the possible emergence of a hybrid swarm. Ecology and Evolution, 2012, 2 (12), pp.3061 - 3078. 10.1002/ece3.406. hal-01922698

HAL Id: hal-01922698

https://hal.science/hal-01922698

Submitted on 4 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Gene flow at major transitional areas in sea bass (*Dicentrarchus labrax*) and the possible emergence of a hybrid swarm

Nolwenn Quéré<fnref n="1" />,<fnref n="2" />, Erick Desmarais<fnref n="1" />,<fnref n="3" />, Costas S. Tsigenopoulos<fnref n="4" />, Khalid Belkhir<fnref n="1" />, François Bonhomme<fnref n="1" />,<fnref n="2" /> & Bruno Guinand<fnref n="1" />,<fnref n="2" />

**Keywords**

Hybrid swarm, introgression, microsatellite, minisatellite, sea bass, tension zones.

**Correspondence**

Bruno Guinand, Institut des Sciences de l'Evolution de Montpellier, CNRS-UMR 5554, Université Montpellier 2, cc63, 34095 Montpellier Cedex 5, France. Tel: +33 467 144 687; Fax: +33 467 144 554; E-mail: bruno.guinand@univ-montp2.fr

**Funding Information**

No funding information provided.

Received: 11 July 2012; Revised: 18 September 2012; Accepted: 25 September 2012

Ecology and Evolution 2012; 2(12): 3061-3078

doi: 10.1002/ece3.406

**Abstract**

The population genetic structure of sea bass (*Dicentrarchus labrax*) along a transect from the Atlantic Ocean (AO) to the Eastern Mediterranean (EM) Sea differs from that of most other marine taxa in this area. Three populations (AO, Western Mediterranean [WM], EM) are recognized today, which were originally two allopatric populations. How two ancestral genetic units have evolved into three distinct units has not been addressed yet. Therefore, to investigate mechanisms that lead to the emergence of the central WM population, its current status, and its connectivity with the two parental populations, we applied 20 nuclear loci that were either gene associated or gene independent. Results confirmed the existence of three distinct gene pools, with higher differentiation at two transitional areas, the Almeria-Oran Front (AOF) and of the Siculo-Tunisian Strait (STS), than within any population. Significant linkage disequilibrium and heterozygote excess indicated that the STS is probably another tension zone, as already described for the AOF. Neutrality tests fail to reveal marker loci that could be driven by selection within or among metapopulations, except for locus *DLA0068*. Collectively, results support that the central WM population arose by trapping two tensions zones at distinct geographic locations of limited connectivity. Population assignment further revealed that WM individuals were more introgressed than individuals from the other two metapopulations. This suggests that this population might result from hybrid swarming, and was or is still seeded by genes received through the filter of each tension zone.

**Introduction**

It is well established that, despite the apparent connectiv-
ity of marine habitats and the high dispersal capabilities
of marine organisms, populations may nevertheless be
organized into well-defined genetic units (Hellberg et al.
2002; Hauser and Carvalho 2008; Hellberg 2009; Sanford
and Kelly 2011). If present-day processes can contribute
to patterns of differentiation (Galindo et al. 2006; Gerlach
et al. 2007; Selkoe et al. 2008, 2010), recognized phyloge-
ographic barriers shared by a wide range of species
represent the main barriers to gene flow in the marine
environment (e.g., Benzie 1999; Barber et al. 2000; Patar-

nello et al. 2007; Galarza et al. 2009; Kelly and Palumbi
2010).

Nevertheless, nearly all marine organisms have sporadic
opportunities to cross almost all physical barriers and
to exchange genes between phylogeographically differenti-
ated populations. Opportunities for any neutral or
unconditionally favorable mutation to invade most of –
or the whole – the species range frequently occur after
secondary contacts, so gene pools should be homogenized
on either side because minimal gene migration easily
overcomes mutation and drift over generations. Genetic
differentiation persisting across phylogeographic barriers
probably indicates nonneutral limitations to gene flow, as

<footnote n="1">¹Institut des Sciences de l'Évolution de Montpellier, CNRS-UMR 5554, Université Montpellier 2, cc63, 34095 Montpellier Cedex 5, France</footnote>

<footnote n="2">²Station Méditerranéenne de l'Environnement Littoral, 2 Avenue des chantiers, 34200 Sète, France</footnote>

<footnote n="3">³LabEx CeMEB, Université Montpellier II, place E. Bataillon, cc63, 34095 Montpellier Cedex 5, France</footnote>

<footnote n="4">⁴Hellenic Center for Marine Research, Institute of Marine Biology and Genetics, PO Box 2214, Gournes Pediados, 71500 Heraklion, Crete, Greece</footnote>
exemplified by numerous marine hybrid zones worldwide (e.g., Gardner 1997; Strelkov et al. 2007; Hobbs et al. 2009; Gagnaire et al. 2011). The classical “porous” (or “semipermeable”) genome model states that differentiation would be wiped out from most genomic regions, except around loci conferring fitness advantage to one of the populations or, on the contrary, loci acquired during isolation and involved in endogenous incompatibilities of the Dobzhansky–Muller type (e.g., Barton and Hewitt 1985; Mallet 2005; Payseur 2010; Bierne et al. 2011). Such postzygotic barriers contribute to natural selection against hybrid F₁, F₂, and backcrosses (Nosil et al. 2005; Burton et al. 2006; Marshall et al. 2010), and may reveal costs of adaptation for immigrants (Hereford 2009). They, hence, contribute to the establishment of so-called tension zones reflecting the conflicting effects of dispersal of parental forms and selection against hybrids (e.g., Barton and Hewitt 1985). As the porous genome model implies wide heterogeneity in the distribution of differentiation along chromosomes, estimation of population divergence will depend on the loci investigated, their position relative to sites under selective pressure, and recombination rate (see, e.g., Bierne 2010; Payseur 2010). This leads to different introgression patterns among (classes of loci) across tension zones (e.g., Barton and Hewitt 1985). Conversely – but considerably less studied – new fit hybrid genotypes may arise from successfully recombinant individuals (Harrison 1993; Edmands et al. 2005; Hwang et al. 2011; Willett 2012). A new population (or taxon), designed as a hybrid swarm, can then emerge and successfully occupy a portion of the available habitat (Moore 1977; Arnold 1993; Jiggins and Mallet 2000).

The European sea bass *Dicentrarchus labrax* L. is a model for evolutionary studies among marine fishes, in both wild and cultured populations (Volckaert et al. 2008; Chatain and Chavanne 2009). Over its distribution area, the species is structured into three main metapopulations (Atlantic, Western, and Eastern Mediterranean [EM]), separated by well-defined barriers to gene flow (Benharrat et al. 1983; Naciri et al. 1999; Bahri-Sfar et al. 2000; Lemaire et al. 2005). Genetic homogeneity predominates in Atlantic and Western Mediterranean (WA) basins (e.g., Garcia de León et al. 1997; Naciri et al. 1999; Fritsch et al. 2007; Coscia and Mariani 2011), while the less well-studied EM basin is more structured (Bahri-Sfar et al. 2000; Castilho and Ciftci 2005). The genetic transitions over the range are each coincident with an oceanic front (the Almeria-Oran front [AOF] and the Siculo-Tunisian Strait [STS]). This may reflect tension zones that have been trapped there, with low dispersal rates and somewhat reduced population densities (Barton’s trapping hypothesis; discussed in Bierne et al. 2011). Evidence that these transitional areas are tension zones in sea bass

has however only been documented for the AOF by Lemaire et al. (2005). Authors demonstrated differential introgression at mtDNA and nuclear markers, but also among anonymous nuclear DNA markers (microsatellites) at the AOF. Lemaire et al. (2005) also reported significant linkage disequilibrium (LD) at population samples located close to the AOF, which are classical features of tension zones (Barton and Hewitt 1985).

All studies on sea bass have been based on a limited number of anonymous microsatellites of unknown genomic location. Therefore, the objectives of this study were to use distinct classes of nuclear markers that were either associated or independent of genes, in order to: (1) provide better estimates on the levels of differentiation at major transitional areas; (2) establish connectivity and introgression patterns for each group of markers or patterns shared by markers irrespective the group they belong to; and (3) improve knowledge on the origin of the central WM population of sea bass, in particular whether it is a hybrid swarm.

# Materials and Methods

## Samples and molecular analysis

Sixteen natural populations (*N* = 474) were sampled among the three phylogeographic units of *D. labrax* distribution area: Atlantic Ocean (AO), WM, and EM basins (Table 1). Most samples have been already analyzed in previous studies (Table 1). Nineteen microsatellite and one minisatellite markers were used, classified into two classes of markers: gene-associated loci (hereafter GAL) and gene-independent loci (GIL). A marker was classified as a GAL when it was situated either inside the span of an annotated gene or at <2 kilobases (kb) from one extremity of a coding sequence (CDS). The limit of 2 kb was retained because it generally encompasses the proximal promoter of most genes in vertebrates. Proximal promoters are known to be under functional constraints and, together with their closeness from CDS, mini- and microsatellite loci they contain may have greater probability to be influenced by selection (e.g., Li et al. 2004; Gemayel et al. 2010). A symmetrical size limit of 2 kb was also retained after the 3'-end of the CDS while – as far as we know – documented functional constraints are less reported on this side of genes. Those limits are, however, operating limits to classify GAL and GIL; they do not mean that loci located farther from CDS cannot be impacted by selection (and vice versa). The GAL comprises one mini- and three microsatellite markers that have been specifically designed for their association with growth hormone (GH: *µGH* and *mGH*), somatolactin (SL: *µSL*), and insulin-like growth factor (IGF-1: *µIGF-1*);
<table><caption>**Table 1.** Samples used in this study.</caption><thead><tr><th>Basin</th><th>Abbreviation</th><th>Population</th><th>Country</th><th>Number of individuals</th><th>Geographic location</th><th>Origin of sample</th></tr></thead><tbody><tr><td>AO (31)</td><td>BB</td><td>Bay of Biscay</td><td>France</td><td>31</td><td>45°44'19"N 1°25'11"O</td><td>Fritsch et al. (2007)</td></tr><tr><td rowspan="6">WM (365)</td><td>GOU</td><td>La Goulette</td><td>Tunisia</td><td>30</td><td>36°50'06"N 10°19'59"E</td><td>Guinand et al. (2008)</td></tr><tr><td>ICH</td><td>Ichkeul</td><td>Tunisia</td><td>49</td><td>37°10'21"N 9°39'55"E</td><td>Bahri-Sfar et al. (2000)</td></tr><tr><td>MAR</td><td>Marsala</td><td>Italia</td><td>24</td><td>37°48'09"N 12°25'11"E</td><td>Naciri et al. (1999)</td></tr><tr><td>ANN</td><td>Annaba</td><td>Algeria</td><td>24</td><td>36°54'47"N 7°46'36"E</td><td>Naciri et al. (1999)</td></tr><tr><td>SBD</td><td>Sabaudia</td><td>Italia</td><td>16</td><td>41°19'36"N 13°59'30"E</td><td>Lemaire et al. (2000)</td></tr><tr><td>FIUM</td><td>Fiumicino</td><td>Italia</td><td>22</td><td>41°51'36"N 12°70'02"E</td><td>Lemaire et al. (2000)</td></tr><tr><td rowspan="7">EM (98)</td><td>OR</td><td>Etang de l'Or</td><td>France</td><td>54</td><td>43°34'35"N 4°01'40"E</td><td>This study</td></tr><tr><td>PER</td><td>Pérols</td><td>France</td><td>30</td><td>43°33'50"N 3°58'10"E</td><td>This study</td></tr><tr><td>LUN</td><td>Lunel</td><td>France</td><td>30</td><td>43°36'05"N 4°05'12"E</td><td>This study</td></tr><tr><td>SET</td><td>Sête</td><td>France</td><td>38</td><td>43°23'06"N 3°42'19"E</td><td>Guinand et al. (2008)</td></tr><tr><td>MUR</td><td>Murcia</td><td>Spain</td><td>48</td><td>37°51'50"N 0°43'31"O</td><td>Lemaire et al. (2005)</td></tr><tr><td>SYR</td><td>Syria</td><td>Syria</td><td>32</td><td>35°34'50"N 35°31'36"E</td><td>This study</td></tr><tr><td>CYP</td><td>Cyprus</td><td>Cyprus</td><td>15</td><td>34°33'03"N 32°58'36"E</td><td>Boutet et al. (2008)</td></tr><tr><td></td><td>SYRA</td><td>Syracuse</td><td>Italia</td><td>26</td><td>37°47'07"N 18°48'27"E</td><td>This study</td></tr><tr><td></td><td>SEL</td><td>Selinunte</td><td>Italia</td><td>25</td><td>36°50'53"N 13°34'07"E</td><td>This study</td></tr></tbody></table>

AO, WM, and EM indicate the basin from where samples originate (AO: Atlantic Ocean; WM: Western Mediterranean; EM: Eastern Mediterranean). Numbers of individuals considered in each basin are given in brackets. References of published work that already used individual samples are reported.

Quéré et al. 2010; description and conditions herein) and seven microsatellite loci from two simple sequence repeat (SSR)-enriched genomic libraries (Tsigenopoulos et al. 2003; Chistiakov et al. 2004) that appeared to be close to CDS afterward. The GIL class contains nine loci originating from the same SSR-enriched genomic libraries than GAL (Table 2).

Markers from SSR libraries were amplified in two multiplexes PCR (Table 2). For the first multiplex, amplification conditions have been reported in Guinand et al. (2008). The second multiplex was optimized for a final volume of 10 µL with 10× Taq buffer, 30 mmol/L MgCl₂, 2.7 mmol/L dNTP, 10 ng of genomic DNA and 4 µmol/L of each primer. Reactions were performed on a PTC-200 (MJ Research, St Bruno, Québec, Canada) as follows: an initial denaturation at 95°C for 3 min, 35 cycles at 94°C for 45 s, annealing at the 58°C for 45 s, and 72°C for 45 s, followed by a final elongation at 72°C for 10 min. Genotyping was performed on ABI PRISM 3130xl or 3700 DNA Analyser (Life Technologies, St Aubin, France), using 5'-labeled reverse primers and the Gene-Scan® TM-500 LIZ Size Standard (Life Technologies) as internal size standard.

## Data analyses

After control for null alleles with MicroChecker (Van Oosterhout et al. 2004), all analyses – except neutrality tests – were first performed with individual samples (Fig. 1; Table 1), then with pooled samples corresponding to the three metapopulations. In order to provide descriptive

summary statistics, the mean number of alleles (Ā), the mean unbiased expected gene diversity (H̄exp; Nei 1987) were computed for each locus and each (meta)population. Significant within or among basins differences in Ā and H̄exp were tested using the nonparametric Mann-Whitney U-test implemented in the Microsoft Excel add-on, StatEL. Complementarily, genetic description of sample was assessed by computing LD. LD between pairs of loci was computed using GENETIX v4.05 (http://genetix.univ-montp2.fr/) in each sample, and tested using permutation tests (5000 permutations). Finally, deviations from Hardy-Weinberg expectations (HWE) within each original sample and within each metapopulation were investigated using f̂ (Weir and Cockerham 1984), also using GENETIX. The null hypothesis (f = 0) of no significant departure from panmixia was tested by randomly permuting alleles (n = 1000) from the original matrix of genotypes.

Among-sample comparisons were assessed by estimating levels of population differentiation using ò, an estimator of FST (Weir and Cockerham 1984). Estimates of ò were independently computed for three data sets (i.e., the full set of 20 loci, then GAL and GIL sets). Their significance was tested according to the permutation procedure available in GENETIX and corrected for multiple tests according to Benjamini and Yukitieli (2001) false discovery rate (FDR; Narum 2006). Distance trees were inferred from Reynolds et al. (1983) coancestry genetic distance matrix [D = 1 - ln(1 - ò)] by the Neighbor procedure of the phylogenetic package Phylip 3.6 (Felsenstein 2005). Trees were visualized using Treeview (Page 1996). Following
<table><caption>**Table 2.** Microsatellite loci used in this study, including size ranges, GenBank accession numbers, fluorochrome labeling, and PCR primers.</caption><thead><tr><th>Locus</th><th>Category</th><th>Annotation</th><th>Motif</th><th>GenBank Accession</th><th>Size range (bp)<sup>1</sup></th><th>Fluorochrome</th><th>Primers (5'→3')</th></tr></thead><tbody><tr><td>DLA0041</td><td>GIL</td><td>—</td><td>(TC)<sub>8</sub>(AC)<sub>21</sub>AA(AC)<sub>3</sub>AA(AC)<sub>3</sub></td><td>DQ363864</td><td>167–201</td><td>FAM</td><td>F AAAAGGAACAGCCCTCCAC<br>R AGCATTGTTCTTCTGAGTGACC</td></tr><tr><td>DLA0044</td><td>GAL</td><td>Unknown</td><td>(TC)<sub>18</sub>TTTT(TC)<sub>6</sub>CTCC</td><td>DQ363867</td><td>105–129</td><td>HEX</td><td>F TCCGCTCCGCACCGAGTGAC<br>R ACCGCCCAAGGGTTGGACTG</td></tr><tr><td>DLA0051</td><td>GAL</td><td>MAPK3<sup>a</sup></td><td>(GT)<sub>16</sub></td><td>DQ363874</td><td>149–181</td><td>ROX</td><td>F AGGTTCTTGGCCTGGGAATC<br>R AGTGACAGCAGCCTCCAGAG</td></tr><tr><td>DLA0060</td><td>GAL</td><td>Bestrophin 3<sup>b</sup></td><td>(CA)<sub>12</sub>(TA)<sub>3</sub>AA(CA)<sub>2</sub></td><td>DQ363883</td><td>119–131</td><td>FAM</td><td>F GAGAGTTCATCCTGTTCGCTC<br>R TGTAGTAATAATGCGCTCTGCAA</td></tr><tr><td>DLA0061</td><td>GIL</td><td>—</td><td>(TG)<sub>14</sub></td><td>DQ363884</td><td>153–167</td><td>FAM</td><td>F AAAGGCCAGTGAAACTCATGT<br>R CTCCCTGTCCATCTGTCCTC</td></tr><tr><td>DLA0066</td><td>GIL*</td><td>—</td><td>(AG)<sub>22</sub></td><td>DQ363887</td><td>135–161</td><td>PET</td><td>F GTTGACCGGAGTCCTAGC<br>R GGCCATATGTGTCTTGCTT</td></tr><tr><td>DLA0068</td><td>GIL</td><td>—</td><td>(CA)<sub>7</sub>CGCACG(CA)<sub>3</sub></td><td>DQ363889</td><td>247–266</td><td>NED</td><td>F CAACACCTGTTCCTCTGAACC<br>R GCATTAGCATTGATTGTCCTG</td></tr><tr><td>DLA0070</td><td>GAL</td><td>Unknown</td><td>(AC)<sub>30</sub></td><td>DQ363891</td><td>126–155</td><td>VIC</td><td>F TCTGCTTGCATCTGTGGAAT<br>R GCCATCTGGCTAGCTTCACT</td></tr><tr><td>DLA0073</td><td>GIL*</td><td>—</td><td>(CT)<sub>36</sub></td><td>DQ363894</td><td>157–181</td><td>NED</td><td>F CATGACTTCATGTGCTAATGTCC<br>R AGTTCAGAGCGGCAACTGT</td></tr><tr><td>DLA0075</td><td>GIL</td><td>—</td><td>(CA)<sub>15</sub></td><td>DQ363896</td><td>180–186</td><td>PET</td><td>F CACATACACAAGCTTAACCC<br>R GGCAGAGATGGGAAATAGACA</td></tr><tr><td>DLA0078</td><td>GAL</td><td>MAN2A1<sup>c</sup></td><td>(AG)<sub>29</sub></td><td>DQ363899</td><td>223–243</td><td>VIC</td><td>F AAGACTGGACCTCTGGAGACC<br>R CACAAGGAACCGAGACAAGA</td></tr><tr><td>DLA0081</td><td>GAL*</td><td>PPP2R5A<sup>d</sup></td><td>(CA)<sub>16</sub></td><td>DQ363902</td><td>202–222</td><td>PET</td><td>F GACGAAGACTTCAGACGAGCTAT<br>R ATACCGAGCGACCATGTTG</td></tr><tr><td>DLA0086</td><td>GIL</td><td>—</td><td>(AC)<sub>26</sub></td><td>DQ363907</td><td>191–205</td><td>FAM</td><td>F GCTAGAGGATTCATGTCGCTT<br>R ACCTGGTGATTGGCAATTCT</td></tr><tr><td>DLA0089</td><td>GAL</td><td>LLGL1<sup>e</sup></td><td>(GT)<sub>15</sub></td><td>DQ363910</td><td>127–135</td><td>NED</td><td>F ACGAGTAATGAGGACCCA<br>R GTCAAAACAGCCCACCTA</td></tr><tr><td>DLA0096</td><td>GIL*</td><td>—</td><td>(TG)<sub>16</sub></td><td>EF471091</td><td>256–270</td><td>ROX</td><td>F AACTTAGTGAAGTAACTTGTGGCAA<br>R TCGATGCATCTAGGACAGGA</td></tr><tr><td>DLA0097</td><td>GIL</td><td>—</td><td>(GT)<sub>4</sub>GC(GT)<sub>2</sub>GC(GT)<sub>13</sub></td><td>EF471092</td><td>216–240</td><td>HEX</td><td>F GCTGCAGGAGTGTGAGAGG<br>R GCGAGAGACTCGAGGAAGA</td></tr></tbody></table>

Categories of loci to which markers belong to are reported, together with their functional annotation for GAL when available. Note that GAL developed by Quéré et al. (2010; µGH, mGH, µSL, and µIGF-1) associated with growth hormone (GH), somatolactin (SL), and insulin-like growth factor 1 (IGF-1) are not reported in this table. Loci marked with an asterisk indicated loci that did not amplify in the EM samples, and were analyzed only in the AO and WM samples. Alphabetical superscripts refer to selected references in which functional role of annotated GAL markers are reported.

O'Reilly et al. (2004), we explored the relationship between $\hat{\theta}$ and $H_{exp}$ to check for effects of size homoplasy. We specifically controlled whether GIL and GAL markers were differently affected. We further explored the possible role of selection in shaping differentiation in the full, GAL, and GIL data sets using neutrality tests developed by Beaumont and Nichols (1996; hereafter B&N and implemented in LOSITAN Antao et al. 2008, http://popgen.eu/soft/lositan/), by Vitalis et al. (2001; implemented

in DETSEL Vitalis et al. 2003, http://www.genetix.univ-montp2.fr/detsel.html), and by Foll and Gaggiotti (2008; hereafter: F&G and implemented in BayeScan, http://www-leca.ujf-grenoble.fr/logiciels.htm). All tests were first performed among populations within each basin, then among the three basins except for DETSEL, specifically designed for pairwise comparisons then between each pair of basins. Outlier loci detected using several methods would be stronger candidate to the

<footnote n="1">¹Size range established before this study.</footnote>

<footnote>*Locus that was not genotyped in the EM basin.</footnote>

<footnote>ªGarcía et al. (1998); Kultz (1998).</footnote>

<footnote>ᵇO'Driscoll et al. (2008); Matchkov et al. (2008).</footnote>

<footnote>ᶜIshida et al. (1999).</footnote>

<footnote>ᵈJanssens and Goris (2001).</footnote>

<footnote>ᵉWodarz (2000), Vasioukhin (2006).</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Map of sampling locations; samples are abbreviated as in Table 1. The colored symbols for samples belonging to each metapopulation will be used in other figures. Location of the Almeria-Oran Front and the Siculo-Tunisian Strait are reported.</figcaption></figure>

action of selection than loci detected once. The B&N and, when considered, the DETSEL tests were both performed according to *P* = 0.95 and *P* = 0.99 successively. F&G test considered a Log₁₀(BF) of 1.3 and 2 (BF = Bayes Factor) corresponding to posterior probabilities of locus effects of 0.95 and 0.99, respectively.

The program MIGRATE-n v3.0 (Beerli 2008) was used to infer *M*, the migration rate, among basins only (*M* = *m*/*μ*, where *m* is the immigration rate per generation). Contrary to most *F*<sub>ST</sub>-based statistics, MIGRATE-n provides estimation of asymmetric migrations among (meta)populations, indicating if a (meta)population is a net donor or net receiver of individuals over evolutionary times. It hence provides a picture how tension zones may function. These calculations used the Brownian mutation model and the mutation rate was considered equal for all loci (μ = 10<sup>-4</sup>) even for GAL. We used coalescent maximum likelihood (ML) based on Markov Chain Monte Carlo with Hastings Metropolis importance sampling to infer the various parameters (Beerli and Felsenstein 1999, 2001). *F*<sub>ST</sub> estimates among basins were used as initial parameters for the estimation of Θ and *M* in MIGRATE-n. For each locus, the ML was run for 10 short and five long chains with 50,000- and 10,000-recorded genealogies, respectively, after discarding the first 1000 genealogies (burn-in) for each chain. One of every 20 reconstructed genealogies was sampled for both the short and long chains. We used an adaptive heating scheme with four concurrent chains. Analyses were all performed in replicates either on the full set of loci or on the GAL and GIL sets independently.

As a complement to estimates of population differenti-
ation and asymmetric migration rates, we computed the
probability of membership of individuals to each
metapopulation (i.e., individuals should have higher
probability of membership in the metapopulation they
were sampled). Assignments of individuals to populations

and associated probabilities were inferred with the soft-
ware STRUCTURE v2.3 (Pritchard et al. 2000; http://
pritch.bsd.uchicago.edu/structure.html) by setting
K = 3 (i.e., three basins) after verification of the most likely
number of independent population clusters was K = 3
using Evanno et al.'s (2005) ΔK method (details not
reported). The software used a Monte Carlo Markov
Chain (MCMC) Bayesian clustering method that maxi-
mizes the within-cluster Hardy-Weinberg and linkage
equilibriums. The admixture model with noncorrelated
allele frequencies was used for the full, GAL, and GIL
data sets. A burn-in length of 50,000 iterations and subse-
quent 500,000 additional MCMC iterations were carried
out. Individuals were assigned to clusters based on the
highest probability of membership (q̂-value). Five repli-
cates were independently performed, giving reproducible
results. Average q̂-values (±1 standard deviation) of indi-
viduals to each metapopulation were hence computed
from those individual assignments.

Analyses using MIGRATE-n and STRUCTURE were first performed using initial sample sizes of the three basins (Table 1), but, as sample sizes were unbalanced, we also performed analyses with 31 individuals in each basin (i.e., the number of individuals from the AO basin). Results were highly reproducible and qualitatively comparable to those obtained with the full set of individuals and will not be reported.

Results

Numerous amplifications failures were observed for fishes
from EM populations, at several loci (GAL: *DLA0081*,
*mGH*; GIL: *DLA0066*, *DLA0073*, *DLA0096*). This was
probably due to loss of sequence homology at PCR prim-
ing sites. Those loci were discarded from some analyses,
as relevant. Therefore, genetic comparisons were per-
formed with 20 loci (11 GAL, nine GIL), except those
involving EM populations, which were only based on 15
loci (nine GAL, six GIL).

Gene and allelic diversities

Allelic and gene diversities are reported in Table 3. There
were no differences among populations within the WM
or the EM basin for the whole set of loci, or for the GAL
and GIL marker sets considered separately (details not
reported), (20 or 15 loci, respectively). Â and H<sub>exp</sub> were
not significantly different between the AO and WM meta-
populations (*P* > 0.05 in each case; 20 loci). In contrast,
Â and H<sub>exp</sub> of the EM metapopulation were significantly
lower than the other metapopulations for the 15 loci
(*P* < 0.001). When considering the GAL and GIL marker
sets independently, Â was significantly larger for GAL
<table><caption>**Table 3.** Gene and allelic diversity indices for each individual samples and metapopulations (AO, EM, WM) of sea bass considered in this study.</caption><thead><tr><th rowspan="2">Basin</th><th rowspan="2">Sample</th><th colspan="3">Multilocus</th><th colspan="3">GAL</th><th colspan="3">GIL</th></tr><tr><th>H<sub>exp</sub></th><th>A<sub>r</sub></th><th>f</th><th>H<sub>exp</sub></th><th>A<sub>r</sub></th><th>f</th><th>H<sub>exp</sub></th><th>A<sub>r</sub></th><th>f</th></tr></thead><tbody><tr><td rowspan="12">AO<br>WM</td><td>GG</td><td>0.77</td><td>8.70</td><td>0.01</td><td>0.78</td><td>9.33</td><td>-0.01</td><td>0.76</td><td>7.75</td><td>0.03</td></tr><tr><td>MUR</td><td>0.79</td><td>9.65</td><td>0.01</td><td>0.79</td><td>10.17</td><td>0.01</td><td>0.78</td><td>8.88</td><td>0.01</td></tr><tr><td>OR</td><td>0.77</td><td>9.20</td><td>0.00</td><td>0.76</td><td>9.48</td><td>0.04**</td><td>0.79</td><td>8.73</td><td>0.06***</td></tr><tr><td>SET</td><td>0.77</td><td>7.84</td><td>0.01</td><td>0.77</td><td>8.64</td><td>-0.01</td><td>0.78</td><td>6.75</td><td>0.03</td></tr><tr><td>PER</td><td>0.77</td><td>8.61</td><td>-0.01</td><td>0.76</td><td>9.17</td><td>-0.01</td><td>0.77</td><td>7.75</td><td>0.00</td></tr><tr><td>LUN</td><td>0.78</td><td>8.90</td><td>0.02</td><td>0.78</td><td>8.58</td><td>0.00</td><td>0.79</td><td>9.38</td><td>0.05</td></tr><tr><td>FIUM</td><td>0.75</td><td>7.65</td><td>0.00</td><td>0.74</td><td>8.25</td><td>-0.04</td><td>0.76</td><td>6.70</td><td>-0.02</td></tr><tr><td>SBD</td><td>0.77</td><td>7.36</td><td>-0.03</td><td>0.74</td><td>7.25</td><td>-0.04</td><td>0.81</td><td>7.53</td><td>-0.02</td></tr><tr><td>MAR</td><td>0.76</td><td>8.13</td><td>-0.01</td><td>0.76</td><td>8.06</td><td>-0.02</td><td>0.76</td><td>8.26</td><td>0.01</td></tr><tr><td>ANN</td><td>0.78</td><td>9.47</td><td>0.00</td><td>0.79</td><td>9.54</td><td>-0.02</td><td>0.77</td><td>9.38</td><td>0.03</td></tr><tr><td>ICH</td><td>0.74</td><td>8.12</td><td>0.02</td><td>0.73</td><td>7.83</td><td>-0.01</td><td>0.76</td><td>8.48</td><td>-0.07***</td></tr><tr><td>GOU</td><td>0.77</td><td>8.54</td><td>-0.03**</td><td>0.76</td><td>8.83</td><td>-0.03**</td><td>0.77</td><td>8.13</td><td>-0.01</td></tr><tr><td>WM metapop.</td><td>0.77</td><td>8.49</td><td>0.01</td><td>0.76</td><td>8.71</td><td>-0.01</td><td>0.78</td><td>8.18</td><td>0.01</td></tr><tr><td rowspan="5">EM</td><td>SEL</td><td>0.68</td><td>5.24</td><td>-0.07</td><td>0.69</td><td>6.10</td><td>-0.08*</td><td>0.66</td><td>4.19</td><td>-0.05</td></tr><tr><td>SYRA</td><td>0.60</td><td>3.93</td><td>-0.22***</td><td>0.61</td><td>4.20</td><td>-0.24***</td><td>0.56</td><td>3.25</td><td>-0.18</td></tr><tr><td>CYP</td><td>0.61</td><td>4.13</td><td>-0.06</td><td>0.67</td><td>4.48</td><td>-0.06</td><td>0.57</td><td>3.62</td><td>-0.07</td></tr><tr><td>SYR</td><td>0.66</td><td>5.21</td><td>-0.08</td><td>0.67</td><td>5.82</td><td>-0.06</td><td>0.62</td><td>4.03</td><td>-0.06</td></tr><tr><td>EM metapop.</td><td>0.64</td><td>4.62</td><td>-0.07***</td><td>0.67</td><td>5.15</td><td>-0.07***</td><td>0.61</td><td>3.76</td><td>-0.06</td></tr></tbody></table>

AO, Atlantic Ocean; EM, Eastern Mediterranean; WM, Western Mediterranean; GAL, gene-associated loci; GIL, gene-independent loci.

Multilocus estimates of *f* are also provided. \* *P* < 0.05; \*\* *P* < 0.01; \*\*\* *P* < 0.001.

than for GIL in the WM metapopulation (Mann–Whitney U-test, *P* = 0.006). No such difference was observed between GAL and GIL in the AO and EM basins for *H̄*<sub>exp</sub> (*P* = 0.055 and *P* = 0.12, respectively).

## Hardy–Weinberg equilibrium and linkage disequilibrium

After corrections for multiple tests, no significant departure from HWE was detected for any individual samples belonging to the AO and WM basins, whatever data were considered (not reported). HWE was also respected after pooling the WM samples into a single metapopulation (all markers: *f* = 0.01, NS; GAL: *f* = -0.01, NS; GIL: *f* = 0.01, NS; Table 3). For the EM populations, there was a significant heterozygote excess over the 15 loci in the SYRA (Syracuse) population (*f* = -0.22; *P* < 0.001), while HWE was respected in other individual samples. This induced a significant heterozygote excess in the EM metapopulation (*f* = -0.07; *P* < 0.001). It was mostly due to GAL markers in the SYRA and SEL (Selinunte) populations (*f* = -0.24, *P* < 0.001; *f* = -0.08, *P* < 0.05, respectively). No significant departure from HWE was detected in the EM metapopulation for the multilocus GIL marker set (*f* = -0.06; NS; Table 3).

Because larger sample sizes prevent detection of spurious genotypic associations among loci, analysis of LD was primarily performed at the metapopulation level for Mediterranean samples (i.e., AO consists of a single sample).

After correction for multiple tests, there were only three significant cases of LD in the WM metapopulation. In the EM metapopulation, 12 significant cases of LD were found, that involved only four loci: two GIL and two GAL markers (DLA0068, DLA0086, DLA0070 and DLA0089).

## Relationships between ˆθ and gene diversity

Relationship between ˆθ and *H*<sub>exp</sub> within WM and EM revealed no significant negative correlations (WM<sub>[20 loci]</sub>: R<sup>2</sup> = 0.062, NS; EM<sub>[15 loci]</sub>: R<sup>2</sup> = 0.129, NS; Fig. S1). There were no significant correlations by pairwise metapopulation comparison, AO–WM, AO–EM and WM–EM, nor comparison among the three metapopulations (AO–WM<sub>[20 loci]</sub>: R<sup>2</sup> = 0.158, NS; AO–EM<sub>[15 loci]</sub>: R<sup>2</sup> = 0.052, NS; WM–EM<sub>[15 loci]</sub>: R<sup>2</sup> = 0.043, NS; three metapopulations<sub>[15 loci]</sub>: R<sup>2</sup> = 0.028, NS). When present, homoplasy did not significantly affect data and these results legitimate further interpretation in terms of differential gene flow and/or selection.

## Genetic differentiation within the WM and EM metapopulations

Within the WM basin, no significant differentiation was detected ( ˆθ = 0.004, NS). There was no differentiation for the GIL or the GAL data separately ( ˆθ = 0.001, NS; ˆθ = 0.004, NS, respectively). Within the EM basin, there
was significant differentiation only for the GAL data ($\hat{\theta} = 0.049$, $P < 0.01$), but not for the full or GIL data ($\hat{\theta} = 0.021$, NS; $\hat{\theta} = -0.006$, NS, respectively).

## Genetic differentiation among basins

Reynolds' distance trees based on the 15 locus data clearly partitioned individual samples into three basins (Fig. 2). Levels of differentiation were generally significant and dependent on the pairwise basin comparison and markers used (GAL or GIL; Fig. 2). The differentiation between AO and WM was greater than between WM and EM, whatever data considered (Fig. 2).

The individual values of $\hat{\theta}$ for each locus in the AO–WM and the WM–EM comparisons are reported in Figure 3A,B, respectively. Between AO and WM (20 loci used), five loci had $\hat{\theta} > 10\%$, including locus *mGH* where $\hat{\theta} \approx 25\%$. The multilocus differentiation estimator $\hat{\theta}$ was 0.056 ($P < 0.001$) with significant differentiation at nine markers (GAL: six markers; GIL: three markers, Fig. 3A). The GAL data had higher values of $\hat{\theta}$ than the GIL data ($\hat{\theta} = 0.073$ vs. $\hat{\theta} = 0.034$; both $P < 0.001$; Fig. 2). Using 15 loci, levels of genetic differentiation were significant between WM and EM ($\hat{\theta} = 0.026$; $P < 0.001$), with significant individual estimates of $\hat{\theta}$ at nine of the 15 loci (GAL: four markers; GIL: five markers; Fig. 3B). Contrary to the AO–WM comparison, differentiation at the GIL markers was greater than at the GAL markers in the WM

–EM comparison (respectively, $\hat{\theta} = 0.036$ and $\hat{\theta} = 0.020$, both $P < 0.001$; Fig. 2).

## Patterns of genetic differentiation

There was wide variability in spatial patterns of allelic frequency changes (Fig. 4), including cases where allele frequencies were (1) almost identical in all metapopulations, reflecting either high gene flow, shared ancestral polymorphisms, or general selective sweeps at these loci (e.g., *DLA0089*, *µGH*, resulting in nonsignificant $\hat{\theta}$ values in Fig. 3B), (2) identical in two adjacent metapopulations (AO–WM or WM–EM) then marking a sharp break at one of the transition zones (e.g., *DLA0051*, *DLA0097*), or (3) in an apparent east–west gradient of allele frequencies (e.g., *DLA0068*). Those latter two classes of loci naturally coincided with loci achieving higher levels of differentiation (Fig. 3). Locus *DLA0060* appeared as a particular case with almost identical allele frequencies for the metapopulations located at each side of the distribution but drastic change in frequency for the most frequent alleles in the central metapopulation.

## Population assignment and migration

Analyses regarding assignment were initially performed for the 15 markers in all metapopulations, then only using GAL ($n = 9$; *DLA0044*, *DLA0051*, *DLA0060*,

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Reynold's distance trees for each marker set used in this study (all markers [15 loci]: GAL [9 loci] and GIL [6 loci]). Figures located at the bottom of each tree summarized patterns of genetic differentiation among samples within metapopulation (no estimate for Atlantic Ocean [AO] as one unique sample was considered in this metapopulation; Table 1), and among metapopulations as estimated by $\hat{\theta}$ for each marker set. In the case of the AO–Western Mediterranean, $\hat{\theta}$ values are also reported in brackets for the 20 loci marker set, then for the set of 11 and 9 GAL and GIL marker sets, respectively. Samples are abbreviated as in Table 1 and colored symbols as in Fig. 1. **$P < 0.01$, ***$P < 0.001$.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Estimates of $\hat{\theta}$ observed at each individual locus for (a) the Atlantic Ocean–Western Mediterranean (WM) comparison (20 loci), and (b) the WM–Eastern Mediterranean (EM) comparison (15 loci). Loci are arranged according to increased values Of $\hat{\theta}$ in each comparison. Horizontal lines indicate mean values of $\hat{\theta}$ observed in each comparison (shaded area = 95% CI). Vertical lines indicate which individual loci were found significant (on the right) compared to those that were not (on the left) for each comparison. Open and black boxes indicate GIL and GAL, respectively. Asterisks indicate loci absent in the WM–EM comparison.</figcaption></figure>

DLA0070, DLA0078, DLA0089, µGH, µIGF1, µSL) or GIL ($n = 6$; DLA0041, DLA0061, DLA0068, DLA0075, DLA0086, DLA0097; Fig. 5). Unfortunately, some loci with a dramatically high $\hat{\theta}$ between AO and WM could not be considered because they failed to amplify in EM (e.g., mGH; Fig. 3A). This lowered the probability of assignment for AO and WM individuals. Nevertheless, better mean average assignment to their population of origin was achieved for AO and EM individuals (mean $\hat{q}$-value = 0.548 ± 0.156 and 0.868 ± 0.151, respectively) than for the WM individuals ($\hat{q}$ = 0.326 ± 0.206; WM individuals were in average better assigned to the EM cluster: $\hat{q}$ = 0.487 ± 0.291 when using all loci). This suggested higher introgression in the WM compared with the other basins (Fig. 5), and correlated well with lower all-marker genetic differentiation estimated at the STS compared to the AOF (Fig. 2). When analyzing the GAL, only EM individuals were better assigned on average to their actual population of origin ($\hat{q}$ = 0.596 ± 0.223), while correct assignment was far lower for other meta-populations (AO: $\hat{q}$ = 0.417 ± 0.117; WM: $\hat{q}$ = 0.370 ± 0.134). When the GIL was used and whatever the

metapopulation considered, probabilities of correct assignment did not differ significantly from $\hat{q}$ = 0.33 (details not reported). It has to be emphasized that results for GIL in this analysis did not indicate lower intrinsic relevance of GIL to assign individual, nor increased gene flow for this class of loci, but were primarily due to the low number of GIL in this particular analysis ($n = 6$) with only one locus demonstrating genetic differentiation among each pair of basins (DLA0068; Fig. 3). At the scale of the three metapopulations, it resulted that those six loci collectively had significant but low average $\hat{\theta}$ = 0.031 ($P$ < 0.001), resulting in poor ability to successful individual assignment. Using simulated data sets in a five deme system, Latch et al. (2006) reported an average value of $\theta$ = 0.030 as a lower bound to consistently assign individuals to a subpopulation (i.e., with a probability distinct from equiprobability). Empirical $\hat{\theta}$ observed in sea bass is very close from this lower bound, easily explaining poor assignment capacities with this data set.

Individual WM samples located close to the STS (e.g., GOU, MAR; Fig. 1) present a compound allele pattern slightly closer to that of individuals from the EM samples
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Variation of allele frequencies observed at each individual locus among sea bass metapopulations. For each marker, colored lines represent variation for a single allele. Because of size of the panels associated to loci, all alleles might not be visible at each locus. Metapopulations are abbreviated as in Table 1.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Summary of results obtained by MIGRATE-n regarding the estimation of number of migrants among sea bass metapopulations for each set of markers considered in this study (15 loci data set; left panel), and STRUCTURE regarding probability of membership of each individual to one of the K = 3 sea bass metapopulations considered in this study. Thickness of arrows is broadly proportional to estimated number of migrants. Samples and metapopulations are abbreviated as in Table 1.</figcaption></figure>

than to WM samples located further west (e.g., MUR) for all assignments, either performed over the full data set or with the random subsets of 31 individuals. Despite this, WM individuals from western and eastern WM samples were not assigned significantly differently from EM samples over five replicates, whatever the number and class of markers considered (details not reported).

Maximum-likelihood estimates of migrations using MIGRATE-n revealed asymmetric exchange from AO to WM with a approximately threefold to fivefold difference in number of successful migrants when using GAL and all markers (GAL: $M_{AO→WM}$ = 43.03 vs. $M_{WM→AO}$ = 12.32; all markers: $M_{AO→WM}$ = 34.47 vs. $M_{WM→AO}$ = 6.22; Fig. 5). This asymmetry of migration between AO and WM was inferior to two for the GIL data set (Fig. 5). Between WM and EM, the asymmetry was also estimated with a higher number of alleles from EM migrating toward WM. The extent of the asymmetrical exchange was however lower when considering all markers ($M_{EM→WM}$ = 15.64 vs. $M_{WM→EM}$ = 9.62), or the GIL and GAL data sets separately (Fig. 5). Results of MIGRATE-n therefore indicate greater introgression in WM than in other metapopulations. That is, WM received more alleles by migration than they provided to the other two metapopulations, suggesting it may represent a hybrid swarm separated by two tension zones at the AOF and STS where

migration of genes are asymmetrically filtered. It was however expected that lower estimates of genetic differentiation should translate in higher estimated number of migrants at the STS. This point has to be investigated further.

## Neutrality tests

Neutrality tests were initially performed on individual samples within the WM and EM metapopulations, respectively. DETSEL was not used at this scale because of too many, not fully independent pairwise comparisons. In WM, no congruence between the B&N or the F&G tests was detected among potentially selected loci when using *P* = 0.95 (or Log₁₀(BF) = 1.3). No locus was significant for a more stringent probability or BF (not reported). Indeed, B&N indicated that three loci were outliers and hence could potentially be under selection (μSL, *DLA0066* and *DLA0070*) in WM, while F&G detected eight different loci (*DLA0044*, *DLA0068*, *DLA0081*, *DLA0089*, *DLA0096*, μIGF-1, μGH, and mGH). In EM, F&G did not detect potentially selected loci whatever the Log₁₀(BF) considered, while B&N detected three loci at *P* = 0.95 (μSL, μGH, and μIGF-1; none at *P* = 0.99). It must be noted that those loci demonstrated the highest levels of differentiation among EM samples (details not shown). As results of neutrality tests failed to
converge, the reality of nonneutrality within metapopulations is open to question.

Results were slightly more conclusive when performing neutrality tests among metapopulations. In the AO–WM comparison, the three tests did not provide similar results, but all indicated locus *DLA0068* (GIL) as a possible outlier (Table 4). For the WM–EM comparison, B&N and DetSel detected the *DLA0060* (GAL), *DLA0068*, and *DLA0097* (GIL) loci as outliers for a minimal probability level of *P* = 0.95. F&G also detected three outlier loci for Log₁₀ (BF) > 1.3: *DLA0068*, *μIGF-1* and *μSL* (Table 4). Locus *DLA0068* was detected by all three tests and was the locus with the sharpest pattern of graded variation across metapopulations (Fig. 4).

## Discussion

Genetic differentiation between Atlantic and Mediterranean populations of marine organisms is usually interpreted as the signature of secondary contact between two formerly allopatric lineages, occurring since the last glacial maximum (Borsa et al. 1997; Patarnello et al. 2007; but see, e.g., Grant 2005; Limborg et al. 2012). Accordingly, support for two distinct mtDNA clades, or nuclear differentiation, between Atlantic and Mediterranean populations at the AOF is numerous in marine organisms (Patarnello et al. 2007). When differentiation at the

<table><caption>**Table 4.** Summary of loci found as being potentially selected at the metapopulation level in each outlier detection selective test used in this study (B&N: Beaumont and Nichols 1996; DetSel: Vitalis et al. 2001; F&G: Foll and Gaggiotti 2008).</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Locus</th><th colspan="3">Neutrality Test</th></tr><tr><th>B&N</th><th>DetSel</th><th>F&G</th></tr></thead><tbody><tr><td rowspan="7">A: AO vs. WM</td><td>DLA0051</td><td>0.993</td><td>&gt;0.99</td><td>—</td></tr><tr><td>DLA0060</td><td>—</td><td>&gt;0.95</td><td>—</td></tr><tr><td>DLA0061</td><td>—</td><td>&gt;0.95</td><td>0.992</td></tr><tr><td><strong>DLA0068</strong></td><td><strong>0.973</strong></td><td><strong>&gt;0.99</strong></td><td><strong>1.0</strong></td></tr><tr><td>DLA0070</td><td>0.995</td><td>&gt;0.95</td><td>—</td></tr><tr><td>DLA0075</td><td>0.999</td><td>—</td><td>0.976</td></tr><tr><td>DLA0096</td><td>0.958</td><td>&gt;0.99</td><td>—</td></tr><tr><td>DLA0097</td><td>0.999</td><td>—</td><td>0.952</td></tr><tr><td rowspan="7">B: WM vs. EM</td><td>mGH</td><td>0.978</td><td>&gt;0.99</td><td>—</td></tr><tr><td>μIGF1</td><td>0.997</td><td>—</td><td>1.0</td></tr><tr><td>μSL</td><td>0.967</td><td>—</td><td>0.983</td></tr><tr><td>DLA0060</td><td>0.993</td><td>&gt;0.95</td><td>—</td></tr><tr><td><strong>DLA0068</strong></td><td><strong>0.964</strong></td><td><strong>&gt;0.99</strong></td><td><strong>0.982</strong></td></tr><tr><td>DLA0097</td><td>0.996</td><td>&gt;0.99</td><td>—</td></tr><tr><td>μIGF1</td><td>—</td><td>—</td><td>1.0</td></tr><tr><td></td><td>μSL</td><td>—</td><td>—</td><td>0.997</td></tr></tbody></table>

AO, Atlantic Ocean; EM, Eastern Mediterranean; WM, Western Mediterranean.

Locus *DLA0068* found as an outlier in each test is indicated in bold.

AOF is not demonstrated (e.g., Bargelloni et al. 2003), differentiation is generally detected in another area of the Mediterranean, including the STS and some other locations (e.g., south of Peloponnesus; Patarnello et al. 2007). The originality of sea bass resides in their partition into three metapopulations, separated by two areas of reduced genetic connectivity: the AOF and the STS (Naciri et al. 1999; Bahri-Sfar et al. 2000; Lemaire et al. 2005). To our knowledge, reports of two locations of reduced connectivity within the Mediterranean are pretty scarce, apart from the Black Sea (but see Rolland et al. 2007; Zulliger et al. 2009; Limborg et al. 2012). In sea bass, the hypothesis of a purely physical barrier shaping variation at the AO–WM transition has been seriously reassessed by Lemaire et al. (2005), who demonstrated a tension zone at this location. By using nuclear markers (GAL and GIL), we expanded on the study by Lemaire et al. (2005) to dissect the sea bass’ patterns of genetic variation further.

## Patterns of genetic differentiation and support for a tension zone at the Siculo-Tunisian strait

This study did not change the perception of sea bass genetic structure within each basin at nuclear loci (i.e., no differentiation inside the WM basin but some differentiation within the EM basin; Garcia de León et al. 1997; Naciri et al. 1999; Bahri-Sfar et al. 2000; Castilho and Ciftci 2005), nor the existence of three metapopulations. Nevertheless, estimates of differentiation among successive pairs of basins were larger than previously estimated with fewer microsatellite loci. The differentiation between the AO and WM basin reported in this study ($\hat{\theta}_{[20 \text{ loci}]}$ = 5.6%) was more than twice that provided by Naciri et al. (1999) over six loci ($\hat{\theta}$ ≈ 2.1%). Estimates of differentiation at the WM–EM transition zone were also higher ($\hat{\theta}_{[15 \text{ loci}]}$ ≈ 2.6%) than those reported by Bahri-Sfar et al. (2000) ($\hat{\theta}$ ≈ 1.4%; six loci). Results hence confirmed a stronger nuclear genetic barrier at the AOF than at the STS. This pattern was also demonstrated for mtDNA in sea bass (Lemaire et al. 2005; R. Rondon, E. Desmarais, F. Bonhomme and B. Guinand, unpubl. results), but also in sprat (Limborg et al. 2012). Using far less nuclear markers, Rolland et al. (2007) suggested on the contrary that differentiation should be stronger at the STS compared with the AOF for the common sole (*Solea solea*).

GAL and GIL markers may both be involved in overall differentiation and/or differentiation among adjacent metapopulations, but the extent of differentiation depended on the particular comparison. Indeed, GAL that are located closer from or even within genes consistently demonstrated in average reduced genetic connectivity compared GIL at the AO–WM transition (20 loci;
$\hat{0}_{\text{GAL}} = 7.3\%; \hat{0}_{\text{GIL}} = 3.4\%$). However, a reverse relationship was found at the WM-EM transition (15 loci; $\hat{0}_{\text{GAL}} = 2.0\%; \hat{0}_{\text{GIL}} = 3.6\%$). This second observation hence partly conflicts with those reported in other studies in which gene-associated markers had, on average, better discriminatory power than independent markers (e.g., Vasemägi et al. 2005; Yatabe et al. 2007; Shikano et al. 2010; Vilas et al. 2010).

Besides differentiation, it is also interesting to note that significant LD and heterozygote excess was concentrated in a few samples. Heterozygote excess was especially marked for GAL at the WM-EM transition on both sides of Sicily: for the SEL and GOU samples, and at the Messina canal for the SYRA sample. LD only occurred in the SEL and SYRA samples from the EM basin. This observation is unlikely to be due to low sizes of SYRA and SEL samples, as others with similar sizes did not demonstrate any heterozygote excess and very few cases of LD. The GOU, SEL, and SYRA samples are all located at the WM-EM transition, suggesting that heterozygote excess and/or LD may be the signature of hybridization among parental pools at the STS, as already demonstrated by Lemaire et al. (2005) at the AOF. There, significant LD was only found in the sample closest to the AOF. In this study, we report a nearly identical situation at the STS (but note that the MAR sample also located close to the front did not show significant LD and/or heterozygote excess, in line with previous results by Naciri et al. (1999) and Bahri-Sfar et al. (2000) who also used this sample). Results are in agreement with the definition of a tension zone rather than with strict isolation by distance or parapatric differentiation. By theoretical modeling, Irwin (2012) has shown that genetically differentiated haploid clades arising from local adaptation along an environmental gradient can arise in a parapatric fashion under certain conditions and persist in a metapopulation even in the presence of moderate gene flow. Nevertheless, extension to multilocus diploid cases in high gene flow species has not been done, so it seems reasonable to stick to the less costly hypothesis of historical differentiation followed by secondary contact. If there was a single reason for this more parsimonious explanation, this would be at least because paleobiogeography indicates that favorable ecological zones have indeed been fragmented several time during the glacial maximums, notwithstanding other theoretical considerations about all the conditions required for the progressive building up of multilocus differentiation in presence of gene flow in highly dispersing species. We hence can reasonably be confident that the encounter of differentiated genetic backgrounds in the sea bass is a likely possibility as demonstrated for other marine organisms with large dispersal capabilities (e.g., bivalve mollusks, Nikula et al. 2008).

## Spatial patterns support tension zones, but rarely selection at target loci

In addition to LD and/or heterozygote excess, spatial patterns of allele frequency changes are in accordance with the prediction of a tension zone at both the AO-WM and WM-EM transitions. The theory predicts that allele frequency changes are more or less pronounced depending on the marker investigated, with steeper transitions for loci closer to “isolation genes” (Barton and Hewitt 1985). In this study, some markers present no change in allele frequency, reflecting either one unconstrained high neutral gene flow, a selective sweep of one uniformly advantageous allele or the inheritance of an ancestral state where allele frequencies were already homogeneous. On the other hand, some loci show relatively sharp changes at AO-WM transition as already shown for mtDNA (Lemaire et al. 2005), or at the WM-EM transition, which has been less extensively investigated (Bahri-Sfar et al. 2000). A vicariant event that produced at least two separate Atlantic and Mediterranean metapopulations of sea bass during one of the quaternary climatic oscillations will have led to the emergence of isolation genes. These endogenous barriers influenced many other loci by indirect (hitchhiking) selection at genomic scales, depending on their strength and the amount of recombination (Barton and Hewitt 1985; Piálek and Barton 1997; Slatkin and Wiehe 1998). According to general hybrid zone models, we postulate that barriers have been trapped at various locations which correspond to areas of reduced density and/or environmental change, such as the AOF and STS. A rationale of this phenomenon is provided in Bierne et al. (2011), including a review of marine examples (see also Luttikhuizen et al. 2012).

Locus *DLA0060* (GAL; located in intron one of the bestrophin three gene) presents a distinct pattern of spatial allele frequency compared to other loci, with the two most frequent alleles at higher frequencies in WM than in the marginal AO and EM metapopulations. In the present case, we cannot distinguish whether locus *DLA0060* was linked to an “endogenous” locus whose variation in allele frequency reflects genomic incompatibilities, or to an “exogenous” locus whose variation reflects local adaptation, both fitting with the coupling hypothesis of Bierne et al. (2011).

In hitchhiking models of genetic variation, direct selection on markers is not needed (Piálek and Barton 1997; Bierne et al. 2003). In sea bass, neutrality tests were unable to identify loci with clear footprint of selection, except *DLA0068* (GIL). This locus presents a steeper gradient of allele frequencies than other loci, with a nearly fixed allele in the AO sample and an alternative one that predominates in the EM sample. This marker was not clearly associated with a gene. It is situated at approximately
4 kb from the 3' end of a Retinal G Protein Coupled Receptor gene (RGR) on one side and approximately 6 kb from the 5' end of the CDS of a solute carrier gene (SLC18A3 also known as VACHT, e.g., Weihe et al. 1996; Karila et al. 1998) on the other. We cannot therefore reach firm conclusions about its functional importance.

## A hybrid swarm in the Western Mediterranean?

The WM contains strongly introgressed individuals, but also receives migrants from other basins. Taken together, this indicates that the WM metapopulation may have a dual origin from both AO and EM nuclear genomes, which have introgressed into each other and led to the emergence of a new reproductive unit, known as a hybrid swarm (Arnold 1993; Jiggins and Mallet 2000). The absence of LD among loci in each WM sample, except the GOU sample located close to the STS, also supports this hypothesis of a sea bass hybrid swarm. Interestingly, Schunter et al. (2011) recently reported a similar result in the comber *Serranus cabrilla* within the WM basin. In their case, combers from locations located between the two recognized fronts for this species (AOF and the Ibiza Channel) were significantly more introgressed than individuals from populations on the external sides of each front, although the presence of tension zones was not investigated. Our results on sea bass are also comparable to experimental studies by Edmands et al. (2005) on the copepod *Tigriopus californicus*. After approximately 15 generations of mixing, Edmands et al. (2005) found more hybrids than expected, together with heterozygote excess. These experimental results conform to the situation observed in wild sea bass near tension zones, as heterozygote excess and LD were recorded at the STS and AOF (this study and Lemaire et al. 2005, respectively). Overall, the results suggest that the WM metapopulation of sea bass could result from hybrid swarming, that is separated from parental forms by two distinct tension zones.

Hybrid swarms are poorly documented in marine organisms (McMillan et al. 1999; Roques et al. 2001; Roberts et al. 2010) by comparison with terrestrial or freshwater taxa (e.g., Arnold 1997; Nolte et al. 2009; Stemshorn et al. 2011). They may be restricted to narrow marine zones (estuaries; see Roberts et al. 2010) or extend over larger areas (Ayres et al. 2004; Strelkov et al. 2007; Nikula et al. 2008). Strelkov et al. (2007) analyzed patterns of genetic variation in a hybrid swarm of *M. balthica* that extends over hundreds of kilometers in the Norwegian, Barents, and Pechora seas. The samples within this swarm also showed a component of regional differentiation that was not observed in sea bass: within WM, no differentiation

was detected and there was no significant indication of differential introgression with AO and EM individuals for samples located to the west or the east.

The random coupling of exo- and endogenous incompatibilities hypothesized by Bierne et al. (2011) can elicit emergence of distinct tension zones between successive fragmented patches of habitat. Using a two-locus stepping-stone model of gene flow with overdominance, Goodisman and Crozier (2001) also demonstrated that two spatially distinct tension zones can emerge between formerly isolated populations, giving rise to an additional central population. Goodisman and Crozier (2001) also demonstrated that LD accumulated at the two tension zones separating the three successive populations, while it disappeared in the central population, as observed in sea bass in the WM basin.

Although sea bass population structure and pattern of population differentiation conform to the expectations of some models, it is essential to acquire further data to investigate the hypothesis of a hybrid swarm in the WM. Such caution particularly derives from statistical considerations. This includes the discrepancy between estimated levels of differentiation and estimated number of migrants at the STS. It was expected that lower genetic differentiation should translate in higher number of estimated migrants, but this trend was not observed possibly casting doubt on this particular result. Furthermore, numerous methods for analysis of hybridization processes have emerged over recent years. They range from identification of F₁ to Fₓ hybrids and backcrosses (e.g., Anderson and Thompson 2002), to evaluating the fitness effects of alleles transferred across different genetic backgrounds (i.e., introgression; Gompert and Buerkle 2009), to genetic diffusion models (Stemshorn et al. 2011). Methods aim to analyze introgression (admixture) patterns together with the dynamics and history of introgression (Nolte and Tautz 2010; Lamaze et al. 2012). To derive consistent outputs, however, they require clear recognition of parental genotypes (Nolte et al. 2006, 2009; Gagnaire et al. 2011; Luttikhuizen et al. 2012). This is typically not observed in sea bass, where parental genotypes (AO, but also EM metapopulations) are both significantly introgressed, hence limiting the capacity for reliable inference. High throughput genome-wide genotyping will probably provide the statistical power to disentangle recent from old and shared from unshared polymorphism across transition zones in sea bass.

## Acknowledgments

Authors wish to thank B. Chatain, M. Vandeputte (Ifremer, Palavas, France), D.M. Rand (Brown University, Providence, USA), and F. Volckaert (Katholieke Universi-
teit Leuven, Belgium) for constructive advices. M. Hassan and M. Arculeo kindly provided, respectively, the Syrian and Sicilian samples used in this study. We are also grateful to R. Reinhardt, H. Kuhl, M. Tine at MPI for Molecular Genetics, Berlin-Dahlem, for supplying unpublished data from their “European sea bass sequencing project,” and to D. McKenzie for English editing. Genetic data were produced at facilities of the HCMR and of the LabEx CeMEB (Centre Méditerranéen Environnement Biodiversité). This is publication ISEM2012-154 from the Institut des Sciences de l'Évolution de Montpellier (CNRS UMR 5554).

## Conflict of Interest

None declared.

## References

Anderson, E. C., and E. A. Thompson. 2002. A model-based method for identifying species hybrids using multilocus genetic data. Genetics 160:1217-1229.

Antao, T., A. Lopes, R. J. Lopes, A. Beja-Pereira, and G. Luikart. 2008. LOSITAN: A workbench to detect molecular adaptation based on a $F_{ST}$-outlier method. BMC Bioinformatics 9:323. doi:10.1186/1471-2105-9-323.

Arnold, J. 1993. Cytonuclear disequilibria in hybrid zones. Annu. Rev. Ecol. Syst. 24:521-554.

Arnold, J. 1997. Natural hybridization and evolution. Oxford Univ. Press, Oxford, U.K. 232 pp.

Ayres, D. R., D. L. Smith, K. Zaremba, S. Klorh, and D. R. Strong. 2004. Spread of exotic cordgrasses and hybrids (*Spartina* sp.) in the tidal marshes of San Francisco Bay, California, USA. Biol. Invasions 6:221-231.

Bahri-Sfar, L., C. Lemaire, O. K. Ben Hassine, and F. Bonhomme. 2000. Fragmentation of seabass populations in the western and eastern Mediterranean as revealed by microsatellite polymorphism. Proc. R. Soc. B Biol. Sci. B267:929-935.

Barber, P. H., S. R. Palumbi, M. V. Erdmann, and M. K. Moosa. 2000. A Marine Wallace's Line? Nature 406:692-693.

Bargelloni, L., J. A. Alarcon, M. C. Alvarez, E. Penzo, A. Magoulas, C. Reis, et al. 2003. Discord in the family Sparidae (Teleostei): divergent phylogeographical patterns across the Atlantic-Mediterranean divide. J. Evol. Biol. 16:1149-1158.

Barton, N. H., and G. M. Hewitt. 1985. Analysis of hybrid zones. Annu. Rev. Ecol. Syst. 16:113-148.

Beaumont, M. A., and R. A. Nichols. 1996. Evaluating loci for use in the genetic analysis of population structure. Proc. R. Soc. B Biol. Sci. B263:1619-1626.

Beerli, P. 2008. MIGRATE documentation (version 3.0). Technical report. Available at http://popgen.sc.fsu.edu. Accessed May 23, 2010.

Beerli, P., and J. Felsenstein. 1999. Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. Genetics 152:763-773.

Beerli, P., and J. Felsenstein. 2001. Maximum likelihood estimation of a migration matrix and effective population sizes in *n* subpopulations by using a coalescent approach. Proc. Natl Acad. Sci. USA 98:4563-4568.

Benharrat, K., N. Pasteur, Y. Siau, and A. Bouain. 1983. Polymorphisme biochimique de loups (*Dicentrarchus labrax*) originaires de quatre populations naturelles et d'un élevage. Pp. 1-17. Recherches Biologiques et Aquaculture. Vol. I. Journées d'études de Montpellier. CNEXO-BNDO, Brest, France (in French).

Benjamini, Y., and D. Yukitieli. 2001. The control of false discovery rate under dependency. Ann. Stat. 29:1165-1188.

Benzie, J. A. H. 1999. Genetic structure of coral reef organisms: ghosts of dispersal past? Am. Zool. 39:131-145.

Bierne, N. 2010. The distinctive footprints of local hitchhiking in a varied environment and global hitchhiking in a subdivided population. Evolution 64:3254-3272.

Bierne, N., C. Daguin, F. Bonhomme, P. David, and P. Borsa. 2003. Direct selection on allozymes is not required to explain heterogeneity among marker loci across a *Mytilus* hybrid zone. Mol. Ecol. 12:2505-2510.

Bierne, N., J. Welch, E. Loire, F. Bonhomme, and P. David. 2011. The coupling hypothesis: why genome scans may fail to map local adaptation genes. Mol. Ecol. 20:2044-2072.

Borsa, P., M. Naciri, L. Bahri, L. Chikhi, F. J. Garcia de León, G. Kotoulas, et al. 1997. Zoogéographie infraspécifique de la mer Méditerranée: analyse des données génétiques populationnelles sur seize espèces Atlanto-Méditerranéennes (poissons et invertébrés). Vie Milieu 47:295-305.

Boutet, I., N. Quéré, F. lecomte, J. F. Agnèse, and B. Guinand. 2008. Putative transcription factor binding sites and polymorphisms in the proximal promoter of the *PRL-A* gene in percomorphs and European sea bass (*Dicentrarchus labrax*). Mar. Ecol. 29:354-364.

Burton, R. S., C. K. Ellison, and J. S. Harrison. 2006. The sorry state of $F_2$ hybrids: consequences of rapid mitochondrial DNA evolution in allopatric populations. Am. Nat. 168:S14-S24.

Castilho, R., and Y. Ciftci. 2005. Genetic differentiation between close eastern Mediterranean *Dicentrarchus labrax* (L.) populations. J. Fish Biol. 67:1746-1752.

Chatain, B., and H. Chavanne. 2009. La génétique du bar (*Dicentrarchus labrax*, L.). Cah. Agric. 18:249-255.

Chistiakov, D. A., B. Hellemans, C. S. Tsigenopoulos, A. S. Law, N. Bartley, D. Bertotto, et al. 2004. Development and
linkage relationships for new microsatellite markers of the
sea bass (*Dicentrarchus labrax* L.). Anim. Genet. 35:53–57.

Coscia, I., and S. Mariani. 2011. Phylogeography and population structure of European sea bass in the north-east Atlantic. Biol. J. Linn. Soc. 104:364–377.

Edmands, S., H. V. Feaman, J. S. Harrison, and C. C. Timmerman. 2005. Genetic consequences of many generations of hybridization between divergent copepod populations. J. Hered. 96:114–123.

Evanno, G., S. Regnaut, and J. Goudet. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14:2611–2620.

Felsenstein, J. 2005. PHYLIP (phylogeny inference package) version 3.6. Department of Genome Sciences, Univ. of Washington, Seattle, WA.

Foll, M., and O. Gaggiotti. 2008. A genome scan method to identify selected loci appropriate for both dominant and co-dominant markers: a Bayesian perspective. Genetics 180:977–993.

Fritsch, M., Y. Morizur, E. Lambert, F. Bonhomme, and B. Guinand. 2007. Assessment of sea bass (*Dicentrarchus labrax*, L.) stock delimitation in the Bay of Biscay and the English Channel based on mark-recapture and genetic data. Fish. Res. 83:123–132.

Gagnaire, P. A., Y. Minegishi, S. Zenboudji, P. Valade, J. Aoyama, and P. Berrebi. 2011. Within-population structure highlighted by differential introgression across semipermeable barriers to gene flow in *Anguilla marmorata*. Evolution 65:3513–3525.

Galarza, J. A., J. Carreras-Carbonell, E. Macpherson, M. Pascual, S. Roques, G. F. Turner, et al. 2009. The influence of oceanographic fronts and early-life history traits on connectivity among littoral fish species. Proc. Natl Acad. Sci. USA 106:1473–1478.

Galindo, H. M., D. B. Olson, and S. R. Palumbi. 2006. Seascape genetics: a coupled oceanographic-genetic model predicts population structure of Caribbean corals. Curr. Biol. 16:1622–1626.

García de León, F. J., L. Chikhi, and F. Bonhomme. 1997. Microsatellite polymorphism and population subdivision in natural populations of European sea bass *Dicentrarchus labrax* (Linnaeus, 1758). Mol. Ecol. 6:51–62.

García, F., G. Zalba, G. Páez, I. Encío, and C. de Miguel. 1998. Molecular cloning and characterization of the human p44 mitogen-activated protein kinase gene. Genomics 50:69–78.

Gardner, J. P. A. 1997. Hybridization in the sea. Adv. Mar. Biol. 31:1–78.

Gemayel, R., M. D. Vinces, M. Legendre, and K. J. Verstrepen. 2010. Variable tandem repeats accelerate evolution of coding and regulatory sequences. Annu. Rev. Genet. 44:445–477.

Gerlach, G., J. Atema, M. Kingsford, K. P. Black, and V. Miller-Sims. 2007. Smelling home can prevent dispersal of reef fish larvae. Proc. Natl Acad. Sci. USA 104:858–863.

Gompert, Z., and C. A. Buerkle. 2009. A powerful regression-based method for admixture mapping of isolation across the genome of hybrids. Mol. Ecol. 18:1207–1224.

Goodisman, M. A., and R. H. Crozier. 2001. Clines maintained by overdominant selection in hybrid zones. Hereditas 134:161–169.

Grant, W. S. 2005. A second look at mitochondrial DNA variability in European anchovy (*Engraulis encrasicolus*): assessing models of population structure and the Black Sea isolation hypothesis. *Genetica* 125:293–309.

Guinand, B., E. Dujardin, V. Dufour, and C. S. Tsigenopoulos. 2008. Characterisation of genetic structure of *Dicentrarchus labrax* larvae in two nurseries of the Gulf of Lions (NW Mediterranean). Aquat. Living Resour. 21:81–87.

Harrison, R. G. (ed.) 1993. Hybrid zones and the evolutionary process. Oxford Univ. Press, New York. 364 pp.

Hauser, L., and G. R. Carvalho. 2008. Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts. Fish. Fish. 9:333–362.

Hellberg, M. E. 2009. Gene flow and isolation among populations of marine animals. Annu. Rev. Ecol. Evol. Syst. 40:291–310.

Hellberg, M. E., R. S. Burton, J. E. Neigel, and S. R. Palumbi. 2002. Genetic assessment of marine population connectivity. Bull. Mar. Sci. 70(Suppl.):273–290.

Hereford, J. 2009. A quantitative survey of local adaptation and fitness trade-offs. Am. Nat. 173:579–588.

Hobbs, J. P. A., A. J. Frisch, G. R. Allen, and L. Van Herwerden. 2009. Marine hybrid hotspot at Indo-Pacific biogeographic border. Biol. Lett. 5:258–261.

Hwang, A. S., S. L. Northrup, J. K. Alexander, K. T. Vo, and S. Edmands. 2011. Long-term experimental hybrid swarms between moderately incompatible *Tigrionus californicus* populations: hybrid inferiority in early generations yields to hybrid superiority in later generations. Conserv. Genet. 12:895–902.

Irwin, D. E. 2012. Local adaptation to smooth ecological gradients causes phylogeographical breaks and phenotypic clustering. Am. Nat. 180:35–49.

Ishida, N., S. Yoshioka, Y. Chiba, M. Takeuchi, and M. Kawakita. 1999. Molecular cloning and functional expression of the human Golgi UDP-N-Acetylglucosamine transporter. J. Biochem. 126:68–77.

Janssens, V., and J. Goris. 2001. Protein phosphatase 2A: a highly regulated family of serine/threonine phosphatases implicated in cell growth and signalling. Biochem. J. 353:417–439.

Jiggins, C. D., and J. Mallet. 2000. Bimodal hybrid zones and speciation. Trends Ecol. Evol. 15:250–255.

Karila, P., F. Shahbazi, J. Jensen, and S. Holmgren. 1998. Projections and actions of tachykininergic, cholinergic, and serotonergic neurones in the intestine of the Atlantic cod. Cell Tissue Res. 291:403–413.
Kelly, R. P., and S. R. Palumbi. 2010. Genetic structure among 50 species of the northeastern Pacific rocky intertidal community. PLoS ONE 5:e8594. doi:10.1371/journal.pone.0008594.

Kultz, D. 1998. Phylogenetic and functional classification of mitogen- and stress-activated protein kinases. J. Mol. Evol. 46:571–588.

Lamaze, F. C., C. Sauvage, A. Marie, D. Garant, and L. Bernatchez. 2012. Dynamics of introgressive hybridization assessed by SNP population genomics of coding genes in stocked brook charr (Salvelinus fontinalis). Mol. Ecol. 21:2877–2895.

Latch, E. K., G. Dharmarajan, J. C. Glaubitz, and O. E. Rhodes Jr. 2006. Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation. Conserv. Genet. 7:295–302.

Lemaire, C., G. Allegrucci, M. Naciri, L. Bahri-Sfar, H. Kara, and F. Bonhomme. 2000. Do discrepancies between microsatellite and allozyme variation reveal differential selection between sea and lagoon in the sea bass (*Dicentrarchus labrax*)? Mol. Ecol. 9:457–467.

Lemaire, C., J. J. Versini, and F. Bonhomme. 2005. Maintenance of genetic differentiation across a transition zone in the sea: discordance between nuclear and cytoplasmic markers. J. Evol. Biol. 18:70–80.

Li, Y. C., A. B. Korol, T. Fahima, and E. Nevo. 2004. Microsatellites within genes: structure, function, and evolution. Mol. Biol. Evol. 21:991–1007.

Limborg, M. T., R. Hanel, P. V. Debes, A. K. Ring, C. André, C. S. Tsigenopoulos, et al. 2012. Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (*Sprattus sprattus*). Heredity 102:96–107.

Luttikhuizen, P. C., J. Dren, K. T. C. A. Peijnenburg, H. C. van der Veer, and K. Johannesson. 2012. Genetic architecture in a marine hybrid zone: comparing outlier detection and genomic clines analysis in the bivalve *Macoma balthica*. Mol. Ecol. 21:3048–3061.

Mallet, J. 2005. Hybridization as an invasion of the genome. Trends Ecol. Evol. 20:229–237.

Marshall, D. J., K. Monro, M. Bode, M. J. Keough, and S. Swearer. 2010. Phenotype–environment mismatches reduce connectivity in the sea. Ecol. Lett. 13:128–140.

Matchkov, V. V., P. Larsen, E. V. Bouzinova, A. Rojek, D. M. B. Boedtkjer, V. Golubinskaya, et al. 2008. Bestrophin-3 (Vitelliform Macular Dystrophy 2-Like 3 Protein) is essential for the cGMP-dependent calcium-activated chloride conductance in vascular smooth muscle cells. Circ. Res. 103:864–872.

McMillan, W. O., L. A. Weigt, and S. R. Palumbi. 1999. Color pattern evolution, assortative mating, and genetic differentiation in brightly colored butterflyfishes (*Chaetodontidae*). Evolution 53:247–260.

Moore, W. S. 1977. An evaluation of narrow hybrid zones in vertebrates. Q. Rev. Biol. 52:263–278.

Naciri, M., C. Lemaire, P. Borsa, and F. Bonhomme. 1999. Genetic study of the Atlantic/Mediterranean transition in seabass (*Dicentrarchus labrax*). J. Hered. 90:591–596.

Narum, S. R. 2006. Beyond Bonferroni: less conservative analyses for conservation genetics. Conserv. Genet. 7: 783–787.

Nei, M. 1987. Molecular evolutionary genetics. Columbia Univ. Press, New York. 512 pp.

Nikula, R., P. Strelkov, and R. Vaïnölä. 2008. A broad transition zone between an inner Baltic hybrid swarm and a pure North Sea subspecies of *Macoma balthica* (Mollusca, Bivalvia). Mol. Ecol. 17:1505–1522.

Nolte, A. W., and D. Tautz. 2010. Understanding the onset of hybrid speciation. Trends Genet. 26:54–58.

Nolte, A. W., J. Freyhof, and D. Tautz. 2006. When invaders meet locally adapted types: rapid moulding of hybrid zones between sculpins (*Cottus*, Pisces) in the Rhine system. Mol. Ecol. 15:1983–1993.

Nolte, A. W., Z. Gompert, and C. A. Buerkle. 2009. Variable patterns of introgression in two sculpin hybrid zones suggest that genomic isolation differs among populations. Mol. Ecol. 18:2615–2627.

Nosil, P., T. H. Vines, and D. J. Funk. 2005. Reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution 59:705–719.

O'Driscoll, K. E., W. J. Hatton, H. R. Burkin, N. Leblanc, and F. C. Britton. 2008. Expression, localization, and functional properties of Bestrophin 3 channel isolated from mouse heart. Am. J. Physiol. Cell Physiol. 295:C1610–C1624.

O'Reilly, P. T. O., M. F. Canino, K. M. Bailey, and P. Bentzen. 2004. Inverse relationship between $F_{ST}$ and microsatellite polymorphism in the marine fish, walleye Pollock (*Theragra chalcogramma*): implications for resolving weak population structure. Mol. Ecol. 13:1799–1814.

Page, R. D. M. 1996. TreeView: an application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12:357–358.

Patarnello, T., F. A. J. M. Volckaert, and R. Castilho. 2007. Pillars of Hercules: is the Atlantic–Mediterranean transition a phylogeographical break? Mol. Ecol. 16:4426–4444.

Payseur, B. A. 2010. Using differential introgression in hybrid zones to identify genomic regions involved in speciation. Mol. Ecol. Resour. 10:806–820.

Piálek, J., and N. H. Barton. 1997. The spread of an advantageous allele across a barrier: the effects of random drift and selection against heterozygotes. Genetics 145:493–504.

Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference of population structure using multilocus genotype data. Genetics 155:945–959.
Quéré, N., B. Guinand, H. Kuhl, R. Reinhardt, F. Bonhomme, and E. Desmarais. 2010. Genomic sequences and genetic differentiation at associated tandem repeat markers in growth hormone, somatolactin and insulin-like growth factor-1 genes of the sea bass, *Dicentrarchus labrax*. Aquat. Living Resour. 23:285–296.

Reynolds, J., B. S. Weir, and C. C. Cockerham. 1983. Estimation of the coancestry coefficient, basis for a short term genetic distance. Genetics 105:767–779.

Roberts, D. G., C. A. Gray, R. J. West, and D. J. Ayre. 2010. Marine genetic swamping: hybrids replace an obligately estuarine fish. Mol. Ecol. 19:508–520.

Rolland, J. L., F. Bonhomme, F. Lagardère, M. Hassan, and B. Guinand. 2007. Population structure of the common sole (*Solea solea*) in the Northeastern Atlantic and the Mediterranean Sea: revisiting the divide with EPIC markers. Mar. Biol. 151:327–341.

Roques, S., J. M. Sévigny, and L. Bernatchez. 2001. Evidence for broadscale introgressive hybridization between two redfish (genus *Sebastes*) in the North-west Atlantic: a rare marine example. Mol. Ecol. 10:149–165.

Sanford, E., and M. W. Kelly. 2011. Local adaptation in the sea. Ann. Rev. Mar. Sci. 3:509–535.

Schunter, C., J. Carreras-Carbonell, E. Macpherson, J. Tintoré, E. Vidal-Vijande, A. Pascual, et al. 2011. Matching genetics with oceanography: directional gene flow in a Mediterranean fish species. Mol. Ecol. 20:5167–5181.

Selkoe, K. A., C. M. Henzler, and S. D. Gaines. 2008. Seascape genetics and the spatial ecology of marine populations. Fish Fish. 9:363–377.

Selkoe, K. A., J. R. Watson, C. White, T. Ben Horin, M. Iacchei, S. Mitarai, et al. 2010. Taking the chaos out of genetic patchiness: seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species. Mol. Ecol. 19:3708–3726.

Shikano, T., J. Ramadevi, and J. Merilä. 2010. Identification of local and habitat-dependent selection: scanning functionally important genes in nine-spined sticklebacks (*Pungitius pungitius*). Mol. Biol. Evol. 27:2775–2789.

Slatkin, M., and T. Wiehe. 1998. Genetic hitch-hiking in a subdivided population. Genet. Res. 71:155–160.

Stemshorn, K. C., F. A. Reed, A. W. Nolte, and D. Tautz. 2011. Rapid formation of distinct hybrid lineages after secondary contact of two fish species (*Cottus* sp.). Mol. Ecol. 20:1475–1491.

Strelkov, P., R. Nikula, and R. Väinölä. 2007. *Macoma balthica* in the White and Barents Seas: properties of a widespread marine hybrid swarm (Mollusca: Bivalvia). Mol. Ecol. 16:4110–4127.

Tsigenopoulos, C. S., B. Hellemans, D. A. Chistiakov, A. Libertini, G. Kotoulas, and F. Volckaert. 2003. Eleven new

microsatellites of the sea bass (*Dicentrarchus labrax* L.). Mol. Ecol. Notes 3:352–354.

Van Oosterhout, C., W. F. Hutchinson, D. P. M. Willis, and P. F. Shipley. 2004. Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4:535–538.

Vasemägi, A., J. Nilsson, and C. R. Primmer. 2005. Expressed sequence tag-linked microsatellites as a source of gene-associated polymorphisms for detecting signatures of divergent selection in Atlantic salmon (*Salmo salar* L.). Mol. Biol. Evol. 22:1067–1076.

Vasioukhin, V. 2006. Lethal giant puzzle of Lgl. Dev. Neurosci. 28:13–24.

Vilas, R., C. Bouza, M. Vera, A. Millán, and P. Martinez. 2010. Variation in anonymous and EST-microsatellites suggest adaptive population divergence in turbot. Mar. Ecol. Prog. Ser. 420:231–239.

Vitalis, R., K. Dawson, and P. Boursot. 2001. Interpretation of variation across marker loci as evidence of selection. Genetics 158:1811–1823.

Vitalis, R., K. Dawson, P. Boursot, and K. Belkhir. 2003. DetSel 1.0: a computer program to detect markers responding to selection. J. Hered. 94:429–431.

Volckaert, F. A. M., C. Batargias, A. Canario, D. Chatziplis, D. Chistiakhov, C. Haley, et al. 2008. European sea bass. Pp. 117–133 in T. D. Kocher and C. Cole, eds. Genome mapping and genomics in animals. Vol. 2: Genome mapping and genomics in fishes and aquatic animals. Springer-Verlag, Berlin.

Weihe, E., J. H. Tao-Cheng, M. K. H. Schaefer, J. D. Erickson, and L. E. Eiden. 1996. Visualization of the vesicular acetylcholine transporter in cholinergic nerve terminals and its targeting to a specific population of small synaptic vesicles. Proc. Natl Acad. Sci. USA 93:3547–3552.

Weir, B. S., and C. C. Cockerham. 1984. Estimating *F*-statistics for the analysis of population structure. Evolution 38:1358–1370.

Willett, C. S. 2012. Hybrid breakdown weakens under thermal stress in population crosses of the copepod *Tigrionus californicus*. J. Hered. 103:103–114.

Wodarz, A. 2000. Tumor suppressors: linking cell polarity and growth control. Curr. Biol. 10:R624–R626.

Yatabe, Y., N. C. Kane, C. Scotti-Saintagne, and L. H. Rieseberg. 2007. Rampant gene exchange across a strong reproductive barrier between the annual sunflowers, *Helianthus annuus* and *H. petiolaris*. Genetics 175:1883–1893.

Zulliger, D. E., S. Tanner, M. Ruch, and G. Rubi. 2009. Genetic structure of the high dispersal Atlanto-Mediterranean sea star *Astropecten aranciacus* revealed by mitochondrial DNA sequences and microsatellite loci. Mar. Biol. 156:597–610.
**Supporting Information**

Additional Supporting Information may be found in the online version of this article at the publisher's website:

**Figure S1.** Single-locus estimates of *F*<sub>st</sub> as estimated by *θ* (mean ± 95% CI) vs. expected heterozygosity for (a) the 20 loci used for the Western Mediterranean samples, and

(b) the 15 loci used for the Eastern Mediterranean samples. Relationship between *θ* and *H*<sub>exp</sub> revealed no significant negative correlations (Western Mediterranean: *R*<sup>2</sup> = 0.062, NS; Eastern Mediterranean: *R*<sup>2</sup> = 0.129, NS). Open and filled diamonds indicate GIL and GAL, respectively. Labels of individual loci are not reported.