Skip to main content

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

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Background

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 [4], 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 [4].

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 [6] and Fryer [7] have shown that the growth and orientation of the anchor are affected by the anatomy of the host. As discussed by Fryer [7], Yashouv [12] 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 [13]. Finally, a recent study [10], 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 [17], has been subsequently reported from more than a dozen fish species [18]. 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 [19], 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 [21], 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 [24] 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.

Methods

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].

Table 1 List of collected Lernaea specimens with collection details

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 [32]. 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 [33]. 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 [28].

Fig. 1
figure1

a A representative Lernaea cyprinacea specimen from goldfish. b A representative Lernaea cruciata specimen from mosquitofish. Abbreviations: aas, anchor arms; aw, anchor width; va, ventral arm; adp, anterior dorsal process; pdp, posterior dorsal process

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 [14]. 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) [36] 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 [37]. 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 [38]. The number of haplotypes for each molecular marker was calculated using DnaSP v6.11.01 [39]. 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 [14] 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 [40], and BI using MrBayes 3.2 [41]. We used the akaikeʼs information criterion (AIC), implemented in ModelFinder [42] in IQ-Tree, to select the best-fit evolutionary models for each dataset (Additional file 4: Table S4).

Results

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).

Fig. 2
figure2

A Lernaea specimen collected from an experimentally infected mosquitofish, and identified as L. cruciata, exhibiting minuscule processes on one of its arms, somewhat resembling those of L. cyprinacea. The arrow highlights the feature somewhat resembling a posterior dorsal process

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).

Table 2 Morphometrics of Lernaea cyprinacea (Lcy) and L. cruciata (Lcr) specimens collected from experimental infections and in the wild

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.

Fig. 3
figure3

Principal components analysis (PCA) ordination showing the difference of morphological features between the two different Lernaea species collected in the wild and from an experimental (laboratory) infection. Abbreviations: Lcy, L. cyprinacea; Lcr, L. cruciate; l, specimens collected from the laboratory infection; w, specimens collected in the wild. Notes: Samples are divided into four groups (1 to 4). For example, 4_Lcr_l indicates group 4, comprising specimens of L. cruciata collected from the laboratory infection

Sequence comparisons

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.

cox1

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.

Table 3 Pairwise identity (below the diagonal) and variable sites (above the diagonal) for the cox1 dataset

18S

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.

Table 4 Pairwise identity (below the diagonal) and variable sites (above the diagonal) for the 18S dataset

28S

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.

Table 5 Pairwise identity (below the diagonal) and variable sites (above the diagonal) for the 28S dataset

Phylogenetic analyses

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, 67) 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.

Fig. 4
figure4

Phylogenetic trees inferred using the 18S dataset and Bayesian inference (BI) (a) and maximum likelihood (ML) (b) methods. Numbers next to nodes indicate bootstrap (ML)/posterior probability (BI) support. Abbreviations: Lcy, Lernaea cyprinacea; Lcr, L. cruciate; w, specimens collected in the wild; l, specimens collected from the laboratory infection

Fig. 5
figure5

Phylogenetic trees inferred using the 28S dataset and Bayesian inference (BI) (a) and maximum likelihood (ML) (b) methods. Numbers next to nodes indicate bootstrap (ML)/posterior probability (BI) support. Abbreviations: Lcy, Lernaea cyprinacea; Lcr, L. cruciate; w, specimens collected in the wild; l, specimens collected from the laboratory infection

Fig. 6
figure6

Phylogenetic trees inferred using the cox1 dataset and Bayesian inference (BI) (a) and maximum likelihood (ML) (b) methods. Numbers next to nodes indicate bootstrap (ML)/posterior probability (BI) support. Abbreviations: Lcy, Lernaea cyprinacea; Lcr, L. cruciate; w, specimens collected in the wild; l, specimens collected from the laboratory infection

Fig. 7
figure7

Phylogenetic trees inferred using the concatenated genes (cox1, 18S and 28S) dataset, and Bayesian inference (BI) (a) and maximum likelihood (ML) (b) methods. Numbers next to nodes indicate bootstrap (ML)/posterior probability (BI) support. Abbreviations: Lcy, Lernaea cyprinacea; Lcr, L. cruciate; w, specimens collected in the wild; l, specimens collected from the laboratory infection

