Phylogeography of the veined squid, Loligo forbesii, in European waters

Scientific Reports volume 12, Article number: 7817 (2022) Cite this article

1398 Accesses

4 Citations

7 Altmetric

Metrics details

The veined squid, Loligo forbesii Steenstrup, 1856, occurs at the European Shelf areas including the Azores and represents a valuable resource for the European commercial fishery in the North East Atlantic. However, very little is known about its population structure and phylogeography. This lack of knowledge also impedes the development of sustainable fishery management for this species. The present study combined the use of two types of markers that retrieve patterns of gene flow in different time spans; the analysis of 16 nuclear microsatellites and sequencing of the mitochondrial cytochrome oxidase subunit I (COI). Whereas the high mutation rate of microsatellites allows the description of recent patterns of connectivity in species, the lower mutation rate of COI provides phylogeographic patterns on a longer timescale. A total of 347 individuals of L. forbesii were investigated from nearly the entire distribution range of the species, including the North East Atlantic Shelf, the Azores and the Mediterranean. Individuals from the Western and Eastern Mediterranean Sea have never been included in a genetic study before. We were able to analyse COI sequences from all 12 sampling areas and define three clades of L. forbesii. Due to our large sampling area, we are presenting 13 COI-haplotypes that were previously unknown. The microsatellite analysis does not include the Azores but three main clades could be identified at the remaining 11 sampling sites. Low FST values indicate gene flow over large geographical distances. However, the genetically significant differences and an additional slight grouping in the microsatellite structure reveal that geographical barriers seem to influence the population structure and reduce gene flow. Furthermore, both markers provide strong evidence that the observed phylogeographic pattern reflects the geographical history of the Azores and the Mediterranean Sea.

Within the last decades biomass of various cephalopod populations and commercial catches increased world-wide1. Especially in the Atlantic Ocean, the commercial importance of the loliginids has grown2,3. Loliginids landings reached about 12,000 t3 in 2017 showing a fivefold increase between 2000 and 2017 in the entire Northeast Atlantic, which illustrates the growing socio-economic importance. Among the loliginids, the veined squid Loligo forbesii Steenstrup, 1856 is one of the most important species for the European fishery3. This species is neritic, associated to the shelf and equally distributed on the lower shelf (80–200 m) and the upper slope (200–500 m) in the northern Mediterranean4. In the North East Atlantic it occurs from the North Sea and UK waters to the Canary Islands and Azores2. This opportunistic predator is characterized by high growth rates. During its 12–16 months of life span, it can reach a dorsal mantle length of up to 900 mm5,6. As a semelparous species L. forbesii spawns only once in its lifetime2.

In most regions L. forbesii is unregulated fished as bycatch, but European target fisheries exist in the English Channel, Scottish and Portuguese waters3 and recreational fishery develops in UK and Norway. First prediction models for L. forbesii have been developed7, but up to now it has been difficult to develop a sustainable fishery management for this valuable resource due to knowledge gaps regarding population structure8, life cycle and spawning areas9. Here, we will focus on the first aspect and will shed new light on the genetic population structure of L. forbesii.

The mitochondrial cytochrome-c-oxidase I (COI) gene is one of the most commonly used markers for molecular systematics and is often used to infer phylogeographic patterns10. However, to elucidate the complex population structure of L. forbesii, allozymes and mitochondrial DNA sequences were found apparently unsuitable because of their extremely low levels of genetic variability11,12. In such cases, a more appropriate approach is the analysis of microsatellites. Microsatellites are short repetitive sequences in the nuclear genome showing high mutation rates and multiple alleles13, allowing a powerful analysis of recent patterns of population connectivity. Shaw et al.12 found a higher genetic variability by using microsatellites for L. forbesii compared with previous studies using allozyme and mitochondrial DNA (mtDNA) markers, and were able to identify subtle genetic differentiation between species occurring at the European shelf seas (Scotland to Northern Spain) and offshore population including the Azores, Rockall and Faroes. They suggest that water depth and isolating current regimes are responsible for the genetic differences between offshore and onshore areas.

To shed new light on the genetic structure and phylogeography of L. forbesii we performed a mitochondrial COI and nuclear microsatellite analysis using the same dataset. In comparison to Shaw et al.12 we were able to geographically expand the study area to almost the entire distribution range of L. forbesii including the North-East Atlantic (from the North Sea to the Gulf of Cadiz, including the Azores) and the Mediterranean Sea (from the Balearic to the Aegean Sea). The new findings in population structure, gene flow and habitat connectivity of L. forbesii might have a substantial influence on the management of this species’ fisheries, as they provide new information regarding the identification of different appropriate management units, so called "stocks"14.

Loligo forbesii was sampled at 12 European areas, including shelf and upper slope areas and the Azores in 2019 (Fig. 1), resulting in a total of 347 analysed individuals (see Supplementary Table S1 online). Squids were trawled during different research cruises and frozen on board for processing on land. Additional individuals were collected at commercial markets.

Sampled areas: A = Aegean Sea (19 individuals), B = Balearic Sea (28), C = Gulf of Cadiz (20), D = South Adriatic Sea (15), I = East Ionian Sea (30), K = Celtic Sea (19), L = English Channel (30), N = North Sea (64), O = east coast of Sardinia (30), S = Bay of Biscay (24), W = west coast of Sardinia (21), Z = Azores (47) (Ocean Data View version 4.6.2).

