Hostname: page-component-7b9c58cd5d-v2ckm Total loading time: 0 Render date: 2025-03-15T22:47:55.089Z Has data issue: false hasContentIssue false

Integrating Anisakis spp. parasites data and host genetic structure in the frame of a holistic approach for stock identification of selected Mediterranean Sea fish species

Published online by Cambridge University Press:  22 August 2014

S. MATTIUCCI*
Affiliation:
Department of Public Health and Infectious Diseases, Section of Parasitology, ‘Sapienza University of Rome’, Rome, Italy
R. CIMMARUTA
Affiliation:
Department of Ecological and Biological Sciences, Tuscia University, Viterbo, Italy
P. CIPRIANI
Affiliation:
Department of Public Health and Infectious Diseases, Section of Parasitology, ‘Sapienza University of Rome’, Rome, Italy Department of Ecological and Biological Sciences, Tuscia University, Viterbo, Italy
P. ABAUNZA
Affiliation:
Instituto Español de Oceaonografía, Oceanographic Centre of Santander, Santander, Spain
B. BELLISARIO
Affiliation:
Department of Ecological and Biological Sciences, Tuscia University, Viterbo, Italy
G. NASCETTI
Affiliation:
Department of Ecological and Biological Sciences, Tuscia University, Viterbo, Italy
*
* Corresponding author: Department of Public Health and Infectious Diseases, Section of Parasitology, Sapienza University of Rome, P.le Aldo Moro, 5 00185 Rome, Italy. E-mail: simonetta.mattiucci@uniroma1.it
Rights & Permissions [Opens in a new window]

Summary

The unique environment of the Mediterranean Sea makes fish stock assessment a major challenge. Stock identification of Mediterranean fisheries has been based mostly from data on biology, morphometrics, artificial tags, otolith shape and fish genetics, with less effort on the use of parasites as biomarkers. Here we use some case studies comparing Mediterranean vs Atlantic fish stocks in a multidisciplinary framework. The generalized Procrustes Rotation (PR) was used to assess the association between host genetics and larval Anisakis spp. datasets on demersal (hake) and pelagic (horse mackerel, swordfish) species. When discordant results emerged, they were due to the different features of the data. While fish population genetics can detect changes over an evolutionary timescale, providing indications on the cohesive action of gene flow, parasites are more suitable biomarkers when considering fish stocks over smaller temporal and spatial scales, hence giving information of fish movements over their lifespan. Future studies on the phylogeographic analysis of parasites suitable as biomarkers, and that of their fish host, performed on the same genes, will represent a further tool to be included in multidisciplinary studies on fish stock structure.

Type
Fisheries
Copyright
Copyright © Cambridge University Press 2014 

INTRODUCTION

A stock can be defined as a group of individuals of a particular fish species whose genetic and life-history characteristics are more similar to each other, rather than to the other subpopulations (or stocks) (Waldman, Reference Waldman, Cadrin, Friedland and Waldman2005). The stock concept has been introduced to better understand the state of over-exploited fish subpopulations, and to highlight which areas within the species range are at higher risk of being over-fished. The correct identification of biological stocks, including the determination of their boundaries between adjacent conspecific subpopulations characterized by differences in abundance, growth, reproduction, spawning and genetic diversity, is a useful approach in the management of economically important fish species, since it is the unit of both management practices and conservation measures (Waples, Reference Waples1998).

Marine fish have been thought to have little genetic subdivision among geographically isolated populations (Ward et al. Reference Ward, Woodwark and Skibinski1994). Indeed the gene flow, through the passive transport of pelagic eggs and larvae, or by active adult migration over large distances, contributes to maintaining the genetic homogeneity of large and stable fish populations of the same species over large geographical scales. Recent data have partially modified this impression, highlighting geographic structure in many different fish species (Cowen and Sponaugle, Reference Cowen and Sponaugle2009). Although the use of genetic/molecular data has become a common basic method to identify fish population structure, allowing the evaluation of genetic differentiation and gene flow between subpopulations across the geographical distribution of a fish species, the assessment of fish stock is based nowadays on a holistic approach (Waldman, Reference Waldman1999). This includes the use of fish morphometrics, artificial tags, fish genetics and other biological tags to clarify stock structure. Among the biological tags, the parasite-based methodology (including parasite community structure of the fish species along its geographical distribution, parasitic infection levels and genetic/molecular characterization of parasite species) represents an important approach in defining a fish stock. Indeed, parasites can provide ecological information on the origin, migration, nursery ground and life history of the fish species (Thomas et al. Reference Thomas, Verneau, de Meeûs and Renaud1996).

The basic principle of the use of parasites as biological tags is that a fish becomes infected by a parasite species only when it is in the endemic area of that parasite. The endemic area of the parasite is the geographic region where the abiotic (temperature, salinity) and biotic (presence of suitable intermediate and definitive hosts) factors are suitable for the transmission and completion of its life cycle. Thus, we can assume that when a fish population is found infected by a parasite species, it means that the fish has spent part of its life history in the endemic area of the parasite, where fish behaviour and feeding habits could result in different infection levels by that parasite species. As a consequence, when the parasite fauna of two populations of the same fish species sampled from two different geographic areas is different, it means that the life history of those fish samples were different (MacKenzie et al. Reference MacKenzie, Campbell, Mattiucci, Ramos, Pinto and Abaunza2008).

Basically, a parasite can be used as a suitable biological tag for fish stock identification when its geographical distribution and life cycle are known, and when the parasite's residence time in the host is long enough compared with the lifespan of the fish host. In this sense, the parasite as biological marker reflects the geographic origin of the fish population on a spatial scale. This is the main reason why concordant and discordant results in the existence of different stocks of a fish species between the genetic/molecular datasets on the fish and the results from the use of some parasites as biological tags, may occur. Indeed, the genetic/molecular markers define the stock on the evolutionary temporal scale, while the parasite taxa characterize the stock on a spatial/geographical scale.

This paper is a review of the research carried out on the use of parasites as biological markers to assess the stock structure of some pelagic and demersal fish species in the Mediterranean Sea. Specifically, it focuses on the research carried out, in recent years, on the simultaneous use of population genetics analysis of the fish species and biological markers (parasites of the genus Anisakis) to assess sub-population structuring of some pelagic and demersal fish species in the Mediterranean Sea, in comparison with populations of the same fish species in the Atlantic Ocean.

THE MEDITERRANEAN SEA AND FISHERIES

The Mediterranean Sea is a semi-closed sea with a surface of about 3 million km2 (including the Black Sea) representing 0·8% of the total world marine surface area. Its range of latitude is from 30°N to 46°N, in a temperate zone of the northern hemisphere. It consists of two connected mid-latitude basins, the Western (WM) and the Eastern Mediterranean (EM) separated by the shallow (400 m) sill of the Sicily Channel, and is characterized by a limited water exchange with the Atlantic Ocean. The thermohaline circulation of the Mediterranean Sea is generally described as an open basin-wide cell, resulting in the transformation of surface water of Atlantic origin, on its eastward progress into Levantine Intermediate Water (LIW), and two closed secondary cells, one in the Western Mediterranean (WM) and the other in the EM, which involve the transformation of the surface and intermediate water into Eastern and Western deep waters (Gačić et al. Reference Gačić, Schroeder, Civitarese, Vetrano and Eusebi Borzelli2012). Water mass is stratified during the summer, however the temperature is around 13 °C in the deep water (below the 400 m) all year round. The absence of precipitation during the summer period is the main characteristic of the Mediterranean Sea. The loss of water through evaporation is partially balanced by the contribution of Atlantic water through the Strait of Gibraltar. In the Mediterranean Sea the evaporation exceeds precipitation and run-off: the salinity increase, due to evaporation over the Mediterranean surface, is compensated by an outflow of salty water into the Atlantic Ocean. Cooling in winter causes convection to intermediate depths forming the LIW. This salty intermediate water returns to the west underneath the AW and becomes the main component of the Mediterranean Outflow to the Atlantic. The LIW is a crucial component of the Mediterranean thermohaline, because it is the warmest and saltiest Mediterranean water mass formed in large amounts, and because it mainly flows westwards along the northern continental slopes of both basins, influencing the formation of North Atlantic Deep water (Gačić et al. Reference Gačić, Schroeder, Civitarese, Vetrano and Eusebi Borzelli2012).

The Mediterranean and Atlantic regions are long-recognized marine zoogeographic provinces. Forty-five per cent of the Mediterranean fish species that co-occur in the Atlantic are considered to have invaded the Mediterranean since the opening of the Gibraltar Strait in the early Pliocene (Quignard, Reference Quignard1978). Dramatic palaeogeographic and palaeoclimatic events have played a crucial role in the evolutionary history of many marine species in the Mediterranean region (Valsecchi et al. Reference Valsecchi, Pasolini, Bertozzi, Garoia, Ungaro, Vacchi, Sabelli and Tinti2005). The closure between the Atlantic Ocean and the Mediterranean Sea, which occurred between about 5·96 and 5·33 mya, determined the Messinian salinity crisis. The re-opening of the Strait of Gibraltar (about 5·33 mya) was characterized by the colonization events.

Mediterranean fisheries actually represent a small proportion of the total European production, only 20% of the total catch; the mean annual landings (mainly sold as fresh) have been around 1·5 million of tons in the last decade. The Mediterranean borders 21 countries (excluding the Black Sea); among them, the European countries (Spain, France, Greece and Italy) represent one third of the Mediterranean Sea coastline. From a socio-economic point of view, fishing activity is highly diverse and varies by geographic area and country. This is mainly related to different environmental conditions and is due to various socio-economic conditions in the countries whose fisheries are based in the Mediterranean Sea. The Mediterranean fisheries are characterized by fragmented fleets, usually composed of relatively small vessels, multi-species catches and a large number of landing sites. In general, the lack of connection between fish stock assessment and management is the main problem for the Mediterranean fisheries. This is because available historical data on landings are scarce and, given the complexity and diversity of the Mediterranean fisheries, insufficient data are available for the assessment of most fish species (Lleonart and Maynou, Reference Lleonart and Maynou2003). Ninety-five per cent of the fish stocks are overfished, which means that catches are higher than they should be in order to make the fish stocks grow to sustainable sizes. The management of the Mediterranean fisheries is based mainly on effort control. No Total Allowable Catches (TAC) are implemented (except for bluefin tuna and swordfish), nor any other type of adaptive management. Other technical measures, such as closed seasons and areas, minimum landing sizes and minimum mesh sizes, are also implemented, but they are fewer than in the Atlantic fisheries (Lleonart and Maynou, Reference Lleonart and Maynou2003).

Large pelagic, small pelagic and demersal fisheries are the three main different groups of fish resources in the Mediterranean Sea. The bluefin tuna (Thunnus thynnus) and the swordfish (Xiphias gladius) are the most important large pelagic species in the Mediterranean Sea; they represent around 4% of the total reported landings. The International Commission of Atlantic Tunas (ICCAT) recognizes the existence of a single stock of bluefin tuna for the North-Eastern Atlantic Ocean. Swordfish is the second large pelagic fish species, and ICCAT recognizes the existence so far of a single stock unit of swordfish in the Mediterranean Sea, distinct from those of the Atlantic Ocean (see paragraph devoted to X. gladius). Small (anchovies, sardine) and medium-size (mackerels, horse mackerel, bogue, etc.) pelagic fish are the main contributors (more than 50%) to the total landings in the Mediterranean. Demersal species in the Mediterranean Sea include more than 100 species; they have commercial value, despite the large number of species, according to the FAO, and they represent only 3% of the total catch (Lleonart and Maynou, Reference Lleonart and Maynou2003).

CURRENT METHODOLOGIES USED FOR THE IDENTIFICATION OF MEDITERRANEAN FISH STOCK UNITS

Definition of fish unit stocks has been largely absent or in some cases is still problematic in the Mediterranean fisheries (Caddy, Reference Caddy1995; Cacaud, Reference Cacaud2005; Palomera et al. Reference Palomera, Olivar, Salat, Sabates, Coll, Garcıa and Morales-Nin2007). The Scientific Advisory Committee of the GFCM (General Fisheries Commission for the Mediterranean) has established management areas based on political and statistical considerations, rather than on biological ones (FAO, 1999). The Mediterranean basin has been split into 27 Geographical Subareas (GSA) for statistical purposes and data on captures and production are reported accordingly (GFCM, 2012). The availability of appropriate data is considered one of the main shortcomings in the assessment of the Mediterranean fisheries. There is a large amount of reliable data but they are not often adequate to fit standard stock assessment methods. However, scientists have tried to assess the fisheries with several tools, using available data. Analytical methods based on population dynamics (long-term series of size/age relationships, growth, etc.) have been applied to different stocks of demersal, small pelagic and large pelagic fishes. However, approaches using a single method are now considered unsatisfactory in fish stock identification (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ).

On the other hand, the marine environment and the historical data from the Mediterranean Sea basin make stock assessment of pelagic and demersal fish species, with a wide geographic distribution, a challenging argument to be investigated using different approaches. Generally, marine fish stock identification has been challenging, due to the few physical barriers and strong ocean currents existing in the marine realm, which facilitate extensive gene flow between fish populations. Even in the absence of physical barriers and within small geographic distances, however, marine fish species can exhibit cryptic population structuring based on differences in their recruitment, spawning behaviour, and mechanisms of retention of larval stages and general life-history strategies and traits (Knutsen et al. Reference Knutsen, Jorde, André and Stenseth2003). However, the Mediterranean and Atlantic populations of marine organisms were isolated due to the sea level changes during the glacial maxima and then reconnected several times in the past. The genetic divergence of populations during isolation followed by secondary contacts after the re-opening of the Strait of Gibraltar, accompanied by a reduced gene flow between the geographically divergent populations, could therefore have led to the present occurrence of population substructuring observed in some fish species across the Mediterranean and Atlantic waters. Further, the Almeria-Oran oceanographic front and the occurrence of a well-defined hydrographic barrier have been suggested to be responsible for its position as the boundary between Mediterranean and Atlantic populations of several marine organisms (Quesada et al. Reference Quesada, Beynon and Skibinski1995; Pannacciulli et al. Reference Pannacciulli, Bishop and Hawkins1997). Thus, the Strait of Gibraltar is not regarded as the main barrier between the Mediterranean and the Atlantic marine populations. The transport of North Atlantic water through a discontinuous one-way surface current flow throughout the Strait of Gibraltar into the Alboran Sea is also documented (Tintore et al. Reference Tintore, La Violette, Blade and Cruzado1988; Patarnello et al. Reference Patarnello, Volckaert and Castilho2007).

