# Implications for management and conservation of the population genetic structure of the wedge clam Donax trunculus across two biogeographic boundaries

Amandine D. Marie, Christophe D Lejeusne, Evgenia A Karapatsiou, José A. Cuesta, Pilar A Drake, Enrique A Macpherson, Louis A Bernatchez, Ciro Rico

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

## To cite this version:

Amandine D. Marie, Christophe D Lejeusne, Evgenia A Karapatsiou, José A. Cuesta, Pilar A Drake, et al.. Implications for management and conservation of the population genetic structure of the wedge clam Donax trunculus across two biogeographic boundaries. Scientific Reports, 2016, 6, pp.39152. 10.1038/srep39152. hal-01436790

HAL Id: hal-01436790

https://hal.sorbonne-universite.fr/hal-01436790

Submitted on 16 Jan 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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
SCIENTIFIC REPORTS

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

OPEN

# Implications for management and conservation of the population genetic structure of the wedge clam *Donax trunculus* across two biogeographic boundaries

Received: 20 May 2016

Accepted: 18 November 2016

Published: 19 December 2016

Amandine D. Marie<fnref n="1" />,*, Christophe Lejeusne²,³,⁴, Evgenia Karapatsiou⁵, José A. Cuesta⁶, Pilar Drake⁶, Enrique Macpherson⁷, Louis Bernatchez⁸ & Ciro Rico<fnref n="1" />,²,*

In a resource management perspective, the understanding of the relative influence of the physical factors on species connectivity remains a major challenge and is also of great ecological and conservation biology interest. Despite the overfishing threat on the wedge clam *Donax trunculus* in Europe, relatively little information is known about its population genetic structure and connectivity and their consequences on conservation policies. We employed 16 microsatellite loci to characterise the genetic diversity and population structure of *D. trunculus*. A total of 514 samples from seven different localities along the Atlantic-Mediterranean transition, from the Atlantic (Gulf of Cádiz) to the north-western Mediterranean were genotyped. The analysis of the population genetic structure displayed a clear distinction along the Atlantic-Mediterranean transition with different clusters in the Atlantic Ocean, the Alboran Sea and the north-western Mediterranean. Consequently, we recommend that these three areas should be considered as different management units. We showed that all populations seem to be at high long-term risk of extinction with the exception of the protected Doñana National Park population which still seems to have evolutionary potential. Therefore, our results emphasized the necessity of protection of this economic resource and the validity of molecular tools to evaluate the population dynamics.

Many soft and hard faunal barriers to dispersal have been identified in marine environments throughout the world, including, for instance, the Isthmus of Panama<fnref n="1" />⁻³ and the Terminal Tethyan Event⁴,⁵ as hard barriers, as well as the East Pacific Barrier⁶,⁷ and the Benguela upwelling⁸ as soft barriers. Some potential physical factors such as ocean currents, habitats discontinuities or temperature gradients are recognised as driving the biogeographic barriers structuring these regions. Thus, in a resource management perspective, the understanding of the relative influence of the physical factors on species connectivity remains a major challenge and is also of great ecological and conservation biology interest.

Genetic studies have demonstrated that gene flow between conspecifics either side of barrier is often limited for marine⁹⁻¹¹ and terrestrial fauna¹²⁻¹⁴. Among the major marine biogeographic boundaries, the Strait of Gibraltar that separates the Mediterranean Sea from the Atlantic Ocean and the Almería Oran Front, which is a semi-permanent dynamic oceanographic front connecting the main jet of incoming Atlantic water and the

<footnote n="1">¹School of Marine Studies, Molecular Diagnostics Laboratory, Faculty of Science Technology and Environment, The University of South Pacific, Suva, Fiji Islands. ²Estación Biológica de Doñana, Consejo Superior de Investigaciones Científicas (EBD, CSIC), Sevilla 41092, Spain. ³Sorbonne Universités, UPMC Université Paris 06, UMR 7144, Station Biologique de Roscoff, 29688 Roscoff, France. ⁴CNRS, UMR 7144, Station Biologique de Roscoff, 29688 Roscoff, France. ⁵University of Aegean, Department of Environmental Sciences, University, Hill, Mytilene 81100, Greece. ⁶Instituto de Ciencias Marinas de Andalucía, Consejo Superior de Investigaciones Científicas (ICMAN, CSIC), Cádiz 11519, Spain. ⁷Centre d'Estudis Avançats de Blanes, Consejo Superior de Investigaciones Científicas (CEAB, CSIC), Blanes 17300, Spain. ⁸Institut de Biologie Intégrative et des Systèmes (IBIS), Département de Biologie, Université Laval, Pavillon Charles-Eugène-Marchand, Québec G1V 0A6, Canada. ⁹These authors contributed equally to this work. Correspondence and requests for materials should be addressed to A.D.M. (email: amandine.marie@usp.ac.fj)</footnote>
Mediterranean Sea, have been confirmed by several studies as effective barriers to dispersal<sup>11,15,16</sup>. The complexity of historical and contemporary physical factors may explain the causes of the reduced dispersal encountered in this region. It is also well established that the Mediterranean Sea was strongly affected by sea level fluctuations during the Messinian Salinity Crisis in the Miocene and subsequently during the Pleistocene glaciations, isolating it several times from the Atlantic<sup>15</sup>. However, although the actual reopening of the Gibraltar Strait allowed free migration across this area, genetic patterns of some marine species reveal reduced gene flow levels between populations at either side, confirming that it represents a major biogeographic barrier for many taxa. Strong genetic differentiation between Atlantic and Mediterranean populations has been shown in fish as well as in shellfish species with contrasting life histories<sup>10,17,18</sup>. To explain this phenomenon, two mains factors have been proposed: 1) the one-way surface current of Atlantic water flowing through the Strait of Gibraltar into the Mediterranean<sup>19</sup> and, 2) the presence of a front of surface waters between Almeria in Southeast of Spain and Oran in Algeria, the so-called Almeria-Oran front (AOF)<sup>16</sup>. These factors could represent a barrier to gene flow between Atlantic and Mediterranean populations and be responsible for the reduced migration observed.

The wedge clam, *Donax trunculus* (Bivalvia: Donacidae), is an Atlanto-Mediterranean warm-temperate species found in the Black Sea, the Mediterranean Sea<sup>20</sup> and from Senegal to the northern Atlantic coast of France<sup>21</sup>. *D. trunculus* constitutes a very important fishing resource due to its high economic value. Only in Europe, the recorded landings over the last 12 years equal 11,202 tons, with a maximum yield of 1,355 tons in 2005 followed by a steady decline reaching only 516 tons in 2011<sup>22</sup>. Although landing data are only available since 1972 and fishermen were not obliged to declare their catches until 1998<sup>23</sup>, the species has been intensively exploited over the last decades using primarily either artisanal hand-operated dredges in shallow water or commercial boat-operated dredges from subtidal to deeper waters<sup>24</sup>.

Until recently<sup>25,26</sup> and despite the overexploitation threat faced by this shellfish resource, *D. trunculus* has not been studied from a conservation biology or genetics perspective. The few studies that have paid any attention to its genetic makeup have focused on its genome structure<sup>27</sup> or its mitochondrial DNA inheritance<sup>25,28</sup> only, but not with the aim to characterise its diversity. However, the availability of highly polymorphic codominant genetic markers, such as microsatellites, allows essential research on the population structure, its effective population size, contemporary phylogeography and conservation genetics of wedge clamps. Nanton *et al.*<sup>26</sup> and Rico *et al.* (GenBank Accession Numbers HG792255 to HG792275) recently characterised respectively the first nineteen and sixteen polymorphic microsatellite markers in *D. trunculus*, providing a first battery of genetic tools for the study of the genetic structure of the species to allow the delineation of conservation units essential in fisheries management<sup>29</sup>.

