Open Access
Issue
Parasite
Volume 33, 2026
Article Number 39
Number of page(s) 10
DOI https://doi.org/10.1051/parasite/2026038
Published online 27 July 2026

© A. Delhomme et al., published by EDP Sciences, 2026

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Introduction

Anisakidae Railliet & Henry, 1912 is a family of zoonotic parasitic nematodes that infect mainly marine animals [28]. This family includes genera such as Anisakis, Contracaecum, and Phocanema (ex-Pseudoterranova [3]), which are the most important in terms of zoonotic significance. All of these have received significant interest because they are ubiquitous parasites of marine fish species and because the consumption of raw or uncooked infected fish may generate hazards for human consumers [28]. The lifecycle of Anisakidae includes several hosts. Eggs are released into the environment with faeces of definitive hosts, cetaceans, pinnipeds, or piscivorous birds and hatch in the water. In the water column, larvae are consumed by zooplanktonic crustaceans. Zooplanktonic crustaceans are considered either paratenic or intermediate hosts (with development of the larva in the host), where larvae occur at L2 or L3 stages [25, 49]. Most frequently, fish and cephalopods are then infested, whether by consuming zooplankton, for zooplankton feeding species [7], or by consuming infested fish of lower trophic levels, for piscivorous fish species. Fish and cephalopods are mainly considered to be paratenic hosts, meaning that no morphological transformation of Anisakidae occurs in these species. The cycle is then completed when infested fish are preyed upon by a suitable definitive host, mainly cetaceans, pinnipeds, or piscivorous birds. Human infestation by Anisakidae occurs through the consumption of raw or undercooked infected seafood. Anisakidosis symptoms are mainly digestive or allergic. Most cases are observed in countries with high levels of fish consumption [45]. There is consequently a human-health concern related to increased consumption of infected and raw or undercooked fish products [9]. Due to this risk and given the need for infection data in fish to support risk assessment, documenting levels of infection by Anisakidae in fish has benefited from intense research activity. Similarly, as marine mammals are protected and, in some cases, endangered species, and as parasitism may be one of the potential causes of mortality (or a factor of morbidity for this group), parasites and specifically Anisakidae have also been screened for in stranded individuals [10, 26, 33], but with a lower effort than in fish.

Nonetheless, addressing the importance of the other members of the lifecycle, i.e. zooplanktonic crustaceans, is needed to better understand the ecological mechanisms underlying infection. In the context of implementation of the “One Health” approach, which calls for more integrated studies of zoonoses, understanding and predicting hazards related to the consumption of infected fish requires further study, and consideration of the other members of the cycle [9]. However, as illustrated by Poulin et al. [43], integrating parasitology into an ecological framework is a major challenge. Most studies are driven by human health related issues, and therefore focus on edible species in order to limit infection hazards for human consumers and overlook other members of the cycle [47]. Levels of Anisakidae infection in zooplankton have consequently been largely neglected, despite the importance of this group as a driver of fish infection and as a potential bottleneck in the cycle of Anisakidae [13, 32, 33, 42]. These studies mostly describe Euphausiacea as potential hosts of Anisakidae [15, 32, 4951], even though Anisakidae larvae have also been observed in other groups, like Mysidacea and Copepoda. Ontogenic variation of the infection in jack mackerel Trachurus trachurus is, for example, directly related to ontogenic changes in diet, and the inclusion of Euphausiacea in the diet of larger fish individuals [7]. In the few studies focused on Anisakidae infection in Euphausiacea, prevalence values are generally low, and lower than 0.005%. Values ranging between 0 and 4% are nonetheless reported ([32], Table 1). A very peculiar case, with an extremely high prevalence of 40%, and “an exceptionally high rate of 78.0%” was recorded for Anisakis simplex larvae occurring in the Euphausiacea Thysanoessa spp. in the North Sea off Scotland and the Faroe Islands [50]. This value should, nevertheless, be considered with caution, particularly as the number of Euphausiacea individuals is not provided. But, as noted by Marcogliese [32], there is a major discrepancy in prevalence levels between zooplankton and fish, with levels being lower in zooplankton by several orders of magnitude. He highlighted this discrepancy as an interesting paradox that has to be elucidated in order to understand the ecological mechanisms of parasite transfer, and proposed hypotheses related either to massive consumption of zooplankton or to areas or periods where infection is favoured. The scarcity of data about Anisakidae infection in zooplankton did not allow any firm conclusions to be drawn.

Table 1

Existing data on the infection of zooplanktonic crustaceans by Anisakidae larvae. nd: no data provided in the original article.

One of the major challenges in determining Anisakidae infection levels in zooplankton is largely technical. The studies previously cited are based on visual observation of freshly collected samples of zooplankton, whether through binoculars or microscopes, and followed by zooplankton dissection and molecular identification of the parasite species, for the most recent paper [15]. This approach is associated with several drawbacks. Firstly, sorting, identifying, and observing numerous zooplankton individuals is time-consuming. This is a major issue if observations are to be generalised, over large spatial or temporal scales, as the low prevalence requires the collection of several hundred or even thousands of zooplankton individuals to gain a robust picture of infection patterns (Table 1). In addition, visual observation has to be performed immediately after collection. Most zooplankton preservation methods, whether freezing or chemical fixation with formalin-derived solutions, make them opaque and the larvae can no longer be seen through them [13, 50]. These types of observations are, therefore, limited to coastal sampling, where it is possible to return to the laboratory on land immediately after sampling. For samples taken during offshore campaigns, microscopic observation on board can be complicated by rough sea conditions. In addition, zooplankton is not the focus group of most surveys, usually dedicated to the collection of data about fish to ensure fisheries sustainability. Furthermore, sea campaigns are rare opportunities. They often involve a great deal of work, carried out by small teams. As a result, it is usually not possible to devote a significant amount of time to ancillary work on zooplankton, as time between two plankton-sampling operations is limited. Progressing towards more exhaustive investigation of Anisakidae levels in zooplankton therefore requires new methodological developments [32].