From each individual, a pea-sized piece of muscle tissue was cut off the mantle and transferred into 96% undenatured ethanol. Ethanol was changed after 2 and 24 h if possible, to extract all water from the tissue and preserve the DNA. Samples, excluding North Sea samples which were prepared at the Thünen Institute of Baltic Sea Fisheries (TI-OF), were sent to the TI-OF, alcohol was changed once more and stored until the beginning of the genetic analysis at University of Rostock, where DNA was extracted using "innuPrepMini Kit" (Analytik Jena, Germany) following the manufacturer's standard protocol.

Loligo forbesii is not protected under any legislation and not considered threatened or endangered. Atlantic samples were collected as bycatch during ICES coordinated international fishing trawl surveys (International Bottom Trawl Survey: North Sea, Bay of Biscay; Spanish Gulf of Cadiz Bottom Trawl Survey: Cadiz) and CEFAS Otter Trawl Survey (Celtic Sea) or at the fish market (English Channel). L. forbesii samples from the Mediterranean were bycatch from the International bottom trawl survey in the Mediterranean (MEDITS) except samples from the Balearic Sea which were collected from a commercial vessel. Azorean samples have been collected by the artisanal jigging fishery. All individuals were already dead when they were handed over to collect the tissue samples for our study and therefore handled in accordance with relevant guidelines and regulations. Fishing licenses for all surveys and trips were available.

Population genetics of 300 individuals from eleven sampling areas were performed via microsatellite analysis at the University of Rostock using 16 specific primers (Lfor1–Lfor16) for L. forbesii15,16. Specimens from the Azores were not included in this analysis, because (1) we received the samples only 1 year after completion of the microsatellite analysis and (2) Shaw et al.12 already showed the uniqueness of L. forbesii in this area based on microsatellites. All primers were marked with a fluorescent dye. The master mix for the PCR included 0.5 µL of both primers (10 µM), 1.0 µL dNTPs (10 mM, 2.5 mM each), 1.0 µL 10 × Buffer, 0.085 µL Taq-Polymerase (5 units/µL) and between 0.8 and 1.2 µL MgCL2 (25 mM), referring to the primer. This was filled up to a volume of 9 µL with water. For the PCR reaction 1 µL DNA isolate was added. The reaction conditions followed the description of Shaw et al.12. Next, the concentration of DNA was estimated using calibrated gel electrophoresis and diluted accordingly. The diluted PCR product was added to 9.8 µL HiDi™ formamide and 0.02 µL size standard GeneScanTM LIZ™ (Applied Biosystems, Foster City, California, USA). The samples were denatured at 96 °C for 4 min and fragments were size-separated in a capillary sequencer (Hitachi 3130xl Genetic Analyzer, Applied Biosystems, Foster City, California, USA). If the signal in the electropherogram was too weak, the whole procedure was repeated up to three times.

The recorded electropherograms were evaluated using "Geneious Prime" (version 2019.04, Biomatters, New Zealand). The evaluation followed the guidelines of Butler17. Thereafter, the software performed the allele binning and the length of the different alleles was determined.