Understanding connectivity among populations through genetic structure has important conservation implications because genetic assessment is a central criterion in determining the appropriate units and spatial scale for conservation and management<sup>30,31</sup>. Thus, in a resource management perspective, the aim of this study was to assess the genetic diversity as well as the genetic structure of population of the wedge clam *Donax trunculus*. For that, 514 individuals from 7 locations along the Atlantic and Mediterranean coasts of Spain were genotyped using 16 microsatellite markers. After evaluating the genetic diversity, the importance of the two well-known characterised biogeographic boundaries (i.e. the Strait of Gibraltar and the Almeria-Oran Front) on the population genetic structure was assessed. Finally, we discuss these findings in a perspective of biodiversity conservation and resource management.

## Results

**Genetic diversity.** A total of 514 individuals from seven different localities were genotyped at 16 microsatellite loci. The genetic diversity found in each locus across samples is summarised in Table 1. Highly significant deviations from Hardy-Weinberg equilibrium were found in all sampling sites and for each locus. Exact tests for genotypic linkage disequilibrium confirmed the absence of linkage among most loci (91.6%; *P* < 0.05 after Bonferroni correction). Genotyping was consistent across the 514 samples as revealed by mis-scoring rates from two independent scorers below 2%. The number of alleles (N<sub>A</sub>) per locus ranged from 12 to 33 (χ = 22.75, σ = 6.8). Observed (H<sub>O</sub>) and expected (H<sub>E</sub>) heterozygosity and polymorphic information content ranged from 0.207 to 0.768 (χ = 0.49, σ = 0.18), 0.557 to 0.947 (χ = 0.81, σ = 0.11) and 0.534 to 0.942 (χ = 0.79, σ = 0.12), respectively. The software Microchecker<sup>32</sup> depicted the occurrence of null alleles in all loci and the presence of large alleles dropout for some loci. Thus, to better estimate the frequency of null alleles, we relied on the software Cervus<sup>33</sup> (Table 1). The null allele frequencies results from Cervus revealed an estimate for each locus ranging from 6.6% to 55.3% (χ = 27.3, σ = 0.16). and all loci, but two (Dtru 16 and Dtru 23), had a null alleles frequency above 10% (Table 1). Table 2 summarises the genetic diversity found for each sampling site. Observed and expected heterozygosity ranged from 0.423 to 0.533 (χ = 0.49, σ = 0.04) and from 0.756 to 0.826 (χ = 0.79, σ = 0.03), respectively. Despite the strong departure from HWE observed, we estimated the effective population size, using NeEstimator to explore current demographic patterns in localities exposed to different fishing pressure. The Ne estimates ranged from 210 for Isla Canela to 1191.7 for Doñana (95% confidence interval of the parametric test ranging from 165.1 to infinite; Table 2).

**Genetic structure and gene flow.** The analysis of molecular variance (AMOVA) revealed that the majority (98.6%) of the observed variation was explained by differences within the sampling sites, however a small but significant portion of the variation (1.1%) occurred among sampling sites between groups (Table 3). Results of the F<sub>ST</sub> analysis revealed a significant genetic differentiation (*P* < 0.05) between the Gulf of Cádiz (Southern Atlantic Spanish coast), the Alboran Sea (Southwestern Mediterranean coast) and the Northwestern Mediterranean coast. This corroborates the presence of two well characterised biogeographical barriers, the Strait of Gibraltar<sup>34</sup> and the AOF (*P* > 0.05; Table 4).
<table><caption>**Table 1.** Details of the 16 microsatellite markers developed and optimized. 514 individuals from 7 sampling sites were used for these analyses. Locus identity, (N<sub>A</sub>) number of alleles found, (Range) range of allele sizes, (H<sub>O</sub>) observed and (H<sub>E</sub>) expected heterozygosity, (PIC) polymorphism information content, (NAF) null-allele frequency and GenBank accession numbers.</caption><thead><tr><th>Locus ID</th><th>NA</th><th>Range</th><th>HO</th><th>HE</th><th>PIC</th><th>NAF</th><th>GenBank Acc. N°</th></tr></thead><tbody><tr><td>D.tru2</td><td>18</td><td>97-148</td><td>0.590</td><td>0.830</td><td>0.809</td><td>0.133</td><td>HG792255</td></tr><tr><td>D.tru4</td><td>14</td><td>151-179</td><td>0.641</td><td>0.750</td><td>0.716</td><td>0.123</td><td>HG792256</td></tr><tr><td>D.tru6</td><td>29</td><td>69-133</td><td>0.408</td><td>0.796</td><td>0.777</td><td>0.327</td><td>HG792257</td></tr><tr><td>D.tru8</td><td>14</td><td>127-169</td><td>0.670</td><td>0.823</td><td>0.796</td><td>0.109</td><td>HG792258</td></tr><tr><td>D.tru11</td><td>26</td><td>97-153</td><td>0.469</td><td>0.715</td><td>0.700</td><td>0.211</td><td>HG792259</td></tr><tr><td>D.tru14</td><td>19</td><td>75-135</td><td>0.207</td><td>0.557</td><td>0.534</td><td>0.540</td><td>HG792260</td></tr><tr><td>D.tru15</td><td>33</td><td>215-279</td><td>0.592</td><td>0.947</td><td>0.942</td><td>0.218</td><td>HG792261</td></tr><tr><td>D.tru16</td><td>26</td><td>70-124</td><td>0.768</td><td>0.912</td><td>0.903</td><td>0.093</td><td>HG792262</td></tr><tr><td>D.tru19</td><td>28</td><td>126-184</td><td>0.542</td><td>0.935</td><td>0.929</td><td>0.223</td><td>HG792263</td></tr><tr><td>D.tru22</td><td>12</td><td>215-249</td><td>0.280</td><td>0.689</td><td>0.639</td><td>0.410</td><td>HG792264</td></tr><tr><td>D.tru23</td><td>33</td><td>101-167</td><td>0.738</td><td>0.835</td><td>0.825</td><td>0.066</td><td>HG792265</td></tr><tr><td>D.tru26</td><td>28</td><td>55-119</td><td>0.590</td><td>0.924</td><td>0.917</td><td>0.223</td><td>HG792266</td></tr><tr><td>D.tru29</td><td>21</td><td>262-322</td><td>0.261</td><td>0.826</td><td>0.808</td><td>0.487</td><td>HG792268</td></tr><tr><td>D.tru32</td><td>25</td><td>196-284</td><td>0.323</td><td>0.674</td><td>0.659</td><td>0.425</td><td>HG792271</td></tr><tr><td>D.tru40</td><td>23</td><td>119-167</td><td>0.578</td><td>0.889</td><td>0.878</td><td>0.229</td><td>HG792273</td></tr><tr><td>D.tru49</td><td>15</td><td>144-182</td><td>0.238</td><td>0.801</td><td>0.771</td><td>0.553</td><td>HG792275</td></tr></tbody></table>