Furthermore, the hazards associated with eating infected fish have led to the development and standardisation of robust methods for collecting, detecting, and identifying Anisakidae larvae in fish. Fish muscle is either crushed under a hydraulic press, and then observed under UV light [20] or dissolved in pepsin solution [30]. Then, after preliminary identification based on morphological features, quantitative or conventional PCR followed by Sanger analysis allows for molecular identification of the larvae to the species level [11, 38, 46].

In this context, the aim of the present study was to adapt the molecular method developed to detect Anisakidae in fish for Euphausiacea. Numerous constraints, such as the relative size between Euphausiacea and the parasite inside, the impossibility of knowing whether there will be UV fluorescence, the opacity of the Euphausiacea or the difference in matrix, preclude direct application of the methods used in fish to Euphausiacea. It was therefore necessary to optimise all the steps to produce a protocol with sufficiently low sensitivity and high specificity for the detection of Anisakidae. The method also has to be simple to use and fast enough to process a large number of samples, as required by the expected low prevalence, and in order to apply the method over large spatial and temporal scales. Development was carried out in two steps. As it was impossible to know a priori whether Euphausiacea samples were infected, samples containing Euphausiacea and Anisakidae were artificially created first, to test and adapt the method originally created for fish. The method was then applied to real samples collected during a survey in the North Sea. This preliminary parasitological study provided some initial results on the levels of infection in Euphausiacea, in the North Sea.

Material and methods

Field sampling

Zooplankton samples were collected at 4 stations in the North Sea during the International Bottom Trawl Survey (IBTS) in January and February 2024 [2] (Fig. 1). IBTS is an international survey aimed at informing fisheries management by producing indices of abundance for the exploited fish species. Additionally, it provides an integrated picture of the ecosystem, by sampling all of its components, from plankton to top predators. Evaluation of commercial fish species also includes collection of their larvae. This sampling was used to opportunistically collect zooplankton occurring with fish larvae using a Methot Isaac Kidd (MIK) net. MIK is a 13 m long, 2 m diameter cylindroconic net, of 1.6 mm mesh size, narrowing to 500 μm in the final metre. The net was hauled at night during a double oblique tow from 5 m above the sea floor up to the surface, for at least 15 min. After hauling, samples were quickly processed on board, to discard all non-target organisms (e.g. gelatinous organisms and fish larvae) then stored frozen (−20 °C).

Thumbnail: Figure 1 Refer to the following caption and surrounding text. Figure 1

Map of stations sampled during the IBTS 2024 survey. The general location of sampling is presented in the bottom left inset, in which the rectangle highlights the more precise location of the sampling stations. Euphausiacea individuals from station C2043 (in italics) were used to calibrate the method, which was then applied to samples from the other three stations. Depth of the stations is illustrated by the colour of the dots.

Sample preparation

Back in the laboratory, and after defrosting, samples were carefully processed to select Euphausiacea individuals only, on the basis on specific morphological parameters [16]. Identification to the species level was not possible, as the conservation of Euphausiacea in sea water and thawing led to the detachment of the thoracic limbs, which are mandatory parts for identification.

Euphausiacea from station C2043 were used for method development because it was thought on board that this was the station where they were most abundant, and considering a priori that the repeated steps of method calibration would require a large number of samples. Then, Euphausiacea individuals from stations C2042, C2046, and C2051 were used for method application and the determination of infection levels. Euphausiacea individuals from all four stations were separated individually and frozen again.

The feasibility of the UV light detection method was tested by applying the same approach as for fish: all Euphausiacea individuals from stations C2042, C2046, and C2051 were observed under UV light (366 nm), following the standardised method ISO 23036-1:2021(F) [17]. Because of their thinness, the Euphausiacea were not crushed by the press. Once the method was calibrated, Euphausiacea were analysed individually with the developed method for the parasitological study. At each station, 198 Euphausiacea were analysed, leading to a total number of samples analysed of 594 over the three stations.

DNA extraction

Prior to the experiment per se, a subset of Euphausiacea and Anisakidae L3 larvae were individually weighed, to estimate the average mass ratio between one Anisakidae individual and one Euphausiacea individual. One Euphausiacea and half of an Anisakidae L3 larvae were weighed using a precision scale, (to the nearest 10−3 mg), considering that in a Euphausiacea, the Anisakidae should be at most half the size of a L3 larvae. Assessing these mean weights was necessary to set up artificial samples with a relative importance of Euphausiacea and Anisakidae consistent with natural mass ratios (Table 2).

Table 2

Mass of Anisakidae and Euphausiacea for each sample used for method calibration. Samples with 0% and 100% proportion of Anisakidae are negative and positive controls, respectively and are presented for information purposes only.

For method development, Euphausiacea from station C2043 and Anisakis individuals (L3 larvae) collected in the visceral cavity of hake Merluccius merluccius caught in the Celtic Sea (FAO Zone 27.VII.h) in a previous project were used. Euphausiacea and Anisakis were separately ground using a Mixer Mill MM 400 (Retsch GmbH, Haan, Germany) at 25 Hz for 2 min, with six 7 mm diameter stainless steel beads in a grinding bowl. Then, shredded Anisakidae and Euphausiacea were artificially mixed. The term “artificially” is used to refer to samples prepared by mixing Euphausiacea and Anisakis larvae in different proportions. The ratio of 1% was considered to be a representative value for the average mass proportion between Anisakidae and Euphausiacea. For the parasitological study, Euphausiacea from stations C2042, C2046, and C2051 were ground with two 3 mm diameter beads, following the same parameters as above.