For the statistical evaluation several programs were applied. The "Microsatellite Toolkit" by Park18 calculated the number of alleles per locus as well as the expected heterozygosity19 and the observed heterozygosity20. Furthermore, statistics on linkage disequilibrium and Hardy–Weinberg Principle21 were performed with "GenePop"22,23. "HP RARE"24 was used to determine allelic richness and private allelic richness, based on 15 individuals. The calculation and the pairwise comparison of the FST values were performed in "Arlequin Ver 3.5"25. Cluster analyses based on the microsatellite data were conducted with "STRUCTURE" version 2.3.4426. All settings were selected according to the standard settings for microsatellite data as recommended by Gilbert et al.27 and Porras-Hurtado et al.28. For a K-value from 1 to 10 a total of 25 iterations were calculated. Next, the best K-value was calculated with the ΔK method29, using the online software "Clumpak" (http://clumpak.tau.ac.il/)30.

To describe the phylogeography of L. forbesii, 218 individuals from all 12 sampling areas were analysed. The primers LCO1490 und HCO219831 were used to amplify a 648 bp long part of the mitochondrial cytochrome oxidase I (COI) gene. For samples with poor DNA quality, internal primers (see Supplementary Table S2 online) were designed to bisect the sequence into two shorter parts. These two internal sequences were then aligned to generate the full sequence. The PCR-Mix contained: 3 µL DNA-isolate, 3 µL dNTPs (10 mM, 2.5 mM each), 3 µL of each primer (10 µM), 3 µL 10 × Buffer, 4.8 µL MgCl2 (25 mM), 0.21 µL Taq-Polymerase (5 units/µL) and 9.99 µL water. The PCR was conducted in 38 cycles: 30 s at 94 °C, 30 s at 55 °C, 60 s at 72 °C and finally 300 s at 72 °C.

For the sequencing, the "BigDye Terminator v1.1 Kit" (Applied Biosystems, Foster City, California, USA) was used. Cycle sequencing products were analysed by using capillary separation on an ABI Genetic Analyzer 3130 xl (Applied Biosystems/Hitachi). All products were sequenced in both directions. The sequences obtained were analysed with the software "CEQ8000 Genetic Analysis System" (Version 9.0.25, Beckman Coulter GmbH, Germany) and an alignment of 249 sequences was created ("BioEdit" version 7.2.5.32) of which 28 were retrieved from GenBank32,33,34,35,36,37,38,39,40. The phylogenetic analyses (Maximum Likelihood, ML) was performed using "MEGA" version 641 with 1000 bootstraps. We determined the best fitted model of the ML method with "MEGA 6" on the basis of the Akaike Information Criterion (AIC), TN93 + G + I (Fig. 5).

"MrBayes 3.2.7"42 was used for the Bayesian Inference (BI) method. Two independent runs with four chains were performed for 2 million generations using the Markov Chain Monte Carlo (MCMC) method. Calculations of the consensus tree, including clade posterior probability (PP), were conducted based on the trees sampled after the chains converged using "Tracer 1.7"43. The first 25% were discarded as burn-in. We determined the best fitted model of the BI method with "Modeltest", implemented in "MEGA 6", on the basis of the Bayesian information criterion (BIC), T92 + G.

Finally, a haplotype network with 249 individuals was calculated (Median-Joining Method44) using "Network 5.0.1.0" and its standard settings (Fluxus Technology, Suffolk, UK). Additionally, a second haplotype network was calculated containing 313 individuals with shorter sequences (432 bp), because some individuals yielded incomplete or shorter sequences (see Supplementary Fig. S1 online). Gene flow (Nm) and FST values for the COI-gene were calculated using "DnaSP v5"45 (see Supplementary Table S3 online).

All maps were created with "Ocean Data View" version 4.6.246.

In total, 16 microsatellite loci were tested for 300 individuals of L. forbesii. The locus Lfor7 showed up to eight peaks per individual. The loci Lfor14 and Lfor15 showed additional polymorphic amplification products outside the expected size range. Consequently, all three loci were excluded from statistical analyses. In addition, 4.92% of the samples could not be analysed.

The remaining loci showed multiple alleles, except Lfor9 which was monomorphic. The number of alleles per locus ranged from 11 (Lfor2) to 49 (Lfor12). The test of linkage disequilibrium showed that all loci are inherited independently. Some deviations of expected and observed heterozygosity were found to be significant. On average across all loci for each sampled area, only non-significant deviations were determined, meaning the Hardy–Weinberg equilibrium could not be rejected.

The average private allelic richness (Table 1) across all loci ranged from 0.06 (West-Sardinia) to 1.15 (Celtic Sea) and revealed high gene flow (for more details see Appendix S4).

The pairwise FST values between the sampled populations (areas) were low but after Bonferroni-correction still some significant differences (p < 0.0009) were found (highlighted in Table 2). This is especially the case for more distant geographical sampling areas, like North Sea and Balearic Sea, or North Sea and west coast of Sardinia. However, within the Mediterranean Sea some significant differences (Ionic Sea – west and east coast of Sardinia), were identified as well.

The STRUCTURE results revealed three clusters as the most likely number of clusters for this dataset (Fig. 2; Supplementary Fig. S2). One dominant cluster occurred in the Atlantic (Fig. 2; blue), while the Mediterranean Sea showed two different clusters (Fig. 2; violet and orange) with differences between the coasts of Sardinia and the remaining areas. The Bay of Biscay, the Gulf of Cadiz and the Balearic Sea represent a transition area with admixed origin from the three clusters in different proportions according to their respective geographic position. Whereas the Bay of Biscay is more genetically related to the East Atlantic samples (North Sea, English Channel and Celtic Sea), the Bay of Cadiz appears to represent equally the Atlantic and the Mediterranean, and the Balearic Sea is clearly more related to the Adriatic, Ionian and Aegean Seas.

Structure analysis based on microsatellite data revealing three genetic clusters (blue = Atlantic cluster, orange & violet = Mediterranean cluster) for individuals from the sampled areas.

The analysis of the mitochondrial COI-sequences identified 16 different haplotypes and illustrates the existence of two major clades, one exclusive of the Azores. The other clade can be split further into two clades, one representing mostly individuals from the East Atlantic, henceforth referred to as the "East Atlantic" clade, and the other representing mostly individuals from the Mediterranean, henceforth referred to as the "Mediterranean" clade (Fig. 3). All sequences were submitted to GenBank (Accession Numbers OK135754-OK135769).

Haplotype network after Median-Joining Method for L. forbesii (haplotypes 1–10, 21–26) representing 249 COI-sequences (561 bp); black coloured GenBank data33,34,35,36,37 by various authors originate all from North East Atlantic individuals, yellow coloured including GenBank data38,39,40 .

The "East Atlantic" clade and "Mediterranean" clade were separated by at least three mutations, while the Azorean clade was separated by five and six mutations from the "Mediterranean" and the "East Atlantic" clades, respectively. Individuals caught in the eastern Mediterranean Sea (Aegean Sea, Ionic Sea, Adriatic Sea) belong exclusively to the "Mediterranean" clade, while individuals from the eastern Atlantic, with few exceptions, are members of the "East Atlantic" clade. The Gulf of Cadiz and the west and east coast of Sardinia represent a transition zone between the North Atlantic and the Mediterranean, with both clades occurring in these areas (Fig. 4).

Geographical distribution of the Azores clade (yellow), "East Atlantic" clade (blue) and the "Mediterranean" clade (orange) of L. forbesii in Europe (sample size depicted inside circles) based on COI-sequences (Ocean Data View version 4.6.2).

The shorter and longer sequences show similar haplotype networks, although some mutations are omitted in the network performed with the shorter sequences length (see Supplementary Fig. S1, Table S5 online).

The ML algorithm identify L. forbesii as a monophyletic species with 100% bootstrap support (Fig. 5). The sister group (outgroup) used consisting of Loligo vulgaris and Loligo reynaudii was identified in a more comprehensive analysis by calculating a ML tree with several loliginid species (Supplement Fig. S3). Within the species L. forbesii, two main lineages can be identified (Fig. 5). One lineage represents exclusively the individuals from the Azores and consists of the Azorean clade, whereas the other lineage consists of two more clades, the "East Atlantic" clade and the "Mediterranean" clade. The additional Bayesian phylogenetic analysis shows similar result with a comparable support of the clades (Supplement Fig. S4).

Molecular phylogenetic analysis after Maximum Likelihood (ML) method (1000 Bootstraps). Labelling next to the species names indicate the GenBank Accession Numbers for L. reynaudii and L. vulgaris.

Based on clades defined by COI-sequence data, populations predominantly belonging to the same clade show low pairwise FST values based on microsatellites as well as relatively high gene flow rates (Table S5). For example, the North Sea individuals show high gene flow rates with other east Atlantic individuals like the fished individuals from the English Channel or the Celtic Sea (Nm = 16.49 and Nm = 12.97, respectively). Conversely, gene flow between the North Sea individuals and those from the eastern Mediterranean Sea, e.g. Ionic Sea and Aegean Sea, is far more restricted (Nm = 0.03 and Nm = 0.03, respectively; Table S5).

Our results based on L. forbesii microsatellites analysis illustrate the existence of three genetic units along the Atlantic and Mediterranean shelf area (North East Atlantic, western and eastern Mediterranean) with high gene flow between populations within and, partially, between the genetic clusters. Due to our substantially expanded sampling area and increased number of microsatellite loci in comparison to previous studies, we were able to identify a higher number of alleles giving us a more robust dataset for gene flow inferences12,15,16,47,48.

Locus Lfor7 showed up to eight peaks per individual, which was never described before12,47. This inconsistency between different studies may be due to different methods of allele scoring. Furthermore, loci Lfor14 and Lfor15 showed additional polymorphic amplification products. As unspecific priming and chain-termination can be excluded, we consider that these amplifications may eventually represent repeated regions in the genome. Not every locus delivers suitable results for every sample. This is not unusual as microsatellites are sensitive and the quality of the tissue sample plays an important role. Our dataset has 4.92% missing data, which is slightly higher than the 3.75% in Shaw and Boyle47. Not every locus was affected equally by missing data. The differences in the present study are most likely due to low DNA quality. This was especially the case for some samples from the Bay of Biscay and Adriatic Sea for which the ethanol may have not been changed adequately.

The observed low values of private allelic richness, as well as the low pairwise FST values are indicative of a considerable gene flow across the entire range of this species. Only sampling areas with a substantial geographic distance show low significant differences. Our STRUCTURE results show clear differences between Atlantic and Mediterranean L. forbesii individuals. The genetic composition of individuals from the Gulf of Cadiz to the western Mediterranean assigning to two or three genetic clusters is an evidence for recent gene flow, or otherwise a short period of genetic isolation. It is well known that the relatively shallow area of the Strait of Gibraltar has the potential to be a barrier for some species49,50. But it is not surprising that the area represents a transition zone for a shelf-living species such as L. forbesii. In addition, admixture within the Atlantic and the Mediterranean sampling sites seems to occur, represented by high gene flow rates (Supplementary Table S3), which is in line with the known migration behaviour of L. forbesii2. Other members of the taxon of Loliginidae also show gene flow over long geographic distances. Populations from Doryteuthis opalescens (formerly known as Loligo opalescens) and Doryteuthis gahi (formerly known as Loligo gahi) show poor genetic differentiation as well51,52. Our results from northern Atlantic sample sites correspond in general to the findings of Shaw et al.12 and Brierly et al.11. Hence, although we were not able to incorporate the Azorean samples in our microsatellite analysis, we anticipate that the Azores would represent a fourth microsatellite cluster as Shaw et al.12 uncovered significant genetic differences between L. forbesii from the archipelago and the North East Atlantic Shelf area.

We found a total of 16 COI-haplotypes for L. forbesii, whereof three have been reported before32,33,34,35,36,37,38,39,40 (H7 and H10, "East Atlantic" clade; H25, Azorean clade, see Fig. 3). We found three different mitochondrial clades in European waters: the exclusive Azorean clade, a clade dominated by East Atlantic individuals and a clade dominated by Mediterranean individuals. The two latter clades are not geographically isolated, because some individuals from the Atlantic and some individuals from the Mediterranean share some haplotypes from each clade. In detail, the "East Atlantic" clade is dominant in the Atlantic, while in the more eastern Mediterranean (Adriatic Sea, Ionian Sea, Aegean Sea) only members of the "Mediterranean" clade were found. Both clades co-occur in the Gulf of Cadiz as well as at the west and east coast of Sardinia, indicating an exchange of individuals between the Atlantic and Mediterranean Sea via the Strait of Gibraltar. Around the Balearic islands only members of the Mediterranean clade were found which may reflect the low number of samples. We found three North Sea individuals with Mediterranean haplotypes that were juvenile of 103–130 days old53 and were caught at a known spawning area during spawning time54. The elektropherogramms revealed that a sibling or half-sibling relationship can be excluded for these three individuals, additionally indicating multiple events of gene flow from the Mediterranean into the Atlantic. It is known that Loligo reynaudii, a closely related species, performs migrations over large geographical distances during spawning time and does not exhibit homing behavior55. When ready to spawn, individuals of this species show an average migration speed of 3 km per day, and exceptional migration events were also reported, as a movement of 207 km in 18 days55. While this information is not available for L. forbesii, it may be hypothesised that the spawning behaviour and the ability to move over long distances of L. reynaudii may be also present in L. forbesii.

The three mitochondrial clades could also be differentiated from the phylograms (Fig. 5, Supplement Fig. S3). The bootstrap values are rather low most likely reflecting the less than 1% overall differences for the COI-sequence data, but eventually also because the data was not partitioned by codon position. This low differentiation let us assume that the separation of the clades occurred in the Pleistocene56 (Supplement Fig. S5). Our phylogeographic results are concordant with COI-data from Sepia officinalis, a species with a similar distribution to L. forbesii. Perez-Losada et al.57 inferred that S. officinalis expanded from Northwestern Europe coast into the Mediterranean, and observed that Aegean and Ionian populations of S. officinalis are genetically clearly separated from the remaining Mediterranean. We also found this differentiation between Eastern and Western Mediterranean but for L. forbesii mitochondrial haplotypes are generally more admixed between the sampling sites. This result may be due to gene flow following secondary contact between these clusters. Differences between local patterns of L. forbesii and S. officinalis may correspond to the different dispersal ability of the two species, with L. forbesii being a way more mobile species. Both the geological history of the Mediterranean and physical barriers between eastern and western Mediterranean Sea seem to have had a huge impact on the phylogeographical patterns of both cephalopod species.

The microsatellite and mitochondrial COI analyses show concurring results in population structure of L. forbesii. Comparing the same geographical range for both markers, we found two mitochondrial clades and three microsatellite clusters, which expectably express the different mutation rates of these markers. As microsatellites are known to better describe recent patterns of connectivity between populations13, our results reveal a recent differentiation into three clades within the European Shelf population.

Shaw et al.12 described the influence of geographic barriers, like deep sea basins, on the migration patterns of L. forbesii. The Italian peninsula as well as the Hellenic trench are potential migration barriers in the Mediterranean Sea. The offshore location of the Azores is very likely the reason for the genetic isolation between L. forbesii from the archipelago and the remaining North East Atlantic shelf area12, as the squids prefer the upper shelf slope2,4. Active migrations over large distances also seem to be important58 regarding the exchange between the Atlantic and Mediterranean over the Strait of Gibraltar. The two-layer water exchange in the Strait of Gibraltar described by Izquierdo et al.59 may influence the distribution of L. forbesii. Less dense Atlantic water forms a current in the upper water layer which flows into the Mediterranean, while the denser Mediterranean water forms the Mediterranean outflow, an underwater current which brings Mediterranean water masses into the Atlantic. This water exchange may affect differently adult L. forbesii, which spend the daytime close to the ground and ascend to surface waters at night for feeding2. This way, adults are directly influenced by the water flow in both directions, which supports the exchange of individuals between Atlantic and Mediterranean triggering the genetic admixture that can be seen in the present results. It is also expectable that the transport of paralarvae is influenced by currents.

With our new findings of the COI-gene analysis, we are able to describe uncovered phylogeographic patterns for L. forbesii. We found a subtle genetic differentiation between the Atlantic and the Mediterranean, but statoliths of Loligo, indistinguishable from those of current L. forbesii and L. vulgaris, were collected in early Miocene deposits at southern French shores60. This means that an ancestral of the current L. forbesii population existed in the Mediterranean sea, but we anticipate that the final establishment of the current Mediterranean population took place only after the Messinian Salinity Event 5.5 million years ago61,62, as the survival of the species during this crisis is very unlikely. We thus assume a repetitive colonialization of the Mediterranean by L. forbesii from the Atlantic across the Strait of Gibraltar. A geographical separation of the species took apparently place during the last ice age in the Late Pleistocene63, and the East Atlantic and Mediterranean lineage could evolve in allopatry. Today, the two previously evolved lineages have come into contact, supported by the combination of gene flow rates and the distribution of clade-specific haplotypes.

Under the assumption of the above mentioned scenario, the Sardinian cluster may have established following a second wave of immigration into the Mediterranean Sea. This cluster is represented in mitochondrial COI diversity by the "East Atlantic" haplotypes 5 and 6, which were only found around Sardinia. However, an alternative explanation for the genetic patterns observed is that a subpopulation of an ancestral population of L. forbesii that inhabited both the Atlantic and Mediterranean, got isolated in the Mediterranean due to sea level change eventually during the Late Pleistocene glaciations, becoming later the common ancestor of the three clades observed today, the "East Atlantic" clade, the "Mediterranean clade " and the endemic clade in the Azores.

Based on our results, we recommend that future fishery management should consider at least three different genetic groups of L. forbesii. Even knowing that the structure analysis illustrates admixture in some regions, a more accurate classification into smaller management units seems not possible using the current genetic markers. Therefore, understanding the complex population structure of L. forbesii and differentiate stocks for a more precise management needs further research and more alternative methods. Some studies indicate that statolith shape analysis and statolith elemental analysis or the combination of both might be potential methods to identify smaller management units64.

Doubleday, Z. A. et al. Global proliferation of cephalopods. Curr. Biol. 26, R406–R407 (2016).

Article CAS PubMed Google Scholar

Jereb, P. et al. Cephalopod biology and fisheries in Europe: II. Species Accounts. ICES Cooperative Research Report No vol. 325 (2015).

ICES. ICES WGCEPH REPORT 2015 Interim Report of the Working Group on Cephalopod Fisheries and Life History (WGCEPH). 8–11 (2019).

Quetglas, A. et al. Long-term spatiotemporal dynamics of cephalopod assemblages in the Mediterranean sea. Sci. Mar. 83, 33–42 (2019).

Article Google Scholar

Martins, H. R. Biological studies of the exploited stock of Loligo forbesi (Mollusca: Cephalopoda) in the Azores. J. Mar. Biol. Assoc. United Kingdom 62, 799–808 (1982).

Article Google Scholar

Guerra, A. & Rocha, F. The life history of Loligo vulgaris and Loligo forbesi (Cephalopoda: Loliginidae) in Galician waters (NW Spain). Fish. Res. 21, 43–69 (1994).

Article Google Scholar

Pierce, G. J. & Boyle, P. R. Empirical modelling of interannual trends in abundance of squid (Loligo forbesi) in Scottish waters. Fish. Res. 59, 305–326 (2003).

Article Google Scholar

Lishchenko, F. et al. A review of recent studies on the life history and ecology of European cephalopods with emphasis on species with the greatest commercial fishery and culture potential. Fish. Res. 236, 105847 (2021).

Article Google Scholar

Laptikhovsky, V. et al. Identification of benthic egg masses and spawning grounds in commercial squid in the English Channel and Celtic Sea: Loligo vulgaris vs L. forbesii. Fish. Res. 241, 106004 (2021).

Article Google Scholar

Souza, H. V. et al. Analysis of the mitochondrial COI gene and its informative potential for evolutionary inferences in the families Coreidae and Pentatomidae (Heteroptera). Genet. Mol. Res. 15, 1–14 (2016).

CAS Google Scholar

Brierley, A. S. et al. Genetic variation in the neritic squid Loligo forbesi (Myopsida: Loliginidae) in the northeast Atlantic Ocean. Mar. Biol. 122, 79–86 (1995).

Article Google Scholar

Shaw, P. W. et al. Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers. Mol. Ecol. 8, 407–417 (1999).

Article CAS Google Scholar

Ellegren, H. Microsatellites: Simple sequences with complex evolution. Nat. Rev. Genet. 5, 435–445 (2004).

Article CAS PubMed Google Scholar

Begg, G. A. & Waldman, J. R. An holistic approach to fish stock identification. Fish. Res. 43, 35–44 (1999).

Article Google Scholar

Shaw, P. W. Polymorphic microsatellite markers in a cephalopod: The veined squid Loligo forbesi. Mol. Ecol. 6, 297–298 (1997).

Article CAS Google Scholar

Emery, A. M. et al. New microsatellite markers for assessment of paternity in the squid Loligo forbesi (Mollusca: Cephalopoda). Mol. Ecol. 9, 110–112 (2000).

Article CAS PubMed Google Scholar

Butler, J. M. Advanced Topics in Forensic DNA Typing: Interpretation (Elsevier Academic Press, 2015).

Google Scholar

Park, S. D. E. Trypanotolerance in West African Cattle and the Population Genetics Effects of Selection. Trinity Coll. (2001).

Nei, M. Molecular Evolutionary Genetics (Columbia University Press, 1987).

Google Scholar

Hedrick, P. W. Genetics of Populations (Science Books International, 1983).

Google Scholar

Weir, B. S. & Cockerham, C. C. Estimating F statistics for Population Structure. Evolution 38, 1358–1370 (1984).

CAS PubMed Google Scholar

Raymond, M. & Rousset, F. GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. J. Hered. 86, 248–249 (1995).

Article Google Scholar

Rousset, F. GENEPOP’007: A complete re-implementation of the GENEPOP software for Windows and Linux. Mol. Ecol. Resour. 8, 103–106 (2008).

Article PubMed Google Scholar

Kalinowski, S. T. HP-RARE 1.0: A computer program for performing rarefaction on measures of allelic richness. Mol. Ecol. Notes 5, 187–189 (2005).

Article CAS Google Scholar

Excoffier, L. et al. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinforma. 1, 117693430500100 (2005).

Article Google Scholar

Pritchard, J. K. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).