Three main hypotheses may explain the subpopulation structuring in marine pelagic and demersal fish of the Mediterranean fish resources: (i) environmental factors, including past sea-level changes, past barriers and present oceanic barriers, responsible for the mixing or interrupted gene flow between populations; (ii) increasing of geographic distances; (iii) life-history traits, including potential for larval dispersal, differential homing and spawning area behaviour, and a possible abiotic factor such as oceanic currents responsible for larval retention in coastal zones.

Generally speaking, transition zones are expected to play an important role in limiting the gene flow among populations, thus promoting their genetic substructuring. The existence of the Atlantic Ocean–Mediterranean Sea transition has been shown for several fish species such as the swordfish, X. gladius (Kotoulas et al. Reference Kotoulas, Magoulas, Tsimenides and Zouros1995; Chow et al. Reference Chow, Okamoto, Uozumi, Takeuchi and Takeyama1997, Reference Chow, Clarke, Nakadate and Okazaki2007; Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998; Alvarado-Bremer et al. Reference Alvarado Bremer, Mejutoc, Gómez-Márquezd, Boánd, Carpinterod, Rodríguez, Viñase, Greig and Elyg2005a , Reference Alvarado Bremer, Viñasa, Mejutoc, Elyd and Plae b ; Kasapidis et al. Reference Kasapidis, Mejuto, Tserpes, Antoniou, Garcia-Cortes, Peristeraki, Oikioonomaki, Kotoulas and Magoulas2007), common dentex, Dentex dentex (Bargelloni et al. Reference Bargelloni, Alarcon, Alvarez, Penzo, Magoulas, Reis and Patarnello2003), mullets Mugil cephalus (Crosetti et al. Reference Crosetti, Nelson and Avise1994), sand goby, Pomatoschistus minutus (Stefanni and Thorley, Reference Stefanni and Thorley2003), Mediterranean sea bass, Dicentrarchus labrax (Naciri et al. Reference Naciri, Lemaire, Borsa and Bonhomme1999; Bahri-Sfar et al. Reference Bahri-Sfar, Lemaire, Hassine and Bonhomme2000) and poor cod, Trisopterus minutus capelanus (Mattiangeli et al. Reference Mattiangeli, Ryan, Galvin, Mork and Cross2003). Moreover, an east-west fragmentation around the Siculo-Tunisian strait has been suggested in D. labrax (Barhi-Sfar et al. Reference Bahri-Sfar, Lemaire, Hassine and Bonhomme2000) and in X. gladius (Kasapidis et al. Reference Kasapidis, Mejuto, Tserpes, Antoniou, Garcia-Cortes, Peristeraki, Oikioonomaki, Kotoulas and Magoulas2007). Similarly, differences in specific hydrological, topographical and oceanographic features (depth, temperature, salinity) between the northern and southern part of the Adriatic Sea seem to be responsible for the genetic isolation of Adriatic sardine populations from the rest of the Mediterranean Sea (Alegria-Hernandez, Reference Alegria-Hernandez1986); these findings do not seem to be supported by the application of some genetic markers (Tinti et al. Reference Tinti, Di Nunno, Guarniero, Talenti, Tommasini, Fabbri and Piccinetti2002), despite a subpopulation structuring suggested recently based on the microsatellites loci analysis between Adriatic and Ionian Sea sardines (Ruggeri et al. Reference Ruggeri, Splendiani, Bonanomi, Arneri, Cingolani, Santojanni, Colella, Donato, Giovannotti, Caputo and Barucchi2013).

Most of the information regarding stock identification has been collected by applying morphological characteristics, morphometric traits, artificial tags, otolith shape and genetic methodologies. Indeed, in the last 15 years, genetic/molecular markers have become the most favoured method to assess fish population stock identity. Fishery scientists have applied many different types of genetic/molecular markers in the study of Mediterranean fish population structure, such as allozymes (Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998, Reference Cimmaruta, Bondanelli and Nascetti2005, Reference Cimmaruta, Bondanelli, Ruggi and Nascetti2008), mitochondrial DNA loci (Tinti et al. Reference Tinti, Di Nunno, Guarniero, Talenti, Tommasini, Fabbri and Piccinetti2002) and, more recently, microsatellites loci (Ruggeri et al. Reference Ruggeri, Splendiani, Bonanomi, Arneri, Cingolani, Santojanni, Colella, Donato, Giovannotti, Caputo and Barucchi2013, and reference therein) and single nucleotide polymorphisms (SNPs) (Milano et al. Reference Milano, Babbucci, Cariani, Atanassova, Bekkevold, Carvalho, Espiñeira, Fiorentino, Garofalo, Geffen, Hansen, Helyar, Nielsen, Ogden, Patarnello, Stagioni, Tinti and Bargelloni2014, and references therein). These are the most favoured markers to establish fish stock structure, to assess their genetic differentiation and to estimate their genetic variability levels. Genetic combined with biological studies have so far been carried out for the stock definition and assessment of anchovy (Bembo et al. Reference Bembo, Carvalho, Snow, Cingolani and Pitcher1995, Reference Bembo, Carvalho, Cingolani, Arneri, Giannetti and Pitcher1996a , Reference Bembo, Carvalho, Cingolani and Pitcher b ; Pla et al. Reference Pla, Tudela, Garcia-Marin, Palomera and Rubies1996; Tudela et al. Reference Tudela, Garcia-Marin and Pla1999), swordfish (Kotoulas et al. Reference Kotoulas, Magoulas, Tsimenides and Zouros1995; Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998; Kasapidis et al. Reference Kasapidis, Mejuto, Tserpes, Antoniou, Garcia-Cortes, Peristeraki, Oikioonomaki, Kotoulas and Magoulas2007), hake (Roldan et al. Reference Roldan, Garcia-Marin, Utter and Pla1998; Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005; Milano et al. Reference Milano, Babbucci, Cariani, Atanassova, Bekkevold, Carvalho, Espiñeira, Fiorentino, Garofalo, Geffen, Hansen, Helyar, Nielsen, Ogden, Patarnello, Stagioni, Tinti and Bargelloni2014); mackerel (Zardoya et al. Reference Zardoya, Castilho, Grande, Favre-Krey, Caetano, Marcato, Krey and Patarnello2004), sardines (Ruggeri et al. Reference Ruggeri, Splendiani, Bonanomi, Arneri, Cingolani, Santojanni, Colella, Donato, Giovannotti, Caputo and Barucchi2013) and mullets (Crosetti et al. Reference Crosetti, Nelson and Avise1994).

The literature concerning the use of parasites as biomarkers in the fish stock assessment of Mediterranean fish resources is not large. Among the methodologies applied in this field, the morphological identification of the parasite species, the use of genetic/molecular markers for the parasites’ detection, combined with the estimation of parasitic infection levels and multivariate analyses of the parasite community composition, are those mainly used so far (MacKenzie et al. Reference MacKenzie, Campbell, Mattiucci, Ramos, Pinto and Abaunza2008).

Power et al. (Reference Power, Balbuena and Raga2005) used parasite infracommunities as predictors of harvest location of bogue, Boops boops, in two fisheries off the Atlantic coast of Spain and one off the Spanish Mediterranean coast. Two parametric methods of classification were compared with other non-parametric methods to assign individual fish to their fishery ground. Five predictor parasite species, including the monogenean Microcotyle erythrini and the digeneans Aphanurus stossichii, Bacciger israelensis, Hemiurus communis and Lecithocladium excisum, present at high abundance in the whole host population, were considered as suitable biomarkers for the classification of bogues collected from different geographical areas of the Atlantic vs the western Mediterranean Sea (Spanish coast).

Similar studies have been carried out on parasite distribution in the striped red mullet, Mullus surmuletus, one of the most abundant and widely distributed fish in the coastal zone of the Mediterranean Sea and along the sublittoral zone of the North East Atlantic Ocean. Prevalence and abundance of three long-lived parasite taxa differed significantly among three localities off the Spanish Mediterranean coast (Ferrer-Castello et al. Reference Ferrer-Castelló, Raga and Aznar2007), indicating their potential as stock indicators of the striped red mullet. A cluster analysis based on prevalence and MANOVA on abundance datasets indicated strong inter-sample variability, even within the same locality, with poor spatial segregation among samples. A linear discriminant analysis (LDA) based on the abundance of 17 parasite taxa correctly assigned over 80% of fish to their capture locality, indicating good and stable predictive power based on the parasitological data (Ferrer-Castello et al. Reference Ferrer-Castelló, Raga and Aznar2007).

Klimpel et al. (Reference Klimpel, Kleinertz and Palm2008) found regional differences in the distribution of the parasite taxa identified in fish host populations of M. surmuletus sampled from the Mediterranean Sea and along the sublittoral zone of the North East Atlantic Ocean, suggesting the possible role of the history of the Mediterranean Sea as an hypothesis to explain the differential distribution of the parasite species detected.

Regional differences in infection levels by larval anisakids of the genus Anisakis from anchovies, Engraulis encrasicolus, fished in different fishing grounds of the Mediterranean Sea have been suggested to be related to the existence of different populations of this fish along its Mediterranean distribution (Rello et al. Reference Rello, Adroher, Benitez and Valero2009). The authors found the prevalence of the infection by Anisakis spp. in anchovies from the Ligurian Sea to be significantly higher than that in three other fishing areas of the southwestern Mediterranean Sea (Gulf of Lion, Catalonia coast and Alboran Sea) (Rello et al. Reference Rello, Adroher, Benitez and Valero2009). In similar studies, no infection by Anisakis was found in anchovies fished from off the Tarragona (NE Spain) (Serracca et al. Reference Serracca, Cencetti, Battistini, Rossini, Prearo, Pavoletti, Fioravanti, Righetti, Di Donfrancesco and Ercolini2013). On the other hand, anchovies caught in the Adriatic Sea off the Croatian coast showed moderately higher infection levels by Anisakis spp. larvae with respect to other Mediterranean and Atlantic waters which have shown strikingly lower prevalences (Mladineo et al. Reference Mladineo, Simat, Miletic, Beck and Poliack2012). However, when considering the host population dynamics as having a possible role in infection by larval Anisakis spp., the fish host size has not always been taken into account to explain regional differences observed in parasitic infection values from the same fish species.

When carrying out a multimethodological approach to fish stock identification, including morphological and biological traits (including parasites) and the genetic structure of the fish host, it is fundamental to include fish samples from the whole range of the species and to use the same fish samples for all the methodologies applied (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ). Furthermore, while fish size is not a confounding variable in detecting genetic population structure, it is in parasite studies, especially in long-lived species such as larval Anisakis, where it could have a confounding effect. This explains why, whenever possible, the sample size of the fish host should be as homogenous as possible and comparable for all the methodological approaches. However, it could also happen that when comparing fish population structure from different basin waters such as the Mediterranean Sea and the Atlantic Ocean, the size of the fish hosts, despite being the same age, may be different due to other environmental variables (abiotic factors, food availability and overfishing). This is for instance the case with the Mediterranean and Atlantic populations of swordfish (Ehrhardt, Reference Ehrhardt1992; Tserpes and Tsimenides, Reference Tserpes and Tsimenides1995).

