Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Mitochondrial DNA Reveals Genetic Structuring of Pinna nobilis across the Mediterranean Sea

  • Daria Sanna,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Piero Cossu,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Gian Luca Dedola,

    Affiliation Dipartimento di Medicina Veterinaria - Sezione di Anatomia, Università di Sassari, Sassari, Italy

  • Fabio Scarpa,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Ferruccio Maltagliati,

    Affiliation Dipartimento di Biologia, Università di Pisa, Pisa, Italy

  • Alberto Castelli,

    Affiliation Dipartimento di Biologia, Università di Pisa, Pisa, Italy

  • Piero Franzoi,

    Affiliation Dipartimento di Scienze Ambientali, Informatica e Statistica, Università Cà Foscari, Venezia, Italy

  • Tiziana Lai,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Benedetto Cristo,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Marco Curini-Galletti,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Paolo Francalacci,

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

  • Marco Casu

    marcasu@uniss.it

    Affiliation Dipartimento di Scienze della Natura e del Territorio - Sezione di Zoologia, Archeozoologia e Genetica, Università di Sassari, Sassari, Italy

Abstract

Pinna nobilis is the largest endemic Mediterranean marine bivalve. During past centuries, various human activities have promoted the regression of its populations. As a consequence of stringent standards of protection, demographic expansions are currently reported in many sites. The aim of this study was to provide the first large broad-scale insight into the genetic variability of P. nobilis in the area that encompasses the western Mediterranean, Ionian Sea, and Adriatic Sea marine ecoregions. To accomplish this objective twenty-five populations from this area were surveyed using two mitochondrial DNA markers (COI and 16S). Our dataset was then merged with those obtained in other studies for the Aegean and Tunisian populations (eastern Mediterranean), and statistical analyses (Bayesian model-based clustering, median-joining network, AMOVA, mismatch distribution, Tajima’s and Fu’s neutrality tests and Bayesian skyline plots) were performed. The results revealed genetic divergence among three distinguishable areas: (1) western Mediterranean and Ionian Sea; (2) Adriatic Sea; and (3) Aegean Sea and Tunisian coastal areas. From a conservational point of view, populations from the three genetically divergent groups found may be considered as different management units.

Introduction

The assessment of genetic variation in marine species and the evaluation of connectivity among populations are crucial components for conservation purposes and resource management [1], [2]. However, direct estimates of dispersal are difficult to obtain for marine organisms [3] because of the multiple factors (e.g., pelagic larval duration, larval phylopatry, regional and local water circulation, and habitat specificity) that influence dispersal capability [4]. Molecular tools are often employed to make indirect inferences about the levels of connectivity among natural populations [5], [6]. Molecular approaches can be particularly useful for determining the influence of such factors as hydrodynamic conditions [7], [8], habitat specificity [9], [10] and various ecological parameters (e.g., [11], [12]) on gene flow and spatial genetic structure. However, the spatial genetic structure of a species may reflect not only its current dispersal capability but also its phylogeographic history. For example, the occurrence of different genetic lineages in Mediterranean species has been explained as the result of either the present dispersal capability or past geomorphologic processes or both (see [13] and references therein). Furthermore, human activities (e.g., intense exploitation, industrialisation, and coastal tourism) often lead to significant reduction of population size, which may result in an increased magnitude of genetic drift.

The fan mussel Pinna nobilis (Linnaeus, 1758) (Mollusca: Bivalvia) is an endemic Mediterranean species of great conservation interest that can be regarded as a flag species sensu Walpole and Leader-Williams [14] and Heywood [15]. As one of the largest bivalves, with maximum lengths in excess of 1 m, this popular species is able to capture the public imagination [16]. Pinna nobilis occurs at depths between 0.5 and 60 m in soft bottoms primarily characterised by seagrass meadows [16]. Its veliger-stage larvae drift in the water column [17] before settling in the sediment and anchoring via a byssus [18]. Although De Gaulejac and Vicente [19] hypothesised a period of 5–10 days, little is known about the length of the larval life cycle of P. nobilis. Moreover, Peharda and Vilibić [20] postulated that P. nobilis veliger might exhibit daily vertical migration, spending the daylight hours in darker, deeper waters and occupying superficial layers at night, as observed for other bivalves [21].

As a result of the population declines associated with human activity (primarily the harvesting of byssus for the manufacture of so-called “sea silk” and shell collecting [22], [23], secondarily the use as food of the abductor muscle in some Mediterranean regions [24], and the collection of the little, with poor commercial value, pearls [25]), P. nobilis has been included in the list of Mediterranean endangered species (Annex IV of the Habitat Directive and Annex II of the Barcelona Convention) since 1995. Remarkably, it is likely that some of the known geographic and ecological barriers [26], [27], [28], which define different marine ecoregions and biogeographic sectors in the Mediterranean [29], [30], may affect the degree of connectivity among P. nobilis populations. In the Mediterranean region, several authors have documented population subdivision in relation to physical barriers in other marine species with a high dispersal potential (e.g., [13], [31], [32]).

Despite the conservation importance of P. nobilis, few studies have investigated the biology and ecology of this species (e.g., [20], [23], [33], [34]), and fewer still have studied its population genetics. To date, only two molecular surveys have been conducted, one in the Aegean Sea [35] and one along the Tunisian coasts [36]. The former [35] surveyed the genetic variability of four Aegean populations by means of the mitochondrial DNA (mtDNA) markers Cytochrome c Oxidase subunit I (COI) and 16S ribosomal DNA (16S). A high level of haplotypic diversity was found for the COI gene, whereas the 16S gene showed a lower level of variability. These results suggested the lack of genetic structuring among the Aegean populations. In the second study [36], the authors used COI sequences to describe the genetic variability of five populations from the northern, eastern, and southern Tunisian coasts. A North-East decreasing gradient of genetic variability was found among these populations, which was explained in light of the variance in the hydrodynamic regime of the areas analysed.

Knowledge of the amount of genetic variability and distribution in space and time is crucial for a correct diagnosis of the conservation status and viability of populations and the threats to them [37], [38]. Accordingly, it is evident that a molecular analysis performed on a larger number of Mediterranean P. nobilis populations could provide a deeper insight into the possible influence of environmental or anthropogenic stress on the population dynamics and evolutionary history of this endangered species, allowing the development of effective conservation measures. Based on this background, the present study was aimed to investigate the large-scale patterns of spatial genetic variation of P. nobilis in the Mediterranean. Our primary focus was the appraisal of the genetic structuring and population connectivity of P. nobilis in the western Mediterranean, Ionian Sea, and Adriatic Sea ecoregions. We surveyed populations from the Sardinian-Corsican region, Elba Island, Sicily, and the Venetian Lagoon and compared these populations with the available data from the Aegean Sea [35] and Tunisian coasts [36]. This large area is characterised by a complex and variable pattern of marine circulation (e.g., [39]) and encompasses several biogeographic sectors [28]. To accomplish the objective of the present study, we employed two mitochondrial markers (COI and16S) previously used for the investigation of Aegean [35] and Tunisian populations [36].

Materials and Methods

Ethics Statement

No field studies involving impacting manipulation, dislocation, or removal of Pinna nobilis individuals were performed. For each location under protection all necessary permits were obtained for the sampling activities by the authority responsible for each protected area:

- Gianfranco Russino, Director of the Marine Protected Area of Capo Caccia-Isola Piana (year 2010); location: Baia di Porto Conte;

- Augusto Navone, Director of the Marine Protected Area of Tavolara-Punta Coda Cavallo (year 2010); locations: Molara, Monte Petrosu;

- Mauro Gargiulo, Director of the National Park Arcipelago di La Maddalena (year 2011); location La Maddalena;

- Lorenzo Mascia, Director of the Marine Protected Area of Penisola del Sinis-Isola di Mal di Ventre (year 2011); locations: Oristano, Isola di Mal di Ventre;

- Maddy Cancemi, Director of the Réserve Naturelle des Bouches de Bonifacio (year 2011); location: Isola Piana.

No specific permissions were required for the locations placed out of protected areas, since those locations are not privately-owned or protected in any way.

Sampling