Article CAS PubMed PubMed Central Google Scholar

Gilbert, K. J. et al. Recommendations for utilizing and reporting population genetic analyses: The reproducibility of genetic clustering using the program structure. Mol. Ecol. 21, 4925–4930 (2012).

Article PubMed Google Scholar

Porras-Hurtado, L. et al. An overview of STRUCTURE: Applications, parameter settings, and supporting software. Front. Genet. 4, 1–13 (2013).

Article CAS Google Scholar

Evanno, G. et al. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).

Article CAS PubMed Google Scholar

Kopelman, N. M. et al. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 15, 1179–1191 (2015).

Article CAS PubMed PubMed Central Google Scholar

Folmer, O. et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 3, 294–299 (1994).

CAS PubMed Google Scholar

Hall, T. A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95–98 (1999).

CAS Google Scholar

Anderson, F. E. Phylogeny and historical biogeography of the loliginid squids (Mollusca: Cephalopoda) based on mitochondrial DNA sequence data. Mol. Phylogenet. Evol. 15, 191–214 (2000).

Article CAS PubMed Google Scholar

Gebhardt, K. & Knebelsberger, T. Identification of cephalopod species from the North and Baltic Seas using morphology, COI and 18S rDNA sequences. Helgol. Mar. Res. 69, 259–271 (2015).

