Three hundred seventy-one
Twenty-two percent of the fish samples were parasitized with Anisakidae and Raphidascarididae.
Phylogenetic analysis using ITS2 and 18S gene markers demonstrated the species identification of
Fishes are important food sources that provide nutrition as well as livelihood to millions of people worldwide. The Food and Agriculture Organization (FAO) of the United Nations has documented that 88% of the global fish production in 2016 was intended for human consumption, surpassing that of the meat from terrestrial animals (FAO, 2018). Moreover, the fisheries sector is recognized as a large part of the trade industry for food commodities, with high export rates from developing countries such as the Philippines (FAO, 2020).
Filipinos are generally fish eaters because of the rich aquatic resources throughout the archipelago. One of the staple sources of proteins are the round scad fish of the family Carangidae (Rafinesque, 1815), genus
According to the studies from American, European, and Mediterranean countries, more than 100,000 species of fish parasites have been estimated to be present globally (Mishra and Chubb, 2009; Gibson et al., 2014; Costa et al., 2017). Parasites of marine pelagic fishes are very diverse and serve a very important role in the ecological perspective (Poulin and Morand, 2004). More importantly, fish parasitological studies shed light on zoonotic parasites that cause potent problems to human health. In the Philippines, studies on marine parasitology are scarce. The only published checklist of Philippine fish parasites by far was made available by Arthur and Lumanlan-Mayo (1997). Velasquez (1972) reported infection in marine fishes sold in Manila and nearby provinces. Furthermore, a survey on the parasites of several fish species from the Visayas region in the early 1990s was conducted. Parasites identified were microsporidia, cestodes, and nematodes, and all of which are found in low prevalence, posing low risks of infection. Reports of
This study aimed to determine the prevalence of ascaridoid nematode infection in
Ethical approval for animal use was obtained from the Institutional Animal Care and Use Committee (IACUC) University of the Philippines, Los Baños (Protocol No: CAS-2019-001). Approval and sampling were coordinated to the National Stock Assessment Program (NSAP) of the Bureau of Fisheries and Aquatic Resources Region IV-A (BFAR-4A) in Lucena City Regional Office.
Balayan Bay (13.8447° N, 120.8499° E) and Tayabas Bay (13.5999° N, 121.7195° E) are among the productive fishing grounds of the Philippines. These bays are adjacent to the West Philippine Sea and are part of the Verde Island Passage, which is considered to be rich in marine biodiversity. These two bays are included in Fisheries Management Area 12 (FMA 12) designated by the Bureau of Fisheries and Aquatic Resources (BFAR) for conservation and management of fisheries in the Philippines.
Fish samples were obtained from two fish landing areas in Batangas and Lucena in May 2018. The Batangas fish landing area located in Lemery and Calatagan has most of their fish catch from the adjacent Balayan Bay. The fish landing in Dalahican Port in Lucena acquires most of the fishes and other marine resources from Tayabas Bay (Fig. 1).
Figure 1
Fish landing area at Balayan Bay and Tayabas Bays, southern Luzon, Philippines, for fish sample collection.