Furthermore, when choosing a parasite as the best suitable biomarker, in addition to the ecological characteristics of the parasite it is also important to apply the most suitable method for its identification, especially when the parasite includes a complex of cryptic (or sibling) species. This is, for instance, the case for larval nematodes of the Anisakis simplex (sensu lato) complex, which comprises so far three biological species, named as Anisakis pegreffii, A. simplex (sensu stricto) and A. berlandi (see Mattiucci et al. Reference Mattiucci, Nascetti, Cianchi, Paggi, Arduino, Margolis, Brattey, Webb, D'Amelio, Orecchia and Bullini1997, Reference Mattiucci, Cipriani, Webb, Paoletti, Marcer, Bellisario, Gibson and Nascetti2014a ). The third-stage larvae of these three species, recovered from the fish hosts, show the morphotype indicated as Anisakis type I or Anisakis type II (Berland, Reference Berland1961), but they cannot be identified to their species level based on traditional morphological analysis. Among the genetic methodologies used for the species identification of Anisakis spp. larvae, the allozyme markers, as based on several diagnostic loci among the species of Anisakis so far genetically detected, are able to recognize the three species at any of their life-history stages. In addition, because of the temporal and geographic stability of the allozyme data, they are used in assessing the distribution of the different recognized species of Anisakis in the same fish host captured from several geographic areas. Finally, in order to have a comprehensive pattern of distribution of a parasite species in a fish species along its geographic range, a large number of parasites (up to several thousand) must be identified by molecular/genetic methodologies to their species level. This was, for instance, the case with the larval Anisakis spp. identified in samples of M. merluccius (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004), Trachurus trachurus (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008) and X. gladius (Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ). In addition, among the other molecular markers used for Anisakis detection (Mattiucci et al. Reference Mattiucci, Cipriani, Webb, Paoletti, Marcer, Bellisario, Gibson and Nascetti2014a ), the high substitution rate of the mtDNA cox2 sequences found so far in Anisakis spp. (Mattiucci et al. Reference Mattiucci, Cipriani, Webb, Paoletti, Marcer, Bellisario, Gibson and Nascetti2014a ) also allows the use of the phylogeography of Anisakis spp. to provide support to the subpopulation structure, as inferred for its fish host.

In the following paragraphs particular attention is focused on reviewing the use of some parasite species (i.e. larval nematode of the genus Anisakis) as biomarkers in three case studies of populations of large, small, pelagic and demersal fish species of the Mediterranean Sea within a multidisciplinary approach, including the use of genetic/molecular markers to detect fish host structure. The parasitological datasets concerning the occurrence of Anisakis spp. larvae are reviewed here in the context of the ‘fish stock identification’ as based on a ‘holistic approach’. Particular attention is then given to the results achieved by us, in recent years, on the population genetic structure of three fish species (European hake M. merluccius, horse mackerel T. trachurus and swordfish X. gladius) as inferred from allozyme markers, and on some parasite species (larval Anisakis spp.) which are the most suitable biological tags to characterize fish stocks.

Meta-analysis of parasites in their fish hosts was used to correlate the genetic datasets on the fish host and larval Anisakis as biological markers, as obtained on the same samples of the pelagic and demersal fish species collected in localities within the range of distribution of the three species in Mediterranean and Atlantic waters.

Congruence and incongruence in depicting the stocks of those fish species according to the two different datasets, i.e. fish genetics and Anisakis spp. larval distribution in fish populations, are analysed by a generalized Procrustes Rotation (PR) (also known as analysis of congruence) and discussed in the light of the biological features of the species. Indeed, Procrustes Analysis has proven to be a useful approach to test the congruence between two or more multivariate datasets, such as in the case of application of co-phylogenetic analysis of host-parasites associations (Balbuena et al. Reference Balbuena, Miguez-Lozano and Blasco-Costa2013), and genetic and morphological datasets achived on parasite species (Mattiucci et al. Reference Mattiucci, Cipriani, Webb, Paoletti, Marcer, Bellisario, Gibson and Nascetti2014a ).

In this paper, a generalized PR, performed by the Software R (R Development Core Team, 2012) was used to compare different ordinations of specimens of the three fish species, object of the review, based on combined genetic data of the fish host and Anisakis spp. larvae distribution and burden distance matrices. The latter was computed by considering the number of individual parasites (fourth-root transformed to improve normality and remove the mean/variance relationship) belonging to different species of Anisakis, genetically identified. A Procrustes test with n = 1000 bootstrap replicates was used to estimate the significance of the Procrustes statistics, assessing similarities between different ordinations, and the Procrustes correlation ‘r ’ was derived from the symmetric Procrustes residual (Peres-Neto and Jackson, Reference Peres-Neto and Jackson2001).

CASE STUDIES OF FISH STOCK IDENTIFICATION IN THE MEDITERRANEAN SEA

Demersal fish – the case of the European hake, Merluccius merluccius

The European hake is a demersal fish living in eastern Atlantic waters from the North Sea to the African coasts, including the Mediterranean and Black Sea. Hake fisheries are well developed throughout the whole species range and the available data suggest that at least the Mediterranean population could be overexploited (GFCM, 2012).

Here we present data acquired on the pattern of distribution and abundance of different species of Anisakis, supporting the existence of different stocks of M. merluccius sampled from several areas covering its whole geographic range. The data were acquired in the framework of a multidisciplinary international project GENHAKE (FAIR CT97-3494), which was aimed at clarifying the stock structure of M. merluccius in European waters. The project included the analysis of biological and morphometric characteristics of the fish host, and the use of different nuclear genetic/molecular markers to define the stock structure. The parasitological examination was carried out on the same individuals used for the genetic analyses and other biological studies.

Among the wide range of parasite taxa found in M. merluccius across its range of distribution, attention was focused on some selected parasites proved to be of particular value in stock definition. These are the larval nematodes of the genus Anisakis: they are long-lived parasites, occurring at high frequency in this fish species throughout its geographic range and they are not pathogenic to the fish host. Therefore, they comply with the main characteristics defining a parasite as a good ‘biomarker’. The relative proportions of larvae of different species of Anisakis in different sampling localities of M. merluccius allowed the recognition of Mediterranean and Atlantic stocks of European hake (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004) (Fig. 1A). Anisakis  pegreffii was identified in all the hakes sampled from the Mediterranean Sea, except in the Levantine Sea (off the Cyprus coast). In the western Mediterranean Sea, A. pegreffii was the most prevalent and abundant species in the Ligurian Sea; this finding seems to be related to the fact that this area is comprised in the Protected Marine Area the Pelagos Sanctuary of Mediterranean Cetaceans, a Mediterranean protected area inhabited by various dolphin species (Notarbartolo di Sciara et al. Reference Notarbartolo di Sciara, Agardy, Hyrenbach, Scovazzi and Van Klaveren2008), which are the main definitive hosts of A. pegreffii (Mattiucci and Nascetti, Reference Mattiucci and Nascetti2008; Mattiucci et al. Reference Mattiucci, Cipriani, Webb, Paoletti, Marcer, Bellisario, Gibson and Nascetti2014a ). This species also occurred at high prevalence and density in fish sampled from the eastern part of the Mediterranean Sea, such as in the Ionian, Aegean and Cretan Seas. Interestingly, hakes from Levantine Sea waters harboured only the species A. typica, a species parasitic in many cetacean species in warm temperatures and tropical waters (Mattiucci et al. Reference Mattiucci, Paggi, Nascetti, Portes Santos, Costa, Di Benedetto, Ramos, Argyrou, Cianchi and Bullini2002). Anisakis physeteris was the most prevalent species in hakes from the western Mediterranean Sea off the Balearic Islands, Alicante and Malaga (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004). The high proportion of this Anisakis species found in the fish from this deep area of the Mediterranean basin around the Balearic Islands could be related to the presence of the main migration routes for the sperm whale Physeter catodon, the definitive host of A. physeteris (Mattiucci and Nascetti, Reference Mattiucci and Nascetti2008). However, to date hakes from the western Mediterranean Sea were not found infected with A. simplex (s.s.) (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004), although this parasite species has been documented in some pelagic fish species, such as mackerel Scomber scombrus (Nascetti et al. Reference Nascetti, Paggi, Orecchia, Smith, Mattiucci and Bullini1986) and horse mackerel T. trachurus (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008; see next paragraph) from some of the collecting areas. This finding supports the hypothesis that there are no migrations of M. merluccius from and into Atlantic waters. Conversely, in the North-East Atlantic hake samples from north of the Strait of Gibraltar, A. simplex (s.s.) was the most prevalent species, while it occurred off the Spanish Atlantic coast in sympatry with A. pegreffii. Indeed, hakes that exhibited mixed infections of A. simplex (s.s.) and A. pegreffii represented >20% of the fish examined from the Atlantic coast of Galicia and 14% from the Bay of Biscay (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004). Moving south in the North-East Atlantic from the Strait of Gibraltar, mixed infections of different species of Anisakis in hakes caught along the Atlantic coast of Morocco were observed. More than 22% of the fish examined were found to be parasitized by five species of Anisakis: the major component species was A. pegreffii, followed by A. physeteris, A. ziphidarum, and, at a lower percentage, A. brevispiculata and A. paggiae (Fig. 1A). Thus, according to the pattern of distribution of larval Anisakis spp. and their abundance in the hake samples, different populations of M. merluccius in European waters (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004) were identified by PCA analysis (Fig. 1B). The larval distribution and abundance of the different species of Anisakis recognized in hakes from the different fishing grounds indicate that: (1) there are two stocks of M. merluccius, from Mediterranean and Atlantic waters respectively; (2) in the North Atlantic area, at least two distinct subpopulations are present, one north of the Strait of Gibraltar (samples 3, 6, 7 Fig. 1B) and another (no. 11, Fig. 1B) from off the Atlantic coast of Morocco; (3) some substructuring of the western and eastern parts of the Mediterranean hake populations seems to be recognized, with the population from Levantine Sea distinct from the other Mediterranean ones (Fig. 1B) (Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004, Reference Mattiucci, Abaunza, Damiano, Garcia, Santos and Nascetti2007). Interestingly, in the PCA analysis, the sample from the Alboran Sea was more similar to Atlantic populations (Fig. 1B).

Fig. 1. The case study of stock identification of a demersal fish species, the European hake Merluccius merluccius, from the Mediterranean Sea and Atlantic Ocean, as inferred from biological data, biomarkers (Anisakis spp. parasites) and the population genetics data (allozymes) of the fish populations (data from Nascetti et al. Reference Nascetti, Abaunza, Mattiucci, Magoulas, Cimmaruta and Bullini2000; Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004; Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005). (A) Relative proportions of the Anisakis spp. larvae genetically identified (N parasites = 1950) in M. merluccius from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Mattiucci et al. Reference Mattiucci, Abaunza, Ramadori and Nascetti2004), plotted with the stock units as identified by the biological and morphometric analyses (data from Nascetti et al. Reference Nascetti, Abaunza, Mattiucci, Magoulas, Cimmaruta and Bullini2000). A different colour identifies a putative stock unit as defined by morphometric datasets of the fish host. Dotted lines identify the boundaries between hake stock units, as identified by genetic/molecular markers; (B) Principal Component Analysis (PCA) inferred from the larval Anisakis spp. bio-markers, identified genetically in the European hake specimens (M. merluccius sampled in the Mediterranean (white circles) and Atlantic (black circles) basins). The first axis of PCA explains about 60% of total ordination, and PCA 2 about 30%; (C) A generalized Procrustes Rotation (PR) performed with R (R Development Core Team, 2012) was used to compare different ordinations of the European hakes M. merluccius (N fish = 1306), based on combined molecular and parasite datasets. The latter was computed by considering the number (N parasites = 1950) of larval parasites belonging to A. pegreffii, A. simplex (s. s.), A. physeteris, A. typica, A. ziphidarum and A. brevispiculata. A Procrustes test (also known as analysis of congruence) with n = 1000 bootstrap replicates was used to estimate the significance of the Procrustes statistics, assessing similarities between different ordinations, and the Procrustes correlation r was derived from the symmetric Procrustes residual (Peres-Neto and Jackson, Reference Peres-Neto and Jackson2001). Analysis showed a significant correlation (Procrustes r = 0·59, P<0·01) between the Anisakis spp. burden and molecular traits, showing the clustering of sampled specimens in two well-defined clusters corresponding to the Mediterranean (white circles) and Atlantic stocks (black circles). Code numbers and coloured ellipses are related to the sampling localities, as appearing in Fig. 1A.

The parasitological findings are in agreement with the bulk of biological data, showing separate Mediterranean and Atlantic stocks of European hake since they are known to live in markedly different environments and to differ in many biological and demographic features such as growth rate, spawning season and recruitment (Froese and Pauly, Reference Froese and Pauly2014).

In addition, the parasitological data were largely in accordance with the genetic results obtained by using different molecular nuclear markers (allozymes, microsatellites, SNPs), showing a main subdivision between Atlantic and Mediterranean stocks (Roldan et al. Reference Roldan, Garcia-Marin, Utter and Pla1998; Lundy et al. Reference Lundy, Moran, Rico, Milner and Hewitt1999, Reference Lundy, Rico and Hewitt2000; Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005; Milano et al. Reference Milano, Babbucci, Cariani, Atanassova, Bekkevold, Carvalho, Espiñeira, Fiorentino, Garofalo, Geffen, Hansen, Helyar, Nielsen, Ogden, Patarnello, Stagioni, Tinti and Bargelloni2014). However, the molecular markers indicated that the boundary between the two stocks is located along the Almeria-Oran Front (AOF), instead of the Strait of Gibraltar. Within the Atlantic basin, the hake genetic structure is coincident with the Anisakis findings. However, the management areas of the International Council for the Exploration of the Sea (ICES) for European hake are not coincident with the parasitological analysis (Fig. 1C). The ICES distinguishes a northern and a southern stock of European hake with a boundary limit in the Cap Breton Canyon (Southeastern of the Bay of Biscay). This definition is based more on administrative reasons rather than on biological supporting, although the available data on morphometrics and biology seem to support a substructure in European hake from Atlantic waters to the North of the Strait of Gibraltar (Fig. 1A). The Moroccan area is out of the ICES area, and based on our parasitological data, the Northwest African waters could be a candidate for a southern stock of European hake, differentiated from the northern stock which inhabits European Atlantic waters (Fig. 1B and C).

Both allozymes and SNPs clearly identified a prominent role of selection in moulding the genetic stock structure between and within sea basins (Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005; Milano et al. Reference Milano, Babbucci, Cariani, Atanassova, Bekkevold, Carvalho, Espiñeira, Fiorentino, Garofalo, Geffen, Hansen, Helyar, Nielsen, Ogden, Patarnello, Stagioni, Tinti and Bargelloni2014), indicating the sea water temperature and salinity as the main driving factors shaping the genetic structure of a number of outlier loci. This finding suggests that the European hake represents a case study showing how marine environmental variables, such as temperature and salinity, are correlated to both the genetic structure of the fish host species, and to the pattern of distribution of the larval parasite Anisakis spp.

The high agreement among the different datasets in identifying the subpopulations of European hake was supplemented with a generalized Procrustes Rotation (Fig. 1C), which compared different ordinations of the European hake specimens as inferred from the combination of the genetic analysis (allozymes, data from Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005) with the larval distribution of the different Anisakis species, as detected in the same subpopulations. The analysis showed a significant correlation (Procrustes r = 0·59, P<0·01) between the distribution and abundance of the larval Anisakis and the pattern of fish-host genetic structure. The sampled specimens clustered in two well-defined groups corresponding to the Mediterranean (white circles) and Atlantic stocks (black circles) (Fig. 1C). Interestingly, in the PR analysis, the population from the Alboran Sea (no. 18, Fig. 1C) showed more similarity to the Atlantic samples than to those from the Mediterranean basin. In addition, the eastern population from the Levantine Sea (no. 22, Fig. 1C) resulted as being distinct from the western Mediterranean populations (Fig. 1C). Both datasets contributed strongly to the ordination recorded, with the molecular distance between samples accounting for about 93% of total ordination and the distance based on the parasitic burden of Anisakis spp. larvae for about 75%.

The results achieved in the case of European hake are evidence of how Anisakis spp. larvae are good biological markers in recognizing fish populations from different fishing grounds, given the assumption that these parasites are related to differences in life-history traits of their individual fish host. Indeed, the ecological differences of Anisakis spp., including their biogeography and life cycles, reflect the environmental conditions and feeding behaviour of adult fish related to each single fishing ground in the Atlantic and Mediterranean Sea waters.

Small pelagic fish – the case of the horse mackerel, Trachurus trachurus

