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

Open Access

Peer-reviewed

Research Article

Octopus vulgaris (Cuvier, 1797) in the Mediterranean Sea: Genetic Diversity and Population Structure

Abstract

The common octopus, Octopus vulgaris Cuvier 1797, is a largely exploited cephalopod species in the Mediterranean Sea and the Atlantic Ocean, as well as along the coasts of Africa, Brazil and Japan, where its taxonomic identity is still debated. The assessment of its genetic structure is a pressing need to correctly manage the resource and to avoid overfishing and collapsing of local stocks. Here we analysed genetic variation and population structure of O. vulgaris using thirteen microsatellite loci in seven sampling localities from the Mediterranean Sea and one from the Atlantic Ocean. We also used a DNA barcoding approach by COI gene fragment to understand the phylogenetic relationships among the specimens here investigated and the ones whose sequences are available in literature. Our results reveal high levels of allelic richness and moderate heterozygosity in all samples investigated, and a pronounced differentiation of the Atlantic and Sicilian specimens. This latter aspect seems to support the isolation of the biota within the Strait of Messina. A certain degree of differentiation was detected among the other geographic samples within the Mediterranean Sea, which is more compatible with an island model than isolation by distance. The occurrence of null alleles affected more genetic diversity indices than population structure estimations. This study provides new insights about the genetic diversity and structure of O. vulgaris in the area of interest, which can be used as guidelines for a fisheries management perspective.

Citation: De Luca D, Catanese G, Procaccini G, Fiorito G (2016) Octopus vulgaris (Cuvier, 1797) in the Mediterranean Sea: Genetic Diversity and Population Structure. PLoS ONE 11(2): e0149496. https://doi.org/10.1371/journal.pone.0149496

Editor: Lorenzo Zane, University of Padova, ITALY

Received: July 5, 2015; Accepted: February 2, 2016; Published: February 16, 2016

Copyright: © 2016 De Luca 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: All relevant data are within the paper and its Supporting Information files.

Funding: Funding provided by SP2_WP4AZ2UO03, www.ritmare.it, Ritmare Italian Flagship Project, GP GF and www.szn.it, Stazione Zoologica Anton Dohrn di Napoli, GP GF. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

Introduction

Marine species show contrasting patterns of population structure due to the elevated number of intervening factors. The common view of marine populations as demographically open units has been questioned over the past decades [1] and evidence is accumulating for more complex scenarios [2]. Interactions among physical and biological factors, such as marine current patterns, sea bottom topology, and dispersal capability of the species at any biological stage (gametes, larvae, juveniles, adults) can account for different levels and patterns of gene flow, requiring ad hoc assessments for each specific case [3–6].

Among invertebrates, molluscs include sessile, slow walking or fast swimming species characterised by internal or external fertilization and direct or indirect development, living in a wide variety of habitats such as open sea, rocky or sandy bottoms. Several patterns of genetic structure have been identified, from genetic homogeneity [7,8] to isolation by distance [9,10], and island models of isolation [11]. Due to their short life-cycles, variable temperature- [12–14] and food availability-dependent growth rate, cephalopod stocks are highly unstable, being susceptible to overfishing and sensitive to environmental changes [15].

The concepts of differentiation, population structure and stock boundaries are strictly linked; a stock is commonly defined as “a group of organisms whose demographic/genetic trajectory is largely independent from other such groups” [1]. Since many mollusc species also represent important fishery resources, assessment of genetic structure can give important insights in stock management [16,17]. The insufficient knowledge on population boundaries does not allow to design adequate strategies for stock management and for improvement of aquaculture techniques (e.g. [18]), which would be essential to preserve the species from local collapse or extinction.

In the last few years, the decline of other fishery resources has led to an increase of cephalopod exploitation, especially in Europe [19], making a better knowledge of cephalopod ecology and life history traits a priority.

Studies dealing with genetic structure of cephalopods have increased over years [20–26], but the available information in many geographic regions is still scarce. In the Mediterranean Sea and the nearby Atlantic Ocean (the Eastern portion interesting European and African coasts), for example, few studies have investigated the genetic structure of cuttlefishes and squids (e.g. [27,28]), while more contributions exist on the common octopus (Octopus vulgaris Cuvier, 1797) ([29–34]; reviewed in [35]). However, the use of different molecular markers in not overlapping sampling sites led to several uncertainties. All these studies, together with others from different geographic regions, agree in finding genetic structure at different scales [33], generally not related to geographic distance among samples [29,36], although a latitudinal gradient along the Atlantic coasts of the Iberian Peninsula and Canary Islands was found [34]. If population structure is more influenced by passive movement of paralarvae through marine currents, or by active migration of adults, is still unclear.

O. vulgaris lifecycle spans from 9 to 15 months [37,38], the 10–15% of which represented by the planktonic phase at paralarval stage [39] and sexual maturity reached in a few months [32]. Since adult dispersal seems to be limited to forage habits (around a 15 m diameter, [40]) and short inshore and offshore migrations due to spawning behaviour [41], the major contributions to gene flow is supposed to depend on currents-driven paralarval movements. However, also climate change can directly affect gene flow and consequently population structure at both juvenile and adult level, determining critic fluctuations in stock abundance, or indirectly, through currents intensity and direction changes [42–44]. In this work we use microsatellite markers [45,46] to assess genetic diversity and structure of the common octopus, O. vulgaris, in the Mediterranean Sea and the nearby Atlantic Ocean. We also use assignment test to assess connectivity among geographic samples and we calculate the effective population size, for fishery and management purposes. Furthermore, we use COI sequences as DNA barcode to confirm the taxonomic assignment of the specimens here investigated and to explore their phylogenetic relationships.

Materials and Methods

Sampling

All the analyses have been carried out using preserved or freshly dead specimens collected from local fishermen. No use of live animals has been required for this study and no specific permissions were needed for the sampling activities in all of the investigated areas because our species of interest is commercially harvested (not endangered nor protected) and it was caught in areas where fishing is allowed. A total of 193 individuals were collected in seven localities across the Mediterranean Sea and in one locality in the nearby Atlantic Ocean (Fig 1; see Table 1 for details).

thumbnail
Download:
Fig 1. Map of sampling localities.

Eight geographic samples were collected across the Mediterranean Sea and the Atlantic Ocean. See Table 1 for geographic coordinates and other information.

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

thumbnail
Download:
Table 1. Characteristics of sampling sites, including geographic information and genetic diversity indices per geographic sample.

N: number of individuals sampled in each site; Npa: number of private alleles; Ar: allelic richness per geographic sample; He and Ho indicate expected and observed heterozygosity respectively; FIS is the inbreeding coefficient calculated in each population sample. In italics values after correction for null alleles.

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

According to [47], we sampled whenever possible about 25 individuals per sampling site (assumed as geographic samples). Samples collected from the Gulf of Naples (central Tyrrhenian Sea, Italy) in two geographical locations distant less than 7 km, were included to test for fine-scale population structuring due to differences in the sea-bottom landscape [48,49]. Adult octopuses were caught by local fishermen not far from the coast. For each individual, a piece of arm (2–3 cm) was excised and stored in 95–100% ethanol until DNA extraction.