Article ADS Google Scholar

Lobo, J. et al. Enhanced primers for amplification of DNA barcodes from a broad range of marine metazoans. BMC Ecol. 13, 1–8 (2013).

Article CAS Google Scholar

de Luna Sales, J. B. et al. New molecular phylogeny of the squids of the family Loliginidae with emphasis on the genus Doryteuthis Naef ,1912: Mitochondrial and nuclear sequences indicate the presence of cryptic species in the southern Atlantic Ocean. Mol. Phylogenet. Evol. 68, 293–299 (2013).

Article Google Scholar

Tatulli, G. et al. A rapid colorimetric assay for on-site authentication of cephalopod species. Biosensors 10, 3–10 (2020).

Article CAS Google Scholar

Velasco, A. et al. A new rapid method for the authentication of common octopus (Octopus vulgaris) in seafood products using recombinase polymerase amplification (rpa) and lateral flow assay (lfa). Foods 10, 1825 (2021).

Article CAS PubMed PubMed Central Google Scholar

Luz, A. & Keskin, E. Building Reference Library for Marine Fish Species of Azores Archipelago and Bio-monitoring via DNA Metabarcoding. https://www.ncbi.nlm.nih.gov/nuccore/MT491734 (2020).

BoldSystems. https://boldsystems.org/index.php/Public_RecordView?processid=AZB030-20 (2018). (Accessed 2 May 2022).