<table><caption>**Table 2.** Characteristics of the sampling sites including the number of clams genotyped (N), observed (H<sub>O</sub>) and expected (H<sub>E</sub>) heterozygosity and the genetic effective population size (Ne estimate) with the lower (lower CI) and the higher (higher CI) confidence intervals. The use of ∞ indicates inestimable upper confidence limits*<sup>56</sup>.</caption><thead><tr><th></th><th>N</th><th>Ho</th><th>He</th><th>Ar</th><th>Ne estimate</th><th>lower CI</th><th>higher CI</th><th>Situation*</th></tr></thead><tbody><tr><td>Isla Canela</td><td>85</td><td>0.533</td><td>0.799</td><td>14.327</td><td>210</td><td>165.1</td><td>238.6</td><td>at high long-term risk of extinction</td></tr><tr><td>Doñana</td><td>86</td><td>0.517</td><td>0.817</td><td>14.689</td><td>1191.7</td><td>526.9</td><td>∞</td><td>at evolutionary potential</td></tr><tr><td>Caleta de Vélez</td><td>56</td><td>0.516</td><td>0.773</td><td>13.109</td><td>782.5</td><td>320</td><td>∞</td><td>at high long-term risk of extinction</td></tr><tr><td>Cabo de Gata</td><td>86</td><td>0.463</td><td>0.826</td><td>13.573</td><td>357.2</td><td>245</td><td>635.6</td><td>at high long-term risk of extinction</td></tr><tr><td>Gandía</td><td>53</td><td>0.468</td><td>0.769</td><td>12.577</td><td>332.5</td><td>203.3</td><td>845.1</td><td>at high long-term risk of extinction</td></tr><tr><td>Sant Carles de la Ràpita</td><td>62</td><td>0.423</td><td>0.758</td><td>11.136</td><td>602.7</td><td>282.9</td><td>∞</td><td>at high long-term risk of extinction</td></tr><tr><td>Roses</td><td>86</td><td>0.486</td><td>0.756</td><td>15.202</td><td>600.3</td><td>340.4</td><td>2210.8</td><td>at high long-term risk of extinction</td></tr></tbody></table>

<table><caption>**Table 3.** Results of an analysis of molecular variance (AMOVA) showing the distribution of the genetic variation among groups, among populations within groups and within populations based on 16 microsatellites. Gulf of Cádiz (Atlantic coast): Isla Canela and Doñana. Alboran Sea (southwestern Mediterranean coast): Caleta de Vélez and Cabo de Gata. Northwestern Mediterranean coast: Gandía, Sant Carles de la Ràpita and Roses. *Significant P-value (<0.05).</caption><thead><tr><th>Source of variation</th><th>Sum of squares</th><th>Variance components</th><th>Percentage variation</th><th></th></tr></thead><tbody><tr><td>Among groups</td><td>21.734</td><td>0.023</td><td>1.08*</td><td>F<sub>CT</sub> = 0.011*</td></tr><tr><td>Among populations within groups</td><td>12.101</td><td>0.007</td><td>0.31</td><td>F<sub>SC</sub> = 0.003</td></tr><tr><td>Within populations</td><td>2127.061</td><td>2.083</td><td>98.60*</td><td>F<sub>ST</sub> = 0.014*</td></tr><tr><td>Total</td><td>2160.896</td><td>2.113</td><td></td><td></td></tr></tbody></table>

STRUCTURE analyses using the RECESSIVEALLELES and LOCPRIOR functions further substantiated the existence of three genetically differentiated clusters: one in the Gulf of Cádiz; a second one for the samples collected east of the Strait of Gibraltar, in the Alboran Sea, but west of AOF in the Mediterranean Sea; and a third one
<table><caption>**Table 4.** Pairwise *F*<sub>ST</sub> (above diagonal) and associated *P*-value (below diagonal) among pooled sampling sites. Gulf of Cadiz (Atlantic coast); Isla Canela and Doñana. Alboran Sea (southwestern Mediterranean coast): Caleta de Vélez and Cabo de Gata. Northwestern Mediterranean Coast: Gandía, Sant Carles de la Ràpita and Roses. *Significant *P*-value.</caption><thead><tr><th></th><th>Gulf of Cádiz</th><th>Alboran Sea</th><th>Northwestern Mediterranean coast</th></tr></thead><tbody><tr><td>Gulf of Cádiz</td><td>-</td><td>0.018</td><td>0.032</td></tr><tr><td>Alboran Sea</td><td>0.017*</td><td>-</td><td>0.028</td></tr><tr><td>Northwestern Mediterranean coast</td><td>0.017*</td><td>0.017*</td><td>-</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Assignment of individuals clams to each pooled sampling site in a conventional bar plot (results obtained from STRUCTURE). The red cluster identifies the population of the Gulf of Cádiz, the green cluster the population of the Alboran Sea and the blue cluster the population of the Northwestern Mediterranean Coast.</figcaption></figure>

for the samples collected northwest of the AOF in the Northwestern Mediterranean coast. A clear discrimination was obtained between these three clusters (Fig. 1). The assigned *q*-value for each pooled sampling site were of 0.806 ± 0.080 (range from 0.462 to 0.920) for the Gulf of Cadiz, 0.542 ± 0.084 (range from 0.357 to 0.736) for the Alboran Sea, and 0.845 ± 0.100 (range from 0.524 to 0.959) for the Northwestern Mediterranean Coast. Some individuals were not included into the mean individual admixture proportions of each pooled sampling site because they clustered into a different group: one individual from Isla Canela (*q*-value = 0.589) clustered into the Gulf of Cadiz and, two individuals from Caleta de Vélez (*q*-value = 0.438 and *q*-value = 0.355) and two from Cabo de Gata (*q*-value = 0.444 and *q*-value = 0.458) into the Northwestern Mediterranean Coast.

Analyses using BayesAss suggest a consistent gene flow among the groups of populations (Table 5). The average migration rate in all pairwise comparisons ranged from 0.003 to 0.034. Each value represents the proportion of individuals that is derived from a corresponding source population for each generation. Altogether, the average migration rate was very low between samples from the Alboran Sea (i.e. source) and the two other groups of populations (i.e. sink), 0.034 and 0.017, respectively for the Gulf of Cadiz and the Northwestern Mediterranean Coast and were higher than any other pairwise comparison.

## Discussion

Several recent studies have clearly shown that the degree of connectivity among marine taxa separated by biogeographic barriers varies from extensive isolation to complete panmixia and thus the actual influence of these marine transitions remains controversial³⁵. Furthermore, they have also shown that the differences observed are largely dependent on early life history traits, such as pelagic larval duration of the species concerned¹¹,³⁶. In fact, the number of studies that have utilised different molecular markers to unravel the contrasting results of the influence of biogeographical barriers to dispersal in marine environments is so vast that a new synthesis is required (e.g.³⁷). The failure to account for differences in heterozygosity among markers and species may explain to some extent for the irreconcilable patterns observed across ecological and evolutionary timescales³⁸. For example, a recent study on European sea bass *Dicentrarchus labrax* confirmed that the Atlantic and Mediterranean basins harbour two distinct lineages³⁹. Similarly, Luttikhuizen *et al.*¹⁰ showed a very clear phylogeographic differentiation among the common shrimp (*Crangon crangon*) populations from the Atlantic Ocean and the Mediterranean Sea. Very high levels of genetic differentiation between Atlantic and Mediterranean populations (pairwise θ<sup>B</sup> values range from 0.395 to 0.510) were also found in the Lusitanian sea star (*Asterina gibbosa*)¹⁸. On the other hand, the analysis of the genetic structure among samples from Atlantic Ocean and Mediterranean Sea of the European
<table><caption>**Table 5.** Estimates of migration rates (proportion of individuals). Values shown are the means of the posterior distributions of *m*, the migration rate into each group of populations, and their respective 95% confidence intervals in parentheses. Values in bold are the proportion of individuals derived from the source population each generation. Column headings are the source population, and row headings are the destination populations. Gulf of Cádiz (Atlantic coast): Isla Canela and Doñana. Alboran Sea (Southwestern Mediterranean Coast): Caleta de Vélez and Cabo de Gata. Northwestern Mediterranean Coast: Gandía, Sant Carles de la Ràpita and Roses.</caption><thead><tr><th rowspan="2">Destination</th><th colspan="3">Source</th></tr><tr><th>Gulf of Cadiz</th><th>Alboran Sea</th><th>Southwestern Mediterranean coast</th></tr></thead><tbody><tr><td>Gulf of Cadiz</td><td><strong>0.985 (0.007)</strong></td><td>0.009 (0.005)</td><td>0.006 (0.004)</td></tr><tr><td>Alboran Sea</td><td>0.032 (0.014)</td><td><strong>0.952 (0.016)</strong></td><td>0.017 (0.009)</td></tr><tr><td>Southwestern Mediterranean coast</td><td>0.003 (0.003)</td><td>0.012 (0.006)</td><td><strong>0.985 (0.006)</strong></td></tr></tbody></table>

