Sphaerospora molnari is a myxozoan parasite causing skin and gill sphaerosporosis in common carp (Cyprinus carpio) in central Europe. For most myxozoans, little is known about the early development and the expansion of the infection in the fish host, prior to spore formation. A major reason for this lack of information is the absence of laboratory model organisms, whose life-cycle stages are available throughout the year. We have established a laboratory infection model for early proliferative stages of myxozoans, based on separation and intraperitoneal injection of motile and dividing S. molnari stages isolated from the blood of carp. In the present study we characterize the kinetics of the presporogonic development of S. molnari, while analyzing cellular host responses, cytokine and systemic immunoglobulin expression, over a 63-day period. Our study shows activation of innate immune responses followed by B cell-mediated immune responses. We observed rapid parasite efflux from the peritoneal cavity (< 40 hours), an initial covert infection period with a moderate proinflammatory response for about 1–2 weeks, followed by a period of parasite multiplication in the blood which peaked at 28 days post-infection (dpi) and was associated with a massive lymphocyte response. Our data further revealed a switch to a massive anti-inflammatory response (up to 1456-fold expression of il-10), a strong increase in the expression of IgM transcripts and increased number of IgM+ B lymphocytes, which produce specific antibodies for the elimination of most of the parasites from the fish at 35 dpi. However, despite the presence of these antibodies, S. molnari invades the liver 42 dpi, where an increase in parasite cell number and indistinguishable outer cell membranes are indicative of effective exploitation and disguise mechanisms. From 49 dpi onwards, the acute infection changes to a chronic one, with low parasite numbers remaining in the fish. To our knowledge, this is the first time myxozoan early development and immune modulation mechanisms have been analyzed along with innate and adaptive immune responses of its fish host, in a controlled laboratory system. Our study adds important information on host–parasite interaction and co-evolutionary adaptation of early metazoans (Cnidaria) with basic vertebrate (fish) immune systems and the evolution of host adaptation and parasite immune evasion strategies.

It is assumed that complex life cycles in cnidarian parasites belonging to the Myxozoa result from incorporation of vertebrates into simple life cycles exploiting aquatic invertebrates. However, nothing is known about the driving forces and implementation of this event, though it fostered massive diversification. We performed a comprehensive search for myxozoans in evolutionary ancient fishes (Chondrichthyes), and more than doubled existing 18S rDNA sequence data, discovering seven independent phylogenetic lineages. We performed cophylogenetic and character mapping methods in the largest monophyletic dataset and demonstrate that host and parasite phylogenies are strongly correlated, and that tectonic changes may explain phylogeographic clustering in recent skates and softnose skates, in the Atlantic. The most basal lineages of myxozoans inhabit the bile of chondrichthyans, an immunologically privileged site and protective niche, easily accessible from the gut via the bile duct. We hypothesize that feed-integration is a likely mechanism of host acquisition, an idea supported by feeding habits of chimaeras and ancient sharks and by multiple entries of different parasite lineages from invertebrates into the new host group. We provide exciting first insights into the early evolutionary history of ancient metazoan parasites in a host group that embodies more evolutionary distinctiveness than most other vertebrates.

Intertidal rock pools where fish and invertebrates are in constant close contact due to limited space and water level fluctuations represent ideal conditions to promote life cycles in parasites using these two alternate hosts and to study speciation processes that could contribute to understanding the roles of parasitic species in such ecosystems. Gall bladder and liver samples from five clinid fish species (Blenniiformes: Clinidae) were morphologically and molecularly examined to determine the diversity, prevalence, distribution and host specificity of Ceratomyxa parasites (Cnidaria: Myxozoa) in intertidal habitats along the coast of South Africa. Phylogenetic relationships of clinid ceratomyxids based on the SSU rDNA, LSU rDNA and ITS regions were assessed additionally to the investigation of population genetic structure of Ceratomyxa cottoidii and subsequent comparison with the data known from type fish host Clinus cottoides. Seven Ceratomyxa species including previously described Ceratomyxa dehoopi and C. cottoidii were recognized in clinids. They represent a diverse group of rapidly evolving, closely related species with a remarkably high prevalence in their hosts, little host specificity and frequent concurrent infections, most probably as a result of parasite radiation after multiple speciation events triggered by limited host dispersal within restricted spaces. C. cottoidii represents the most common clinid parasite with a population structure characterized by young expanding populations in the south west and south east coast and by older populations in equilibrium on the west coast of its distribution. Parasite and fish host population structures show overlapping patterns and are very likely affected by similar oceanographic barriers possibly due to reduced host dispersal enhancing parasite community differentiation. While fish host specificity had little impact on parasite population structure, the habitat preference of the alternate invertebrate host as well as tidal water exchange may be additional crucial variables affecting the dispersal and associated population structure of C. cottoidii.