Tamura, K. et al. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013).

Article CAS PubMed PubMed Central Google Scholar

Ronquist, F. & Huelsenbeck, J. P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 1572–1574 (2003).

Article CAS PubMed Google Scholar

Rambaut, A. et al. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 67, 901–904 (2018).

Article CAS PubMed PubMed Central Google Scholar

Bandelt, H.-J. et al. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (2009).

Article Google Scholar

Librado, P. & Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).

Article CAS PubMed Google Scholar

Schlitzer, R. Ocean Data View. (2013).

Shaw, P. W. & Boyle, P. R. Multiple paternity within the brood of single females of Loligo forbesi (Cephalopoda: Loliginidae), demonstrated with microsatellite DNA markers. Mar. Ecol. Prog. Ser. 160, 279–282 (1997).

Article ADS Google Scholar

Emery, A. M. et al. Assignment of paternity groups without access to parental genotypes: Multiple mating and developmental plasticity in squid. Mol. Ecol. 10, 1265–1278 (2001).

Article CAS PubMed Google Scholar

Catarino, D. et al. The role of the Strait of Gibraltar in shaping the genetic structure of the Mediterranean Grenadier, Coryphaenoides mediterraneus, between the Atlantic and Mediterranean Sea. PLoS ONE 12, 1–24 (2017).