Overall, we sampled 236 specimens of P. nobilis from 25 locations (Table 1 and Figure 1), sited in the following ecoregions according to Spalding et al. [30]: (1) the central part of the western Mediterranean (the Sardinian Sea, the North Tyrrhenian Sea, the South Tyrrhenian Sea, and the Strait of Sicily); (2) the Ionian Sea (the South-East of Sicily); and (3) the Adriatic Sea (the Venetian Lagoon). Two individuals from the Levantine Sea (Cyprus) were also analysed as the easternmost Mediterranean outliers. We developed a specific non-lethal sampling method, performed by SCUBA divers, which does not cause significant damages to the shell and soft tissues of P. nobilis. The valves of a given individual were held open with a wooden stick (diameter = 0.5 cm), put in proximity (4–5 cm) of the hinge ligament, and a 20–50 mg sample of mantle tissue was excised using Hartman Alligator Forceps 3.5. The stick was then removed, and the tissue sample stored in a 5 ml tube. The tissue sample was then transferred to the laboratory in a refrigerated box and there preserved in 75% ethanol. This method ensures the survival of the sampled individuals.

thumbnail
Figure 1. Mediterranean sampling plan.

Map of the Mediterranean indicating the localisation of the sampling sites. The relative geographic coordinates are reported in Table 1. The populations are labelled as reported in Table 1.

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

DNA Extraction and PCR

DNA was isolated using the Qiagen DNeasy tissue kit. Mitochondrial regions were amplified using specific primers for COI (L: 5′-GGTTGAACTATHTATCCNCC-3′ and H: 5′-GAAATCATYCCAAAAGC-3′) and 16S (L: 5′-TGCTCAATGCCCAAGGGGTAAAT-3′ and H: 5′-AACTCAGATCACGTAGGG-3′) designed by the authors, since those provided by Katsares et al. [35] did not give satisfactory results.

Each 25 µl PCR mixture contained approximately 100 ng of total genomic DNA, 0.32 µM of each primer, 2.5 U of EuroTaq DNA Polymerase (Euroclone), 1× reaction buffer and 200 µM of each dNTP. MgCl2 concentration was set at 3 mM, and 12.5 µg of BSA were added to the reaction mixture. PCR amplifications were performed according to the following steps: 1 cycle of 2 min at 94°C, 35 cycles of 1 min at 94°C, 1 min at 46°C and 1 min 30 s at 72°C. A post-treatment of 5 min at 72°C was applied. After electrophoresis on 2% agarose gels, the PCR products were purified using ExoSAP-IT (USB Corporation) and sequenced using an external sequencing core service (Macrogen Europe). Dual peaks of similar height, which could be interpreted as evidence of mitochondrial pseudogenes in the nucleus (Numts) or heteroplasmy, were not observed in any of the electropherograms. The PCR products did not show any occurrence of aspecificity, excluding the possibility of multiple nuclear mtDNA-like sequences.

Statistical Analysis

Sequences were aligned using the program Clustal W [40] and deposited in GenBank (Accession Nos.: COI, JX854788-JX855023; 16S, JX854562-JX854787). Estimates of the number of polymorphic sites (S), the number of haplotypes (H), haplotype diversity (h), nucleotide diversity (π), and the mean number of pairwise differences (d) were obtained using the software package DnaSP 5.10 [41].

After having added to our COI and 16S sequences the COI [35], [36] and 16S [35] sequences available in the literature, we constructed two different datasets. One dataset included 236 COI sequences obtained in the present study along with 26 and 49 sequences from the Aegean Sea [35] and Tunisian coasts [36], respectively (see Table 1 for details). The second dataset included the concatenated COI and 16S regions and consisted of the 219 sequences of the present study along with 25 from the Aegean Sea [35]. Concatenating these two mitochondrial genes was appropriate because COI and 16S did not show significant heterogeneity by partition-homogeneity test (P = 0.06) performed with PAUP* 4.0b10 [42].

a) Analysis of the COI dataset.

The presence of population genetic structure was assayed by the Bayesian model-based clustering algorithm implemented in BAPS 5.2 [43]. Clustering was performed using the module for linked molecular data and by applying the codon linkage model, which is appropriate for sequence data. Each analysis was run ten times with a vector of K-values = 1 to 22, each with six replicates. Haplotypes were organised into haplogroups following the partition of sequences into the distinct genetic clusters evidenced by Bayesian clustering. The Bayesian clustering algorithms probabilistically assign individuals to groups based on nucleotide frequencies of DNA sequence data without presuming pre-defined populations [44]. For each sampling location, we computed the proportion of a given haplogroup to build maps of haplogroup frequency in the Mediterranean.

Genetic relationships among haplotypes were investigated by a median-joining network [45] using the software package Network 4.5.0.1 (www.fluxus-engineering.com). Mutations were inversely weighted according to the number of times they originated. Therefore, different weights were assigned to the point mutations: sites involved in a single mutational event = 90; sites involved in repeated mutational events = 90/number of repetitions.

Patterns of genetic differentiation at the population level were assessed by Arlequin 3.5.1.3 [46] using a matrix of Tamura and Nei’s [47] genetic distances with a gamma correction according to the best-fitting model of sequence evolution obtained with jModeltest 2.1.1 [48]. First, we calculated pairwise ΦST values between sampling localities. Samples with less than five individuals were excluded from population level analysis due to a lack of statistical power. We next conducted a hierarchical analysis of molecular variance (AMOVA) [49]. Alternative groupings of populations were taken into account. Initially, we ran analyses in which populations were grouped a posteriori according to the pattern of genetic structure suggested by the pairwise ΦST values. Significance of both pairwise ΦST values and AMOVA was assessed by a permutation test (with 10,000 random replicates). Then, we defined a priori regional groups that corresponded to the main biogeographic subdivisions of the Mediterranean [29]: (1) western vs. eastern Mediterranean; (2) western vs. eastern Mediterranean vs. Adriatic. Sicilian populations located close to the border separating the two Mediterranean basins were twofold tested to be either entirely part of western Mediterranean, or part of western (San Vito lo Capo-SVC, Mondello-MON, Milazzo-MLZ) and eastern (Pachino-PAC, Ognina-OGN) Mediterranean. Where necessary, probability values were corrected for multiple testing by applying the False Discovery Rate method [50].

Aspects of historical demography were inferred by examining the fit of the observed mismatch distribution of DNA pairwise differences to a model of demographic population expansion [51], [52]. Such distributions are unimodal when populations have experienced recent expansion and multimodal at demographic equilibrium or when populations are significantly subdivided. Mismatch distributions and relative Rogers’ parameters [53] (τ, θ0, θ1) were assessed using Monte-Carlo simulations of 10,000 random samples as implemented in Arlequin 3.5.1.3. The fit of the observed mismatch distribution to the expected distribution under a model of demographic expansion was assessed by Monte-Carlo simulations of 10,000 random samples. The sum of squared deviations (SSD) between observed and expected mismatch distributions was used to test the probability of obtaining a simulated SSD that was larger than or equal to the observed one.

Tajima’s D [54] and Fu’s Fs [55] neutrality tests were used to infer departures from population equilibrium models. Significant negative values of Tajima’s D are expected to occur in cases of recent population expansion or after a selective sweep [55]. Positive values are expected in the case of balancing selection, population subdivision or recent bottlenecks [56]. Similarly, significant negative Fu’s Fs values indicate an excess of rare haplotypes, which might be caused by recent population expansion [55], whereas positive values indicate balancing selection, population structure or moderate bottlenecks [57]. Combining different neutrality tests can help to distinguish among the different evolutionary processes responsible for departures from equilibrium; Fu’s Fs can better detect demographic expansions, whereas Tajima’s D can better detect bottlenecks and populations contractions [58].

b) Analysis of the concatenated COI and 16S dataset.

Statistical treatments were conducted as with the COI dataset. Note that for this dataset, redefined population groupings were required for some analyses (e.g., AMOVA) due to the lack of 16S sequences from Tunisian coasts.

c) Estimates of divergence time.

The modal value (τ) of the mismatch distribution of pairwise genetic differences was also used to estimate time since expansion using the formula t = τ/2u [51]. Here, t is the time since expansion in generations and u is the mutation rate of the entire DNA fragment, equal to u = µs, where µ is the mutation rate per nucleotide per generation, and s is the sequence length. Uncertainty in the estimates of expansion time since expansion was accounted for computing the 95% confidence interval (CI) of τ from the parametric bootstrap approach with 10,000 pseudo-replicates [52]. Time since expansion expressed in years was estimated assuming a constant molecular clock, a 2-years generation time [17], and a divergence rate per nucleotide site per million years of 0.52% as estimated in Luttikhuizen et al. [59] for bivalves.

