# Absence of spatial genetic structure in common dentex (Dentex dentex Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers

Alexiane Viret, Dimitris Tsaparis, Costas Tsigenopoulos, Patrick Berrebi, Andrea Sabatini, Marco Arculeo, Chiheb Fassatoui, Antonios Magoulas, Michel Marengo, Beatriz Morales-Nin, et al.

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

## ► To cite this version:

Alexiane Viret, Dimitris Tsaparis, Costas Tsigenopoulos, Patrick Berrebi, Andrea Sabatini, et al. Absence of spatial genetic structure in common dentex (Dentex dentex Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers. PLoS ONE, 2018, 13 (9), pp.e0203866. 10.1371/journal.pone.0203866. hal-03117680

HAL Id: hal-03117680

https://hal.umontpellier.fr/hal-03117680

Submitted on 21 Jan 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
RESEARCH ARTICLE

# Absence of spatial genetic structure in common dentex (*Dentex dentex* Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers

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

Alexiane Viret¹, Dimitris Tsaparis², Costas S. Tsigenopoulos², Patrick Berrebi¹,
Andrea Sabatini³, Marco Arculeo⁴, Chiheb Fassatoui⁵, Antonios Magoulas²,
Michel Marengo⁶,⁷, Beatriz Morales-Nin⁸, Nathalie Caill-Milly⁹, Eric D. H. Durieux⁶,⁷*

1 ISEM, Université de Montpellier, CNRS, IRD, EPHE, Montpellier, France, 2 Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research (HCMR), Thalassocosmos, Heraklion, Crete, Greece, 3 University of Cagliari, Biologia Animale ed Ecologia, Cagliari, Italy, 4 Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche, Università di Palermo, Palermo, Italy, 5 Institut National Agronomique de Tunisie, UR03AGRO1 Écosystèmes & Ressources Aquatiques, Tunisia, 6 UMR 6134 CNRS-UCPP Sciences pour l'Environnement, Université de Corse Pasquale Paoli, Corte, France, 7 UMS 3514, CNRS-UCPP Plateforme marine Stella Mare, Université de Corse Pasquale Paoli, Biguglia, France, 8 Department of Natural Resources, Mediterranean Institute for Advanced Studies (IMEDEA CSIC-UIB), Esporlas, Balearic Islands, Spain, 9 Ifremer, Laboratory Environnement Resources of Arcachon, UFR Sciences and Technics, 1 allée du Parc Montaury, Anglet, France

* durieux@univ-corse.fr

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

OPEN ACCESS

**Citation:** Viret A, Tsaparis D, Tsigenopoulos CS, Berrebi P, Sabatini A, Arculeo M, et al. (2018) Absence of spatial genetic structure in common dentex (*Dentex dentex* Linnaeus, 1758) in the Mediterranean Sea as evidenced by nuclear and mitochondrial molecular markers. PLoS ONE 13 (9): e0203866. https://doi.org/10.1371/journal.pone.0203866

**Editor:** Tzen-Yuh Chiang, National Cheng Kung University, TAIWAN

**Received:** December 4, 2017

**Accepted:** August 29, 2018

**Published:** September 12, 2018

**Copyright:** © 2018 Viret et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** MtDNA sequences are all submitted to Genbank and accession numbers given in the text. All sequences (COI and D-loop) and microsatellite genotypes are provided as Supporting Information.

**Funding:** This work was supported by Région Corse, Union Européenne.

**Competing interests:** The authors have declared that no competing interests exist.

## Abstract

The common dentex, *Dentex dentex*, is a fish species which inhabits marine environments in the Mediterranean and Northeast Atlantic regions. This is an important species from an ecological, economic and conservation perspective, however critical information on its population genetic structure is lacking. Most samples were obtained from the Mediterranean Sea (17 sites) with an emphasis around Corsica (5 sites), plus one Atlantic Ocean site. This provided an opportunity to examine genetic structuring at local and broader scales to provide science based data for the management of fishing stocks in the region. Two mitochondrial regions were examined (D-loop and COI) along with eight microsatellite loci. The COI data was combined with publicly available sequences and demonstrated past misidentification of common dentex. All markers indicated the absence of population genetic structure from the Bay of Biscay to the eastern Mediterranean Sea. Bayesian approaches, as well as the statistical tests performed on the allelic frequencies from microsatellite loci, indicated low differentiation between samples; there was only a slight (p = 0.05) indication of isolation by distance. Common dentex is a marine fish species with a unique panmictic population in the Mediterranean and likely in the Atlantic Ocean as well.
## Introduction

Patterns of geographic genetic structure in marine fishes rely mainly on life histories and environmental drivers over space and time [1, 2, 3]. Compared to terrestrial organisms, marine species, and marine fishes in particular, generally show a high genetic diversity and a weak spatial genetic differentiation [1, 2, 4, 5, 6]. This is mainly attributed the high prevalence of external fertilization and the production of a huge amount of larvae with high dispersive capability [7] associated with the lack of barriers in marine waters which in turn should facilitate high gene flow between remote populations [6]. However, the development and use of polymorphic genetic markers allow for the describing of a certain degree of spatial genetic structure in marine fishes due to specific connectivity in both adults' movement and early life stage dispersal [8, 9, 10]. Knowledge of the genetic structure of exploited populations is therefore essential to conservation and sustainable management of fishery resources [8, 9, 10, 11]. This type of research is especially necessary to identify populations, delimitate stocks, and determine the connectivity and potential resilience within populations.

The Sparidae family (Rafinesque, 1810) belongs to teleost fishes comprising more than 118 species distributed among 35 genera found in all tropical and temperate oceans [12, 13]. These species are of significant economic interest since the majority are exploited by commercial fishing activities; some have also become important for the aquaculture industry, such as the very popular gilthead sea bream (*Sparus aurata*) in the Mediterranean Sea and the Japanese red sea bream (*Pagrus major*) in East Asia. The common dentex, *Dentex dentex* (Linnaeus, 1758), is a sparid living along the Mediterranean, and the Atlantic Ocean coasts from the British Islands to Senegal, and occasionally in the Black Sea [14]. The species lives near the bottom of sea, from a few meters to 200 m depth, preferably on a rocky substrate. Adults are usually solitary and congregate only in the spring for reproduction [15]. In Corsica, fishermen observed that they congregate between 40 and 100 m depth on hard substrate (rocks, wrecks) [14]. The common dentex is gonochoristic, reaching sexual maturity between 2 and 4 years old [16, 17, 18]. It grows to a maximum length of 100 cm and a weight of 13 kg, with a relatively long life span (up to 36 years) [19]. Due to its large size, flesh quality, and high commercial value, the species is of great interest to both artisanal and recreational fisheries [20]. Furthermore, the common dentex is classified by the International Union for the Conservation of Nature (IUCN) as "vulnerable" in the Red List of Threatened Species. However, little is known about the population genetic structure of the common dentex. Previous studies sampling partly the same populations as in this survey and based either on allozymes and partial sequences of the mitochondrial D-loop (control) region [21] or morphological biometrics [22] have proposed a strong divergence between the Atlantic and Mediterranean populations. In fact, the Strait of Gibraltar is considered a biogeographical breakpoint [23] between the Mediterranean Sea and the Atlantic Ocean both for pelagic/epipelagic and demersal fish species. Species such as the European sea bass *Dicentrarchus labrax*, the common dolphinfish *Coryphaena hippurus* and other sparids like the striped *Lithognathus mormyrus*, and the black sea bream *Spondyliosoma cantharus* show significant genetic differentiation in this area (reflected by high Fst values between populations both in allozymes and mtDNA) [23, 24].

The aim of the present study was to assess the species identification and then the genetic structure and phylogeography of the common dentex by combining a representative population sampling scheme in the Mediterranean Sea and multiple markers. A denser population sampling scheme includes samples already analyzed by Bargelloni et al. [21] and adds new samples covering the Mediterranean Sea from the east (N. Aegean Sea, Greece), to the west (Alicante, Spain), and up to the Bay of Biscay in the Atlantic Ocean (France). First, sequences from the mitochondrial DNA (mtDNA) cytochrome oxidase (COI) locus were used for
barcoding purposes to check for any taxonomic inconsistencies when an individual's species is determined; subsequently, mtDNA D-loop sequences and multi-locus microsatellite genotype data were employed to refine the intra-specific common dentex population genetic structure within the study area. A recent study [25] established that there is no genetic structure around Corsica. To a lesser extent, we aim to define the species structure in the Mediterranean Sea and one site in the Atlantic Ocean.

## Materials and methods

### Sampling

Common dentex tissues (fin clips), sampled in 2002 and kept in ethanol in IMBBC for Bargelloni et al. [21], were re-used in the present study. More recent samples (fin clips and muscle tissues) were obtained by professional fishermen, and fish were caught between the years 2012 and 2016 with longlines, gillnets, trammel nets and trolling. Sampling localities and sizes are reported in Fig 1 and Table 1. These individuals originated from 19 different fishing locations distributed across the Northeast Atlantic Ocean (France and Portugal), the Western basin of the Mediterranean Sea (Spain: two sites, France: seven sites, Sardinia, Sicily and Tunisia) and the Eastern basin of the Mediterranean Sea (Lampedusa, southern Tunisia, the Adriatic Sea, Crete and N. Aegean Sea). At a finer scale, six samples were located all around the Corsican coasts: Bonifacio in the South, Ajaccio, Galeria and Saint Florent in the West, Giraglia in the North and Bastia in the East (Fig 1).

### DNA extraction

DNA was extracted from fin clips and muscle tissues. For each sample, a piece of tissue about 2 mm x 4 mm was taken from ethanol-preserved specimens, dried and dipped into 200 µL of Chelex previously heated to 56°C. The DNA extraction was performed according to a Chelex protocol [26] with proteinase K with a final concentration of 0.198 mg. mL-1. The mix was placed overnight in an oven at 56°C and then at 95°C for 5 min. Finally, samples were centrifuged at 3,000 rotations per min for 2 min. The supernatant was diluted 1:20 with sterile water. The samples were then frozen (-20°C) until use.

### MtDNA sequencing analyses

In order to resolve any taxonomic uncertainties concerning either old samples [21] or new samples, we performed a DNA barcoding analysis using the cytochrome c oxidase subunit 1 (COI) mitochondrial gene as a marker for species-level identification [27]. A 612 bp fragment at the 5' end of COI was amplified with primers FishCoxI-F2 (5' -TCGACTAATCATAAAGATATCGG CAC-3') and FishCoxI-R2 (5' -ACTTCAGGGTGACCGAAGAATCAGAA-3') [28]. The 10 µL PCR mix included 20-50 ng of template DNA, 1x Taq buffer, 0.2 µM of each primer, 0.2 mM dNTP mix, 1 U of Taq polymerase, 2.5 mM MgCl2 and ultrapure water. The PCR cycling protocol consisted of an initial step of 2 min at 95°C, followed by 35 cycles of 0.5 min at 95°C, 0.5 min at 50°C, and 1 min at 72°C, followed by a final extension step at 72°C for 10 min.

A 444-456 bp long fragment of the D-loop region was also amplified in order to study the mtDNA variability and geographic structure at the species level. Primers tRNA-Pro-L (5' -AC CATTGGCTCCCAAAGCTA-3') and H16498 (5' -CCTGAAGTAGGAACCAGATG-3') [29] were used for this amplification. PCR mix concentrations and cycling protocol were slightly different from COI with an MgCl2 concentration of 3 mM and an annealing temperature of 49°C.

Purification of PCR products in both cases was performed according to a standard ethanol precipitation protocol. Sequencing reactions were carried out using the BigDye Terminator
<figure><img src="image_5.png" /><figcaption>Fig 1. Map of sampling localities of common dentex (*Dentex dentex*) with respective years (blue circle: 2002 Bargelloni [21] samples; yellow circle: 2012–2014 samples; green circle: 2016 sample).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203866.g001

version 3.1 Cycle Sequencing Kit (Applied Biosystems, Inc.). Products were sequenced in both directions on an ABI 3730 capillary sequencer following the manufacturer's instructions. Individual sequences were edited with MEGA 6.06 [30], re-examined manually by visual inspection of raw Fluorigram data and then aligned with ClustalW (as implemented in MEGA 6.06).

## MtDNA sequence diversity and phylogeographic analyses

COI sequences from our samples were aligned to those from *Dentex* species and other closely related Sparidae species available in GenBank, using MEGA 6.06 [30]. All haplotypes from the final multiple alignment were used for the construction of a haplotype network with the median-joining network method [31] and default settings as implemented in the program Network 5.0.0.0 (Fluxus Technology Ltd., http://www.fluxus-engineering.com/). Networks are
<table><caption>Table 1. Collection date and number of analyzed samples (N) for each of the 19 localities and for each of the three markers used.</caption><thead><tr><th>Locality</th><th>Date</th><th>N microsatellites</th><th>N COI</th><th>N D-loop</th></tr></thead><tbody><tr><td>Heraklion (North Crete, Greece)</td><td>2002</td><td>24</td><td>5</td><td>5</td></tr><tr><td>North Aegean Sea (Greece)</td><td>2002</td><td>13</td><td>5</td><td>6</td></tr><tr><td>Otranto (South Adriatic, Italy)</td><td>2002</td><td>23</td><td>5</td><td>5</td></tr><tr><td>Sicily (Italy)</td><td>2014</td><td>13</td><td>-</td><td>-</td></tr><tr><td>Lampedusa (Italy)</td><td>2013</td><td>13</td><td>5</td><td>5</td></tr><tr><td>North Tunisia</td><td>2013</td><td>14</td><td>2</td><td>2</td></tr><tr><td>South Tunisia</td><td>2013</td><td>14</td><td>1</td><td>-</td></tr><tr><td>Sardinia (Italy)</td><td>2013</td><td>23</td><td>6</td><td>5</td></tr><tr><td>Giraglia (Corsica, France)</td><td>2012–2013</td><td>25</td><td>5</td><td>5</td></tr><tr><td>Bastia (Corsica, France)</td><td>2012–2013</td><td>-</td><td>6</td><td>5</td></tr><tr><td>St Florent (Corsica, France)</td><td>2012</td><td>18</td><td>5</td><td>3</td></tr><tr><td>Galeria (Corsica, France)</td><td>2012–2013</td><td>31</td><td>4</td><td>3</td></tr><tr><td>Ajaccio (Corsica, France)</td><td>2012–2013</td><td>15</td><td>5</td><td>5</td></tr><tr><td>Bonifacio (Corsica, France)</td><td>2012</td><td>24</td><td>5</td><td>5</td></tr><tr><td>St Raphael (France)</td><td>2013–2014</td><td>14</td><td>8</td><td>6</td></tr><tr><td>Baleares (Spain)</td><td>2013</td><td>19</td><td>7</td><td>5</td></tr><tr><td>Alicante (Spain)</td><td>2002</td><td>21</td><td>5</td><td>6</td></tr><tr><td>Faro (Portugal)</td><td>2002</td><td>24</td><td>20</td><td>-</td></tr><tr><td>Bay of Biscay (France)</td><td>2016</td><td>8</td><td>8</td><td>8</td></tr><tr><td>SUM</td><td></td><td>336</td><td>107</td><td>79</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0203866.t001

very helpful for the inferring of haplotype relationships and are used here in order to graphically represent species delimitation. A Maximum Likelihood phylogenetic tree was also constructed with the online software PHYML 3 [32] as implemented in the platform ATGC (http://www.atgc-montpellier.fr/). Model selection was made automatically by the program with the Smart Model Selection application [33], while branching support was surveyed with the aBayes method [34] and bootstrapping (10000 bootstraps). Basic diversity indices as the number of haplotypes (nh), nucleotide diversity (π), and haplotype (gene) diversity (h) for each sampled population were estimated using DnaSP 5 [35].

Using D-loop sequences, a hierarchical analysis of molecular variance (AMOVA) implemented in ARLEQUIN 3.5 was used to check for hypothesized patterns of spatial genetic structure assessing variance components between areas, between populations within areas and among individuals within populations [36]. Phylogeographic history of *D. dentex* was investigated using NETWORK 5.0.0.0 with default settings to reconstruct a median-joining network of haplotypes [31]. In order to explore demographic history, a mismatch distribution analysis in the whole data set of D-loop sequences (sum of sampled populations) was performed with ARLEQUIN 3.5. The simulated distribution of pairwise nucleotide differences (under the sudden expansion model and 10000 bootstrap replications) was compared with the observed distribution using the sum of square deviations (SSD) and the raggedness index (r) [36, 37].

## Microsatellite marker analysis

Microsatellite analyses were limited to the *D. dentex* samples alone, determined according to the COI marker. Eight microsatellite hypervariable markers established from other spaid species were applied to *D. dentex* (S1 Table data provides all technical information). For each pair of markers, one of the 5' ends of the two primers was end-labelled with a fluorescent dye, either 6-FAM, HEX or NED. Polymerase chain reactions (PCR) were performed using the Qiagen
multiplex PCR kit (Qiagen, Courtaboeuf, France) in a final volume of 10 µL containing 3 µL of genomic DNA diluted at 10ng/µL, 5 µL of Qiagen PCR Master Mix, 1 µL of Qiagen Q-solution, and 1 µL of primers mix at 2 µM each (Eurofins MWG Operon, Ebersberg, Germany). Amplifications were carried out in a GeneAmp PCR System 2700 thermal cycler (Applied Biosystems) according to the supplier's instructions (Qiagen multiplex PCR kit): initial denaturation at 95°C for 15 min; followed by 35 cycles of denaturation at 94°C (30 s), annealing (56 and 58°C, as indicated in S1 Table, during 90 s) and extension (72°C, 60 s); with a final extension step at 60°C for 30 min. Amplified PCR fragments were then diluted and separated on an ABIPRISM 3130/xl/ sequencer (Applied Biosystems) with GeneScan 500 Rox dye size standards. Allele sizes were determined using the GeneMapper v4.1 software system (Applied Biosystems, Life Technologies). The genotype matrix was then constructed and used as a basis for all of the following statistical analyses.

A first selection of genotypes was processed; missing data were limited to 2 loci among 8 for each individual in the matrix to reduce statistical bias. The software GENETIX [38] was used, as for most of the following tests (diversity, F statistics, IBD, FCA).

Samples' polymorphism were estimated from the expected unbiased heterozygosity (Hnb) [39], observed heterozygosity (Ho) and average number of alleles per locus (A) parameters estimate the polymorphism of each sample.

The intra-population fixation index, Fis, was assessed with the Weir and Cockerham's f esti-
mator [40], in order to check the Hardy-Weinberg equilibrium. Its significance was tested by
comparing the natural Fis to 5,000 within the sample allele's permutated matrices. The
Wright's index Fst, applied to nuclear genotypes, allowed for evaluating the differentiation
between samples. Fst was calculated through the Weir and Cockerham's θ estimator; signifi-
cance was tested for 5,000 permutations of genotypes between samples.

Isolation by distance (IBD) was tested by comparing a first matrix of genetic distances and a second matrix of geographical distances (obtained from Google Earth, tools path). Then the Mantel's test [41] was applied using the genetic distance based on Fst/(1-Fst) as recommended by Rousset [42]. The value of Z, the Mantel's coefficient, between the two matrices of distances was calculated with the true data, and then the significance of each test was assessed by comparing the true value with the series of pseudo-values produced by 5,000 permutations of the populations' order of one of the two matrices of distances.

Multivariate analyses display graphical representations of correlations between individuals
depending on different variables. The factorial correspondence analysis, FCA, was based on an
allelic disruptive matrix.

The existence and the number of sub-groups contained in the analyzed *D. dentex* sampling
can be estimated by assignment method using the Bayesian clustering program STRUCTURE
[43, 44, 45]. In our analysis, we selected the free clustering option with no prior population infor-
mation and the admixture model. For the parameter settings, we set the number of MCMC repe-
titions at 200,000 after an initial burn-in of 100,000 repetitions. Values of K were set from 1 to 17
with 10 iterations for each value of K. Last, the ΔK Evanno's method [46] was used as an aid to
define the most informative partition by calculating the variation of the rate of likelihood between
K and K+1. These values correspond to local maximums of the curve ΔK function of K. The best
K value is automatically calculated by the online program STRUCTURE HARVESTER [47].

Results

COI DNA barcoding

In total, 107 individuals from 19 locations (Table 1) were sequenced for the COI gene. Haplo-
type and nucleotide diversity values overall and for each sampling site are presented in
<table><caption>Table 2. Review of basic genetic diversity indices estimated for *Dentex dentex* populations in 16 of the 19 sampling localities and for both mitochondrial markers.</caption><thead><tr><th rowspan="2">LOCALITY</th><th colspan="4">COI 612bp</th><th colspan="4">D-loop 444-455bp</th></tr><tr><th>N</th><th>nh</th><th>h</th><th>π</th><th>N</th><th>nh</th><th>h</th><th>π</th></tr></thead><tbody><tr><td>North Aegean Sea (Greece)</td><td>5</td><td>2</td><td>0.4</td><td>0.00065</td><td>6</td><td>5</td><td>0.9</td><td>0.020089</td></tr><tr><td>Otranto (South Adriatic, Italy)</td><td>5</td><td>3</td><td>0.7</td><td>0.00131</td><td>5</td><td>4</td><td>0.9</td><td>0.021875</td></tr><tr><td>Lampedusa (Italy)</td><td>5</td><td>2</td><td>0.6</td><td>0.00078</td><td>5</td><td>5</td><td>1.0</td><td>0.017857</td></tr><tr><td>North Tunisia</td><td>2</td><td>2</td><td>1.0</td><td>0.00163</td><td>2</td><td>2</td><td>1.0</td><td>0.020089</td></tr><tr><td>South Tunisia</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Sardinia (Italy)</td><td>6</td><td>4</td><td>0.8</td><td>0.00163</td><td>5</td><td>5</td><td>1.0</td><td>0.023214</td></tr><tr><td>Giraglia (Corsica, France)</td><td>5</td><td>1</td><td>0.0</td><td>-</td><td>5</td><td>5</td><td>1.0</td><td>0.022098</td></tr><tr><td>Bastia (Corsica, France)</td><td>6</td><td>3</td><td>0.6</td><td>0.00109</td><td>5</td><td>4</td><td>0.9</td><td>0.016295</td></tr><tr><td>St Florent (Corsica, France)</td><td>5</td><td>3</td><td>0.7</td><td>0.00196</td><td>3</td><td>3</td><td>1.0</td><td>0.039474</td></tr><tr><td>Galeria (Corsica, France)</td><td>4</td><td>3</td><td>0.8</td><td>0.00245</td><td>3</td><td>3</td><td>1.0</td><td>0.017857</td></tr><tr><td>Ajaccio (Corsica, France)</td><td>5</td><td>2</td><td>0.4</td><td>0.00131</td><td>5</td><td>5</td><td>1.0</td><td>0.028571</td></tr><tr><td>Bonifacio (Corsica, France)</td><td>5</td><td>2</td><td>0.4</td><td>0.00065</td><td>5</td><td>4</td><td>0.9</td><td>0.016923</td></tr><tr><td><strong>Corsica all sites</strong></td><td><strong>30</strong></td><td><strong>7</strong></td><td><strong>0.462</strong></td><td><strong>0.00371</strong></td><td><strong>26</strong></td><td><strong>22</strong></td><td><strong>0.987</strong></td><td><strong>0.022794</strong></td></tr><tr><td>St Raphael (France)</td><td>8</td><td>4</td><td>0.6</td><td>0.00123</td><td>6</td><td>6</td><td>1.0</td><td>0.022173</td></tr><tr><td>Baleares (Spain)</td><td>7</td><td>1</td><td>0.0</td><td>-</td><td>5</td><td>5</td><td>1.0</td><td>0.021875</td></tr><tr><td>Alicante (Spain)</td><td>5</td><td>3</td><td>0.7</td><td>0.00131</td><td>6</td><td>6</td><td>1.0</td><td>0.014435</td></tr><tr><td>Bay of Biscay (France)</td><td>8</td><td>2</td><td>0.3</td><td>0.00041</td><td>8</td><td>7</td><td>1.0</td><td>0.020727</td></tr><tr><td><strong>ALL</strong></td><td><strong>87</strong></td><td><strong>14</strong></td><td><strong>0.465</strong></td><td><strong>0.00097</strong></td><td><strong>79</strong></td><td><strong>57</strong></td><td><strong>0.985</strong></td><td><strong>0.021041</strong></td></tr></tbody></table>

N = number of individuals analyzed in each sampling locality; nh = number of haplotypes analyzed, h = haplotype diversity and π = nucleotide diversity.

https://doi.org/10.1371/journal.pone.0203866.t002

Table 2. In addition, we included all the available COI sequences from *D. dentex*, as well as some sequences from other *Dentex* species (*D. gibbosus*, *D. maroccanus*, *D. macrophthalmus*, *D. angolensis*, *D. canariensis*) and other Sparidae species (*Pagrus caeruleostictus*, *Cheimerius nufar*, *Viridentex acromegalus*) already published in GenBank. A total of 66 sequences were downloaded from GenBank (S2 Table) and added to the final alignment/data set.

The construction of the COI median-joining network (Fig 2) was based on 51 different haplotypes (S2 and S4 Tables). All COI sequences of *D. dentex* published in GenBank and all the present study's sequences from 17 locations (Faro excluded) were clustered together in one major haplogroup in dark blue (Fig 2A). On the contrary, Bargelloni et al.'s [21] sampled sequences from Faro Portugal first attributed to *D. dentex* were significantly divergent and most of them (17 out of 20) constituted a distinct haplogroup (haplotypes HC21-HC27) in green. This Portuguese group, in which the published sequences for *D. gibbosus* are also included (HC21), probably corresponds to this last species. Haplotypes HC28 and HC29 were obtained from three samples of the same location [21] (Portugal) but seem to belong to a different haplogroup (in grey) which is significantly differentiated from *D. dentex* and *D. gibbosus* as well. Actually, these specific haplotypes (HC28, HC29) are closer to the *Chemerius nufar* haplogroup (in brown) than any other species.

The Maximum Likelihood phylogenetic analysis is based on the same marker, and the same samples gave similar results (Fig 3). For this tree construction, four extra sequences of *Sparus aurata* and *Diplodus puntazzo* were added as an outgroup (S2 Table). The best substitution model was selected by PHYML was HKY85 +G. Intraspecific genetic distances in *D. dentex* were low (less than 2%) and significantly lower than interspecific distances with all other putative species. DNA barcoding successfully discriminated between most of the putative species with the exception of *D. maroccanus* and *D. angolensis* (Fig 3). Moreover this analysis shows that specimens captured in Faro, Portugal in 2002 belong to one or probably two different
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Median-joining network for the COI haplotypes identified in the common dentex DNA barcoding analysis (S2 Table). Each circle represents a haplotype and its size is proportional to haplotype frequency. Colors indicate different species in A, sampling regions in B for common dentex only, and source of all sequences used in C. Small red nodes represent possible median vectors while numbers indicate the number of nucleotide differences.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203866.g002

species (*D. gibbosus* and a "new" species with no COI sequence available in GenBank) and not to *D. dentex*.

The samples confirmed to belong to *D. dentex* by the COI marker include mainly Mediterranean specimens but also samples from the Atlantic Ocean, Bay of Biscay, and France collected in the present study (Fig 2B). Only these samples are involved in the following D-loop analysis.

**D-loop diversity in *D. dentex***

Mitochondrial DNA polymorphism analysis was based on D-loop sequences of 79 common dentex individuals from 16 sampling locations (Table 1 and S4 Table), i.e. only those confirmed to belong to *D. dentex* by COI barcoding analysis. The sequence length ranged from
<figure><img src="image_7.png" /></figure>
**Fig 3.** Maximum Likelihood phylogenetic tree constructed from COI haplotypes of *Dentex dentex* and 10 other species (S2 Table) with PHYML. Numbers on branches indicate aBayes / bootstrap support values (only values above 50% are shown).

https://doi.org/10.1371/journal.pone.0203866.g003

444 to 456 bp, a difference due to a duplication event. In total, 69 variable sites were found, with one 4-base and one 8-base duplication. Both global nucleotide diversity $\pi$ and haplotype diversity $h$ were high ($\pi$ = 0.02104, $h$ = 0.9854), indicating high levels of genetic variation in *D. dentex*. Estimations for each sampled population are summarized in Table 2. Among the 57 haplotypes identified, only few were distributed in more than one location (HD15, HD33, HD35, HD38, HD39, HD43, HD45, HD47, HD51), while most were sample-specific (see the median-joining D-loop network, Fig 4).

Nonetheless, differentiation between sampled populations was in most cases low and statistically not significant, as indicated by pairwise Fst estimations (S1 Fig). No significant spatial genetic structure was revealed by AMOVA analysis and testing of different groupings (Tamura & Nei distance matrix, 100,000 permutations). According to the AMOVA, the most probable scenario includes 3 areas/groups separated by the Gibraltar and the Siculo-Tunisian Straits: Eastern Mediterranean (Aegean, Crete, Adriatic, Lampedusa), Western Mediterranean (Tunisia, Sardinia, Corsica locations, France, Baleares, Spain) and Atlantic Ocean (Bay of Biscay). The largest component of variation (96.86%) was due to variation among individuals within populations (Fst = 0.031, $p$ = 0.11), while variation among areas represents only 1.10% of total variation (FCT = 0.011, $p$ = 0.24).

Phylogeographic analysis with median-joining networks also indicates a lack of spatial structure in *D. dentex*. In order to simplify the full median network initially produced which contained hypercubes (homoplasies), an analysis with MP option [48] was performed with NETWORK 5.0.0. for a final interpretation. The network (Fig 4) does not support the existence of any geographical structuring of haplotypes. Mismatch distribution showed a unimodal distribution of pairwise differences (S2 Fig) indicating that *D. dentex* has undergone periods of population growth (SSD = 0.00076, $p$ = 0.67). A pattern of demographic expansion was also confirmed by the raggedness index, which failed to reject the null hypothesis of sudden population growth ($r$ = 0.0044, $p$ = 0.58). The existence of distinct clades with different demographic history would have resulted in multimodal distribution.

## Microsatellites diversity in *D. dentex*

Data consists of the genotyping of 312 *D. dentex* individuals for 8 microsatellite loci (described in S1 Table). General parameters of polymorphism and Fis parameter are given in Table 3. Significant intra-specific fixation indices are generally low but can be locally high (from 0.05526 to 0.30670).

Inter-sample differentiation (Fst) values are rarely significant (11 out of 128 tests before correction) and low (< 0.032) (Table 4). No geographic organization of the genetic diversity is observed since the seven inter-Mediterranean significant tests became non-significant after Bonferroni correction [49]. Moreover, organization among samples has also been investigated through multidimensional analyses, FCA (S3 Fig). This second analysis does not allow for the distinguishing of any structure among groups except for the additional fifteen Bargelloni et al. 2003 [21] samples confirming their misidentification.

Additionally, assignment Bayesian analysis led to the same result, i.e. that no structure is evidenced (Fig 5). The $\Delta K$ Evanno's method indicates a possible cluster number for K = 2 (Ln'(K) = -41.8). For this value of K, the plot from the STRUCTURE software does not allow determining a structure; each individual is assigned to both groups in equal proportion (±
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Median-joining network of 57 D-loop haplotypes of *Dentex dentex* (S3 Table). Each circle represents a haplotype and its size is proportional to haplotype frequency. Colors indicate different geographic origin (red: Eastern Mediterranean, yellow: Western Mediterranean, blue: Atlantic). Small red nodes represent possible median vectors while numbers indicate nucleotide differences greater than one between neighboring haplotypes.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203866.g004

50%). Thus all the individuals would form the same population. Likewise, for other bigger val-
ues of K (3 to 15), the same pattern is observed, cutting in rather than between individuals.
Finally, the Mantel's test indicates a p = 0.05 probability of being right by rejecting the null
hypothesis of independence of the two matrices. This limit of significance indicates a slight
structure probably due to the geographical distance between the 17 stations.
# Discussion

## COI sequences: Species misidentification

The only genetic analysis available so far in *D. dentex* is that of Bargelloni et al. [21] using allozymes and mtDNA D-loop sequences. They also compared four other spaid species' genetic structures to that of *D. dentex*. This last species displayed the highest degree of genetic differentiation between Mediterranean and Atlantic populations. More specifically, in this study *D. dentex* presented a repeated motif in the 3' end of the sequenced region of all Mediterranean individuals that made the sequence alignment with other sequences difficult. This perturbing region was deleted and all further analyses were performed on a reduced data set of only 98 bp. This removal of unalignable regions could lead to a bias in the estimation of the divergence between the *D. dentex* haplotype groups [21]. In the same way, a morphological study using the same samples [22] have indicated a geographical gradient for the common dentex only, from Atlantic Ocean to the east of the Mediterranean.

In order to reach a rational phylogenetic picture, the COI sequences provided in the present analysis (107 sequences) were augmented by GenBank data (66 sequences). This was necessary for a comparison with other studies, but also to detect species misidentification considering COI sequences as a barcoding marker. Bargelloni et al. [21], using D-loop sequences from the same specimens from Faro, Portugal, have proposed the existence of two Mediterranean/ Atlantic clades in *D. dentex*, but this finding is strongly questioned by our results. The median-joining network presented in Fig 2, based on COI sequences, clearly points out one major haplogroup in dark blue (Fig 2A) representing the species *Dentex dentex* according to numerous GenBank sequences, and a distinct haplogroup (haplotypes HC21-HC27) in green, in which most of Bargelloni's Portuguese group and the GenBank *D. gibbosus* sequences are included.

<table><caption>Table 3. Estimated (He), unbiased (Hnb) and observed (Ho) heterozygosities and Fis values for the 17 samples of *Dentex dentex* genotyped with eight microsatellite loci.</caption><thead><tr><th></th><th>N</th><th>He</th><th>Hnb</th><th>Ho</th><th>A</th><th>F<sub>IS</sub></th><th>F<sub>IS</sub> signif.</th></tr></thead><tbody><tr><td>Heraklion (North Crete, Greece)</td><td>24</td><td>0.723</td><td>0.739</td><td>0.740</td><td>10</td><td>-0.002</td><td>ns</td></tr><tr><td>North Aegean Sea (Greece)</td><td>13</td><td>0.690</td><td>0.721</td><td>0.523</td><td>6.6</td><td>0.284</td><td>***</td></tr><tr><td>Otranto (South Adriatic, Italy)</td><td>23</td><td>0.681</td><td>0.696</td><td>0.641</td><td>9.2</td><td>0.081</td><td>*</td></tr><tr><td>Sicily (Italy)</td><td>13</td><td>0.656</td><td>0.682</td><td>0.611</td><td>5.9</td><td>0.109</td><td>*</td></tr><tr><td>Lampedusa (Italy)</td><td>13</td><td>0.689</td><td>0.716</td><td>0.721</td><td>7.5</td><td>-0.007</td><td>ns</td></tr><tr><td>NorthTunisia</td><td>14</td><td>0.649</td><td>0.677</td><td>0.476</td><td>6.6</td><td>0.307</td><td>***</td></tr><tr><td>South Tunisia</td><td>14</td><td>0.677</td><td>0.705</td><td>0.558</td><td>7.4</td><td>0.216</td><td>***</td></tr><tr><td>Sardinia (Italy)</td><td>23</td><td>0.666</td><td>0.681</td><td>0.683</td><td>7.5</td><td>-0.002</td><td>ns</td></tr><tr><td>Giraglia (Corsica, France)</td><td>25</td><td>0.699</td><td>0.714</td><td>0.675</td><td>9.7</td><td>0.055</td><td>ns</td></tr><tr><td>St Florent (Corsica, France)</td><td>18</td><td>0.681</td><td>0.702</td><td>0.590</td><td>8.0</td><td>0.165</td><td>***</td></tr><tr><td>Galeria (Corsica, France)</td><td>31</td><td>0.690</td><td>0.703</td><td>0.631</td><td>10.2</td><td>0.104</td><td>**</td></tr><tr><td>Ajaccio (Corsica, France)</td><td>15</td><td>0.687</td><td>0.712</td><td>0.568</td><td>8.2</td><td>0.208</td><td>***</td></tr><tr><td>Bonifacio (Corsica, France)</td><td>24</td><td>0.673</td><td>0.689</td><td>0.618</td><td>8.9</td><td>0.105</td><td>**</td></tr><tr><td>St Raphael (France)</td><td>14</td><td>0.699</td><td>0.725</td><td>0.741</td><td>7.4</td><td>-0.023</td><td>ns</td></tr><tr><td>Baleares (Spain)</td><td>19</td><td>0.677</td><td>0.698</td><td>0.642</td><td>8.4</td><td>0.083</td><td>*</td></tr><tr><td>Alicante (Spain)</td><td>21</td><td>0.694</td><td>0.711</td><td>0.683</td><td>9.2</td><td>0.041</td><td>ns</td></tr><tr><td>Bay of Biscay (France)</td><td>8</td><td>0.658</td><td>0.703</td><td>0.734</td><td>5.9</td><td>-0.048</td><td>ns</td></tr></tbody></table>

Ns = not significant

* = p<0.05

** = p<0.01

*** = p<0.001.

https://doi.org/10.1371/journal.pone.0203866.t003
<table><caption>Table 4. Calculation of inter-samples microsatellite differentiation for *Dentex dentex* (Weir and Cockerham 1984 estimator θ) and significance of each test.</caption><thead><tr><th></th><th>HE</th><th>NA</th><th>OT</th><th>SI</th><th>LA</th><th>NT</th><th>ST</th><th>SA</th><th>GI</th><th>SF</th><th>GA</th><th>AJ</th><th>BO</th><th>SR</th><th>BA</th><th>AL</th><th>BB</th></tr></thead><tbody><tr><td>HE</td><td>0.000</td><td>0.000</td><td>0.007</td><td>0.025</td><td>-0.008</td><td>-0.013</td><td>0.011</td><td>-0.007</td><td>-0.004</td><td>0.009</td><td>0.008</td><td>-0.017</td><td>-0.005</td><td>0.003</td><td>0.001</td><td>0.002</td><td>0.015</td></tr><tr><td>NA</td><td></td><td>0.000</td><td>0.002</td><td>*0.016</td><td>0.005</td><td>-0.007</td><td>0.005</td><td>0.005</td><td>0.000</td><td>0.001</td><td>0.003</td><td>0.000</td><td>0.000</td><td>-0.007</td><td>-0.004</td><td>0.002</td><td>0.013</td></tr><tr><td>OT</td><td></td><td></td><td>0.000</td><td>0.009</td><td>-0.003</td><td>-0.001</td><td>0.001</td><td>*0.012</td><td>-0.001</td><td>-0.001</td><td>-0.002</td><td>0.001</td><td>0.004</td><td>-0.002</td><td>-0.006</td><td>0.001</td><td>0.012</td></tr><tr><td>SI</td><td></td><td></td><td></td><td>0.000</td><td>0.010</td><td>0.005</td><td>-0.007</td><td>*0.022</td><td>0.005</td><td>-0.001</td><td>0.005</td><td>0.012</td><td>0.004</td><td>0.001</td><td>0.001</td><td>0.011</td><td>0.025</td></tr><tr><td>LA</td><td></td><td></td><td></td><td></td><td>0.000</td><td>0.006</td><td>0.003</td><td>0.013</td><td>0.001</td><td>0.006</td><td>-0.002</td><td>-0.012</td><td>0.008</td><td>0.006</td><td>0.000</td><td>0.005</td><td>*0.020</td></tr><tr><td>NT</td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>0.007</td><td>-0.004</td><td>-0.008</td><td>-0.003</td><td>-0.002</td><td>-0.013</td><td>-0.004</td><td>-0.004</td><td>-0.014</td><td>-0.007</td><td>0.011</td></tr><tr><td>ST</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>*0.026</td><td>0.003</td><td>-0.012</td><td>-0.001</td><td>0.008</td><td>0.005</td><td>-0.006</td><td>-0.003</td><td>0.004</td><td>0.004</td></tr><tr><td>SA</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.004</td><td>*0.017</td><td>*0.013</td><td>-0.003</td><td>0.008</td><td>*0.016</td><td>0.008</td><td>0.009</td><td>*0.032</td></tr><tr><td>GI</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.004</td><td>0.000</td><td>-0.005</td><td>-0.002</td><td>0.005</td><td>-0.008</td><td>-0.005</td><td>0.013</td></tr><tr><td>SF</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.005</td><td>0.004</td><td>-0.007</td><td>-0.000</td><td>-0.007</td><td>0.003</td><td>0.008</td></tr><tr><td>GA</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.004</td><td>0.006</td><td>-0.001</td><td>-0.006</td><td>-0.003</td><td>*0.019</td></tr><tr><td>AJ</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>0.010</td><td>0.006</td><td>-0.003</td><td>-0.009</td><td>0.022</td></tr><tr><td>BO</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>0.000</td><td>-0.002</td><td>0.003</td><td>0.003</td></tr><tr><td>SR</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.002</td><td>0.001</td><td>0.018</td></tr><tr><td>BA</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>-0.004</td><td>0.012</td></tr><tr><td>AL</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td><td>*0.025</td></tr><tr><td>BB</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>0.000</td></tr></tbody></table>

Significance:

* = p<0.05

** = p<0.01

*** = p<0.001.

HE = Heraklion (North Crete, Greece); NA = North Aegean Sea (Greece); OT = Otranto (South Adriatic, Italy); SI = Sicily (Italy); LA = Lampedusa (Italy);

NT = NorthTunisia; ST = South Tunisia; SA = Sardinia (Italy); GI = Giraglia (Corsica, France); SF = St Florent (Corsica, France); GA = Galeria (Corsica, France);

AJ = Ajaccio (Corsica, France); BO = Bonifacio (Corsica, France); SR = St Raphael (France); BA = Balearic Islands (Spain); AL = Alicante (Spain); BB = Bay of Biscay (France).

https://doi.org/10.1371/journal.pone.0203866.t004

The sequences used both in Bargelloni et al. [21] and Palma and Andrade [22] included a third group (in grey in Fig 2A) which is likely another unknown species. These results are further confirmed by nuclear (microsatellite) markers (S3 Fig). This type of misidentification is in fact quite common and has also occurred in other marine fishes [50, 51, 52]. As a consequence, the morphological gradient observed between Atlantic and Mediterranean seas by Palma and Andrade [22] is likely biased by the misidentification demonstrated here. The large unalignable region of D-loop reported by Bargelloni et al. [21] was likely a result of having multiple species in their data set.

In the present study, COI barcode analysis placed most Portuguese samples of Bargelloni et al. [21] within a distinct species that is probably *D. gibbosus*. Moreover, the Atlantic *D. dentex* sampled in the Bay of Biscay (close to Hendaye, France) for the present study are genetically undifferentiated from the Mediterranean populations. By its global coherence, the COI network seems very efficient in attributing the true species to each published and new sequences. Based on the current available sequences generated in the present study and from GenBank, a set of *D. dentex* samples were identified for further analyses. Support for monophyly of *D. dentex* cluster is high as shown Fig 3.

## MtDNA D-loop sequences

After attribution of the true species to each sequence, 57 *D. dentex* D-loop haplotypes were analyzed for a structure description. The D-loop median-joining network presented Fig 4 is a
<figure><img src="image_9.png" /><figcaption>Fig 5. Assignment tests for 2, 4, 6 and 15 subdivisions for *Dentex dentex* population structure analysis.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203866.g005

first method to find subgroups. However the three geographic categories (East-Mediterranean, West-Mediterranean and Atlantic populations) were randomly dispatched in the entire network evidencing no structure. Although most minority haplotypes are found in only one region, most frequent haplotypes were shared in two (HD38, HD43, HD47...) or three regions (HD51).

Such absence of structure in a very large range, from the Aegean Sea to French Atlantic regions, is quite surprising for a coastal species. However similar homogeneity has been observed in coastal-lagoonal fish species [53], jellyfish [54] and in crustacean species [55, 56], highlighting the probable importance of larval dispersal in a species range.

## Microsatellites loci

Microsatellites are generally highly variable, valuable nuclear markers. Their allele diversity generally guarantees a good description of the structure, even at a very local scale [57, 58]. Here, genotyping of 8 microsatellite loci on *D. dentex* populations from East-Mediterranean, West-Mediterranean on French and Spanish coast, Tyrrhenian islands and North African
coasts, and finally in the Bay of Biscay on the Atlantic coast, i.e. about 4–5,000 km long, showed no indication of genetic structure at all. Inter samples Fst were not significant and assignment analysis was unable to distinguish subgroups, cutting inside individuals rather than between them, even for K = 2 (S4 and S5 Figs).

Mainly, nuclear markers did not allow for the detecting of any structure in Mediterranean and Atlantic populations of *D. dentex*, i.e. in the whole sampled distribution range, confirming the mtDNA results. This absence of structure has been observed in other species [59, 60, 61, 62]. These markers also highlighted the species misidentification involved in the Bargelloni et al. study [21]. However, the Mantel test is at the limit of significance. This suggests a slight east-west Mediterranean structure reminiscent of the observed Greece-to-Spain gradient thanks on morphometric measurements [22].

## High gene flow between Mediterranean common dentex populations and between the Mediterranean and Atlantic regions

In most genetic analyses of marine species, the use of variable markers allowed the description of subgroups in the Mediterranean and Atlantic range (e.g. *Platichthys flesus* [63]; *Scomber scombrus* [60]; *Lithognathus mormyrus* [64]; *Symphodus tinca* [65]). In fact, there are only few cases showing a total lack of genetic differentiation (e.g. *Sebastes mentella* [59]; *Scomber japonicus* [60]; *Epinephelus marginatus* [61]; *Conger conger* [62], and now *D. dentex*. According to Patarnello et al. [23], comparing the life history of several species of the same family with comparable biological traits can provide new understanding. Thus, the Sparidae family presents two patterns: *Lithognathus mormyrus*, *Spondyliosoma cantharus* [21] and *Diplodus puntazzo* [66] showing a clear Atlantic-Mediterranean differentiation, unlike *Pagellus bogaraveo*, *Pagrus pagrus* [21], *Diplodus sargus* [66] and *D. dentex* (the present study).

The different analyses performed based on microsatellites and mitochondrial DNA data all lead to the same observation: the absence of genetic structure in populations of common dentex. The inter-population fixation index and AMOVA's results are not significant. The FCA does not allow for the determination of groups among individuals, except the misidentified ones. The method of assignment does not give structured results. Only the IBD would reveal a slight structure separating the Mediterranean and Atlantic individuals.

It is remarkable that the sampling, all around Mediterranean, stretches along four years (2012–2014) for the new samples and along fifteen years if we consider the comparison with the Bargelloni et al. sampling [21]. This could be a bias if the Mediterranean populations were structured and dynamic with high level exchanges. However, the picture given by several markers is a very stable structure, without genetic contrast between regions. This diachronic sampling has so no effect on the general results.

The same pattern is observed around Corsica. Individuals all seem to belong to the same panmictic population. No Fst value is significant between the different stations and the other methods employed do not show any differentiation among populations. This absence of structure at this local Corsican scale has also been investigated with a multi-method approach [25] using a combination of markers that have different spatial and temporal scales of integration: microsatellite DNA markers, otolith shape analysis and parasites communities. However, although no genetic structure was found using microsatellites, a complex population structure is suggested in this study with the other marker distinguishing different “ecological population units” around Corsica. These markers, otoliths shape and parasitological tags, revealed an ecological timescale on finer temporal and spatial scales probably without genetic consequence. The holistic approach makes it possible to take into account phenotypic characteristics and thus to define more precisely some structures depending on local ecology and migrations [67].
The genetic homogeneity of many marine species is thought commonly to be due to two factors that minimize accumulation of genetic differences among populations: a large effective population size that limits genetic drift and life history characteristics that favor dispersal in continuous dynamic oceanic environments [68, 69]. This high connectivity between sampling sites may be due to the migration of adults, which are capable of large movements and/or larval dispersal. This larval dispersal may be, however, more limited than expected [70] by distance, for example. The IBD at the limit of significance and the holistic approach of Marengo et al. [25] would support this hypothesis. Among the three markers chosen for the holistic method, only the genetic nuclear markers did not detect any structure.

## Perspectives on the species distribution in the Atlantic Ocean

The established distribution range of *D. dentex* will probably change in the future for two reasons. First, in the Atlantic Ocean, the common dentex is considered present from the British Isles to Senegal [14]. But the misidentification described above could call into question this vast range. In the present survey, this question has not been fully developed (only one sample); more samples have to be done in the Atlantic Ocean, including the real distribution of the common dentex. The individuals observed at the extremities of the zone could belong to similar species. The species distribution in the Atlantic Ocean is scarce and more limited than previously thought. A detailed description of the sparis present along the European and African Atlantic coasts is necessary, sustained by DNA barcoding. Secondly, in a few years, global warming could again challenge the area of distribution. Indeed, the common dentex prefers a warm water environment [71]. Global warming is gradually impacting the temperature of Mediterranean waters. The seasonal stratification is still normal, but it is subject to oscillations. This phenomenon therefore has a positive impact on the common dentex, which can then extend its range to the north, along the Atlantic Ocean coasts.

## Supporting information

**S1 Table.** Information on the eight microsatellite loci used in the present study and the PCR conditions. (CL168 and CL1014 [72]; Ds33 and Dxd16 [73]; SaGT41b [74]; SaI19 [75]; SauE82INRA and SauI41INRA [76]).
(DOCX)

**S2 Table.** All information used for COI gene sequence analyses in Figs 2 and 3. np: not published.
(DOCX)

**S3 Table.** D-loop haplotypes for each individual.
(DOCX)

**S4 Table.** COI and D-loop sequences included in mtDNA analyses.
(DOCX)

**S5 Table.** *Dentex dentex* genotypes at 8 microsatellites. (0 = no amplification).
(DOCX)

**S1 Fig.** Heat-map of D-loop pairwise Fst's between all pairs of localities. Asterisk indicates $P \le 0.05$ in terms of statistical significance.
(TIFF)

**S2 Fig.** Mismatch distribution for D-loop haplotypes of *Dentex dentex*.
(TIFF)
**S3 Fig. Factorial Correspondences Analysis (FCA) of the whole *Dentex dentex* specimens.**

In the blue envelope are the very homogeneous samples of Mediterranean and Bay of Biscay origin. The red envelope gathers the Bargelloni et al. 2003 Portugal samples. Clearly, Bargelloni's Portuguese fish are heterogeneous (at least three clouds) and do not belong to *D. dentex* (except two individuals: see red arrows). Consequently, the Bay of Biscay samples are common dentex similar to those of Mediterranean populations.
(TIFF)

**S4 Fig. Diagram from Structure Harvester expressing Delta K according to K.**

(TIFF)

**S5 Fig. Diagram from Structure Harvester expressing likelihood according to K. Standard**

error bars are also indicated.

(TIFF)

**Acknowledgments**

The authors are grateful to Martin Laporte (University Laval, Quebec, Canada) for his help in autocorrelation calculations.

**Author Contributions**

**Conceptualization:** Patrick Berrebi, Chiheb Fassatoui, Eric D. H. Durieux.

**Funding acquisition:** Patrick Berrebi, Eric D. H. Durieux.

**Investigation:** Alexiane Viret, Dimitris Tsaparis, Costas S. Tsigenopoulos, Patrick Berrebi, Eric D. H. Durieux.

**Methodology:** Alexiane Viret, Dimitris Tsaparis, Costas S. Tsigenopoulos, Patrick Berrebi, Eric D. H. Durieux.

**Project administration:** Eric D. H. Durieux.

**Resources:** Eric D. H. Durieux.

**Supervision:** Eric D. H. Durieux.

**Writing – original draft:** Alexiane Viret, Dimitris Tsaparis, Costas S. Tsigenopoulos, Patrick Berrebi, Eric D. H. Durieux.

**Writing – review & editing:** Alexiane Viret, Dimitris Tsaparis, Costas S. Tsigenopoulos, Patrick Berrebi, Andrea Sabatini, Marco Arculeo, Chiheb Fassatoui, Antonios Magoulas, Michel Marengo, Beatriz Morales-Nin, Nathalie Caill-Milly.

**References**

1. Hauser L, Carvalho GR. Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts. Fish and Fisheries 2008; 9(4): 333–62.

2. Selkoe KA, Henzler CM, Gaines SD. Seascape genetics and the spatial ecology of marine populations. Fish and Fisheries 2008; 9(4): 363–77.

3. Chopelet J, Waples RS, Mariani S. Sex change and the genetic structure of marine fish populations. Fish and Fisheries 2009; 10(3): 329–43.

4. Nesbø CL, Rueness EK, Iversen SA, Skagen DW, Jakobsen KS. Phylogeography and population history of Atlantic mackerel (*Scomber scombrus* L.): a genealogical approach reveals genetic structuring among the eastern Atlantic stocks. Proceedings of the Royal Society of London Series B 2000; 267 (1440): 281–92. https://doi.org/10.1098/rspb.2000.0998 PMID: 10714883
5. Beheregaray LB, Sunnucks P. Fine-scale genetic structure, estuarine colonization and incipient speciation in the marine silverside fish *Odontesthes argentinensis*. Molecular Ecology 2001; 10(12): 2849–66. PMID: 11903897

6. Plough LV. Genetic load in marine animals: a review. Current Zoology 2016; 62(6): 567–79. https://doi.org/10.1093/cz/zow096 PMID: 29491946

7. Calò A, Félix-Hackradt FC, Garcia J, Hackradt CW, Rocklin D, Otón JT, et al. A review of methods to assess connectivity and dispersal between fish populations in the Mediterranean Sea. Advances in Oceanography and Limnology 2013; 4(2): 150–75.

8. Ward RD. Genetics in fisheries management. In: Marine Genetics. Springer, Dordrecht; 2000. p. 191–201.

9. Cowen RK, Paris CB, Srinivasan A. Scaling of connectivity in marine populations. Science 2006; 311 (5760): 522–7. https://doi.org/10.1126/science.1122039 PMID: 16357224

10. Waples RS, Punt AE, Cope JM. Integrating genetic data into management of marine resources: how can we do it better? Fish and Fisheries 2008; 9(4): 423–49.

11. Laikre L, Palm S, Ryman N. Genetic population structure of fishes: implications for coastal zone management. AMBIO: A Journal of the Human Environment 2005; 34(2): 111–9.

12. Hanel R, Tsigenopoulos CS. Phylogeny, evolution and taxonomy of sparids with some notes on their ecology and biology. In: Pavlidis MA, Mylonas CC, editors. Sparidae. Wiley-Blackwell; 2011. p. 51–73.

13. Froese R, Pauly D. Editors. FishBase. World Wide Web electronic publication. www.fishbase.org, version (02/2017)

14. Marengo M, Durieux EDH, Marchand B, Francour P. A review of biology, fisheries and population structure of *Dentex dentex* (Sparidae). Reviews in Fish Biology and Fisheries 2014; 24(4): 1065–88.

15. Gimenez G, Estévez A, Lahnsteiner F, Zecevic B, Bell JG, Henderson RJ, et al. Egg quality criteria in common dentex (*Dentex dentex*). Aquaculture 2006; 260(1–4): 232–43.

16. Morales-Nin B, Moranta J. Life history and fishery of the common dentex (*Dentex dentex*) in Mallorca (Balearic Islands, western Mediterranean). Fisheries Research 1997; 30(1–2): 67–76.

17. Loir M, Le Gac F, Somarakis S, Pavlidis M. Sexuality and gonadal cycle of the common dentex (*Dentex dentex*) in intensive culture. Aquaculture 2001; 194(3–4): 363–81.

18. Grau A, Saborido-Rey F, Pastor E, Palmer M, Massutí-Pascual E, Quetglas A, et al. Reproductive strategy of common dentex *Dentex dentex*: management implications. Mediterranean Marine Science 2016; 17(2): 552–66.

19. Baudouin M, Marengo M, Pere A, Culioli J-M, Santoni M-C, Marchand B, et al. Comparison of otolith and scale readings for age and growth estimation of common dentex *Dentex dentex*. Journal of Fish Biology 2016; 88(2): 760–6. https://doi.org/10.1111/jfb.12816 PMID: 26563912

20. Marengo M, Culioli J-M, Santoni M-C, Marchand B, Durieux EDH. Comparative analysis of artisanal and recreational fisheries for *Dentex dentex* in a Marine Protected Area. Fisheries Management and Ecology 2015; 22(3): 249–60.

21. Bargelloni L, Alarcon JA, Alvarez MC, Penzo E, Magoulas A, Reis C, et al. Discord in the family Sparidae (Teleostei): divergent phylogeographical patterns across the Atlantic-Mediterranean divide. Journal of Evolutionary Biology 2003; 16(6): 1149–58. PMID: 14640406

22. Palma J, Andrade JP. Morphological study of *Pagrus pagrus*, *Pagellus bogaraveo*, and *Dentex dentex* (Sparidae) in the eastern Atlantic and the Mediterranean Sea. Journal of the Marine Biological Association of the United Kingdom 2004; 84(2): 449–54.

23. Patarnello T, Volckaert FAM, Castilho R. Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break? Molecular Ecology 2007; 16(21): 4426–44. https://doi.org/10.1111/j.1365-294X.2007.03477.x PMID: 17908222

24. Sacco F, Marrone F, Lo Brutto S, Besbes A, Nfati A, Gatt M, et al. The Mediterranean Sea hosts endemic haplotypes and a distinct population of the dolphinfish *Coryphaena hippurus* Linnaeus, 1758 (Perciformes, Coryphaenidae). Fisheries Research 2017; 186(Part 1): 151–8.

25. Marengo M, Baudouin M, Viret A, Laporte M, Berrebi P, Vignon M, et al. Combining microsatellite, otolith shape and parasites community analyses as a holistic approach to assess population structure of *Dentex dentex*. Journal of Sea Research 2017; 128(Supplement C): 1–14.

26. Estoup A, Largiader CR, Perrot E, Chourrout D. Rapid one-tube DNA extraction for reliable PCR detection of fish polymorphic markers and transgenes. Molecular Marine Biology and Biotechnology 1996; 5 (4): 295–8.

27. Zemlak TS, Ward RD, Connell AD, Holmes BH, Hebert PDN. DNA barcoding reveals overlooked marine fishes. Molecular Ecology Resources 2009; 9: 237–42. https://doi.org/10.1111/j.1755-0998.2009.02649.x PMID: 21564983
28. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PD. DNA barcoding Australia's fish species. Philosophical transactions of the Royal Society of London B 2005; 360(1462): 1847-57.

29. Meyer A, Kocher TD, Basasibwaki P, Wilson AC. Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial-DNA Sequences. Nature 1990; 347: 550-553. https://doi.org/10.1038/347550a0 PMID: 2215680

30. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution 2013; 30(12): 2725-9. https://doi.org/10.1093/molbev/mst197 PMID: 24132122

31. Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution 1999; 16(1): 37-48. https://doi.org/10.1093/oxfordjournals.molbev.a026036 PMID: 10331250

32. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biology 2010; 59(3): 307-21. https://doi.org/10.1093/sysbio/syq010 PMID: 20525638

33. Lefort V, Longueville J-E, Gascuel O. SMS: Smart model selection in PhyML. Molecular Biology and Evolution 2017; 34(9): 2422-2424. https://doi.org/10.1093/molbev/msx149 PMID: 28472384

34. Anisimova M, Gil M, Dufayard JF, Dessimoz C, Gascuel O. Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Systematic Biology 2011; 60(5): 685-99. https://doi.org/10.1093/sysbio/syr041 PMID: 21540409

35. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009; 25(11): 1451-2. https://doi.org/10.1093/bioinformatics/btp187 PMID: 19346325

36. Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 1992; 131(2): 479-91. PMID: 1644282

37. Harpending RC. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Human Biology 1994; 66: 591-600. PMID: 8088750

38. Belkhir K, Borsa P, Goudet J, and Bonhomme F. 2004. GENETIX 4.05: logiciel sous Windows pour la génétique des populations. Laboratoire Génome et Population, CNRS-UPR, Université de Montpellier II, Montpellier, France. 2004

39. Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 1978; 89(3): 583-90. PMID: 17248844

40. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution 1984; 38(6): 1358-70. https://doi.org/10.1111/j.1558-5646.1984.tb05657.x PMID: 28563791

41. Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Research 1967; 27(2 Part 1): 209-20. PMID: 6018555

42. Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 1997; 145(4): 1219-28. PMID: 9093870

43. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 2003; 164(4): 1567-87. PMID: 12930761

44. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics 2000; 155(2): 945-59. PMID: 10835412

45. Pritchard J. Structure. 2012. http://pritchardlab.stanford.edu/structure.html

46. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Molecular Ecology 2005; 14(8): 2611-20. https://doi.org/10.1111/j.1365-294X.2005.02553.x PMID: 15969739

47. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources 2012; 4(2): 359-61.

48. Polzin T, Daneshmand SV. On Steiner trees and minimum spanning trees in hypergraphs. Operations Research Letters 2003; 31(1): 12-20.

49. Rice WR. Analyzing tables of statistical tests. Evolution 1989; 43(1): 223-5. https://doi.org/10.1111/j.1558-5646.1989.tb04220.x PMID: 28568501

50. Rock J, Costa FO, Walker DI, North AW, Hutchinson WF, Carvalho GR. DNA barcodes of fish of the Scotia Sea, Antarctica indicate priority groups for taxonomic and systematics focus. Antarctic Science 2008; 20(3): 253-62.

51. Ward RD, Hanner R, Hebert PDN. The campaign to DNA barcode all fishes, FISH-BOL. Journal of Fish Biology 2009; 74(2): 329-56. https://doi.org/10.1111/j.1095-8649.2008.02080.x PMID: 20735564
52. Garcia-Vazquez E, Machado-Schiaffino G, Campo D, Juanes F. Species misidentification in mixed hake fisheries may lead to overexploitation and population bottlenecks. Fisheries Research 2012; 114 (Supplement C): 52–5.

53. Boissin E, Micu D, Janczyszyn-Le Goff M, Neglia V, Bat L, Todorova V, et al. Contemporary genetic structure and postglacial demographic history of the black scorpionfish, *Scorpaena porcus*, in the Mediterranean and the Black Seas. Molecular Ecology 2016; 25(10): 2195–209. https://doi.org/10.1111/mec.13616 PMID: 26989881

54. Glynn F, Houghton JDR, Bastian T, Doyle TK, Fuentes V, Lilley MKS, et al. High-resolution genetic analysis reveals extensive gene flow within the jellyfish *Pelagia noctiluca* (Scyphozoa) in the North Atlantic and Mediterranean Sea. Biological Journal of the Linnean Society 2016; 117(2): 252–63.

55. García-Merchán VH, Robainas-Barcia A, Abelló P, Macpherson E, Palero F, García-Rodríguez M, et al. Phylogeographic patterns of decapod crustaceans at the Atlantic-Mediterranean transition. Molecular Phylogenetics and Evolution 2012; 62(2): 664–72. https://doi.org/10.1016/j.ympev.2011.11.009 PMID: 22138160

56. Marra A, Mona S, Sà RM, D'Onghia G, Maiorano P. Population genetic history of *Aristeus antennatus* (Crustacea: Decapoda) in the western and central Mediterranean Sea. PLOS ONE 2015; 10(3): e0117272. https://doi.org/10.1371/journal.pone.0117272 PMID: 25775363

57. Bryja J, Smith C, Konečný A, Reichard M. Range-wide population genetic structure of the European bit-terling (*Rhodeus amarus*) based on microsatellite and mitochondrial DNA analysis. Molecular Ecology 2010; 19(21): 4708–22. https://doi.org/10.1111/j.1365-294X.2010.04844.x PMID: 20958813

58. Gagnaire PA, Minegishi Y, Zenboudji S, Valade P, Aoyama J, Berrebi P. The genetic footprint of allopat-ric isolation in *Anguilla marmorata*: introgression across a semi-permeable barrier to gene flow high-lights fine population structure. Evolution 2011; 65: 3413–3427. https://doi.org/10.1111/j.1558-5646.2011.01404.x PMID: 22133215

59. Roques S, Sévigny J-M, Bernatchez L. Genetic structure of deep-water redfish, *Sebastes mentella*, populations across the North Atlantic. Marine Biology 2002; 140(2): 297–307.

60. Zardoya R, Castilho R, Grande C, Favre-Krey L, Caetano S, Marcato S, et al. Differential population structuring of two closely related fish species, the mackerel (*Scomber scombrus*) and the chub mack-erel (*Scomber japonicus*), in the Mediterranean Sea. Molecular Ecology 2004: 1785–98. https://doi.org/10.1111/j.1365-294X.2004.02198.x PMID: 15189203

61. Schunter C, Carreras-Carbonell J, Planes S, Sala E, Ballesteros E, Zabala M, et al. Genetic connectivity patterns in an endangered species: The dusky grouper (*Epinephelus marginatus*). Journal of Experi-mental Marine Biology and Ecology 2011; 401(1): 126–33.

62. Miralles L, Lee LS, Borrell Y, Garcia-Vazquez E. Population genetic structure of the European conger (*Conger conger*) in North East Atlantic and West Mediterranean Sea. Fisheries Research 2016; 174 (Supplement C): 245–9.

63. Borsa P, Blanquer A, Berrebi P. Genetic structure of the flounders *Platichthys flesus* and *P. stellatus* at different geographic scales. Marine Biology 1997; 129(2): 233–46.

64. Sala-Bozano M, Ketmaier V, Mariani S. Contrasting signals from multiple markers illuminate population connectivity in a marine fish. Molecular Ecology 2009; 18(23): 4811–26. https://doi.org/10.1111/j.1365-294X.2009.04404.x PMID: 19863720

65. Carreras C, Ordóñez V, Zane L, Kruschel C, Nasto I, Macpherson E, et al. Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation. Scientific Reports 2017; 7.

66. Bargelloni L, Alarcon JA, Alvarez MC, Penzo E, Magoulas A, Palma J, et al. The Atlantic-Mediterranean transition: Discordant genetic patterns in two seabream species, *Diplodus puntazzo* (Cetti) and *Diplo-dus sargus* (L.). Molecular Phylogenetics and Evolution 2005; 36(3): 523–35. https://doi.org/10.1016/j.ympev.2005.04.017 PMID: 15936957

67. Begg GA, Waldman JR. An holistic approach to fish stock identification. Fisheries Research 1999; 43 (1): 35–44.

68. Waples RS. Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. Journal of Heredity 1998; 89(5): 438–50.

69. Ward RD, Woodwark M, Skibinski DOF. A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes. Journal of Fish Biology 1994; 44(2): 213–32.

70. Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Annual Review of Marine Science 2009; 1(1): 443–66.

71. Aguzzi J, Doya C, Tecchio S, Leo FCD, Azzurro E, Costa C, et al. Coastal observatories for monitoring of fish behaviour and their responses to environmental changes. Reviews in Fish Biology and Fisheries 2015; 25(3): 463–83.
72. Vogiatzi E, Lagnel J, Pakaki V, Louro B, Canario AVM, Reinhardt R, et al. In silico mining and characterization of simple sequence repeats from gilthead sea bream (*Sparus aurata*) expressed sequence tags (EST-SSRs); PCR amplification, polymorphism evaluation and multiplexing and cross-species assays. Marine Genomics 2011; 4(2): 83–91. https://doi.org/10.1016/j.margen.2011.01.003 PMID: 21620329

73. De la Herrán R, Magoulas A, Garrido-Ramos MA, Ruiz-Rejón C, Ruiz-Rejón M, Zouros E. Desarrollo de microsatélites en tres especies de espáridos con interés en acuicultura. In: La acuicultura como actividad económica de las zonas costeras. IX Congreso Nacional de Acuicultura. Ed. Consejería de Agricultura y Pesca. Junta de Andalucía. Sevilla. Spain, 2005. pp133-136.

74. Batargias C, Dermitzakis E, Magoulas A, Zouros E. Characterization of six polymorphic microsatellite markers in gilthead seabream, *Sparus aurata* (Linnaeus 1758). Molecular Ecology 1999; 8: 897–8. PMID: 10368971

75. Brown RC, Tsalavouta M, Terzoglou V, Magoulas A, Mcandrew BJ. Additional microsatellites for *Sparus aurata* and cross-species amplification within the Sparidae family. Molecular Ecology Notes 2005; 5(3): 605–7.

76. Launey S, Krieg F, Haffray P, Bruant JS, Vannier A, Guyomard R. Twelve new microsatellite markers for gilted seabream (*Sparus aurata* L.): characterization, polymorphism and linkage. Molecular Ecology Notes 2003; 3(3): 457–9.