The environmental DNA (eDNA) metabarcoding approach has become a useful tool for detecting the species diversity of different animal groups, including parasites. Myxozoa (Malacosporea and Myxosporea) represent a unique group of morphologically simplified endoparasites that mainly infest fish and whose diversity remains largely unexplored. Metabarcoding of DNA from the aquatic environment is a promising non‐invasive method that allows us to assess myxozoan biodiversity at a given site. This is essential not only for describing myxozoan communities, but also for the development of disease control methods. Using an alignment of 330 sequences, we employed in silico PCR to score primer pairs, designed to amplify the V4 region of the SSU rDNA of different myxosporean sublineages comprising the entire diversity of oligochaete‐infecting (freshwater) myxosporeans. We selected eight clade‐specific primer sets for metabarcoding, avoiding amplification of DNA from other organisms present in eutrophic freshwaters. The metabarcoding approach used in the analysis of eDNA sediment samples detected a high myxosporean diversity even in small sample volumes (in total 44 OTUs). Furthermore, metabarcoding analysis of myxosporeans in fish tissue samples selected for primer testing revealed 91 different myxosporean OTUs, more than double the number obtained by classical PCR screening and Sanger sequencing with general myxozoan primers and almost seven times higher detection than by microscopic examination. Our results further suggest quantitative sampling requirements for realistic future diversity estimates by comparing OTUs from fish tissue metabarcoding and eDNA samples. The use of specific primer sets enabled the detection of a high proportion of myxosporean reads (63–100%) in all datasets, even in highly eutrophic habitats. This shows our metabarcoding approach as an excellent tool for non‐invasive and sensitive detection of myxosporean biodiversity in aquatic sediments, potentially useful for monitoring myxozoan disease agents that threaten economically important fish in aquaculture.

Myxobolus pseudodispar Gorbunova, 1936 (Myxozoa) was originally described as a parasite of common roach, Rutilus rutilus (Linnaeus), with developing stages in muscles and spores disseminated in macrophage centres of different organs and tissues. Later, this parasite was described from several other cyprinids, but with relatively large intraspecific differences based on SSU rDNA gene sequences. Within our long-term study on myxozoan biodiversity, we performed a broad microscopic and molecular screening of various freshwater fish species (over 450 specimens, 36 species) from different localities. We investigated the cryptic species status of M. pseudodispar. Our analysis revealed four new unique SSU rDNA sequences of M. pseudodispar as well as an infection in new fish host species. Myxobolus pseudodispar sequence analysis showed clear phylogenetic grouping according to fish host criterion forming 13 well-recognised clades. Using 1% SSU rDNA-based genetic distance criterion, at least ten new species of Myxobolus B?tschli, 1882 may be recognised in the group of M. pseudodispar sequences. Our analysis showed the paraphyletic character of M. pseudodispar sequences and the statistical tests rejected hypothetical tree topology with the monophyletic status of the M. pseudodispar group. Myxobolus pseudodispar represents a species complex and it is a typical example of myxozoan hidden diversity phenomenon confirming myxozoans as an evolutionary very successful group of parasites with a great ability to adapt to a new hosts with subsequent speciation events.

Myxidium rhodei Léger, 1905 (Cnidaria: Myxozoa) is a kidney-infecting myxosporean that was originally described from the European bitterling Rhodeus amarus. Subsequently, it has been documented based on spore morphology in more than 40 other cypriniform species, with the roach Rutilus rutilus being the most commonly reported host. This study introduces the first comprehensive data assessment of M. rhodei, conducted through morphological, ecological and molecular methods. The morphological and phylogenetic analyses of SSU rDNA sequences of Myxidium isolates obtained from European bitterling and roach did not support parasite conspecificity from these fish. In fact, the roach-infecting isolates represent three distinct parasite species. The first two, M. rutili n. sp. and M. rutilusi n. sp., are closely related cryptic species clustering with other myxosporeans in the freshwater urinary clade, sharing the same tissue tropism. The third one, M. batuevae n. sp., previously assigned to M. cf. rhodei, clustered in the hepatic biliary clade sister to bitterling-infecting M. rhodei. Our examination of diverse cypriniform fishes, coupled with molecular and morphological analyses, allowed us to untangle the cryptic species nature of M. rhodei and discover the existence of novel species. This underscores the largely undiscovered range of myxozoan diversity and highlights the need to incorporate sequence data in diagnosing novel species.