We also used the Bayesian coalescent approach implemented in BEAST 1.7.4 [60] to infer historical demography, as it allows improved recovery of the historical signal within DNA sequences [61]. As there was a lack of 16S sequences for Tunisian populations and of an appropriate divergence rate for 16S region for the species, analysis was carried out on the COI dataset only.

We used the Bayesian skyline plot model [62] assuming a piecewise constant model with ten coalescent intervals; Markov-Chain-Monte-Carlo simulations were run under the GTR+G+I model with four gamma categories using a relaxed uncorrelated molecular clock [62]. Divergence rate of the COI region was modeled using a uniform distribution that ranged from 0.14% to 0.52% (divergence rate per site per million years) according to Luttikhuizen et al. [59]. For each simulation three independent runs were carried out; depending on the size (number of sequences in the dataset) 100 million up to 400 million iterations were used (thinning was adjusted to obtain a final sample of 10,000 records). Tracer 1.5 [63] was used to assess convergence of runs and the Effective Sample Size (ESS) of each parameter. In order to obtain an adequate effective sample size (ESS ≥200), the three independent runs performed for each simulation were combined using the Logcombiner utility of the BEAST 1.7.4 package. The resulting file was used to estimate population size change through time, that was visualized by the Bayesian skyline plot computed with Tracer 1.5.

Results

COI Dataset

Specific primers amplified an internal portion of COI 338 base pair (bp) long. Therefore, after the addition to our data of those from Katsares et al. [35] and Rabaoui et al. [36], we obtained a 338 bp sequence alignment from 311 Pinna nobilis individuals. Sixty-two different haplotypes (Table S1), defined by 45 polymorphic sites, were found. In our dataset (236 individuals), many haplotypes were found only at a single locality (63%). This percentage increased to 79% when the samples from the Aegean Sea [35] and Tunisian coasts [36] were included. Total mean haplotype and nucleotide diversity, calculated on 311 individuals, were h = 0.910 and π = 0.007, respectively. The lowest values of haplotype and nucleotide diversity were found among Aegean and Tunisian populations (h = 0.720 and π = 0.004 for Aegean Sea and h = 0.605 and π = 0.002 for Tunisian coasts); overall higher values of genetic diversity were found in the other sites. Estimates of COI genetic diversity for each site are reported in Table S2.

a) Genetic structuring.

The Bayesian assignment analysis identified four distinct haplotype groups, hereafter denominated P1, P2, P3 and P4 (Table S3). The groups P1, P3 and P4 includeed samples from Sardinia, Corsica, Elba Island, Sicily, the Venetian Lagoon and Cyprus and lack internal geographic structuring, whereas the Aegean and Tunisian samples all exhibited P2 haplotypes, except for one Aegean individual from Epanomi-EP and one Tunisian specimen from Monastir-M, both belonging to group P3. The geographic distributions of haplogroups over the Mediterranean map are provided in Figure 2. The group P1 was the most widespread (41.8%), being present in Sardinia (with the exception of all samples from Isola di Mal di Ventre-IMV), Corsica, Elba Island, Sicily and the Venetian Lagoon (two individuals only). The group P2 was found in 23.5% of individuals and characterised populations from the Aegean Sea and Tunisian coasts. Although P3 was the group spread most widely across the Mediterranean (occurring in all populations except for the majority of the Aegean and Tunisian samples), it occurred at low frequency (26.7%), reaching higher proportion in areas as distant as the northern Adriatic and western Sardinia. The group P4 was characterised by the lowest average frequency (8%) and occurred in north-eastern Sardinia, southern Corsica, at Elba Island, northern Sicily (plus one individual from south-eastern Sicily, Pachino-PAC) and in the Venetian Lagoon (two individuals only).

thumbnail
Figure 2. COI dataset: Bayesian cluster distribution.

Frequency distribution of the four groups of haplotypes P1, P2, P3, and P4, as evidenced by Bayesian analysis over the Mediterranean map. Due to the high number of populations, the size of the pie charts in the Sardinian-Corsican region was reduced, and the relative populations were merged in three distinct groups as reported below. North-western Sardinia group includes Baia di Porto Conte-BPC, Torre del Porticciolo-POR, Lazzaretto-LAZ, Ospedale Marino-OSM; central-eastern Sardinia group includes Molara-MOL, Capo Ceraso-CCE, Le Saline-SAL, Monte Petrosu-MPE, Porto Ottiolu-OTT; central-western Sardinia group includes Oristano-ORI, Marceddì-MAR, Isola di Mal di Ventre-IMV.

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

The median-joining network analysis identified a low level of divergence among haplotypes, as, with a few exception, connections were characterised by one or two point mutations, (Figure S1, and Figure 3). Three haplotypes (PN1, PN6, PN7; Table S1), present in 134 of 311 (43.1%) individuals, were the most common among all samples, with the exception of the Aegean and Tunisian samples, of which 85.3% (64 of 75 individuals) showed two main haplotypes (PN52, PN54; Table S1). No evidence of geographic structuring among regions (all characterised by many instances of haplotype sharing) was detected for the entire sample, with the exception of populations from the Aegean Sea and Tunisian coasts that did not share haplotypes with the remaining populations (Figure S1). When haplotypes were grouped on the network according to the results of Bayesian assignment, each cluster (P1, P2, P3 or P4) was generally characterised by one or two widely-dispersed haplotypes, with radiating star-like branching patterns (Figure 3).

thumbnail
Figure 3. COI dataset: network analysis.

Median-joining network (on the left) with haplotypes grouped according to the results of the Bayesian assignment. Small red plots on the nodes, labelled as “mv”, show median vectors representing hypothetic connecting sequences, calculated with a maximum parsimony method. The long branches leading to isolated haplotypes were shortened and indicated with “\\”. The distribution of haplogroups within populations (on the right) are displayed in the histogram. The x axis reports populations and the y axis the absolute frequency of distribution. The populations are labelled as reported in Table 1. The number of mutations on the network branches are reported in the Figure S1.

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

The pairwise ΦST estimates indicated significant genetic differentiation between the Aegean and Tunisian samples and the other sites (Table S4). The Venetian Lagoon was also significantly differentiated in most comparisons, with the only exceptions being three Sardinian populations (Le Saline-SAL, Porto Ottiolu-OTT, Marceddì-MAR) and one Sicilian (Milazzo-MLZ). The AMOVA (Table 2) maximized the largest differences among groups defined a posteriori when the samples from the Aegean Sea and Tunisian coasts were grouped together and the samples from Sardinia, Corsica, Elba Island, Sicily and the Venetian Lagoon were considered a unique separate group (ΦCT = 0.369, P<0.001) (Table 2, A). However, a similar value of molecular variance (ΦCT = 0.367, P<0.001) was obtained when treating the Venetian Lagoon as a separate third group (Table 2, B). When alternative a priori groupings of samples were tested, based on biogeographic criteria, the AMOVA showed a decrease in the proportion of ΦCT variance (Table 2, C and D).

b) Historical demography.

The mismatch distribution carried out on the entire dataset showed a unimodal distribution of pairwise DNA differences (Figure 4A), fitting the Rogers and Harpending’s model [51] of demographic expansion (SSD = 0.002, P>0.05). When assuming the occurrence of discrete groups of populations that, according to AMOVA, represent three geographic areas of (1) Sardinia, Corsica, Elba Island and Sicily, (2) the Venetian Lagoon (Adriatic Sea), and (3) Aegean Sea and Tunisian coasts, we found that the two former groups displayed a unimodal mismatch distribution fitting the demographic expansion model (Figures 4B and 4C), whereas the mismatch distribution of the Aegean and Tunisian samples (Figure 4D) did not correspond to the expected distribution under a demographic expansion model (SSD = 0.014, P<0.05). Neutrality tests showed a significant departure from equilibrium, with the exception of Tajima’s D for the Venetian Lagoon (Adriatic) population (Figure 4C). The negative values resulting from the neutrality tests were consistent with population expansions in all three groups. The non-significant value of D encountered in the Venetian Lagoon may reflect the reduced power of this test compared to the Fu’s Fs to detect demographic expansions [58].