spiny lobster (*Palinurus elephas*) showed small but significant differences ($F_{ST}$ for pairwise comparisons within basins = 0.003 ± 0.004 and $F_{ST}$ among basins = 0.011 ± 0.005) indicating the presence of two partially-overlapping groups (overlapping of some localities from the Atlantic and Mediterranean basins)$^{40}$. In contrast, a recent study on the genetic structure the holoplanktonic jellyfish, *Pelagia noctiluca*, across the Atlantic Ocean and Mediterranean Sea indicated a high degree of connectivity with little evidence of population subdivision between basins$^{41}$. A very small but significant differentiation ($F_{ST}$ = 0.005) among Atlantic and Mediterranean stocks was detected when using microsatellite allele frequencies but no evidence of differentiation was observed using mtDNA haplotype diversity. The magnitude of genetic differentiation observed between the different pooled sampling sites in the wedge clam can be considered as large values for marine taxa (from 2 to 3%).

In this study, the level of genetic differentiation observed between the main oceanographic discontinuities along the Gulf of Cádiz to the Western Mediterranean Sea (the Gibraltar Strait and the AOF) confirms that these biogeographic barriers to dispersal effectively restrict gene flow between *D. trunculus* populations as supported by our Bayesass analysis and also as reported in many other taxa$^{11}$. Before the entry of the Atlantic waters throughout the Gibraltar Strait a branch of these waters recirculates near the Strait, in front of the Cape of Trafalgar, towards the northwest along the coast of Cádiz. This area is also influenced by the intense tidal-current regime of the Strait of Gibraltar and the strong topographic interaction between the swift along-shore tidal flow and a submerged ridge running perpendicular to the shoreline$^{15,42}$. These processes originate persistently a patch of cold water that can affect the connectivity between populations at both sides of the Gibraltar Strait$^{11,43}$. On the other hand, the AOF is a semi-permanent dynamic oceanographic front connecting the main jet of incoming Atlantic water and the Mediterranean Sea and it is situated at approximately 400 km east of the Strait of Gibraltar$^{16}$. This front is formed by the convergence between two different water masses (of Atlantic and Mediterranean origin) present in the upper 300 m, and it is known to act as a barrier to gene flow in numerous species$^{11,15}$. This phenomenon is clearly demonstrated through our results (Table 4 and Fig. 1), which displayed a significant genetic differentiation between the populations from the Gulf of Cádiz, the Alboran Sea and the Northwestern Mediterranean Coast. As such, these three areas should be considered as different management units. However, the levels of genetic differentiation obtained in this study might represent a conservative estimate as levels of heterozygosity may be artificially decreased due to the presence of null alleles. However, despite the heterozygous deficits of the loci used in this study, they were all sufficient for clearly differentiating the populations between the Atlantic Ocean and the Southeast and Northwest Mediterranean Coast of Spain.

The high incidence of null alleles found in this study is consistent with many reports in natural populations of bivalve species analysed with DNA markers such as microsatellites$^{26,44,45}$. In a previous study on *D. trunculus*, Nanton *et al.*$^{26}$ showed the presence of null alleles (range from 0.109 to 0.277) in 10 loci out of 19. Another study aiming to examine the genetic variability and relationships between two shellfish (*Ruditapes philippinarum*) populations in Korea also reported the occurrence of null alleles in 7 loci out of 10$^{46}$. The presence of null alleles has also been reported in a study that assessed the genetic diversity of four wild and six hatchery stocks of the hard clam (*Mercenaria mercenaria*) in Florida. Null alleles were present in all loci (range from 0.015 to 0.296), but for 4 loci out of 7, the frequency of null alleles was particularly high (>0.146)$^{47}$. Although the causes of these heterozygote deficits remain uncertain, several hypotheses have been suggested to explain them including the Wahlund effect, selection, inbreeding and the presence of null alleles$^{48}$. Inbreeding seems unlikely given the large population size of bivalves. Moreover, we cannot demonstrate from these results, selection against heterozygotes. Indeed, although microsatellite markers are supposed to be neutral, it may be possible that one or more loci are linked to genes under selection$^{49}$ or are under selective constrains due to DNA structure$^{50}$. Another reason that may explain the presence of heterozygote deficits would be scoring errors. However, we can confidently discard scoring errors as a cause of heterozygous deficits because the robustness of the genotyping was thoroughly assessed by double checking independently all genotypes by two experienced scorers and secondly by using the software Microchecker (unpublished data).

Alternatively and certainly a far more interesting result would be that the deviations from HWE are due to the known severe exploitation that the species has been subjected over the last decades. For example, Hoarau *et al.*$^{51}$ used DNA from archived otoliths collected between 1924 and 1972 together with 2002 juvenile's tissue to estimate effective population size (Ne) in plaice (*Pleuronectes platessa*). They found that populations examined between 1924 and 1960 were in Hardy-Weinberg equilibrium, whereas populations examined after approximately 1970
were not. They performed extensive testing to rule out that these deviations for HWE were not due to genotyping artefacts and Wahlund effects. Subsequently, the authors were able to attribute the significant heterozygote deficiencies found from 1970 onward to inbreeding and overexploitation. Unfortunately, we do not have access to archived samples to verify if the populations where in HWE in the past. However, we attempted to rule out genotyping errors (see above) and PCR artefacts using the following approaches. First, we redesigned primers whenever the flanking sequences where adequate (D.tru2, D.tru8, D.tru14, D.tru28, and D.tru32). Secondly, we re-amplified 10 loci including the loci for which we redesigned primers and 5 more which produced codominant amplicons for a subsample of individuals (approx. 200) under relaxed annealing temperatures (4°C lower) to verify that additional alleles had not been missed. None of these strategies resulted in a reduction of the HWE deviations and thus, it was concluded that PCR artefacts were minimal and insufficient to explain the large heterozygous deficits observed. While it is impossible with the current dataset and with the funding that as available for the project to demonstrate that overexploitation may play and important part in the lack of conformance with HWE, this explanation remains plausible and interesting.

In this study, the microsatellite loci characterized for *D. trunculus* showed variable levels of genetic variation than those previously found²⁶. The allele number of these polymorphic markers ranged from 12 to 33, which is higher than previous results on *D. trunculus*²⁶ as well as the most part of the studies on other bivalve species such as *Sinonovacula constricta*⁵² and *Mactra chinensis*⁵³. A heterozygosity deficit has been observed for each locus and all sampling sites (Table 1; Table 2). Although these results are higher than those reported by Nanton *et al.*²⁶, and are consistent with the levels of variability found in both for *D. trunculus* and other clam species⁵⁴.