Google Scholar

Gonzalez, E. G. & Zardoya, R. Relative role of life-history traits and historical factors in shaping genetic population structure of sardines (Sardina pilchardus). BMC Evol. Biol. 7, 1–12 (2007).

Article CAS Google Scholar

Reichow, D. & Smith, M. J. Microsatellites reveal high levels of gene flow among populations of the California squid Loligo opalescens. Mol. Ecol. 10, 1101–1109 (2001).

Article CAS PubMed Google Scholar

Shaw, P. W. et al. DNA markers indicate that distinct spawning cohorts and aggregations of Patagonian squid, Loligo gahi, do not represent genetically discrete subpopulations. Mar. Biol. 144, 961–970 (2004).

Article CAS Google Scholar

Göpel, A. Populationsgenetik und Phylogeographie des Nordischen Kalmars Loligo forbesii Steenstrup, 1856 in Europäischen Gewässern. Masterthesis, Univ. Rostock in German, 76pp (2020).

Oesterwind, D. et al. Biology and meso-scale distribution patterns of North Sea cephalopods. Fish. Res. 106, 141–150 (2010).

Article Google Scholar

Sauer, W. H. H. et al. Tag recapture studies of the chokka squid Loligo vulgaris reynaudii d’Orbigny, 1845 on inshore spawning grounds on the south-east coast of South Africa. Fish. Res. 45, 283–289 (2000).

Article ADS Google Scholar

Knowlton, N. & Weigt, L. A. New dates and new rates for divergence across the Isthmus of Panama. Proc. R. Soc. B Biol. Sci. 265, 2257–2263 (1998).

Article Google Scholar

Pérez-Losada, M. et al. Testing hypotheses of population structuring in the Northeast Atlantic Ocean and Mediterranean Sea using the common cuttlefish Sepia officinalis. Mol. Ecol. 16, 2667–2679 (2007).

Article PubMed Google Scholar

O’Dor, R. K. Can understanding squid life-history strategies and recruitment improve management?. South African J. Mar. Sci. 7615, 193–206 (1998).

Article Google Scholar

Izquierdo, A. et al. Modelling in the Strait of Gibraltar: From operational oceanography to scale interactions. Fundam. i Prikl. Gidrofiz. 9, 15–24 (2016).