thumbnail
Figure 4. COI dataset: mismatch analysis.

Graphs of the mismatch distributions of (A) the entire sample, (B) Sardinia, Corsica, Elba Island, and Sicily, (C) the Venetian Lagoon, and (D) the Aegean Sea and Tunisian coasts. The x axis reports the observed distribution of pairwise nucleotide differences, and the y axis reports the frequencies. N, sample sizes; SSD, sum of squared deviations; D, Tajima’s D value; Fs, Fu’s F value. The significant values are given in bold.

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

The time since expansion calculated according to Rogers and Harpending [51] for P. nobilis in Mediterranean was estimated to be the Pleistocene. In particular, for populations from Sardinia, Corsica, Elba Island and Sicily, we found that the demographic expansion started around 1.6 Mya (95% CI: 0.512–2.643 Mya). A similar estimate was calculated for the Venetian Lagoon population, with time since expansion dating back to about 1.3 Mya (95% CI: 0.643–1.956 Mya). A more recent demographic expansion was inferred for the Aegean-Tunisian group starting around 0.5 Mya (95% CI: 0.319–0.848 Mya).

The demographic history inferred by the Bayesian skyline plot revealed a signature of demographic expansion starting approximately 0.270 Mya and ending around 0.100 Mya (Figure 5A). However, only samples from Sardinia, Corsica, Elba Island and Sicily showed a demographic expansion that corresponded to this time interval (Figure 5B). Conversely, the Venetian Lagoon (Figure 5C) and Aegean-Tunisian samples (Figure 5D) had a virtually constant size throughout their demographic history.

thumbnail
Figure 5. COI dataset: Bayesian skyline plots.

Bayesian skyline plots depicting the demographic history of Pinna nobilis. The effective population size (y-axis) is shown as a function of time (x-axis), with the units expressed as KYR (1,000 years). (A) entire dataset, (B) the regional group encompassing the Sardinian-Corsican region, Elba Island, and Sicily, (C) the Venetian lagoon and (D) the Aegean Sea and Tunisian coasts. The black lines denote the median population estimates, and the grey lines denote the 95% high posterior density (HPD) confidence limits.

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

Concatenated COI and 16S Dataset

The concatenated dataset comprised 338 bp for COI and 450 bp for 16S. The sequences of 219 specimens were determined here. Twenty-five COI and 16S sequences from Katsares et al. [35] were added. Among these 244 individuals, 103 haplotypes (Table S5), defined by 78 polymorphic sites (S) were found. Total mean haplotype diversity and nucleotide diversity were h = 0.961 and π = 0.005, respectively. A large number of haplotypes (77%) was found only at a single locality. The Aegean populations showed the lowest average values of haplotype diversity and nucleotide diversity (h = 0.333 and π = 0.001), whereas high levels of genetic diversity were found for almost all other sites. Estimates of genetic diversity for the concatenated COI and 16S dataset are provided in Table S2. Furthermore, estimates of genetic diversity and haplotype frequencies for the 16S region are given in Table S2, and Table S6, respectively.

a) Genetic structuring.

The Bayesian analysis identified four haplotype groups, hereafter denominated N1, N2, N3 and N4 (Table S7). The most frequent group (N2), representing half of the individuals analysed, included the two most common haplotypes (PNCS 1, PNCS 6; Table S5), found in all populations, with a few exceptions in Sardinia (Isola di Mal di Ventre-IMV), Cyprus and the Aegean Sea (Epanomi-EP, Aggeloxori-AG, and Korinthiakos Gulf-KO). The second group (N3) (32.8% of individuals) included samples from Sardinia, Corsica, Elba Island, Sicily, the Venetian Lagoon and Cyprus. The groups N1 and N4 were less frequent, occurring in 9% and 8.2% of individuals only, respectively. N1 occurred in individuals from Sardinia, Corsica, Elba Island, Sicily and the Venetian Lagoon, while N4 included individuals from all Aegean populations with the exception of those from Xios-XI.

The haplogroups frequency distribution largely resembled that obtained by the COI analysis (Figure 6). In particular, groups N1 and N3 approximately mirrored the distributions of P4 and P3, respectively. A striking exception is represented by haplogroup N2, whose geographic distribution overlapped that of the COI haplogroup P1, spreading as far as the Aegean Sea. At Xios-XI, the N2 group reached frequencies as high as those found in western Mediterranean and Ionian Sea populations. In contrast, the geographic distribution of haplogroup N4 was restricted to the Aegean Sea.

thumbnail
Figure 6. COI-16S dataset: Bayesian cluster distribution.

Frequency distribution of the four groups of haplotypes N1, N2, N3 and N4, as evidenced by Bayesian analysis over the Mediterranean map. Due to the high number of populations, the size of the pie charts in the Sardinian-Corsican region was reduced, and the relative populations were merged in three distinct groups as reported below. North-western Sardinia group includes Baia di Porto Conte-BPC, Torre del Porticciolo-POR, Lazzaretto-LAZ, Ospedale Marino-OSM; central-eastern Sardinia group includes Molara-MOL, Capo Ceraso-CCE, Le Saline-SAL, Monte Petrosu-MPE, Porto Ottiolu-OTT; central-western Sardinia group includes Oristano-ORI, Marceddì-MAR, Isola di Mal di Ventre-IMV.

https://doi.org/10.1371/journal.pone.0067372.g006

The median-joining network of haplotypes (Figure S2 and Figure 7) showed a pattern of high genetic variability distributed across the Mediterranean. The two most frequent haplotypes (PNCS 1, PNCS 6; Table S5) were found in 23.8% of individuals and were shared only among individuals from Sardinia, Elba Island and Sicily. Different star-like branching patterns radiating out from the inferred root haplotypes occurred in the network (Figure S2). A few point mutations (1 to 6) distinguished haplotypes from each other, and a diffused unresolved reticulation of branches was found among individuals. Haplotypes belonging to the Venetian Lagoon and Aegean Sea populations occupied a peripheral position in the network. The two Cyprian haplotypes were very divergent from one another (Figure S2). No evidence of spatial structuring was apparent at either the local or large geographic scales. However, an exception was represented by Aegean individuals, which did not share haplotypes with other samples. Considering groupings detected by Bayesian assignment analysis (Figure 7), the network showed that the most diffused group (N2) was characterised by two main haplotypes and many derived ones, that diverged for a single point mutation. The group N3 was also characterised by a high level of haplotype variation, with unique haplotypes radiating from the central ones. The group N1 exhibited a similar pattern, although a reduced level of polymorphism was found. The group N4 showed the lowest level of diversity.

thumbnail
Figure 7. COI-16S dataset: network analysis.

Median-joining network (on the left) with haplotypes grouped according to the results of the Bayesian assignment. Small red plots on the nodes, labelled as “mv”, show median vectors, representing hypothetic connecting sequences, calculated with a maximum parsimony method. The long branches leading to isolated haplotypes were shortened and indicated with “\\”.The x axis reports populations and the y axis the absolute frequency of distribution. The populations are labelled as reported in Table 1. The number of mutations on the network branches are reported in the Figure S2.

https://doi.org/10.1371/journal.pone.0067372.g007

The pairwise ΦST estimates obtained for the concatenated dataset (Table S8) were similar to those reported for COI. A significant genetic differentiation occurred between the Aegean samples of Epanomi-EP and Aggeloyesori-AG, and the other samples. As previously noticed, the Venetian Lagoon sample was significantly differentiated from the others in most comparisons with the only exceptions being two Sardinian populations (Le Saline-SAL, Porto Ottiolu-OTT). Nonetheless, combining data from different genes resulted in a larger number of significant comparisons (Table S8). In particular, samples from Corsica appeared to be genetically differentiated not only from samples outside the Sardinian-Corsican region, but also from their northern Sardinian counterparts, and samples from Xios-XI were significantly different from those from the Aegean Sea.

Consistently with results of the COI analysis, AMOVA (Table 3) detected maximum genetic differentiation among groups defined a posteriori on the pairwise ΦST values, when samples from the Aegean Sea were grouped separately from the others (ΦCT = 0.396, P<0.001) (Table 3, A). A similarly high level of molecular variance was found when the Venetian Lagoon was considered as a third group (ΦCT = 0.356, P<0.001) (Table 3, B). AMOVA tested on a priori groupings of samples also showed a decrease in the proportion of variance (Table 3, C and D); particularly no significant differences were found between groups corresponding to the western and eastern Mediterranean (ΦCT = 0.045, P>0.05) (Table 3, C and D).