Discussion

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 [12] observed specimens with different anchor shapes, some of which corresponded morphologically to other described Lernaea species. However, Yashouv [12] 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 [7] 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 [17] and the white bass Morone chrysops [18]. 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 [43]. Lernaea spp. infestations of mosquitofish had been reported from Bangladesh [44], Turkey [43, 45] and China [24]. Lernaea cyprinacea specimens from Bangladesh [44] were found to exhibit a typical L. cyprinacea morphology, whereas specimens from Turkey [45] morphologically better corresponded to L. cruciata. Razavi et al. [46] 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 [24], 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 [12], specimens from mosquitofish infected by Lernaea larvae were found to be quite different morphologically from L. cyprinacea, so Fryer [7] 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 [12], so we had to rely on Fryer’s review of it [7]. 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 [44], 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 [47], 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 [63]. 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 [55]. Although molecular studies also have their limitations [64], 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 [65], and the complete mitochondrial genome [66], 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 [16]. 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.

Conclusions

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.

Abbreviations

cox1:

cytochrome c oxidase subunit 1 gene

ML:

maximum likelihood

BI:

Bayesian inference

AIC:

akaikeʼs information criterion

References

  1. 1.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Tamuli KK, Shanbhogue SL. Incidence and intensity of anchor worm Lernaea bhadraensis infection on cultivated carps. Environ Ecol. 1996;282:8.

    Google Scholar 

  3. 3.

    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.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  4. 4.

    WoRMS: World Register of Marine Species. 2018. http://marinespecies.org/aphia.php?p=taxdetails&id=149776. Accessed 6 Oct 2018.

  5. 5.

    Amin OM. On the crustacean ectoparasites of fishes from Southeast Wisconsin. Trans Am Microsc Soc. 1981;100:142–50.

    Article  Google Scholar 

  6. 6.

    Harding JP. On some species of Lernaea, crustacea, copepoda, parasites of freshwater fish. Bull Brit Mus (Nat Hist). 1950;1:1–27.

    Article  Google Scholar 

  7. 7.

    Fryer G. Variation and systematic problems in a group of lernaeid copepods. Crustaceana. 1961;2:275–85.

    Article  Google Scholar 

  8. 8.

    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.

    Article  Google Scholar 

  9. 9.

    Shariff M, Sommerville C. Morphometrics of the larval stages of Lernaea polymorpha Yu and Lernaea cyprinacea Linnaeus (Copepoda). Crustaceana. 1989;57:136–9.

    Google Scholar 

  10. 10.

    Pallavi B, Shankar KM, Abhiman PB, Ahmed I. Molecular identification of the fish parasite Lernaea. Indian J Fish. 2017;64:76–82.

    Article  Google Scholar 

  11. 11.

    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.

    PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Yashouv A. On the biology of Lernaea in fish ponds. Bamidgeh. 1959;11:80–9.

    Google Scholar 

  13. 13.

    Poddubnaya A. Variability and specificity of Lernaea parasite on pond fishes. Trudy Vses Nauch Inst Prud Ryb Khoz (in Russian). 1973;22:159–73.

    Google Scholar 

  14. 14.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    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.

    Article  Google Scholar 

  16. 16.

    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.

  17. 17.

    Lesueur CA. On three new species of parasitic vermes belonging to the Linnaean genus Lernaea. J Acad Nat Sci Phila. 1824;3:286–93.

    Google Scholar 

  18. 18.

    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.

    Article  Google Scholar 

  19. 19.

    Piasecki W, Goodwin AE, Eiras JC, Nowak BF. Importance of Copepoda in freshwater aquaculture. Zool Stud. 2004;43:193–205.

    Google Scholar 

  20. 20.

    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.

    Article  Google Scholar 

  21. 21.

    Grabda J. Life cycle and morphogenesis of Lernaea cyprinacea. Acta Parasitol Pol. 1963;11:169–98.

    Google Scholar 

  22. 22.

    Kabata Z. Parasites and diseases of fish cultured in the tropics. Abingdon: Taylor & Francis Ltd.; 1985.

    Google Scholar 

  23. 23.

    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.

    Google Scholar 

  24. 24.

    Chen R. A report on the Lernaea parasites of freshwater fishes in Pearl River Delta. Jinan Liyi Xuebao. 1981;1:64–7.

    Google Scholar 

  25. 25.

    USGS Parasitic Copepods database. https://www.glsc.usgs.gov/greatlakescopepods/Image.php?PATH=Parasite/Introduction/&SRC=42Lcrubody.jpg. Accessed 15 Aug 2018.

  26. 26.

    Shariff M. Studies on the biology, host-parasite interactions and distribution of Lernaea spp. in West Malaysia. PhD thesis. Stirling: University of Stirling; 1985.

  27. 27.

    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.

    Article  Google Scholar 

  28. 28.

    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.

    Article  Google Scholar 

  29. 29.

    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.

    Article  Google Scholar 

  30. 30.

    Kuang P. Parasitic Crustacea of freshwater fishes. Beijing: Science Press; 1991.

    Google Scholar 

  31. 31.

    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.

    Google Scholar 

  32. 32.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Shields RJ. Procedures for the laboratory rearing of Lernaea cyprinacea L. (Copepoda). Crustaceana. 1978;35:259–64.

    Article  Google Scholar 

  34. 34.

    Anderson M. A new method for non-parametric multivariate analysis of variance. Aus Ecol. 2001;26:32–46.

    Google Scholar 

  35. 35.

    Mcardle BH, Anderson MJ. Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology. 2001;82:290–7.

    Article  Google Scholar 

  36. 36.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    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.

    PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    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.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  41. 41.

    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.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    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.

    Google Scholar 

  44. 44.

    Hossain MMM, Ferdoushi J, Haque A. Biology of anchor worms (Lernaea cyprinacea). J Entomol Zool Stud. 2018;6:910–7.

    Google Scholar 

  45. 45.

    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.

    Google Scholar 

  46. 46.

    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.

    Google Scholar 

  47. 47.

    Nadler SA, León DE. Integrating molecular and morphological approaches for characterizing parasite cryptic species: implications for parasitology. Parasitology. 2011;138:1688–709.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    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.

    PubMed  Article  Google Scholar 

  51. 51.

    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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    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.

    PubMed  Article  CAS  Google Scholar 

  53. 53.

    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.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    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.

    PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    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.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  56. 56.

    Poulin R, Morand S. The diversity of parasites. Quart Rev Biol. 2000;75:277–93.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    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.

    PubMed  Article  PubMed Central  Google Scholar 

  58. 58.

    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.

    PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    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.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  60. 60.

    Dantas-Torres F. Species concepts: what about ticks? Trends Parasitol. 2018;34:1017–26.

    PubMed  Article  Google Scholar 

  61. 61.

    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.

    CAS  Article  Google Scholar 

  62. 62.

    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.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  63. 63.

    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.

    PubMed  Article  Google Scholar 

  64. 64.

    Degnan JH, Rosenberg NA. Discordance of species trees with their most likely gene trees. PLoS Genet. 2006;2:e68.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  65. 65.

    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.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Su YB, Wang LX, Kong SC, Chen L, Fang R. Complete mitochondrial genome of Lernaea cyprinacea (Copepoda: Cyclopoida). Mito DNA. 2016;27:1503–4.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr Rui Song and Min Xie for their assistance in field sampling.

Funding

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).

Author information

Affiliations

Authors

Contributions

WXL, GTW and CJH designed the study. CJH and HZ conducted the experiments. CJH and DZ conducted the data analysis. CJH and IJ wrote the paper. All authors have revised the manuscript critically for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wen X. Li.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

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.

Additional file 2: Table S2.

Sequences used for the comparative analyses in this study.

Additional file 3: Table S3.

Sequences used for the phylogenetic analyses in this study.

Additional file 4: Table S4.

Best-fit models for four datasets (cox1, 18S, 28S and concatenated), selected based on the Akaikeʼs information criterion using ModelFinder software.

Additional file 5: Table S5.

The identity values between partial cox1, 18S and 28S genes obtained in this study and Lernaea cyprinacea homologues available on GenBank.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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

Download citation

Keywords

  • Lernaea
  • Parasitic copepods
  • Phenotypic plasticity
  • Molecular taxonomy