Google Scholar

Clarke, M. & Hart, M. Treatise Online no. 102: Part M, Chapter 11: Statoliths and coleoid evolution. Treatise Online (2018).

Hsü, K. J. et al. Late Miocene desiccation of the mediterranean. Nature 242, 240–244 (1973).

Article ADS Google Scholar

Garcia-Castellanos, D. et al. Catastrophic flood of the Mediterranean after the Messinian salinity crisis. Nature 462, 778–781 (2009).

Article ADS CAS PubMed Google Scholar

Thunell, R. C. et al. Atlantic-mediterranean water exchange during the late neocene. Paleoceanography 2(6), 661 (1987).

Article ADS Google Scholar

Green, C. P. et al. Combining statolith element composition and fourier shape data allows discrimination of spatial and temporal stock structure of arrow squid (Nototodarus gouldi). Can. J. Fish. Aquat. Sci. 72, 1609–1618 (2015).

Article Google Scholar

Download references

This project has been partly funded by the EU through the European Maritime and Fisheries Fund (EMFF) within the Spanish National Program of collection, management and use of data in the fisheries sector and support for scientific advice regarding the Common Fisheries Policy and genetic sampling was supported by the Cephs and Chefs INTERREG project. Sampling in the Azores was funded by CESAM (UIDP/50017/2020+UIDB/50017/2020+LA/P/0094/2020) that is financed by FCT/MCTES through national funds. We thank Christopher Zimmermann for the funding support of the study, Matthias Kloppmann and his crew for his considerate cruise lead during sampling in the North Sea and Nicholas Badouvas, Nikolaos Fotiadis, for their sample preparation at HCMR. We also thank Lukas Krebes and Sören Möller from the University of Rostock for sharing their experiences with the microsatellite analysis. Last but not least, we would like to thank the students Vivian Fischbach and Chantal Petong for their valuable lab work, two anonymous reviewers for their very constructive comments and especially the editor (Raquel Godinho) for her patient management of the review process and her huge support to get this manuscript published.

Open Access funding enabled and organized by Projekt DEAL.

Thünen Institute of Baltic Sea Fisheries, Alter Hafen Süd 2, 18069, Rostock, Germany

Anika Göpel & Daniel Oesterwind

Institute of Biological Sciences, University of Rostock, Albert-Einstein-Str. 3, 18059, Rostock, Germany

Anika Göpel & Ralf Bastrop

Cefas Laboratory, Pakefield Rd, Lowestoft, NR33 0HT, UK

Christopher Barrett & Vladimir Laptikhovsky

Hellenic Centre for Marine Research, Institute of Marine Biological Resources and Inland Waters, 576 SideRD Vouliagmenis Ave, 16452, Athens, Greece

Evgenia Lefkaditou

Technical University of Denmark, National Institute of Aquatic Resources, Nordsøen Forskerpark, Willemoesvej 2, 9850, Hirtshals, Denmark

Kai Wieland

University of Caen Normandy, CS 14032, 14032, Caen Cedex 05, France

Angela Larivain & Jean-Paul Robin

Centre Oceanográfic de les Balears s/n, Instituto Español de Oceanografía (IEO), 07015, Palma, Spain

Maria Valls

Department of Life and Environmental Sciences, University of Cagliari, Cagliari, Italy

Rita Cannas & Maria Cristina Follesa

COISPA Tecnologia & Ricerca, Via dei Trulli, 18-20, Bari, Italy

Pierluigi Carbonara & Marilena Donnaloia

Centro Oceanográfico de Cádiz, Instituto Español de Oceanografía, Puerto Pesquero, Muelle de Levante S/N, 11006, Cádiz, Spain

Luis Silva Caparro & Ignacio Sobrino

Centro Oceanográfico de Vigo, Instituto Español de Oceanografía (IEO), Subida a Radio Faro, 50, 36390, Vigo, Spain

Maria Begoña Santos & Julio Valeiras

CESAM - Centre for Environmental and Marine Studies, Department of Biology, University of Aveiro, Campus de Santiago, 3810-193, Aveiro, Portugal

Hugo C. Vieira

You can also search for this author in PubMed Google Scholar

A.G.: Performed the laboratory analysis and statistics and wrote the draft manuscript and led the revision; D.O.: conceived the study project, organized samples and contributed to the drafting of the manuscript and led the revisions; R.B.: supervised lab work and statistics and contributed to the drafting of the manuscript and to the revisions; C.B., R.C., L.S.C., P.C., M.D., M.C.F., A.L., V.L., E.L., J.-P.R., M.B.S., I.S., J.V., M.V., K.W.: prepared samples and reviewed the draft manuscript and the revisions. H.C.V. prepared Azorean samples and contributed to the revisions of the manuscript. All authors read and approved the manuscript.

Correspondence to Daniel Oesterwind.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Göpel, A., Oesterwind, D., Barrett, C. et al. Phylogeography of the veined squid, Loligo forbesii, in European waters. Sci Rep 12, 7817 (2022). https://doi.org/10.1038/s41598-022-11530-z

Download citation

Received: 10 January 2021

Accepted: 21 April 2022

Published: 12 May 2022

DOI: https://doi.org/10.1038/s41598-022-11530-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Marine Biology (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

News

Phylogeography of the veined squid, Loligo forbesii, in European waters