- Open Access
Morphology is not a reliable taxonomic tool for the genus Lernaea: molecular data and experimental infection reveal that L. cyprinacea and L. cruciata are conspecific
Parasites & Vectors volume 12, Article number: 579 (2019)
Species belonging to the genus Lernaea are cosmopolitan parasites that can infect many different freshwater fish hosts. Due to a high degree of morphological intraspecific variability and high levels of interspecific similarities, their classification is extremely difficult and controversial. Although the suitability of the shape of cephalic horns has been questioned decades ago by some experimental infection studies, this character still plays the central role in the identification of Lernaea spp.
We used the nominal species Lernaea cyprinacea and Lernaea cruciata to test the hypothesis that the shape of the anchor can exhibit host-induced morphological variability, and that the two taxa may be synonymous.
We examined 517 wild or farmed specimens of five host fish species (four cyprinids and a mosquitofish), and found that all 16 parasite specimens collected from mosquitofish could be morphologically identified as L. cruciata, whereas the remaining 25 parasite specimens were all identified as L. cyprinacea. We experimentally infected goldfish and mosquitofish specimens with offspring (copepodids) of a single L. cyprinacea specimen: the adult parasites from goldfish were morphologically identified as L. cyprinacea, and those from mosquitofish as L. cruciata. We then used molecular data to corroborate that all these specimens are conspecific.
Our results suggest that L. cyprinacea and L. cruciata may be synonyms, misidentified as different species as a result of host-induced morphological variation. Given the current shortage of molecular data for the genus Lernaea, in order to resolve the taxonomy of this genus (determine the exact number of species), future studies should aim to sequence as much molecular data as possible, and conduct further experimental infections.
Cosmopolitan parasitic copepods belonging to the genus Lernaea Linnaeus, 1758 (Cyclopoida: Lernaeidae) can infect many different freshwater fish species [1,2,3], causing lernaeosis, a disease that can cause serious pathogenic effects on their hosts. The taxonomy of this genus is still largely unresolved due to the existence of an exceptionally large number (109) of recorded nominal species , a small number of morphological traits useful for species identification, a high degree of intraspecific morphological variability, and a lack of clear morphological distinction between some species [5,6,7,8]. As a result, almost half (48) of the nominal species are believed to be synonymous .
The shape of the cephalic anchors, also known as the “horns” or “antlers”, and their processes, has been traditionally used as the most reliable characteristics for taxonomic identification of Lernaea spp. [9, 10]. However, experimental infection studies put a major question mark over the validity of the prevailing taxonomy of the genus by showing that the anchor exhibits high intraspecific variation [5, 11]. For example, Harding  and Fryer  have shown that the growth and orientation of the anchor are affected by the anatomy of the host. As discussed by Fryer , Yashouv  collected larvae from adults settled on one host and then successfully infected another fish species (host not specified by Fryer) and tadpoles with these larvae; most of the specimens from carp and buffalo fish were assignable to L. cyprinacea, but specimens from Gambusia were morphologically notably different from L. cyprinacea. A subsequent study of Poddubnaya (1973) even further casted the shadow of doubt over the usefulness of the anchor for species identification: when different hosts were infected by larvae from a single batch of L. elegans (a synonym of L. cyprinacea) eggs, adult parasites exhibited different anchor characteristics, some of which corresponded to other described species of Lernaea . Finally, a recent study , used 18S and 28S gene fragments to identify four Lernaea specimens, which could be assigned to different species on the basis of their morphology (the authors did not indicate which species), and found that the specimens shared a similar genotype, so all were identified as L. cyprinacea. These results clearly demonstrate that the anchor characteristics are often merely structural adaptations of a single parasitic specimen to a different host species, and therefore have very limited reliability for taxonomic purposes.
Although molecular data (partial sequences of 18S and 28S rRNA) have been used in the identification of Lernaea species [10, 11, 14,15,16], their availability remains extremely limited. To illustrate this, in December 2018 there were only 21 sequences for just two genes (18S and 28S), and a mitochondrial genome sequence, all belonging to a single species, L. cyprinacea, available in the GenBank database (Additional file 1: Table S1).
The reliability of the two single gene-based molecular markers (18S and 28S) that have been traditionally used in the identification and phylogenetic studies of Lernaea species remains questionable, or at least unconfirmed, largely due to a limited number of studies and limited amount of molecular data publicly available [10, 16]. The DNA barcode marker, cox1, has not been used in the identification or phylogenetic studies. Therefore, we can conclude that due to the shortage of molecular data and the shortcomings of the anchor as a tool for taxonomy and identification, scientists currently do not have a single sufficiently reliable tool for the identification of Lernaea species at disposal.
Lernaea cruciata Lesueur, 1824, first reported from the body surface of the rock bass Cichla aenea , has been subsequently reported from more than a dozen fish species . The anchor of this species is X-shaped, with four simple and short arms approximately equal in size. Lernaea cyprinacea is the most widely distributed species of the genus , which also exhibits very low host specificity, infecting a wide range of freshwater fishes, as well as some amphibians [19, 20], but its anchoring apparatus is much more complex. It also has two pairs of arms , but the dorsal pair is larger than the ventral pair, it branches out into pairs of processes at the tip; hence, some studies refer to it as the T- or Y-shaped dorsal ramified pair, or as “antlers” [5, 21,22,23]. For instance, experimental infection trials on Gambusia sp. infected with L. cyprinacea showed a significant change in the shape of the anchor. In this case, the anchor of L. cyprinacea specimens developed into the adults which were similar to those described for L. gobioptera  and/or L. cruciata [17, 25], and not the typical form of L. cyprinacea. However, despite the fact that the evidence for its extreme variability was presented several decades ago, perhaps due to the absence of other morphological features that could be used for this purpose, the shape of the anchor remains commonly used as the most useful morphological trait, for the identification of Lernaea species [10, 11, 16, 26,27,28,29]. The issue of reliability of this trait has not been revisited in decades, and it has never been studied using a combination of experimental infection and molecular tools. Our working hypothesis was that L. cyprinacea and L. cruciata are synonyms, but the shape of the anchor is prone to host-induced morphological variability, which results in taxonomic misidentification. This would also prove that this morphological trait is not a reliable tool for the identification of (at least these two) Lernaea species. In a preliminary survey, specimens morphologically (anchor) corresponding to L. cyprinacea were found on four cyprinid fish species, whereas specimens morphologically corresponding to L. cruciata were found only on Gambusia holbrooki (the eastern mosquitofish, referred to as mosquitofish henceforth). To test our hypothesis, and the taxonomic validity of L. cruciata, we infected different host species with these specimens, and used molecular data to corroborate the identity of specimens collected from different hosts.
Sample collection and identification of Lernaea cyprinacea and L. cruciata from wild fishes
In 2016 and 2017, we collected 517 fish specimens belonging to 5 species (Carassius auratus, Cyprinus carpio, Erythroculter ilishaeformis, Gambusia holbrooki, Squaliobarbus curriculus) from six locations, comprising wild, farmed and captive fish populations (Table 1). Except for the mosquitofish, which were captured using dip nets with 5 mm stretched mesh size, fishes were captured using trawl nets. Immediately after sampling, the fish specimens were visually examined (oral cavity, branchial cavities, head, fins and skin) for the presence of copepod parasites. Uninfected fish were returned to the water unharmed, and fishes infected with copepod parasites were placed in aerated tanks and transported to the laboratory alive. Parasites were collected from hosts using forceps and a dissecting needle, and preserved in 70% ethanol for detailed identification by light microscopy. Species identification was conducted according to the available literature [16, 17, 25, 30, 31].
Experimental infections of goldfish and mosquitofish with L. cyprinacea
Goldfish (Carassius auratus) specimens with a mean body weight of 6.24 ± 0.78 g were obtained from a local fish farm in Wuhan City, China. The population had no history of infection with L. cyprinacea (i.e. no previous reports of visible signs of infection). Mosquitofish were collected from the Donghu Lake (Table 1), and (apparently) gravid female specimens were selected to be taken to the laboratory and used to produce broodstock via induced breeding (mosquitofish is ovoviviparous). The fingerlings used for the experimental infection were reared under controlled nursery conditions; they were apparently healthy and of same age. The two experimental populations were kept in two separate, aerated tanks (30% water exchanged fortnightly) for 30 days: goldfish (n = 125) in a 500 litre tank, and mosquitofish (n =73) in a 30 litre tank. To remove all ectoparasites, all experimental fish were treated with three consecutive baths in 1:10,000 formalin solution for 12 h at 48 h intervals . Treated fish were then examined under a dissection microscope to ensure that they are free of Lernaea parasites. A subset of specimens of both species (n = 30 for each species) was then selected for the experimental inspection. These specimens were moved to small (5 litre) plastic aquaria (1 specimen/aquarium) filled with static dechlorinated tap water and equipped with aerators and acclimatized and quarantined for at least seven days before being used in experiments. The light:dark cycle was 12:12 h, water pH ranged between 7.4–8.0, and temperature ranged between 24–28 °C. Fish were fed twice daily with commercial pelleted feed, at 2% of the estimated total fish biomass. Egg-sacs of a Lernaea specimen, obtained from Erythroculter ilishaeformis obtained from the Yangluo farm (Table 1), morphologically and genetically (cox1, 18S and 28S; GenBank accession numbers MH982220, MH982197 and MH982208, respectively) identified as L. cyprinacea, were removed and cultured under laboratory conditions as previously described . This laboratory stock of L. cyprinacea (offspring of a single parasite) used for all experimental infestations was maintained on a population (n = 20) of laboratory-reared goldfish. Copepodids-I (the infective stages) of L. cyprinacea were obtained by culturing parasite eggs extracted from the egg-sacs of gravid female parasites under controlled laboratory conditions. Viability of the hatched copepodid-I specimens was inferred by counting the active swimming specimens in diluted stock suspensions under the microscope (4× magnification). Only active parasites were taken into consideration for calculating the number of copepodid-I specimens used to infect each group.
Fifteen goldfish and 15 mosquitofish specimens were isolated and exposed to infective stages of L. cyprinacea at the rates of 70 and 30 copepodids, respectively. Each fish specimen remained isolated in a container with 5 litres of water to prevent the other fish from removing their parasites, and examined daily using a stereomicroscope to determine the success of the infection. Two parallel control groups (same species, same size, n =15), unexposed to parasites, were maintained in an identical environment (isolated in aquaria).
Water exchange and removal of faecal material were not carried out during the initial seven days post-infection to prevent accidental removal of free-living infective stages. After this initial period, faecal and uneaten feed material settled at the bottom of the pools were removed by siphoning with plastic tubes regularly. Regular water exchange (25%) was performed on alternate days. Observations were conducted daily from the onset of infection. Parasites with egg-sacs were removed when they were fully developed. Lernaea parasites were carefully detached using forceps and dissecting needle, and examined immediately. Species identification was conducted as described above.
Morphological and statistical analysis
All of the examined parasite specimens were photographed, and subjected to detailed morphological and morphometric analyses. The specimens were visualized and photographed under a stereomicroscope (Leica S8AP0, Wetzlar, Germany), equipped with a SPOT Insight 2.0 Mp digital camera (Sterling Heights, MI, USA). The length of the whole body, arms and processes, and anchor width were measured in mm (Fig. 1), with allowance being made for the bends and curves of the body, as described by Robinson & Avenant-Oldewage .
Principal components analysis (PCA) was performed to determine the difference between different groups (i.e. L. cyprinacea and L. cruciata samples collected in the wild and from experimental infections). All physical variables were standardized to zero mean and unit variance to make them dimensionless. Results were considered significant at the 95% level (P < 0.05). Ordination and multivariate analysis of variance (MANOVA) analyses were performed using the vegan package in R version 3.6.1 [34, 35].
Molecular data amplification and sequencing
Molecular analysis was performed for 7 L. cyprinacea and 5 L. cruciata specimens collected from host fish in the wild and from those collected during the laboratory infection experiment. Total genomic DNA was extracted from the posterior part of the parasite body (to preserve the anchor intact) using TIANamp Micro DNA Kit (Tiangen Biotech, Beijing, China), according to the manufacturer’s instructions. DNA was stored at − 20 °C for subsequent molecular analysis. Three molecular markers were amplified by PCR and sequenced: a fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1), a partial sequence of the 18S rRNA gene, and a partial sequence of the 28S rRNA gene. cox1 was amplified using newly designed primers (forward: 5′-TAG TTG GAA TTT GGG CTG GC-3v and reverse: 5′-ATT AGG GGC CTT GTT GGG AAG-3′). The PCR reaction mix (20 µl) was comprised of 1 µl genomic DNA, 0.6 µl of both primers, 7.4 µl ddH2O, 10 µl 2× PCR buffer (Mg2+, dNTP plus, Takara, Dalian, China), 0.4 µl rTaq polymerase (250 U, Takara); conditions were: initial denaturation at 98 °C for 2 min, 40 cycles of 98 °C for 10 s, 50 °C for 15 s, 68 °C for 1 min/kb, and the final extension at 68 °C for 10 min. Primers and PCR parameters for 18S and 28S were previously published . Amplicons were subjected to electrophoresis on a 1.2% agarose gel stained with GoldView dye (Solarbio, Beijing, China), and then purified and sequenced by the Sangon Co. (Shanghai, China). The 36 newly generated sequences were compared with similar sequences available in the GenBank database using the Basic Local Alignment Search Tool (BLAST)  to confirm their identity and deposited in GenBank under the accession numbers MH982192-MH982227 (Additional file 2: Table S2).
Comparative and phylogenetic analyses
Sequences obtained in this study (cox1, 18S, 28S genes of the nominal L. cyprinacea and L. cruciata specimens) were aligned with all Lernaea spp. (all belonging to L. cyprinacea) available on GenBank (Additional file 1: Table S1) using MAFFT 7.149 . For comparative analyses, 18S and 28S sequences obtained from GenBank that did not exhibit a full overlap with the sequences obtained in this study were removed from the alignment. Multiple alignment, pairwise identity and the number of variable sites among the sequences of cox1, 18S and 28S were conducted using Geneious 8.1.3 . The number of haplotypes for each molecular marker was calculated using DnaSP v6.11.01 . Sites with gaps and missing data were not considered.
For phylogenetic analyses, we used all available Lernaea sequences and added outgroups. As the only cox1 gene sequence belonging to the family Lernaeidae available on GenBank was that from the mitogenome of L. cyprinacea (KM235194), for the outgroup we used two sequences for Sinergasilus polycolpus belonging to another family of the order Cyclopoida (Ergasilidae) (Additional file 3: Table S3). For the 18S data, we selected all available homologues belonging to the closest-related  genus in the family Lernaeidae, Lamproglena, to root the tree (Additional file 3: Table S3). For the 28S dataset, five sequences belonging to two Lamproglena spp. were chosen as the outgroup (Additional file 3: Table S3). For the fourth dataset, to maximize the phylogenetic resolution, we concatenated these three genes into a single alignment (cox1, 18S and 28S). Due to lack of suitable data, we did not include an outgroup in this analysis. Phylogenetic analyses were conducted using four datasets (cox1, 18S, 28S and concatenated) and two methods: maximum likelihood (ML) and Bayesian inference (BI). ML was conducted using IQTree v1.6.3 , and BI using MrBayes 3.2 . We used the akaikeʼs information criterion (AIC), implemented in ModelFinder  in IQ-Tree, to select the best-fit evolutionary models for each dataset (Additional file 4: Table S4).
Morphology of Lernaea spp. collected in the wild
After carefully inspecting 517 fish specimens belonging to five species (Table 1) collected in the wild, and discarding all parasite specimens which could not be identified with confidence (not fully developed or with incomplete anchor arms, usually damaged while removing from the body surface of fish), we morphologically identified 25 specimens as L. cyprinacea and 16 as L. cruciata. Whereas all L. cruciata specimens were collected from the mosquitofish, L. cyprinacea specimens were collected from hosts belonging to four different species: Cyprinus carpio, C. auratus, Squaliobarbus curriculus and Culter alburnus.
Experimental infection results and comparative morphology of L. cyprinacea and L. cruciata
Specimens found on fish after the experimental infection were regarded as mature (fully developed) when they were ovigerous. In the C. auratus infection experiment, 15 Lernaea specimens were collected from six fish specimens and morphologically identified as L. cyprinacea. In the mosquitofish infection experiment, only 5 adult Lernaea specimens were collected from 5 fish and all morphologically identified as L. cruciata. One of these five parasites exhibited somewhat mixed morphological traits, with one of its arms exhibiting minuscule processes, but it predominantly exhibited a morphology corresponding to L. cruciata (Fig. 2).
Among the parasites obtained in experimental infection, the average body length and anchor width of L. cyprinacea (8.86 ± 1.71 and 2.31 ± 1.37 mm, respectively) were greater than those of L. cruciata specimens (7.10 ± 1.00 and 1.19 ± 0.29 mm, respectively) (Table 2). The average ventral anchor arms of L. cyprinacea (0.55 ± 0.21 mm) were also longer than those of L. cruciata (0.49 ± 0.10 mm). However, none of these differences between samples of L. cyprinacea and L. cruciata were statistically significant (t-test: body length, t(18) = 2.159, P = 0.45, anchor width, t(18) = 1.78, P = 0.92, ventral arms, t(18) = 0.654, P = 0.522).
Multivariate morphometric comparison of the wild and experimental L. cyprinacea and L. cruciata specimens
PCA (Fig. 3) revealed that PC1, mainly determined by the width of the anchor and the length of the ventral arms, explained 50.6% of the variation contained within all morphological characteristics. PC2, mainly associated with the length of the body and ventral arms, explained 38.1% of the total variation. Therefore, PC1 and PC2 combined explained almost 90% (88.7%) of the variance contained within the dataset. Apart from a single outlier (Fig. 3), L. cruciata specimens collected in the wild and from the laboratory infection clustered together. The distribution of L. cyprinacea specimens (both wild and laboratory-reared) was scattered, but several laboratory-reared specimens overlapped with the L. cruciata cluster.
Partial cox1, 18S and 28S genes (1241, 1389 and 706 bp, respectively) all exhibited very high identity (BLAST) values to the corresponding available L. cyprinacea sequences in the GenBank database: 99%, 99–100% and 99–100%, respectively (Additional file 5: Table S5). We aligned our sequences with all fully overlapping L. cyprinacea homologs available in the GenBank.
All cox1 sequences obtained in this study and one cox1 sequence of L. cyprinacea downloaded from GenBank were used to compare the sequence similarity. High pairwise identity (98.2–100%) was found among the cox1 sequences of L. cyprinacea (Lcy) and putative L. cruciata (Lcr) (Table 3). The highest number of variable sites among the sequences was 22 (between Lcy2 and Lcy6, Lcy2 and Lcy7). There were 9 haplotypes and 45 variable sites among these sequences.
Twelve 18S sequences of L. cyprinacea (seven sequenced in this study and five from GenBank) and five 18S sequences of L. cruciata (sequenced in this study) also exhibited very high pairwise identity (99.8–100%) (Table 4). There were 7 haplotypes and only 7 variable sites. It is also noteworthy that all of the newly generated 18S sequences of L. cyprinacea, two 18S sequences of L. cyprinacea downloaded from GenBank (DQ107555, KP235363), and two 18S sequences of L. cruciata, shared the same haplotype.
Twelve 28S sequences of L. cyprinacea (seven sequenced in this study and four downloaded from GenBank) and five 28S sequences of L. cruciata (all sequenced in this study) were compared (Table 5). Due to the poorly conserved 3′-end of the LcyGB KM281817 sequence (No. 17), which might be a sequencing artefact, pairwise identity between this sequence and other sequences (98.9–99.1%) was marginally lower than pairwise identity among all remaining sequences (99.4–100%). Variable sites between KM281817 and other sequences ranged from 6 to 8, while variable sites among the remaining sequences ranged from 0 to 4. There were 10 haplotypes and 15 variable sites among all sequences. It is also noteworthy that three 28S sequences of L. cyprinacea (Lcy1, Lcy3, Lcy6) and all 28S sequences of L. cruciata shared the same haplotype.
Despite minor variations in the topology between different algorithms (ML and BI) and datasets (18S, 28S, cox1, concatenated genes), all eight obtained phylograms (Figs. 4, 5, 6, 7) produced identical results in two important aspects: monophyletic Lernaea clade and highly admixed (paraphyletic) L. cyprinacea and L. cruciata sequences. Sequences from the two morphospecies did not exhibit clear clustering according to morphotypes or the geographical origin. In contrast, they appeared to be randomly distributed within the cluster.
To date, species of Lernaea are still mainly differentiated on the basis of a combination of morphological characters, primarily the shape of the anchors, which is believed to be the most reliable characteristic for identification. However, it has been shown more than half a century ago that anchor morphology to some extent depends on the thickness of the tissue to which the parasite is attached [6, 7]. The age of the parasite, host species, and the site of infestation, can all affect the shape of anchors in Lernaea spp. [5, 6]. Most notably, after infecting different fish hosts with larvae from a single batch of eggs of L. cyprinacea, Yashouv  observed specimens with different anchor shapes, some of which corresponded morphologically to other described Lernaea species. However, Yashouv  did not have molecular tools at disposal to corroborate this observation. Herein, we compared morphological characteristics of specimens collected in the wild, conducted a laboratory infection experiment wherein different fish hosts were infected with the offspring of a single specimen of L. cyprinacea, and sequenced nuclear and mitochondrial genes. Our results showed that the offspring of L. cyprinacea collected from mosquitofish could be morphologically identified as L. cruciata, while genetically they corresponded to the parental species, L. cyprinacea. The observed morphological variation must therefore be host-induced phenotypic variation, affecting the diagnostic morphological features. Despite the lack of sequences for other Lernaea species on GenBank, very high identity values among the gene fragments sequenced in this study indicate that specimens from our study genetically correspond to L. cyprinacea. Therefore, our results show that the shape of the anchors cannot be considered as a reliable character for the identification of Lernaea spp. As there may exist some reproductive barriers between populations parasitizing on different host species, or putatively a genetic preference for a certain host species, further studies should target samples from other populations and a broad range of hosts, and attempt to precisely (if possible) identify the genetic boundaries of this species and the amount of gene flow between its different morphotypes. Importantly, the use of molecular data should be considered a prerequisite for any study that requires taxonomic identification of species belonging to this genus.
In terms of underlying reasons for this host-induced morphological variation, Fryer  suggested that the anchor of L. cyprinacea is reduced, and that branching of the dorsal anchors (processes) tends to be suppressed, when it parasitizes on a small host, where a limited volume of tissue would limit their full development. However, we found that some of these specimens with suppressed dorsal anchors were attached on the abdomen and the base of the dorsal fin of the mosquitofish, where they appeared to have sufficient space to develop their anchors. Furthermore, L. cruciata was also reported on some lager-sized fish species, such as the rock bass Cichla aenea  and the white bass Morone chrysops . We can therefore only speculate that a number of host characteristics may affect the morphology of the parasite; it could be the size of the host, the structure of the host’s muscle and the hardness of the tissue, nutritional substances that the parasite obtains from the host, skin (or mucus) chemistry, etc.
Mosquitofish (Gambusia affinis and G. holbrooki) are native to the southeastern USA, but now occur on every continent except Antarctica due to introduction programmes implemented since the early 20th century . Lernaea spp. infestations of mosquitofish had been reported from Bangladesh , Turkey [43, 45] and China . Lernaea cyprinacea specimens from Bangladesh  were found to exhibit a typical L. cyprinacea morphology, whereas specimens from Turkey  morphologically better corresponded to L. cruciata. Razavi et al.  investigated an infection on the Farsi toothcarp (Aphanius farsicus) and found that the collected Lernaea specimens morphologically best corresponded to L. cruciata, which led them to speculate that L. cruciata may have been translocated into the Maharlou Lake Basin by mosquitofish. Lernaea gobioptera, another species first described on the body surface of Gobiopterus macrolepis and mosquitofish in the delta of the Pearl River in China , is morphologically very similar to L. cruciata [24, 25, 46]. Both species share a typical trait: a single pair of branched holdfasts, and dorsal and ventral arms of the anchor of about equal size, which makes it look X-shaped. Therefore, we suggest that L. cyprinacea tends to develop a specific shape of anchors when it infects specific fish taxa, e.g. Gobiopterus sp. and mosquitofish, which results in misidentification of these specimens as different species.
In the experiment by Yashouv , specimens from mosquitofish infected by Lernaea larvae were found to be quite different morphologically from L. cyprinacea, so Fryer  suggested that these specimens look like a new species. Having reviewed the images in Fryerʼs paper, we argue that they resemble L. cruciata. Unfortunately, we could not access the original study of Yashouv , so we had to rely on Fryer’s review of it . Regardless of this limitation, our experimental infection results correspond to those by Yashouv, which indicates that the observed anchor shape change of L. cyprinacea specimens when parasitizing on Gambusia was not an outlier. However, the specimen which had intermediate morphological characteristics between L. cyprinacea and L. cruciata in our study, and a L. cyprinacea specimen obtained from mosquitofish in the field , suggest that L. cyprinacea specimens infecting mosquitofish may exhibit three morphotypes: typical morphological features of L. cruciata; intermediate morphological characteristics; and typical morphological features of L. cyprinacea.
Although morphology still plays an important role in species description, identification and taxonomy , there is a growing amount of evidence, from a broad range of taxa, that morphology alone often does not provide adequate taxonomic resolution and that it often leads to erroneous conclusions in some taxa [48,49,50,51,52,53,54,55]. Intriguingly, this appears to be particularly widespread in parasitic taxa [54,55,56,57,58,59,60,61,62], and taxonomic artefacts caused by host-induced morphological variability have been reported in several other parasitic taxa. For example, a molecular study revealed that Caryophyllaeus laticeps tapeworms that parasitize breams are merely a morphotype of Caryophyllaeus brachycollis Janiszewska, 1953, which parasitizes other cyprinid fishes . Also, morphological traits of Isthmiophora melis are highly variable and host-dependent, and without the support of molecular data they could easily lead to a misidentification of several apparently distinct species, or even genera . Although molecular studies also have their limitations , owing to the wide use of molecular techniques in taxonomic studies, the number of valid taxa has changed in most major taxonomic groups during the last few decades. We expect that a major re-evaluation based on molecular data would also result in a notably reduced number of valid species in the genus Lernaea.
Although molecular data have been used in the identification of Lernaea species [10, 11, 14,15,16, 65, 66], previous molecular studies appear to have focused exclusively on L. cyprinacea, which is therefore the only species for which there are molecular data currently (October 2018) available in public databases. Although these also include a transcriptome , and the complete mitochondrial genome , identification of L. cyprinacea based on molecular data has been focused principally on partial sequences of two nuclear rRNA genes: 18S and 28S [11, 14,15,16]. This unavailability of molecular data for the remaining Lernaea species presents a major obstacle to their application as a tool for species identification . Although our findings indicate that L. cruciata is a synonym of L. cyprinacea, due to this shortage of molecular data we cannot make this claim with confidence. We therefore urge scientists to sequence relevant genes of other Lernaea morphospecies and re-examine the status of the species currently recognized in this genus.
The results of our experimental infections indicate that L. cyprinacea sometimes exhibits different morphological features when parasitizing on different hosts. Considering its wide host range, we suspect that this host-specific morphology of conspecific parasites has resulted in numerous taxonomic artefacts, i.e. misidentifications of morphotypes as new species, and that many of the described species of Lernaea are actually one and the same species, L. cyprinacea. This hypothesis should be tested and validated using both molecular and morphological data, as well as experimental laboratory infections.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files. The newly generated nucleotide sequences were submitted to the GenBank database under the accession numbers MH982192–MH982227.
cytochrome c oxidase subunit 1 gene
akaikeʼs information criterion
Hemaprasanth KP, Raghavendra A, Singh R, Sridhar N, Raghunath MR. Efficacy of doramectin against natural and experimental infections of Lernaea cyprinacea in carps. Vet Parasitol. 2008;156:261–9.
Tamuli KK, Shanbhogue SL. Incidence and intensity of anchor worm Lernaea bhadraensis infection on cultivated carps. Environ Ecol. 1996;282:8.
Singh R, Raghavendra A, Sridhar N, Raghunath MR, Eknath AE. Comparative susceptibility of carp fingerlings to Lernaea cyprinacea infection. Vet Parasitol. 2011;178:156–62.
WoRMS: World Register of Marine Species. 2018. http://marinespecies.org/aphia.php?p=taxdetails&id=149776. Accessed 6 Oct 2018.
Amin OM. On the crustacean ectoparasites of fishes from Southeast Wisconsin. Trans Am Microsc Soc. 1981;100:142–50.
Harding JP. On some species of Lernaea, crustacea, copepoda, parasites of freshwater fish. Bull Brit Mus (Nat Hist). 1950;1:1–27.
Fryer G. Variation and systematic problems in a group of lernaeid copepods. Crustaceana. 1961;2:275–85.
Kularatne M, Shariff M, Subasinghe RP. Comparison of larval morphometrics of Lernaea minuta, a copepod parasite of Puntius gonionotus from Malaysia, with those of L. cyprinacea and L. polymorpha. Crustaceana. 1994;67:288–95.
Shariff M, Sommerville C. Morphometrics of the larval stages of Lernaea polymorpha Yu and Lernaea cyprinacea Linnaeus (Copepoda). Crustaceana. 1989;57:136–9.
Pallavi B, Shankar KM, Abhiman PB, Ahmed I. Molecular identification of the fish parasite Lernaea. Indian J Fish. 2017;64:76–82.
Soares IA, Salinas V, Ponti OD, Mancini MA, Luque JL. First molecular data for Lernaea cyprinacea (Copepoda: Cyclopoida) infesting Odontesthes bonariensis, a commercially important freshwater fish in Argentina. Rev Bras Parasitol Vet. 2018;27:106–9.
Yashouv A. On the biology of Lernaea in fish ponds. Bamidgeh. 1959;11:80–9.
Poddubnaya A. Variability and specificity of Lernaea parasite on pond fishes. Trudy Vses Nauch Inst Prud Ryb Khoz (in Russian). 1973;22:159–73.
Song Y, Wang GT, Yao WJ, Gao Q, Nie P. Phylogeny of freshwater parasitic copepods in the Ergasilidae (Copepoda: Poecilostomatoida) based on 18S and 28S rDNA sequences. Parasitol Res. 2008;102:299–306.
Stavrescu-Bedivan MM, Popa OP, Popa LO. Infestation of Lernaea cyprinacea (Copepoda: Lernaeidae) in two invasive fish species in Romania, Lepomis gibbosus and Pseudorasbora parva. Knowl Manag Aquat Ecosyst. 2014;414:12.
McCredden M. Anchors away: the susceptibility and response to infection between native and co-introduced fishes to the alien anchor worm Lernaea cyprinacea. PhD thesis, Murdoch University, Perth; 2016.
Lesueur CA. On three new species of parasitic vermes belonging to the Linnaean genus Lernaea. J Acad Nat Sci Phila. 1824;3:286–93.
Joy JE, Jones LP. Observations on the inflammatory response within the dermis of a white bass, Morone chrysops (Rafinesque), infected with Lernaea cruciata (Copepoda: Caligidea). J Fish Biol. 1973;5:21–3.
Piasecki W, Goodwin AE, Eiras JC, Nowak BF. Importance of Copepoda in freshwater aquaculture. Zool Stud. 2004;43:193–205.
Salinas ZA, Biolé FG, Grenat PR, Pollo FE, Sala NE, Martino AL. First report of Lernaea cyprinacea (Copepoda: Lernaeidae) in tadpoles and newly-metamorphosed frogs in wild populations of Lithobates catesbeianus (Anura: Ranidae) in Argentina. Phyllomedusa. 2016;15:43–50.
Grabda J. Life cycle and morphogenesis of Lernaea cyprinacea. Acta Parasitol Pol. 1963;11:169–98.
Kabata Z. Parasites and diseases of fish cultured in the tropics. Abingdon: Taylor & Francis Ltd.; 1985.
Shariff M, Kabata Z, Sommerville C. Host susceptibility to Lernaea cyprinacea L. and its treatment in a large aquarium system. J Fish Biol. 1986;9:393–401.
Chen R. A report on the Lernaea parasites of freshwater fishes in Pearl River Delta. Jinan Liyi Xuebao. 1981;1:64–7.
USGS Parasitic Copepods database. https://www.glsc.usgs.gov/greatlakescopepods/Image.php?PATH=Parasite/Introduction/&SRC=42Lcrubody.jpg. Accessed 15 Aug 2018.
Shariff M. Studies on the biology, host-parasite interactions and distribution of Lernaea spp. in West Malaysia. PhD thesis. Stirling: University of Stirling; 1985.
Ahnelt H, Konecny R, Gabriel A, Bauer A, Pompei L, Lorenzoni M, et al. First report of the parasitic copepod Lernaea cyprinacea (Copepoda: Lernaeidae) on gobioid fishes (Teleostei: Gobonellidae) in southern Europe. Knowl Manag Aquat Ecosyst. 2018;419:34.
Robinson J, Avenant-Oldewage A. Aspects of the morphology of the parasitic copepod Lernaea cyprinacea Linnaeus, 1758 and notes on its distribution in Africa. Crustaceana. 1996;69:610–26.
El-Mansy AI. On the occurrence of adult females of Lernaea species (Crustacea: Copepoda) parasitic on goldfish Carassius auratus (Linnaeus) in some commercial aquaria in Egypt. Egypt J Aquat Biol Fish. 2009;13:7–36.
Kuang P. Parasitic Crustacea of freshwater fishes. Beijing: Science Press; 1991.
Yin W, Lin M, Xu G, Chen Y, Kuang P, Zhu X. Studies on the lernaeosos (Lernaea, Copepoda parasite) of fresh-water fishes of China. Shui Sheng Sheng Wu Xue Ji Kan. 1963;11:48–117.
Zhou S, Zou H, Wu SG, Wang GT, Marcogliese DJ, Li WX. Effects of goldfish (Carassius auratus) population size and body condition on the transmission of Gyrodactylus kobayashii (Monogenea). Parasitology. 2017;144:1221–8.
Shields RJ. Procedures for the laboratory rearing of Lernaea cyprinacea L. (Copepoda). Crustaceana. 1978;35:259–64.
Anderson M. A new method for non-parametric multivariate analysis of variance. Aus Ecol. 2001;26:32–46.
Mcardle BH, Anderson MJ. Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology. 2001;82:290–7.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.
Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2014;32:268–74.
Ronquist F, Teslenko M, Mark PVD, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Kalyaanamoorthy S, Minh BQ, Wong TK, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587.
Innal D, Avenant-Oldewage A, Dogangil B, Stavrescu-Bedivan MM, Ozmen O, Mavruk S. Susceptibility of endemic and non-indigenous fish to Lernaea cyprinacea (Copepoda: Lernaeidae): a case study from Duger Spring Creek (Burdur-Turkey). Bull Eur Ass Fish Pathol. 2017;37:100–9.
Hossain MMM, Ferdoushi J, Haque A. Biology of anchor worms (Lernaea cyprinacea). J Entomol Zool Stud. 2018;6:910–7.
Innal D, Avenant-Oldewage A. Occurrence of Lernaea cyprinacea on mosquito fish (Gambusia affinis) from Kundu Estuary (Antalya-Turkey). Bull Eur Ass Fish Pathol. 2012;32:140–7.
Razavi S, Alavi-Yeganeh M, Peyghan R. Report of infection to Lernaea cruciata (Le Sueur, 1824) in Farsi toothcarp Aphanius farsicus. Fish Sci Technol. 2017;5:128–33.
Nadler SA, León DE. Integrating molecular and morphological approaches for characterizing parasite cryptic species: implications for parasitology. Parasitology. 2011;138:1688–709.
Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci USA. 2004;101:14812–7.
Urantowka AD, Strzała T, Grabowski KA. Complete mitochondrial genome of endangered Maroon-fronted Parrot (Rhynchopsitta terrisi)—conspecific relation of the species with thick-billed Parrot (Rhynchopsitta pachyrhyncha). Mitochondrial DNA. 2014;25:424–6.
Táncsics A, Benedek T, Farkas M, Máthé I, Márialigeti K, Szoboszlay S, et al. Sequence analysis of 16S rRNA, gyrB and catA genes and DNA-DNA hybridization reveal that Rhodococcus jialingiae is a later synonym of Rhodococcus qingshengii. Int J Syst Evol Microbiol. 2014;64:298–301.
Peng HS, Yuan QJ, Li QQ, Huang LQ. Molecular systematics of genus Atractylodes (Compositae, Cardueae): evidence from internal transcribed spacer (ITS) and trnL-F sequences. Int J Mol Sci. 2012;13:14623–33.
Gómez F, Qiu D, Lin S. The synonymy of the toxic dinoflagellates Prorocentrum mexicanum and P. rhathymum and the description of P. steidingerae sp. nov. (Prorocentrales, Dinophyceae). J Eukaryot Microbiol. 2017;64:668–77.
Li Y, Zhou X, Feng G, Hu H, Niu L, Hebert PD, et al. COI and ITS2 sequences delimit species, reveal cryptic taxa and host specificity of fig-associated Sycophila (Hymenoptera, Eurytomidae). Mol Ecol Resour. 2010;10:31–40.
Hernández-Orts JS, Scholz T, Brabec J, Kuzmina T, Kuchta R. High morphological plasticity and global geographical distribution of the Pacific broad tapeworm Adenocephalus pacificus (syn. Diphyllobothrium pacificum): molecular and morphological survey. Acta Trop. 2015;149:168–78.
Hildebrand J, Adamczyk M, Laskowski Z, Zaleśny G. Host-dependent morphology of Isthmiophora melis (Schrank, 1788) Luhe, 1909 (Digenea, Echinostomatinae)—morphological variation vs molecular stability. Parasit Vectors. 2015;8:481.
Poulin R, Morand S. The diversity of parasites. Quart Rev Biol. 2000;75:277–93.
Brazenor AK, Saunders RJ, Miller TL, Hutson KS. Morphological variation in the cosmopolitan fish parasite Neobenedenia girellae (Capsalidae: Monogenea). Int J Parasitol. 2018;48:125–34.
Miller TL, Cribb TH. Dramatic phenotypic plasticity within species of Siphomutabilus n. g. (Digenea: Cryptogonimidae) from Indo-Pacific caesionines (Perciformes: Lutjanidae). Syst Parasitol. 2013;86:101–12.
Hogans W. Review of Pennella Oken, 1816 (Copepoda: Pennellidae) with a description of Pennella benzi sp. nov., a parasite of escolar, Lepidocybium flavobrunneum (Pisces) in the northwest Atlantic Ocean. Zootaxa. 2017;4244:1–38.
Dantas-Torres F. Species concepts: what about ticks? Trends Parasitol. 2018;34:1017–26.
Perkins EM, Donnellan SC, Bertozzi T, Chisholm LA, Whittington ID. Looks can deceive: molecular phylogeny of a family of flatworm ectoparasites (Monogenea: Capsalidae) does not reflect current morphological classification. Mol Phylogen Evol. 2009;52:705–14.
Zhang D, Zou H, Wu SG, Li M, Jakovlić I, Zhang J, et al. Sequencing of the complete mitochondrial genome of a fish-parasitic flatworm Paratetraonchoides inermis (Platyhelminthes: Monogenea): tRNA gene arrangement reshuffling and implications for phylogeny. Parasit Vectors. 2017;10:462.
Barčák D, Oros M, Hanzelová V, Scholz T. Phenotypic plasticity in Caryophyllaeus brachycollis Janiszewska, 1953 (Cestoda: Caryophyllidea): does fish host play a role? Syst Parasitol. 2014;88:153–66.
Degnan JH, Rosenberg NA. Discordance of species trees with their most likely gene trees. PLoS Genet. 2006;2:e68.
Pallavi B, Shankar KM, Abhiman PB, Iqlas A. Identification of putative genes involved in parasitism in the anchor worm, Lernaea cyprinacea by de novo transcriptome analysis. Exp Parasitol. 2015;153:191–7.
Su YB, Wang LX, Kong SC, Chen L, Fang R. Complete mitochondrial genome of Lernaea cyprinacea (Copepoda: Cyclopoida). Mito DNA. 2016;27:1503–4.
The authors would like to thank Dr Rui Song and Min Xie for their assistance in field sampling.
This study was funded by the Earmarked Fund for China Agriculture Research System (CARS-45-15), the National Natural Science Foundation of China (31872604, 31572658), and the Major Scientific and Technological Innovation Project of Hubei Province (2015ABA045).
Ethics approval and consent to participate
All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
All molecular sequences belonging to the genus Lernaea currently (December 2018) available on GenBank. Mt.gen. is a complete mitochondrial genome, sequence size is given in bp, C column indicates whether the sequence was used for our comparative analyses, and Ref. is the associated reference.
Sequences used for the comparative analyses in this study.
Sequences used for the phylogenetic analyses in this study.
Best-fit models for four datasets (cox1, 18S, 28S and concatenated), selected based on the Akaikeʼs information criterion using ModelFinder software.
The identity values between partial cox1, 18S and 28S genes obtained in this study and Lernaea cyprinacea homologues available on GenBank.
About this article
Cite this article
Hua, C.J., Zhang, D., Zou, H. et al. Morphology is not a reliable taxonomic tool for the genus Lernaea: molecular data and experimental infection reveal that L. cyprinacea and L. cruciata are conspecific. Parasites Vectors 12, 579 (2019). https://doi.org/10.1186/s13071-019-3831-y
- Parasitic copepods
- Phenotypic plasticity
- Molecular taxonomy