b) Historical demography.

The results of demographic analyses of the concatenated dataset partly agreed with those of the COI dataset. The mismatch distribution carried out on the entire dataset fitted the model of demographic expansion (SSD = 0.001, P>0.05), with a unimodal distribution of pairwise DNA differences (Figure 8A). In contrast to the results obtained for the COI dataset, all the three groups, (1) Sardinia, Corsica, Elba Island and Sicily (Figure 8B); (2) the Venetian Lagoon (Figure 8C); and (3) the Aegean Sea (Figure 8D), displayed unimodal mismatch distributions fitting the model of demographic expansion. Nonetheless, only the group including samples from Sardinia, Corsica, Elba Island and Sicily (Figure 8B) showed significant negative values for Tajima’s D and Fu’s Fs, which can be consistent with a population expansion model. The non-significant values of neutrality tests obtained for the populations from the Venetian Lagoon and the Aegean Sea (Figure 8C, 8D) were likely to reflect the reduced power of the tests to detect departures from equilibrium, due to small sample size. Indeed, when we used only samples from Aegean Sea, Fu's Fs statistic did not displayed a significant departure from equilibrium for the COI dataset too (data not shown).

thumbnail
Figure 8. COI-16S dataset: mismatch analysis.

Graphs of the mismatch distributions of (A) the entire sample, (B) Sardinia, Corsica, Elba Island, and Sicily, (C) the Venetian Lagoon and (D) the Aegean Sea. The x axis reports the observed distribution of pairwise nucleotide differences, and the y axis reports the frequencies. N, sample sizes; SSD, sum of squared deviations; D, Tajima’s D value; Fs, Fu’s F value. The significant values are given in bold.

https://doi.org/10.1371/journal.pone.0067372.g008

Discussion

This study provides the first analysis on mtDNA genetic variation of Pinna nobilis populations from a wide study area, roughly corresponding to the central part of the western Mediterranean, Ionian Sea and Adriatic Sea ecoregions (sensu Spalding et al. [30]). Supplementing the sequences obtained in the present study with those from previous investigations conducted in the Aegean Sea and along the Tunisian coasts provided a deep insight into (1) the large-scale patterns of spatial genetic variation of P. nobilis and (2) the role of Mediterranean biogeographic boundaries in shaping the species’ genetic structuring.

The high variability of the mtDNA COI region of P. nobilis further supports the effectiveness of this marker in investigating the genetic architecture of marine species (see [64] and references therein). In addition, the sample sizes adopted in the present survey allowed us to identify 16S haplotypes undetected in the previous study by Katsares et al. [35]. As the results obtained from the analysis of the COI sequences and of the merged COI and 16S sequences were largely consistent, we discuss them jointly in the following sections.

Genetic Structure of Pinna nobilis in the Mediterranean

Our study showed the occurrence of a moderate level of genetic structuring in P. nobilis populations in the study area. The low genetic divergence among the haplogroups and haplotypes, as evinced by the reduced number of point mutations separating the sequences from each other, suggested a common origin of the P. nobilis populations. However, Bayesian assignment analysis, pairwise ΦST values and AMOVA noted the occurrence of at least two main genetically divergent groups of populations: (1) the Sardinian-Corsican region, Elba Island and Sicily (western Mediterranean and Ionian Sea) and (2) the Aegean Sea and Tunisian coasts (eastern Mediterranean).

We detected a primary, significant pattern of genetic structuring along the West-East direction, with samples from the Aegean Sea (excluding the Xios-XI sample for the concatenated COI and 16S dataset) and Tunisian coasts, which were genetically differentiated from the remaining samples. This pattern could be explained as a possible consequence of a reduced gene flow between the western and eastern Mediterranean. Several studies have shown that a hindrance to gene flow exists between the western and eastern Mediterranean, even for species with high dispersal capabilities and large population sizes (e.g., [65], [66], [67]). This result was partially consistent with findings reported for other bivalves with a high potential for dispersal, along with significant genetic differentiation and a pattern of isolation by distance (e.g., [68], [69], [70]).

Other studies have also shown that the genetic break between the western and eastern Mediterranean corresponds to the Siculo-Tunisian Strait (e.g., [71]); nonetheless, for different marine organisms the geographic position of the genetic break may shift eastward with respect to the strait (e.g., [29], [72] and references therein). In our case, the pattern of genetic structuring between the Aegean with Tunisian samples and other populations (Sardinian-Corsican region, Elba Island and Sicily) did not mirror the West-East separation of the Mediterranean along the Siculo-Tunisian Strait (see [73]). Moreover, the sample from the southern coast of the strait (Tunisian populations), which was homogeneous with the Aegean populations, was genetically divergent from the samples of the northern coast of the strait (Ionian Sicily), suggesting the occurrence of a further genetic break in a North-South direction. This finding could be due to the hydrodynamic characteristics of this area: the Algerian Current, which flows eastward from the Strait of Gibraltar, along the northern coasts of Africa and passes over the Siculo-Tunisian Strait [39], represents a hydrodynamic barrier between the facing coasts of Tunisia and Sicily. Such a barrier, which reduces the connectivity among populations via marine currents, has also been invoked for the identification of distinct mitochondrial lineages in the vermetid gastropod Dendropoma petraeum in adjacent locations in Tunisia and Sicily [74]. Furthermore, the genetic affinity of the Sicilian samples with the western Mediterranean ones, led us to hypothesise a shift of the boundary of the genetic break eastward from the Siculo-Tunisian Strait, as proposed by Pérès and Picard [75] and Bianchi [29]. Although few studies on population genetics involve this area, a similar picture has been reported for Patella rustica [76], whereby the western mtDNA lineage reaches as far as the southern coast of the Italian Peninsula, well beyond the Siculo-Tunisian Strait. A genetic break located at the level of the southern tip of Calabria has also been described for Posidonia oceanica [77]. However, at this stage of knowledge, we cannot rule out that the genetic pattern observed around the western-eastern edge of the Mediterranean may reflect a primary genetic break in a North-South direction along the Siculo-Tunisian Strait. Furthermore, as has already been reported for other species [77], the genetic homogeneity found between the samples from the Tunisian coasts and the Aegean Sea could be ascribed to the persistence of a relict population in the eastern, semi-enclosed basin, rather than to an effective gene flow between the two Mediterranean basins. More individuals from unsampled regions of both western and eastern Mediterranean are needed to disentangle this point.

Our results also suggest the occurrence of a weak but significant genetic divergence of the Venetian Lagoon population from the above-mentioned two groups of populations. This finding is further evidence of the semi-enclosed nature of the Adriatic Sea, which represents a well-defined phylogeographic region within the Mediterranean [13]. A genetic discontinuity between the Adriatic Sea and the rest of Mediterranean has also been found in other species, such as the common cuttlefish Sepia officinalis [78], and the sea urchin Paracentrotus lividus [79]. Furthermore, several studies have revealed the occurrence of a genetic break among the Adriatic and Aegean seas, that was explained by the hydrographical isolation of these two basins [80], [81], [82], [83].

Interestingly, samples from Cyprus, at the easternmost part of the Mediterranean, exhibited haplotypes typical of the westernmost populations. The human-mediated introduction of planktonic larvae passively transported in the ballast water of commercial ships is a common finding in marine molluscs invasion [84], [85], [86], [87] and may account for this outcome. However, the small number of samples from Cyprus prevented any statistically supported inference.

Historical Demography

The results of historical demographic analysis and the estimates of time since expansion lead to the hypothesis of a Pleistocene scenario for P. nobilis, with expansions/contractions of its populations due to sea-level fluctuations of this period, repeatedly leading to the isolation, or partial isolation, of the Mediterranean basins [88]. Conversely, vicariant events, leading to ancient divergence of populations, are unlikely to have occurred in P. nobilis because we would expect a bimodal mismatch distribution of pairwise DNA differences for the entire dataset, which may occur when different genetic lineages are present [11], [89]. It should be noted that the estimates of time since expansion may be biased by the lack of a mutation rate specific for P. nobilis and the time dependency of molecular rate estimates [90]. In the present investigation, we used a range of inter-specific divergence rates estimated for bivalves [59]. However, the estimated molecular rates were much higher over short timescales at the intra-specific level than between species ([90], [91] and references therein but see also [92]). Despite these caveats, our results, nonetheless, provide initial insights into the dynamics of the colonisation of P. nobilis populations in the Mediterranean.