The horse mackerel, T. trachurus is a small carangid fish living on the continental shelf of the North-Eastern Atlantic Ocean, from Norway to the Cape Verde Islands, and in the Mediterranean Sea. Horse mackerel catches have a high relevance in fisheries, with landings from the Atlantic area strongly decreasing after 2000 and now partially recovering (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ; Fernandes and Cook, Reference Fernandes and Cook2013; ICES, 2013). The horse mackerel is characterized by a high migratory ability and by a reproductive strategy with a rather extensive spawning season (4–8 months, according to latitude) and an asynchronous maturation of oocytes that mature and are spawned in a number of consecutive clutches (Abaunza et al. Reference Abaunza, Gordo, Karlou-Riga, Murta, Eltink, Garcia Santamaria, Zimmermann, Hammer, Lucio, Iversen, Molloy and Gallo2003; Murta et al. Reference Murta, Pinto and Abaunza2008). The combination of high vagility and multiple spawning makes the boundaries of horse mackerel stocks very difficult to define, while its economic relevance makes this issue a crucial one.

In recent years, the horse mackerel has been the target of a wide multi-methodological research, aimed at identifying the biological stocks of horse mackerel by simultaneously studying subpopulations by means of different genetic/molecular markers, morphometry, life-history traits and biological parameters, including parasites (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a , Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Gallo, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Ruggi, Sanjuan, Santos, Stransky and Zimmerman b ). All the genetic/molecular datasets were concordant in describing a low level of differentiation among the Mediterranean and Atlantic populations, whilst a high genetic variation was found at the intra-population level. Both nuclear (allozymes and DNA microsatellites loci) and mitochondrial (ND4L and ND2 sequences analysis of mtDNA) markers found the larger part of the variation to be distributed within the horse mackerel populations sampled from the different fishing grounds rather than between them (Cimmaruta et al. Reference Cimmaruta, Bondanelli, Ruggi and Nascetti2008; Comesaña et al. Reference Comesaña, Martinez-Areal and Sanjuan2008; Kasapidis and Magoulas, Reference Kasapidis and Magoulas2008).

The allozyme analysis suggested an association between the levels of genetic differentiation and the geographic distance in the Atlantic Ocean (Isolation By Distance, IBD), although it was impossible to distinguish discrete stocks on the basis of genetic data alone. These findings are in agreement with the biological features of the species, a vagile batch spawner, supporting the idea that the low genetic structuring of the fish species is related to a high gene flow due to the fish behaviour and mobility, rather than to the passive dispersion driven by marine gyres during the larval stage (Cimmaruta et al. Reference Cimmaruta, Bondanelli, Ruggi and Nascetti2008).

Despite the low resolution power of genetic/molecular markers to clearly distinguish the stocks of horse mackerel in European waters, several biological approaches (morphometrics, life-history traits, otoliths and parasites) were able to detect discrete populations (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ). Indeed, the above-mentioned techniques separated a Mediterranean stock from the Atlantic ones. A further subdivision in the NE Atlantic was observed: a ‘North Sea’, a ‘Western’ and a ‘Southern’ stock were identified. Also in the Mediterranean Sea, a western, a central and an eastern subpopulation were outlined (Fig. 2A). In addition, the biological data acquired indicated an influence of the Atlantic populations on the western extreme (Alboran Sea) of the Mediterranean Sea (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a , Reference Abaunza, Gordo, García Santamaría, Iversen, Murta and Gallo c ; Gordo et al. Reference Gordo, Costa, Abaunza, Lucio, Eltink and Figueiredo2008; MacKenzie et al. Reference MacKenzie, Campbell, Mattiucci, Ramos, Pinto and Abaunza2008; Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008; Murta et al. Reference Murta, Pinto and Abaunza2008; Stransky et al. Reference Stransky, Murta, Schlickeisen and Zimmermann2008).

Fig. 2. The case study of stock assessment of a small pelagic fish species, the horse mackerel Trachurus trachurus, from the Mediterranean Sea and Atlantic Ocean waters, as inferred from biological data, biomarkers (parasites), and the population genetics dataset (allozymes) of the fish populations (data from Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ; Cimmaruta et al. Reference Cimmaruta, Bondanelli, Ruggi and Nascetti2008; Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008; Murta et al. Reference Murta, Pinto and Abaunza2008). (A) Relative proportions of the Anisakis spp. identified genetically (N parasites = 2200) in T. trachurus from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008), plotted with the stock units identified by the biological and morphometric analyses (data from Murta et al. Reference Murta, Pinto and Abaunza2008). Different colours identify each putative stock unit as defined by morphometric datasets of the fish host; (B) Principal Component Analysis (PCA) inferred from larval Anisakis spp. as biomarkers, identified genetically in the horse-mackerel specimens (T. trachurus) sampled in the Mediterranean (white circles) and Atlantic (black circles). The first axis of PCA explains about 92% of total ordination, whilst PCA 2 about 7%; (C) A generalized Procrustes Rotation (PR), performed with R (R Development Core Team, 2012) used to compare different ordinations of the Atlantic horse mackerels (T. trachurus) analysed (N fish = 2241), based on combined molecular and Anisakis spp. burden distance matrices. The latter was computed by considering the number of individual parasites (N parasites = 2200) (fourth-root transformed to improve normality and remove the mean/variance relationship) belonging to A. pegreffii, A. simplex (s.s.), A. physeteris and A. typica. A Procrustes test (also known as analysis of congruence) with n = 1000 bootstrap replicates was used to estimate the significance of the Procrustes statistics, assessing similarities between different ordinations, and the Procrustes correlation r was derived from the symmetric Procrustes residual (Peres-Neto and Jackson, Reference Peres-Neto and Jackson2001). Analysis showed significant correlation (Procrustes r = 0·39, P = 0·05) between the parasitic burden and molecular traits, as shown by the clustering of sampled specimens corresponding to the Mediterranean (white circles) and Atlantic stocks (black circles). Molecular distance between specimens accounted for about 37% of total ordination, whilst the distance based on the parasitic burden with Anisakis for about 95%. Code numbers and coloured ellipses are related to the sampling localities, as appearing in Fig. 2A.

The use of Anisakis parasites as biological tags gave concordant results with otoliths and morphometric data analysis (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008), rather than with the genetic datasets.

Of the five species of Anisakis identified in horse mackerel, A. pegreffii and A. simplex (s.s.) were found to be the dominant species in all the fish sampled along its geographic distribution, whereas the other three species, i.e. A. physeteris, A. typica and A. nascettii (the last indicated as Anisakis sp. in Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008), were identified at very low percentages in the central Tyrrhenian Sea, the Cretan Sea and in the NE Atlantic fish samples (Fig. 2A). Anisakis pegreffii and A. simplex (s.s.) showed statistically significant differences in their relative proportions in the horse mackerel samples throughout the NE Atlantic and the Mediterranean Sea. Indeed, A. pegreffii was identified as the main species parasitizing horse mackerel in the Mediterranean Sea; it was the only species found in the Adriatic and Ionian Seas samples, while in the horse mackerels from the Aegean Sea, A. pegreffii occurred in mixed infections with A. typica (1%). Similarly, in the western Mediterranean Sea, A. pegreffii was the dominant species, even though A. physeteris occurred in mixed infections in horse mackerel fished along the Tyrrhenian Sea coast (Fig. 2A). Mattiucci et al. (Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008) found that the samples of horse mackerel fished in the Alboran Sea showed mixed infections by both the sibling species, A. pegreffii and A. simplex (s.s.) in almost equal proportions (Fig. 2A). Moving to the NE Atlantic Sea, the occurrence of A. pegreffii progressively decreased from almost 87% off the Portuguese coast of Algarve, to 30% in the horse mackerels from off the Spanish Galician coast. Anisakis pegreffii was also rarely identified from horse mackerel from off the south coast of Ireland and from the Bay of Biscay (Fig. 2A). The opposite trend was found in the occurrence of A. simplex (s.s.), which showed increasing proportions in the horse mackerels fished from off the Portuguese coast (almost 13%) to the Spanish Atlantic coast of Galicia (71%). It became the only species present in the fish species in the Norwegian Sea and along the western English Channel (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008) (Fig. 2A).

In horse mackerel fished along the Portuguese and Spanish Atlantic coast, a contact zone between A. pegreffii and A. simplex (s.s.), the two parasite species occurred at high percentages in mixed infections, with the highest level in 43% of the fish captured from off the North Portuguese coast (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008). Similarly, in the Alboran Sea, the two species were found in syntopy in almost 60% of the infected fish, while most (35%) of the remaining fish examined from this basin were found infected only by the species A. pegreffii (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008). Further, few (N = 15 out of 2200 larvae analysed) F1 hybrid individuals between the two species A. pegreffii and A. simplex (s.s.) – i.e. showing a heterozygote genotype at all the diagnostic allozyme loci between the two species – were recognized in the sympatric areas of the two parasite species (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008 and unpublished data).

According to the genetically recognized species of Anisakis in the horse mackerel (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008) and from abundance levels of infection, different populations of the fish species can be distinguished, as shown in Fig 2B. Indeed, a discrete structuring of subpopulations of horse mackerel from the Mediterranean Sea (excluding the Alboran Sea) and those from the Atlantic can be recognized according to PCA analysis (Fig. 2B). Moreover, the same analysis seems to support the existence of a ‘southern’ stock of horse mackerel along the Spanish and Portuguese coasts (including the samples no. 7, 8, 9, 10), distinct from the ‘northern’ stock, as represented by the areas of the Northern North Sea and Southern Ireland coast (samples no. 1, 2, 3, 4, 5, 6) (Fig. 2A and B). Finally, it seems that the sample from the Alboran Sea (no. 18, Fig. 2A and B), according to its larval Anisakis spp. composition, clusters with the other Atlantic samples rather than with those from the Mediterranean Sea (Fig. 2B). This finding, and the presence of a low percentage of fish infected only by A. simplex (s.s.) in the Alboran Sea (Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008), suggests that this population of horse mackerel could be the result of a migration of this small pelagic fish species from Atlantic waters into this extreme western area of the Mediterranean Sea, possibly mixing with the other Mediterranean Sea populations. The migrations of Atlantic horse mackerel into the western part of the Mediterranean Sea is also supported by the finding of the monogenean Heteraxinoides atlanticus in the same fish sample examined from the Alboran Sea (MacKenzie et al. Reference MacKenzie, Campbell, Mattiucci, Ramos, Pinto and Abaunza2008), a monogenean species recognized in the Atlantic horse mackerel. Thus, in turn, the finding that the fish samples from the Alboran Sea are more similar, in their parasite fauna, to the Atlantic populations rather than to the Mediterranean ones, is supported by the existence of the oceanographic front, a particular water circulation, called the AOF. This is an oceanographic front located from Almeria (Spain) to Oran (Morocco), consisting of inflowing Atlantic water which, entering at the surface through the Gibraltar Strait, describes in the Alboran Sea a quasi-permanent anticyclonic gyre in the west with less variation in the eastern part. This oceanographic front exhibits a pronounced step temperature (1·4 °C) and salinity gradient (2 ppt) over a distance of 2 km with an average current speed of 40 cm s−1 flowing southeastward from the Spanish coast to the coast of North Africa (Tintore et al. Reference Tintore, La Violette, Blade and Cruzado1988). This water circulation could explain the finding of A. simplex (s.s.), a parasite species widespread between 35°N and the Arctic Circle, in horse mackerel fished from the Alboran Sea. These horse mackerels could have been infected by preying upon infected invertebrates and/or smaller pelagic fish species in the Alboran basin.

Thus, horse mackerel shows a low level of genetic differentiation between the Mediterranean and Atlantic populations, a stable genetic structure and a high level of genetic variability without apparent geographic structuring (Cimmaruta et al. Reference Cimmaruta, Bondanelli, Ruggi and Nascetti2008; Comesaña et al. Reference Comesaña, Martinez-Areal and Sanjuan2008; Kasapidis and Magoulas, Reference Kasapidis and Magoulas2008). However, several other approaches (morphometrics and parasites) support the separation between the Mediterranean and Atlantic populations of horse mackerel whilst parasites, morphometrics and life-history traits support the existence of the ‘southern’ and ‘northern’ stocks in Atlantic waters (Murta et al. Reference Murta, Pinto and Abaunza2008), with the possible migration of adult fish following the west coast in the North-East Atlantic between the Celtic and Northern North Seas.

Finally, parasites and morphometric analyses were congruent in depicting the existence of sub-structured areas of horse mackerel in the Mediterranean Sea (western, central and eastern) (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a , Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Gallo, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Ruggi, Sanjuan, Santos, Stransky and Zimmerman b , Reference Abaunza, Gordo, García Santamaría, Iversen, Murta and Gallo c ; MacKenzie et al. Reference MacKenzie, Campbell, Mattiucci, Ramos, Pinto and Abaunza2008; Mattiucci et al. Reference Mattiucci, Farina, Campbell, Mackenzie, Ramos, Pinto, Abaunza and Nascetti2008; Murta et al. Reference Murta, Pinto and Abaunza2008) (Fig. 2C).

The Generalized PR comparing different ordinations of the horse mackerel populations based on combined molecular and parasitic burden datasets showed slightly significant correlation (Procrustes r = 0·39, P = 0·05) between the parasitic burden and molecular traits, as shown by the clustering of sampled specimens corresponding to the Mediterranean and Atlantic horse mackerel (Fig. 2C). Indeed, despite the fact that a trend similar to that shown by PCA of parasites (Fig. 2B) in depicting the existence of a stock structure of horse mackerel in European waters, fish populations genetic data accounted for only about 37% of total ordination at the PR analysis (Fig. 2C), whilst the differences in the parasitic burden by larval Anisakis spp. accounted for about 95%. This finding further supports the idea that parasites show structuring of their fish host populations at regional and temporal scales. Indeed, the latitudinal patterns of the larval Anisakis spp. distribution delineated the horse mackerel putative stocks in the Mediterranean and Atlantic waters according to the different regions of basin waters inhabited by different host meta-populations of the pelagic fish species.

Large pelagic fish – the case of the swordfish, Xiphias gladius