Artificial samples of Euphausiacea with Anisakidae were used to compare the performance of four DNA extraction methods: proteinase K internal protocol (named proteinase K), a Promega Wizard Genomic DNA purification kit (Promega, Madison, WI, USA), a QIAGEN Genomic-tip kit (QIAGEN, Hilden, Germany), and a Nucleobond Column AXG kit (Macherey-Nagel, Düren, Germany).

Efficiency of the four extraction methods was tested using artificial Anisakidae/Euphausiacea samples of 1% and 0.5% proportions. Then, the proportion of Anisakidae in the artificial sample was decreased (0.25% and 0.1%) to test the selected method’s detection limits.

For the first method, samples were lysed using proteinase K at 20 mg mL−1 (QIAGEN) and 10× lysis buffer (60 mg L−1 Tween 20, 0.5 M EDTA and Trisbase) for 6 h at 55 °C, at 600 rpm in a Thermomixer C. Then, the proteinase K was inactivated for 10 min at 95 °C, at 300 rpm. After brief centrifugation, supernatant was collected and preserved. All three commercial kits were applied, following the manufacturer’s instructions. Two control solutions, one positive with only Anisakidae tissue, one negative with only Euphausiacea tissues, were also tested. DNA extractions were performed in triplicate.

PCR amplification and electrophoresis

For all samples (whether artificially created or originating from natural samples), a fragment of 530 bp of the mitochondrial cytochrome c oxidase subunit II (mtDNA Cox2) gene was amplified using the primers F-univ-nem (5′– GGT GTT CTT TCT TTT GTT TCT G –3′) and R-univ-nem (5′– ATA AAA CTA TGG TTA GCC CCA C –3′), according to Seesao et al. [46]. These primers are specific to mtDNA Cox2 from Anisakidae and Hysterothylacium, formerly belonging to the Anisakidae family but now included in the Raphidascarididae family [52]. PCR reactions were carried out in a total volume of 50 μL with 2 μL of extracted DNA and 48 μL of PCR mix consisting of 10× reaction buffer, 25 mM MgCl2, 10 mM dNTPs, 50 μM of each primer (Eurofins, Nantes, France) and 5 U/μL of HotStarTaq DNA polymerase (QIAGEN, Germany).

Amplifications were carried out in a 2720 Therma Cycler (Applied Biosystems, Foster City, CA, USA) with a preheating step at 95 °C for 15 min, 50 cycles of amplification (95 °C for 30 s, 58 °C for 30 s, 72 °C for 30 s) and a final extension step at 72 °C for 7 min. All PCR products were visualised on a 1% agarose gel with 3.5‰ of ethidium bromide for visualisation.

Result analyses and parasitological descriptors

Each DNA extract was tested in triplicate by PCR. The method was validated (i.e. considered reproducible) if, for each proportion, (i) at least two amplifications were positive, (ii) if controls displayed expected results, i.e. amplification for positive control and no amplification for negative control, and (iii) in the absence of false positive or negative. For the parasitological study, each extracted sample was analysed three times by PCR. It was considered positive if DNA amplification was observed after PCR on the electrophoresis gel for at least two out of three times for one sample, and if controls displayed expected results as mentioned above.

Euphausiacea infection was expressed as prevalence, i.e. as the numbers of individuals where Anisakidae were detected divided by the total number of individuals analysed [6].

In the parasitological study, in order to determine the parasite species observed in Euphausiacea, PCR products of positive samples and 530 bp were sequenced from both sides (forward and reverse) using Sanger sequencing (Eurofins). The obtained sequences were analysed by GenBank BLAST software (version 2023.2.1.0).

Results

Selection of the DNA extraction method

Using solutions with 0.5% and 1% Anisakidae/Euphausiacea mass ratios, reproductive results were obtained with a Wizard Genomic DNA purification kit only, without any false positives or negatives (Fig. 2). Results obtained with the other three methods were inconclusive, and provided inconsistent results, such as no effect of increased concentration.

Thumbnail: Figure 2 Refer to the following caption and surrounding text. Figure 2

Agarose gel electrophoresis of PCR products from DNA extracted from 0.5 and 1% Anisakidae/Euphausiacea mass ratio samples after (A) proteinase K protocol, (B) Wizard Genomic DNA purification kit, (C) QIAGEN Genomic-tip kit, and (D) Nucleobond Column AXG kit. Each condition (method/mass ratio) was extracted three times (numbered 1, 2, and 3). Each sample was amplified and deposited in triplicate. Water was used as the negative control (*).

Consequently, the Promega kit was used on artificial samples of Euphausiacea with proportions of 0.25% and 0.1% of Anisakidae for further validation. Again, each proportion was extracted three times and each sample was amplified three times. At least two replicates out of three for each proportion were positive.

Parasitological study

Application of DNA extraction with the Wizard Genomic DNA purification kit to Euphausiacea from stations C2042, C2046, and C2051 detected 30 positive samples out of 594 Euphausiacea individuals.