A total of 371 fish samples were collected from Tayabas Bay and Balayan Bay and were identified based on key morphological characteristics found in FishBase (Froese and Pauly, 2019). Three species of
The fish samples were maintained in cold temperature and immediately transported to the Parasitology Research Laboratory, Institute of Biological Sciences, University of the Philippines Los Baños (UPLB) for dissection and processing. Before dissection, the total body length and weight of individual fish samples were measured. Sex of fish samples was also recorded based on observation of the gonads following dissection (Ohshimo et al., 2014).
The body cavity, stomach, and intestine of each dissected fish were carefully examined with the naked eye and under a stereomicroscope (Leica EZ4; Leica Microsystems Pte Ltd, Singapore) for the presence of helminth larvae (third-stage larvae of nematodes). The helminths were removed, counted, and preserved in 70% ethanol in separate Eppendorf tubes. These larvae were then cleared in a mixture of 5% glycerin in 70% ethanol for further examination under a compound microscope (Nikon, Tokyo, Japan). Parasites were identified based on morphological characteristics (Mattiucci et al., 2005; Borges et al., 2012; Hien et al., 2021). The specimens were photographed using digital microscope cameras OptixCam Summit Series SK2-5.2Xx USB (Microscope LLC, VA, USA) and Optika version 2.1 (Optika Microscopes, Bergamo, Italy). All parasites were removed and preserved in 70% ethanol, but only representative individuals of those initially identified as marine ascaridoid parasites through light microscopy were subjected to molecular processing for species identification.
Fifteen (five per fish species) morphologically similar ascaridoid larvae were randomly selected for molecular identification. Genomic DNA was extracted from ethanol-fixed midsections of individual larvae using a DNeasyTM Tissue Kit (Qiagen, Hilden, Germany). Identification of the parasites was carried out through analysis of the small subunit ribosomal ribonucleic acid sequences. Primers were designed using the Primer-BLAST (Primer3web v. 4.1.0) tool (Ye et al., 2012) to target the 18S rRNA gene and the internal transcribed spacer 2 (ITS2) region, flanked by partial sequences of 5.8S and 28S genes. The 18S gene marker was amplified using the primer pair SSU-F (5'-ATGAGAGGGCAAGTCTGGTG-3') and SSU-R (5'-CTGTCAATCCTCACGGTGTC-3'). The amplification reactions for this gene were performed in a 25 ml volume mastermix containing 2.5 ml 10X PCR buffer, 1.875 ml 10 mM dNTPs, 2 ml 25mM MgCl2, 1 ml 10 mM each of forward and reverse primers, 0.125 U Taq polymerase, and 2 ml 40ng/ml DNA template. The PCR conditions included an initial denaturation at 94°C for 4 min, followed by 40 cycles of denaturation at 94°C for 1 min, annealing at 55°C for 30 s, extension at 65°C for 1 min, and a final extension step at 72°C for 10 min. The 5.8S-ITS2-28S region was amplified using forward primer AscF (5'-GATCGATGAAGAACGCAGCC-3') and reverse primer AscR (5'-TTTGCAACTTTCCCTCACGG-3') with the following PCR conditions: initial denaturation of 94°C for 4 min, followed by 40 cycles of denaturation at 94°C for 1 min, annealing at 60°C for 30 s, and extension at 65°C for 1 min. A final extension step was added at 72°C for 10 min. The PCR was performed in a 25 ml volume with a similar composition of the PCR mastermix for the 18S gene amplification. The PCR products were resolved using a 1% agarose gel stained with 1% EtBr run in a horizontal gel electrophoresis system at 100 V for 30 min. Bands were visualized using a UV transilluminator. Distinct bands were sliced out from the gel and were stored in a properly labeled 2-ml microcentrifuge tube ready for gel extraction.
From the PCR amplification, 14 amplicons with a molecular weight of about 1 kb were purified using a QIAquick Gel Extraction Kit (Qiagen, Germany) and were sent to Macrogen, South Korea, for standard Sanger sequencing. Sequence assembly was performed using the Staden Package 4.0 (Staden et al., 2000). The generated consensus sequences were subjected to local alignment using the Basic Local Alignment Search Tool (
GenBank sequences of close species matches used for phylogeny construction.
Ascaridoid species | Host | Habitat | Geographical Origin | ITS2 | 18S rRNA |
---|---|---|---|---|---|
Marine | Tunisia | MT820022 | - | ||
Marine | China | - | MF072697 | ||
Marine | Italy | - | EF180082 | ||
Marine | Norway | AY826723 | - | ||
Marine | Japan | LC621351 | - | ||
Marine | China | - | MF072711 | ||
Marine | Egypt | HF911524 | - | ||
Bottlenose dolphin | Marine | China | KF673776 | - | |
Marine | Thailand | AB432909 | - | ||
Marine | Brazil | AY826724 | - | ||
Marine | Indonesia | KC928262 | - | ||
Marine | Brazil | EU327686 | - | ||
Marine | Portugal | JN005767 | - | ||
Marine | Mauritius | EU718473 | - | ||
Marine | South Africa | AY826725 | - | ||
Marine | Australia | - | EF180072 | ||
NM | NM | NM | - | AY702702 | |
NM | NM | NM | - | U94370 | |
Marine | Denmark | JX845137 | - | ||
Dolphin fish | Marine | South Korea | HQ702733 | - | |
Marine | China | - | MF072693 | ||
Marine | Brazil | - | JF718550 | ||
NM | NM | NM | - | U94374 | |
NM | NM | NM | - | U94376 | |
Marine | China | - | MF072705 | ||
Marine | China | - | MF072702 | ||
Marine | China | - | MF072707 | ||
Marine | Denmark | KM273088 | - | ||
NM | NM | NM | - | U94380 | |
Freshwater | Brazil | - | KP726276 | ||
Freshwater | Brazil | - | KP726278 | ||
Freshwater; brackish water | Finland | - | DQ503460 | ||
Freshwater | Brazil | - | KX859077 | ||
NM | Marine | South China Sea | JN102362 | - | |
Marine | China | KP326546 | MF072704 | ||
NM | NM | China | MF422212 | - | |
Marine | China | MH211584 | MF072692 | ||
NM | Marine | China | JF809816 | - | |
NM | NM | NM | - | U94381 | |
a | Freshwater | Brazil | MK141032 | - | |
a | NM | NM | NM | U94366 |
aOutgroup species
ITS2, internal transcribed spacer 2; NM, not mentioned.
The prevalence of parasite infection was calculated by dividing the number of infected fish by the total number of fish examined, multiplied by 100, and expressed as percentage. Mean intensity was calculated by dividing the total number of parasites by the total number of infected fish hosts examined. On the other hand, chi-square tests were performed to determine variations in helminth prevalence and mean intensities among fish host species and between sampling sites. This was performed using IBM SPSS Statistics v. 25 with 95% confidence interval and 0.05 level of significance. Parasite aggregation in the fish host population was measured using Poulin’s discrepancy index (
Three hundred seventy-one
Figure 2
Light microscopic images of representative marine ascaridoid larvae collected from

Figure 3
Prevalence (%) and mean intensity (worms/host ± SD) of marine ascaridoid larvae across the three

The variation in the prevalence and mean intensity of marine ascaridoid larvae with respect to the sampling areas is shown in Figure 4. A total of 185 individuals were collected from Tayabas Bay fish landing area in Quezon Province, while 186 individuals were from Balayan Bay in Batangas. Of the 371 fish samples, 51 (27.57%) were observed as infected from Tayabas Bay and 29 (15.59%) from Balayan Bay. While a higher prevalence is observed in Tayabas Bay than in Balayan Bay, there was a difference in the mean intensity of infection of fish samples collected from the two sites, 3 ± 5.33 in Tayabas Bay and 2 ± 1.49 in Balayan Bay. Statistical analysis revealed that prevalence (x2 [5] = 6.0,
Figure 4
Prevalence (%) and mean intensity (worms/host ± SD) of marine ascaridoid larvae across the three

The marine ascaridoid parasites showed an aggregated distribution in
Figure 5
Aggregated distribution of marine ascaridoid larval parasites from (A) pooled

Evaluation of the identity of the marine ascaridoid larvae was performed using molecular analysis based on the 18S rRNA and the 5.8S-ITS2-28S markers. Twelve of the 15 selected larval samples generated DNA sequences that were subjected to subsequent molecular analyses. Single-gene alignment of both 18S rRNA and 5.8S-ITS2-28S regions identified two genera,
BLAST results for the 18S rRNA and 5.8S-ITS2-28S sequences of ascaridoid samples.
Specimen | 18S rRNA | 5.8S-ITS2-28S | Family | ||
---|---|---|---|---|---|
accession (18S/ITS2) | GenBank accession | % Match | GenBank accession | % Match | |
OK659774/ OK659786 | MF072704 ( | 99.84 | MF422212 ( | 99.18 | Raphidascarididae |
OK659775/ OK659787 | MF072711 ( | 99.52 | HF911524 ( | 99.88 | Anisakidae |
OK659776/ OK659788 | MF072704 ( | 99.84 | MF422212 ( | 99.79 | Raphidascarididae |
OK659777/ OK659789 | MF072711 ( | 99.68 | HF911524 ( | 99.88 | Anisakidae |
OK659778/ OK659790 | MF072704 ( | 99.52 | MF422212 ( | 99.79 | Raphidascarididae |
OK659779/ OK659791 | MF072704 ( | 99.52 | MF422212 ( | 99.59 | Raphidascarididae |
OK659780/ OK659792 | MF072711 ( | 99.20 | HF911524 ( | 98.66 | Anisakidae |
OK659781/ OK659793 | MF072711 ( | 99.68 | HF911524 ( | 98.79 | Anisakidae |
OK659782/ OK659794 | MF072711( | 99.52 | HF911524 ( | 99.87 | Anisakidae |
OK659783/ OK659795 | MF072711( | 99.68 | HF911524 ( | 98.86 | Anisakidae |
OK659784/ OK659796 | MF072711 ( | 99.68 | HF911524 ( | 98.61 | Anisakidae |
OK659785/ OK659797 | MF072711 ( | 99.68 | HF911524 ( | 99.87 | Anisakidae |
ITS2, internal transcribed spacer 2.
Pairwise distance comparison (Table 3) indicates low genetic variation within species for both 18S and 5.8S-ITS2-28S markers. However, sequence variation between species and between genera is at lowest for 18S rRNA gene, with
Figure 6
Intra- and interspecific sequence variations in Anisakidae and Raphidascarididae obtained in the GenBank and sample isolate sequences. (A) 18S rRNA and (B) 5.8S-ITS2-28S regions.

Phylogenetic reconstruction using the ML and NJ methods based on the 18S rRNA gene showed two major clusters of the currently studied ascaridoid larvae. Four samples were resolved in close association with
Figure 7
ML phylogenetic tree of Anisakidae and Raphidascarididae inferred from 628 nucleotide positions of the 18S rRNA gene using the TrNef + G model of DNA substitution and rooted on

Corrected (TrNef + G) genetic distances of 18S rRNA gene and corrected (HKY + G) genetic distances of the 5.8S-ITS2-28S region of Anisakidae and Raphidascarididae nematodes.
Ascaridoid comparisons | Number of pairwise comparisons | 18S rRNA | 5.8S-ITS2-28S region | ||
---|---|---|---|---|---|
Minimum | Maximum | Minimum | Maximum | ||
67 (18S), 118 (ITS2) | 0 | 0.0182098 | 0 | 0.0077422 | |
494 (18S), 410 (ITS2) | 0 | 0.04259699 | 0.00493606 | 1.8776778 |
ITS2, internal transcribed spacer 2.
The phylogeny based on the partial sequence of the 5.8S-ITS2-28S rRNA region (Fig. 8) showed similar grouping of the four
Figure 8
ML phylogenetic tree of Anisakidae and Raphidascarididae inferred from 460 nucleotide positions from the partial sequence of the 5.8S, whole sequence of the ITS2, and partial sequence of 28S rRNA gene using the HKY + G model of DNA substitution and rooted on

This study surveyed
Early reports on anisakid infection in marine fish showed that the parasites infect muscles around the belly area of fish (Petersen et al., 1993). This study also harvested the parasites from the abdominal cavity, stomach, and intestine of host fish. This is similar to the findings of various studies showing the gastrointestinal tract to be the preferred site among parasitic helminths (Gosselle et al., 2008; Onyedineke et al., 2010; Dauda et al., 2016). It is important to note that the method implemented for isolation and examination of the parasites was the common routine visual observation coupled with identification based on taxonomic descriptions from previous studies. However, the major disadvantage of this method is the tendency of missing out parasites with predilection sites other than the gastrointestinal tract, especially those that are embedded within tissues, which were not covered in this study. Another simple method was explored by various researchers, such as Shamsi and Suthar (2016), which will be beneficial if considered. Here, it was initially assumed that the helminth parasite larvae were identified as marine ascaridoid parasites based on the morphological characteristics. The results revealed that 22% of the fishes were detected with the parasite, with infection rates not significant between male and female fish.
Aggregation analysis revealed a highly aggregated distribution of parasites in the host population, with
Only a limited number of larvae were subjected to molecular analysis. We speculated that these larvae belong to the same taxonomic group due to gross morphological similarities. However, morphological characterization of ascaridoid larvae poses difficulty in identifying species, particularly from families with cryptic larval features. Although morphological characteristics are useful for species description and identification, the characterization may be limited because many macroscopic structures are produced infrequently and microscopic features of ascaridoid larvae are similar across different species (Mattiucci and Nascetti, 2006; Quiazon et al., 2008; Cavallero et al., 2014). Identifying which parasite species infects the fish populations is imperative since there are occurrence of zoonotic parasites that may affect the health of fish-eating populations.
Different genera of Anisakidae have been described using molecular tools more accurately than using classification based on their morphological characteristics and can be classified into either family Anisakidae or family Raphidascarididae (Kalay et al., 2009). Phylogenetic analyses based on the 5.8S-ITS2-28S DNA region confirmed close association of
The study revealed the presence of
Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

Figure 7

Figure 8

GenBank sequences of close species matches used for phylogeny construction.
Ascaridoid species | Host | Habitat | Geographical Origin | ITS2 | 18S rRNA |
---|---|---|---|---|---|
Marine | Tunisia | MT820022 | - | ||
Marine | China | - | MF072697 | ||
Marine | Italy | - | EF180082 | ||
Marine | Norway | AY826723 | - | ||
Marine | Japan | LC621351 | - | ||
Marine | China | - | MF072711 | ||
Marine | Egypt | HF911524 | - | ||
Bottlenose dolphin | Marine | China | KF673776 | - | |
Marine | Thailand | AB432909 | - | ||
Marine | Brazil | AY826724 | - | ||
Marine | Indonesia | KC928262 | - | ||
Marine | Brazil | EU327686 | - | ||
Marine | Portugal | JN005767 | - | ||
Marine | Mauritius | EU718473 | - | ||
Marine | South Africa | AY826725 | - | ||
Marine | Australia | - | EF180072 | ||
NM | NM | NM | - | AY702702 | |
NM | NM | NM | - | U94370 | |
Marine | Denmark | JX845137 | - | ||
Dolphin fish | Marine | South Korea | HQ702733 | - | |
Marine | China | - | MF072693 | ||
Marine | Brazil | - | JF718550 | ||
NM | NM | NM | - | U94374 | |
NM | NM | NM | - | U94376 | |
Marine | China | - | MF072705 | ||
Marine | China | - | MF072702 | ||
Marine | China | - | MF072707 | ||
Marine | Denmark | KM273088 | - | ||
NM | NM | NM | - | U94380 | |
Freshwater | Brazil | - | KP726276 | ||
Freshwater | Brazil | - | KP726278 | ||
Freshwater; brackish water | Finland | - | DQ503460 | ||
Freshwater | Brazil | - | KX859077 | ||
NM | Marine | South China Sea | JN102362 | - | |
Marine | China | KP326546 | MF072704 | ||
NM | NM | China | MF422212 | - | |
Marine | China | MH211584 | MF072692 | ||
NM | Marine | China | JF809816 | - | |
NM | NM | NM | - | U94381 | |
a |
Freshwater | Brazil | MK141032 | - | |
a |
NM | NM | NM | U94366 |
BLAST results for the 18S rRNA and 5.8S-ITS2-28S sequences of ascaridoid samples.
Specimen | 18S rRNA | 5.8S-ITS2-28S | Family | ||
---|---|---|---|---|---|
accession (18S/ITS2) | GenBank accession | % Match | GenBank accession | % Match | |
OK659774/ OK659786 | MF072704 ( |
99.84 | MF422212 ( |
99.18 | Raphidascarididae |
OK659775/ OK659787 | MF072711 ( |
99.52 | HF911524 ( |
99.88 | Anisakidae |
OK659776/ OK659788 | MF072704 ( |
99.84 | MF422212 ( |
99.79 | Raphidascarididae |
OK659777/ OK659789 | MF072711 ( |
99.68 | HF911524 ( |
99.88 | Anisakidae |
OK659778/ OK659790 | MF072704 ( |
99.52 | MF422212 ( |
99.79 | Raphidascarididae |
OK659779/ OK659791 | MF072704 ( |
99.52 | MF422212 ( |
99.59 | Raphidascarididae |
OK659780/ OK659792 | MF072711 ( |
99.20 | HF911524 ( |
98.66 | Anisakidae |
OK659781/ OK659793 | MF072711 ( |
99.68 | HF911524 ( |
98.79 | Anisakidae |
OK659782/ OK659794 | MF072711( |
99.52 | HF911524 ( |
99.87 | Anisakidae |
OK659783/ OK659795 | MF072711( |
99.68 | HF911524 ( |
98.86 | Anisakidae |
OK659784/ OK659796 | MF072711 ( |
99.68 | HF911524 ( |
98.61 | Anisakidae |
OK659785/ OK659797 | MF072711 ( |
99.68 | HF911524 ( |
99.87 | Anisakidae |
Corrected (TrNef + G) genetic distances of 18S rRNA gene and corrected (HKY + G) genetic distances of the 5.8S-ITS2-28S region of Anisakidae and Raphidascarididae nematodes.
Ascaridoid comparisons | Number of pairwise comparisons | 18S rRNA | 5.8S-ITS2-28S region | ||
---|---|---|---|---|---|
Minimum | Maximum | Minimum | Maximum | ||
67 (18S), 118 (ITS2) | 0 | 0.0182098 | 0 | 0.0077422 | |
494 (18S), 410 (ITS2) | 0 | 0.04259699 | 0.00493606 | 1.8776778 |
A Seinhorst Model Determined the Host-Parasite Relationships of Meloidogyne javanica Infecting Fenugreek cv. UM202Deladenus bonabensis n. sp. (Tylenchomorpha: Neotylenchidae) from East Azarbaijan province, northwestern Iran: A morphological and molecular phylogenetic studyEvaluation of Meloidogyne incognita andRotylenchulus reniformis nematode-resistant cotton cultivars with supplemental Corteva Agriscience nematicidesExternal and internal microbiomes of Antarctic nematodes are distinct, but more similar to each other than the surrounding environment Molecular and morphological characterization of Tylenchus zeae n. sp. (Nematoda: Tylenchida) from Corn (Zea mays ) in South Carolina