All populations analysed in this study are likely to have shared a common origin. Our data suggest an eastward expansion across the Mediterranean, likely mediated by marine currents (e.g., the Algerian Current), and followed by one or more founder events, that lead to lower diversity and private haplotypes in the Aegean and Tunisian populations. A similar trend was also reported for populations of the sea cucumber Holothuria mammata from the Aegean Sea [65]. In particular, the group of populations including the Sardinian-Corsican region, Elba Island, Sicily and the Venetian Lagoon samples survived the Late Pleistocene climatic cycles. Our data also suggest that the Sardinian-Corsican region, Elba Island and Sicily populations further experienced their last demographic expansion during the Pleistocene glaciation. Conversely, no evidence of demographic expansion was detected for the Venetian Lagoon, Aegean and Tunisian populations, with the exception of a weak trace of expansion that may have occurred later in the Early-Middle Pleistocene.

Final Remarks

Rabaoui et al. [36] suggested that the conservation and management of P. nobilis should benefit from improvements in the sampling plan (with the inclusion of samples from the western Mediterranean). Although part of the range of the species remains unsampled, our study represents an effort in this direction, as it triples the existing genetic information on P. nobilis with new data from a wide area of the western Mediterranean.

The general picture revealed by our study is the presence of a large, homogeneous group of populations of P. nobilis spanning two Mediterranean marine ecoregions (the western Mediterranean and Ionian Sea), which genetically diverges from the Adriatic population and those from the Aegean Sea and Tunisian coasts. The semi-enclosed nature of the Adriatic Sea may explain the genetic divergence of the Venetian Lagoon population. Hindrances to gene flow, that are related to biogeographic boundaries among these Mediterranean sectors, and the Pleistocene changes in sea levels, can be invoked to explain the pattern of genetic structuring shown in this study. Our results further suggest the occurrence of a West-East genetic break located eastward of the Siculo-Tunisian Strait, and also provide evidence of a further genetic break in this area which is in the North-South direction. In this context, it would be worth noting that other samples from the western (e.g. Spain and Balearic Islands) and eastern Mediterranean (e.g. from Egypt and Israel) would be very informative for the understanding of the species’ evolutionary history.

Interestingly, although the data for the local or regional levels of genetic erosion due either to human disturbance or to bio-ecological changes are unavailable, the higher levels of mitochondrial variability found in the Sardinian-Corsican region, Elba Island, Sicily and the Venetian Lagoon relative to the populations from the Aegean Sea and Tunisian coasts support the occurrence of an eastward demographic expansion. The lack of strong genetic divergence among these three geographic areas (e.g., only one point mutation separates the Aegean and Tunisian samples from those from Sardinia, Corsica, Elba Island, Sicily and the Venetian Lagoon) is evidence of a common origin of P. nobilis populations in the Mediterranean.

From a conservation point of view, the three genetically divergent groups, (1) the samples from Sardinia, Corsica, Elba Island and Sicily, (2) the samples from the Aegean Sea and Tunisian coasts and (3) the Venetian Lagoon sample, should be considered as different management units. Moreover, the samples from the Venetian Lagoon may be considered a “peripheral isolate” sensu Frey [93], with particular conservation relevance, given its peculiar geographic position.

Supporting Information

Figure S1.

COI dataset: network analysis. Median-joining network showing the haplotypes relationships among Pinna nobilis populations. Small red plots on the nodes, labelled as “mv”, show median vectors, representing hypothetic connecting sequences, calculated with a maximum parsimony method. Haplotypes diverge each other for a single mutation except where Arabic numbers on network branches indicate the occurrence of a higher number of point mutations. Populations are labelled as reported in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s001

(TIF)

Figure S2.

COI-16S dataset: network analysis. Median-joining network showing the haplotypes relationships among Pinna nobilis populations. Small red plots on the nodes, labelled as “mv”, show median vectors, representing hypothetic connecting sequences, calculated with a maximum parsimony method. Haplotypes diverge each other for a single mutation except where Arabic numbers on network branches indicate the occurrence of a higher number of point mutations. Populations are labelled as reported in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s002

(TIF)

Table S1.

COI dataset: haplotype frequencies. Frequency distribution of COI haplotypes in 311 individuals from 34 populations of Pinna nobilis. N: absolute frequency; %: relative frequency within Mediterranean populations. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s003

(DOC)

Table S2.

COI, COI-16S, 16S datasets: genetic divergence estimates. Sample sizes and genetic diversity estimates obtained for the mitochondrial regions analysed in Pinna nobilis. N: sample sizes; S: number of polymorphic sites; H: number of haplotypes; h: haplotype diversity; π: nucleotide diversity; d: mean of pairwise nucleotide differences. Populations are labelled as in Table 1. Sites with gaps were not considered.

https://doi.org/10.1371/journal.pone.0067372.s004

(DOC)

Table S3.

COI dataset: Bayesian COI haplogroup frequencies. N: absolute frequency; %: relative frequency within Mediterranean populations of Pinna nobilis. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s005

(DOC)

Table S4.

COI dataset: pairwise ΦST values among sampling localities. Pairwise ΦST values between Pinna nobilis populations with at least five individuals. Significance was assessed by permutation test. Significant values after correction for multiple testing are reported in bold. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s006

(DOC)

Table S5.

COI-16S dataset: haplotype frequencies. Frequency distribution of mitochondrial haplotypes in 244 individuals from 29 populations of Pinna nobilis. N: absolute frequency; %: relative frequency within Mediterranean populations. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s007

(DOC)

Table S6.

16S dataset: haplotype frequencies. Frequency distribution of haplotypes in 251 individuals from 29 populations of Pinna nobilis. N: absolute frequency; %: relative frequency within Mediterranean populations. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s008

(DOC)

Table S7.

COI-16S dataset: Bayesian COI-16S mitochondrial region haplogroup frequencies. N: absolute frequency; %: relative frequency within Mediterranean populations of Pinna nobilis. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s009

(DOC)

Table S8.

COI-16S dataset: pairwise ΦST values among samples localities. Pairwise ΦST values between Pinna nobilis populations with at least five individuals. Significance was assessed by permutation test. Significant values after correction for multiple testing are reported in bold. Populations are labelled as in Table 1.

https://doi.org/10.1371/journal.pone.0067372.s010

(DOC)

Acknowledgments

We are indebted with the entire staff of the Marine Protected Areas of Capo Caccia-Isola Piana, Tavolara-Punta Coda Cavallo, Penisola del Sinis-Isola di Mal di Ventre, the National Park Arcipelago di La Maddalena, and the Réserve Naturelle des Bouches de Bonifacio, and in particular with Dr. Sarah Caronni, Dr. Francesco Mura, Dr. Pier Panzalis and Dr. Alberto Ruiu for their precious collaboration during the sampling campaigns.

Author Contributions

Conceived and designed the experiments: DS MC. Performed the experiments: DS GLD FS TL BC. Analyzed the data: DS PC FM MC. Contributed reagents/materials/analysis tools: MCG AC P. Franzoi P. Francalacci MC. Wrote the paper: DS PC MC.