Sequencing of positive samples made it possible to identify 22 samples as Anisakis simplex, A. pegreffii, and Hysterothylacium aduncum (Table 3). Of note, A. simplex was the most abundant parasite, with a proportion of 82% among all parasites. Seven samples from station C2046 and one from station C2051 could not be identified with sequencing. Either sequencing failed or obtained sequences were not analysable. Thus, prevalences were calculated from the number of infected Euphausiacea confirmed by Sanger sequencing. They ranged from 2.5 to 4.5%.

Table 3

Anisakidae prevalence at the three stations C2042, C2046, and C2051 of the IBTS 2024 survey and results of species identification by sequencing. Number of positive Euphausiacea: number of Euphausiacea samples with positive results on PCR detection of Anisakidae. Parasite identification: results of species identification by Sanger sequencing.

The standard UV light method provided inconclusive results: fluorescence was observed on the cephalothorax of some Euphausiacea; PCR was negative for these samples. In addition, no fluorescence was observed for samples with positive PCR results.

Discussion

Method development

The first objective of this study was to develop and validate a fast and reliable method to detect Anisakidae larvae in zooplankton and more specifically in Euphausiacea. Importantly, former studies on the presence of Anisakidae in zooplanktonic crustaceans relied on visual observations. However, this method is time-consuming and requires working on fresh material since crustaceans rapidly become opaque after death, reducing the ability to see the Anisakidae larvae within the crustaceans.

Following former studies on the fluorescence of Anisakidae larvae in fish [20, 40], a similar approach was adopted. However, the results based on fluorescence and molecular analysis were inconsistent, both because the fluorescence observed in some individuals could not be correlated with the detection of Anisakidae, and because PCR-positive individuals did not fluoresce (data not shown). Nonetheless, even though inconclusive, these results must be published and discussed, to avoid repeating them, and to challenge the limitations of the UV light method, consistently with the importance of publishing negative results [4]. Fluorescence-based detection of Anisakidae relies on the presence of lipofuscin in the parasite, a molecule that accumulates in tissues with parasite’s age and becomes fluorescent under UV light [40]. It can, therefore, be hypothesised that lipofuscin is absent or present in low concentrations in L2 larvae, which prevents the observation of fluorescence at this stage, and accumulates afterwards in the cycle when the parasite enters the fish. On the contrary, the observation of fluorescence in individuals for which PCR results were negative can be explained by the herbivorous diet of Euphausiacea, and consumption of microalgae whose pigments fluoresce when exposed to low wavelengths [41]. A molecular approach was thus assessed for the fast and reliable detection of Anisakidae larvae in zooplanktonic crustacean.

The first step was to constitute infected material. Since there were no natural infected crustaceans available, artificial samples were constituted. To mimic the natural infection of Euphausiacea by Anisakidae, mean weights of both organisms were obtained. The weight ratio between Anisakidae larvae (extracted from fish) and Euphausiacea individuals was between 0.8 and 2% (data not shown). Thus, artificial mixtures of 0.5 and 1% Anisakidae weight/Euphausiacea weight were built, so as to reproduce the presence of one Anisakidae larvae in one Euphausiacea. However, since no data were available about the weight of Anisakidae larvae in the crustacean host, more stringent conditions were also assessed with 0.25 and 0.1% weight ratios. Additionally, only the Wizard Genomic DNA purification kit allowed for amplification of Anisakidae DNA, even with these stringent and limited conditions.

Moreover, complex matrices such as whole individuals might lead to the presence of PCR inhibitors or nonspecific PCR amplification. DNA extraction from whole Euphausiacea leads to a mixture of DNA coming from Euphausiacea, its stomach content, organisms living on its carapace and potential Anisakidae. Greenstone et al. [14] and O’Rorke et al. [36], while studying stomach contents of crustaceans, described trouble shooting of PCR amplification from complex mixtures and/or different organisms, whereas Itoïz et al. [19] described the presence of PCR inhibitors in clams depending on the tissue used for the extraction. To counter these drawbacks, confirmation steps were added to the present protocol. Like in non-specific amplifications, Sanger sequencing of PCR products made it possible to confirm the specificity of the analysis. Among 30 Euphausiacea positive for Anisakidae in PCR, the presence of Anisakidae or Raphidascaridae DNA was confirmed in 22 individuals. For the eight remaining individuals, either no sequencing or non-relevant sequences were obtained, reinforcing the hypothesis of non-specific amplification. This result confirms that, at this step of method development, and before further analyses of the potential mechanisms originating the non-specific amplification, sequencing of positive results is required to ensure robust calculation of prevalence data on the basis of Sanger sequencing results. Like for PCR inhibitors, supplementary steps were carried out during protocol optimisation and/or during the sample analyses. The Wizard Genomic purification kit allowed for repeatable and systematic amplification of all tested artificial samples, even when Anisakidae DNA was present at very low concentrations in the samples (0.1% weight ratio), validating the elimination of most PCR inhibitors during the DNA extraction step. Moreover, as a supplementary precaution, during the sample analyses, each Euphausiacea DNA was tested in triplicate. Each PCR product reporting a positive result was then sequenced to confirm the results. This method was optimised and validated on Euphausiacea. It yielded very conclusive and promising results. However, even though other zooplanktonic hosts of Anisakidae are mainly crustaceans, such as Mysidacea or Copepoda, a validation step will be mandatory before extending this method to study the distribution of Anisakidae in these taxonomic groups. Naturally, to have an overall view of the circulation of Anisakidae among the zooplanktonic community, every member will have to be considered and carefully assessed for its role in this parasite lifecycle.

Ecology of the parasite