The effective population size (Ne) is an important parameter of numerically small populations, which is positively correlated with their ability to persist in a changing environment and to evolve. The 50/500 rule often cited by conservation practitioners⁵⁵ has been recently revised to ≥100/1000 by Frankham *et al.*⁵⁶. They estimated that a minimum of 100 for Ne was necessary to avoid immediate risk of extinction due to inbreeding depression. Populations with Ne less than 500–1000 were at high long-term risk of extinction due to the loss of ecologically relevant genetic diversity by drift and populations with Ne above 1000 were considered as being able to retain their evolutionary potential. *D. trunculus* populations have been intensively exploited over the last decades, and particularly in Europe²⁴, and thus, all populations seem to be at high long-term risk of extinction except in Doñana, a managed area located in a National Park where effective population sizes appear to be larger (Table 2)⁵⁷. However, although we used the $P_{crit}$ value recommended for microsatellite markers ($P_{crit} = 0.02_{58}$), these results should be taken with caution given the presence of null alleles, and as a consequence of the heterozygote deficits. Interestingly however, Isla Canela suffers from severe overexploitation from both, registered and licensed fishermen and from Illegal, Unregulated and Unreported (IUU) extraction and has the lowest Ne estimate. This result supports to some extent the robustness of the analyses⁵⁷.

Both populations, Isla Canela and Doñana, are geographically very close to each other (85 km) and share a similar habitat structure: long sandy beaches close to a river mouth (Guadiana and Guadalquivir Rivers, respectively) (Fig. 2). However, they also present clear differences in respect to the management of the *Donax trunculus* fishery. The Doñana beach belongs to the Doñana National Park, and since 1996 has a regulation of the *D. trunculus* fishery: only 160 shellfish collectors are allowed (have a license). These shellfish collectors have to work using only an artisanal rake with a long net attached to the end of the drag (fishing from boats is forbidden) on the intertidal zone. During five days a week, they can collect a maximum of 25 kg of *D. trunculus* individuals per day and these must be larger than 25 mm. Furthermore, they have to respect the breeding season (two months a year) without fishing, and they cannot collect when toxic algae are detected in the wedge clams. In the case of Isla Canela, in addition to the shellfish collectors with license that work in the *D. trunculus* fishery, there are numerous illegal collectors who operate in the area throughout the year and sell their catches in the streets of Huelva and Seville and other smaller villages in both Provinces of Andalusia and, in general, visitors to the beach also harvest clams mostly during the summer months (CR & JAC personal observations). Fishing from boats is also allowed, so the area of collection is wider, covering both the subtidal and intertidal zones. Moreover, there is neither a strict and respected reproductive period without fishery, nor a strict control over the *D. trunculus* minimum sizes. All these differences in the exploitation of *D. trunculus* fishery are reflected in the size structure of both populations. For example, in July of 2011, 2012 and 2014, although significant differences in *D. trunculus* density between the two beaches were only observed in 2014, individuals from Doñana were always significantly larger than those from Isla Canela (Supplementary Methods and Results; Figs 1 and 2), possibly reflecting different age structures resulting from over-harvesting in Isla Canela. This observation is also in line with the observed differences in effective population size between these locations.

In conclusion, for the first time, the genetic diversity and population structure of *Donax trunculus* is described along the Atlantic-Mediterranean transition. Using 16 microsatellite markers, we identified the presence of a large proportion of null alleles in this species. However, the presence of null alleles in the loci appears to have no effect on the population genetic structure uncovered. Our results showed that all populations seem to be at high long-term risk of extinction with the exception of the Doñana National Park population which still seems to have evolutionary potential. Interestingly, this population is rigorously managed by the authorities of the National Park with restrictions in the number of fishing licences, the amount of allowed harvest and the use of hand operated dredges. This work represents a contribution for the understanding of the genetic diversity and population structure of this species, and de facto in terms of biodiversity conservation and resource management.

## Methods