References

  1. 1. Palumbi SR (2004) Marine reserves and ocean neighborhoods: The spatial scale of marine populations and their management. Annu Rev Environ Resour 29: 31–68.
  2. 2. OSPAR (2006) Report on North Sea pilot project on ecological quality objectives. OSPAR Biodiversity Series 2006/239. 126 p.
  3. 3. Kinlan BP, Gaines SD (2003) Propagule dispersal in marine and terrestrial environments: a community perspective. Ecology 84: 2007–2020.
  4. 4. Kelly RP, Palumbi SR (2010) Genetic structure among 50 species of the northeastern Pacific rocky intertidal community. PLoS One 5: e8594.
  5. 5. Jones GP, Almany GR, Russ GR, Sale PF, Steneck RS, et al. (2009) Larval retention and connectivity among populations of corals and reef fishes: history, advances and challenges. Coral Reefs 28: 307–325.
  6. 6. Horne JB, Momigliano P, Welch DJ, Newman SJ, van Herwerden L (2011) Limited ecological population connectivity suggests low demands on self-recruitment in a tropical inshore marine fish (Eleutheronema tetradactylum: Polynemidae). Mol Ecol 20: 2291–2306.
  7. 7. Wing SR, Botsford LW, Largier JL, Ralston S (1998) Meroplankton distribution and circulation associated with a coastal retention zone in the northern California upwelling zone. Limnol Oceanogr 43: 1710–1721.
  8. 8. Knutsen H, Jorde PE, André C, Stenseth NC (2003) Fine scaled geographic population structuring in a highly mobile marine species: the Atlantic cod. Mol Ecol 12: 385–394.
  9. 9. Pitchford JW, Brindley J (2001) Prey patchiness predator survival and fish recruitment. Bull Math Biol 63: 527–546.
  10. 10. Biktashev VN, Brindley J (2004) Phytoplankton blooms and fish recruitment rate: effects of spatial distribution. Bull Math Biol 66: 233–259.
  11. 11. Petersen CH (2007) Historical demography and contemporary spatial genetic structure of an estuarine crab in the northeast Pacific (Hemigrapsus oregonensis). Mar Biol 150: 1289–1300.
  12. 12. Chen G, Hare MP (2008) Cryptic ecological diversification of a planktonic estuarine copepod, Acartia tonsa. Mol Ecol 17: 1451–1468.
  13. 13. Patarnello T, Volckaert FAMJ, Castilho R (2007) Pillars of Hercules: is the Atlantic–Mediterranean transition a phylogeographical break? Mol Ecol 16: 4426–4444.
  14. 14. Walpole MJ, Leader-Williams N (2002) Tourism and flagship species in conservation. Biodivers Conserv 11: 543–547.
  15. 15. Heywood VH (1995) Global biodiversity assessment. Cambridge: University Press. 1440 p.
  16. 16. Zavodnik D, Hrs-Braenko M, Legac M (1991) Synopsis on the fan shell Pinna nobilis L. in the eastern Adriatic Sea. In: Boudouresque CF, Avon M, Gravez V, editors. Les Espèces Marines à Protéger en Méditerranée. Marseille: GIS Posidonie publications. 169–178.
  17. 17. Katsanevakis S (2005) Population ecology of the endangered fan mussel Pinna nobilis in a marine lake. Endang Species Res 1: 51–59.
  18. 18. Butler A, Vicente N, De Gaulejac B (1993) Ecology of the pterioid bivalves Pinna bicolor Gmelin and P. nobilis L. Marine Life. 3: 37–45.
  19. 19. De Gaulejac B, Vicente N (1990) Ecologie de Pinna nobilis (L.), mollusque bivalve sur les côtes de Corse. Essais de transplantation et expériences en milieu contrôlé. Haliotis 10: 83–100.
  20. 20. Peharda M, Vilibić I (2008) Modelling the recruitment effect in a small marine protected area: the example of saltwater lakes on the Island of Mljet (Adriatic Sea). Acta Adriat 49: 25–35.
  21. 21. Gosling E (2003) Bivalve molluscs: biology, ecology and culture. Oxford: Fishing News Books. 443 p.
  22. 22. Richardson CA, Peharda M, Kennedy H, Kennedy P, Onofri V (2004) Age, growth rate and season of recruitment of Pinna nobilis (L) in the Croatian Adriatic determined from Mg: Ca and Sr: Ca shell profiles. J Exp Mar Biol Ecol 299: 1–16.
  23. 23. Rabaoui L, Tlig-Zouari S, Katsanevakis S, Hassine OKB (2010) Modelling population density of Pinna nobilis (Bivalvia) on the eastern and southeastern coast of Tunisia. J Molluscan Stud. 76: 340–347.
  24. 24. Greenwald M (1996) The cruising chef cookbook. Arcata, California (USA): Paradise Cay publication. 300 p.
  25. 25. Kunz GF, Stevenson CH (1993) The book of the pearl; the history, art, science, and industry of the queen of gems. Mineola, New York (USA): Dover Publications Inc. 548 p.
  26. 26. Astraldi M, Balopoulos S, Candela J, Font J, Gacic M, et al. (1999) The role of straits and channels in understanding the characteristics of Mediterranean circulation. Prog Oceanogr 44: 65–108.
  27. 27. Béranger K, Mortier L, Gasparini GP, Gervasio L, Astraldi M, et al. (2004) The dynamics of the Sicily Strait: a comprehensive study from observations and models. Deep-Sea Res PT II 51: 411–440.
  28. 28. Robinson AR, Leslie WG, Theocharis A, Lascaratos A (2001) Ocean circulation currents: Mediterranean Sea circulation. In: Turekian KK, Thorpe SA, editors. Encyclopedia of Ocean Sciences. London: Academic Press. 1689–1706.
  29. 29. Bianchi CN (2007) Biodiversity issues for the forthcoming tropical Mediterranean Sea. Hydrobiologia 580: 7–21.
  30. 30. Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, et al. (2007) Marine ecoregions of the World: a bioregionalization of coastal and shelf areas. Bioscience 57: 573–583.
  31. 31. Magoulas A, Tsimenides N, Zouros E (1996) Mitochondrial DNA phylogeny and the reconstruction of the population history of a species: the case of the European anchovy (Engraulis encrasicolus). Mol Biol Evol 13: 178–190.
  32. 32. Viñas J, Alvarado Bremer J, Pla C (2004) Phylogeography of the Atlantic bonito (Sarda sarda) in the northern Mediterranean: the combined effects of historical vicariance, population expansion, secondary invasion, and isolation by distance. Mol Phylogenet Evol 33: 32–42.
  33. 33. Rabaoui L, Tlig-Zouari S, Ben Hassine OK (2008) Distribution and habitat of the fan mussel Pinna nobilis Linnaeus, 1758 (Mollusca: Bivalvia) along the northern and eastern Tunisian coasts. Cah Biol Mar 49: 67–78.
  34. 34. Galinou-Mitsoudi S, Vlahavas G, Papoutsi O (2006) Population study of the protected bivalve Pinna nobilis (Linnaeus, 1758) in Thermaikos Gulf (North Aegean Sea). J Molluscan Stud 5: 47–53.
  35. 35. Katsares V, Tsiora A, Galinou-Mitsoudi S, Imsiridou A (2008) Genetic structure of the endangered species Pinna nobilis (Mollusca: Bivalvia) inferred from mtDNA sequences. Biologia 63: 412–417.
  36. 36. Rabaoui L, Mejri R, Tlig-Zouari S, Bahri L, Hassine OKB (2011) Genetic variation among populations of the endangered fan mussel Pinna nobilis (Mollusca: Bivalvia) along the Tunisian coastline. Hydrobiologia 678: 99–111.
  37. 37. Frankham R (1995) Effective population size/adult population size ratios in wild life review. Genet Res 66: 95–107.
  38. 38. Escudero A, Iriondo MJ, Torres ME (2003) Spatial analysis of genetic diversity as a tool for plant conservation. Biol Conserv 113: 351–365.
  39. 39. Rio MH, Poulain PM, Pascual A, Mauri E, Larnicol G, et al. (2007) A mean dynamic topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model. J Mar Syst 65: 484–508.
  40. 40. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.
  41. 41. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452.
  42. 42. Swofford DL (2003) PAUP* Phylogenetic Analysis Using Parsimony (*and other methods), Version 4. Sinauer Associates editors. Sunderland Massachusetts.
  43. 43. Corander J, Tang J (2007) Bayesian analysis of population structure based on linked molecular information. Math Biosci 205: 19–31.
  44. 44. Pearse DE, Crandall KA (2004) Beyond FST: analysis of population genetic data for conservation. Conserv Genet 5: 585–602.
  45. 45. Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48.
  46. 46. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567.
  47. 47. Tamura K, Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: 512–526.
  48. 48. Posada D (2008) jModelTest: phylogenetic model averaging. Mol Biol Evol 25: 1253–1256.
  49. 49. Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction Data. Genetics 131: 479–491.
  50. 50. Benjamini Y, Hochberg Y (1995) Controlling the false dicovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B-Stat Methodol. 57: 289–300.
  51. 51. Rogers AR, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9: 552–569.
  52. 52. Schneider S, Excoffier L (1999) Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics 152: 1079–1089.
  53. 53. Rogers AR (1995) Genetic evidence for a Pleistocene population explosion. Evolution 49: 608–615.
  54. 54. Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA Polymorphism. Genetics 123: 585–589.
  55. 55. Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and backgroud selection. Genetics 147: 915–925.
  56. 56. Hein J, Schierup MH, Wiuf C (2005) Gene genealogies, variation and evolution. Oxford: University Press. 276 p.
  57. 57. Depaulis F, Mousset S, Veuille M (2003) Power of neutrality tests to detect bottlenecks and hitchhiking. J Mol Evol 57: 190–200.
  58. 58. Soriano AR, Ramos-Onsins SE, Rozas J, Calafell F, Navarro A (2008) Statistical power analysis of neutrality tests under demographic expansions, contractions and bottlenecks with recombination. Genetics 179: 555–567.
  59. 59. Luttikhuizen PC, Drent J, Baker AJ (2003) Disjunct distribution of highly diverged mitochondrial lineage clade and population subdivision in a marine bivalve with pelagic larval dispersal. Mol Ecol 12: 2215–2229.
  60. 60. Drummond AJ, Suchard MA, Xie D, Rambaut A (2012) Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 29: 1969–1973.
  61. 61. Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22: 1185–1192.
  62. 62. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A (2006) Relaxed phylogenetics and dating with confidence. PLoS Biol 4: e88.
  63. 63. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology 7: 214.
  64. 64. Casu M, Sanna D, Cossu P, Lai T, Francalacci P, et al. (2011) Molecular phylogeography of the microturbellarian Monocelis lineata (Platyhelminthes: Proseriata) in the North-East Atlantic. Biol J Linnean Soc 103: 117–135.
  65. 65. Borrero-Pérez GH, González-Wangüemert M, Marcos C, Pérez-Ruzafa A (2011) Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern. Mol Ecol 20: 1964–1975.
  66. 66. Dailianis T, Tsigenopoulos CS, Dounas C, Voultsiadou E (2011) Genetic diversity of the imperilled bath sponge Spongia officinalis Linnaeus, 1759 across the Mediterranean Sea: patterns of population differentiation and implications for taxonomy and conservation. Mol Ecol 20: 3757–3772.
  67. 67. Mokhtar-Jamaï K, Pascual M, Ledoux JB, Coma R, Féral JP, et al. (2011) From global to local genetic structuring in the red gorgonian Paramuricea clavata: the interplay between oceanographic conditions and limited larval dispersal. Mol Ecol 20: 3291–3305.
  68. 68. Varney RL, Galindo-Sanchez CE, Cruz P, Gaffney PM (2009) Population genetics of the eastern oyster Crassostrea virginica (Gmelin, 1791) in the Gulf of Mexico. J Shellfish Res 28: 855–864.
  69. 69. Xiao J, Cordes JF, Wang H, Guo X, Reece KS (2010) Population genetics of Crassostrea ariakensis in Asia inferred from microsatellite markers. Mar Biol 157: 1767–1781.
  70. 70. Lazoski C, Gusmão J, Boudry P, Solé-Cava AM (2011) Phylogeny and phylogeography of Atlantic oyster species: evolutionary history, limited genetic connectivity and isolation by distance. Mar Ecol-Prog Ser 426: 197–212.
  71. 71. Tarnowska K, Chenuil A, Nikula R, Féral J-P, Wolowicz M (2010) Complex genetic population structure of the bivalve Cerastoderma glaucum in a highly fragmented lagoon habitat. Mar Ecol-Prog Ser 406: 173–184.
  72. 72. Fernández MV, Heras S, Maltagliati F, Roldán MI (2013) Deep genetic divergence in giant red shrimp Aristaeomorpha foliacea (Risso, 1827) across a wide distributional range. J Sea Res. In press.
  73. 73. Bianchi CN, Morri C (2000) Marine Biodiversity of the Mediterranean Sea: situation, problems and prospects for future research. Mar Pollut Bull 40: 367–376.
  74. 74. Calvo M, Templado J, Oliverio M, Machordom A (2009) Hidden Mediterranean biodiversity: molecular evidence for a cryptic species complex within the reef building vermetid gastropod Dendropoma petraeum (Mollusca: Caenogastropoda). Biol J Linnean Soc 96: 898–912.
  75. 75. Pérès JM, Picard J (1964) Nouveau manuel de bionomie benthique de la mer Méditerranée. Recl Trav Stn mar Endoume 31: 1–137.
  76. 76. Sá-Pinto A, Baird SJE, Pinho C, Alexandrino P, Branco M (2010) A three-way contact zone between forms of Patella rustica (Mollusca: Patellidae) in the central Mediterranean Sea. Biol J Linnean Soc 100: 154–169.
  77. 77. Arnaud-Haond S, Migliaccio M, Diaz-Almela E, Teixeira S, Van De Vliet MS, et al. (2007) Vicariance patterns in the Mediterranean Sea: East-West cleavage and low dispersal in the endemic seagrass Posidonia oceanica. J Biogeogr 34: 963–976.
  78. 78. Pérez-Losada M, Nolte MJ, Crandall KA, Shaw PW (2007) Testing hypotheses of population structuring in the Northeast Atlantic Ocean and Mediterranean Sea using the common cuttlefish Sepia officinalis. Mol Ecol 16: 2667–2679.
  79. 79. Maltagliati F, Di Giuseppe G, Barbieri M, Castelli A, Dini F (2010) Phylogeography and genetic structure of the edible sea urchin Paracentrotus lividus (Echinodermata: Echinoidea) inferred from the mitochondrial cytochrome b gene. Biol J Linnean Soc 100: 910–923.
  80. 80. Nikula R, Väinölä R (2003) Phylogeography of Cerastoderma glaucum (Bivalvia: Cardiidae) across Europe: a major break in the Eastern Mediterranean. Mar Biol 143: 339–350.
  81. 81. Domingues VS, Bucciarelli G, Almada VC, Bernardi G (2005) Historical colonization and demography of the Mediterranean damselfish, Chromis chromis. Mol Ecol 14: 4051–4063.
  82. 82. Peijnenburg KTCA, Fauvelot C, Breeuwer AJ, Menken BJ (2006) Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers. Mol Ecol 15: 3319–3338.
  83. 83. Zulliger DE, Tanner S, Ruch M, Ribi G (2009) Genetic structure of the high dispersal Atlanto-Mediterranean sea star Astropecten aranciacus revealed by mitochondrial DNA sequences and microsatellite loci. Mar Biol 156: 597–610.
  84. 84. Carlton JT, Hodder J (1995) Biogeography and dispersal of coastal marine organisms: experimental studies on a replica of a 16th-century sailing vessel. Mar Biol 121: 721–730.
  85. 85. Johnson LE, Carlton JT (1996) Post-establishment spread in large-scale invasions: dispersal mechanisms of the zebra mussel Dreissena polymorpha. Ecology 77: 1686–1690.
  86. 86. Apte S, Holland BS, Godwin LS, Gardner JPA (2000) Jumping ship: a stepping stone event mediating transfer of nonindigenous species via a potentially unsuitable environment. Biol Invasions 2: 75–79.
  87. 87. Gollasch S (2007) Is ballast water a major dispersal mechanism for marine organisms? In: Nentwig W, editors. Biological Invasion, vol. 193. Berlin: Springer. 49–57.
  88. 88. Svitoch AA, Selivanov AO, Yanina TA (2000) The Pont-Caspian and Mediterranean basins in the Pleistocene (paleogeography and correlation). Oceanology 40: 868–881.
  89. 89. Larmuseau MHD, Van Houdt JKJ, Guelinckx J, Hellemans B, Volckaert FAM (2009) Distributional and demographic consequences of Pleistocene climate fluctuations for a marine demersal fish in the north-eastern Atlantic. J Biogeogr 36: 1138–1151.
  90. 90. Ho SYW, Phillips MJ, Cooper A, Drummond AJ (2005) Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol 22: 1561–1568.
  91. 91. Ho SYW (2007) Calibrating molecular estimates of substitution rates and divergence times in birds. J Avian Biol 38: 409–414.
  92. 92. Emerson BC (2007) Alarm bells for the molecular clock? No support for Ho et al.'s model of time-dependent molecular rate estimates. Syst Biol 56: 337–345.
  93. 93. Frey JK (1993) Modes of peripheral isolate formation and speciation. Syst Biol 42: 373–381.