Prevalence values observed here (2.5–4.5%) are markedly higher than those reported so far in the literature (0–4%, but usually less than 0.01%, Table 1). The higher sensitivity of the molecular method, when compared with visual observation, can be proposed as the first cause of these higher values. This would of course be consistent with the technical limitations of the visual method already mentioned and would be an advantage in favour of the PCR method. In addition, an increase in infection would also be in line with the general trend of an increase of Anisakidae infection, related to several factors, such as global warming, increases in cetaceans populations, and an increase in cetacean infection by Anisakidae [10, 33]. An increased abundance of the harbour porpoise (Phocoena phocoena), the most abundant cetacean species in the North Sea, has notably been observed since the early 1990s in the south of the North Sea [5]. Harbour porpoise is known to host both Anisakis and Hysterothylacium genera [26]. A recent (2009–2015) increase in Anisakis infection was also observed for wild salmon in Scotland [22]. Simultaneously, global warming drove a decrease in the abundance of Euphausiacea in the entire North Atlantic and also in the area studied here [8]. Euphausiacea remained, nonetheless, important in the western part of the North Sea, where it locally drives the functioning of the planktonic ecosystem. It is notably an indicator taxon of this area, and explains the difference in plankton assemblages between the west of the North Sea and the other areas (A. Delegrange unpubl. results; [35]). These simultaneous trends for both first and final hosts of Anisakidae are a reasonable explanation for increased prevalence. If the number of Anisakidae eggs emitted in the water is on the rise, as the number of cetaceans is increasing, while the number of Euphausiacea is going down, the probability of an Anisakidae larva encountering an Euphausiacea increases mathematically. Similarly, Anisakidae eggs hatch faster in warmer waters. In the context of global climate change, this also increases infection probabilities. All of these factors seem to confirm that an increase in Anisakidae infection in Euphausiacea is plausible, and that the pattern detected by the method is consistent with ecological patterns.

The species diversity of nematodes identified with sequencing is also consistent with the pattern observed for other members of the cycle, and may illustrate global climate change effects. Euphausiacea were mostly infected with A. simplex and, to a lesser extent, with H. aduncum, a pattern already described by Smith [50]. Regarding fish, A. simplex is commonly observed as the dominant species of Anisakidae in this area [1, 7, 11]. It dominates over H. aduncum in several fish species in the area, e.g. grey gurnard Eutrigla gurnardus [21], beaked redfish Sebastes mentella [23], Atlantic salmon Salmo salar [22], or cod Gadus morhua [39]. Lastly, A. simplex is the main Anisakidae species infecting P. phocoena in the North Sea [26]. The observation of one A. pegreffii sequence can be considered unexpected, as this species is considered to be restricted to the western Mediterranean Sea and the southern coast of the North Atlantic, i.e. off Spain and Portugal [34]. Nonetheless, some recent occurrences of A. pegreffii in cod [11] or Atlantic mackerel Scomber scombrus [27] in the North Sea, along with the observation of an A. simplex × A. pegreffii hybrid in hake Merluccius merluccius off Southwest Ireland [1] may challenge our understanding of the current geographical distribution of A. pegreffii. In these studies, low prevalence (1–9 infected fish) is interpreted as the transport of A. pegreffii away from its usual area of presence through host feeding migrations. However, increasing observations might illustrate a shift in the geographical distribution of this species. A. simplex and A. pegreffii have different thermic preferenda, as A. simplex eggs and L1 larvae are more adapted to colder waters than A. pegreffii’s [12]. Warming of North Sea waters might favour the achievement of A. pegreffii cycle, and allow some eggs of this species to hatch in these waters, a phenomenon that is likely to intensify as warming continues. More repeated analyses on all members of the life cycle would be needed to ascertain a shift northward in A. pegreffii distribution, and potentially relate it to warming waters. The present method would be appropriate to do so in zooplanktonic hosts.

Eventually, this preliminary sampling will provide an opportunity to continue exploring the impact of depth on Anisakidae infection. Anisakidae are well-known to be pelagic parasites. Their entire cycle implies pelagic organisms, and they have developed adaptation to the pelagic realm, like cuticles that increase the buoyancy of the second stage larvae [25]. Previous studies on fish have demonstrated higher infection at the deepest stations for individuals of the same species [7, 18, 44] or, when comparing certain species, for those with a higher ability to forage in the pelagic realm [29]. Similarly, for Euphausiacea, Smith [50] confirmed a depth effect, as he observed no infected Euphausiacea from shallow stations, but infected ones in 10 out of the 31 deep stations (>100 m depth), and with a prevalence of 0.54% (15 individuals infected out of 2 687 observed). Results obtained here are consistent with this pattern, as the prevalence is the highest for the deepest station. The work performed here was not dedicated to an actual analysis of real data, explaining why so few stations were used. However, the data appear promising and will need to be confirmed by applying this method to a larger sample. For example, spatial variation of Euphausiacea infection in the English Channel would be needed to better support the role of Euphausiacea in the infection pattern observed for horse mackerel [7]. Similarly, investigating seasonal variations in Euphausiacea infection would be required to improve our understanding of seasonal patterns in fish.

Acknowledgments

We are grateful to the crew of the R/V Thalassa for their help during sample collection as part of the IBTS 2024 cruise. We would also like to thank the French Parasitological Society for allocating a master’s grant to Anais Delhomme that allowed this work to be performed. This work benefited from funding from the Hauts-de-France Region, through the STIMULE project ANIMODE “Circulation of Anisakidae in marine food webs”.

Conflicts of interest

The authors declare no conflicts of interest.