The correct identification of swordfish stocks is a ‘hot topic’ in fisheries science, being the rationale for the assessment and management of this important resource. As a consequence, different genetic/molecular techniques have been fruitfully used to study the stock structure of this cosmopolitan large pelagic fish that is characterized by a philopatric behaviour and a limited gene exchange between different parts of its range, in contrast to its remarkable migration ability. As a consequence of these characteristics, different stocks of swordfish corresponding to many reproductive areas are now recognized. In the Pacific Ocean, swordfish from north-eastern, north-central, south-eastern and south-western areas could be genetically differentiated (Alvarado-Bremer et al. Reference Alvarado Bremer, Hinton and Grieg2006). In the Atlantic Ocean, a ‘northern’ and a ‘southern’ Atlantic stock are now recognized, with their boundary still questioned based on genetic data: some markers confirm the 5°N position officially set by the International Commission for the Conservation of Atlantic Tunas (ICCAT) (Kasapidis et al. Reference Kasapidis, Mejuto, Tserpes, Antoniou, Garcia-Cortes, Peristeraki, Oikioonomaki, Kotoulas and Magoulas2007), while other data move it northward (Chow et al. Reference Chow, Clarke, Nakadate and Okazaki2007) or suggest an intergradation zone in the equatorial area (Smith and Alvarado-Bremer, Reference Smith and Alvarado-Bremer2010; Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011) (Fig. 3A). As to the Mediterranean population, it is considered as a further stock, highly isolated and divergent from the Atlantic one and characterized by small population size and low genetic variability (Smith and Kotoulas et al. Reference Kotoulas, Magoulas, Tsimenides and Zouros1995; Chow et al. Reference Chow, Okamoto, Uozumi, Takeuchi and Takeyama1997; Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998; Chow and Takeyama, Reference Chow and Takeyama2000; Vinas et al. Reference Viñas, Pérez-Serra, Vidal, Alvarado-Bremer and Pla2010). These features are explained by the history of this population, which was likely isolated and restricted in the eastern Mediterranean basin during the last two glacial periods and is still exchanging genes very little with the Atlantic populations (Alvarado-Bremer et al. Reference Alvarado Bremer, Viñasa, Mejutoc, Elyd and Plae2005b ). Indeed, Mediterranean swordfish are known to move to Atlantic waters via the Strait of Gibraltar for feeding, while Atlantic specimens may enter the Mediterranean Sea. However, the spawning grounds of the Mediterranean and Atlantic populations, and thus their gene pools, remain substantially separated due to the fish's philopatric behaviour (Magoulas et al. Reference Magoulas, Kotoulas, de la Serna, De Metrio, Tsimenides and Zouros1992; Chow et al. Reference Chow, Okamoto, Uozumi, Takeuchi and Takeyama1997). In the Mediterranean Sea, a genetic sub-structuring of swordfish has been revealed by phylogeographic analysis of mitochondrial control gene, supporting the existence of ‘eastern’ and ‘western’ subpopulation units (Viñas et al. Reference Viñas, Pérez-Serra, Vidal, Alvarado-Bremer and Pla2010).

Fig. 3. The case study of stock assessment of a large pelagic fish species, the swordfish Xiphias gladius from the Mediterranean Sea and Atlantic Ocean, as inferred from biological data, biomarkers (parasites) and the genetic datasets (allozymes) of the same fish populations (data from Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011; Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ). (A) Relative proportions of the Anisakis spp. identified genetically (N parasites = 405) in X. gladius from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011; Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ), plotted with the stock units as suggested by the genetic/molecular results (Kotoulas et al. Reference Kotoulas, Magoulas, Tsimenides and Zouros1995; Alvarado-Bremer et al. Reference Alvarado Bremer, Mejutoc, Gómez-Márquezd, Boánd, Carpinterod, Rodríguez, Viñase, Greig and Elyg2005a , Reference Alvarado Bremer, Viñasa, Mejutoc, Elyd and Plae b ; Chow et al. Reference Chow, Clarke, Nakadate and Okazaki2007). Dotted lines identify the boundaries between swordfish stock units; (B) Principal Component Analysis (PCA) inferred from the parasitic burden of the larval Anisakis spp. identified genetically and used as biomarkers, in the swordfish X. gladius specimens sampled in the Mediterranean (white circles) and Atlantic (black circles) basins. PCA analysis was computed by considering the number of individual parasites (fourth-root transformed to improve normality and remove the mean/variance relationship) belonging to A. pegreffii, A. simplex (s.s.), A. physeteris and A. brevispiculata. The first axis of PCA explains about 54% of total ordination, and PCA 2 about 22%; (C) Principal Component Analysis (PCA) inferred from Nei's distance matrix based on allozyme data (source: Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998) performed on swordfish specimens analysed genetically (N fish = 323). The first and second axis accounted for 96 and 2% of total ordination, respectively.

The phylogeography and the behaviour of swordfish accounts for its well-defined pattern of genetic variation, making the distinct populations/stocks well identifiable by using genetic/molecular markers on the fish.

According to the results achieved on the different species of larval Anisakis recognized in swordfish from the Mediterranean Sea and Atlantic Ocean fishing grounds, different populations of swordfish have been distinguished from the two basin waters (Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011; Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ). Indeed, among all the metazoan taxa of swordfish identified in several samples of the fish captured in Mediterranean Sea and Atlantic Ocean waters, the larval stages of Anisakis were among those mainly accounting for the differences detected between the Mediterranean and Atlantic populations (Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ). In fact, A. pegreffii was identified only in swordfish from the Mediterranean Sea fishing grounds, whereas the larvae of A. brevispiculata, A. paggiae, Anisakis sp. 2 and A. simplex (s.s.), were absent in X. gladius sampled from the same Mediterranean waters (Fig. 3A). Anisakis pegreffii was also found in co-infection with A. physeteris genetically identified in samples from the central and the southern Tyrrhenian Sea. However, the significantly different infection pattern with A. pegreffii and the absence of A. physeteris in the swordfish sampled from the Ionian Sea compared with the central and southern Tyrrhenian Sea (Fig. 3A) suggests a further investigation of swordfish from other parts of the eastern Mediterranean Sea to search for a possible accordance as suggested by genetic/molecular data (Viñas et al. Reference Viñas, Pérez-Serra, Vidal, Alvarado-Bremer and Pla2010). On the other hand, a high proportion of A. physeteris was found and showed higher infection levels in swordfish fished in east-tropical and tropical-equatorial Atlantic waters than in those observed in Mediterranean samples (Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011; Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ). Indeed, in the Atlantic swordfish, A. physeteris, A. brevispiculata, A. paggiae and Anisakis sp. 2 are the dominant taxa. These findings are in accordance with the geographic distribution reported so far for those species (Fig. 3A and B). The PCA analysis based on the distribution pattern and infection levels with Anisakis spp. from swordfish showed the Mediterranean samples clustering separately from the Atlantic ones (Fig. 3B). In addition, in the Atlantic waters, a cluster formed by the ‘northern’ sample distinct from a cluster including the ‘southern’ sample was observed (Mattiucci et al. Reference Mattiucci, Garcia, Cipriani, Santos, Nascetti and Cimmaruta2014b ) (Fig. 3B).

The use of Anisakis spp. as biomarkers agrees with the genetic results gathered on the fish host when describing the existence of a Mediterranean Sea stock and two discrete Atlantic stocks, i.e. a ‘northern’ and a ‘southern’ Atlantic stock of the fish species (Fig. 3C). Indeed, a population genetic analysis based on allozyme loci of different swordfish populations from Mediterranean and Atlantic waters clearly distinguished the Mediterranean stock as distinct from the Atlantic one, as summarized by the PCA analysis (Fig. 3C, data from Cimmaruta et al. Reference Cimmaruta, Iaconelli, Nascetti and Bullini1998 and unpublished). The first axis of PCA alone explained 96% of total ordination (another 2% was explained by the second axis) illustrating that for this fish species even conservative markers, such as allozymes, have a high-resolution power.

However, when defining the boundary between the two Atlantic stocks, the results obtained by different genetic/molecular markers on X. gladius are discordant with each other, either confirming the 5°N position, moving it north of 5°N, or suggesting an intergradation zone in the equatorial area. According to the Anisakis spp. used as markers, the distinction between the NW swordfish stock from the remaining Atlantic populations is in agreement with the bulk of genetic data, while the boundary between the NW and the ‘southern’ stock seems to be less certain. This suggests a possible boundary between the northern and southern stocks, north of 5°N, or at least the existence of a mixing area in the tropical-equatorial Atlantic waters (Garcia et al. Reference Garcia, Mattiucci, Damiano, Santos and Nascetti2011; Mattiucci, personal communication) (Fig. 3A).

The North Atlantic swordfish recently benefited from the management actions that rebuilt its stock consistency (Neilson et al. Reference Neilson, Arocha, Cass-Calay, Mejuto, Ortiz, Scott, Smith, Travassos, Tserpes and Andrushchenko2013). On the contrary, the Mediterranean swordfish is still in need of effective management to counteract the loss of its spawning stock biomass and the ever-increasing percentage of juveniles in its landings, signalling a possible stock collapse into the next generation time of 7–10 years (Tserpes et al. Reference Tserpes, Tzanatos and Peristeraki2011).

DISCUSSION AND FUTURE PERSPECTIVES

A high genetic homogeneity of marine fish species has been postulated for years, based on the idea that in the marine environment there are few physical barriers limiting gene flow between fish populations (Palumbi, Reference Palumbi1994). However, in the last decade data have increasingly shown that marine species may be spatially structured even on a very small scale due to a number of factors ranging from local self-recruitment to historical vicariance (Cowen et al. Reference Cowen, Lwiza, Sponaugle, Paris and Olson2000, Reference Cowen, Paris and Srinivasan2006). Also, the marine environment proved to have geographic barriers, habitat discontinuities and transition zones which were able to exert different selective pressures on geographically discrete populations of marine species (Cimmaruta et al. Reference Cimmaruta, Bondanelli and Nascetti2005; Lemaire et al. Reference Lemaire, Versini and Bonhomme2005; Riginos et al. Reference Riginos, Douglas, Jin, Shanahan and Treml2011; Sà-Pinto et al. Reference Sà-Pinto, Branco, Alexandrino, Fontaine and Baird2012; Milano et al. Reference Milano, Babbucci, Cariani, Atanassova, Bekkevold, Carvalho, Espiñeira, Fiorentino, Garofalo, Geffen, Hansen, Helyar, Nielsen, Ogden, Patarnello, Stagioni, Tinti and Bargelloni2014). The historical and oceanographic characteristics of the Mediterranean–Atlantic transition are among the main factors in designating distinctive and geographically separated stocks of several fish species, as shown by an increasing number of investigations carried out at both large and fine geographic scales and providing phylogeographical scenarios for several fish species from these basins. Surveys of the genetic structure of fish populations over a large geographic range, spanning the North-East Atlantic Ocean, the Mediterranean Sea (west and east) and the Black Seas, included many fish species, among them Engraulis encrasicolus (Magoulas et al. Reference Magoulas, Tsimenides and Zouros1996), sparid species (Bargelloni et al. Reference Bargelloni, Alarcon, Alvarez, Penzo, Magoulas, Reis and Patarnello2003, Reference Bargelloni, Alarcon, Alvarez, Penzo, Magoulas, Reis and Patarnello2005), the anglerfish Lophius budegassa and Lophius piscatorius (Charrier et al. Reference Charrier, Chenel, Durand, Girard, Quiniou and Laroche2006), the Atherinidae Atherina spp. (Pujolar et al. Reference Pujolar, Zane and Congiu2012) and the mackerels Scomber scombrus and Scomber japonicus (Zardoya et al. Reference Zardoya, Castilho, Grande, Favre-Krey, Caetano, Marcato, Krey and Patarnello2004). In spite of the high number of studies, very few contributions report a multidisciplinary approach to fish stock assessment in Mediterranean vs Atlantic waters, which combines genetic data, biological (parasites) and morphological surveys of the fish host. To date, the only representative cases are those concerning fish stock assessment of horse mackerel and European hake, performed in the framework of the European projects HOMSIR (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ) and GENHAKE (Nascetti et al. Reference Nascetti, Abaunza, Mattiucci, Magoulas, Cimmaruta and Bullini2000) respectively. Some other authors have used a meta-analysis approach on biological (not including parasites) and genetic datasets acquired for various fish species to provide an integrated picture of Mediterranean–Atlantic population structure and phylogeography (Patarnello et al. Reference Patarnello, Volckaert and Castilho2007), also in the context of global patterns (Galarza et al. Reference Galarza, Carreras-Carbonell, Macpherson, Pascual, Roque, Turner and Rico2009; Riginos et al. Reference Riginos, Douglas, Jin, Shanahan and Treml2011). Their findings suggested an overall poor ability of biological factors as predictors of population/stock genetic structure (Galarza et al. Reference Galarza, Carreras-Carbonell, Macpherson, Pascual, Roque, Turner and Rico2009) or failed to show any particular association between biological traits and genetic variables (Patarnello et al. Reference Patarnello, Volckaert and Castilho2007), but all highlighted the role of biogegraphic-environmental factors in shaping the genetic structure of fish species.

In the present review we show that Anisakis parasites, used as biological tags, do provide a biological dataset able to identify fish stock structure generally in agreement with the findings of molecular approaches on the fish species. The meta-analysis approach performed by generalized Procrustes Rotation analysis on the genetic (allozyme markers) and biological (Anisakis spp. parasites) datasets acquired from the same fish populations of the Mediterranean Sea and Atlantic Ocean, in the framework of a multidisciplinary approach, produced a significant association between parasites as biomarkers, and the genetic variation detected at the allozyme level and agreed in delineating the stock structure of the demersal and pelagic fish studied (Figs 1C and 2C). The results obtained also suggest that the combined datasets produce stronger information on the population structure of a fish species rather than that indicated by a single dataset (Fig. 2C). This is due to the different features of the two kinds of data, where the parasites (biomarkers) provide knowledge of the movements of fish specimens for feeding purposes during their lifespan, while the fish population genetic data (molecular/genetic markers) provide those for gene flow averaged over many generations. Thus, genetic studies on fish host populations can detect changes over an evolutionary timescale. These differences may cause discordance in the results achieved by the two different methodological approaches, as shown for example in the case of T. trachurus stock assessment, where fish genetic data provided less structuring power than the parasites which, on the contrary, accounted for most of the differences found at the meta-population level. It is also appropriate to mention here that in stock identification studies, it may be a mistake to conclude that there is only one stock when differences are not observed, given the fact that failing to reject H0 (=Hypothesis ‘zero’, according to which, one single stock exists when differences are not observed) does not mean that the hypothesis has been proven to hold true (Gelman and Hill, Reference Gelman and Hill2007). On the other hand, parasites as biomarkers are generally more concordant with the biological and life-history traits of the fish host, as demonstrated in the case of T. trachurus (Abaunza et al. Reference Abaunza, Murta, Campbell, Cimmaruta, Comesana, Dahle, Garcia Santamaria, Gordo, Iversen, MacKenzie, Magoulas, Mattiucci, Molloy, Nascetti, Pinto, Quinta, Ramos, Sanjuan, Santos, Stransky and Zimmermann2008a ) and M. merluccius. The occurrence and composition of the parasite fauna in fish populations, as for other fish biological parameters, are very good and suitable biomarkers when considering a fish stock in a temporal (as related to the fish life-history) and spatial (as related to the geographic region) scale, because they reflect the environmental variables in which the fish species spends its life history. Thus, their potential application as biomarkers ranges from providing knowledge on fish migrations across different marine ecosystems to the identification of the marine grounds of origin for a given fish species.

