The phylogenetics of Anguillicolidae (Nematoda: Anguillicoloidea), swimbladder parasites of eels
© Laetsch et al.; licensee BioMed Central Ltd. 2012
Received: 27 October 2011
Accepted: 19 April 2012
Published: 4 May 2012
Anguillicolidae Yamaguti, 1935 is a family of parasitic nematode infecting fresh-water eels of the genus Anguilla, comprising five species in the genera Anguillicola and Anguillicoloides. Anguillicoloides crassus is of particular importance, as it has recently spread from its endemic range in the Eastern Pacific to Europe and North America, where it poses a significant threat to new, naïve hosts such as the economic important eel species Anguilla anguilla and Anguilla rostrata. The Anguillicolidae are therefore all potentially invasive taxa, but the relationships of the described species remain unclear. Anguillicolidae is part of Spirurina, a diverse clade made up of only animal parasites, but placement of the family within Spirurina is based on limited data.
We generated an extensive DNA sequence dataset from three loci (the 5' one-third of the nuclear small subunit ribosomal RNA, the D2-D3 region of the nuclear large subunit ribosomal RNA and the 5' half of the mitochondrial cytochrome c oxidase I gene) for the five species of Anguillicolidae and used this to investigate specific and generic boundaries within the family, and the relationship of Anguillicolidae to other spirurine nematodes. Neither nuclear nor mitochondrial sequences supported monophyly of Anguillicoloides. Genetic diversity within the African species Anguillicoloides papernai was suggestive of cryptic taxa, as was the finding of distinct lineages of Anguillicoloides novaezelandiae in New Zealand and Tasmania. Phylogenetic analysis of the Spirurina grouped the Anguillicolidae together with members of the Gnathostomatidae and Seuratidae.
The Anguillicolidae is part of a complex radiation of parasitic nematodes of vertebrates with wide host diversity (chondrichthyes, teleosts, squamates and mammals), most closely related to other marine vertebrate parasites that also have complex life cycles. Molecular analyses do not support the recent division of Anguillicolidae into two genera. The described species may hide cryptic taxa, identified here by DNA taxonomy, and this DNA barcoding approach may assist in tracking species invasions. The propensity for host switching, and thus the potential for invasive behaviour, is found in A. crassus, A. novaezelandiae and A. papernai, and thus may be common to the group.
653 - 663
51.43 - 52.22
30.55 - 31.64
651 - 663
51.89 - 52.53
30.55 - 34.91
30.55 - 31.27
In East Asia, A. crassus parasitises the native Japanese An. japonica as well as introduced and cultured An. anguilla and Anguilla rostrata (the American eel) [16, 17]. In its native host this parasite is minimally pathogenic, with a small adult body mass and low infection intensity . In the 1980s the parasite was introduced from Taiwan into Europe as a result of the live eel trade [18, 19], and the parasite colonised wild European eels, An. anguilla. In this new host the parasite attains much higher intensities. Anguillicoloides crassus infections have since spread through wild and farmed populations of An. anguilla in Europe and North Africa and are associated with cases of mass mortality when paired with environmental stressors such as high temperatures and low dissolved oxygen levels . A two-stage colonisation pattern has been described, consisting of rapid spread upon introduction into a water system followed by equilibration at ceiling levels . A. crassus was subsequently, and likely independently, introduced into populations of An. rostrata in North America. The introduction into North America is considered to have been from Japan [19, 22–24]. A. crassus has also recently been reported from the island of Réunion near Madagascar, where it was found in three indigenous Anguilla species . Being a “global invader” (sensu and ), it is considered an important pathogen of the economically relevant Atlantic eel species An. anguilla and An. rostrata[15, 27].
A. novaezelandiae Moravec and Taraschewski 1988 was first described from the short-fin eel Anguilla australis from New Zealand . However, an explant population of this species was recorded in Lake Bracciano in Italy following stocking with An. australis in 1975 . After the introduction of A. crassus into Lake Bracciano in 1993, the A. novaezelandiae population appears to have disappeared . A. australiensis Johnston and Mawston 1940 is recorded only from North-East Australia where it parasitises the native long-fin eel Anguilla reinhardtii, and no evidence for pathological effects on the host have been found . A. papernai Moravec and Taraschewski 1988 is endemic to South Africa and Madagascar where it parasitises the African long-fin eel Anguilla mossambica. A. papernai has been shown to be able to complete its life cycle in European intermediate and final hosts in laboratory infections . A. globiceps Yamaguti 1935 is known from certain prefectures of Japan and provinces of China [31, 32] where it infects wild populations of the Japanese eel An. japonica.
A phylogenetically correct and robust taxonomy of Nematoda is critical to understanding of biodiversity, biogeography and host-parasite coevolution . In the case of the Anguillicolidae, it is important to understand the phylogenetic distribution of traits associated with colonisation of and pathogenicity in new hosts, such as the plasticity of life-cycle traits and the increased per-host intensity for invasive A. crassus. DNA sequence data are good characters for taxonomic inference (i.e. DNA taxonomy) and for analysis of the deeper phylogenetic history of organisms. However, loci with different rates of evolution are typically required to resolve these different types of questions.
The nuclear small subunit ribosomal RNA gene (nSSU or 18S) has been used extensively for analysis of nematode phylogenetics [34, 35], but has limited resolution at the congeneric level, including within Anguillicolidae . In analyses of Nematoda using nSSU, A. crassus is placed within the entirely animal parasitic Clade III (suborder Spirurina in order Rhabditida sensu De Ley and Blaxter [34, 35, 37, 38]), consisting of the Gnathostomatomorpha, Ascaridomorpha, Oxyuridomorpha, Rhigonematomorpha and Spiruromorpha. De Ley and Blaxter [37, 38] placed Anguillicolidae, along with other taxa in Dracunculoidea, as incertae sedis within Spirurina. Focussed analyses of Spirurina using nSSU have identified significant conflicts with both classical systematics and the revisions of De Ley and Blaxter [39, 40], in particular the non-monophyly of 'Ascaridomorpha' and the placement of A. crassus together with the vertebrate-parasitic genus Gnathostoma, and not allied to other dracunculoids. The ichthyoparasitic cucullanid Truttaedactinis truttae was identified as sister to other analysed Spirurina . Moravec in his synoptic revision  removed the Anguillicolidae from the superfamily Dracunculoidea and erected the superfamily Anguillicoloidea with the family Anguillicolidae as its only member. In addition, based on morphological characters, the sub-genera Anguillicoloides and Anguillicola were proposed to be promoted to the taxonomic rank of genera .
The nuclear large subunit ribosomal RNA gene (nLSU or 28S), in particular the section spanning diversity loops D1 and D3, is an attractive alternative to nSSU, as there are both highly conserved parts and regions of more rapid evolution. nLSU D1-D3 sequences have also been used for analyses of Spirurina , but the available data are much more sparse than for nSSU. A third candidate locus is the mitochondrial gene cytochrome oxidase I (COX1, specifically the 5' half or 'Folmer region'), which has been proposed as a universal DNA taxonomy and DNA barcoding target for Metazoa [43–47]. This mitochondrial locus evolves very rapidly compared to nSSU, such that there is appreciable variation within species. COX1 is beginning to be used widely in analyses of nematode population structure and phylogeography , including Anguillicola/Anguillicoloides, generating extensive datasets with, however, low taxonomic coverage.
The assumption underpinning DNA taxonomy is the existence of a “barcoding gap”, resulting from the non-overlapping, discrete distribution of intra- and inter-specific variation of the DNA fragment analysed . Given this, unidentified specimens that differ by less than a threshold sequence divergence from reference voucher sequences (from specimens reliably identified to species, for example) can be assigned to that species. Symmetrically, specimens yielding sequences that are more divergent than the specified threshold can be allocated to a different taxon. For Metazoa, a general COX1 DNA taxonomy threshold of 2% has been proposed , but molecular operational taxonomic units (MOTU)  can be defined at any cutoff, and exploration of the pattern of MOTU count and cutoff is warranted for previously understudied groups . This approach can simplify the identification of potentially morphologically cryptic taxa and encourage revisionary taxonomy .
To investigate boundaries of species and genera in Anguillicolidae, and assess placement of this family in the diversity of Spirurina, we have generated data from nSSU, nLSU and COX1 for a large sample of Anguillicolidae encompassing all five nominal species. We assessed the empirical support for barcoding gaps in the three genes by analysing the impact of increasing sequence divergence thresholds, and the congruence between MOTU derived from different loci and species assignments. We analyse and compare the ability of the three genes to resolve the evolutionary relationships between the Anguillicolidae and test whether the morphological taxonomy of this family is supported by molecular evidence. Additionally, we infer a phylogenetic tree for Spirurina based on published and new nSSU sequences, to further refine our understanding of the relationships of the Anguillicolidae, and thus better inform hypotheses as to their origins.
Multi-locus marker development from Anguillicolidae
Location of sampling sites
Collector* and date**
Townsville, Queensland (Australia)
Brisbane, Queensland (Australia)
Deloraine, Tasmania (Australia)
Guangzhou, Guangdong (China)
Zhuhai, Guangdong (China)
Albufera de Valencia (Spain)
Rußheimer Altrhein (Germany)
EH, UW 2009
Natural water system, Wakayama (Japan)
HS 2006 2007
Ambatondrazaka (Madagascar) ***
Sniardwy Lake, Mikolajki (Poland)
Ribeira das Lampreias (Portugal)
Farm Dam, Fort Hare (South Africa)
Koonap River (South Africa)
Sunday’s River, Slagboom Dam (South Africa)
Sunday’s River, Darlington Dam (South Africa)
Great Fish River (South Africa)
Nahoon River (South Africa)
Asi River, Hatay (Turkey)
Sinyuan, Kaoping River (Taiwan)
Eel culturing pond, Budai, (Taiwan)
Species boundaries and MOTU in Anguillicolidae
Phylogenetic analyses of Spirurina
Phylogenetic relationships within the Anguillicolidae
Analysis of COX1 (Figure 5) divided the Anguillicolidae into two lineages with respect to morphospecies: a weakly supported clade formed by A. globiceps and A. papernai (with BPP of only 0.83 for COX1, but 1.00 for nLSU), and a robustly supported clade consisting of A. crassus, A. novaezelandiae and A. australiensis. A close relationship of the Oceania species A. novaezelandiae and A. australiensis was strongly supported (BPP for nLSU = 1.00 and for COX1 = 0.97). In the COX1 analyses, the A. novaezelandiae population sampled in Tasmania is more closely related to A. australiensis than to the A. novaezelandiae population from New Zealand. Mirroring the MOTU analyses, there is significant divergence within A. papernai, with three strongly supported clades (COX1 BPP = 1.0) as distinct from each other as are Tasmanian A. novaezelandiae from A. australensis.
Population diversity in Anguillicola crassus
Anguillicolidae species diversity analysed using DNA taxonomy
Primers that selectively amplify DNA from certain taxa are crucial to surveys in which parasite DNA has to be amplified over the background of host DNA. Studies on Anguillicolidae clearly fall into this category due to the haematophagous life-style of the pre-adult (larval) stages, and thus host blood cells in their digestive system. The primers used here [34, 53, 54] were broadly successful in amplification of the desired nematode fragments for subsequent sequencing. Thus DNA barcoding approaches are readily applicable to these nematodes.
Although there were five distinct nSSU sequences, corresponding to the five morphological species, these sequences differed by no more than 3 bases over the 867 bases sequenced. This close similarity raises the issue of potential misidentification due to sequencing errors, and thus nSSU may not be useful for species-level identification surveys of Anguillicolids. As expected the nLSU D2-D3 region displayed higher intra- and inter-specific variation among morphospecies than nSSU, with 9 unique sequences. This variation, along with the presence of indels, should also enable the design of species-specific amplification primers. This nLSU fragment has been successfully amplified from partially degraded specimens of A. crassus encapsulated within host tissue , and from morphologically indistinguishable life stages such as larvae. nLSU could be used for PCR-based diagnostics of anguillicoloidosis in eels and can serve as a reliable tool in life-cycle studies.
The methodology of MOTU definition, as implemented in jMOTU , explores a defined range of sequence divergence thresholds upon which a plateau of MOTU richness, i.e. potential barcoding gaps, can be identified. For both the nuclear loci studied, no clear evidence of a barcoding gap was found. The lack of interspecific diversity for the nSSU locus limits the utility of this marker at congeneric levels. Multiple within-(morpho)species nLSU sequence types were only observed in A. crassus (two sequence types) and A. papernai (four sequence types). One A. crassus sequence type is found globally, while the second is confined to one Japanese population, supporting the hypothesis of an Asian origin of this species .
The highest number of sequence types (54) was observed for the COX1 fragment, which has previously been used to investigate global population structure in A. crassus. We were able to obtain only a single COX1 fragment from A. globiceps that contained a correct open reading frame, despite multiple trials. This phenomenon has been observed in other surveys of this species (S. Wielgoss, pers. comm.), and may be due to preferential amplification of nuclear copies of mitochondrial DNA  or the existence of RNA editing of this mitochondrial transcript, as has been observed in other nematodes . Sequencing of COX1 transcripts and the mitochondrial genome may be informative.
Our A. crassus COX1 data, when added to those previously determined, reinforces the view that the invasion of West European and North American hosts has only happened a few times (likely once for each location) and that the origins of the invading parasite differ. North American A. crassus are robustly linked to Japanese specimens, and European to Taiwanese. By sampling additional locations across China and Japan, we identified many additional, unique COX1 haplotypes limited to East Asia, affirming this as the area of highest diversity, and thus the likely origin of diversity of this species. Interestingly, Turkish A. crassus may have distinct origins, as some haplotypes group with Chinese and Taiwanese sequences not previously observed in Western Europe.
The COX1 locus of the Anguillicolidae exhibits a likely barcoding gap with over 5.6% sequence divergence, exceeding by more than twofold the proposed threshold (2%) for Metazoa . Although the species A. crassus, A. globiceps and A. australiensis were stably grouped into individual MOTUs within this plateau, sequences of both A. novaezelandiae and A. papernai were found in multiple, distinct COX1 MOTUs. These distinct MOTUs within a single nominal species are possible evidence of cryptic speciation. The Madagascar populations of A. papernai may be the source of introduction into South Africa, since all sequence types from the five sampling sites in South Africa occur in, or are closely related to, those from Madagascar. This inference is supported by the pattern observed for the four nLSU sequence types of A. papernai, which are all present in Madagascar. We note that the mitochondrial locus COX1 has uniparental inheritance and thus may not be an unbiased reporter of phylogeographic history, but propose these hypotheses as testable inferences from our data.
Specimen SNR118 was a single, unidentified larva from the swim bladder serosa of a South African A. mossambica, and the whole specimen was used for DNA extraction. Its nSSU sequence showed that it is a member of Spirurina B. Sampling site, host species and phylogenetic position in the Spirurina suggest this specimen belongs to or is closely affiliated to Paraquimperia africana (Seuratoidea) .
DNA taxonomy and molecular phylogenetic systematics of the Anguillicolidae
We found no support for the division of the Anguillicolidae into the two genera Anguillicola and Anguillicoloides as proposed by Moravec , as Anguillicola globiceps was recovered within the radiation of Anguillicoloides species, suggesting that that the morphological criteria used to erect these putative genera (spinosity of cuticular ornamentation and the structure of the oesophagus) may not be phylogenetically informative. As noted previously , modifications of the structure of the oesophagus, a trait related to the mode of nutrition, may occur independently in the course of trophic adaptations. Our data support a single generic division for the species in Anguillicolidae, which by priority should be called Anguillicola, restoring Anguillicoloides crassus to Anguillicola crassus and ensuring continuity with historical literature on this important species.
Based on the analysis of the DNA sequences of three genes, a reevaluation of the taxonomy of the Anguillicolidea may be necessary. We identified eight discrete COX1 MOTUs over a wide range of divergence, up to 5%, suggesting that there may be eight species-level taxa represented if a barcoding cutoff of 2% is accepted. Three MOTUs are congruent with morphological species identifications: A. crassus, A. globiceps and A. australiensis. However specimens unambiguously assigned to the morphological species A. novaezelandiae and A. papernai are found in multiple MOTUs (two MOTUs for A. novaezelandiae and three for A. papernai). These additional MOTUs could represent morphologically cryptic, or previously unrecognised species.
In the case of A. novaezelandiae phylogenetic analyses of the COX1 locus suggest the paraphyly of this species, as the specimens from New Zealand (that form one MOTU) are robustly placed as the sister group to the clade consisting of A. novaezelandiae specimens from Tasmania (the second MOTU) plus A. australiensis. Since nLSU D2-D3 sequences are not available for A. novaezelandiae from New Zealand, additional representatives should be sampled to more thoroughly test the monophyly of this species.
A. papernai displays the greatest number of nLSU D2-D3 sequence types among the Anguillicolidae and its COX1 locus displays a clear division into three clades and is represented by three MOTUs. Whether these are distinct species level taxa, or merely distinct diverse populations of a widespread species will require analysis of additional loci. A. papernai is found in well separated watersheds, with no linking waterways, and an absence of human-induced admixture through stocking. The diversity within this species, together with the observation that it is able to complete its life cycle in European eel hosts , supports the need for further research.
A similar branching pattern is observed for both nLSU and COX1 for all species analysed except A. crassus. Analyses of the two nuclear loci revealed A. crassus to be sister to the other sampled species within Anguillicolidae. However the mitochondrial COX1 gene placed it as the sister taxon to the clade comprising A. australiensis and A. novaezelandiae. Despite the low number of informative sites in the nLSU and nSSU loci, we favour the hypothesis based on the nuclear loci, as no close outgroup sequences were available for COX1 and the sequences used (Toxocara cati and Strongylida sp.) may have impacted on proper rooting of the COX1 tree. Reanalyses with COX1 data from Seuratoidea or Gnathostomatinae would address this issue.
These data suggest a scenario of Asian origin for the Anguillicolidae, since A. papernai is the only endemic species west of the 80 degree line of longitude. Speciation of A. crassus could have taken place in East Asian waters, followed by the cladogeneses of the Oceanian species and the clade comprising A. papernai and A. globiceps. The biogeographic distribution of A. papernai observed today could thus be explained by a host-related dispersal of its ancestor.
Phylogenetic relationships within Spirurina
Spirurina includes parasites with a direct lifecyle as well as those that include biological vector and transport (paratenic) hosts in complex, multi-species systems . Although some applications of parsimony principles might suggest that the simpler, direct life cycle should be ancestral for the group, molecular phylogenetic analyses robustly place two clades of parasites (Spirurina A and B) that have complex life histories as successive sister groups, at the base of the Spirurina, to clades that include direct life cycle parasites. This phylogenetic hypothesis favors intepretation of complex life histories as the ancestral state for the group. We found that members of the ichthyoparasitic family Cucullanidae (Seuratoidea: C. robustus, D. mexicanus and T. truttae) form a well-supported clade (Spirurina A; Figure 3) that is sister to the remaining Spirurina [40, 41]. Anguillicolidae are placed in Spirurina B along with certain members of the superfamily Seuratoidea and all sampled Gnathostomatoidea. Spirurina C includes the remaining taxa including the abundantly sampled Ascaridomorpha and Spiruromorpha. Previous analyses of taxa in Spirurina C have suggested that many classical groups are non-monophyletic, in particular "Ascaridomorpha" and "Spiruromorpha" [39, 40, 42], a result echoed here. The position of the three Gnathostoma species (G. binucleatum, G. neoprocyonis and G. turgidum) outside of Spirurina C has been observed previously  and was recovered consistently in all analyses. Members of the other two genera of the Gnathostomatinae, i.e. T. tiara and E. overstreeti, are also displaced, and share direct ancestors with Linstowinema sp. (Seuratoidae) and members of Anguillicolidae, respectively. This result indicates the paraphyly of both Gnathostomatinae and Seuratoidea. The Dracunculoidea is also rendered paraphyletic, as members are found in Spirurina B and Spirurina C. This phylogenetic hypothesis reveals an enormous definitive host diversity within the Spirurina B, comprising fresh-water teleosts (Anguillicolidae), chondrichthyes (Echinocephalus), mammals (Gnathostoma and Linstowinema) and squamates (Tanqua), and an even wider diversity of vector and paratenic hosts. Denser taxon sampling in the Spirurina, with special emphasis on the morphologically diverse Spirurina B and C, is required to fully explore this fascinating group.
We have investigated the genetic diversity within, and genetic distinctiveness of the five described species of the family Anguillicolidae. We found no support for the erection of two genera within the family, and identified two species (A. novaezelandiae and A. papernai) where within-species divergences and phylogenetic tree topology suggests the presence of distinct cryptic taxa. The role of host vicariance in the speciation of A. globiceps and A. papernai is an intriguing topic for future study, as is the evolutionary history of the Oceanian species (A. australiensis and the two MOTUs within A. novaezelandiae). We revisited the phylogeny of the Spirurina and identified a clade of vertebrate-parasitic taxa that includes Anguillicoloidae and members of other families, which must therefore be suspected of being paraphyletic. Our analyses highlight the possibility that A. papernai might transfer to and be pathogenic in new, economically important, eel hosts, and the data generated here will, we hope, act as reference for future DNA barcoding surveys of eel swim bladder parasites worldwide.
DNA extraction and sequencing
Sampling and identification of nematode specimens was performed by José Lino Costa, Kerstin Geiss, Ercüment Genç, Emanuel Heitlinger, Albert Keim, Lea Perseke, Pilar Muñoz Ruíz, Hiroshi Sato, Björn Schäffner, Horst Taraschewski, Urszula Weclawski and Olaf Weyl. A total of 150 anguillicolid nematodes were extracted from the swim bladders of their respective hosts from 21 different locations (see Table 2). In addition, an unidentified nematode larva (SNR118) was extracted from the serosa of the swim bladder of a specimen of A. mossambica from South Africa. DNA was prepared from single nematodes as described in . Lysates were used directly as templates in PCR reactions. For all loci, 25 μl polymerase chain reactions (PCRs) were carried out using 0.1 μl Taq DNA Polymerase (5 units/μl) (Quiagen, Hilden, Germany), 2.7 μl 10x PCR Buffer (containing 15 mM MgCl2) (Quiagen), 2.7 μl 2 mM dNTP (2 mM dATP, dTTP, dGTP, dCTP), 0.4 μl of each PCR primer (10 μM), 2 μl template DNA and 17.1 μl ddH2O (Milli-Q). For partial nSSU amplification the forward primer SSU_F04 5’-GCTTGTCTCAAAGATTAAGCC-3’ and the reverse primer SSU_R26 5’-CATTCTTGGCAAATGCTTTCG-3’  were used. Amplification of the nLSU D2-D3 region was carried out using forward primer D2A 5’-ACAAGTACCGTGAGGGAAAGT-3’ and the reverse primer D3B 5’-TGCGAAGGAACCAGCTACTA-3’ . The cytochrome c oxidase subunit I (COX1) was amplified using the forward primer LCO1490 5’-GGTCAACAAATCATAAAGATATTGG-3’ and the reverse primer HCO2198 5’-TAAACTTCAGGGTGACCAAAAAAT-3’ . PCR products were purified using shrimp alkaline phosphatase and Escherichia coli exonuclease I (USB Corporation, USA) as described in  and sequenced on an automated ABI Prism 3730 Genetic Analyzer using ABI BigDye v3.1 Terminator sequencing chemistry (Applied Biosystems, Foster City, CA) in the GenePool Genomics Facility, Edinburgh (http://genepool.bio.ed.ac.uk). Sequencing of each PCR product was carried out in both directions to minimise PCR artefacts, ambiguities and base-calling errors. Directly sequenced COX1 PCR products of A. globiceps displaying indels disrupting the open reading frame were cloned (PCRII Topo TA cloning kit, Invitrogen) and only specimens for which a single, correct open reading frame (ORF) bearing sequence was identified have been included in the analyses. Raw ABI chromatograph files of the sequences were processed using trace2seq.pl (a perl program that uses phred [63, 64] to identify high-quality base calls; A. Anthony and M. Blaxter, unpublished). After screening for contaminants using NCBI BLAST , the high quality forward and reverse sequences of each gene from each sample were aligned and a consensus sequence was inferred using a phred score of 30, i.e. a 99.9% accuracy of a base-call, as the detection threshold, resulting in an increase in sequence lengths and an improvement in sequence reliability. The 352 sequences were deposited in the EMBL database (accession numbers JF805371 – JF805722). Aligned sequence datasets for all analyses are available as Additional file 1 and also at datadryad.org under doi:10.5061/dryad.8h5p7p00.
Construction of the datasets
For each of the three genes, an individual dataset was created (nSSU, nLSU and COX1) containing the respective consensus sequences from different individual nematode specimens. Trimmed datasets containing only sequence types derived from the subset of 79 specimens for which all three genes were successfully sequenced (termed nSSU*, nLSU* and COX1*). Additional COX1 sequences , obtained using the same primers, were retrieved from GenBank and included in the COX1 dataset. Sequences contained in the nSSU and COX1 datasets were of uniform length (787 bp and 550 bp, respectively) and could be aligned unambiguously using CLUSTALW 2.0.9  with the penalties for gap opening and extension set to 10 and 0.2, respectively. Sequences in the nLSU alignment showed several insertion/deletion events (indels), which were binary coded for phylogenetic analyses  implemented in Seqstate 1.4 . The data was included in the phylogenetic analysis as a binary partition under the restriction site model with the ascertainment bias set to variable, as suggested by .
jMOTU 1.0.7  is a Java application that defines molecular operational taxonomic units (MOTUs) based on nucleotide sequences and a user defined range of cutoff values, the maximum number of base pair differences between two sequences, by using global alignment. MOTUs were defined at cut-off values ranging from 0 to 14% sequence divergence in intervals of 1 base pair, on sequences from specimens for which all three genes were sequenced (nSSU*, nLSU* and COX1*). To investigate the degree of influence sampling depth has on MOTU richness and membership, the results were compared to those of the MOTU analysis performed on the COX1 dataset.
For phylogenetic analysis of the Spirurina, additional sequences were retrieved from GenBank [39–41] and aligned to the anguillicolid sequences. New sequences from spirurine nematodes from a dataset provided by S. Nadler were also included [accession numbers JF934725–JF934737]. Regions within the resulting Spirurina alignment in which determination of homology was ambiguous (i.e. many apparently independent insertion-deletion events over 10 contiguous bases of the alignment) were excluded. For both the nLSU and COX1 data, outgroups were chosen based on closest sequence matches in the public databases, while for nSSU outgroups were chosen based on previous analyses. Bayesian phylogenetic inference was carried out using MrBayes 3.1.2 . Phylogenies for individual datasets were inferred under the GTR + I + Γ model of sequence evolution, partitioned by codon position for COX1 and partitioned by datatype for nLSU (i.e. nucleotide data and binary insertion-deletion data). For each analysis, two independent Markov chain Monte Carlo (MCMC) runs of four Metropolis-coupled chains were performed with the gamma shape parameter, the proportion of invariable sites, base frequencies and substitution rates unlinked across partitions and assuming default priors. Chains were sampled every 1,000 generations for 7.5x106 (nSSU), 7.5x106 (nLSU) and 5x106 (COX1) generations. Convergence of Markov chains was assessed using Tracer 1.4 (Rambaut A, Drummond AJ (2007); available from http://beast.bio.ed.ac.uk/Tracer) and saved trees from the first 750,000 (nSSU, nLSU), 500,000 (Spirurina) and 300,000 (COX1) generations were discarded as burn-in.
Statistical parsimony network analysis
Population designations for Anguillicola crassus population structure analyses with COX1
Åland Islands (Finland)
Kullen, Øresund/Kattegat (Sweden)
Slapton Ley, Cornwall (Great Britain)
Lake Neagh (Great Britain)
Lough Dergh, Shannon (Ireland)
Rußheimer Altrhein (Germany)
Sniardwy Lake, Mikolajki (Poland)
Bois Joli, Frémur (France)
Angers, Loire (France)
Camargue, Rhône (France)
Albufera de Valencia (Spain)
Ribeira das Lampreias (Portugal)
Roma, Tiber (Italy)
Lake Bracchiano (Italy)
Asi River, Hatay (Turkey)
St. Jones River (New Jersey, USA)
Tung-Chiang, Kao-Ping River (Taiwan)
Sinyuan, Kaoping River (Taiwan)
Eel culturing pond, Budai, (Taiwan)
Zhuhai, Guangdong (China)
Guangzhou, Guangdong (China)
Mikawa Bay (Japan)
Yamaguchi, Fushino (Japan)
Natural water system, Wakayama (Japan)
We gratefully thank our collegues José Lino Costa, Kerstin Geiss, Ercüment Genç, Albert Keim, Lea Perseke, Pilar Muñoz Ruíz, Hiroshi Sato, Björn Schäffner, Urszula Weclawski and Olaf Weyl for collection and identification of samples. We thank Lin Hao-Ran and Han Yu-San for support during collections. Collection of South African and Madagascar samples was financially and logistically supported by the South African Institute for Aquatic Biodiversity (SAIAB) and the National Research Foundation of South Africa (NRF). HT was supported by the DFG and Karlsruher Universitätsgesellschaft. We thank the GenePool Genomics Facility, Edinburgh (http://genepool.bio.ed.ac.uk) for sequencing. Thanks are due to Sebastien Wielgoss and Francois Lefebrve for comments on the manuscript.
- Kuwahara A, Itagaki H, Itagaki H: Studies of a nematode parasitic in the air bladder of the eel: 1. Description of Anguillicola crassus n. sp. (Philometridae, Anguillicolidae). Japanese, Journal of Parasitology. 1974, 23: 275-279.Google Scholar
- Moravec F, Taraschewski H: Revision of the genus Anguillicola Yamaguti, 1935 (Nematoda: Anguillicolidae) of the swimbladder of eels, including descriptions of two new species, A. novaezelandiae sp. n. and A. papernai sp. n. Folia Parasitol. 1988, 35: 125-146.PubMedGoogle Scholar
- Moravec F: Dracunculoid and anguillicoloid nematodes parasitic in vertebrates. 2006, Academia, PragueGoogle Scholar
- De Charleroy D, Grisez L, Thomas K, Belpaire C, Ollevier FP: The life cycle of Anguillicola crassus. Dis Aquat Organ. 1990, 8: 77-84.View ArticleGoogle Scholar
- Wang P, Zhao Y: Observations on the life history of Anguillicola globiceps (Nematoda: Anguillicolidae). Acta Zoologica Sinica. 1980, 26: 243-249.Google Scholar
- Moravec F, Konecny R: Some new data on the intermediate and paratenic host of the Nematode Anguillicola crassus Kuwahara, Niimi et Itagaki, 1974 (Dracunculoidea), a swimbladder parasite of eels. Folia Parasitol. 1994, 41: 65-70.Google Scholar
- Taraschewski H, Boomker J, Moravec F: Studies on the morphology and ecology of Anguillicola papernai (Nematoda: Anguillicolidae) parasitizing the swimbladder of the African longfin eel, Anguilla mossambica, and on the helminth community in this eel. Dis Aquat Organ. 2005, 63: 185-195.PubMedView ArticleGoogle Scholar
- Moravec F, Nagasawa K, Miyakawa M: First record of ostracods (Ostracoda) as natural intermediate hosts of Anguillicola crassus (Nematoda: Anguillicolidae), a pathogenic swimbladder parasite of eels (Anguilla spp.). Dis Aquat Organ. 2005, 66: 171-173.PubMedView ArticleGoogle Scholar
- Székely C: Experimental studies on the infectivity of Anguillicola crassus third-stage larvae (Nematoda) from paratenic hosts. Folia Parasitol. 1996, 4: 183-186.Google Scholar
- Moravec F: Aquatic invertebrates (snails) as new paratenic hosts of Anguillicola crassus (Nematoda: Dracunculoidea) and the role of paratenic hosts in the life cycle of this parasite. Dis Aquat Organ. 1996, 27: 237-239.View ArticleGoogle Scholar
- Moravec F, Skoríková B: Amphibians and larvae of aquatic insects as new paratenic hosts of Anguillicola crassus (Nematoda: Dracunculoidea), a swimbladder parasite of eels. Dis Aquat Organ. 1988, 34: 217-222.View ArticleGoogle Scholar
- Kirk RS: The impact of Anguillicola crassus on European eels. Fisheries Management and Ecology. 2003, 10: 385-394. 10.1111/j.1365-2400.2003.00355.x.View ArticleGoogle Scholar
- Taraschewski H: Host and parasites as aliens. J Helminthol. 2006, 80: 99-128. 10.1079/JOH2006364.PubMedView ArticleGoogle Scholar
- Kennedy CR: The pathogenic helminth parasites of eels. J Fish Dis. 2007, 30: 319-334. 10.1111/j.1365-2761.2007.00821.x.PubMedView ArticleGoogle Scholar
- Székely C, Palstra A, Molnár K, Thillart G: Impact of the Swim-Bladder Parasite on the Health and Performance of European Eels. In: Spawning Migration of the European Eel. 2009, Springer, Netherlands, 201-226. vol. 30View ArticleGoogle Scholar
- Egusa S: Changes of epizootics in pond eels in Japan. Bulletin-Office international des épizooties. 1968, 69: 1035-1044.PubMedGoogle Scholar
- Han YS, Chang YT, Taraschewski H, Chang SL, Chen CC, Tzeng WNT: The swimbladder parasite Anguillicola crassus in native Japanese eels and exotic American eels in Taiwan. Zoological Studies. 2008, 47: 667-675.Google Scholar
- Neumann W: Schwimmblasenparasit Anguillicola bei Aalen. Fischer und Teichwirt. 1985, 36: 322-Google Scholar
- Wielgoss S, Taraschewski H, Meyer A, Wirth T: Population structure of the parasitic nematode Anguillicola crassus, an invader of declining North Atlantic eel stocks. Mol Ecol. 2008, 17: 3478-3495. 10.1111/j.1365-294X.2008.03855.x.PubMedView ArticleGoogle Scholar
- Molnár K, Baska F, Csaba G, Glávits R, Székely C: Pathological and histopathological studies of the swimbladder of eels Anguilla anguilla infected by Anguillicola crassus (Nematoda: Dracunculoidea). Dis Aquat Organ. 1993, 15: 41-50.View ArticleGoogle Scholar
- Lefebvre F, Crivelli AJ: Anguillicolosis: dynamics of the infection over two decades. Dis Aquat Organ. 2004, 62: 227-232.PubMedView ArticleGoogle Scholar
- Fries LT, Williams DJ, Johnson SKEN: Occurrence of Anguillicola crassus, an Exotic Parasitic Swim Bladder Nematode of Eels, in the Southeastern United States. Trans Am Fish Soc. 1996, 125: 794-797. 10.1577/1548-8659(1996)125<0794:NOOCAE>2.3.CO;2.View ArticleGoogle Scholar
- Barse AM, Secor DH: An exotic nematode parasite of the Amercian eel. Fisheries. 1999, 24: 6-10.View ArticleGoogle Scholar
- Rockwell LS, Jones KMM, Cone DK: First reocrd of Anguillicoloides crassus (Nematoda) in American Eels (Anguilla rostrata) in Canadian Estuaries, Cape Breton, Nova Scotia. J Parasitol. 2009, 96: 483-486.View ArticleGoogle Scholar
- Sasal P, Taraschewski H, Valade P, Grondin H, Wielgoss S, Moravec F: Parasite communities in eels of the Island of Reunion (Indian Ocean): a lesson in parasite introduction. Parasitol Res. 2008, 102: 1343-1350. 10.1007/s00436-008-0916-5.PubMedView ArticleGoogle Scholar
- Colautti RI, MacIsaac HJ: A neutral terminology to define invasive species. Diversity and Distributions. 2004, 10: 135-141. 10.1111/j.1366-9516.2004.00061.x.View ArticleGoogle Scholar
- Vogel G: Europe Tries to Save Its Eels. Science. 2010, 329: 505-10.1126/science.329.5991.505.PubMedView ArticleGoogle Scholar
- Paggi L, Orecchia P, Minervini R, Mattiucci S: Sulla comparsa di Anguillicola australiensis Johnston e Mawson, 1940 (Dracunculoidea: Anguillicolidae) in Anguilla anguilla del Lago di Bracciano. Parasitologia. 1982, 24: 139-144.Google Scholar
- Moravec F, Di Cave D, Orecchia P, Paggi L: Present occurrence of Anguillicola novaezelandiae (Nematoda: Dracunculoidea) in Europe and its development in the intermediate host. Folia Parasitol. 1994, 41: 203-208.PubMedGoogle Scholar
- Kennedy CR: The distribution and abundance of the nematode Anguillicola australiensis in eels Anguilla reinhardtii in Queensland, Australia. Folia Parasitol. 1994, 41: 279-285.Google Scholar
- Yamaguti S: Studies on the helmith fauna of Japan, part 9. Nematodes of fishes. Japanese Journal of Zoology. 1933, 6: 8-Google Scholar
- Nagasawa K, Kim Y-G, Hirose H: Anguillicola crassus and A. globiceps (Nematoda: Dracunculoidea) parasitic in the swimbladder of eels (Anguilla japonica and A. anguilla) in East Asia: a review. Folia Parasitol. 1994, 41: 127-137.Google Scholar
- Adams BJ: Species concepts and the evolutionary paradigm in modern nematology. Journal of Nematology. 1998, 30: 1-21.PubMedPubMed CentralGoogle Scholar
- Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, Vanfleteren JR, Mackey LY, Dorris M, Frisse LM, et al: A molecular evolutionary framework for the phylum Nematoda. Nature. 1998, 392: 71-75. 10.1038/32160.PubMedView ArticleGoogle Scholar
- Holterman M, van der Wurff A, van den Elsen S, van Megen H, Bongers T, Holovachov O, Bakker J, Helder J: Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown Clades. Mol Biol Evol. 2006, 23: 1792-1800. 10.1093/molbev/msl044.PubMedView ArticleGoogle Scholar
- Hirose H, Yabu T, Hirono I, Aoki T: The phylogeny of Anguillicola crassus and A. globiceps based on partial 18 S ribosomal RNA sequences. J Fish Dis. 1998, 21: 265-271. 10.1046/j.1365-2761.1998.00103.x.View ArticleGoogle Scholar
- De Ley P, Blaxter ML: A new system for the Nematoda: combining morphological characters with molecular trees, and translating clades into ranks and taxa. Proceedings of the fourth International Congress of Nematology. Edited by: Cook R, Hunt DJ. 2004, Brill (Leiden, Netherlands), Tenerife, 865-Google Scholar
- De Ley P, Blaxter ML: Systematic position and phylogeny. The Biology of Nematodes. Edited by: Lee D. 2002, Taylor & Francis, London, 1-30.View ArticleGoogle Scholar
- Wijová M, Moravec F, Horák A, Lukes J: Evolutionary relationships of Spirurina (Nematoda: Chromadorea: Rhabditida) with special emphasis on dracunculoid nematodes inferred from SSU rRNA gene sequences. Int J Parasitol. 2006, 36: 1067-1075. 10.1016/j.ijpara.2006.04.005.PubMedView ArticleGoogle Scholar
- Nadler SA, Carreno RA, Mejia-Madrid H, Ullberg J, Pagan C, Houston R, Hugot JP: Molecular phylogeny of clade III nematodes reveals multiple origins of tissue parasitism. Parasitology. 2007, 134: 1421-1442. 10.1017/S0031182007002880.PubMedView ArticleGoogle Scholar
- van Megen H, van den Elsen S, Holterman M, Karssen G, Mooyman P, Bongers T, Holovachov O, Bakker J, Helder J: A phylogenetic tree of nematodes based on about 1200 full-length small subunit ribosomal DNA sequences. Nematology. 2009, 11: 927-950. 10.1163/156854109X456862.View ArticleGoogle Scholar
- Kuzmin Y, Tkach VV, Snyder SD, Bell JA: Camallanus Railliet et Henry, 1915 (Nematoda, Camallanidae) from Australian freshwater turtles with descriptions of two new species and molecular differentiation of known taxa. Acta Parasitol. 2011, 56: 213-226. 10.2478/s11686-011-0015-0.View ArticleGoogle Scholar
- Dorris M, De Ley P, Blaxter ML: Molecular analysis of nematode diversity and the evolution of parasitism. Parasitol Today. 1999, 15: 188-193. 10.1016/S0169-4758(99)01439-8.PubMedView ArticleGoogle Scholar
- Floyd R, Abebe E, Papert A, Blaxter M: Molecular barcodes for soil nematode identification. Mol Ecol. 2002, 11: 839-850. 10.1046/j.1365-294X.2002.01485.x.PubMedView ArticleGoogle Scholar
- Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Series B. 2003, 270: 313-321.Google Scholar
- De Ley P, De Ley IT, Morris K, Abebe E, Mundo-Ocampo M, Yoder M, Heras J, Waumann D, Rocha-Olivares A, Jay Burr AH, et al: An integrated approach to fast and informative morphological vouchering of nematodes for applications in molecular barcoding. Philos Trans R Soc Lond B Biol Sci. 2005, 360: 1945-1958. 10.1098/rstb.2005.1726.PubMedPubMed CentralView ArticleGoogle Scholar
- Bhadury P, Austen MC, Bilton DT, Lambshead PJD, Rogers AD, Smerdon GR: Development and evaluation of a DNA-barcoding approach for the rapid identification of nematodes. Mar Ecol Prog Ser. 2006, 320: 1-9.View ArticleGoogle Scholar
- Derycke S, Vanaverbeke J, Rigaux A, Backeljau T, Moens T: Exploring the use of cytochrome oxidase c subunit 1 (COI) for DNA barcoding of free-living marine nematodes. PLoS One. 2010, 5: e13716-10.1371/journal.pone.0013716.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyer CP, Paulay G: DNA barcoding: error rates based on comprehensive sampling. PLoS Biol. 2005, 3: e422-10.1371/journal.pbio.0030422.PubMedPubMed CentralView ArticleGoogle Scholar
- Hebert PD, Ratnasingham S, deWaard JR: Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc Lond B Biol Sci. 2003, 270 (Suppl 1): S96-S99.View ArticleGoogle Scholar
- Jones M, Ghoorah A, Blaxter M: jMOTU and Taxonerator: turning DNA Barcode sequences into annotated operational taxonomic units. PLoS One. 2011, 6: e19259-10.1371/journal.pone.0019259.PubMedPubMed CentralView ArticleGoogle Scholar
- de Perez-Ponce León G, Nadler SA: What we don't recognize can hurt us: a plea for awareness about cryptic species. J Parasitol. 2010, 96: 453-464. 10.1645/GE-2260.1.View ArticleGoogle Scholar
- Nunn GB: Nematode Molecular Evolution: An Investigation of Evolutionary Patterns Among Nematodes Based on DNA Sequences. PhD thesis. 1992, University of Nottingham, NottinghamGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994, 3: 294-PubMedGoogle Scholar
- Heitlinger EH, Laetsch DR, Weclawski U, Han Y, Taraschewski H: Massive encapsulation of larval Anguillicoloides crassus in the intestinal wall of Japanese eels. Parasites and Vectors. 2009, 2: 48-10.1186/1756-3305-2-48.PubMedPubMed CentralView ArticleGoogle Scholar
- Vanfleteren JR, Vierstraete AR: Insertional RNA editing in metazoan mitochondria: the cytochrome b gene in the nematode Teratocephalus lirellus. RNA. 1999, 5: 622-624. 10.1017/S135583829999009X.PubMedPubMed CentralView ArticleGoogle Scholar
- Hebert PDN, Ratnasingham S, de Waard JR: Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society of London Series B. Biological Sciences. 2003, 270: S96-S99.PubMedPubMed CentralView ArticleGoogle Scholar
- Cernotíková E, Horák A, Moravec F: Phylogenetic relationships of some spirurine nematodes (Nematoda: Chromadorea: Rhabditida: Spirurina) parasitic in fishes inferred from SSU rRNA gene. Folia Parasitol. 2011, 58: 135-148.PubMedView ArticleGoogle Scholar
- Havermans C, Nagy ZT, Sonet G, De Broyer C, Martin P: Incongruence between molecular phylogeny and morphological classification in amphipod crustaceans: A case study of Antarctic lysianassoids. Mol Phylogenet Evol. 2010, 55: 202-209. 10.1016/j.ympev.2009.10.025.PubMedView ArticleGoogle Scholar
- Minegishi Y, Aoyama J, Inoue JG, Miya M, Nishida M, Tsukamoto K: Molecular phylogeny and evolution of the freshwater eels genus Anguilla based on the whole mitochondrial genome sequences. Mol Phylogenet Evol. 2005, 34: 134-146. 10.1016/j.ympev.2004.09.003.PubMedView ArticleGoogle Scholar
- Anderson RC: Nematode Parasites of Vertebrates. Their development and transmission. 2000, C.A.B. International, Wallingford, 2View ArticleGoogle Scholar
- Hanke M, Wink M: Direct DNA sequencing of PCR-amplified vector inserts following enzymatic degradation of primer and dNTPs. Biotechniques. 1994, 17 (5): 858-PubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X Version 2.0. Bioinformatics. 2007, 23: 2947-10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Simmons MP, Ochoterena H: Gaps as characters in sequence-based phylogenetic analyses. Syst Biol. 2000, 49: 369-381.PubMedView ArticleGoogle Scholar
- Müller K: SeqState - primer design and sequence statistics for phylogenetic DNA data sets. Appl Bioinformatics. 2005, 4: 65-69. 10.2165/00822942-200504010-00008.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.PubMedView ArticleGoogle Scholar
- Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.