**Sample collection.** With the assistance of a licensed fisherman and using both artisanal hand- and boat-operated dredges, *Donax trunculus* individuals were obtained in seven different locations in the Month of July 2011 in Spain (Fig. 2): Gulf of Cádiz/Atlantic coast (Isla Canela and Doñana National Park), Alboran Sea/Southwestern Mediterranean Coast (Caleta de Vélez and Cabo de Gata) and Northwestern Mediterranean Coast (Gandía, Sant
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Geographical location of sampling sites on the East Atlantic Ocean and Western Mediterranean Sea Coast in Spain (Google Earth, 2016, version 7.1.5.1557, https://www.earth.google.com).</figcaption></figure>

Carles de la Ràpita and Roses). All the samples used in this study belong to species targeted by the commercial fishing industry and are neither protected, nor require approval of any animal ethics committee since they are invertebrates. Permit of sampling in the Doñana National Park was granted for the study.

**DNA extraction.** DNA was extracted from different wedge clam's tissues using a slightly modified version of Aljanabi and Martinez⁵⁹ salting out method. Namely, after addition of the saline solution, the mixture was centrifuged for 30 min at 10,000 g. Also, DNA precipitation was achieved by incubation at -20 °C with 600 µL isopropanol for 30 min. then, after washing out the pellet with 70% ethanol, we centrifuged at 10,000 g for 10 min.

**Microsatellite analysis.** A total of 514 of *Donax trunculus* individuals were genotyped with 16 microsatellite markers (D.tru2, D.tru4, D.tru6, D.tru8, D.tru11, D.tru14, D.tru15, D.tru16, D.tru19, D.tru22, D.tru23, D.tru26, D.tru29, D.tru32, D.tru40, D.tru49). The PCR amplifications were done on an Applied Biosystems 9700 DNA thermal cycler. The PCR products were then separated on a capillary sequencer (ABI 3130x Genetic Analyzer, Applied Biosystems, USA) using GeneScan™ 500 LIZ® Size Standard. Allele sizes were determined with the Gene Mapper®V4.0 program (Applied Biosystems, USA).

**Genetic diversity analyses.** Deviation from Hardy-Weinberg equilibrium and linkage disequilibrium were tested using Genepop (version 4.0⁶⁰ for each marker. We used the software MICROCHECKER³² to test the presence of genotyping errors due to stuttering, allele drop-out or null alleles. Results showed the presence of null alleles for almost each marker, as well as stuttering and large allele drop-out for some of them (results not shown), but no genotyping errors from allele dropouts or stuttering affected allele scoring. Additionally, a deficit of heterozygosity for almost each marker was observed. Consequently, all samples were re-amplified for each of the 10 loci under relaxed annealing temperature (4 °C lower) to verify that additional alleles had not been missed. After scoring all individuals, the same analyses were re-conducted for each marker. Furthermore, primers were redesigned for the loci where sufficient and adequate flanking regions were available. Primer sequences for this new primers are; D.tru2 F 5'FAMACCACCAATTCTCCTACGG3', D.tru2 R 5'ATGTGGGCGATGATTTCCT3'; D.tru8 F 5'PET TAAAATTGCCATGCGTCAG3', D.tru8 R 5'GAAATATATTGCAGGCTGGTAGG3'; D.tru14 F 5'VIC TTTTTGTTCTTCTGAATAGTGCAA3', D.tru14 R 5'TTTACAGCGTGTCGCCATCT3'; D.tru28 F 5'FAM CTGAGAAGTAATAAAACGTGAATTGTTG3', D.tru28 R 5'CAACAGAGCCACTGATGACAA3'; D.tru32 F 5'NED TGGGTCCTGGAGGGTAAAAT3', D.tru32 R 5'TTAGTCAGAACCTCAACTCTCAAA3'.

The total number of alleles (N<sub>A</sub), the observed (H<sub>O</sub>) and expected (H<sub>E</sub>) heterozygosity, the polymorphism information content<sup>61</sup> and the null-allele frequency (NAF) were also determined for each locus and each sampling site using the program CERVUS 3.0<sup>33</sup>. Effective population size for all the populations were also estimated using the linkage disequilibrium method<sup>62</sup> as implemented in NeEstimator (version 2.0<sup>58</sup>).

**Analyses of population genetic structure.** In order to assess the significance of the level of genetic differentiation, an analysis of molecular variance (AMOVA) was implemented in Arlequin (version 3.5.2.1)⁶³. Specifically, we evaluated the amount of variance attributable to groups of sampling sites (Gulf of Cádiz: Isla Canela and Doñana; Alboran Sea: Caleta de Vélez and Cabo de Gata; Northwestern Mediterranean coast: Gandía,
San Carles de la Ràpita and Roses), sampling sites within the groups and within sampling sites. The extent of genetic differentiation among the group of sampling locations was assessed using pairwise $F_{ST}$ estimates⁶⁴ obtained in FSTAT version 2.9.3⁶⁵.

The individual admixture proportions (*q*-values) for each individual in each sampling site, considering each one like a population, was performed using STRUCTURE v.2.2⁶⁶. Because of the presence of null alleles, analyses were done incorporating RECESSIVEALLELES option. In this case, the program assumes that the recessive allele is never observed in homozygous state but it might be present (e.g. when there might be null alleles). Moreover, the genetic differentiation of population being weak, the LOCPRIOR model has been used. The program assumes an initially probability that each individual is a resident of its sampling locality calculating posterior probabilities that individuals belong to their sampled locality. The number of clusters for each analysis was therefore assessed using the ad hoc statistic $\Delta K$⁶⁷. An admixture model with correlated allele frequency was used for each analysis with 250,000 steps of the Markov-Chain preceded by a burnin-period of 100,000 steps.

BAYESASS version 3.0⁶⁸ was used to estimate direction and rates of migration that has occurred more recently (i.e., within the last several generations) between the different groups of populations (Gulf of Cadiz, Alboran Sea and Northwestern Mediterranean coast). Analysis were performed pooling locations as recommended by Meirmans⁶⁹ (i.e. few populations with many individuals). This method allows for deviations from HWE but assumes linkage disequilibrium and constant migration rates for two generations prior to sampling. Migrants are defined by hybrid genotypes and the program determines a migration rate *m*, which is the proportion of first generation migrants. Parameters were as follows: number of iterations was 3,000,000, of which $10^6$ are burn in, and the sampling frequency was 2,000.

## References

1. Farris, D. W. et al. Fracturing of the Panamanian Isthmus during initial collision with South America. *Geology* **39**, 1007–1010, doi: 10.1130/g32237.1 (2011).

2. O'Dea, A. et al. Environmental change preceded Caribbean extinction by 2 million years. *PNAS* **104**, 5501–5506, doi: 10.1073/ pnas.0610947104 (2007).

3. Knowlton, N. & Weigt, L. A. New dates and new rates for divergence across the Isthmus of Panama. *Philos. Trans. R. Soc. B-Biol. Sci.* **265**, 2257–2263, doi: 10.1098/rspb.1998.0568 (1998).

4. Rögl, F. Palaeogeographic Considerations for Mediterranean and Paratethys Seaways (Oligocene to Miocene). *Annalen des Naturhistorischen Museums in Wien. Serie A, Für Mineralogie und Petrographie, Geologie und Paläontologie, Anthropologie und Prähistorie* **99**, 279–310 (1997).

5. Cowman, P. F., Bellwood, D. R. & Rocha, L. A. The historical biogeography of coral reef fishes: global patterns of origination and dispersal. *J. Biogeogr.* **40**, 209–224, doi: 10.1111/jbi.12003 (2013).

6. Lessios, H. A., Kessing, B. D. & Robertson, D. R. Massive gene flow across the world's most potent marine biogeographic barrier. *Proc. R. Soc. Lond. Ser. B-Biol. Sci.* **265**, 583–588, doi: 10.1098/rspb.1998.0334 (1998).

7. Lessios, H. A. & Robertson, D. R. Crossing the impassable: genetic connections in 20 reef fishes across the eastern Pacific barrier. *Proc. R. Soc. Lond. Ser. B-Biol. Sci.* **273**, 2201–2208, doi: 10.1098/rspb.2006.3543 (2006).

8. Gibbons, M. J., Yuecher, E., Thibault-Botha, D. & Helm, R. R. Patterns in marine hydrozoan richness and biogeography around southern Africa: implications of life cycle strategy. *J. Biogeogr.* **37**, 606–616, doi: 10.1111/j.1365-2699.2009.02237.x (2010).

9. Quesada, H., Zapata, C. & Alvarez, G. A multilocus allozyme discontinuity in the mussel *Mytilus galloprovincialis* - the interaction of ecological and life-history factors. *Mar. Ecol.-Prog. Ser.* **116**, 99–115, doi: 10.3354/meps116099 (1995).

10. Luttikhuizen, P. C., Campos, J., Bleijswijk, J. v., Peijnenburg, K. T. C. A. & van der Veer, H. W. Phylogeography of the common shrimp, *Crangon crangon* (L.) across its distribution range. *Mol. Phylogenet. Evol.* **46**, 1015–1030, doi: http://dx.doi.org/10.1016/j. ympev.2007.11.011 (2008).

11. Galarza, J. A., Turner, G. F., Macpherson, E. & Rico, C. Patterns of genetic differentiation between two co-occurring demersal species: the red mullet (*Mullus barbatus*) and the striped red mullet (*Mullus surmuletus*). *Can. J. Fish. Aqua. Sci.* **66**, 1478–1490, doi: 10.1139/f09-098 (2009).

12. Gómez-Díaz, E., González-Solís, J., Peinado, M. A. & Page, R. D. M. Phylogeography of the *Calonectris* shearwaters using molecular and morphometric data. *Mol. Phylogenet. Evol.* **41**, 322–332, doi: http://dx.doi.org/10.1016/j.ympev.2006.05.006 (2006).

13. Jaramillo-Correa, J. P. et al. The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey. *Mol. Ecol.* **19**, 5452–5468, doi: 10.1111/j.1365-294X.2010.04912.x (2010).

14. Mata, V. A., da Silva, L. P., Lopes, R. J. & Drovetski, S. V. The Strait of Gibraltar poses an effective barrier to host-specialised but not to host-generalised lineages of avian *Haemosporidia*. *Int. J. Parasit.* **45**, 711–719, doi: http://dx.doi.org/10.1016/j.ijpara.2015.04.006 (2015).

15. Patarnello, T., Volckaert, F. A. & Castilho, R. Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break? *Mol. Ecol.* **16**, 4426–4444, doi: 10.1111/j.1365-294X.2007.03477.x (2007).

16. Tintore, J., La Violette, P. E., Blade, I. & Cruzado, A. A Study of an Intense Density Front in the Eastern Alboran Sea: The Almeria-Oran Front. *J. Phys. Oceanogr.* **18**, 1384–1397, doi: 10.1175/1520-0485(1988)018<1384:ASOAID>2.0.CO;2 (1988).

17. Lowe, C. D., Martin, L. E., Montagnes, D. J. S. & Watts, P. C. A legacy of contrasting spatial genetic structure on either side of the Atlantic-Mediterranean transition zone in a marine protist. *PNAS* **109**, 20998–21003 (2012).

18. Baus, E., Darrock, D. J. & Bruford, M. W. Gene-flow patterns in Atlantic and Mediterranean populations of the Lusitanian sea star *Asterina gibbosa*. *Mol. Ecol.* **14**, 3373–3382, doi: 10.1111/j.1365-294X.2005.02681.x (2005).

19. Bryden, H. L., Candela, J. & Kinder, T. H. Exchange through the Strait of Gibraltar. *Prog. Oceanogr.* **33**, 201–248, doi: http://dx.doi.org/10.1016/0079-6611(94)90028-0 (1994).

20. Bayed, A. & Guillou, J. A contribution to the study of populations of the *Donax* genus - The population of *Donax trunculus* L. (Mollusca, Bivalvia). *Annales De L Institut Oceanographique* **61**, 139–147 (1985).

21. Tebble, N. British bivalve shells. *A handbook for identification* (1966).

22. FAO-FIGIS. In *Fisheries Global Information System* (http://www.fao.org/fishery/statistics/global-capture-production/query/en/, 2013).

23. Gaspar, M. B., Ferreira, R. & Monteiro, C. C. Growth and reproductive cycle of *Donax trunculus* L., (Mollusca: Bivalvia) off Faro, southern Portugal. *Fish Res.* **41**, 309–316, doi: 10.1016/s0165-7836(99)00017-x (1999).

24. Gaspar, M. B. et al. Depth segregation phenomenon in *Donax trunculus* (Bivalvia: Donacidae) populations of the Algarve coast (southern Portugal). *Sci. Mar.* **66**, 111–121, doi: 10.3989/scimar.2002.66n2111 (2002).

25. Nantón, A., Freire, R., Arias-Pérez, A., Gaspar, M. B. & Méndez, J. Identification of four *Donax* species by PCR-RFLP analysis of cytochrome c oxidase subunit I (COI). *Eur. Food Res. Technol.* **240**, 1129–1133, doi: 10.1007/s00217-015-2416-z (2015).
26. Nantón, A., Arias-Pérez, A., Méndez, J. & Freire, R. Characterization of nineteen microsatellite markers and development of multiplex PCRs for the wedge clam *Donax trunculus* (Mollusca: Bivalvia). *Mol. Biol. Rep.* **41**, 5351–5357, doi: 10.1007/s11033-014-3406-0 (2014).

27. Plohl, M. & Cornudella, L. Characterization of Interrelated Sequence Motifs in Four Satellite DNAs and Their Distribution in the Genome of the Mollusc *Donax trunculus*. *J. Mol. Evol.* **44**, 189–198, doi: 10.1007/pl00006135 (1997).

28. Theologidis, I., Fodelianakis, S., Gaspar, M. B. & Zouros, E. Doubly uniparental inheritance (DUI) of mitochondrial dna in *Donax trunculus* (Bivalvia: Donacidae) and the problem of its sporadic detection in Bivalvia. *Evolution* **62**, 959–970, doi: 10.1111/j.1558-5646.2008.00329.x (2008).

29. Hauser, L., Waples, R. S. & Carvalho, G. R. Special Issue: Advances in Marine Fish and Fisheries Genetics. *Fish. Fish.* **9**, 331–332, doi: 10.1111/j.1467-2979.2008.00307.x (2008).

30. Waples, R. S. & Gaggiotti, O. What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. *Mol Ecol* **15**, 1419–1439, doi: 10.1111/j.1365-294X.2006.02890.x (2006).

31. Funk, W. C., McKay, J. K., Hohenlohe, P. A. & Allendorf, F. W. Harnessing genomics for delineating conservation units. *Trends Ecol Evol* **27**, 489–496, doi: 10.1016/j.tree.2012.05.012 (2012).

32. van Oosterhout, C., Hutchinson, W. F., Wills, P. M. & Shipley, P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. *Mol. Ecol. Notes* **4**, 535–538, doi: 10.1111/j.1471-8286.2004.00684.x (2004).

33. Kalinowski, S. T., Taper, M. L. & Marshall, T. C. Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment. *Mol. Ecol.* **16**, 1099–1106, doi: 10.1111/j.1365-294X.2007.03089.x (2007).

34. Riesgo, A. et al. Comparative description of ten transcriptomes of newly sequenced invertebrates and efficiency estimation of genomic sampling in non-model taxa. *Front. Zool.* **9**, 33, doi: 10.1186/1742-9994-9-33 (2012).

35. García-Cisneros, A., Palacín, C., Ben Khadra, Y. & Pérez-Portela, R. Low genetic diversity and recent demographic expansion in the red starfish *Echinaster sepositus* (Retzius 1816). *Sci. Rep.* **6**, doi: 10.1038/srep33269 (2016).

36. Wangensteen, O. S., Turon, X., Pérez-Portela, R. & Palacín, C. Natural or Naturalized? Phylogeography Suggests That the Abundant Sea Urchin *Arbacia lixula* Is a Recent Colonizer of the Mediterranean. *PLoS One* **7**, doi: 10.1371/journal.pone.0045067 (2012).

37. Bowen, B. W. et al. Comparative phylogeography of the ocean planet. *PNAS* **113**, 7962–7969, doi: 10.1073/pnas.1602404113 (2016).

38. Hellberg, M. E. Gene Flow and Isolation among Populations of Marine Animals. *Annu. Rev. Ecol. Evol. Syst.* **40**, 291–310, doi: 10.1146/annurev.ecolsys.110308.120223 (2009).

39. Souche, E. L. et al. Range-wide population structure of European sea bass *Dicentrarchus labrax*. *Biol. J. Linnean Soc.* **116**, 86–105, doi: 10.1111/bij.12572 (2015).

40. Palero, F., Abelló, P., Macpherson, E., Beaumont, M. & Pascual, M. Effect of oceanographic barriers and overfishing on the population genetic structure of the European spiny lobster (*Palinurus elephas*). *Biol. J. Linnean Soc.* **104**, 407–418, doi: 10.1111/j.1095-8312.2011.01728.x (2011).

41. Glynn, F. et al. High-resolution genetic analysis reveals extensive gene flow within the jellyfish *Pelayia noctiluca* (Scyphozoa) in the North Atlantic and Mediterranean Sea. *Biol. J. Linnean Soc.* **117**, 252–263, doi: 10.1111/bij.12654 (2016).

42. García Lafuente, J. & Ruiz, J. The Gulf of Cádiz pelagic ecosystem: A review. *Prog. Oceanogr.* **74**, 228–251, doi: 10.1016/j. pocean.2007.04.001 (2007).

43. Sala-Bozano, M., Ketmaier, V. & Mariani, S. Contrasting signals from multiple markers illuminate population connectivity in a marine fish. *Mol. Ecol.* **18**, 4811–4826, doi: 10.1111/j.1365-294X.2009.04404.x (2009).

44. Launey, S., Ledu, C., Boudry, P., Bonhomme, F. & Naciri-Graven, Y. Geographic Structure in the European Flat Oyster (*Ostrea edulis* L.) as Revealed by Microsatellite Polymorphism. *J. Hered.* **93**, 331–351, doi: 10.1093/jhered/93.5.331 (2002).

45. Chiesa, S. et al. Null alleles of microsatellites for Manila clam *Ruditapes philippinarum*. *Anim. Genet.* **47**, 135–136, doi: 10.1111/ age.12382 (2016).

46. Kim, E. M. et al. New polymorphic microsatellite markers for the Korean manila clam (*Ruditapes philippinarum*) and their application to wild populations. *Genet. Mol. Res.* **13**, 8163–8173, doi: 10.4238/2014.October.7.11 (2014).

47. Hargrove, J. S., Sturmer, L., Scarpa, J. & Austin, J. D. Assessment of Genetic Diversity in Wild and Aquaculture Stocks of *Mercenaria mercenariensis* Florida. *J. Shellfish Res.* **34**, 355–365, doi: 10.2983/035.034.0218 (2015).

48. Zouros, E. & Foltz, D. W. Possible explanations of the heterozygote deficiency in Bivalve mollusks. *Malacologia* **25**, 583–591 (1984).

49. Astanei, I., Gosling, E., Wilson, J. I. M. & Powell, E. Genetic variability and phylogeography of the invasive zebra mussel, *Dreissena polymorpha* (Pallas). *Mol. Ecol.* **14**, 1655–1666, doi: 10.1111/j.1365-294X.2005.02530.x (2005).

50. Rico, C., Rico, I. & Hewitt, G. 470 Million Years of Conservation of Microsatellite Loci among Fish Species. *Proc. R. Soc. Lond. Ser. B-Biol. Sci.* **263**, 549–557, doi: 10.1098/rspb.1996.0083 (1996).

51. Hoarau, G. et al. Low effective population size and evidence for inbreeding in an overexploited flatfish, plaice (*Pleuronectes platessa* L.). *Proc. R. Soc. Lond. Ser. B-Biol. Sci.* **272**, 497–503, doi: 10.1098/rspb.2004.2963 (2005).

52. Ma, H. T., Jiang, H. B., Liu, X. Q., Wu, X. P. & Wei, X. M. Polymorphic microsatellite loci for the razor clam, *Sinonovacula constricta*. *Genet. Mol. Res.* **14**, 145–148, doi: 10.4238/2015.January.15.17 (2015).

53. Li, H. J., Zhang, J. J., Li, H. B., Gao, X. G. & He, C. B. Development and characterization of 13 polymorphic microsatellite markers for the Chinese surf clam (*Mactra chinensis*) through Illumina paired-end sequencing. *Conserv. Genet. Resour.* **6**, 877–879, doi: 10.1007/s12686-014-0229-1 (2014).

54. Pereira, S., Arias, A., Méndez, J., Insua, A. & Freire, R. Isolation of twelve microsatellite markers in the pullet carpet shell *Venerupis pullastra* (Bivalvia: Veneridae). *Conserv. Genet. Resour.* **2**, 201–203, doi: 10.1007/s12686-009-9134-4 (2010).

55. Jamieson, I. G. & Allendorf, F. W. How does the 50/500 rule apply to MVPs? *Trends Ecol. Evol.* **27**, 578–584, doi: 10.1016/j. tree.2012.07.001 (2012).

56. Frankham, R., Bradshaw, C. J. A. & Brook, B. W. Genetics in conservation management: Revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. *Biol. Conserv.* **170**, 56–63, doi: 10.1016/j.biocen.2013.12.036 (2014).

57. Rico, C., Drake, P., Macpherson, E. & Cuesta, J. A. Estimación de la diversidad genética y del tamaño efectivo de la población de coquina *Donax trunculus* del Parque Nacional de Doñana y su contribución a áreas no protegidas., (Internal Report. Organismo Autónomo Parques Nacionales, Ministerio de Agricultura, Alimentación y Medio Ambiente, Madrid Spain, 2015).

58. Do, C. et al. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. *Mol. Ecol. Resour.* **14**, 209–214, doi: 10.1111/1755-0998.12157 (2014).

59. Aljanabi, S. M. & Martinez, I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. *Nucleic Acids Res.* **25**, 4692–4693, doi: 10.1093/nar/25.22.4692 (1997).

60. Rousset, F. *Genepop*007: a complete re-implementation of the genepop software for Windows and Linux. *Mol. Ecol. Res.* **8**, 103–106, doi: 10.1111/j.1471-8286.2007.01931.x (2008).

61. Huang, D. W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nat. Protocols* **4**, 44–57, doi: 10.1038/nprot.2008.211 (2008).

62. Waples, R. S. & Do, C. H. I. Ldne: a program for estimating effective population size from data on linkage disequilibrium. *Mol. Ecol.* **8**, 753–756, doi: 10.1111/j.1755-0998.2007.02061.x (2008).

63. Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. *Mol. Ecol. Res.* **10**, 564–567, doi: 10.1111/j.1755-0998.2010.02847.x (2010).

64. Weir, B. S. & Cockerham, C. C. Estimating F-Statistics for the Analysis of Population Structure. *Evolution* **38**, 1358–1370, doi: 10.2307/2408641 (1984).
65. Goudet, J. F. S. T. A. T. (Version 1.2): A Computer Program to Calculate F-Statistics. *J. Hered.* **86**, 485-486 (1995).

66. Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of Population Structure Using Multilocus Genotype Data. *Genetics* **155**, 945-959, doi: 10.3410/f.1015548.197423 (2000).

67. Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software structure: a simulation study. *Mol. Ecol.* **14**, 2611-2620, doi: 10.1111/j.1365-294X.2005.02553.x (2005).

68. Wilson, G. A. & Rannala, B. Bayesian Inference of Recent Migration Rates Using Multilocus Genotypes. *Genetics* **163**, 1177-1191, doi: 10.3410/f.1012487.188326 (2003).

69. Meirmans, P. G. Nonconvergence in Bayesian estimation of migration rates. *Mol. Ecol. Res.* **14**, 726-733, doi: 10.1111/1755-0998.12216 (2014).

## Acknowledgements

C.R. gratefully acknowledges the Spanish Ministry of Education (Grant N° PR2010-0601) and the Spanish Scientific Research Council for the financial support received during his stay at L.B. laboratory. C.R. also acknowledges F. Hiraldo Cano, J.J. Negro Balmaseda, Eric Normandeau and Guillaume Côté for their support during his stay in Canada. This research was financially supported by Organismo Autónomo de Parques Nacionales Grant OAPN122_2010 to C.R., E.M., J.A.C. and P.D. We are also grateful for the support given, throughout the project, by the Coordination Office of Doñana National Park and the Consejería de Agricultura, Pesca y Desarrollo Rural of the Andalusian Government. We also thank Dr Serrão who handled the MS as EBM and two anonymous referees for their constructive comments on a previous version.

## Author Contributions

C.R., E.M., J.A.C. and P.D. designed the project. P.D. and J.C. organised and sampled all specimens. C.R. conducted the laboratory work, the microsatellite development and the first round of genotyping. E.K. conducted the second round of genotyping and all the genotyping of the second round of amplifications. C.L. contributed with microsatellite experimental design and analyses. A.D.M. analysed the results. A.D.M and C.R. wrote the manuscript. L.B. provided advice in the experimental design, execution of the project and laboratory facilities. All authors reviewed and contribute to the writing of the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Marie, A. D. *et al*. Implications for management and conservation of the population genetic structure of the wedge clam *Donax trunculus* across two biogeographic boundaries. *Sci. Rep.* **6**, 39152; doi: 10.1038/srep39152 (2016).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2016