Studies of Myxidium giardi Cépède, 1906 infections in Icelandic eels identifies a genetically diverse clade of myxosporeans that represents the Paramyxidium n. g. (Myxosporea: Myxidiidae)
Parasites & Vectorsvolume 11, Article number: 551 (2018)
The myxosporean Myxidium giardi Cépède, 1906 was described infecting the kidney of the European eel, Anguilla anguilla (L.), having spindle-shaped myxospores and terminal sub-spherical polar capsules. Since then, numerous anguillid eels globally have been documented to have similar Myxidium infections. Many of these have been identified using the morphological features of myxospores or by the location of infection in the host, and some have been subsequently synonymised with M. giardi. Therefore, it is not clear whether M. giardi is a widely distributed parasite, infecting numerous species of eels, in multiple organs, or whether some infections represent other, morphologically similar but different species of myxosporeans. The aim of the present study was to assess the status of M. giardi infections in Icelandic eels, and related fish hosts in Malaysia and to use spore morphology and molecular techniques to evaluate the diversity of myxosporeans present.
The morphologies of the myxospores from Icelandic eels were very similar but the overall dimensions were significantly different from the various tissue locations. Myxospores from the kidney of the Malaysian tarpon, Megalops cyprinoides (Broussonet), were noticeably smaller. However, the SSU rDNA sequences from the different tissues locations in eels, were all very distinct, with percentage similarities ranging from 92.93% to as low as 89.8%, with the sequence from Malaysia being even more dissimilar. Molecular phylogenies consistently placed these sequences together in a clade that we refer to as the Paramyxidium clade that is strongly associated with the Myxidium clade (sensu stricto). We erect the genus Paramyxidium n. g. (Myxidiidae) to accommodate these histozoic taxa, and transfer Myxidium giardi as Paramyxidium giardi Cépède, 1906 n. comb. as the type-species.
There is not a single species of Myxidium (M. giardi) causing systemic infections in eels in Iceland. There are three species, confirmed with a robust phylogeny, one of which represents Paramyxidium giardi n. comb. Additional species probably exist that infect different tissues in the eel and the site of infection in the host fish is an important diagnostic feature for this group (Paramyxidium n. g. clade). Myxospore morphology is generally conserved in the Paramyxidium clade, although actual spore dimensions can vary between some species. Paramyxidium spp. are currently only known to infect fishes from the Elopomorpha.
The importance of the shape of myxospores and its relative usefulness in myxosporean taxonomy has been scrutinised in recent years, and it has been unambiguously demonstrated that numerous genera are now polyphyletic in molecular phylogenetic analyses due to the use of spore morphology as the principal taxonomic measure [1,2,3]. However, myxospores remain the primary diagnostic feature of myxosporean infections in fish and other vertebrates  and their morphology, especially when combined with other important information from the host, can sometimes be sufficient to make preliminary identifications or diagnoses without the use of DNA analysis. The polyphyletic nature of some myxosporean genera in molecular analyses has occurred, in part, as our systematic framework for deciding which family and genus to place novel species has become less stringent over time, and gives priority to basic myxospore morphology over other characters . This is combined with the fact that there also seems to be a general reluctance to establish new genera when necessary, with the historical preference being to further loosen the descriptive boundaries of existing genera . This has led to artificially high numbers of species being placed in certain genera, which in some cases are clearly not closely related to each other, and has hence contributed to polyphyly in subsequent molecular phylogenies [2, 4].
Myxidium giardi Cépède, 1906 was originally described from the kidney of the European eel, Anguilla anguilla (L.), in France, with myxospores described as spindle-shaped with terminal sub-spherical polar capsules , but have since been reported in numerous sizes and described as both wide and slender, some with almost spherical polar capsules [6, 7]. Globally, numerous species of anguillid eels have been reported to have similar Myxidium spp. infections, many of which have been distinguished using morphological features such as spore size and number or lack of valvular striations [8, 9] or the site of infection in the eel (skin, gills, kidney, stomach, etc.) . However, many descriptions of Myxidium spp. infecting eels have subsequently been synonymized with M. giardi, as less emphasis was placed on details such as the site of infection  and M. giardi was considered to have an almost worldwide distribution . Currently, it is not clear whether M. giardi is a widely distributed parasite infecting numerous species of eels, in multiple organs, or whether some infections represent other, morphologically similar, species of myxosporeans. One short rDNA sequence exists for M. giardi from the urinary system of A. anguilla from Scotland, and phylogenetic analyses place this in the freshwater urinary clade, with numerous other myxosporeans found infecting the urinary systems of fish, but none with a similar spore morphology .
Anguillid eels (Order Anguilliformes) are catadromous fish consisting of 17 species in a single genus, Anguilla Garsault, distributed in all the world’s oceans. Two species inhabit the Atlantic Ocean, the European eel, A. anguilla, and the American eel, Anguilla rostrata (Lesueur) . A dramatic decline has been experienced in populations of the European eel and diseases/parasites are considered to be one of several factors responsible. Consequently, since 2009, the European eel has been on the red list of the International Union for the Conservation of Nature (IUCN) and hence is listed as a ‘critically endangered’ species . The Indo-Pacific tarpon, Megalops cyprinoides (Broussonet) (Order Elopiformes), inhabits tropical coastal and brackish waters of the Indo-Pacific, migrating between the open sea and inland rivers and mangroves . Elopiformes are related to, but do not resemble anguillid eels, and like eels they spawn at sea producing leptocephalic larvae that later migrate inland . Tarpon (Elopiformes) together with their sister group the eels (Anguilliformes) form the Elopomorpha which are one of the oldest major extant teleost lineages .
In the present study, we aimed to examine eels from Iceland and genetically characterise Myxidium giardi-like parasites from different infection sites in the hosts. In addition, we aim to compare these with morphologically similar parasites previously observed by one of us (MAF) in tarpon from Malaysia.
Thirty-one live European eels, A. anguilla, (growing yellow eels; total length range 42–68 cm) were collected from Lake Vifilsstadavatn in Iceland during May 2013 using fyke nets. The fish were transported to the laboratory, placed in tanks with freshwater and kept alive until examination. Prior to examination, the fish were euthanized in MS-222. All eels were dissected and examined for myxosporeans in the gills, kidney and stomach wall using dissecting and compound microscopes. Images were taken of twenty spores from each tissue location for size estimations using ImageJ. In addition, fourteen mature Pacific tarpon, Megalops cyprinoides total length range 22–34 cm, were caught with gill nets in Kilim mangroves on Langkawi Island, Malaysia. Tarpon were examined for gill and kidney myxosporeans, using dissecting and compound microscopes. In Malaysia, myxospores were photographed in the field using a portable compound microscope and a Dino-Eye eyepiece camera. Infected tissues from eels and tarpon were fixed in 10% buffered formalin and processed for standard histology, 3 μm thick sections prepared, stained with Giemsa and mounted in resin-based medium.
Oligochaetes (n = 100) were taken from sediment samples from the lake Vifilsstadavatn in Iceland, isolated and incubated at 10 °C in filtered (pore size 0.45 μm) lake water in 24-well tissue culture plates. Oligochaetes were examined daily using a dissecting microscope for the presence of actinospore production, typically characterized by a cloudy secretion from the anus of the worms, and confirmed using a compound microscope. Every other day, the water in the wells was exchanged. The oligochaetes that produced actinospores were used for DNA analysis and images of the actinospores taken as detailed above.
Infected tissue (kidney, stomach wall) or spore-filled cysts (gills), equal to approximately 20–40 mg, were placed directly into DNA lysis buffer for molecular analysis. Total DNA was extracted using a GeneMATRIX DNA isolation kit (EURx, Gdansk, Poland) following the tissue protocol and used as templates in subsequent PCR reactions. Parasite small subunit ribosomal DNA (SSU rDNA) was amplified using the general myxosporean primers and methodology previously described  with a new reverse primer M-790r (5'-ACG ACC AAT TAA GGC TAT GC-3'), paired with 18e, utilising the PCR conditions as in . In addition, the primer pair Mg-50f (5'-ACT AAG CCA TGC ATG TCT ATG T-3') and Mg-1170r (5'-TGA TCA ATC GAA ACG GTC TAG G-3') were designed to assist in further studies of these myxosporeans from anguillid eels and amplified a 1150 bp region of the SSU rDNA including the phylogenetically informative variable V1-V5 regions; PCR conditions were the same as in , except for using an annealing temperature of 58 °C and extension time of 1 min. PCRs were conducted on parasite DNA, from each tissue, from 4 different fish and amplified in triplicate. PCR products of the expected sizes were recovered using a GeneMATRIX PCR products extraction kit (EURx Poland) and sequencing reactions were performed using BigDyeTM Terminator cycle sequencing chemistry utilising the same oligonucleotide primers that were used for the original PCRs. DNA sequencing was performed in both forward and reverse directions for all PCR products and nucleotide BLAST searches performed for each sequence read to confirm a myxosporean origin . The contiguous sequences were obtained manually using CLUSTAL X and BioEdit [20, 21]. CLUSTAL X was used for the initial SSU rDNA sequence alignments of the novel sequences to other related myxosporean sequences identified from the BLAST searches and from relevant literature [1, 13, 22, 23]. Percentage divergence matrices were constructed from selected aligned taxa in CLUSTAL X using the neighbour-joining method based on the Kimura 2-parameter model .
Phylogenetic analyses were performed using the maximum likelihood methodology in PhyML  with the automatic smart model selection (selection criterion: Akaike Information Criterion), running the general time-reversible substitution model (GTR+G6+I) with 1000 bootstrap repeats. Bayesian inference (BI) analysis was performed using MrBayes v. 3.2 . For the BI analysis, models of nucleotide substitution were first evaluated for the alignment using MrModeltest v. 2.2 . The most parameter-rich evolutionary model based on the AIC was the general time-reversible, GTR+I+G model of evolution. Therefore, the settings used for the analysis were nst = 6, with the gamma-distributed rate variation across sites and a proportion of invariable sites (rates = invgamma). The priors on state frequency were left at the default setting (Prset statefreqpr = dirichlet (1,1,1,1)), and Chloromyxum leydigi (AY604199) was set as the outgroup. Posterior probability distributions were generated using the Markov Chain Monte Carlo (MCMC) method with four chains being run simultaneously for 1,000,000 generations. ‘Burn-in’ was set at 2500 and trees were sampled every 100 generations making a total of 7500 trees used to compile the majority rule consensus trees.
Measurements of myxospores from each organ (kidney, gills and stomach wall) were tested for normality using the Shapiro-Wilk test. To test whether there was a significant difference in analogous measurement between the different myxospores, a Student’s t-test was applied. Both tests were performed in RStudio (version 0.98.1062).
All the fish examined in this study appeared in good condition and no gross signs of a disease were present. From the 31 eels examined, small myxosporean cysts were observed on the gills of 25 fish (prevalence of 83%) using a dissection stereoscope, ranging between 2–10 cysts per gill arch. Infections in the kidney and the stomach wall were only detected using a compound microscope. Spores were detected in the kidney of 10 fish (32%) and from stomach scrapings from 6 eels (19%). Infections were light in most cases, being more apparent in the gills. All myxospores detected had a similar morphology that conformed to the currently accepted shape and size range for Myxidium giardi (Table 1) . No gill cysts were observed from 14 tarpon, examined in Malaysia, but the kidney of 4 fish (29%) were heavily infected with a myxosporean, observed in fresh tissue preparations. The myxospores from tarpon kidney had a similar morphology to the ones from the eels. However, they were noticeably, and significantly, smaller than all three species observed from the Icelandic eels (Tables 1 and 2).
The overall morphology and dimensions of the myxospores from Icelandic eels appeared similar (Fig. 1). However, taking average dimensions from 20 measured spores, they were significantly different with regard to spore length and/or width (Tables 1 and 2). Most fresh spores were typically observed in the sutural plane and had an oval to bluntly-rounded shape with almost spherical polar capsules opening sub-terminally (Fig. 1a-e). The sutural line, when visible, was sigmoidal and striations were visible on the surface of the spore valves (Fig. 1c, f), although these structures were not always clearly visible. In the valvular view, spores were more lemon-shaped than a true spindle or fusiform shape (Fig. 1e), but clearly resembled the original drawings by Cépède , however they did not have terminal capsular foramina. All myxospores in this study were between 62–70% as wide as they were long, which is similar to the original description of Myxidium giardi from eels in France (Table 1) . The myxospores from tarpon in Malaysia, had a remarkably similar overall morphology to those from Icelandic eels (Fig. 1g-k), but they were significantly smaller in size (Tables 1 and 2). When myxospores were stained, it was clear that the sporoplasm contained two nuclei, in all cases, as originally described for M. giardi . Only one oligochaete, morphologically identified as Lumbriculus variegatus, produced actinospores with an Aurantactinomyxon morphology (Fig. 1l, Table 3).
Paramyxidium giardi infecting the kidney of eels forms polysporous plasmodia which vary greatly in size (c.55–120 μm in this study). The plasmodia were most commonly found as intratubular, both in the proximal and distal kidney tubules, although they were also, to a lesser extent, observed in the renal interstitium. The associated histopathology was mostly due to mechanical disruption of the tissue as the large plasmodia often caused extensive widening of the kidney tubules followed by atrophy and degeneration of the tubular epithelial cells (Fig. 2a) which is in agreement with the pathology described by Copland . Renal infections described in his study  were more extensive than in this study. Consequently, he observed varying degrees of pathology in other parts of the kidney, e.g. the bowman capsule and the kidney interstitium.
Paramyxidium magi n. sp. is histozoic in the stomach wall. It forms polysporous plasmodia of variable sizes (c.60–95 μm) in the gastric gland interstitium (Fig. 2b). Large spore masses in the gastric glands cause focal compaction and disruption of the adjacent gastric glands similar to previous reports .
Parayxidium branchialis n. sp. is histozoic in the gills. As observed by Copland , the great majority of infections were observed in the secondary lamellae. The plasmodia are polysporous and differ significantly in size (c.35–142 μm). Infections cause separation of the lamellar epithelial cells from the basement membrane and disruption of pillar cells (Fig. 2c) and in some cases oedema in the basal part of the secondary lamellae.
Paramyxidum bulani n. sp. infects the kidney. The plasmodia are polysporous, somewhat elongated (size range: length, 25–75 μm; width, 10–22 μm) and develop as intratubular, both in proximal and distal tubules. They are seen at various stages in development, some of which are fully mature, and commonly attached to the brush border of the proximal kidney tubule epithelium. In the distal kidney tubules, the proportion of mature plasmodia seem higher (Fig 2d). No pathology was associated with infections.
No host response, such as encapsulation by fibroblasts or infiltration of immune cells, was detected in association with any of these four different myxosporeans.
SSU rDNA sequences were successfully sequenced for P. giardi (2076 bp) and for the three species recognised here as new, i.e. P. magi n. sp. (2070 bp), P. branchialis n. sp. (2082 bp), P. bulani n. sp. (2056 bp) and the Aurantiactinomyxon (865 bp), which have been assigned accession numbers in GenBank: MH414925-MH414929. BLAST searches confirmed a myxosporean origin for all sequences with closest matches in the database being to numerous actinosporean sequences and also Myxidium lieberkuehni Bütschli, 1882, the type-species of the genus Myxidium. The SSU rDNA sequences from infected tissues from eels, were all very distinct, with percentage similarities ranging from 92.93% to as low as 89.8%, over > 2000 bp of sequence data (Table 4). The SSU rDNA sequence from P. bulani n. sp. was even more dissimilar, being < 84% similar to any of the other myxospore sequences, and only 81.75% similar to the kidney-infecting P. giardi from eels (Table 4); only M. lieberkuehni was more dissimilar, being < 79% similar to any of the sequences generated from myxospore infections (Table 4). The Aurantiactinomyxon sequence did not match any of the sequences generated from the myxospores in this study and was most similar (96.07%) to P. magi n. sp.
Molecular phylogenetic analyses produced trees with very similar topologies irrespective of the methodology used (Fig. 3). All Paramyxidium spp. were fully supported in a clade with numerous actinosporean sequences (Paramyxidium clade) shown as the lilac box in Fig. 3. Paramyxidium bulani n. sp. from a tarpon in Malaysia, represented the basal sequence of the clade. The Myxidium clade (sensu stricto), that includes the type-species M. lieberkuehni, was also fully supported (red box, Fig. 3) and formed a fully supported sister group to the clade containing Chloromyxum spp. (green box, Fig. 3). Together these formed a robustly supported sister clade to the Paramyxidium clade. The Platysporina (blue box, Fig. 3) formed a fully supported sister clade to this Paramyxidium / Myxidium (s.s.) / Chloromyxum grouping. Basal to this is the hepatic biliary clade, a well-supported clade which contains numerous taxa currently assigned to the Myxidiidae, some of which infect vertebrates other than fish .
Class Myxosporea Bütschli, 1881
Suborder Variisporina Lom & Noble, 1984
Family Myxidiidae Thélohan, 1892
Genus Paramyxidium n. g.
Myxospore length width ratio between 1.4–1.6:1 resulting in a lemon-shaped form in valvular aspect and oval to bluntly-rounded form in sutural view. Valvular striations present, with sigmoidal sutural line. Polar capsules two, almost spherical. Two nuclei in the sporoplasm positioned between polar capsules. Capsular foramina open sub-terminally. Histozoic polysporous plasmodia develop in various host tissues with different species infecting different host tissues.
Type-species: Paramyxidium giardi (Cépède, 1906) n. comb.
Other species: Paramyxidium branchialis n. sp.; Paramyxidium bulani n. sp.; Paramyxidium magi n. sp.
Etymology: The new genus is named Paramyxidium n. g. as it is placed next to the type-species of the genus Myxidium (s.s.) in phylogenetic analyses.
ZooBank registration: To comply with the regulations set out in article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new genus have been submitted to ZooBank. The Life Science Identifier (LSID) of the article is: urn:lsid:zoobank.org:pub:716F81BB-DAFC-4B33-89ES-A57684065765. The LSID for the new genus name Paramyxidium is: urn:lsid:zoobank.org:act:27B9FBC6-345A-4DD9-93D6-A13EACB82A18.
Paramyxidium spp. are sufficiently different with respect to either myxospore morphology/ characteristics, tissue location, host type or DNA sequences to necessitate the erection of a novel genus within the family Myxidiidae Thélohan, 1892 to accommodate them. The family Myxidiidae currently contains five genera: Myxidium Bütschli, 1882; Coccomyxa Léger & Hesse, 1907; Cystodiscus Lutz, 1889; Enteromyxum Palenzuela, Redondo & Alvarez-Pellitero, 2002; Soricimyxum Prunescu, Prunescu, Pucek & Lom, 2007; and Zschokkella Auerbach, 1910. Species of Cystodiscus and Soricimyxum are only found in non-fish hosts . Enteromyxum spp. have very large and elongated polar capsules, spores that lack valvular striations, and are phylogenetically associated with morphologically similar forms from the gastrointestinal tract of marine fishes . Coccomyxa spp. are typically found in the gall-bladder of marine fishes, are more ellipsoidal in morphology, and are phylogenetically distantly placed in the marine gall-bladder (marine Myxidium) clade . The genera Myxidium and Zschokkella are more complicated to decipher for as both are polyphyletic, with many species likely being wrongly attributed to each genus [1, 2, 4]. However, DNA sequence data convincingly placed all Paramyxidium taxa in a robust clade which is sister to the Myxidium clade (s.s.). In addition, host differences and tissue tropism are apparent, with Paramyxidium taxa infecting the Elopomorpha in numerous organs and Myxidium spp. (s.s.) infecting the renal system of members of the Protacanthopterygii.
Therefore, we erect the genus Paramyxidium to accommodate these myxosporeans, and transfer Myxidium giardi as Paramyxidium giardi Cépède, 1906 n. comb. as the type-species. The annelid hosts for Paramyxidium spp. are likely to be oligochaetes. Currently known species are restricted to freshwater environments and fish hosts from the Elopomorpha, presumably becoming infected during migratory periods into fresh water.
Paramyxidium magi n. sp.
Type-host: Anguilla anguilla (L.) (Anguilliformes: Anguillidae).
Type-locality: Lake Vifilsstadavatn (64°4'47.39"N, 21°52'26.45"W), Iceland.
Type-specimens: Hapantotypes (histological sections and phototypes) were deposited at the Natural History Museum, London, UK, under the accession number NHMUK 2018.235.
Site in host: Stomach.
Prevalence: 19% (6/31).
Representative DNA sequences: A partial SSU rDNA sequence (2070 bp) was deposited in the GenBank database under the accession number MH414927.
ZooBank registration: To comply with the regulations set out in article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new species have been submitted to ZooBank. The LSID for the new species name Paramyxidium magi is: urn:lsid:zoobank.org:act:FEC23203-00E5-4862-8D24-FA5D6D8FFB0C.
Etymology: Specific name derived from Icelandic for stomach ‘magi’
Spore. Mature spores lemon-shaped in valvular aspect (Fig. 1a), oval to bluntly-rounded in sutural view with sigmoidal sutural line (Fig. 1b), measuring 10.8–12.9 × 7.0–8.4 (11.6 × 7.6) (L × W) (n = 20). Polar capsules almost spherical, 3.6–4.6 × 3.1–4.0 (4.0 × 3.6) (n = 20), opening sub-terminally (Fig. 1a-b). Fine spore striations visible across both spore valves (Fig. 1c).
Histozoic in the stomach wall forming polysporous plasmodia of variable sizes (c.60–95 μm) in the gastric gland interstitium (Fig. 2b). Differs from other Paramyxidium spp. primarily on the site of infection and genetic distance.
Paramyxidium bulani n. sp.
Type-host: Megalops cyprinoides Broussonet (Elopiformes: Megalopidae).
Type-locality: Kilim mangroves, Langkawi Island, Malaysia.
Other localities: Pangkor Island, Malaysia.
Type-specimens: Hapantotypes (histological sections and phototypes) were deposited at the Natural History Museum, London, UK, under the accession number NHMUK 2018.237.
Site in host: Kidney.
Prevalence: 29% (4/14).
Representative DNA sequences: A partial SSU rDNA sequence (2059 bp) was deposited in the GenBank database under the accession number MH414929.
ZooBank registration: To comply with the regulations set out in article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new species have been submitted to ZooBank. The LSID for the new species name Paramyxidium bulani is: urn:lsid:zoobank.org:act:DA7B51B0-590F-45A5-88CF-32A015B99593.
Etymology: Specific name refers to the local Malay name for the fish host (ikan bulan).
Spore. Mature spores oval to bluntly-rounded in sutural view, lemon-shaped in valvular aspect, 6.1–6.9 × 4.2–5.1 (6.7 × 4.7) (L × W) (n = 20). Polar capsules almost spherical, 2.1–2.5 × 1.5–2.0 (2.3 × 1.8) (n = 20), opening sub-terminally. Sutural line sigmoidal and inconspicuous (Fig. 1g-k), spore valves with fine striations.
Plasmodia are elongated and polysporous (size range: length, 25–75 μm; width, 10–22 μm) and develop as intratubular, in both proximal and distal kidney tubules. Differs from other Paramyxidium spp. by having much smaller myxospores and infecting the kidney of Pacific tarpon. Stained spores have a sporoplasm containing two nuclei (Figs. 1i and 2d)
Paramyxidium branchialis n. sp.
Type-host: Anguilla anguilla (L.) (Anguilliformes: Anguillidae).
Type-locality: Lake Vifilsstadavatn (64°4'47.39"N, 21°52'26.45"W), Iceland.
Type-specimens: Hapantotypes (histological sections and phototypes) were deposited at the Natural History Museum, London, UK, under the accession number NHMUK 2018.236.
Site in host: Gills.
Prevalence: 83% (25/31).
Representative DNA sequence: A partial SSU rDNA sequence (2082 bp) was deposited in the GenBank database under the accession number MH414926.
ZooBank registration: To comply with the regulations set out in article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new species have been submitted to ZooBank. The LSID for the new species name Paramyxidium branchialis is: urn:lsid:zoobank.org:act:2481F7F4-5B9F-42C8-B2EB-FED264FA7283.
Etymology: Specific name refers to site of infection in the host.
Spore. Mature spores oval to bluntly-rounded in sutural view, lemon-shaped in valvular aspect, 10.7–12.3 × 6.6–7.8 (11.6 × 7.3) (L × W) (n = 20). Polar capsules almost spherical, 3.9-4.5 × 3.2-4.2 (4.2 × 3.8) (n = 20), opening sub-terminally. Sutural line sigmoidal and inconspicuous. Fine spore striations visible across both spore valves.
Plasmodia are polysporous, differ significantly in size (c.35–142 μm), and develop in the gills (Fig. 2c), whereas other species in the genus develop in different host tissues. Differs from other Paramyxidium spp. primarily on the site of infection and genetic distance.
Morphological features, such as myxospore size and number of valvular striation, have been used to identify or differentiate between species of Myxidium infecting eels in the past [6,7,8,9,10]. However, with such conserved myxospore morphology within the Paramyxidium n. g. clade, we do not find this to be a reliable feature to successfully differentiate between species. The location of the parasite within certain specific host tissues would appear to be a far more useful characteristic, with similar tissue tropism being reported for other myxosporean clades in the past [1, 2, 4]. The European eel and its parasites have been extensively studied [5,6,7, 14] and there have been numerous reports of Myxidium giardi infecting the European eel [7, 12, 30]. However, many of these reports of M. giardi infections are from tissues other than kidney and it is therefore highly likely that many of these diagnoses are in fact detailing a morphologically similar parasite from the genus Paramyxidium but not M. giardi. Numerous taxa that fit within the Paramyxidium may already be described [8,9,10,11,12], or have been synonymised with other species due to similar spore morphologies [7, 12]. Given the clear differences in SSU rDNA data that exist between the known members of this new genus, it is imperative that such data is supplied when validating or describing novel species within this group. The Paramyxidium clade currently contains sequences from myxospores that infect various tissues of fish from the superorder Elopomorpha (tarpon/ladyfish and eels). The clade is also comprised of numerous sequences generated from actinospores that probably represent a hidden diversity of unidentified myxosporeans infecting various tissues from anguillid eels; as the sequences have all been generated from European freshwater systems that lack elopiform fish other than eels and Myxidium-like infections have been noted in numerous tissues from the European eel [31, 32].
Our phylogenies show that the Paramyxidium group is robustly supported as a sister clade to the Myxidium clade (s.s.) / Chloromyxum grouping that form as fully supported subclades (Fig. 3). These three clades all demonstrate levels of fish host association and degrees of tissue tropism. All currently known members of the Myxidium clade (s.s.) are described from the renal systems of either pike or salmonids (Protacanthopterygii), whilst the adjacent Chloromyxum sister group are all described from Cypriniformes, when infecting the gall-bladder (dashed green box, Fig. 3), with the basal species from the clade, Chloromyxum kurisis (GenBank: KJ526212), being from the urinary tract of an Atheriniformes fish. This suggests that these two clades have evolved from a common ancestor that infected the renal system of fish potentially having a Chloromyxum-like form. The Paramyxidium, as a sister to these two groups, have only been described from the Elopomorpha, with the basal species in the clade again coming from a renal infection, further indicating that the common ancestor for the whole grouping was likely to have been a renal myxosporean. A deeper extension of this phenomenon can also be applied to our tree topology. The node from which the Platysporina clade forms as a sister to the Paramyxidium/Myxidium (s.s.) clade, could also represent a common ancestor that infected the renal system of fish. Recently the Myxobilatidae (all renal-infecting) have been shown to be phylogenetically placed within the platysporinids , therefore, a renal myxosporean ancestry for the Platysporina is plausible, which could also be associated with a Chloromyxum-form. This outcome is not unexpected, as it has been previously demonstrated that the Chloromyxum morphotype is important in phylogenetic studies and appears to represent a basal morphotype for numerous myxosporean clades [22, 23, 34]. Therefore, additional sequence data from this basal morphotype, in particular when infecting more ancient fish lineages (Atheriniformes, Cypriniformes, Elopomorpha and Protacanthopterygii), will assist in future phylogenetic studies of the Myxosporea and help to better resolve the relationships between certain reproducible clades.
Certain anomalies do exist in our phylogenetic tree, which can be explained as follows. The Myxidium clade (s.s.) currently contains myxosporeans that infect the urinary system of fish from the Protacanthopterygii, a superorder of more basal teleosts, which includes the Salmoniformes (salmonids) and Esociformes (pikes). The sequence in the tree representing Sphaerospora oncorhynchi Kent, Whitaker & Margolis, 1993 (GenBank: AF201373) is, therefore, more likely to be from a fish that was co-infected with a Myxidium, likely related to the Myxidium sp. described as the CKX myxosporean from coho salmon in Canada . Indeed, studies have found both S. oncorhynchi and Myxidium salvelini Konovalov & Shul’man, 1966, in the renal systems from the same Oncorhynchus spp. [36, 37] and other striated Myxidium spp. have been reported from the urinary systems of salmonids: Myxidium minteri Yasutake & Wood, 1957; M. salvelini Konovalov et Shul’man, 1966; Myxidium noblei Konovalov, 1966; Myxidium sp. of Mavor & Strasser, 1918.
The positioning of the sequences Zschokkella sp. (GenBank: AJ581918), Chloromyxum schurovi (GenBank: AJ581917) and Myxidium giardi (GenBank: AJ582213) (dashed black box, Fig. 3), all generated from the same study , within a sub-clade of the Myxobilatidae (Platysporina) remains puzzling, especially the Myxidium/Zschokkella spore forms that exclusively belong to the Variisporina and should not therefore be present in the Platysporina . This concern is reinforced, as the sequence given for M. giardi in eels from Scotland is almost identical to that of C. schurovi from salmonids . In his description of the genus Acauda, Whipps  demonstrated the clustering of the genera Acauda, Myxobilatus and Hoferellus in phylogenetic analyses and noted that the myxospores shared similar morphological features, all having longitudinally striated spore valves and a sutural plane that bisects a pair of polar capsules. This combination of features justified the reuse of the family Myxobilatidae Shul’man, 1953, to accommodate these three genera. In addition, it has been demonstrated that Ortholinea should be included in the family Myxobilatidae, as Ortholinea spp. share similar morphological features, group within the same clade, and collectively they are all parasitic in the renal system of fish . Chloromyxum myxospores also share some synapomorphies with this group, having subspherical or elongated spores that may have caudal appendages, with a straight central suture. Chloromyxum is known to be polyphyletic genus and the retention of some features (4 polar capsules) in some derived species could explain why evolved chloromyxid forms are present in some clades, such as the Myxobilatidae clade (C. schurovi) and the hepatic biliary clade (Chloromyxum trijugum). Indeed, Holzer et. al. , reported that C. schurovi, whilst being genetically similar to members of the clade it occupied (Myxobilatidae), it shares less than 70% identity to other more basal Chloromyxum species and concluded it had evolved further than all other sequenced chloromyxids . Therefore, the inclusion of C. schurovi with the Myxobilatidae can potentially be explained in terms of DNA sequence, spore morphology and tissue tropism. However, the presence of Myxidium giardi (GenBank: AJ582213) and Zschokkella sp. (GenBank: AJ581918) still cannot, as spores with these morphologies share no synapomorphies with any members of the clade or the predicted ancestral morphotypes, which are assumed to be chloromyxids. One possibility is that the sequences generated for M. giardi and Zschokkella sp. either do not correspond to the myxospores that were observed in the tissue samples, or that they were misidentified. They may represent variants of a sequence for Zschokkella stettinensis Wierzbicka, 1987, described from the urinary bladder of the European eel . Indeed, the author of the sequences in question also indicated, in their PhD thesis, that Z. stettinensis was a possible identity for this myxosporean . The myxospores of Z. stettinensis look superficially like those of M. giardi but are also very similar to the Neomyxobolus morphotype (known to infect the renal system of freshwater fishes), being described as spores having two delicate streaks running parallel to the suture , a feature almost unique to Neomyxobolus spp. . Therefore, we suggest that these sequences represent variations of Z. stettinensis or genetically related forms, which are better placed in the genus Neomyxobolus Chen et Hsieh, 1960. Myxidium streisingeri Whipps, Murray & Kent, 2015 also groups in the Myxobilatidae and is another species showing several traits in common with the genus Neomyxobolus , indicating that it has been mistakenly placed as Myxidium. Furthermore, these species also have the tissue location (renal system) and spore characteristics to allow placement in the Myxobilatidae.
In the present study, coinfections with other myxosporeans were not observed in Icelandic eels, only Paramyxidium morphologies were observed. However, the Pacific tarpon is known to have other myxosporean parasites . Therefore, it is possible when sampling tissues such as the kidney, that blood-borne stages of non-target myxosporean could be amplified in error. This is considered to be unlikely in this case, as most of the myxosporean taxa from the Pacific tarpon from the type location have been identified [4, 43] (MAF, unpublished data). In addition, the fact that these distantly related parasites form a robust monophyletic clade, support the fact that the sequences were derived from the myxospore forms observed.
Paramyxidium giardi n. comb. (syn. Myxidium giardi) does not cause systemic infections in eels. In Iceland there are three species now confirmed, using morphological and molecular data. Additional species probably exist that infect different tissues, such as the skin, and a further member of this clade was identified here from an Aurantiactinomyxon sequence. Therefore, the site of infection in the host eels is an important diagnostic feature for this group (Paramyxidium clade). Myxospore morphology is highly conserved in the Paramyxidium clade, and although some myxospore dimensions are noted as significantly different between isolates, differentiating species based on spore dimension alone is not recommended. Despite very similar myxospore morphologies being present across the group there are relatively low genetic similarities between SSU rDNA sequences. The Paramyxidium clade is a well-supported sister to the Myxidium (s.s.) and Chloromyxum clades, each demonstrating a degree of fish host group specificity.
Basic Local Alignment Search Tool
Life Science Identifier
- s.s :
Scanning electron microscopy
Fiala I. The phylogeny of Myxosporea (Myxozoa) based on small subunit ribosomal RNA gene analysis. Int J Parasitol. 2006;36:1521–34.
Freeman MA, Yokoyama H, Ogawa K. Description and phylogeny of Ceratomyxa anko sp. n. and Zschokkella lophii sp. n. from the Japanese anglerfish, Lophius litulon (Jordan). J Fish Dis. 2008;31:921–30.
Bartošová P, Freeman MA, Yokoyama H, Caffara M, Fiala I. Phylogenetic position of Sphaerospora testicularis and Latyspora scomberomori n. gen. n. sp. (Myxozoa) within the marine urinary clade. Parasitology. 2011;138:381–93.
Freeman MA, Kristmundsson Á. Histozoic myxosporeans infecting the stomach wall of elopiform fishes represent a novel lineage, the Gastromyxidae. Parasit Vectors. 2015;8:517.
Cépède C. Myxidium giardi Cépède, et la prétendue immunité des Anguilles a l’ègard des infections myxosporidiennes. Compt Rend Séances Soc Biol Ses Fil. 1906;6:170–3.
Treasurer JW, Cox D. The occurrence of Myxidium giardi Cépède, in cultured eels, Anguilla anguilla L., in West Scotland. Bull Eur Ass. Fish Pathol. 1997;17:171–3.
Copland JW. The occurrence and distribution of Myxidium giardi Cépède, 1906 in wild and cultured European eels, Anguilla anguilla L., in England. J Fish Dis. 1981;4:231–42.
Fujita T. Studies on Myxosporidia of Japan. J Fac Agric Hokkaido Imp Univ. 1927;16:229–47.
Hine PM. Three new species of Myxidium (Protozoa: Myxosporidia) parasitic in Anguilla australis Richardson, 1848 and A. dieffenbachii Gray, 1842 in New Zealand. J R Soc N Z. 1975;5:153–61.
Ishii S. Myxosporidiosis of the Japanese eel. Zool Magazine Tokyo. 1915;27:372–82.
Hine PM. Variations in the spores of Myxidium zealandicum Hine, 1975 (Protozoa: Myxosporidia). N Z J Mar Freshwater Res. 1978;12:189–95.
Hine PM. A review of some species of Myxidium Bütschli, 1882 (Myxosporea) from eels (Anguilla spp.). J Protozool. 1980;27:260–7.
Holzer AS, Sommerville C, Wootten R. Molecular relationships and phylogeny in a community of myxosporeans and actinosporeans based on their 18S rDNA sequences. Int J Parasitol. 2004;34:1099–111.
Tesch FW. The Eel. 5th ed. Edited by Thorpe JE. Oxford, UK: Blackwell Publishing; 2003. p. 73–119.
Jacoby D, Gollock M. Anguilla anguilla. The IUCN Red List of Threatened Species 2014: e.T60344A45833138. https://doi.org/10.2305/IUCN.UK.2014-1.RLTS.T60344A45833138.en.
Shen KN, Chang CW, Iizuka Y, Tzeng WN. Facultative habitat selection in Pacific tarpon Megalops cyprinoides as revealed by otolith Sr:Ca ratios. Mar Ecol Prog Ser. 2009;387:255–63.
Tsukamoto Y, Okiyama M. Growth during the early life-history of the Pacific tarpon Megalops cyprinoides. Japan J Ichthyol. 1993;39:379–86.
Chen JN, López JA, Lavoué S, Miya M, Chen WJ. Phylogeny of the Elopomorpha (Teleostei): evidence from six nuclear and mitochondrial markers. Mol Phylogenet Evol. 2014;70:152–61.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL-X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res. 1997;24:4876–82.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Kristmundsson A, Freeman MA. Sphaeromyxids form part of a diverse group of myxosporeans infecting the hepatic biliary systems of a wide range of host organisms. Parasit Vectors. 2013;6:13.
Sanders JL, Jaramillo AG, Ashford JE, Feist SW, Lafferty KD, Kent ML. Two myxozoans from the urinary tract of topsmelt, Atherinops affinis. J Parasitol. 2015;101:577–86.
Saitou N, Nei M. The neighbour-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Guindon S, Dufayard JF, Lefort V, Anisimov M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL. Bayesian phylogenetic analysis of combined data. Syst Biol. 2004;53:47–67.
Copland JW. The pathology of Myxidium giardi Cépède, 1906 infections in wild and cultured eels, Anguilla anguilla L. J Fish Dis. 1983;6:451–60.
ICZN. International Commission on Zoological Nomenclature: Amendment of articles 8, 9, 10, 21 and 78 of the International Code of Zoological Nomenclature to expand and refine methods of publication. Bull Zool Nomencl 2012;69:161–9.
Jakob E, Walter T, Hanel R. A checklist of the protozoan and metazoan parasites of European eel (Anguilla anguilla): checklist of Anguilla anguilla parasites. J Appl Ichthyol. 2016;32:757–804.
Kristmundsson Á, Helgason S. Parasite communities of eels Anguilla anguilla in freshwater and marine habitats in Iceland in comparison with other parasite communities of eels in Europe. Folia Parasitol. 2007;54:141–53.
Kristmundsson Á. Sjúkdómar/sníkjudýr í villtum álum, Anguilla spp. á Íslandi [Diseases/parasites of wild eels, Anguilla spp., in Iceland]. MSc Thesis, University of Iceland, Reykjavík, Iceland; 2003, 1–78 (In Icelandic).
Karlsbakk E, Kristmundsson Á, Albano M, Brown P, Freeman MA. Redescription and phylogenetic position of Myxobolus aeglefini and Myxobolus platessae n. comb. (Myxosporea), parasites in the cartilage of some North Atlantic marine fishes, with notes on the phylogeny and classification of the Platysporina. Parasitol Int. 2017;66:952–9.
Fiala I, Bartošová P. History of myxozoan character evolution on the basis of rDNA and EF-2 data. BMC Evol Biol. 2010;10:228.
Jones S, Prosperi-Porta G, Dawe S, Blackbourn J, Taylor K, Lowe G, Osborn A. Proliferative renal myxosporidiosis in spawning coho salmon (Oncorhynchus kisutch) in British Columbia and Washington. Folia Parasitol. 2004;51:221–7.
Kent ML, Khattra J, Hervio DML, Devlin RH. Ribosomal DNA sequence analysis of isolates of the PKX myxosporean and their relationship to members of the genus Sphaerospora. J Aquat Anim Health. 1998;10:12–21.
Kent ML, Margolis L, Whitaker DJ, Hoskins GE, McDonald TE. Review of Myxosporea of importance in salmonid fisheries and aquaculture in British Columbia. Folia Parasitol. 1994;41:27–37.
Whipps CM. Interrenal disease in bluegills (Lepomis macrochirus) caused by a new genus and species of myxozoan. J Parasitol. 2011;97:1159–65.
Holzer AS, Sommerville C, Wootten R. Molecular identity, phylogeny and life cycle of Chloromyxum schurovi Shul’man & Ieshko 2003. Parasitol Res. 2006;99:90–6.
Wierzbicka J. Zschokkella stettinensis sp. n. (Myxospora, Bivalvulida) a parasite of eel, Anguilla anguilla (L.). Acta Protozool. 1987;26:79–82.
Holzer AS. Molecular studies on Sphaerospora truttae and other freshwater myxozoans. PhD Thesis, Stirling University, Stirling, UK, 2001. http://hdl.handle.net/1893/21826.
Chen CL, Hsieh SR. Studies on sporozoa from the freshwater fishes Ophiocephalus maculatus and O. argus of China. Acta Hydrobiol. 1960;2:171–96.
Fiala I, Hlavničková M, Kodádková A, Freeman MA, Bartošová-Sojková P, Atkinson SD. Evolutionary origin of Ceratonova shasta and phylogeny of the marine myxosporean lineage. Mol Phylogenet Evol. 2015;86:75–89.
The authors would like to thank the technical services for assistance with histology and SEM.
Funding for the study was provided by a University of Malaya Research Grant (UMRG) No: RG201-12SUS. Funding for publication costs was provided by the Centre for Conservation Medicine and Ecosystem Health, Ross University School of Veterinary Medicine.
Availability of data and materials
The data supporting the conclusions of this article are included within the article. The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The newly generated sequences were submitted to the GenBank database under the accession numbers MH414925-MH 414929. Hapantotypes, represented as histological sections and phototypes of the new species were deposited in the Zoological collections, Natural History Museum, London, UK, under the accession numbers NHMUK 2018.234-NHMUK 2018.237.
Ethics approval and consent to participate
This research was approved under The Icelandic Food and Veterinary Authority (MAST).
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.