References

  1. Atroch F, Rangel LF, Ramos P, Ayra-Pardo C, Santos MJ. 2025. Preliminary observations on muscle distribution of Anisakis sp. L3 in European hake off Southwest Ireland. Food and Waterborne Parasitology, 41, e00303. [Google Scholar]
  2. Auber A, Caboche J, Bled Defruit G. 2024. IBTS 2024 cruise, RV Thalassa. https://doi.org/10.17600/18002960. [Google Scholar]
  3. Bao M, Giulietti L, Levsen A, Karlsbakk E. 2023. Resurrection of genus Phocanema Myers, 1959, as a genus independent from Pseudoterranova Mozgovoĭ, 1953, for nematode species (Anisakidae) parasitic in pinnipeds and cetaceans, respectively. Parasitology International, 97, 102794. [CrossRef] [PubMed] [Google Scholar]
  4. Bik EM. 2024. Publishing negative results is good for science. Access Microbiology, 6(4), 000792. [Google Scholar]
  5. Bouveroux T, Waggitt JJ, Belhadjer A, Cazenave PW, Evans PGH, Kiszka JJ. 2020. Modelling fine-scale distribution and relative abundance of harbour porpoises in the Southern Bight of the North Sea using platform-of-opportunity data. Journal of the Marine Biological Association of the United Kingdom, 100, 481–489. [Google Scholar]
  6. Bush AO, Lafferty KD, Lotz JM, Shostak AW. 1997. Parasitology meets ecology on its own terms: Margolis et al. revisited. Journal of Parasitology, 83, 575–583. [CrossRef] [Google Scholar]
  7. Cresson P, Bourgau O, Cordier R, Couvreur C, Rouquette M, Gay M. 2023. Fish length, diet, and depth drive Anisakis levels in a zooplankton-feeding fish. Canadian Journal of Fisheries and Aquatic Sciences, 80, 1495–1508. [Google Scholar]
  8. Edwards M, Hélaouët P, Goberville E, Lindley A, Tarling GA, Burrows MT, Atkinson A. 2021. North Atlantic warming over six decades drives decreases in krill abundance with no associated range shift. Communications Biology, 4, 644. [Google Scholar]
  9. European Food Safety Authority (EFSA) Panel on Biological Hazards (BIOHAZ) 2010. Scientific opinion on risk assessment of parasites in fishery products. EFSA Journal, 8, 1543. [Google Scholar]
  10. Fiorenza EA, Wendt CA, Dobkowski KA, King TL, Pappaionou M, Rabinowitz P, Samhouri JF, Wood CL. 2020. It’s a wormy world: Meta-analysis reveals several decades of change in the global abundance of the parasitic nematodes Anisakis spp. and Pseudoterranova spp. in marine fishes and invertebrates. Global Change Biology, 26, 2854–2866. [Google Scholar]
  11. Gay M, Bao M, MacKenzie K, Pascual S, Buchmann K, Bourgau O, Couvreur C, Mattiucci S, Paoletti M, Hastie LC, Levsen A, Pierce GJ. 2018. Infection levels and species diversity of ascaridoid nematodes in Atlantic cod, Gadus morhua, are correlated with geographic area and fish size. Fisheries Research, 202, 90–102. [Google Scholar]
  12. Gomes TL, Quiazon KM, Itoh N, Fujise Y, Yoshinaga T. 2023. Effects of temperature on eggs and larvae of Anisakis simplex sensu stricto and Anisakis pegreffii (Nematoda: Anisakidae) and its possible role on their geographic distributions. Parasitology International, 92, 102684. [Google Scholar]
  13. Gómez-Gutiérrez J, Robinson CJ, Kawaguchi S, Nicol S. 2010. Parasite diversity of Nyctiphanes simplex and Nematoscelis difficilis (Crustacea: Euphausiacea) along the northwestern coast of Mexico. Diseases of Aquatic Organisms, 88, 249–266. [Google Scholar]
  14. Greenstone MH, Weber DC, Coudron TA, Payton ME, Hu JS. 2012. Removing external DNA contamination from arthropod predators destined for molecular gut-content analysis. Molecular Ecology Resources, 12, 464–469. [Google Scholar]
  15. Gregori M, Roura Á, Abollo E, González ÁF, Pascual S. 2015. Anisakis simplex complex (Nematoda: Anisakidae) in zooplankton communities from temperate NE Atlantic waters. Journal of Natural History, 49, 755–773. [Google Scholar]
  16. Guglielmo L, Granata A, Guglielmo R. 2015. Order Euphausiacea. Revista IDE@-SEA B, 86, 1–20. [Google Scholar]
  17. International Organization for Standardization. 2021. ISO 23036-1:2021(F). [Google Scholar]
  18. Irigoitia MM, Palomba M, Braicovich PE, Gutiérrez MP, Canel D, Mattiucci S, Timi JT. 2025. Anisakid nematodes in skates (Chondrichthyes: Rajiformes) along a depth gradient: Do bathymetry and host size influence the presence of zoonotic parasites? International Journal for Parasitology, 56, 104747. [Google Scholar]
  19. Itoïz S, Perennou M, Mouronvalle C, Derelle E, Le Goïc N, Bidault A, De Montaudouin X, Arzul I, Soudant P, Chambouvet A. 2021. Development of duplex TaqMan-based real-time PCR assay for the simultaneous detection of Perkinsus olseni and P. chesapeaki in host Manila clam tissue samples. Journal of Invertebrate Pathology, 184, 107603. [Google Scholar]
  20. Karl H, Leinemann M. 1993. A fast and quantitative detection method for nematodes in fish fillets and fishery products. Archiv für Lebensmittelhygiene, 44, 124–125. [Google Scholar]
  21. Karl H, Levsen A. 2011. Occurrence and distribution of anisakid nematodes in Grey gurnard (Eutrigla gurnardus L.) from the North Sea. Food Control, 22, 1634–1638. [Google Scholar]
  22. Kent AJ, Pert CC, Briers RA, Diele K, Rueckert S. 2020. Increasing intensities of Anisakis simplex third-stage larvae (L3) in Atlantic salmon of coastal waters of Scotland. Parasites & Vectors, 13, 62. [Google Scholar]
  23. Klapper R, Kuhn T, Münster J, Levsen A, Karl H, Klimpel S. 2015. Anisakid nematodes in beaked redfish (Sebastes mentella) from three fishing grounds in the North Atlantic, with special notes on distribution in the fish musculature. Veterinary Parasitology, 207, 72–80. [Google Scholar]
  24. Klimpel S, Palm HW, Rückert S, Piatkowski U. 2004. The life cycle of Anisakis simplex in the Norwegian Deep (northern North Sea). Parasitology Research, 94, 1–9. [Google Scholar]
  25. Køie M, Berland B, Burt MDB. 1995. Development to third-stage larvae occurs in the eggs of Anisakis simplex and Pseudotetranova decipiens (Nematoda, Ascaridoidea, Anisakidae). Canadian Journal of Fisheries and Aquatic Sciences, 52, 134–139. [Google Scholar]
  26. Lakemeyer J, Siebert U, Abdulmawjood A, Ryeng KA, IJsseldijk LL, Lehnert K. 2020. Anisakid nematode species identification in harbour porpoises (Phocoena phocoena) from the North Sea, Baltic Sea and North Atlantic using RFLP analysis. International Journal for Parasitology: Parasites and Wildlife, 12, 93–98. [Google Scholar]
  27. Levsen A, Cipriani P, Mattiucci S, Gay M, Hastie LC, MacKenzie K, Pierce GJ, Svanevik CS, Højgaard DP, Nascetti G, González AF, Pascual S. 2018. Anisakis species composition and infection characteristics in Atlantic mackerel, Scomber scombrus, from major European fishing grounds — reflecting changing fish host distribution and migration pattern. Fisheries Research, 202, 112–121. [Google Scholar]
  28. Levsen A, Svanevik CS, Cipriani P, Mattiucci S, Gay M, Hastie LC, Bušelić I, Mladineo I, Karl H, Ostermeyer U, Buchmann K, Højgaard DP, González ÁF, Pascual S, Pierce GJ. 2018. A survey of zoonotic nematodes of commercial key fish species from major European fishing grounds—Introducing the FP7 PARASITE exposure assessment study. Fisheries Research, 202, 4–21. [CrossRef] [Google Scholar]
  29. Lile NK. 1998. Alimentary tract helminths of four pleuronectid flatfish in relation to host phylogeny and ecology. Journal of Fish Biology, 53, 945–953. [Google Scholar]
  30. Llarena-Reino M, Piñeiro C, Antonio J, Outeriño L, Vello C, González ÁF, Pascual S. 2013. Optimization of the pepsin digestion method for anisakids inspection in the fishing industry. Veterinary Parasitology, 191, 276–283. [Google Scholar]
  31. Makings P. 1981. Mesopodopsis slabberi (Mysidacea) at Millport, W. Scotland, with the parasitic nematode Anisakis simplex. Crustaceana, 41, 310–312. [Google Scholar]
  32. Marcogliese DJ. 1995. The role of zooplankton in the transmission of helminth parasites to fish. Reviews in Fish Biology and Fisheries, 5, 336–371. [CrossRef] [Google Scholar]
  33. Mastick NC, Fiorenza E, Wood CL. 2024. Meta-analysis suggests that, for marine mammals, the risk of parasitism by anisakids changed between 1978 and 2015. Ecosphere, 15, e4781. [Google Scholar]
  34. Mattiucci S, Farina V, Campbell N, MacKenzie K, Ramos P, Pinto AL, Abaunza P, Nascetti G. 2008. Anisakis spp. larvae (Nematoda: Anisakidae) from Atlantic horse mackerel: Their genetic identification and use as biological tags for host stock characterization. Fisheries Research, 89, 146–151. [Google Scholar]
  35. Neven CJ, Giraldo C, Girardin R, Lefebvre A, Lefebvre S, Loots C, Meunier CL, Marchal P. 2024. Winter distribution of zooplankton and ichthyoplankton assemblages in the North Sea and the English Channel. PLoS ONE, 19, e0308803. [Google Scholar]
  36. O’Rorke R, Jeffs AG, Fitzgibbon Q, Chow S, Lavery S. 2013. Extracting DNA from whole organism homogenates and the risk of false positives in PCR based diet studies: A case study using spiny lobster larvae. Journal of Experimental Marine Biology and Ecology, 441, 1–6. [Google Scholar]
  37. Oshima T, Shimazu T, Koyama H, Akahane H. 1969. On the larvae of the genus Anisakis (Nematoda: Anisakinae) from euphausiids. Japanese Journal of Parasitology, 18, 241–248. [Google Scholar]
  38. Paoletti M, Mattiucci S, Colantoni A, Levsen A, Gay M, Nascetti G. 2018. Species-specific real time-PCR primers/probe systems to identify fish parasites of the genera Anisakis, Pseudoterranova and Hysterothylacium (Nematoda: Ascaridoidea). Fisheries Research, 202, 38–48. [Google Scholar]
  39. Perdiguero-Alonso D, Montero FE, Raga JA, Kostadinova A. 2008. Composition and structure of the parasite faunas of cod, Gadus morhua L. (Teleostei: Gadidae), in the North East Atlantic. Parasites & Vectors, 1, 23. [Google Scholar]
  40. Pippy JHC. 1970. Use of ultraviolet light to find parasitic nematodes in situ. Journal of the Fisheries Research Board of Canada, 27, 963–965. [Google Scholar]
  41. Poryvkina L, Babichenko S, Leeben A. 2000. Analysis of phytoplankton pigments by excitation spectra of fluorescence. EARSeL-SIG-Workshop LIDAR. Institute of Ecology/LDI: Tallinn, Estonia. pp. 224–232. [Google Scholar]
  42. Poulin R. 2025. Crustaceans as intermediate hosts in the life cycle of parasites, in Ecology and Conservation, Gutow L, Thiel M, Editors. Oxford University Press. p. 126–145. [Google Scholar]
  43. Poulin R, Blasco-Costa I, Randhawa HS. 2016. Integrating parasitology and marine ecology: Seven challenges towards greater synergy. Journal of Sea Research, 113, 3–10. [Google Scholar]
  44. Pulleiro-Potel L, Barcala E, Mayo-Hernández E, Muñoz P. 2015. Survey of anisakids in commercial teleosts from the western Mediterranean Sea: Infection rates and possible effects of environmental and ecological factors. Food Control, 55, 12–17. [Google Scholar]
  45. Rahmati AR, Kiani B, Afshari A, Moghaddas E, Williams M, Shamsi S. 2020. World-wide prevalence of Anisakis larvae in fish and its relationship to human allergic anisakiasis: a systematic review. Parasitology Research, 119, 3585–3594. [Google Scholar]
  46. Seesao Y, Audebert C, Verrez-Bagnis V, Merlin S, Jérôme M, Viscogliosi E, Dei-Cas E, Aliouat-Denis CM, Gay M. 2014. Monitoring of four DNA extraction methods upstream of high-throughput sequencing of Anisakidae nematodes. Journal of Microbiological Methods, 102, 69–72. [Google Scholar]
  47. Shamsi S. 2021. The occurrence of Anisakis spp. in Australian waters: past, present, and future trends. Parasitology Research, 120, 3007–3033. [Google Scholar]
  48. Smith JW. 1971. Thysanoessa inermis and T. longicaudata (Euphausiidae) as first intermediate hosts of Anisakis sp. (Nematoda: Ascaridata) in the northern North Sea, to the North of Scotland and at Faroe. Nature, 234, 478–478. [Google Scholar]
  49. Smith JW. 1983. Anisakis simplex (Rudolphi, 1809, det. Krabbe, 1878)(Nematoda: Ascaridoidea): morphology and morphometry of larvae from euphausiids and fish, and a review of the life-history and ecology. Journal of Helminthology, 57, 205–224. [Google Scholar]
  50. Smith JW. 1983. Larval Anisakis simplex (Rudolphi, 1809, det. Krabbe, 1878) and larval Hysterothylacium sp. (Nematoda: Ascaridoidea) in euphausiids (Crustacea: Malacostraca) in the North-East Atlantic and northern North Sea. Journal of Helminthology, 57, 167–177. [Google Scholar]
  51. Smith JW, Snyder JM. 2005. New locality records for third-stage larvae of Anisakis simplex (sensu lato) (Nematoda: Ascaridoidea) in euphausiids Euphausia pacifica and Thysanoessa raschii from Prince William Sound, Alaska. Parasitology Research, 97, 539–542. [Google Scholar]
  52. WoRMS Editorial Board. 2025. Raphidascarididae Hartwich, 1954. Available from: https://marinespecies.org/aphia.php?p=taxdetails&id=731849. [Google Scholar]

