- Open Access
Influence of confluent marine currents in an ecotonal region of the South-West Atlantic on the distribution of larval anisakids (Nematoda: Anisakidae)
Parasites & Vectors volume 11, Article number: 583 (2018)
In the marine environment, transitional zones between major water masses harbour high biodiversity, mostly due to their productivity and by containing representatives of species characteristic of adjacent communities. With the aim of assessing the value of larval Anisakis as zoogeographical indicators in a transitional zone between subtropical and sub-Antarctic marine currents, larvae obtained from Zenopsis conchifer were genetically identified. Larvae from Pagrus pagrus and Merluccius hubbsi from two adjacent zoogeographical provinces were also sequenced.
Four species were genetically identified in the whole sample, including Anisakis typica, A. pegreffii, A. berlandi and a probably new species related to A. paggiae. Anisakis typica and A. pegreffii were identified as indicators of tropical/subtropical and sub-Antarctic waters, respectively, and their presence evidenced the transitional conditions of the region. Multivariate analyses on prevalence and mean abundance of Anisakis spp. of 18 samples represented by 9 fish species caught south of 35°S determined that host trophic level and locality of capture were the main drivers of the distribution of parasites across zoogeographical units in the South-West Atlantic.
Most samples followed a clear zoogeographical pattern, but the sample of Z. conchifer, composed mostly of A. typica, was an exception. This finding suggests that population parameters of A. typica and A. pegreffii could differ enough to be considered as a surrogates of the identity of larvae parasitizing a given host population and, therefore, a step forward the validation of the use of larval Anisakis as biological indicators for studies on host zoogeography.
Members of the genus Anisakis are known worldwide because of their implication in human health as the causative agents of anisakiosis, resulting from the ingestion of infective third-stage larvae in raw or undercooked marine fish products [1,2,3,4] and considered as one of the most significant emerging food-borne zoonoses . Nevertheless, the biological relevance of anisakids in general goes beyond their epidemiological transcendence. Indeed, larval anisakids have been identified among the most suitable biological tags for stock discrimination because they have a lifespan or remain in an identifiable form in the host long enough to cover the timescale of such investigations [6,7,8,9]. However, a limitation to their effectiveness as markers is imposed by difficulties in their identification, since third-stage larvae of several species cannot be identified to species based on traditional morphological analyses, except to the level of the morphotypes of Berland , Anisakis Type I and Type II . Some of these cryptic species often occur in sympatry and syntopy in fish hosts and the lack of taxonomic resolution can affect comparative studies. This limitation in specific reconnaissance for larval stages of Anisakis can, however, be solved with the application of molecular tools, which have recently proved to be of value when parasites are used as biological tags, especially for studies carried out at large geographical scales [12,13,14,15].
The success of larval Anisakis as tags to discriminate host populations at large spatial scales relies on the fact that members of this genus display species-specific distribution patterns within different climate zones and oceans which, in turn, are congruent with those of their respective final hosts . For this reason, the species composition of these parasites in fish can reveal the transitional nature of ecotonal zones between zoogeographical marine regions or interface areas between masses of water, such that observed for hake and blue whiting between the cold Atlantic and the warm Mediterranean waters [12, 17].
In the Argentine Sea, larvae of Anisakis are commonly reported in fish hosts (see  and references therein) as Anisakis sp. or A. simplex (s.l.). The only published reports of genetically identified species of Anisakis are that of Anisakis pegreffii in Merluccius hubbsi  and the skates Sympterygia bonapartii and Zearaja chilensis , and that of a single specimen of Anisakis berlandi in S. bonapartii , highlighting the considerable uncertainty existing in the species composition of this genus in this region. Adult Anisakis have been also reported in several species of cetaceans in the Argentine Sea , all of them based on morphological identifications, and most reported as Anisakis sp. or A. simplex (s.l.), although Berón-Vera et al.  also reported A. physeteris. A similar situation occurs in Brazilian waters, where some species have been morphologically identified in several species of marine cetaceans . However, recent papers, based on genetic identification of larvae and adults, have recorded a higher diversity, A. typica being the most abundant and widely reported species of the genus, occurring in both cetaceans and fish hosts along Brazilian coasts [22, 23].
A promissory couple of species to evaluate the relative influence of confluent marine currents in the South-West Atlantic is represented by A. typica and A. pegreffii. According to data based on genetic identification, the former occurs in warmer temperate and tropical waters between 35°N and 30°S, whereas the distribution range of the later in the Southern Hemisphere extends in temperate to colder regions from 30°S to 60°S [4, 24]. However, in the South-West Atlantic, the border between the distributions of both species could be displaced to higher latitudes due to the influence of the Brazil Current, which flows southwards carrying subtropical waters to collide with the northward flowing Malvinas Current, composed of sub-Antarctic waters, on the continental slope around 38°S in the Argentine-Uruguayan Common Fishing Zone . In the South-West Atlantic, cetacean species distribute differentially along a latitudinal-temperature gradient [26,27,28], and a contribution of different Anisakis species typical for warmer and colder regions should be expected in the confluence region. In a recent paper, Lanfranchi et al.  evaluated the utility of parasites as indicators of marine ecotones by analyzing data on the assemblages of long-lived larval parasites of Zenopsis conchifer inhabiting deep waters in the region of convergence between the Brazil and Malvinas currents, the southernmost limit of its distribution in the South American Atlantic. The ecology of Z. conchifer is little known; however, there is no evidence of migratory movements in the South-West Atlantic, except a shift towards deeper waters as fish grow . Indeed, this fish is considered as a poor swimmer with restricted mobility [31, 32] and consequently constitutes a suitable model to evaluate the presence of infective stages of anisakids in their habitat, by acting as a passive sampler of the available larvae in their prey.
Lanfranchi et al.  included data on other host species recognized as harbouring parasite assemblages representative of neighbouring zoogeographical regions, characterized by these masses of water [9, 33, 34]. These waters, with subtropical and sub-Antarctic origins, affected the structure of parasite communities in the ecotone by acting as sources of infective stages of helminth species (acanthocephalans, nematodes, cestodes) typical of adjacent zoogeographical units, which were considered as reliable indicators to define such transitional regions. Among suitable markers, Lanfranchi et al.  reported larval Anisakis in Z. conchifer, identifying most of them as A. simplex (s.l.) based on morphology. These parasites were found at a prevalence of 77.3%, unexpectedly higher than the prevalences reported in more coastal fishes at similar latitudes. In the present paper, a genetic identification of a subsample of larval Anisakis from the same samples of Z. conchifer were carried out to assess the relative influence of sub-Antarctic and subtropical waters on the specific composition of this genus. Indeed, this region contained the distributional range of A. pegreffii; their occurrence in sympatry with A. typica should confirm the influence of the Brazil Current, as postulated by Lanfranchi et al. . In fact, the only Anisakis larva found in a sample of ten Z. conchifer from Rio de Janeiro, Brazil was genetically identified as A. typica .
Prevalence and mean abundance of larval Anisakis in the Argentine Sea, presumed to be mostly A. pegreffii, follow a latitudinal pattern increasing southwards irrespective of the host species harbouring them [9, 18]. The locality of capture of Z. conchifer in the ecotonal region also provided the opportunity for evaluating whether population attributes of larval anisakids follow a general distribution pattern across fish species with a similar trophic level in the region, or if they depart from it. Such a departure could be considered not only as evidence of the presence of different species of Anisakis with their own distribution patterns, but also indicative of the influence of warmer waters on the distribution of Anisakis spp. in the subtropical-sub-Antarctic convergence region.
The aim of this paper is, therefore, threefold: (i) to unequivocally identify larval Anisakis in the Argentine Sea based on molecular techniques; (ii) to assess their value as zoogeographical indicators in an ecotonal zone; and (iii) to determine the possible drivers of the distribution of Anisakis spp. across fish species of similar trophic levels in the Argentine Sea. Our results will also contribute to the knowledge and inventory of this genus in the World Ocean, filling a gap on the extant knowledge on the distribution and global diversity of Anisakis.
Fish sampling and parasite inventories
A total of 46 specimens of Z. conchifer were examined for larval Anisakis. Fish were caught by trawl during a research cruise at the Argentine-Uruguayan Common Fishing Zone (35°32'–35°35'S, 53°06'–53°25'W) at depths between 94 and 117 m, in October 2011 (Fig. 1). Data of most of these fish correspond to a previously published paper . Fish were either kept fresh or deep-frozen in plastic bags at -18 °C until examination. Females (n = 31) measured on average 28.4 cm (range 16.5–46.5 cm), males (n = 15) measured on average 24.5 cm (range 17.0–35.0 cm). After defrosting, larval Anisakis were recovered from the mesenteries, body cavity and liver after examination under a stereoscopic microscope. Prevalence and mean abundance were calculated following Bush et al. .
The identification at the species level was obtained by direct sequences analysis of mitochondrial (mtDNA cox2) and nuclear (nDNA EF1 α-1) genes.
A subsample of 19 Type I larvae of Anisakis spp. (equivalent to 18.6% of the total collected) and a unique Type II larva, both from Z. conchifer, along with 16 preserved Type I larvae from another two fish hosts caught during previous studies (Fig. 1), all from the Magellanic zoogeographical province (9 worms from Merluccius hubbsi; 44°4'S, 63°29'W; January 2009), and Argentine zoogeographical province (7 worms from Pagrus pagrus; 36°S, 55°W; March 2016), were identified to the species level by genetic analyses.
DNA extraction was carried out using whole specimens with the DNeasy Blood and Tissue® Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The mtDNA cox2 gene was amplified using the primers 210R (5'-CAC CAA CTC TTA AAA TTA TC-3') and 211F (5'-TTT TCT AGT TAT ATA GAT TGR TTY AT-3') . Additionally, to confirm their identity, the nDNA EF1 α-1 was amplified in larvae of A. pegreffii selected among those identified by the mtDNA cox2 gene using the primers EF-F (5'-TCC TCA AGC GTT GTT ATC TGT T-3') and EF-R (5'-AGT TTT GCC ACT AGC GGT TCC-3') . All PCR reactions were set up in a 25 μl reaction volume using 12.5 μl of HotStarTaq Master Mix (Qiagen), 0.5 μM of each primer and 5 μl of DNA (≥ 10 ng) as a template. The PCR was carried out using the following conditions: 94 °C for 15 min (initial heat activation) followed by 35 cycles at 94 °C for 30 s (denaturation), 50 °C for 100 s (annealing) and 72 °C for 1 min (extension), followed by a final extension step at 72 °C for 10 min. Each PCR product was purified using QIAquick spin columns (QIAquick Gel Extraction Kit, Qiagen). The fragments were sequenced for both DNA strands using the PCR primers. Sequencing was performed using Big Dye Terminator v.3.1 and 3130xl Genetic analyzer (Applied Biosystems, Foster City, CA, USA) at the Genomic Unit, IB-INTA.
Sequences were edited manually in Proseq v.3.5  and deposited in the GenBank database. Generated sequences at the mtDNA cox2 were compared against the NCBI database using the BLAST algorithm . The same sequences were aligned by ClustalW  implemented in the MEGA 7.0 software package  and compared with subsets of ten sequences available in GenBank of each known species of Anisakis, except for A. physeteris and A. cf. paggiae, of which only four and one sequences, respectively, were available for comparison. In the case of A. paggiae, those sequences deposited by Quiazon et al. [43, 44], considered as A. paggiae-related by the authors, were excluded due to specific uncertainty. The interspecific and intraspecific genetic distance between currently described Anisakis species and sequences obtained in the present study were calculated using a Kimura-2-Parameters (K2P) model in MEGA. The sequences obtained at the EF1 α-1 gene (409 bp fragment) of the nDNA for the specimens of A. pegreffii were compared at the diagnostic positions 186 and 286 as previously detailed in .
Distribution of Anisakis spp. across fish species in the Argentine Sea
Data on prevalence and mean abundance of larval Anisakis in other fish species of the region were obtained from previous publications by the research team during the last 20 years (Table 1). Since host trophic level has been recognized as an important determinant of the abundance of long-lived parasites , and because the interest of the analyses relies on the geographical distribution of Anisakis spp., only those fish species with a high trophic level (> 3.5, obtained from Froese & Pauly ) and ichthyophagous diet  were retained (Table 1). In this sense, by having similar diets and due to the non-specificity of larval Anisakis in previous paratenic fish hosts, all host species are expected to be exposed to the same pool of infective stages, diminishing the effect of their trophic level, and therefore allowing other variables to arise as determinants of parasite burdens.
Including data for Z. conchifer, a total of 18 samples corresponding to 9 fish species were analyzed. Each sample was assigned to a region following pre-established zoogeographical schemes ( and references therein); these regions were the Bonaerensean District of the Argentine Province and the Magellanic Province, both displaying characteristic parasite faunas . Samples caught at transitional areas between these two regions , namely San Matías Gulf and the outer shelf of the Bonaerensean region, influenced by sub-Antarctic waters which at these latitudes flow northwards along the slope bordering shelf waters, were assigned to a third region, defined as ecotone.
To analyze the relative contribution of host/abiotic variables on parasite distribution, Euclidean distance matrices of both prevalence and mean abundance were analyzed by distance-based multiple linear regressions (DistLM)  with significance testing based on 9999 permutations. As host-related predictor variables, the host species and their trophic level and mean total length were included in the models due to their known influence on parasite burdens [45, 50, 51]. Abiotic predictor variables were latitude, longitude and year of capture. The central year of the study period (2006) was adopted as the date of the sample Po. Draftsman plots and correlation matrices were used to check for multicollinearity in the predictor variables. Latitude and longitude were highly correlated each other (R = 0.93), due to the northeast to southwest orientation of the Argentine continental shelf, but were combined in a single predictor (locality).
Models including all possible combinations of predictor variables were generated using the Best procedure within the DistLM routine. An information theoretic approach based on modified Akaike’s information criterion (AICc) was used to identify the best model; models with the lowest AICc were considered the most parsimonious . The AICc was devised to handle situations where the number of samples (n) is small relative to the number (v) of predictor variables (n/v < 40) . The difference (Δi) between the AICc value of the best model and the AICc value for each of the other models was calculated; models with Δi between 0 and 2 are considered as having a substantial level of empirical support of the model being therefore as good as the best model ; however, as suggested by Richards , models with Δi up to 6 should not be discounted, thus all models with Δi ≤ 6 were retained. For each selected model, the Akaike weights (wi) were calculated following Burnham & Anderson  to identify and quantify the uncertainty in model selection and further used to estimate the relative importance of each predictor variable (predictor weight). For each predictor, the Akaike weights of all the models (with Δi < 6) that contained that predictor were summed and these values were interpreted as the relative importance of that predictor. Indeed, those predictors occurring consistently in the most likely models have a wi close to 1, whereas variables that are absent from or are only present in poorly fitting models (high AICc values) have a wi close to 0 . Additionally, the relative strength of each candidate model was assessed by calculating the evidence ratio (ER), which provides a measure of how much more likely the best model is than alternative models .
To visualize possible geographical patterns in the prevalence of larval Anisakis across the 18 samples, non-metric multidimensional scaling (nMDS) was carried out using Euclidean distances. A hierarchical agglomerative clustering was applied to the component communities using complete linkage, and resemblance levels were overlaid on the nMDS plot .
A total of 103 larval Anisakis were found parasitizing Z. conchifer. All of them were identified as Type I larva (prevalence: 78.3%; mean abundance: 2.2; range: 0–13), except for one specimen classified as Type II.
The identification through BLAST showed that sequences from the mtDNA cox2 gene of larvae Anisakis Type I exhibited a similarity of 99–100% with sequences for A. typica available on GenBank (n = 20; 16 from Z. conchifer, accession numbers MH443102-MH443117 and 4 from P. pagrus, accession numbers MH443118-MH443121); of 99% with those for A. berlandi (n = 3; 2 from Z. conchifer, accession numbers MH443122-MH443123 and 1 from P. pagrus, accession number MH443124); and of 99–100% with sequences for A. pegreffii (n = 12, 1 from Z. conchifer, accession number MH443127; 2 from P. pagrus, accession numbers MH443126-MH443123; and 9 from M. hubbsi, accession numbers MH443128-MH443136). The only larva Anisakis Type II from Z. conchifer (accession number MH443137) showed a similarity of 95% with the sequences for A. paggiae available on GenBank and of 99% with an undetermined species sequence, A. cf. paggiae.
Interspecific genetic distances (Table 2) ranged between 0.10–0.20, except for those corresponding to the A. simplex complex and the pair Anisakis sp.-A. cf. paggiae, which ranged between 0.05–0.06. The distances between the sequences from the present study and those from GenBank to which they showed the maximum similarity were all 0.01–0.02, the range of most of the observed intraspecific distance, thereby confirming their specific status.
In addition, the presence of a T at the position 186 and of a C at the position 286 of the partial sequences of the EF1 α-1 region of the nDNA obtained in 8 specimens, confirmed the identification of A. pegreffii (GenBank: numbers MH443138-MH443136).
Distribution of larval Anisakis across fish species in the Argentine Sea
The results of the DistLM on the prevalence data showed that the best model included host trophic level and locality as predictor variables (explaining 69.8% of the total variation of the data) (Table 3). The wi indicated that it has 69.0% chance of being the best model and ER showed that it was nearly five times more likely to be the best approximating model than the subsequent one. Indeed, trophic level and locality were included in all and most models with Δi < 6, reaching a predictor weight of 1 and 0.9, respectively, indicating that both variables had the highest probabilities of being a component of the best model (Fig. 2). Regarding mean abundance, a higher number (10) of alternative models were obtained, the best one composed only by locality as a predictor variable (explaining 38% of the total variation of the data) (Table 3). The wi indicated that the first model has a 22.8% chance of being the best one, a value very similar to that of the subsequent (composed by locality and trophic level). Evidence ratios showed that the first two models had similar chances of being the best one, but both were more than one and a half times more likely to be the best approximating models than the subsequent one. The predictor weights indicated that locality had the highest relative importance, followed by trophic level and year with considerably lower importance, whereas the mean host length had quite negligible relevance (Fig. 2). None of the models with Δi < 6 included host species as explanatory variable for the prevalence or mean abundance.
nMDS and cluster analyses on the prevalence data revealed apparent patterns of separation between samples following a zoogeographical pattern, which was substantially different from random as shown by it low stress level (0.01). Indeed, two main groups were evident (Fig. 3), one composed by most Bonaerensean samples and including samples Pb1 and Mh4 from the ecotone between Argentine and Magellanic Provinces and being more heterogeneous (branching at higher distances) than the second, which included the remaining Magellanic samples, but also the southernmost Bonaerensean Pb2 and Z. conchifer. A better picture of samples distribution is obtained by a three-dimensional nMDS (stress level = 0.01) (Additional file 1: Video S1).
Molecular analyses identified four Anisakis species in the whole sample, including A. typica, A. pegreffii, A. berlandi and an unidentified species, Anisakis sp., which seems to be conspecific with A. cf. paggiae from Kogia sima from off north-east Brazil , indicating the possible presence of a new species of Anisakis in the South-West Atlantic. Indeed, the K2P distance between A. paggiae and Anisakis sp. was similar to those between the sibling species of the A. simplex complex. Similar distances have been reported for A. paggiae-related species in Japanese , Philippine  and Brazilian  waters. Our data support the existence of an A. paggiae species complex as suggested by these authors. It is noteworthy that this specimen of Anisakis sp. is the only larva Type II so far recorded by the authors in the region, not only in the present samples, but in many other fish species previously surveyed (personal observation).
Knowledge of geographical distribution patterns of Anisakis spp. is scarce for many species in the genus . Indeed, most species have been reported from temperate, subtropical and tropical waters between the equator and 35° north and south, being apparently more common in the boreal region . The present findings increase the list of Anisakis species in marine fish from the West Atlantic, south to 35°S, by reporting for the first time A. typica and Anisakis sp. Together with the previously reported A. pegreffii and A. berlandi, this results in a regional richness of four species, a richness smaller than that recorded in Brazilian marine waters, where six of the nine known species of the genus [i.e. excepting A. simplex (s.s.), A. pegreffii and A. berlandi] have been recorded . The record of larvae of A. pegreffii in a specimen of Thunnus thynnus from the region of Rio de Janeiro by Mattiucci et al.  is clearly a typing mistake (Table 3 of Mattiucci et al. ), which was subsequently replicated .
The composition of species observed in the study region reflects, furthermore, the transitional nature of its oceanographic conditions. Anisakis typica is known to be restricted to subtropical and tropical waters, and is the most common species in Brazilian waters , whereas A. pegreffii displays a discontinuous distribution, being known from the Mediterranean Sea, the Central East Atlantic and waters of China and Japan in the Northern Hemisphere and displaying a circumpolar distribution in the Southern Hemisphere [4, 24, 59]. Anisakis pegreffii is also the most reported species in the southern Argentine Sea, dominated by sub-Antarctic waters where 121 of 122 previously genetically identified worms were assigned to this species and only one to A. berlandi [18, 60, 61]. Therefore, A. typica and A. pegreffii can be considered as representatives of northern and southern regions of the South-West Atlantic, respectively. Despite being little represented in the examined samples, A. berlandi and Anisakis sp. also seem to have opposite geographical origins in the study region. Anisakis berlandi is typical from southern waters in the Southern Hemisphere [59, 61], although is also present in Pacific Canadian and Californian waters . On the other hand, Anisakis sp. is apparently conspecific with A. cf. paggiae from Kogia sima from off north-east Brazil , a cetacean mainly distributed in deep-water habitats of tropical and temperate zones in the Central Atlantic Ocean . In the marine environment, transitional zones between major water masses harbour high biodiversity, mostly due to their productivity , but also to an “edge effect”, as defined by Odum . Indeed, ecotones contain representatives of species characteristic of adjacent communities , which seems to explain the specific composition of Anisakis in the convergence regions where masses of water supply “northern” and “southern” species. The thermal gradient produced by these marine currents in the South-West Atlantic  explains the dominance of A. pegreffii in M. hubbsi off Patagonia (southern waters) and of A. typica in Z. conchifer. The presence of A. typica in P. pagrus caught at similar latitudes, but in coastal waters, evidences that the influence of the Brazil Current extends to shallower areas surrounding the convergence. This sparid is associated to hard substrates and exhibits considerable site fidelity, remaining in the same patch after recruitment , being therefore only exposed to infective stages present in its habitat.
In addition to the environmental conditions as determinants of distribution on Anisakis species, their distribution patterns at large spatial scales are congruent with those of their respective final hosts . Anisakis typica has been recorded in subtropical and tropical waters from several species of dolphins (Delphinidae), but also in the harbour porpoise, Phocoena phocoena (Phocoenidae) and the franciscana dolphin, Pontoporia blainvillei (Pontoporidae) [4, 56, 58, 68]. Most of these cetaceans are typical of warmer temperate and tropical seas  and consequently the distribution of the parasite can be largely determined by that of their definitive hosts. However, some of these dolphins, namely Tursiops truncatus and Delphinus delphis are also distributed in higher latitudes, reaching the colder Patagonian waters, where larvae of A. typica have not been recorded yet in fishes, and where the adults paraizing them have been morphologically identified as A. simplex (sensu lato) [19, 69, 70]. Similarly, P. blainvillei is infested by A. typica in Uruguayan waters, but harbours A. simplex (s.l.) in central Argentine waters . On the other hand, A. pegreffii has been never reported in Brazilian waters. Therefore, the presence of suitable definitive hosts is a necessary but not sufficient condition to explain the geographical patterns of Anisakis spp., evidencing that environmental conditions play a major role in such a distribution. In agreement with most records in fish and cetacean hosts from lower latitudes , A. typica dominated the anisakid community parasitizing Z. conchifer in waters on the continental slope at depths between 94 and 117 m and 5° south of its known distribution limit, evidencing a higher effect of the warmer Brazil current relative to other masses of water on the local distribution of larval Anisakis. The influence of subtropical waters could, therefore, occur through the transport of infective stages in previous intermediate hosts or be related to the migratory behaviour of the definitive hosts or infested fish from lower latitudes. However, the possibility of A. typica arriving to the study zone with migratory Z. conchifer from northern regions seems unlikely. Although this species seems to perform ontogenetic movements to deeper waters where large fish concentrate for reproduction  the occurrence of latitudinal migrations has not been reported. Furthermore, the size-at-maturity of females is 311 mm (that of males is unknown) , which indicates that most of the examined specimens were juvenile and consequently, based on the available information, they are most probably non-migratory individuals. The presence of A. typica in P. pagrus, as mentioned above, supports the transport of infective stages by the Brazil Current as the main cause of the southern extension in the range of this parasite.
Regarding the distribution of larval Anisakis across fish species south of 35°S, the explanatory variables mean host length and host species had little value as drivers of the prevalence and mean abundance of parasites. The lack of relationship between host species and parasite burdens is relevant since it indicates that, given the extremely low specificity of larval Anisakis, fish hosts act as passive samplers of infective stages through their diets. The identity and load of larvae in fish, therefore, depend on the trophic behaviour of fish hosts, as well as on the availability of larvae in the food web, modelled by environmental conditions, in each region. Indeed, despite only fish species with a high trophic level were considered, this variable was, together with locality, the main determinant of prevalence, being also highly relevant in explaining the mean abundance of larval Anisakis.
Larval anisakids are long-lived in fish hosts, and can be transmitted from one paratenic host to other, persisting in the food web and being potentially available for any host, independent of its trophic level. However, trophic level determines the range of preys a predator can consume, having a direct influence on the abundance and composition of parasite communities [73, 74]. In fact, ichthyophagous fish tend to accumulate higher numbers of larval parasites because they acquire packets of helminth species that travel together in paratenic hosts along food chains . Therefore, the higher their trophic level, the higher the likelihood of becoming infected and the larger the number of infective stages acquired with each individual prey , determining the prevalence and the mean abundance of Anisakis spp. regardless of the locality of capture. The year of capture also had certain relevance in determining mean abundance, indicating that some changes could have occurred in the region during recent decades; however, a reliable explanation to these variations would require additional studies, beyond the scope of the present work. Finally, and as observed for larval Anisakis in skates of the Argentine Sea , the geographical origin of the samples was the main determinant of prevalence and mean abundance of parasites. The environmental conditions in the study region are determined by the subtropical and sub-Antarctic currents flowing along the continental slope . The Malvinas Current dominates adjacent shelf waters, producing a latitudinal cline of temperature which decreases southwards, whereas the effect of the warm Brazil Current is marked at the northern limit of the Argentine sea . This temperature gradient, characteristic of the study area , provides a series of environments in which Anisakis larvae are differentially distributed in terms of prevalence and abundance.
DistLM analyses were carried out irrespective of the Anisakis species comprising the assemblages of each sample. However, nMDS analyses showed that Z. conchifer departed from the general zoogeographical pattern displayed by the other samples, which in turn, agreed with biogeographical schemes recognized in the South-West Atlantic . Indeed, the apparent separation of samples in two groups corresponding to the Bonaerensean District of the Argentine Province and to the Magellanic Province, characterized by low and high values of prevalence, respectively, confirms the value of larval Anisakis as zoogeographical indicators [9, 34]. On the other hand, the two ecotonal hake samples between these two zoogeographical units (Mh3 and Mh4) clustered with Magellanic and Bonaerensean samples, respectively, indicating a higher influence of each mass of water on the populations of Anisakis in fish at these transitional zones.
The exception to this pattern was the assignation of the southernmost Bonaerensean sample of P. brasiliensis (Pb2) to the Magellanic group of samples. This sample, together with that from the ecotonal Pb1 assigned to the Bonaerensean group, displayed the higher prevalence of larval Anisakis for this host. Although they were considered as different stocks based on their parasite assemblages , Avigliano et al.  determined, based on otholith microchemistry, that P. brasilisensis from these two localities represent a single stock. In the light of the present results, migratory movements of these P. brasiliensis between these latitudes could not be disregarded. Finally, Z. conchifer from deep waters at the convergence between subtropical and sub-Antarctic currents grouped with the distant Magellanic samples as a consequence of the unexpectedly high prevalence regarding other northern samples, even some very close ones, but caught on shelf waters. Taking into account that A. typica was the dominant species in Z. conchifer, and assuming that other samples are mostly composed of A. pegreffii, these results could indicate that interspecific differences of populations distribution exist between these congeners. This hypothesis requires further research because no other host species with a high trophic level from a similar latitude and depth have yet been analysed. In case these differences are confirmed in other fish samples, they may be useful tools as surrogates of the identity of larvae parasitizing a given host population.
The specific identification of larval Anisakis in a region where they are considered among the best biological indicators for studies on zoogeography and population distribution of their hosts is undoubtedly a step forward towards the validation of the use of this methodology. Furthermore, considering the diversity of species found, new perspectives for future studies also arise, since members of Anisakis have proven to be excellent tools for studying host-parasite cophylogeny, as indicators of trophic web stability and indicators of habitat disturbance of marine ecosystems [4, 24].
Knowledge on the distribution of Anisakis spp. in the World Ocean is indispensable for future studies on the epidemiology and pathogenicity of anisakiosis, as well as on the possibility of a changing risk of this zoonosis in the time of climate change .
The genetic identification of four species of larval Anisakis in fishes from the southern South-West Atlantic (three known and a probably new species related to A. paggiae) fill a gap in the knowledge of the global distribution of these zoonotic parasites. It was confirmed that, in the study region, larval Anisakis follow a clear zoogeographical pattern, being suitable indicators of tropical/subtropical and sub-Antarctic waters; the main drivers of that pattern were the host trophic level and locality of capture. This information is relevant for both human health and fishery industry, since the species in the genus exhibit differential levels of pathogenicity. The higher taxonomic resolution reached by molecular techniques represents a step forward the validation of the use of larval Anisakis as biological indicators for studies on host zoogeography.
Modified Akaike’s information criterion
- cox2 :
Cytochrome c oxidase subunit 2
Distance-based linear model
- EF1 α-1:
Elongation factor α-1
Kimura 2-parameter model
Non-metric multidimensional scaling
Smith JW, Wootten R. Anisakis and anisakiasis. Adv Parasitol. 1978;16:93–163.
Audicana MT, Ansotegui IJ, de Corres LF, Kennedy MW. Anisakis simplex: Dangerous-dead and alive? Trends Parasitol. 2002;18:20–5.
Audicana MT, Kennedy MW. Anisakis simplex: from obscure infectious worm to inducer of immune hypersensitivity. Clin Microbiol Rev. 2008;21:360–79.
Mattiucci S, Nascetti G. Advances and trends in the molecular systematics of anisakid nematodes, with implications for their evolutionary ecology and host-parasite co-evolutionary processes. Adv Parasitol. 2008;66:47–148.
McCarthy J, Moore TA. Emerging helminth zoonoses. Int J Parasitol. 2000;30:1351–60.
Timi JT. Parasites as biological tags for stock discrimination in marine fish from South American Atlantic waters. J Helminthol. 2007;81:107–11.
Lester RJG, MacKenzie K. The use and abuse of parasites as stock markers for fish. Fish Res. 2009;97:1–2.
Catalano SR, Whittington ID, Donnellan SC, Gillanders BM. Parasites as biological tags to assess host population structure: guidelines, recent genetic advances and comments on a holistic approach. Int J Parasitol Parasites Wildl. 2013;3:220–6.
Cantatore DMP, Timi JT. Marine parasites as biological tags in South American Atlantic waters, current status and perspectives. Parasitology. 2015;142:5–24.
Berland B. Nematodes from some Norwegian marine fishes. Sarsia. 1961;2:1–50.
Cipriani P, Smaldone G, Acerra V, D’Angelo L, Anastasio A, Bellisario B, et al. Genetic identification and distribution of the parasitic larvae of Anisakis pegreffii and Anisakis simplex (s.s.) in European hake Merluccius merluccius from the Tyrrhenian Sea and Spanish Atlantic coast: implications for food safety. Int J Food Microbiol. 2015;198:1–8.
Mattiucci S, Abaunza P, Ramadori L, Nascetti G. Genetic identification of Anisakis larvae in European hake from Atlantic and Mediterranean waters for stock recognition. J Fish Biol. 2004;65:495–510.
Mattiucci S, Garcia A, Cipriani P, Santos MN, Nascetti G, Cimmaruta R. Metazoan parasite infection in the swordfish, Xiphias gladius, from the Mediterranean Sea and comparison with Atlantic populations: implications for its stock characterization. Parasite. 2014;21:35.
Mattiucci S, Cimmaruta R, Cipriani P, Abaunza P, Bellisario B, Nascetti G. 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. Parasitology. 2015;142:90–108.
Abaunza P, Murta AG, Campbell N, Cimmaruta R, Comesana AS, Dahle G, et al. Stock identity of horse mackerel (Trachurus trachurus) in the Northeast Atlantic and Mediterranean Sea: integrating the results from different stock identification approaches. Fish Res. 2008;89:196–209.
Kuhn T, García-Márquez J, Klimpel S. Adaptive radiation within marine anisakid nematodes: a zoogeographical modelling of cosmopolitan, zoonotic parasites. PLoS One. 2011;6:e28642.
Gómez-Mateos M, Valero A, Morales-Yuste M, Martín-Sánchez J. Molecular epidemiology and risk factors for Anisakis simplex s.l. infection in blue whiting (Micromesistius poutassou) in a confluence zone of the Atlantic and Mediterranean: differences between A. simplex s.s. and A. pegreffii. Int J Food Microbiol. 2016;232:111–6.
Irigoitia MM, Braicovich PE, Lanfranchi AL, Farber MD, Timi JT. Distribution of anisakid nematodes parasitizing rajiform skates under commercial exploitation in the Southwestern Atlantic. Int J Food Microbiol. 2018;267:20–8.
Hernández-Orts JS, Paso Viola MN, García NA, Crespo EA, González R, García-Varela M, Kuchta R. A checklist of the helminth parasites of marine mammals from Argentina. Zootaxa. 2015;3936:301–34.
Berón-Vera B, Crespo EA, Raga JA. Parasites in stranded cetaceans of Patagonia. J Parasitol. 2008;94:946–8.
Luque JL, Muniz-Pereira LC, Siciliano S, Siqueira LR, Oliveira MS, Vieira FM. Checklist of helminth parasites of cetaceans from Brazil. Zootaxa. 2010;2548:57–68.
Fonseca MCG, Knoff M, Felizardo NN, Di Azevedo MIN, Torres EJL, Gomes DC, et al. Integrative taxonomy of Anisakidae and Raphidascarididae (Nematoda) in Paralichthys patagonicus and Xystreurys rasile (Pisces: Teleostei) from Brazil. Int J Food Microbiol. 2016;235:113–24.
Di Azevedo MIN, Iñiguez AM. Nematode parasites of commercially important fish from the southeast coast of Brazil: morphological and genetic insight. Int J Food Microbiol. 2018;267:29–41.
Mattiucci S, Paoletti M, Cipriani P, Webb SC, Timi JT, Nascetti G. Inventorying biodiversity of anisakid nematodes from the Austral Region: a hotspot of genetic diversity? In: Klimpel S, Kuhn T, Mehlhorn H, editors. Biodiversity and Evolution of Parasitic Life in the Southern Ocean. Parasitology Research Monographs, vol. 9. Cham: Springer; 2017. p. 109–40.
Piola AR, Martínez Avellaneda N, Guerrero RA, Jardón FP, Palma ED, Romero SI. Malvinas-slope water intrusions on the northern Patagonia continental shelf. Ocean Sci. 2010;6:345–59.
Moreno IB, Zerbini AN, Danilewicz D, de Oliveira Santos MC, Simões-Lopes PC, Lailson-Brito J Jr, Azevedo AF. Distribution and habitat characteristics of dolphins of the genus Stenella (Cetacea: Delphinidae) in the southwest Atlantic Ocean. Mar Ecol Prog Ser. 2005;300:229–40.
Mandiola MA, Giardino GV, Bastida J, Rodríguez DH, Bastida RO. Marine mammal occurrence in deep waters of the Brazil-Malvinas Confluence off Argentina during summer. Mastozool Neotrop. 2015;22:397–402.
Di Tullio JC, Gandra TBR, Zerbini AN, Secchi ER. Diversity and distribution patterns of cetaceans in the subtropical southwestern Atlantic outer continental shelf and slope. PLoS One. 2016;11:e0155841.
Lanfranchi AL, Braicovich PE, Cantatore DMP, Alarcos AJ, Luque JL, Timi JT. Ecotonal marine regions - ecotonal parasite communities: helminth assemblages in the convergence of masses of water in the southwestern Atlantic. Int J Parasitol. 2016;46:809–18.
Haimovici M, Martins AS, Figueiredo JL, Vieira PC. Demersal bony fish of the outer shelf and upper slope of the southern Brazil Subtropical Convergence ecosystem. Mar Ecol Prog Ser. 1994;108:59–77.
Quéro JC, Du Buit MH, Vayne JJ. Les observations de poisons tropicaux et le réchauffement des eaux dans l’Atlantique européen. Oceanol Acta. 1998;21:345–51.
Ragonese S, Giusto GB. Zenopsis conchifera (Lowe, 1852) (Pisces, Actinopterygii, Zeidae): a new alien fish in the Mediterranean Sea. J Fish Biol. 2007;21:345–51.
Pereira AN, Pantoja C, Luque JL, Timi JT. Parasites of Urophycis brasiliensis (Gadiformes: Phycidae) as indicators of marine ecoregions in coastal areas of the South American Atlantic with the assessment of their stocks. Parasitol Res. 2014;113:4281–92.
Braicovich PE, Pantoja C, Pereira AN, Luque JL, Timi JT. Parasites of the Brazilian flathead Percophis brasiliensis reflect West Atlantic biogeograhic regions. Parasitology. 2017;144:169–78.
dos Reis Sardella CJ, Luque JL. Diagnóstico morfológico e molecular de larvas de Anisakis typica e Anisakis brevispiculata em peixes do litoral do Rio de Janeiro. Braz J Vet Med. 2016;38(Supl. 3):124–30.
Bush AO, Lafferty KD, Lotz JM, Shostak AW. Parasitology meets ecology on its own terms: Margolis et al. revisited. J Parasitol. 1997;83:575–83.
Nadler SA, Hudspeth DSS. Phylogeny of the Ascaridoidea (Nematoda: Ascaridida) based on three genes and morphology: hypothesis of structural and sequence evolution. J Parasitol. 2000;86:380–93.
Mattiucci S, Acerra V, Paoletti M, Cipriani P, Levsen A, Webb SC, et al. No more time to stay ‘single’ in the detection of Anisakis pegreffii, A. simplex (s.s.) and hybridization events between them: a multi-marker nuclear genotyping approach. Parasitology. 2016;143:998–1011.
Filatov DA. ProSeq: a software for preparation and evolutionary analysis of DNA sequence data sets. Mol Ecol Notes. 2002;2:621–4.
Altschul FS, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Quiazon K, Yoshinaga T, Santos M, Ogawa K. Identification of larval Anisakis spp. (Nematoda: Anisakidae) in Alaska pollock (Theragra chalcogramma) in northern Japan using morphological and molecular markers. J Parasitol. 2009;95:1227–32.
Quiazon K, Santos M, Yoshinaga T. Anisakis species (Nematoda: Anisakidae) of dwarf sperm whale Kogia sima (Owen, 1866) stranded off the Pacific coast of southern Philippine archipelago. Vet Parasitol. 2013;197:221–30.
Timi JT, Rossin MA, Alarcos AJ, Braicovich PE, Cantatore DMP, Lanfranchi AL. Fish trophic level and the similarity of larval parasite assemblages. Int J Parasitol. 2011;41:309–16.
Froese R, Pauly D. Fish Base. In: World Wide Web Electronic Publication; 2018. www.fishbase.org. Accessed 28 Feb 2018.
Cousseau MB, Perrotta RG. Peces Marinos de Argentina: Biología, Distribución, Pesca. 4th ed. Mar del Plata: Instituto Nacional de Investigación y Desarrollo Pesquero; 2013.
Menni RC, Jaureguizar AJ, Stehmann MFW, Lucifora LO. Marine biodiversity at the community level: zoogeography of sharks, skates, rays and chimaeras in the southwestern Atlantic. Biodivers Conserv. 2010;19:775–96.
Anderson MJ, Gorley RN, Clarke KR. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods. Plymouth: PRIMER-E; 2008.
Timi JT, Lanfranchi AL. Ontogenetic changes in heterogeneity of parasite communities of fish: disentangling the relative role of compositional versus abundance variability. Parasitology. 2013;140:309–17.
Braicovich PE, Ieno EN, Sáez M, Despos J, Timi JT. Assessing the role of host traits as drivers of the abundance of long-lived parasites in fish stock assessment studies. J Fish Biol. 2016;89:2419–33.
Symonds MRE, Moussalli A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav Ecol Sociobiol. 2011;65:13–21.
Burnham KP, Anderson DR. Model Selection and Multimodel Inference. 2nd ed. New York: Springer; 2002.
Richards SA. Testing ecological theory using the information theoretic approach: examples and cautionary results. Ecology. 2005;86:2805–14.
Clarke KR, Gorley RN. PRIMER v7: User Manual/Tutorial. Plymouth: PRIMER-E; 2015.
Di Azevedo MIN, Carvalho VL, Iñiguez AM. Integrative taxonomy of anisakid nematodes in stranded cetaceans from Brazilian waters: an update on parasite’s hosts and geographical records. Parasitol Res. 2017;116:3105–16.
Klimpel S, Kuhn T, Busch MW, Horst K, Palm HW. Deep water life-cycle of Anisakis paggiae (Nematoda: Anisakidae) in the Irminger Sea indicates kogiid whale distribution in North Atlantic waters. Polar Biol. 2011;34:899–906.
Mattiucci S, Paggi L, Nascetti G, Santos CP, Costa G, Di Beneditto AP, et al. Genetic markers in the study of Anisakis typica (Diesing, 1860): larval identification and genetic relationships with other species of Anisakis Dujardin, 1845 (Nematoda: Anisakidae). Syst Parasitol. 2002;51:159–70.
Klimpel S, Palm HW. Anisakid nematode (Ascaridoidea) life cycles and distribution: increasing zoonotic potential in the time of climate change? In: Mehlhorn H, editor. Progress in Parasitology. Parasitology Research Monographs. Heidelberg: Springer Verlag; 2011;2:201–22.
Mattiucci S, Nascetti G, Cianchi R, Paggi L, Arduino P, Margolis L, et al. Genetic and ecological data on the Anisakis simplex complex with evidence for a new species (Nematoda, Ascaridoidea, Anisakidae). J Parasitol. 1997;83:401–16.
Mattiucci S, Nascetti G. Genetic diversity and infection levels of anisakid nematodes parasitic in fish and marine mammals from Boreal and Austral hemispheres. Vet Parasitol. 2007;66:47–148.
Klimpel S, Busch M, Khun T, Rohde A, Palm HW. The Anisakis simplex complex off the South Shetland Islands (Antarctica): endemic populations versus introduction through migratory hosts. Mar Ecol Prog Ser. 2010;34:899–906.
Scales KL, Miller PI, Hawkes LA, Ingram SN, Sims DW, Votier SC. On the front line: frontal zones as priority at-sea conservation areas for mobile marine vertebrates. J Appl Ecol. 2014;51:1575–83.
Odum EP. Fundamentals of Ecology. 2nd ed. Philadelphia: Saunders; 1959.
Baker J, French K, Whelan RJ. The edge effect and ecotonal species: bird communities across a natural edge in southeastern Australia. Ecology. 2002;83:3048–59.
Hoffmann J, Núñez M, Piccolo M. Características climáticas del océano Atlántico sudoccidental. In: Boschi EE, editor. El Mar Argentino y sus Recursos Pesqueros Tomo I. Antecedentes históricos de las exploraciones en el mar y las características ambientales. Mar del Plata: Instituto Nacional de Investigación y Desarrollo Pesquero; 1997. p. 163–93.
De Vries AD. The life history, reproductive ecology and demography of the red porgy, Pagrus pagrus, in the northeastern Gulf of Mexico. PhD Dissertation. Tallahassee: Florida State University; 2006.
Iñiguez AM, Carvalho VL, Motta MR, Pinheiro DC, Vicente AC. Genetic analysis of Anisakis typica (Nematoda: Anisakidae) from cetaceans of the northeast coast of Brazil: new data on its definitive hosts. Vet Parasitol. 2011;178:293–9.
Berón-Vera B, Crespo EA, Raga JA, Fernández M. Parasites communities of common dolphins (Delphinus delphis) from Patagonia: the relation with host distribution and diet and comparison with sympatric host. J Parasitol. 2007;93:1056–60.
Romero MA, Fernández M, Dans SL, García NA, González R, Crespo EA. Gastrointestinal parasites of bottlenose dolphins Tursiops truncatus from the extreme south western Atlantic, with notes on diet composition. Dis Aquat Org. 2014;108:61–70.
Aznar FJ, Raga JA, Corcuera J, Monzón F. Helminths as biological tags for franciscana (Pontoporia blainvillei) (Cetacea, Pontoporiidae) in Argentinian and Uruguayan waters. Mammalia. 1995;59:427–36.
Martins RS, Schwingel PR. Biological aspects of the sailfin dory Zenopsis conchifer (Lowe, 1852) caught by deep-sea trawling fishery off southern Brazil. Braz J Oceanogr. 2012;60:171–9.
Marcogliese DJ. Food webs and the transmission of parasites to marine fish. Parasitology. 2002;124:83–99.
Chen H-W, Liu W-C, Davis AJ, Jordan F, Hwang MJ, Shao K-T. Network position of hosts in food webs and their parasite diversity. Oikos. 2008;117:1847–55.
Piola AR, Rivas AL. Corrientes en la plataforma continental. In: Boschi EE, editor. El Mar Argentino y Sus Recursos Pesqueros Tomo I. Antecedentes históricos de las exploraciones en el mar y las características ambientales. Mar del Plata: Instituto Nacional de Investigación y Desarrollo Pesquero; 1997. p. 119–32.
Guerrero RA, Piola AR. Masas de agua en la plataforma continental. In: Boschi EE, editor. El Mar Argentino y sus Recursos Pesqueros Tomo I. Antecedentes históricos de las exploraciones en el mar y las características ambientales. Mar del Plata: Instituto Nacional de Investigación y Desarrollo Pesquero; 1997. p. 107–18.
Braicovich PE, Timi JT. Parasites as biological tags for stock discrimination of the Brazilian flathead Percophis brasiliensis in the south-west Atlantic. J Fish Biol. 2008;73:557–71.
Avigliano E, Saez MB, Rico R, Volpedo AV. Use of otolith strontium: calcium and zinc: calcium ratios as an indicator of the habitat of Percophis brasiliensis Quoy & Gaimard, 1825 in the southwestern Atlantic Ocean. Neotrop Ichthyol. 2015;13:187–94.
Sardella NH, Timi JT. Parasites of Argentine hake in the Argentine Sea: population and infracommunity structure as evidence for host stock discrimination. J Fish Biol. 2004;65:1472–88.
Alarcos AJ, Timi JT. Parasite communities in three sympatric flounder species (Pleuronectiformes: Paralichthyidae). Similar ecological filters driving toward repeatable assemblages. Parasitol Res. 2012;110:2155–66.
Alarcos AJ, Etchegoin JA. Parasite assemblages of estuarine-dependent marine fishes from Mar Chiquita coastal lagoon (Buenos Aires Province, Argentina). Parasitol Res. 2010;107:1083–91.
Timi JT, Lanfranchi AL. The metazoan parasite communities of the Argentinean sandperch Pseudopercis semifasciata (Pisces: Perciformes) and their use to elucidate the stock structure of the host. Parasitology. 2009;136:1083–91.
We thank the Instituto Nacional de Investigación y Desarrollo Pesquero (INIDEP), Mar del Plata, Argentina, for providing fish samples.
Financial support was provided by grants from CONICET (PIP # 112-201501-009763), ANPCYT (PICT 2015 # 2013) and UNMdP (Exa 915/18).
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional file. The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
All applicable international, national and/or institutional guidelines for the care and use of animals were followed.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lanfranchi, A.L., Braicovich, P.E., Cantatore, D.M.P. et al. Influence of confluent marine currents in an ecotonal region of the South-West Atlantic on the distribution of larval anisakids (Nematoda: Anisakidae). Parasites Vectors 11, 583 (2018). https://doi.org/10.1186/s13071-018-3119-7
- Anisakis pegreffii
- Anisakis berlandi
- Anisakis typica
- Zenopsis conchifer
- Zoogeographical indicators