This, in turn, means that the phenotypic markers (including parasites) of a fish species have greater application in studies on the short-term variation caused by the environment and other abiotic and biotic factors of a marine ecosystem (Begg and Waldman, Reference Begg and Waldman1999). Since the genetic characterization of stocks arises over evolutionary timescales, the comparison of incongruences between genetic and parasitological patterns provides information on both the mechanisms and the timescales involved in stock differentiation. This kind of information can be obtained only within the frame of a multidisciplinary approach to fish stock identification, demonstrating the high power of this approach.

Based on their features as biological tags, parasites can be used as good indicators of ecosystems in those regions inhabited by the fish species. Indeed, the parasites and their density in a region could reflect the temporal and spatial population dynamics of the fish species in a given geographical region. However, in turn, the demography and population dynamics of definitive and intermediate host populations could also influence the demography and density of a given host-specific parasite (Mattiucci and Nascetti, Reference Mattiucci and Nascetti2007, Reference Mattiucci and Nascetti2008). Thus, monitoring a fish subpopulation or stock at the regional scale level, including its specific parasite fauna and their abundance levels as biomarkers, could have implications for fisheries management and conservation measures.

In addition, future comparative co-phylogeographic studies between parasites and their hosts may be useful in determining fish host populations over different geographical areas. Indeed, parasite phylogeography could provide information to infer host population structure (Criscione and Blouin, Reference Criscione and Blouin2007). The relative mutation rates of homologous genes between parasite DNA and that of their hosts suggests that some parasites have a faster mutation rate than their hosts (Whiteman and Parker, Reference Whiteman and Parker2005). This could lead to a faster accumulation of genetic variation in the parasite genome allowing us to infer host population structure from the parasite population structure. In other words, the divergence time separating parasite populations could be shorter than that recorded in the host population. Different inferences could be made by the different genetic markers used, and this applies to both the fish hosts and their parasite species. In this contest, the Procrustes Approach to Co-Phylogeny (PACo) has been proposed as a statistical tool to test the congruence for co-phylogenetic aspects of host-parasite associations (Balbuena et al. Reference Balbuena, Miguez-Lozano and Blasco-Costa2013).

Mitochondrial DNA genes, due to the faster accumulation of base pairs substitution and maternal inheritance, are suitable genetic markers to infer population structure in both fish and parasite species. For instance, Criscione et al. (Reference Criscione, Cooper and Blouin2006) indicated that parasite trematode genotypes are more sensitive for indicating the natal stream of an individual of Oncorhynchus mykiss, rather than the host genotypes. Similarly, Baldwin et al. (Reference Baldwin, Rew, Johansson, Banks and Jacobson2011) used population genetic structure based on the cytochrome c oxidase 2 (cox2) mitochondrial DNA gene of larval nematodes of the genus Anisakis to infer the stock structure of Pacific sardine (Sardinops sagax) in the California current. Mitochondrial haplotypes distribution of larval Anisakis spp. obtained by those authors, from both parsimony networks and analyses of molecular variance, revealed a panmictic distribution of these parasites, which infect sardines throughout the California current ecosystem. Based on the results achieved, Baldwin et al. (Reference Baldwin, Rew, Johansson, Banks and Jacobson2011) hypothesized that the panmictic distribution of the larval Anisakis spp. populations may be a result of the presumed migratory pathways of the intermediate host (the Pacific sardine), moving into the northern portion of the California current in summer and returning to the southern portion to overwinter and spawn in spring. However, the wider geographic range of paratenic (large predatory fish), and final (cetaceans) hosts has also been postulated to explain the observed distribution pattern.

On the other hand, other mitochondrial markers, such as the mtDNA cox1 have been suggested as a suitable indicator of nematode population structure of A. simplex (s.s.), potentially providing knowledge about its fish host, the Atlantic herring Clupea harengus (Cross et al. Reference Cross, Collins, Campbell, Watts, Chubb, Cunningham, Hatfield and MacKenzie2007). This study indicated that mtDNA cox1 gene sequences of A. simplex (s.s.) analysed between and within spawning seasons of C. harengus, showed a high degree of temporal stability, highlighting the potential suitability of this molecular marker in this parasite species as a biological tag of herring stocks (Cross et al. Reference Cross, Collins, Campbell, Watts, Chubb, Cunningham, Hatfield and MacKenzie2007). Future studies on phylogeographic analysis, including population genetic structure of parasites as biomarkers and that of their host evolutionary units, inferred from highly polymorphic genes (such as DNA microsatellites, mitochondrial DNA genes, Single Nucleotide Polymorphisms (SNPs) markers) could represent a further tool to be included in an holistic approach to infer fine-scale genetic structuring of a fish species over its range of distribution.

ACKNOWLEDGEMENTS

We dedicate this review to the memory of our dear friend, Dr Alberto Murta. He left us, prematurely, during the preparation of this article. He was a passionate scientist; we have shared with him the challenge of using the multidisciplinary approach in the framework of fish stock identification, particularly during the HOMSIR Project. We dedicate this paper to his skill and expertise in fish population dynamics, to his work performed on the use of parasites as biological tags, to his enthusiasm for science and for life.

FINANCIAL SUPPORT

Parts of the studies reviewed in this paper were funded by the European Community: GENHAKE [FAIR CT97-3494]; HOMSIR [QLK5-CT99-01438], by the Italian Ministry for Aquaculture and Fisheries (MIPAF), and by the Italian Ministry for the Research and University (MIUR).

References

REFERENCES