DNA extraction

DNA was extracted using the NucleoSpin® Tissue kit (Macherey-Nagel). Starting from the tissue stored in alcohol, about 25 mg were cut, freed of outer skin and put in distilled water for a while to account for rehydration and alcohol removal. Then, the tissue was chopped in small pieces, put in 1.5 mL tubes and processed following the manufacturer’s instructions. Resulting DNA was quantified by NanoDrop Spectrophotometer (ND-1000, Thermo Scientific) and diluted to the optimal concentration for each downstream application.

DNA Barcoding with COI fragment

Despite O. vulgaris is easily recognisable from other octopus species of the Mediterranean Sea and the nearby Atlantic Ocean, a DNA barcode with a fragment of the COI gene was obtained for five individuals from each geographic sampling locality, using the primers developed by [50]. DNA amplification was carried out at the condition specified in [50]. PCR products were run on a 1.5% agarose gel, extracted in slices from the gel and purified using the GenElute™ Gel Extraction Kit (Sigma-Aldrich). DNA fragments were sequenced using LCO1490 primer [50] by Big-Dye Terminator Cycle sequencing kit (Applied Biosystems). The sequences obtained were compared with the ones available on public repositories such as GenBank and BOLD (Barcode Of Life Database) for taxonomic identification, and then aligned using the ClustalW algorithm [51] implemented in BioEdit v7.2.5 (http://www.mbio.ncsu.edu/bioedit/bioedit.html), together with other sequences of O. vulgaris with known geographical provenance [31, 52–60] (S1 Dataset). This dataset was used to build a Maximum Likelihood tree in MEGA6 [61] at the following parameters: GTR + G + I model with 4 gamma categories, 1000 bootstrap replications. The COI sequence of O. insularis was used as outgroup, due to its phylogenetic closeness to O. vulgaris [59].

Microsatellite amplification

Microsatellite amplification was carried out combining 13 loci (previously tested in simplex, [45,46]) in two multiplex PCR reactions (Table 2), using the Type-it® Microsatellite PCR kit (Qiagen) at the following conditions: 15 min at 95°C, 35 cycles at 94°C for 30 s, 59°C for 1 min and 30 s, 72°C for 1 min, and a final extension at 60°C for 30 min. Each reaction was conducted in a final volume of 12.5 μL. Forward primers were 5’ labelled with different fluorescent dyes (Life Technologies, see Table 2). PCR products were run on the Automated Capillary Electrophoresis Sequencer 3730 DNA Analyzer (Applied Biosystems) after a 1:300 dilution.

thumbnail
Download:
Table 2. Description of microsatellite loci utilised in this work.

Abbreviations: Na (number of alleles per locus), Ar (allelic richness), He (expected heterozygosity), Ho (observed heterozygosity). In italics values corrected for null alleles. The column Multiplex refers to the combinations of loci for PCR amplification.

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

Microsatellite analysis and genetic diversity

Microsatellite peaks were analysed with Peak Scanner™ Software v1.0 (Life Technologies) and scored as PCR products; in case of ambiguity in assigning sizes, samples were run at different dilutions or re-amplified. Micro-Checker v2.2.3 [62] was used to detect typographic errors, whilst MicroDrop v1.01 [63] to estimate and correct for allelic dropout rate. The software MicroDrop first calculates the correlation between the amount of missing data and homozygotes among individuals and loci, then fits the results of this analysis into a model and finally creates a new dataset drawing genotypes according to the posterior probabilities of the model. The resulting dataset was used jointly with the original one to evaluate the impact of null alleles on genetic diversity and population structure estimations.

Putative loci under selection were detected with LOSITAN [64], a selection detection workbench based on the FST-outlier method. The options “Neutral mean FST” and “Force mean FST” were used, to increase the reliability of the estimate of the initial mean FST computed from the empirical dataset. A total of 50,000 simulations were computed, plus 5,000 additional, as suggested by the default settings of the software. To test the power of our microsatellite markers on discriminating individuals using multilocus genotypes, we used the software GENCLONE v2.0 [65]. Number of alleles (Na) per locus and inbreeding coefficient (FIS) per population were calculated in Arlequin v3.5.1.3 [66], while FIS per locus [67] was calculated in Genepop on the Web [68,69]. Private alleles were detected using GenAlEx v6.501 [70,71]. Expected (He) and observed heterozygosity (Ho) per locus and population were calculated using GENETIX v4.05.2 [72]. Linkage disequilibrium was tested for each pair of loci in each population using Genepop on the Web [68,69]; the statistical significance of the pair-wise comparisons was determined using the Bonferroni correction for multiple tests.

Hardy-Weinberg equilibrium (HWE) probabilities per locus and population were computed in Genepop on the Web [68,69] both on the original dataset and the one corrected for null alleles; parameters utilised: 10,000 dememorisations, 100 batches and 10,000 iterations per batch. Allelic richness per population, based on minimum sample size (17 individuals) was calculated using standArich v1.00 [73], a package running in R [74].

Relationships among individuals within each sampling locality were calculated using ML-Relate [75] both on the original dataset and the one corrected for null alleles. A two tailed t-test was applied to verify if the results obtained without correcting for null alleles were significantly different from those calculated using correction.

Genetic differentiation and population structure

Genetic differentiation among geographic samples was calculated using Weir and Cockerham’s FST [67] as implemented in Arlequin v3.5.1.3 [66], and Jost’s D [76], using the package DEMEtics [77] implemented in R [74]. In general, FST measures deviation from a condition of panmixia and is calculated on allele frequencies, whilst D measures deviation from total differentiation and is more related to genetic distances among populations [78]. We set up the parameters for Jost’s D as follows: pair-wise comparisons and p-values as statistics to test against the null hypothesis of no genetic differentiation (based on 1,000 bootstrap resamplings). Bonferroni correction for multiple tests was performed for both FST and Jost’s D estimators. To test if the observed patterns of genetic structure conformed to the isolation by distance model, we correlated population pair-wise’s FST previously calculated with the geographic distance (km) between each pair of sampling locality, using the Mantel test function implemented in XLSTAT (Addinsoft, trial version; as add-ins running in Microsoft Excel). Geographical distances were calculated using Google Earth v7.0 (https://www.google.com/earth/), drawing a curve line moving along the coasts according to the circulation of water masses in the Mediterranean Sea [79,80]. We decided to follow this criterion to calculate geographic distances since seasonal adult octopus migrations, where studied, resulted to be very limited [38] and few data about paralarvae migrations are available [81].

After preliminary results, we decided to perform a discriminant analysis of principal components (DAPC) using the R working package ADEGENET v2.0.0 [82] to cluster our samples in groups, due to the peculiarity of multivariate methods of working in a context free of assumptions of Hardy-Weinberg equilibrium and linkage disequilibrium. We tested the reliability of our results using the α score function.

We also used GENECLASS2 [83] to assign individuals to sampling localities into a framework that does not require populations being in HWE, which best fits our data. We used the frequencies-based method [84], with the default frequency for missing alleles (0.01) as suggested by [85], and enabling computation probability under the simulation algorithm described in [85] (parameters applied: 1,000,000 simulated individuals, type I error set at 0.01).

Contemporary effective population size (Ne) from multilocus diploid genotypes was estimated using both the method based on linkage disequilibrium data [86] and the one based on molecular co-ancestry [87], as implemented in the software NeEstimator v2.01 [88]. Both methods calculate the number of effective breeders (Neb) of a cohort from which a sample is obtained that, in populations with non-overlapping generations (such as ours, with rare exceptions), is nearly equal to Ne. For the linkage disequilibrium method, a value of 0.020 was set as the lowest allele frequency used. We used preliminary results of population structure analyses to pool together homogeneous samples before calculating Ne, identifying as homogeneous only the geographic samples that were significantly different from all the others and pooling together the remaining.

Results

DNA Barcoding

All subsamples analysed were unambiguously attributed to O. vulgaris. The 482 bp DNA sequences of the COI region considered in the analysis were identical among all samples from Mediterranean Sea and nearby Atlantic Ocean collected in the present study. Three main clades are present in the ML phylogenetic tree obtained (S1 Fig). The first clade (93% bootstrap) harbours most of the Mediterranean samples and few samples from the nearby Atlantic Ocean (Portugal and Spain), Western and Southern Africa, the Southern Atlantic Ocean (Tristan da Cunha), and the Southern Indian Ocean (Amsterdam and Saint Paul Islands), a second clade (100% bootstrap) harbours samples from South-East Asia and a third one includes samples from Brasil. Samples from Turkey are present both in the Mediterranean and in the Asian clades.

Genetic diversity

All microsatellite loci were successfully amplified in our specimens, except for the locus oct03 in two individuals from the Gulf of Oristano (Italy; Mediterranean Sea) and one from the Balearic Sea (Spain) which also did not result to be successful at locus oct08. We decided to treat these as missing data and included the samples in subsequent analyses. Genotypes for each individual are provided in S2 Dataset. Only locus Vulg10 was found to be a putative “borderline” candidate under positive selection (S2 Fig); after checking that the results of downstream analyses were not affected by the presence of this locus, it was maintained in the final dataset.

A total of 193 unique multilocus genotypes were found over 193 individuals analysed, even after removing loci with missing data.

No correlation (Pearson’s) between amount of homozygotes and amount of missing data was found neither across individuals (r = 0.024, p = 0.36) nor across loci (r = -0.072, p = 0.57), excluding the occurrence of allelic dropout in the dataset. On the contrary, null alleles were detected at all loci, the most at very low frequencies (< 1%), and few around 10–15% (S1 Table). The impact of null alleles is evident in the significant deviations from Hardy-Weinberg equilibrium that were mostly reduced after correction. Indeed, 35 pair-wise locus comparisons that were not in HWE were recovered (Table 3), with locus Vulg06 being in HWE in all of the investigated geographic samples (Table 3). However, null alleles seem unlikely to be the only factor responsible for such pattern of disequilibrium, that otherwise should have been restored after correction.

thumbnail
Download:
Table 3. P values of Hardy-Weinberg equilibrium.

Values in italics refer to the dataset corrected for null alleles. See Table 1 for abbreviations.

https://doi.org/10.1371/journal.pone.0149496.t003

Linkage disequilibrium between loci was only found in two population samples after Bonferroni correction: Croatia, with two couples of loci in linkage (Vulg06 –Ov12; oct08 –Ov12), and Strait of Messina, where each locus was in linkage with at least another one (S2 Table). The high number of loci in linkage disequilibrium in Sicilian specimens is likely due to the occurrence of population structure and inbred individuals (see below) in this sample. Due to this “population” specific pattern, we decided to keep all loci for further analyses.

Microsatellite loci used in our dataset appeared to be moderately to highly polymorphic, with a number of alleles per locus ranging from 6 (Vulg10) to 32 (oct08) and observed heterozygosity values spanning from 0.187 to 0.731 (Table 2). However, it seems that Ho values per locus are affected by the occurrence of null alleles, as attested by the increasing of values after correction (from 0.202 to 0.871, Table 2).

Fifty-five private alleles were detected (Table 1), the most belonging to the Atlantic grouping (23), with a frequency ranging from 0.019 to 0.176 (S3 Table). Genetic diversity indices were quite homogeneous among all geographic samples (Table 1). Observed heterozygosity values per population were always slightly lower than expected, indicating excess of homozygosity, which correlate with the significant and comparable positive values of FIS per locus (Table 4) and geographic samples (Table 1). Significant differences were found in Ho estimates per locus and sample comparing the dataset corrected and not corrected for null alleles (p < 0.01, two-tailed Wilcoxon test).

thumbnail
Download:
Table 4. FIS estimations per locus with and without correction for null alleles.

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

Similarly to the case of HWE, inbreeding estimates corrected for null alleles were still positive (average 0.100 per locus and 0.062 per geographic sample), although lower than those calculated with the original dataset (average 0.275 per locus and 0.244 per geographic sample).

Allelic richness (Ar) per population ranged from 5.615 to 8.115, with the lowest value observed in the samples from the Strait of Messina (Table 1).

Related individuals using ML-Relate were found in all samples (S4 Table). The highest frequencies of full-sibs (FS) pairs were present in the geographic sample from the Strait of Messina (0.22059), whilst comparable frequencies of half-sibs (HS) were detected in all the other samples (S4 Table). The high frequency of HS in the geographic sample of Ria Formosa was reduced after correction for null alleles (S4 Table). Relationship pairs calculated with and without accounting for null alleles were not significantly different (p > 0.50), indicating that these results are not affected by the non-amplification of alleles.

Population differentiation and genetic structure

The estimation of genetic differentiation between geographic samples gave slightly different results according to the coefficient used, with Jost’s D being more conservative than Weir and Cockerham’s FST. Most of the statistically significant pair-wise comparisons were shared between the two coefficients, with the exception of the comparison Porto Cesareo—Spain (Table 5). Samples showing the highest differentiation, with statistically significant values for both indices in all pair-wise comparison with all the others, were Ria Formosa (PTG) and Strait of Messina (STM). Conversely, not significant genetic differentiation was found with both methods in the pair-wise comparison of the two geographic samples from the Gulf of Naples (FST = -0.00390, p > 0.05; Dest = 0.00177, p > 0.05) and between the Adriatic and Ionian specimens, Porto Cesareo and Croatia (FST = 0.00540, p > 0.05; Dest = 0.03085, p > 0.05) (Table 5). Significant differentiation was also found between the Adriatic specimens from Croatia (eastern Mediterranean Sea) and the ones from the Balearic Sea (Spain, western Mediterranean Sea) after Bonferroni correction using both indices (Table 5). No significant differences were found among geographic samples from Spain (SPN), the Gulf of Naples (NA1 and NA2), the Gulf of Oristano and Croatia with both indices (Table 5).

thumbnail
Download:
Table 5. Matrix of pair-wise FST (below the diagonal) and pair-wise Jost’s D (above the diagonal).

https://doi.org/10.1371/journal.pone.0149496.t005

Population differentiation estimations inferred using Weir and Cockerham’s FST did not change significantly using the dataset corrected for null alleles (S5 Table); on the contrary, only differentiation among Ria Formosa and Strait of Messina in respect to all the other sampling localities was retained using Jost’s D.

Despite several significant pair-wise comparisons, no isolation by distance was found using the Mantel test (R2 = 0.010, p = 0.623, see S3 Fig).

The DAPC plot confirmed a strong distinction between samples from Ria Formosa and the Strait of Messina, and an overlapping pattern for the other geographic specimens from the Mediterranean Sea (Fig 2). Loci Vulg15, Ov10 and Ov12 were identified as the ones mostly contributing to the observed variation. No differences in the distribution of individuals within clusters were detected using the α-score function. A similar output was obtained with the dataset corrected for null alleles (S4 Fig).

thumbnail
Download:
Fig 2. Discriminant analysis of principal components (DAPC).

Geographic samples are encompassed by ellipses. See Table 1 for abbreviations.

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

The assignment test performed in GENECLASS2 (Fig 3) allowed to correctly assign 37.8% of individuals (73 out of 193) to their sampling locality, the most belonging to the Atlantic and Sicilian groupings. Mediterranean specimens (except for some from the Strait of Messina, in red) have admixed genotypes and can only be distinguished on the different proportions of membership to each sampling locality: the two population samples from the Gulf of Naples show only slight differences between each other.

thumbnail
Download:
Fig 3. GENECLASS2 assignment test.

In the pie graphs, each value indicates the percentage of individuals correctly assigned to the source sampling locality. nas indicates individuals that were not unambiguously assigned. See Table 1 for abbreviations.

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

The analysis conducted with NeEstimator showed differences in the estimation of contemporary genetic effective population size (Ne) with the two methods used. Estimations with the LD method varied greatly considering or not the correction for null alleles, with overall wide confidence intervals (S6 Table). The method based on molecular co-ancestry showed that the samples with the lowest number of breeders are Ria Formosa (8) and Strait of Messina (6–8): these are also the ones with the highest frequency of related individuals and showing strong genetic structure (S6 Table). Furthermore, the confidence intervals obtained with this approach are narrower than the ones with the LD method, and no large differences were detected with and without considering null alleles (S6 Table).

Discussion

We used a population genetics approach by microsatellite markers to define population structure and genetic diversity of the common octopus Octopus vulgaris in the Mediterranean Sea and nearby Atlantic Ocean. Our results identify specimens from Ria Formosa (Atlantic Ocean) and the Strait of Messina (Mediterranean Sea) as distinct populations, as confirmed by the strong and significant differentiation detected with different methods of analysis. Conversely, the genetic differentiation among all the other specimens from the Mediterranean Sea is not straightforward, reflecting in the inability to correctly assign all the individuals to their locality of origin. Samples from the Adriatic and the Ionian Sea are homogeneous and no small-scale differentiation was detected in the Gulf of Naples. Our study also highlights higher levels of allelic richness but lower heterozygosity than other studies on the same species.

Genetic diversity

Significant deviations from Hardy-Weinberg equilibrium (HWE) were detected in all samples analysed. Correcting the dataset for null alleles, 35 locus by locus evaluations that were not in HWE were recovered, with a general reduction of the inbreeding coefficient FIS and two loci (Vulg06 and Vulg11) showing homozygotes deficit. Homozygosity excess due to null alleles is generally locus-specific [89], while in our analysis it was present across most of the loci (11 out of 13, considering null alleles correction), suggesting other factors being responsible for such pattern. The lowering of homozygosis after correction for null alleles is a clear signal of the impact of non-amplifying alleles, which hamper the estimates of genetic diversity.

Our FIS values per population corrected for null alleles agree with expectations, according to a review on 124 invertebrate species, where reproductive traits were related with FIS value [90]. Authors found that values of FIS < 0.1 are observed in marine invertebrates with internal copulation or sperm transfer such as cephalopods, whilst higher values (FIS > 0.3) in species with planktonic sperm (e.g. anthozoans, bivalves, gastropods, polychaetes and tunicates). In general, our geographic populations showed higher allelic richness but lower heterozygosity than other studies on the same species [33,34,36], while low values of heterozygosity were observed in other cephalopod species [27,91,92]. The moderate heterozygosity and positive FIS values found in our sampling sites could be related to non-random mating within geographic populations, as shown by kinship analysis. On the contrary, an excess of heterozygosis was previously reported for Brazilian populations of O. vulgaris [36], indicating that reproductive behaviour can be different in different geographic areas, although authors suggest a bottleneck event as possible cause of the heterozygosity excess among most of loci [36].

High allelic diversity and low to moderate heterozygosity are contrasting patterns and are difficult to explain. It is possible that the fortuitous migration of individuals among populations (especially within the Mediterranean Sea) is responsible for high allelic richness and low number of private alleles, while the biology of the species accounts for the degree of homozygosity recorded. Small and intense bottleneck events cause substantial loss of allelic diversity, while diffuse bottlenecks spread over generations provoke the same loss of heterozygosity but a smaller reduction of allelic diversity [89]. On the other hand, overfishing causes reduction in effective population size and consequent decrease of allelic diversity and heterozygosity [93,94]. In our study, effective population size (Ne) was particularly low in Atlantic and Sicilian localities, followed by the group encompassing octopuses from Sardinia, Spain and the Gulf of Naples.

Population structure

Our analysis of population structure in the Mediterranean Sea based on microsatellite markers agrees with the analysis conducted by Maltagliati and co-workers [29] using allozymes, in detecting sub-structuring among Mediterranean populations. The authors found populations from Sicily and Crete genetically distant from the ones in the western basin, with no clear structure among the other Mediterranean populations included in the analysis, that were all from the North-Western Mediterranean basin [29]. A further analysis conducted on a single microsatellite locus, which included the populations analysed by Maltagliati et al. [29], found significant structure within the whole Mediterranean basin [30].

Here we detected, in the Mediterranean, a strong isolation of the specimens from Strait of Messina, where the low values of Ho, Na and Ar are indicative of the occurrence of genetic drift. Furthermore, the high frequency of half-sibs detected in this sample strengthen its genetic structure. We are so confident in considering it as a distinct population. Moreover, we suggest the presence of a genetic subgrouping in the basin, supported by both differentiation indices used. Octopuses from the Adriatic and Ionian Sea are not genetically differentiated as well as the ones within the Gulf of Naples. We also found a partial isolation of the central and western basin, with the specimens from Croatia that are significantly different from the ones from Spain and the Gulf of Naples (NA2).

East and West Mediterranean populations appear separated in several marine species, as results of the glacial/post-glacial evolutionary history of the basin and present-day circulation patterns (e.g. [95] for seagrasses; [96] for fishes). We are not able to confirm this assumption; furthermore, according to our results, we can suppose that an island model of migration seems to be more likely to occur in the Mediterranean Sea for O. vulgaris, in respect to isolation by distance, as also suggested by [29].

As common for the most of marine species (e.g. [80]), we also found that the Atlantic samples (Ria Formosa) are well differentiated from the Mediterranean ones. In particular, in the Atlantic population, we found 23 private alleles, against the 2–7 private alleles found among the Mediterranean populations. It is likely that this skewed distribution of different alleles contributes to the observed differentiation. Assignment test was particularly effective in distinguishing Atlantic from Mediterranean specimens, but also highlighted the occurrence, within Mediterranean samples, of genotypes that were likely to belong to the Atlantic ones. Whether this pattern is the result of a recent migration mediated by ocean currents from the Atlantic Ocean to the Mediterranean Sea, or it is related to a past colonization from the Atlantic due to changes in sea level associated with cyclic ice ages during the Quaternary, is still to be assessed. The adult motility seems to be quite limited [32,38] and the dispersal of O. vulgaris paralarvae through marine currents is the most probable vehicle of long distance dispersal, although further investigations of the biology of the common octopus will help in discriminating between these scenarios. Mitochondrial data show differentiation between specimens from central and eastern Mediterranean [32], and within the Eastern basin [31], supporting the hypothesis that the roots of such genetic structure are historical. Our analysis showed that the COI gene is not a good mitochondrial marker to detect genetic structure, at least in the investigated area, contrasting with the findings of [31] along the coasts of Turkey. However, the ambiguous and multiple position of Turkish specimens along our phylogenetic tree highlighted the needing of further studies in that area.

The occurrence of genotypes that can only be assigned to the population from the Strait of Messina can be probably explained with the peculiar characteristics of the site. This is an area of strong currents generated by the opposition of tidal phases between the Tyrrhenian and Ionian basins [97]. These currents and vortices might isolate the biota of the Strait of Messina from the one of the Western basin or connect it with the Eastern basin through a unidirectional water mass flow. However, few studies are available in literature supporting this scenario [95]. It is now important to understand if this pattern is shared by other populations within the Strait and downstream, which would mean that Sicilian octopuses constitute a unique genetic unit in the whole basin, or if they are genetically closer to other specimens in areas of the Eastern Mediterranean. The genetic distinctiveness of the South Sicilian population of Porto Palo in [29], seems to support the first hypothesis.

Population genetic structure has been found in the common octopus also in other geographic areas at different spatial scales [31–34,36]. In some coastal areas, environmental parameters might give rise to barriers to gene flow, promoting differentiation. Along the coasts of Africa, strong seasonal upwelling phenomena are probably responsible for the differentiation of the Northwest African octopus fisheries [33,44]. On the contrary, along the southern coasts of Brazil, the structuring of 11 fishery sites into four genetic populations in absence of apparent oceanographic and geographic barriers is explained by bottleneck events [36]. In our study, we do not find significant differentiation at small scale (about 7 km; Gulf of Naples, Italy). The circulation of the water masses in the Gulf of Naples and relative migration patterns of the common octopus, may contribute to the observed homogeneity at small scale but not allowing connectivity at larger scales.

In the context of fishery management, particular attention should be paid on the octopus population of the Strait of Messina. The high number of loci in linkage disequilibrium, the small effective population size and the relatively low allelic richness are signs of genetic erosion in this population. Extensive fishing may have caused a reduction in effective population size through genetic drift and linkage disequilibrium among loci. Further analyses in this area using more individuals and extending the spatial scale could be of help to investigate the extent of such phenomenon.

The other homogenous groups from the Mediterranean Sea and the Atlantic Ocean derived by our analyses, show contrasting values of Ne according to the method used for estimation, making difficult to draw inferences.

The increasing trend in octopus catches in Europe but also in the world should focus the attention to this issue. According to FAO [98], European production of octopuses for 2010 was 42945 t, the most constituted by Octopus vulgaris, with some Eledone cirrhosa and E. moschata, from the western and central Mediterranean Sea and the Atlantic coasts of Iberian Peninsula. In particular, the five largest producers are Spain (16470 t), Portugal (10934 t), Italy (9884 t), Greece (2676 t), and France (1744 t). In Italy, 3295 t of O. vulgaris have been caught in 2012 (FAO FishStatJ, http://www.fao.org/fishery/statistics/software/fishstatj/en), showing a decline if compared with the previous years. In this scenario, understanding the genetic population structure of O. vulgaris as well as of other exploited species over time become fundamental in order to avoid a reduction of effective population size and consequently erosion of genetic diversity. We hope that the findings of our study will be the starting point for further investigations of genetic structure of the common octopus along its area of distribution.

Supporting Information

S1 Dataset. List of aligned COI sequences.

Files collected from public repositories are identified with the relative code in parenthesis.

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

(DOCX)

S2 Dataset. Microsatellite multilocus genotypes.

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

(XLSX)

S1 Fig. Maximum Likelihood tree based on COI sequences illustrating the evolutionary relationships among O. vulgaris specimens.

Values at the basis of each branch (node) indicate the percentage of occurrence of such specimens together with the other ones in the bootstrap test (1000 replicates).

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

(PDF)

S2 Fig. Outlier loci analysis conducted in LOSITAN.

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

(TIF)

S3 Fig. Mantel test.

Correlation between geographic distance (km) and pairwise genetic differentiation (FST).

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

(TIF)

S4 Fig. Discriminant analysis of principal components (DAPC) using the dataset corrected for null alleles.

Geographic samples are encompassed by ellipses. See Table 1 for abbreviations.

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

(TIF)

S1 Table. Null alleles rate per locus.

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

(DOCX)

S2 Table. Significant pair-wise linkage disequilibrium p-values.

Significance determined after Bonferroni adjustment at p < 0.000641.

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

(DOCX)

S3 Table. Private alleles and their frequency in each geographic sample.

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

(DOCX)

S4 Table. Relationships among individuals.

Abbreviations are as follows: U = unrelated, HS = half- sibs, FS = full-sibs, PO = parent/offspring. In italics are values corrected for null alleles.

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

(DOCX)

S5 Table. Matrix of pair-wise FST (below the diagonal) and Jost’s D (above the diagonal).

Asterisks indicate the p value of each pair-wise comparison: * p < 0.05; ** p < 0.01; *** p < 0.001; ns = not significant. In bold values significant after Bonferroni correction at p < 0.0017672. See Table 1 for abbreviations.

https://doi.org/10.1371/journal.pone.0149496.s011

(DOCX)

S6 Table. Contemporary effective population size (Ne) based on linkage disequilibrium and molecular co-ancestry models, with relative confidence intervals (CIs).

Values in italics refer to the dataset corrected for null alleles.

https://doi.org/10.1371/journal.pone.0149496.s012

(DOCX)

Acknowledgments

We thank Rui Santos (Portugal), Antonio Esteban (Spain), Ivan Guala, Giusy Scarfì, Antonio Terlizzi (Italy), and Ante Zuljevic (Croatia) for having contributed in sampling of octopus specimens; Marlene Jahnke for her precious suggestions on the analysis of the data, and the academic editor Lorenzo Zane and the three anonymous reviewers for the constructive comments on the manuscript. O. vulgaris specimens from Spain were collected during the campaign organised by the Spanish government MEDITS_ES_2014 (GSA 06) in the Western Mediterranean.

Author Contributions

Conceived and designed the experiments: DDL GC GP GF. Performed the experiments: DDL. Analyzed the data: DDL GC. Contributed reagents/materials/analysis tools: GP GF. Wrote the paper: DDL GC GP GF. Supervised the work: GP GF.

References

  1. 1. Waples R. Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. J Hered. 1998;89: 438–450.
  2. 2. Hellberg ME. Gene Flow and Isolation among Populations of Marine Animals. Annu Rev Ecol Evol S. 2009;40: 291–310.
  3. 3. Weersing K, Toonen RJ. Population genetics, larval dispersal, and connectivity in marine systems. Mar Ecol Prog Ser. 2009;393: 1–12.
  4. 4. Hoffman JI, Clarke A, Clark MS, Peck LS. Hierarchical Population Genetic Structure in a Direct Developing Antarctic Marine Invertebrate. PloS One. 2013;8: e63954. pmid:23691125
  5. 5. Villamor A, Costantini F, Abbiati M. Genetic Structuring across Marine Biogeographic Boundaries in Rocky Shore Invertebrates. PloS One. 2014;9: e101135. pmid:24983738
  6. 6. Haye PA, Segovia NI, Muñoz-Herrera NC, Gálvez FE, Martínez A, Meynard A et al. Phylogeographic Structure in Benthic Marine Invertebrates of the Southeast Pacific Coast of Chile with Differing Dispersal Potential. PloS One. 2014;9: e88613 pmid:24586356
  7. 7. Addison JA, Ort BS, Mesa KA, Pogson GH. Range-wide genetic homogeneity in the California sea mussel (Mytilus californianus): a comparison of allozymes, nuclear DNA markers, and mitochondrial DNA sequences. Mol Ecol. 2008;17: 4222–4232. pmid:19378402
  8. 8. Giantsis IA, Mucci N, Randi E, Abatzopoulos TJ, Apostolidis AP. Microsatellite variation of mussels (Mytilus galloprovincialis) in central and eastern Mediterranean: genetic panmixia in the Aegean and the Ionian Seas. J Mar Biol Assoc UK. 2014;94: 797–809.
  9. 9. Rose CG, Paynter KT, Hare MP. Isolation by Distance in the Eastern Oyster, Crassostrea virginica, in Chesapeake Bay. J Hered. 2006;97: 158–170. pmid:16489144
  10. 10. Kenchington EL, Patwary MU, Zouros E, Bird CJ. Genetic differentiation in relation to marine landscape in a broadcast-spawning bivalve mollusc (Placopecten magellanicus). Mol Ecol. 2006;15: 1781–1796. pmid:16689898
  11. 11. Hoffman JI, Clarke A, Clark MS, Fretwell P, Peck LS. Unexpected Fine-Scale Population Structure in a Broadcast-Spawning Antarctic Marine Mollusc. PloS One. 2012;7: e32415. pmid:22403655
  12. 12. Forsythe JW, Van Heukelem WF. Growth in cephalopod life cycles. Comparative Reviews. Boyle PR, editor. London: Academic Press; 1987. pp. 135–136.
    • 13. Forsythe JW. A working hypothesis on how seasonal temperature change may impact the field growth of young cephalopods. In: Recent Advances in Cephalopod Fisheries Biology Okutani T, O’Dor RK, Kubodera T, editors.Tokyo: Tokai University Press; 1993. pp. 133–143.
      • 14. Forsythe JW. Accounting for the effect of temperature on squid growth in nature: from hypothesis to practice. Mar Freshwater Res. 2004;55: 331–339.
      • 15. Pierce GJ, Valavanis VD, Guerra A, Jereb P, Orsi-Relini L, Bellido JM et al. A review of cephalopod–environment interactions in European Seas. Hydrobiologia. 2008;612: 49–70.
      • 16. Shaw PW, Arkhipkin AI, Adcock GJ, Burnett WJ, Carvalho GR, Scherbich JN et al. DNA markers indicate that distinct spawning cohorts and aggregations of Patagonian squid, Loligo gahi, do not represent genetically discrete subpopulations. Mar Biol. 2004;144: 961–970.
      • 17. Buresch KC, Gerlach C, Hanlon RT. Multiple genetic stocks of longfin squid Loligo pealeii in the NWAtlantic: stocks segregate inshore in summer, but aggregate offshore in winter. Mar Ecol Progr Ser. 2006;310: 263–270.
      • 18. Villanueva R, Sykes AV, Vidal EAG, Rosas C, Nabhitabhata J, Fuentes L et al. Current Status and Future Challenges in Cephalopod Culture. In: Iglesias J, Fuentes L, Villanueva R editors. Cephalopod Culture. Springer Netherlands; 2014. pp. 479–489.
        • 19. Xavier JC, Allcock AL, Cherel Y, Lipinski MR, Pierce GJ, Rodhouse PGK et al. Future challenges in cephalopod research. J Mar Biol Assoc UK. 2014;
        • 20. Reichow D, Smith RJ. Microsatellites reveal high levels of gene flow among populations of the California squid Loligo opalescens. Mol Ecol. 2001;10: 1101–1109. pmid:11380869
        • 21. Herke S, Foltz D. Phylogeography of two squid (Loligo pealei and L. plei) in the Gulf of Mexico and northwestern Atlantic Ocean. Mar Biol. 2002;140: 103–115.
        • 22. Ibáñez CM, Cubillos LA, Tafur R, Argüelles J, Yamashiro C, Poulin E. Genetic diversity and demographic history of Dosidicus gigas (Cephalopoda: Ommastrephidae) in the Humboldt Current System. Mar Ecol Progr Ser. 2011;431: 163–171.
        • 23. Ibáñez CM, Argüelles J, Yamashiro C, Adasme L, Cespedes R, Poulin E. Spatial genetic structure and demographic inference of the Patagonian squid Doryteuthis gahi in the south-eastern Pacific Ocean. J Mar Biol Assoc UK. 2012;92: 197–203.
        • 24. Higgins KL, Semmens JM, Doubleday ZA, Burridge CP. Life history matters: comparisons of population structuring in sympatric octopus species that differ in the presence of a pelagic larval stage. Mar Ecol Progr Ser. 2013;486: 203–212
        • 25. Li J, Ye Y, Wu C, Guo B, Gul Y. Genetic diversity and population structure of Sepiella japonica (Mollusca: Cephalopoda: Decapoda) inferred by mitochondrial DNA (COI) variations. Biochem Syst Ecol. 2014;56: 8–15
        • 26. Williams RC, Jackson BC, Duvaux L, Dawson DA, Burke T, Sinclair W. The genetic structure of Nautilus pompilius populations surrounding Australia and the Philippines. Mol Ecol. 2015;24: 3316–3328. pmid:26033519
        • 27. Garoia F, Guarniero I, Ramsak A, Ungaro N, Landi M, Piccinetti et al. Microsatellite DNA variation reveals high gene flow and panmictic populations in the Adriatic shared stocks of the European squid and cuttlefish (Cephalopoda). Heredity. 2004;93: 166–174. pmid:15150540
        • 28. Meriam T, Wafa T, Khawla T, Tarek H, Abdeljelil G, Mhamed E. Genetic diversity and population structure of Sepia officinalis from the Tunisian cost revealed by mitochondrial COI sequences. Mol Biol Rep. 2014.
        • 29. Maltagliati F, Belcari P, Casu D, Casu M, Sartor P, Vargiu G et al. Allozyme genetic variability and gene flow in Octopus vulgaris (Cephalopoda, Octopodidae) from the Mediterranean Sea. B Mar Sci. 2002;71: 473–486.
        • 30. Casu M, Maltagliati F, Meloni M, Casu D, Cossu P, Binelli G et al. Genetic structure of Octopus vulgaris (Mollusca, Cephalopoda) from the Mediterranean Sea as revealed by a microsatellite locus. Ital J Zool. 2002;69: 295–300.
        • 31. Keskin E, Atar HH. Genetic divergence of Octopus vulgaris species in the eastern Mediterranean. Biochem Syst Ecol. 2011;39: 277–282.
        • 32. Fadhlaoui-Zid K, Knittweis L, Aurelle D, Nafkha C, Ezzeddine S, Fiorentino F et al. Genetic structure of Octopus vulgaris (Cephalopoda, Octopodidae) in the central Mediterranean Sea inferred from the mitochondrial COIII gene. CR Biol. 2012;335: 625–636.
        • 33. Murphy JM, Balguerías E, Key LN, Boyle PR. Microsatellite DNA markers discriminate between two Octopus vulgaris (Cephalopoda: Octopoda) fisheries along the northwest African coast. Bull Mar Sci. 2002;71: 545–553.
        • 34. Cabranes C, Fernandez-Rueda P, Martínez JL. Genetic structure of Octopus vulgaris around the Iberian Peninsula and Canary Islands as indicated by microsatellite DNA variation. ICES J Mar Sci. 2008;65: 12–16.
        • 35. De Luca D, Catanese G, Procaccini G, Fiorito G. An integration of historical records and genetic data to the assessment of global distribution and population structure in Octopus vulgaris. Front Ecol Evol. 2014;
        • 36. Moreira AA, Tomás ARG, Hilsdorf AWS. Evidence for genetic differentiation of Octopus vulgaris (Mollusca, Cephalopoda) fishery populations from the southern coast of Brazil as revealed by microsatellites. J Exp Mar Biol Ecol. 2011;407: 34–40.
        • 37. Hernandez-Lopez JL, Castro-Hernandez JJ. Age determined from the daily deposition of concentric rings on common octopus (Octopus vulgaris) beaks. Fish Bull. 2001;99:679–684.
        • 38. Smale MJ, Buchan PR. Biology of Octopus vulgaris off the east coast of South Africa. Mar Biol. 1981;65: 1–12.
        • 39. Villanueva R, Nozais C, Boletzky SV. Swimming behaviour and food searching in planktonic Octopus vulgaris Cuvier from hatching to settlement. J Exp Mar Biol Ecol. 1996;208: 169–184.
        • 40. Mather JA. ‘Home’ choice and modification by juvenile Octopus vulgaris (Mollusca: Cephalopoda): specialized intelligence and tool use?. J Zool. 1994;233: 359–368.
        • 41. Mangold K. Octopus vulgaris. In: Boyle PR editor. Cephalopod Life Cycles. Species Accounts, vol. I. Academic Press: London; 1983. pp. 335–364.
          • 42. Caddy JF. The Cephalopods: Factors Relevant to their Populations Dynamics and to Assessment and Management of Stocks. Advances in Assessment of World Cephalopods Resources. FAO Fish Tech Paper. 1983;231: 416–452.
          • 43. Boyle PR, Boletzky SV. Cephalopod populations: definition and dynamics. Philos T Roy Soc B. 1996;351: 985–1002.
          • 44. Faure V, Inejih CA, Demarcq H, Cury P. The importance of retention process in upwelling areas for recruitment of Octopus vulgaris: the example of the Arguin Bank (Mauritania). Fish Oceanogr. 2000;9: 343–355.
          • 45. Greatorex EC, Jones CS, Murphy J, Key LN, Emery AM, Boyle PR. Microsatellite markers for investigating population structure in Octopus vulgaris (Mollusca: Cephalopoda). Mol Ecol. 2000;9: 641–642. pmid:10792712
          • 46. De Luca D, Catanese G, Fiorito G, Procaccini G. A new set of pure microsatellite loci in the common octopus Octopus vulgaris Cuvier, 1797 for multiplex PCR assay and their cross-amplification in O. maya Voss & Solís Ramírez, 1966. Conserv Genet Resour. 2015;7: 299–301.
          • 47. Hale ML, Burg TM, Steeves TE. Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PloS One. 2012;7: e45170 pmid:22984627
          • 48. Lo Bianco S. Notizie biologiche riguardanti specialmente il periodo di maturità sessuale degli animali del Golfo di Napoli. Mitt Zool Sta Neapel. 1909;19: 513–763.
          • 49. Ranzi S. La distribuzione della vita nel Golfo di Napoli. Atti del XI Congresso geografico italiano. 1930;2: 1–4.
          • 50. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3: 294–299. pmid:7881515
          • 51. Thompson JD, Higgins DG, Gibson TJ. ClustalW: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994;22: 4673–4680. pmid:7984417
          • 52. Takumiya M, Kobayashi M, Tsuneki K, Furuya H. Phylogenetic relationships among major species of Japanese coleoid cephalopods (Mollusca: Cephalopoda) using three mitochondrial DNA sequences. Zoolog Sci. 2005;22: 147–155. pmid:15738635
          • 53. Teske PR, Oosthuizen A, Papadopoulos I, Barker NP. Phylogeographic structure of Octopus vulgaris in South Africa revisited: identification of a second lineage near Durban harbour. Mar Biol. 2007;151: 2119–2122.
          • 54. Guerra Á, Roura Á, González ÁF, Pascual S, Cherel Y, Pérez-Losada M. Morphological and genetic evidence that Octopus vulgaris Cuvier, 1797 inhabits Amsterdam and Saint Paul Islands (southern Indian Ocean). ICES J Mar Sci. 2010;67: 1401–1407.
          • 55. Kaneko N, Kubodera T, Iguchis A. Taxonomic study of shallow-water octopuses (Cephalopoda: Octopodidae) in Japan and adjacent waters using mitochondrial genes with perspectives on Octopus DNA barcoding. Malacologia. 2011;54: 97–108.
          • 56. Dai L, Zheng X, Kong L, Li Q. DNA barcoding analysis of Coleoidea (Mollusca: Cephalopoda) from Chinese waters. Mol Ecol Resour. 2012;12: 437–447. pmid:22348331
          • 57. Keskin E, Atar HH. DNA barcoding commercially important aquatic invertebrates of Turkey. Mitochondrial DNA. 2013;24: 440–450. pmid:23387292
          • 58. Lü ZM, Cui WT, Liu LQ, Li HM, Wu CW. Phylogenetic relationships among Octopodidae species in coastal waters of China inferred from two mitochondrial DNA gene sequences. Genet Mol Res. 2013;12: 3755–65. pmid:24085437
          • 59. Sales JBDL, Rego PSD, Hilsdorf AWS, Moreira AA, Haimovici M, Tomás AR et al. Phylogeographical features of Octopus vulgaris and Octopus insularis in the Southeastern Atlantic Based on the Analysis of Mitochondrial Markers. J Shellfish Res. 2013;32: 325–339.
          • 60. Strugnell JM, Norman MD, Vecchione M, Guzik M, Allcock AL. The ink sac clouds octopod evolutionary history. Hydrobiologia. 2014;725: 215–235.
          • 61. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30: 2725–2729. pmid:24132122
          • 62. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4: 535–538.
          • 63. Wang C, Schroeder KB, Rosenberg NA. A maximum-likelihood method to correct for allelic dropout in microsatellite data with no replicate genotypes. Genetics. 2012;192: 651–669. pmid:22851645
          • 64. Antao T, Lopes A, Lopes RJ, Beja-Pereira A, Luikart G. LOSITAN: a workbench to detect molecular adaptation based on a Fst-outlier method. BMC bioinformatics. 2008;9: 323. pmid:18662398
          • 65. Arnaud-Haond S, Belkhir K. genclone: a computer program to analyse genotypic data, test for clonality and describe spatial clonal organization. Mol Ecol Notes. 2007;7: 15–17.
          • 66. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10: 564–567. pmid:21565059
          • 67. Weir BS, Cockerham CC. Estimating F-Statistics for the Analysis of Population Structure. Evolution. 1984;38: 1358–1370.
          • 68. Raymond M, Rousset F. GENEPOP (Version 1.2): Population Genetics Software for Exact Tests and Ecumenicism. J Hered. 1995;86: 248–249.
          • 69. Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8: 103–106. pmid:21585727
          • 70. Peakall ROD, Smouse PE. genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6: 288–295.
          • 71. Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28: 2537–2539. pmid:22820204
          • 72. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier (France). (1996–2004); http://kimura.univ-montp2.fr/genetix/.
            • 73. Alberto F. standArich_v1.00: an R package to estimate population allelic richness using standardized sample size. 2006; http://www.ccmar.ualg.pt/maree/software.php?soft=sarich.
              • 74. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2015; URL https://www.R-project.org/.
                • 75. Kalinowski ST, Wagner AP, Taper ML. ML-Relate: a computer program for maximum likelihood estimation of relatedness and relationship. Mol Ecol Notes. 2006;6: 576–579.
                • 76. Jost LOU. GST and its relatives do not measure differentiation. Mol Ecol. 2008;17: 4015–4026. pmid:19238703
                • 77. Gerlach G, Jueterbock A, Kraemer P, Deppermann J, Harmand P. Calculations of population differentiation based on GST and D: forget GST but not all of statistics! Mol Ecol. 2010;19: 3845–3852. pmid:20735737
                • 78. Whitlock MC. G'ST and D do not replace FST. Mol Ecol. 2011;20: 1083–1091. pmid:21375616
                • 79. Robinson AR, Leslie WG, Theocharis A, Lascaratos A. Mediterranean Sea Circulation. In: Ocean currents: a derivative of the Encyclopedia of Ocean Sciences. 2001. pp. 1689–1705.
                • 80. Patarnello T, Volckaert FA, Castilho R. Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break? Mol Ecol. 2007;16: 4426–4444. pmid:17908222
                • 81. Moreno A, Dos Santos A, Piatkowski U, Santos AMP, Cabral H. Distribution of cephalopod paralarvae in relation to the regional oceanography of the western Iberia. J Plankton Res. 2009;31: 73–91.
                • 82. Jombart T. Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24: 1403–1405. pmid:18397895
                • 83. Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A. GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered. 2004;95: 536–539. pmid:15475402
                • 84. Paetkau D, Calvert W, Stirling I, Strobeck C. Microsatellite analysis of population structure in Canadian polar bears. Mol Ecol. 1995;4: 347–354. pmid:7663752
                • 85. Paetkau D, Slade R, Burden M, Estoup A. Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol. 2004;13: 55–65. pmid:14653788
                • 86. Waples RS, Do C. Linkage disequilibrium estimates of contemporary Ne using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl. 2010;3: 244–262. pmid:25567922
                • 87. Nomura T. Estimation of effective number of breeders from molecular coancestry of single cohort sample. Evol Appl. 2008;1: 462–474. pmid:25567728
                • 88. Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol Ecol Resour. 2014;14: 209–214. pmid:23992227
                • 89. Allendorf FW, Luikart G, Aitken SN. Conservation and the genetics of populations. 2nd ed. Oxford: John Wiley & Sons; 2013.
                  • 90. Addison JA, Hart MW. Spawning, copulation and inbreeding coefficients in marine invertebrates. Biol Lett. 2005;1: 450–453. pmid:17148230
                  • 91. Shaw PW, Pérez-Losada M. Polymorphic microsatellites in the common cuttlefish Sepia officinalis (Cephalopoda). Mol Ecol. 2000;9: 237–238. pmid:10672169
                  • 92. Pérez-Losada M, Guerra A, Carvalho GR, Sanjuan A, Shaw PW. Extensive population subdivision of the cuttlefish Sepia officinalis (Mollusca: Cephalopoda) around the Iberian Peninsula indicated by microsatellite DNA variation. Heredity. 2002;89: 417–424. pmid:12466983
                  • 93. Hauser L, Adcock GJ, Smith PJ, Bernal Ramírez JH, Carvalho GR. Loss of microsatellite diversity and low effective population size in an overexploited population of New Zealand snapper (Pagrus auratus). P Natl Acad Sci USA. 2002;99: 11742–11747.
                  • 94. Pérez-Ruzafa Á, González-Wangüemert M, Lenfant P, Marcos C, García-Charton JA. Effects of fishing protection on the genetic structure of fish populations. Biol Conserv. 2006;129: 244–255.
                  • 95. Serra IA, Innocenti AM, Di Maida G, Calvo S, Migliaccio M, Zambianchi E et al. Genetic structure in the Mediterranean seagrass Posidonia oceanica: disentangling past vicariance events from contemporary patterns of gene flow. Mol Ecol. 2010;19: 557–568. pmid:20051010
                  • 96. Bahri-Sfar L, Lemaire C, Hassine OKB, Bonhomme F. Fragmentation of sea bass populations in the western and eastern Mediterranean as revealed by microsatellite polymorphism. Proc R Soc B. 2000;267: 929–935. pmid:10853737
                  • 97. Bignami F, Salusti E. Tidal Currents and Transient Phenomena in the Strait of Messina: a Review. In: Pratt LJ, editor. The Physical Oceanography of Sea Straits. Springer Netherlands; 1990 pp. 95–124.
                    • 98. Jereb P, Roper CFE, Norman MD, Finn JK. Cephalopods of the world. An annotated and illustrated catalogue of cephalopod species known to date. Volume 3. Octopods and Vampire Squids. Rome: FAO; 2013.