Cite this article as: Delhomme A, Gay M, Delegrange A, Bourgau O, Duflot M & Cresson P. 2026. A fast and reliable molecular method to detect Anisakidae in zooplanktonic Euphausiacea. Parasite 33, 39. https://doi.org/10.1051/parasite/2026038.

All Tables

Table 1

Existing data on the infection of zooplanktonic crustaceans by Anisakidae larvae. nd: no data provided in the original article.

Table 2

Mass of Anisakidae and Euphausiacea for each sample used for method calibration. Samples with 0% and 100% proportion of Anisakidae are negative and positive controls, respectively and are presented for information purposes only.

Table 3

Anisakidae prevalence at the three stations C2042, C2046, and C2051 of the IBTS 2024 survey and results of species identification by sequencing. Number of positive Euphausiacea: number of Euphausiacea samples with positive results on PCR detection of Anisakidae. Parasite identification: results of species identification by Sanger sequencing.

All Figures

Thumbnail: Figure 1 Refer to the following caption and surrounding text. Figure 1

Map of stations sampled during the IBTS 2024 survey. The general location of sampling is presented in the bottom left inset, in which the rectangle highlights the more precise location of the sampling stations. Euphausiacea individuals from station C2043 (in italics) were used to calibrate the method, which was then applied to samples from the other three stations. Depth of the stations is illustrated by the colour of the dots.

In the text
Thumbnail: Figure 2 Refer to the following caption and surrounding text. Figure 2

Agarose gel electrophoresis of PCR products from DNA extracted from 0.5 and 1% Anisakidae/Euphausiacea mass ratio samples after (A) proteinase K protocol, (B) Wizard Genomic DNA purification kit, (C) QIAGEN Genomic-tip kit, and (D) Nucleobond Column AXG kit. Each condition (method/mass ratio) was extracted three times (numbered 1, 2, and 3). Each sample was amplified and deposited in triplicate. Water was used as the negative control (*).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.