Abaunza, P., Gordo, L., Karlou-Riga, C., Murta, A., Eltink, A. T. G. W., Garcia Santamaria, M. T., Zimmermann, C., Hammer, C., Lucio, P., Iversen, S. A., Molloy, J. and Gallo, E. (2003). Growth and reproduction of horse mackerel, Trachurus trachurus (Carangidae). Review of Fish Biology and Fisheries 13, 2761.CrossRefGoogle Scholar
Abaunza, P., Murta, A. G., Campbell, N., Cimmaruta, R., Comesana, A. S., Dahle, G., Garcia Santamaria, M. T., Gordo, L. S., Iversen, S. A., MacKenzie, K., Magoulas, A., Mattiucci, S., Molloy, J., Nascetti, G., Pinto, A. L., Quinta, R., Ramos, P., Sanjuan, A., Santos, A. T., Stransky, C. and Zimmermann, C. (2008 a). Stock identity of horse mackerel (Trachurus trachurus) in the Northeast Atlantic and Mediterranean Sea: integrating the results from different stock identification approaches. Fisheries Research 89, 196209.CrossRefGoogle Scholar
Abaunza, P., Murta, A. G., Campbell, N., Cimmaruta, R., Comesana, A. S., Dahle, G., Gallo, E., Garcia Santamaria, M. T., Gordo, L. S., Iversen, S. A., MacKenzie, K., Magoulas, A., Mattiucci, S., Molloy, J., Nascetti, G., Pinto, A. L., Quinta, R., Ramos, P., Ruggi, A., Sanjuan, A., Santos, A. T., Stransky, C. and Zimmerman, C. (2008 b). Considerations on sampling strategies for an holistic approach to stock identification: the example of the HOMSIR project. Fisheries Research 89, 104113.Google Scholar
Abaunza, P., Gordo, L., García Santamaría, M. T., Iversen, S. A., Murta, A. G. and Gallo, E. (2008 c). Life history parameters as an important basis for the initial recognition of stock management units in horse mackerel (Trachurus trachurus). Fisheries Research 89, 167180.Google Scholar
Alegria-Hernandez, V. (1986). Catch per unit effort as an index of sardine, Sardina pilchardus (Walb.) distribution changes along the eastern Adriatic. FAO Fisheries Report 345, 127136.Google Scholar
Alvarado Bremer, J. R., Mejutoc, J., Gómez-Márquezd, J., Boánd, F., Carpinterod, P., Rodríguez, J. M., Viñase, J., Greig, T. W. and Elyg, B. (2005 a). Hierarchical analyses of genetic variation of samples from breeding and feeding grounds confirm the genetic partitioning of northwest Atlantic and South Atlantic populations of swordfish (Xiphias gladius L.). Journal of Experimental Marine Biology and Ecology 327, 167182.Google Scholar
Alvarado Bremer, J. R., Viñasa, J., Mejutoc, J., Elyd, B. and Plae, C. (2005 b). Comparative phylogeography of Atlantic bluefin tuna and swordfish: the combined effects of vicariance, secondary contact, introgression, and population expansion on the regional phylogenies of two highly migratory pelagic fishes. Molecular Phylogenetics and Evolution 36, 169187.Google Scholar
Alvarado Bremer, J. R., Hinton, M. G. and Grieg, T. W. (2006). Evidence of spatial genetic heterogeneity in Pacific swordfish (Xiphias gladius) revealed by the analysis of ldh-A sequences. Bulletin of Marine Science 79, 493503.Google Scholar
Bahri-Sfar, L., Lemaire, C., Hassine, B. and Bonhomme, F. (2000). Fragmentation of sea bass populations in the western and eastern Mediterranean as revealed by microsatellite polymorphism. Proceedings of the Royal Society of London Series B – Biological Sciences 267, 929935.Google Scholar
Balbuena, J. A., Miguez-Lozano, R. and Blasco-Costa, I. (2013). PACo: a novel Procrustes application to cophylogenetic analysis. PLOS ONE 8, e61048.CrossRefGoogle Scholar
Baldwin, R. E., Rew, M. B., Johansson, M. L., Banks, M. A. and Jacobson, K. C. (2011). Population structure of three species of Anisakis nematodes recovered from Pacific sardines (Sardinopsis sagax) distributed throughout the California Current System. Journal of Parasitology 97, 545554.Google Scholar
Bargelloni, L., Alarcon, J. A., Alvarez, M. C., Penzo, E., Magoulas, A., Reis, C. and Patarnello, T. (2003). Discord in the family Sparidae (Teleostei): divergent phylogeographical patterns across the Atlantic–Mediterranean divide. Journal of Evolutionary Biology 16, 11491158.Google Scholar
Bargelloni, L., Alarcon, J. A., Alvarez, M. C., Penzo, E., Magoulas, A., Reis, C. and Patarnello, T. (2005). The Atlantic-Mediterranean transition: discordant genetic patterns in two seabream species, Diplodus puntazzo (Cetti) and Diplodus sargus (L.). Molecular Phylogenetic and Evolution 36, 523535.Google Scholar
Begg, G. and Waldman, J. R. (1999). An holistic approach to fish stock identification. Fisheries Research 43, 3544.Google Scholar
Bembo, D. G., Carvalho, G. R., Snow, M., Cingolani, N. and Pitcher, T. J. (1995). Stock discrimination among European anchovies, Engraulis encrasicolus, by means of PCR-amplified mitochondrial DNA analysis. Fishery Bulletin 94, 3140.Google Scholar
Bembo, D. G., Carvalho, G. R., Cingolani, N., Arneri, E., Giannetti, G. and Pitcher, T. J. (1996 a). Allozymic and morphometric evidence for two stocks of the European anchovy Engraulis encrasicolus in Adriatic waters. Marine Biology 126, 529538.Google Scholar
Bembo, D. G., Carvalho, G. R., Cingolani, N. and Pitcher, T. J. (1996 b). Electrophoretic analysis of stock structure in northern Mediterranean anchovies, Engraulis encrasicolus . ICES Journal of Marine Science 53, 115128.Google Scholar
Berland, B. (1961). Nematodes from some Norwegian marine fishes. Sarsia 2, 150.Google Scholar
Cacaud, P. (2005). Fisheries laws and regulations in the Mediterranean: a comparative study. Studies and Reviews/General Fisheries Commission for the Mediterranean 7, 40.Google Scholar
Caddy, J. F. (1995). Cephalopod and demersal finfish stocks: some statistical trends and biological interactions. In Squid 94 Venice, Proceedings of the Third International Cephalopod Trade Conference, Agra Europe, London, UK.Google Scholar
Charrier, G., Chenel, T., Durand, J. D., Girard, M., Quiniou, L. and Laroche, J. (2006). Discrepancies in phylogeographical patterns of two European anglerfishes (Lophius budegassa and Lophius piscatorius). Molecular Phylogenetics and Evolution 38, 742754.Google Scholar
Chow, S. and Takeyama, H. (2000). Nuclear and mitochondrial DNA analyses reveal four genetically separated breeding units of the swordfish (Xiphias gladius). Journal of Fish Biology 56, 10871098.Google Scholar
Chow, S., Okamoto, H., Uozumi, Y., Takeuchi, Y. and Takeyama, H. (1997). Genetic stock structure of the swordfish (Xiphias gladius) inferred by PCR–RFLP analysis of the mitochondrial DNA control region. Marine Biology 127, 359367.Google Scholar
Chow, S., Clarke, S., Nakadate, M. and Okazaki, M. (2007). Boundary between the north and south Atlantic populations of the swordfish (Xiphias gladius) inferred by a single nucleotide polymorphism at calmodulin gene intron. Marine Biology 152, 8793.Google Scholar
Cimmaruta, R., Iaconelli, M., Nascetti, G. and Bullini, L. (1998). Diversità genetica in popolazioni di grandi pelagici del Mediterraneo. Biologia Marina Mediterranea 5, 300310.Google Scholar
Cimmaruta, R., Bondanelli, P. and Nascetti, G. (2005). Genetic structure and environmental heterogeneity in the European hake (Merluccius merluccius). Molecular Ecology 14, 25772591.Google Scholar
Cimmaruta, R., Bondanelli, P., Ruggi, A. and Nascetti, G. (2008). Genetic structure and temporal stability in the horse mackerel (Trachurus trachurus). Fisheries Research 89, 114121.CrossRefGoogle Scholar
Comesaña, A. S., Martinez-Areal, M. T. and Sanjuan, A. (2008). Genetic variation in the mitochondrial DNA control region among horse mackerel (Trachurus trachurus) from the Atlantic and Mediterranean areas. Fisheries Research 89, 122131.Google Scholar
Cowen, R. K. and Sponaugle, S. (2009). Larval dispersal and marine population connectivity. Annual Review of Marine Science 1, 443466.Google Scholar
Cowen, R. K., Lwiza, K. M. M., Sponaugle, S., Paris, C. B. and Olson, D. B. (2000). Connectivity of marine populations: open or closed? Science 287, 857859.Google Scholar
Cowen, R. K., Paris, C. B. and Srinivasan, A. (2006). Scaling of connectivity in marine populations. Science 311, 522527.Google Scholar
Criscione, C. D. and Blouin, M. S. (2007). Parasite phylogeographical congruence with salmon host evolutionary units: implications for salmon conservation. Molecular Ecology 16, 9931005.Google Scholar
Criscione, C. D., Cooper, B. and Blouin, M. S. (2006). Parasite genotypes identify source populations of migratory fish more accurately than fish genotypes. Ecology 87, 823828.Google Scholar
Crosetti, D., Nelson, W. S. and Avise, J. C. (1994). Pronounced genetic structure of mitochondrial DNA among populations of the circumglobally distributed grey mullet (Mugil cephalus). Journal of Fish Biology 44, 4758.Google Scholar
Cross, M. A., Collins, C., Campbell, N., Watts, P. C., Chubb, J. C., Cunningham, C. O., Hatfield, E. M. C. and MacKenzie, K. (2007). Levels of intra-host and temporal sequence variation in a large CO1 subunit from Anisakis simplex sensu stricto (Rudolphi, 1809) (Nematoda: Anisakidae): implications for fisheries managements. Molecular Biology 151, 695702.Google Scholar
Ehrhardt, N. M. (1992). Age and growth of swordfish, Xiphias gladius, in the Northwestern Atlantic. Bulletin of Marine Science 50, 292301.Google Scholar
Fernandes, P. G. and Cook, R. M. (2013). Reversal of fish stock decline in the Northeast Atlantic. Current Biology 23, 14321437.Google Scholar
Ferrer-Castelló, E., Raga, J. A. and Aznar, F. J. (2007). Parasites as fish population tags and pseudoreplication problems: the case of striped red mullet Mullus surmuletus in the Spanish Mediterranean. Journal of Helminthology 81, 169178.Google Scholar
Food and Agriculture Organization. (1999). Aquaculture Production Statistics 1988–1997. Fisheries Circular No. 815 Rev. 11, FAO, Rome, Italy.Google Scholar
Froese, R. and Pauly, D. (eds.) (2014). FishBase. World Wide Web Electronic Publication. Available at: www.fishbase.org, version (01/2014).Google Scholar
Gačić, M., Schroeder, K., Civitarese, G., Vetrano, A. and Eusebi Borzelli, G. L. (2012). On the relationship among the Adriatic-Ionian bimodal Oscillating System (BiOS), the Eastern Mediterranean salinity variations and Western Mediterranean thermohaline cell. Ocean Science Discussion 9, 25612580.Google Scholar
Galarza, J. A., Carreras-Carbonell, J., Macpherson, E., Pascual, M., Roque, S., Turner, G. F. and Rico, C. (2009). The influence of oceanographic fronts and early-life-history traits on connectivity among littoral fish species. Proceedings of the National Academy of Sciences USA 106, 14731478.Google Scholar
Garcia, A., Mattiucci, S., Damiano, S., Santos, M. N. and Nascetti, G. (2011). Metazoan parasites of swordfish, Xiphias gladius (Pisces: Xiphiidae) from the Atlantic Ocean: implications for host stock identification. ICES Journal of Marine Science 68, 175182.Google Scholar
Gelman, A. and Hill, J. (2007). Data Analysis Using Regression and Multilevel/Hierarchical Models. Cambridge University Press, Cambridge, UK.Google Scholar
GFCM (2012). General Fisheries Commission for the Mediterranean (GFCM). Scientific Advisory Committee on Stock Assessment.Google Scholar
Gordo, L. S., Costa, A., Abaunza, P., Lucio, P., Eltink, A. T. G. W. and Figueiredo, I. (2008). Determinate versus indeterminate fecundity in horse mackerel. Fisheries Research 89, 181185.Google Scholar
ICES (2013). Advice 2013. North Sea. Horse mackerel (Trachurus trachurus) in Divisions IIIa, IVb,c, and VIId (North Sea stock). http://www.ices.dk/sites/pub/Publication%20Reports/Advice/2013/2013/hom-nsea.pdf.Google Scholar
Kasapidis, P. and Magoulas, A. (2008). Development and application of microsatellite markers to address the population structure of the horse mackerel Trachurus trachurus . Fisheries Research 89, 132134.Google Scholar
Kasapidis, P., Mejuto, J., Tserpes, G., Antoniou, A., Garcia-Cortes, B., Peristeraki, P., Oikioonomaki, K., Kotoulas, G. and Magoulas, A. (2007). Genetic structure of swordfish (Xiphias gladius) stocks in the Atlantic using microsatellite DNA analysis. Collective Volume of Scientific Papers ICCAT 61, 8998.Google Scholar
Klimpel, S., Kleinertz, S. and Palm, H. W. (2008). Distribution of parasites from red mullets (Mullus surmuletus L., Mullidae) in the North Sea and Mediterranean Sea. Bulletin of Fish Biology 10, 2538.Google Scholar
Knutsen, H., Jorde, P. E., André, C. and Stenseth, N. C. (2003). Fine-scaled geographical population structuring in a highly mobile marine species: the Atlantic cod. Molecular Ecology 12, 385394.Google Scholar
Kotoulas, G., Magoulas, A., Tsimenides, N. and Zouros, E. (1995). Marked mitochondrial DNA differences between Mediterranean and Atlantic populations of the swordfish, Xiphias gladius . Molecular Ecology 4, 473481.Google Scholar
Lemaire, C., Versini, J. J. and Bonhomme, F. (2005). Maintenance of genetic differentiation across a transition zone in the sea: discordance between nuclear and cytoplasmic markers. Journal of Evolutionary Biology 18, 7080.Google Scholar
Lleonart, J. and Maynou, F. (2003). Fish stock assessments in the Mediterranean: state of art. Scientia Marina 67, 3749.Google Scholar
Lundy, C., Moran, P., Rico, C., Milner, R. S. and Hewitt, G. M. (1999). Macrogeographical population differentiation in oceanic environments: a case study of European hake (Merluccius merluccius) a commercially important fish. Molecular Ecology 8, 18891898.Google Scholar
Lundy, C. J., Rico, C. and Hewitt, G. M. (2000). Temporal and spatial genetic variation in spawning grounds of European hake (Merluccius merluccius) in the Bay of Biscay. Molecular Ecology 9, 20672079.Google Scholar
MacKenzie, K., Campbell, N., Mattiucci, S., Ramos, P., Pinto, A. L. and Abaunza, P. (2008). Parasites as biological tags for stock identification of Atlantic horse mackerel Trachurus trachurus L. Fisheries Research 89, 136145.Google Scholar
Magoulas, A., Kotoulas, G., de la Serna, J. M., De Metrio, G., Tsimenides, N. and Zouros, E. (1992). Genetic structure of swordfish (Xiphias gladius) populations of the Mediterranean and the eastern side of the Atlantic: analysis by mitochondrial DNA markers. Collective Volume of Scientific Papers ICCAT XL, 126136.Google Scholar
Magoulas, A., Tsimenides, N. and Zouros, E. (1996). Mitochondrial DNA phylogeny and the reconstruction of the population history of a species: the case of the European anchovy (Engraulis encrasicolus). Molecular Biology and Evolution 13, 178190.Google Scholar
Mattiangeli, V., Ryan, A. W., Galvin, P., Mork, J. and Cross, T. F. (2003). Eastern and western poor cod (Trisopterus minutus capelanus) populations in the Mediterranean Sea: evidence from allozyme and minisatellite loci. Marine Ecology 24, 247258.Google Scholar
Mattiucci, S. and Nascetti, G. (2007). Genetic diversity and infection levels in anisakid nematodes parasitic in fish and marine mammals from Boreal and Austral hemispheres. Veterinary Parasitology 148, 4357.Google Scholar
Mattiucci, S. and Nascetti, G. (2008). Advances and trends in the molecular systematics of Anisakid nematodes, with implications for their evolutionary ecology and host-parasite co-evolutionary processes. Advances in Parasitology 66, 47148.CrossRefGoogle ScholarPubMed
Mattiucci, S., Nascetti, G., Cianchi, R., Paggi, L., Arduino, P., Margolis, L., Brattey, J., Webb, S. C., D'Amelio, S., Orecchia, P. and Bullini, L. (1997). Genetic and ecological data on the Anisakis simplex complex, with evidence for a new species (Nematoda: Ascaridoidea, Anisakidae). Journal of Parasitology 83, 401416.Google Scholar
Mattiucci, S., Paggi, L., Nascetti, G., Portes Santos, C., Costa, G., Di Benedetto, A. P., Ramos, R., Argyrou, M., Cianchi, R. and Bullini, L. (2002). Genetic markers in the study of Anisakis typica (Diesing, 1860): larval identification and genetic relationships with other species of Anisakis Dujardin, 1845 (Nematoda: Anisakidae). Systematic Parasitology 51, 159170.Google Scholar
Mattiucci, S., Abaunza, P., Ramadori, L. and Nascetti, G. (2004). Genetic identification of Anisakis larvae in European hake from Atlantic and Mediterranean waters for stock recognition. Journal of Fish Biology 65, 495510.Google Scholar
Mattiucci, S., Abaunza, P., Damiano, S., Garcia, A., Santos, M. N. and Nascetti, G. (2007). Distribution of Anisakis larvae, identified by genetic markers, and their use for stock characterization of demersal and pelagic fish from European waters: an update. Journal of Helminthology 81, 117127.Google Scholar
Mattiucci, S., Farina, V., Campbell, N., Mackenzie, K., Ramos, P., Pinto, A. L., Abaunza, P. and Nascetti, G. (2008). Anisakis spp. larvae (Nematoda: Anisakidae) from Atlantic horse mackerel: their genetic identification and use as biological tags for host stock identification. Fisheries Research 89, 146151.Google Scholar
Mattiucci, S., Cipriani, P., Webb, S. C., Paoletti, M., Marcer, F., Bellisario, B., Gibson, D. I. and Nascetti, G. (2014 a). Genetic and morphological approaches distinguish the three sibling species of the Anisakis simplex species complex, with a species designation as Anisakis berlandi n. sp. for A. simplex sp. C (Nematoda: Anisakidae). Journal of Parasitology 100, 199214.Google Scholar
Mattiucci, S., Garcia, A., Cipriani, P., Santos, M. N., Nascetti, G. and Cimmaruta, R. (2014 b). Metazoan parasite infection in the swordfish, Xiphias gladius, from the Mediterranean Sea and comparison with Atlantic populations: implications for its stock characterization. Parasite 21, 35. doi: 10.1051/parasite/2014036.Google Scholar
Milano, I., Babbucci, M., Cariani, A., Atanassova, M., Bekkevold, D., Carvalho, G. R., Espiñeira, M., Fiorentino, F., Garofalo, G., Geffen, A. J., Hansen, J. H., Helyar, S. J., Nielsen, E. E., Ogden, R., Patarnello, T., Stagioni, M., FishPopTrace Consortium, Tinti, F. and Bargelloni, L. (2014). Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius). Molecular Ecology 23, 118135.Google Scholar
Mladineo, I., Simat, V., Miletic, J., Beck, R. and Poliack, V. (2012). Molecular identification and population dynamics of Anisakis pegreffii (Nematoda: Anisakidae Dujardin, 1845) isolated from the European anchovy (Engraulis encrasicolus L.) in the Adriatic Sea. International Journal of Food Microbiology 157, 224227.Google Scholar
Murta, A. G., Pinto, A. L. and Abaunza, P. (2008). Stock identification of horse mackerel (Trachurus trachurus) through the analysis of body shape. Fisheries Research 89, 152158.Google Scholar
Naciri, M., Lemaire, C., Borsa, P. and Bonhomme, F. (1999). Genetic study of the Atlantic/Mediterranean transition in sea bass (Dicentrarchus labrax). Journal of Heredity 90, 591596.Google Scholar
Nascetti, G., Paggi, L., Orecchia, P., Smith, J. W., Mattiucci, S. and Bullini, L. (1986). Electrophoretic studies on the Anisakis simplex complex (Ascaridida: Anisakidae) from the Mediterranean and North-East Atlantic. International Journal for Parasitology 16, 633640.Google Scholar
Nascetti, G., Abaunza, P., Mattiucci, S., Magoulas, A., Cimmaruta, R. and Bullini, L. (2000). Genetic Diversity and Parasitic Infections of European Hake (Merluccius merluccius) as Descriptors of Fishing Impact on Mediterranean and Atlantic Stocks. Final Report of the European Project GENHAKE (FAIR CT 97–3494).Google Scholar
Neilson, J., Arocha, F., Cass-Calay, S., Mejuto, J., Ortiz, M., Scott, G., Smith, C., Travassos, P., Tserpes, G. and Andrushchenko, I. (2013). The recovery of Atlantic swordfish: the comparative roles of the regional fisheries management organization and species biology. Reviews in Fisheries Science 21, 5997.Google Scholar
Notarbartolo di Sciara, G., Agardy, T., Hyrenbach, D., Scovazzi, T. and Van Klaveren, P. (2008). The Pelagos Sanctuary for Mediterranean marine mammals. Aquatic Conservation 18, 367391.Google Scholar
Palomera, I., Olivar, M. P., Salat, J., Sabates, A., Coll, A., Garcıa, A. and Morales-Nin, B. (2007). Small pelagic fish in the NW Mediterranean Sea: an ecological review. Progress in Oceanography 74, 377396.Google Scholar
Palumbi, S. R. (1994). Genetic divergence, reproductive isolation and marine speciation. Annual Review in Ecology and Systematics 25, 547572.Google Scholar
Pannacciulli, F. G., Bishop, J. D. D. and Hawkins, S. J. (1997). Genetic structure of populations of two species of Chthamalus (Crustacea: Cirripedia) in the North-East Atlantic and Mediterranean. Marine Biology 128, 7382.Google Scholar
Patarnello, T., Volckaert, F. M. J. and Castilho, R. (2007). Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break? Molecular Ecology 16, 44264444.Google Scholar
Peres-Neto, P. R. and Jackson, D. A. (2001). How well do multivariate data sets match? The advantages of a Procrustean superimposition approach over the Mantel test. Oecologia 129, 169178.Google Scholar
Pla, C., Tudela, S. and Garcia-Marin, J. L. (1996). Diversity and population structure of the European anchovy in the Mediterranean Sea. In The European Anchovy and its Environment (ed. Palomera, I. and Rubies, P.). Sciencia Marina 60, 285286.Google Scholar
Power, A. M., Balbuena, J. A. and Raga, J. A. (2005). Parasite infracommunities as predictors of harvest location of bogue (Boops boops L.): a pilot study using statistical classifiers. Fisheries Research 72, 229239.Google Scholar
Pujolar, J. M., Zane, L. and Congiu, L. (2012). Phylogenetic relationships and demographic histories of the Atherinidae in the Eastern Atlantic and Mediterranean Sea re-examined by Bayesian inference. Molecular Phylogenetic and Evolution 63, 857865.Google Scholar
Quesada, H., Beynon, C. M. and Skibinski, D. O. F. (1995). A mitochondrial DNA discontinuity in the Mussel Mytilus galloprovincialis Lmk: Pleistocene Vicariance Biogeography and secondary Intergradation. Molecular Biology and Evolution 12, 521524.Google Scholar
Quignard, J. P. (1978). La Méditerranée: creuset ichtyologique. Bollettino di Zoologia 45, 2336.Google Scholar
R Development Core Team (2012). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.Google Scholar
Rello, F. J., Adroher, F. J., Benitez, R. and Valero, A. (2009). The fishing area as a possible indicator of the infection by anisakids in anchovies (Engraulis encrasicolus) from southwestern Europe. International Journal for food Microbiology 129, 277281.Google Scholar
Riginos, C., Douglas, K. E., Jin, Y., Shanahan, D. F. and Treml, E. A. (2011). Effects of geography and life history traits on genetic differentiation in benthic marine fishes. Ecography 34, 566575.Google Scholar
Roldan, M. I., Garcia-Marin, J. L., Utter, F. and Pla, C. (1998). Population genetic structure of European hake, Merluccius merluccius . Heredity 81, 327334.Google Scholar
Ruggeri, P., Splendiani, A., Bonanomi, S., Arneri, E., Cingolani, N., Santojanni, A., Colella, S., Donato, F., Giovannotti, M., Caputo, V. and Barucchi, N. (2013). Searching for a stock structure in Sardina pilchardus from the Adriatic and Ionian seas using a microsatellite DNA-based approach. Scientia Marina 77, 565574.Google Scholar
Sà-Pinto, A., Branco, M. S., Alexandrino, P. B., Fontaine, M. C. and Baird, S. J. E. (2012). Barriers to gene flow in the marine environment: insights from two common intertidal limpet species of the Atlantic and Mediterranean. PLOS ONE 7, e50330.Google Scholar
Serracca, L., Cencetti, E., Battistini, R., Rossini, I., Prearo, M., Pavoletti, E., Fioravanti, M. L., Righetti, M., Di Donfrancesco, B. and Ercolini, C. (2013). Survey on the presence of Anisakis and Hysterothylacium larvae from fishes and squids caught in Ligurian Sea. Veterinary Parasitology 196, 547551.Google Scholar
Smith, B. L. and Alvarado-Bremer, J. R. (2010). Inferring population admixture with multiple nuclear genetic markers and Bayesian genetic clustering in Atlantic swordfish (Xiphias gladius). Collective Volume of Scientific Papers ICCAT 65, 185190.Google Scholar
Stefanni, S. and Thorley, J. L. (2003). Mitochondrial DNA phylogeography reveals the existence of an evolutionary significant unit of the sand goby Pomatoschistus minutus in the Adriatic (Eastern Mediterranean). Molecular Phylogenetics and Evolution 28, 601609.Google Scholar
Stransky, C., Murta, A. G., Schlickeisen, J. and Zimmermann, C. (2008). Otolith shape analysis as a tool for stock separation of horse mackerel (Trachurus trachurus) in the Northeast Atlantic and Mediterranean. Fisheries Research 89, 159168.Google Scholar
Thomas, F., Verneau, O., de Meeûs, T. and Renaud, F. (1996). Parasites as to host evolutionary prints: Insights into host evolution from parasitological data. International Journal for Parasitology 26, 677686.Google Scholar
Tinti, S., Di Nunno, C., Guarniero, I., Talenti, M., Tommasini, S., Fabbri, E. and Piccinetti, C. (2002). Mitochondrial DNA sequence variation suggests the lack of genetic heterogeneity in the Adriatic and Ionian stocks of Sardina pilchardus . Marine Biotechnology 4, 163172.Google Scholar
Tintore, J., La Violette, P. E., Blade, I. and Cruzado, G. (1988). A study of an intense density front in the eastern Alboran sea: the Almería-Oran front. Journal of Physical Oceanography 18, 13841397.Google Scholar
Tserpes, G. and Tsimenides, N. (1995). Determination of age and growth of swordfish, Xiphias gladius L, 1758, in the Eastern Mediterranean using anal-fin spines. Fishery Bulletin 93, 594602.Google Scholar
Tserpes, G., Tzanatos, E. and Peristeraki, P. (2011). Use of risk analysis for the evaluation of different management strategies for the Mediterranean swordfish. Collective Volume of Scientific Papers ICCAT 66, 15061514.Google Scholar
Tudela, S., Garcia-Marin, J. L. and Pla, C. (1999). Genetic structure of the European anchovy, Engraulis encrasicolus, in the North-West Mediterranean. Journal of Experimental Marine Biology and Ecology 234, 95109.Google Scholar
Valsecchi, E., Pasolini, P., Bertozzi, M., Garoia, F., Ungaro, N., Vacchi, M., Sabelli, B. and Tinti, F. (2005). Rapid Miocene–Pliocene dispersal and evolution of Mediterranean rajid fauna as inferred by mitochondrial gene variation. Journal of Evolutionary Biology 18, 436446.Google Scholar
Viñas, J., Pérez-Serra, A., Vidal, O., Alvarado-Bremer, J. R. and Pla, C. (2010). Genetic differentiation between eastern and western Mediterranean swordfish revealed by phylogeographic analysis of the mitochondrial DNA control region. ICES Journal of Marine Science 67, 12221229.Google Scholar
Waldman, J. R. (1999). The importance of comparative studies in stock analysis. Fisheries Research 43, 237246.Google Scholar
Waldman, J. R. (2005). Definition of stocks: an evolving concept. In Stock Identification Methods – Application in Fisheries Science (ed. Cadrin, S. X., Friedland, K. D. and Waldman, J. R.) pp. 716. Elsevier Academic Press, Burlington, MA, USA.Google Scholar
Waples, R. S. (1998). Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. Journal of Heredity 89, 438450.Google Scholar
Ward, R. D., Woodwark, M. and Skibinski, D. O. F. (1994). A comparison of genetic diversity levels in marine, freshwater and anadromous fish. Journal of Fish Biology 44, 213232.Google Scholar
Whiteman, N. K. and Parker, P. G. (2005). Using parasites to infer host population history: a new rationale for parasite conservation. Animal Conservation 8, 175181.Google Scholar
Zardoya, R., Castilho, R., Grande, C., Favre-Krey, L., Caetano, S., Marcato, S., Krey, G. and Patarnello, T. (2004). Differential population structuring of two closely related fish species, the mackerel (Scomber scombrus) and the chub mackerel (Scomber japonicus), in the Mediterranean Sea. Molecular Ecology 13, 17851798.Google Scholar
Figure 0

