The susceptibility of diverse species of cultured oligochaetes to the fish parasite Myxobolus pseudodispar Gorbunova (Myxozoa)

Serine protease inhibitors of the whirling disease parasite Myxobolus cerebralis (Cnidaria, Myxozoa): Expression profiling and functional predictions

Here, we studied the expression pattern and putative function of four, previously identified serine protease inhibitors (serpins) of Myxobolus cerebralis, a pathogenic myxozoan species (Cnidaria: Myxozoa) causing whirling disease of salmonid fishes. The relative expression profiles of serpins were determined at different developmental stages both in fish and in annelid hosts using serpin-specific qPCR assays. The expression of serpin Mc-S1 was similar throughout the life cycle, whereas a significant decrease was detected in the relative expression of Mc-S3 and Mc-S5 during the development in fish, and then in the sporogonic stage in the worm host. A decreasing tendency could also be observed in the expression of Mc-S4 in fish, which was, however, upregulated in the worm host. For the first time, we predicted the function of M. cerebralis serpins by the use of several bioinformatics-based applications. Mc-S1 is putatively a chymotrypsin-like inhibitor that locates extracellularly and is capable of heparin binding. The other three serpins are caspase-like inhibitors, and they are probably involved in protease and cell degradation processes during the early stage of fish invasion.

Myxozoan hidden diversity: the case of Myxobolus pseudodispar Gorbunova, 1936

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.

Transcriptomic analysis of oligochaete immune responses to myxosporeans infection: Branchiura sowerbyi infected with Myxobolus cultus

The Myxozoa are endoparasites characterized by a two-host life cycle that typically involves invertebrates and vertebrates as definitive and intermediate hosts, respectively. However, little is known about invertebrate-myxosporean interactions, particularly about patterns of host immune defense. We used RNA-sequencing to identify genes that are possibly involved in the immune responses of the oligochaete Branchiura sowerbyi naturally infected with Myxobolus cultus. De novo assembly of the B. sowerbyi transcriptome yielded 119,031 unigenes, with an average length of 896 bp and an N50 length of 1754 bp. Comparative transcriptome analysis revealed 4059 differentially expressed genes (DEGs) between M. cultus-infected and uninfected B. sowerbyi groups, including 3802 upregulated genes and 257 downregulated genes. Among the B. sowerbyi immune factors implicated in the responses to M. cultus infection, DEGs related to lectins, ubiquitin-mediated proteolysis, phagocytosis, oxidative-antioxidative responses, proteases, and protease inhibitors were upregulated. The expression of some immune-related molecules such as calmodulin, heat shock proteins, antimicrobial peptides, lysenin, and serum amyoid A protein were also significantly upregulated. The expression patterns of 14 immune-related DEGs identified by RNA-seq were validated by quantitative real-time polymerase chain reaction. This study is the first attempt to characterize the B. sowerbyi transcriptome and identify immune-related molecules possibly associated with M. cultus infection. It is also the first report of invertebrate host-myxosporean interactions at the transcriptomic level. Our results will facilitate the elucidation of adaptive evolution mechanisms of myxosporean parasites in the definitive host and the genetic basis for differences in resistance of invertebrate hosts of different genotypes to a myxosporean species.

Selection of suitable reference genes for gene expression studies in myxosporean (Myxozoa, Cnidaria) parasites

Myxozoans (Cnidaria: Myxozoa) are an extremely diversified group of endoparasites some of which are causative agents of serious diseases in fish. New methods involving gene expression studies have emerged over the last years to better understand and control myxozoan diseases. Quantitative RT-PCR is the most extensively used approach for gene expression studies. However, the accuracy of the results depends on the normalization of the data to reference genes. We studied the expression of eight commonly used reference genes, adenosylhomocysteinase (AHC1), beta actin (ACTB), eukaryotic translation elongation factor 2 (EF2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1), DNA-directed RNA polymerase II (RPB2), 18S ribosomal RNA (18S), 28S ribosomal RNA (28S) across different developmental stages of three myxozoan species, Sphaerospora molnari, Myxobolus cerebralis and Ceratonova shasta, representing the three major myxozoan linages from the largest class Myxosporea. The stable reference genes were identified using four algorithms: geNorm, NormFinder, Bestkeeper and ΔCq method. Additionally, we analyzed transcriptomic data from S. molnari proliferative and spore-forming stages to compare the relative amount of expressed transcripts with the most stable reference genes suggested by RT-qPCR. Our results revealed that GAPDH and EF2 are the most uniformly expressed genes across the different developmental stages of the studied myxozoan species.

Susceptibility-related differences in the quantity of developmental stages of Myxobolus spp. (Myxozoa) in fish blood

Here, we investigated the early development of two closely related myxozoan parasites, the highly pathogenic Myxobolus cerebralis, the causative agent of the whirling disease in salmonids, and Myxobolus pseudodispar, a common, non-pathogenic parasite of cyprinids. The aim of our study was to examine under in vivo laboratory conditions whether fish blood is involved in the intrapiscine development of the two parasite species and investigate if there is dissimilarity between the parasite infection intensity in blood and if it varies in terms of host susceptibility and parasite pathogenicity. Highly susceptible, less susceptible and non-susceptible hosts were involved. Blood samples were taken 1 day, 1 week and 1 month post exposure to M. cerebralis and M. pseudodispar, respectively. The prevalence and infection intensity was estimated by parasite-specific quantitative real-time PCR. Although previous findings assumed that M. cerebralis might escape from host immune system by migrating via peripheral nerves, our experimental results demonstrated that M. cerebralis is present in blood during the early stage of intrapiscine development. For the non-pathogenic M. pseudodispar, the highest infection prevalence was found in the original host, common roach Rutilus rutilus, whereas the highest infection intensity was detected in rudd Scardinius erythrophthalmus, a “dead-end” host of the parasite. The presence of M. pseudodispar developmental stages in the blood of both susceptible and non-susceptible cyprinids suggests that the susceptibility differences remain hidden during the early stage of infection. Our findings supply further evidence that host specificity is not determined during the early, intrapiscine development involving the vascular system. Furthermore, we found remarkable differences in the infection dynamics of the two parasite species examined, possibly due to their distinct pathogenicity or variations in adaptive capabilities to immune components in host blood.