Fig. 1. The case study of stock identification of a demersal fish species, the European hake Merluccius merluccius, from the Mediterranean Sea and Atlantic Ocean, as inferred from biological data, biomarkers (Anisakis spp. parasites) and the population genetics data (allozymes) of the fish populations (data from Nascetti et al.2000; Mattiucci et al.2004; Cimmaruta et al.2005). (A) Relative proportions of the Anisakis spp. larvae genetically identified (Nparasites = 1950) in M. merluccius from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Mattiucci et al.2004), plotted with the stock units as identified by the biological and morphometric analyses (data from Nascetti et al.2000). A different colour identifies a putative stock unit as defined by morphometric datasets of the fish host. Dotted lines identify the boundaries between hake stock units, as identified by genetic/molecular markers; (B) Principal Component Analysis (PCA) inferred from the larval Anisakis spp. bio-markers, identified genetically in the European hake specimens (M. merluccius sampled in the Mediterranean (white circles) and Atlantic (black circles) basins). The first axis of PCA explains about 60% of total ordination, and PCA 2 about 30%; (C) A generalized Procrustes Rotation (PR) performed with R (R Development Core Team, 2012) was used to compare different ordinations of the European hakes M. merluccius (Nfish = 1306), based on combined molecular and parasite datasets. The latter was computed by considering the number (Nparasites = 1950) of larval parasites belonging to A. pegreffii, A. simplex (s. s.), A. physeteris, A. typica, A. ziphidarum and A. brevispiculata. A Procrustes test (also known as analysis of congruence) with n = 1000 bootstrap replicates was used to estimate the significance of the Procrustes statistics, assessing similarities between different ordinations, and the Procrustes correlation r was derived from the symmetric Procrustes residual (Peres-Neto and Jackson, 2001). Analysis showed a significant correlation (Procrustes r = 0·59, P<0·01) between the Anisakis spp. burden and molecular traits, showing the clustering of sampled specimens in two well-defined clusters corresponding to the Mediterranean (white circles) and Atlantic stocks (black circles). Code numbers and coloured ellipses are related to the sampling localities, as appearing in Fig. 1A.

Figure 1

Fig. 2. The case study of stock assessment of a small pelagic fish species, the horse mackerel Trachurus trachurus, from the Mediterranean Sea and Atlantic Ocean waters, as inferred from biological data, biomarkers (parasites), and the population genetics dataset (allozymes) of the fish populations (data from Abaunza et al.2008a; Cimmaruta et al.2008; Mattiucci et al.2008; Murta et al.2008). (A) Relative proportions of the Anisakis spp. identified genetically (Nparasites = 2200) in T. trachurus from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Mattiucci et al.2008), plotted with the stock units identified by the biological and morphometric analyses (data from Murta et al.2008). Different colours identify each putative stock unit as defined by morphometric datasets of the fish host; (B) Principal Component Analysis (PCA) inferred from larval Anisakis spp. as biomarkers, identified genetically in the horse-mackerel specimens (T. trachurus) sampled in the Mediterranean (white circles) and Atlantic (black circles). The first axis of PCA explains about 92% of total ordination, whilst PCA 2 about 7%; (C) A generalized Procrustes Rotation (PR), performed with R (R Development Core Team, 2012) used to compare different ordinations of the Atlantic horse mackerels (T. trachurus) analysed (Nfish = 2241), based on combined molecular and Anisakis spp. burden distance matrices. The latter was computed by considering the number of individual parasites (Nparasites = 2200) (fourth-root transformed to improve normality and remove the mean/variance relationship) belonging to A. pegreffii, A. simplex (s.s.), A. physeteris and A. typica. A Procrustes test (also known as analysis of congruence) with n = 1000 bootstrap replicates was used to estimate the significance of the Procrustes statistics, assessing similarities between different ordinations, and the Procrustes correlation r was derived from the symmetric Procrustes residual (Peres-Neto and Jackson, 2001). Analysis showed significant correlation (Procrustes r = 0·39, P = 0·05) between the parasitic burden and molecular traits, as shown by the clustering of sampled specimens corresponding to the Mediterranean (white circles) and Atlantic stocks (black circles). Molecular distance between specimens accounted for about 37% of total ordination, whilst the distance based on the parasitic burden with Anisakis for about 95%. Code numbers and coloured ellipses are related to the sampling localities, as appearing in Fig. 2A.

Figure 2

Fig. 3. The case study of stock assessment of a large pelagic fish species, the swordfish Xiphias gladius from the Mediterranean Sea and Atlantic Ocean, as inferred from biological data, biomarkers (parasites) and the genetic datasets (allozymes) of the same fish populations (data from Garcia et al.2011; Mattiucci et al.2014b). (A) Relative proportions of the Anisakis spp. identified genetically (Nparasites = 405) in X. gladius from the sampling localities from the Mediterranean Sea and the Atlantic Ocean (data from Garcia et al.2011; Mattiucci et al.2014b), plotted with the stock units as suggested by the genetic/molecular results (Kotoulas et al.1995; Alvarado-Bremer et al.2005a, b; Chow et al.2007). Dotted lines identify the boundaries between swordfish stock units; (B) Principal Component Analysis (PCA) inferred from the parasitic burden of the larval Anisakis spp. identified genetically and used as biomarkers, in the swordfish X. gladius specimens sampled in the Mediterranean (white circles) and Atlantic (black circles) basins. PCA analysis was computed by considering the number of individual parasites (fourth-root transformed to improve normality and remove the mean/variance relationship) belonging to A. pegreffii, A. simplex (s.s.), A. physeteris and A. brevispiculata. The first axis of PCA explains about 54% of total ordination, and PCA 2 about 22%; (C) Principal Component Analysis (PCA) inferred from Nei's distance matrix based on allozyme data (source: Cimmaruta et al.1998) performed on swordfish specimens analysed genetically (Nfish = 323). The first and second axis accounted for 96 and 2% of total ordination, respectively.