The kinetics of cellular and humoral immune responses of common carp to presporogonic development of the myxozoan Sphaerospora molnari
Parasites & Vectorsvolume 12, Article number: 208 (2019)
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.
Myxozoans are possibly the oldest metazoan parasites to diverge early in the evolution of free-living cnidarians (650 Mya), initially as parasites of annelid and bryozoans . The acquisition of fish as secondary hosts later in evolution facilitated alternative transmission and dispersion strategies resulting in a distinct success. The conquest of a more evolved intermediate host, however, imposed important challenges on myxozoans, especially because of the hosts’ ability to induce parasite-specific adaptive immune responses [2, 3]. While the innate immunity of annelids has a limited capacity to “remember” pathogens, fish are vertebrates and possess a sophisticated system of antigen-specific B and T cells capable of developing memory with the aim to respond more vigorously and effectively to repeated exposures . Since the initial encounter of myxozoans with fish approximately 440 Mya, these two groups co-diversified and co-evolved by continuous adaptation of immune responses and evasion strategies. This possibly contributed to the reduced antagonism of the parasite , explaining why out of the approximately 2500 described myxozoan species, only a small percentage are represented by serious pathogens . This small percentage, however, causes severe diseases and important economic losses in aquaculture and wild fish stocks. Notably, the highest severity of the diseases is observed in newly acquired species, which do not have a common history with the parasite. With this regard, it is noteworthy that the proliferation and disease severity of myxozoans has been linked to increasing water temperature and epidemical models predict spread of disease and major outbreaks as a result of climate change [7, 8]. Resulting (re-)emergence and dissemination of myxozoan diseases into new areas therefore represents a major threat. Although concentrated efforts have resulted in significant advances in our understanding of the biology, transmission and pathogenicity of these parasites, detailed knowledge about host–parasite interaction throughout disease progression, and the ability of the host’s immune system to eliminate the parasite and protect from re-infection is still limited .
A major reason for limited information on host–parasite interactions is the absence of available in vivo models. Only 55 life-cycles have been elucidated to date  and very few (3–4) are continuously perpetuated in research laboratories, as their maintenance is laborious and time-consuming . We have recently established the first continuous in vivo model system for myxozoan pre-sporogonic proliferative stages. The model employs Sphaerospora molnari, a common parasite of carp in pond cultures and natural freshwater habitats in central and eastern Europe [11, 12]. In natural infections, S. molnari invades the epithelia of skin and gills where it forms spores intracellularly . Importantly, prior to spore formation, S. molnari and other species from the same phylogenetic clade multiply in fish blood [14, 15]. The pluricellular blood stages, initially termed unidentified blood objects (UBOs)  proliferate rapidly and show a unique motility which may help them to evade contact with host immune cells . Our research model takes advantage of these blood stages and bypasses problems related to the development and spore maturation in the invertebrate host, by transferring the parasite via intraperitoneal injection from fish to fish, in a similar fashion as performed by Ibarra et al. [18, 19] for the myxozoan Ceratonova shasta, using parasite stages obtained from ascites of infected fish. This allows for functional analyses of early parasite establishment, migration, development, multiplication and a better understanding of the cross-talk between host and pathogen, under controlled laboratory conditions.
Previously, a number of studies have aimed to elucidate various aspects of the teleost immune response towards myxozoans. Although direct comparison is complicated, they have provided valuable information on the activation of both innate and adaptive immune responses [9, 20, 21]. As with other pathogens, myxozoans avoid the immune system by intracellular proliferation or by exploiting immune privileged sites, e.g. the nervous tissues, where the presence of leukocytes and immune mediators is limited [22,23,24,25]. Despite these efforts, infections lead to the recruitment of macrophages, neutrophils and other granulocytes to infected tissues [26,27,28] and increased production of innate humoral factors including peroxidases, lysozyme or complement, which are involved in killing these parasites . The changes in gene expression induced by myxozoan infections are often only weak or transient with downregulation of pro-inflammatory cytokines and a dominant anti-inflammatory phenotype [29,30,31,32]. With regard to adaptive immunity, a number of reports demonstrated an increase in the number of B cells in infected tissues and presence of myxozoan-specific antibodies in the sera or mucus of survivor fish (summarized in [9, 21, 33]). Although the development of these antibody responses is slower than to other antigens, some antibodies showed partial to complete protection to re-infection [34, 35]. Most of the available data, however, focuses on myxozoans within tissues such as intestine or kidney, and the information on developmental stages in the blood is scarce. Undoubtedly, living in the blood stream, exposed to the full array of humoral and cellular immune effector mechanism, must represent a considerable challenge for the parasite survival, and while motility may explain its escape from effector leukocytes , it is certainly not the only evasion strategy S. molnari exploits. To date, information on parasite replication and changes in the innate and adaptive immunity of the host is not known. The key aims of the presented study thus were to broaden our understanding of fish–myxozoan interaction by analyzing simultaneously and over time, the kinetics of S. molnari presporogonic infection in the host’s body tissues together with an overall activation of the immune system on a cellular and molecular level, over a 63 day period. To shed some light on these questions, we analyzed the expression of selected pro- and anti-inflammatory cytokines in the head kidney, namely il-1β, il-6, tnfα, ifnγ, il-11 and il-10, and investigated changes in the composition of blood leukocytes, as a first indicator of activation of the immune system. Furthermore, to elucidate the development of adaptive immunity, B cell kinetics and the secretion of specific antibodies in the blood were analyzed throughout the experiment.
Parasite kinetics over time: Acute vs chronic infection and liver as an unknown reservoir
The combination of a controlled infection study using our S. molnari laboratory model and a specific qPCR assay which is sensitive enough to consistently detect two parasite stages in a 4 µl blood DNA extraction (data not shown), allowed us to study the number and distribution of S. molnari in different organs during 63 days post-intraperitoneal injection (Fig. 1b; Additional file 1: Figure S1), thereby elucidating the kinetics of infection and identifying previously unknown parasite locations. Sphaerospora molnari was first detected after a latency period of two to three weeks. At 14 days post-infection (dpi), S. molnari was quantifiable by qPCR in one out of five experimentally infected fish. At 21 dpi, all fish tested S. molnari-positive by qPCR and by microscopy (parasites detected in the leukocyte fraction of blood spun in hematocrit capillaries). The highest peak of parasite numbers in the blood was determined at 28 dpi (Fig. 1a). At this time point, parasite SSU rDNA copy numbers in different organs were similar or lower than those in the blood and it can be assumed that S. molnari is represented mainly by proliferative blood stages in the vascular system of these organs. After a decline in all organs but the gills, 35 dpi, parasite numbers showed a second peak of abundance in the blood, 42 dpi, coincident with higher parasite numbers in kidney, muscle, gills and the liver. Notably, in the liver, the maximum S. molnari concentrations were determined, approx. 7 times higher than those detected in the blood on the same day, and 5 times higher than during peak blood proliferation, 28 dpi. The stages in the liver represent a previously unknown parasite location in the host. At 56 and 63 dpi, parasite numbers showed a significant decline (Kruskal-Wallis H-test, χ2 = 6.731, df = 8, P < 0.001) and the acute stage infection changes to a chronic type. Despite noticeable S. molnari gill invasion, 35–42 dpi, we did not observe spore formation in the gills at this time or at any other time point during the experiment.
In Kwik–Diff-stained blood smears and imprints, the cellular composition of S. molnari life-cycle stages was studied. Blood stages (BS) were composed predominantly of small bicellular units (one primary cell harboring one secondary cell); only 17% were larger with 3–20 cells (primary cell containing various secondary cells and sometimes secondary-tertiary cell doublets; Fig. 1b). Until 28 dpi, BS were composed of 1–9 cells; the first fish to show elevated numbers of SSU rDNA copies in the liver (relative copy number of 11.9, Fig. 1) had larger BS stages with up to 20 cells, which also occurred on most of the subsequent days (49, 56 and 63 dpi). However, no statistically significant changes in the cellular composition of S. molnari BS over time were determined (Kruskal-Wallis H-test, χ2 = 3.527, df = 4, P = 0.563), similar to a previous study  (referred to as ‘Csaba cells’ but likely representing S. molnari). BS were often found in contact with host erythrocytes (Fig. 1b, arrowheads) but never within them. In Kwik–Diff-stained imprints of the liver, liver stages (LS) were found to be considerably larger in size than BS and composed of an average of 7 cells (Fig. 1c), but with a range in cellularity identical to that of blood stages (2–20 cells per stage). In LS, secondary cells often occurred in pairs in the cytoplasm of primary cells. Tertiary cells were present only in stages with 6 or more secondary cells and only single tertiary cells were present inside individual secondary cells, the same as in BS. Some parasites showed a separation in a refractile ecto- and more granular endoplasm, usually related to migratory cell motility via a protruding edge (Fig. 1c; [37, 38]). The number of parasite stages detected on Kwik–Diff-stained slides of blood smears correlated significantly with the values obtained by qPCR (linear regression, r2 = 0.2282, P = 0.0182) and allowed the calculation of absolute numbers: BS on 28 dpi reached an average density of 9,390 and a maximum density of 11,350 BS/µl of blood. Based on Ct values and average cellular composition up to 29,480 LS/µg liver tissue were estimated.
BS of S. molnari stained with carboxyfluorescein succinimidyl ester (CSFE) were detectable in peritoneal lavages 16 h post-infection (hpi) (3.16 ± 0.58% of cells) to 40 hpi (0.07 ± 0.11 of cells), but not at 7 or 21 dpi. BS were also detected in the spleen at very low levels (0–12 cells/100,000 cells) but not in head kidney or blood, 16 hpi. At 21 dpi, one out of three fish was microscopically positive for S. molnari BS in the blood but all parasites were CSFE-negative. This indicates that, at this time point, that the parasites had undergone more than 7 cell divisions resulting in 1/128th of their original fluorescence and limiting their cytometric detection .
Host reaction to infection: macroscopic changes and massive lymphocyte response
Previously, systemic pathological changes causing massive swelling of hematopoietic organs of common carp were reported for S. dykovae infections [11, 40, 41]; however, the presence of other myxozoans in these cases cannot be excluded . In the isolated S. molnari infection model, we observed macroscopic changes from 42 dpi onwards and until the end of the experiment, with head kidneys, trunk kidneys and livers being visibly enlarged in all experimental fish. Additionally, and exclusively on 21 dpi, the first day where BS were present in all experimental fish, we noted a blood color change from dark red to brown, without detecting hemolysis in the plasma of centrifuged blood. Furthermore at 21 and 28 dpi, a large percentage of erythrocytes showed cytoplasmic vacuolization and an unusually high number of immature erythrocytes were present. However, erythrocyte number, hematocrit and hemoglobin concentration differed significantly from control fish only at three time points: at 4 and 7 dpi an unexpected decrease of these values was observed in the control group (LM erythrocytes: F(23,68) = 3.05, P = 0.06, 4 and 7 dpi; LM hematocrit: F(23,65) = 4.17, P = 0.01, 4 dpi; LM hemoglobin: F(23,68) = 3.82, P = 0.03, 4 and 7 dpi) and at 56 dpi hematocrit and erythrocyte numbers were meaningfully higher in infected fish than in the control group (Fig. 2, erythrocytes; Additional file 1: Figure S2, hemoglobin and hematocrit).
In the blood, thrombocytes and all leukocyte types showed a substantial increase over time in response to S. molnari infection in carp (Fig. 2) with a maximum increase in lymphocytes (up to 13-fold), especially B cells and monocytes (up to 13-fold), followed by neutrophils (up to 10-fold) and thrombocytes (up to 3-fold), in individual fish. Thrombocyte numbers were higher in infected fish on all days except 1, 2, 14 and 42 dpi (Fig. 2). From 35 dpi onwards, we observed an increase in the number of thrombocytes, with the exception of 42 dpi, when S. molnari resided in the liver and fewer thrombocytes were detected in the blood. This decrease at 42 dpi and during the occurrence of LS was also noticed in other immune cells such as lymphocytes and monocytes, while less obvious in neutrophils. Neutrophil and monocyte numbers (Fig. 2) increased from 14 dpi onwards and remained elevated until the end of the experiment. Neutrophils first peaked at 14 dpi and then reached their maximum at 35 dpi, when S. molnari blood stages are reduced in the blood. Monocytes also first increased 14 dpi and were considerably elevated at 21, 28, 35 and 56 dpi, while showing a decrease towards the end of the experiment (63 dpi). As the most obvious cellular response to S. molnari infection, the number of lymphocytes in infected fish increased massively, from 21 dpi, when S. molnari appeared for the first time in the blood. Numbers thereafter increased considerably, with a significant difference to the control group at all time points except for 42 dpi, and a slight decrease noticeable at 63 dpi, similar to monocytes and coincident with the initiation of a chronic low-level infection by S. molnari at 49 dpi.
IgM expression and B cell responses indicate specific acquired immunity but point to a polyclonal B cell activation
We focused on the changes in the expression of systemic IgM, the kinetic of IgM+ B cells and the presence of specific IgM in the serum of infected fish. IgM gene expression was analyzed in the head kidney, with differential TaqMan assays for IgMmem and IgMsec. We demonstrated that IgMmem expression (Fig. 3; green box plots) was higher than in the control group, from 21 dpi onwards and increased steadily until 42 dpi (significantly higher when compared from 1 to 7 dpi, LM: F(11,43) = 4.9, P < 0.01), where it peaked at > 14-fold expression in individual fish and decreased considerably to close to normal levels thereafter. The prominent role of IgM+ B cells was further supported by a growing number of the IgM+ B cells in peripheral blood. Thus, while the number of B cells remained stable in the blood of fish injected with phosphate-buffered saline (PBS) and ranged between 1.3 and 3 × 103 IgM+ lymphocytes/µl of blood, the infection induced an increase from 21 dpi onwards, peaking at 13 × 104 IgM+ B cells/µl of blood (approx. 60-fold increase), at 56 dpi (Fig. 2).
We further investigated whether the increased number of B cells was linked to an elevated secretion of IgM into the serum. The results of the qPCR analysis indicated that IgMsec expression in the head kidney (Fig. 3; purple box plots) increased significantly from 14 dpi onwards, with an initial peak at 14 dpi (approx. 4-fold expression) and highest levels detected at 42 dpi (up to 10-fold expression in individual fish) (LM: F(11,44) = 8.2, P < 0.001 to 0.049). Thereafter IgMsec expression decreased slightly but was maintained at elevated levels, until the end of the experiment. Although these data suggest an involvement of humoral responses in the infection with S. molnari, it remains to be determined whether the produced antibodies are specific or if the infection merely induces a polyclonal IgM hypersecretion with no effective binding of the parasite. Polyclonal activation is a strategy used by parasites to escape the host-specific immune response by means of diluting pathogen-specific antibodies while increasing irrelevant antibodies (e.g. [43,44,45,46,47]). To elucidate this question, BS of S. molnari were used as antigen in western blotting with sera of experimentally infected fish. The parasite lysate yielded 6 apparent protein bands in sizes of 13, 16, 28, 48, 55 and 70 kDa (Fig. 3). However, after membrane transfer, we were able to demonstrate specific IgM binding to a parasite protein of 16 kDa, from 28 dpi onwards, with increasing intensity over time, peaking at 63 dpi (Fig. 3, blot and red line in figure). No signal was detected in the serum of fish injected with PBS. Taken together, these data provide strong evidence for the induction of a specific carp IgM antibody response to S. molnari.
Cytokine expression in response to S. molnari shows strong involvement of il-10
The expression of six pro- and anti-inflammatory cytokines was studied in the head kidney of experimental and control fish (Fig. 4). Within the first week after the injection of S. molnari into the peritoneal cavity (latent period of infection), we observed only a mild increase in the expression of the pro-inflammatory cytokines driving the acute phase of the immune response. After this initial stage, expression followed these patterns: (i) il-1β, ifnγ and il-11 expression was similar or lower than in the control; (ii) il-6 and tnfα expression was elevated throughout the experiment; and (iii) expression of il-10 increased gradually to the highest levels at day 56. More specifically, the expression of ifnγ, which orchestrates innate and adaptive responses against intracellular pathogens and supports the development of pro-inflammatory macrophages , was highest at 1 dpi with 2.8-fold upregulation, but decreased to levels below those observed in control fish as early as 7 dpi. Similarly, the expression of il-1β, a key regulator of inflammation, increased only slightly during first three days with a maximum 1.5-fold increase at 2 dpi and a subsequent decrease below the expression level observed in control fish, likely indicating inhibition, from 7 dpi onward and for the remainder of the experimental infection (lower expression over time when compared to initial dpi, LM: F(11,44) = 9.1, P < 0.001 to 0.014). The expression of il-11 has previously been shown to be involved in the induction of acute phase proteins and in the maturation of thrombocytes . This member of the il-6 family reached its highest (12-fold) expression at 1 dpi but decreased to control levels by 2 dpi and remained only slightly elevated, with an increase on day 28, 49 and day 63. Conversely, another mediator of acute phase response, il-6 exhibited increased expression from 1 dpi and throughout the whole study, with peaks occurring at 4 dpi (25-fold expression), 21 dpi (19-fold expression) and 56 dpi (21-fold expression). The infection also induced a 4- to 5-fold increase in the expression of tnfα with occasional decreases at 4 and 49 dpi, indicating a prolonged activation of innate arm of the immune response . In contrast to the above mentioned cytokines, the expression of the anti-inflammatory il-10 remained unchanged during the first week, but increased gradually with a 3-fold upregulation at 7 dpi, a 17-fold upregulation at 14 dpi and a subsequent increase over time, reaching up to 1456-fold expression in individual fish (average 921-fold, 56 dpi; Fig. 4; LM: F(11,43) = 39.4, P < 0.001 to 0.047).
For myxozoans, knowledge of the presporogonic intrapiscine development is scarce and detailed information on the course and location of all stages is limited to a few host–parasite systems. The unavailability of such data is explained by problems related to running experimental studies, i.e. lack of known myxozoan life-cycles, limits regarding the production of infective stages in invertebrate hosts, or lacking knowledge of the seasonality of highly infective spore concentrations in natural habitats. In this study, we bypassed these problems as well as potential co-infections by producing S. molnari BS in an in vivo infection system in laboratory-reared specific pathogen-free (SPF) carp.
Using qPCR, S. molnari was first detected in the blood at 14 and 21 dpi, indicating a prolonged period of covert existence in the host. In natural infections of other myxozoans, invasion is fast and initial cycles of parasite multiplication close to the invasion site are initiated immediately [22, 25, 51, 52]. It may be argued that, due to the artificial intraperitoneal introduction route, S. molnari requires additional time to find its entrance to the vascular system. However, we demonstrated that the majority of CSFE-labeled parasites had departed the peritoneal cavity within 40 hpi, possibly via the peritoneal efflux , or surface receptors, which allow it to actively navigate through host tissues to a desired location . We suggest that S. molnari enters specific host cells and starts intracellular proliferation, prior to extracellular proliferation in the blood. Our hypothesis is partially supported by a 2-week delay in the detection of first leukocyte responses, yet the increased expression of pro-inflammatory cytokines in the head kidney during the first four days of infection suggest an activation of systemic immune response on first contact, albeit a mild one. Despite our efforts, the exact location of the parasite before its occurrence in the blood at 14 to 21 dpi remains enigmatic, as parasite DNA was not found in any of the tested tissues.
Once in the blood, cell doublets, the main presporogonic proliferation stage in myxozoans  are predominant. Parasite appearance in the blood coincides with the increase in the number of neutrophils and monocytes, which are essential for the host responses via production of reactive oxygen species, netosis, phagocytosis, activation of the immune system and antigen presentation [56,57,58]. Such changes in the mobilization of leukocytes were described previously in infection experiments with other myxozoan parasites, although a direct comparison of these findings is rather complicated, as the parasites differ in their organ tropism, the timing of the infection and chosen sampling points. In gilthead sea bream, a mobilization of mast cells and Ig+ cells to the intestinal infection site, together with a decrease in the number of granulocytes in the head kidney was observed 40 days after exposure to Enteromyxum leei, the agent of myxozoan emaciation disease . Similarly, infection with Kudoa thrysites in the muscle of Atlantic salmon leads to the activation of macrophages and dendritic cells, i.e. phagocytes and antigen presenting cells . In our study, the elevated expression of cytokines and recruitment of monocytes and neutrophils to blood indicates successful recognition of the parasite and activation of inflammatory responses, yet, the killing capacity of these cells proved insufficient. Partially, this could be attributed to motility, which S. molnari uses to avoid cellular contact with monocytes and neutrophils  but also to additional strategies used by other myxozoan parasites .
We identified 35 dpi as a key point in the course of the S. molnari infection. At this time, the number of parasites was reduced in all tissues examined, with an exception of the gills. Most likely, this is mediated by the presence of S. molnari-specific antibodies, first detectable in the serum of infected fish at 28 dpi. The observed reduction in parasite numbers can likely be attributed to opsonization of the parasite with the specific antibodies, activation of complement and cell-mediated killing by leukocytes . Interestingly, these effector functions did not clear the parasites in the gill, which might be explained by the partially intracellular nature of S. molnari spore-forming stages in the gills, with host epithelial cells serving as a protective niche, where the parasite remains unharmed [13, 62, 63].
Despite the presence of specific antibodies in the serum, parasite numbers recovered during the following days of infection and S. molnari invaded the liver, where the largest parasite quantities were observed. Liver stages are larger, and we demonstrated that blood stages with more than nine cells only occur once the liver becomes infected and likely due to release of these large stages from this host reservoir. After re-entry into the blood, it is unclear if further parasite proliferation occurs or if LS just use the blood stream for transport to the gills. However, at this stage, a second peak of BS occurrence was observed at 42 dpi, which is lower than the previous one but mirrored in all organs due to blood circulation within them, the same as at 28 dpi. Complex migration patterns with portals of entry different from the target tissue and several distinct developmental stages in the fish host are common to the Myxozoa, but proliferation in the blood is a characteristic specific only to sphaerosporids [64,65,66,67]. In S. dykovae, further development of BS occurs in the swim bladder, with secondary cells containing two tertiary cells, a stage identical with the earliest stages of sporogony in the kidney tubules [55, 64], and these triplets were also detected in Sphaerospora truttae . In S. molnari, the typical cell triplets were not observed in BS or LS, and it is assumed that they will only appear later, during sporogonic development in the gills. This suggests that the LS is a larger but not different stage of parasite development. It is, however, likely that nutrients are being taken up from the liver, as the primary cell membrane is almost indistinguishable (Fig. 1), indicating close contact with host cells and potentially active feeding mechanisms. Myxozoan plasmodia have been shown to exhibit a distinct periphery with high pinocytotic activity [69,70,71]. The survival and growth of S. molnari in the liver implies an effective use of exploitation and disguise mechanisms, similar to Plasmodium multiplication in the liver which is clinically silent because the parasite synthesizes host bulk phospholipids  and hijacks specific host pathways to promote its own development , while a number of mechanisms prevent host immune response and phagocytosis .
Between 42 and 49 dpi, the number of parasites decreased in all studied organs and remained low for the rest of the experiment, marking the onset of a chronic infection. However, despite the reduced number of parasites, we still observed increased numbers of thrombocytes and monocytes, possibly involved in parasite clearance by phagocytosis . Alternatively, the increased number of thrombocytes might be related to the resolution of inflammation and wound healing, as activated platelets can be potent drivers of anti-inflammatory responses .
The changes in leukocyte composition were complemented by changes in the expression of selected cytokines with known function in the acute phase of the immune response. Peritoneal injection of the parasite induced an immediate increase in the expression level of some cytokines, suggesting the recognition of S. molnari by pattern recognition receptors, such as toll-like receptors (TLRs). Although the interaction of TLRs with myxozoans has not been studied, their importance in the recognition of fish parasites has previously been demonstrated [27, 77]. Notably, the expression of il-1β and ifnγ increased only to a limited extent early after the injection of S. molnari and thereafter was inhibited (il-1β) or close to the levels observed in control fish (ifnγ). Similarly, inhibition of il-1β has been previously reported in the kidney of brown trout infected with Tetracapsuloides bryosalmonae, the agent of proliferative kidney disease (PKD) , in the intestine infected with E. leei  or C. shasta, the agent of enteronecrosis or gut-rot disease in salmonids . Low expression of ifnγ was previously reported in the kidney during the early stages of PKD (slightly increased by 50 dpi; ), while other myxozoan models suggest its increase and possible contribution to the cytotoxic activities [81, 82]. The timing of the downregulation of these genes is in accordance with the increasing expression of il-10 . Additionally, the transition from low pro-inflammatory to strong anti-inflammatory reaction has been described not only in myxozoans [29, 32, 84] but also in other parasitic infections . These findings suggest that despite initial activation of the pro-inflammatory response, the extent of the inflammation is limited shortly after.
Meanwhile, il-6 and tnfα were upregulated throughout the study indicating a prolonged activation of innate immunity. The increase of tnfα has been reported in other models including carp infected with Trypanosoma borreli  and in all examined organs of fish infected by E. leei . Notably, tnfα does not only play the classical role by stimulating the production of a number of genes associated with inflammation, enhancing the phagocytic activity of leukocytes, regulating homing, proliferation and migration , but it also exhibits lectin-like activity towards parasites such as Trypanosoma brucei . However, the effect of tnfα is not always beneficial, as the excessive production of tnfα was shown to be associated with the development of lesions and loss of barrier function during enteromyxosis in turbot . il-6 is known for its pleitropic effect playing a pivotal role during the transition from innate to acquired immunity. In fish, it has been reported to promote macrophage growth and production of antimicrobial proteins , while at the same time supporting antibody production [88, 89]. Possibly, in the S. molnari infection it might be related to the attraction of monocytes, secretion of IgM from 21 dpi onward and the resolution of the inflammation in the final stages. Involvement of il-6 has also been described in E. leei infected gilthead sea bream  and C. shasta infected Chinook salmon .
From the selected cytokines, the anti-inflammatory il-10 underwent the strongest upregulation, reaching up to 1456-fold increase in expression compared to control fish, at 56 dpi. As carp has two paralogues of il-10 which have almost identical gene structure, synteny, protein sequence and exert identical biological activities in vitro, it is worth noting that the primers used in the present study were targeting the il-10b paralogue, which is preferentially upregulated upon infection . Traditionally, il-10 is not only associated with the control of excessive immune responses but in cell cultures obtained from T. borreli-infected carp, il-10 promotes the proliferation, differentiation and antibody secretion by the IgM+ B cells . Importantly, a growing evidence from other disease models suggest that the protective function of il-10 can be exploited by pathogens . In infection with other Myxozoa, high levels of il-10 were induced by C. shasta , T. bryosalmonae  and also in the intestine, but not in head kidney of gilthead sea bream infected E. leei . These results point to il-10 functioning as a surrogate of myxozoan infections and question its role in host-pathogen interaction. Possibly, the expression of il-10 can be explained by two scenarios. In the traditional understanding, il-10 is induced by the immune system itself in an effort to minimize the collateral damage caused by neutrophils and monocytes fighting the growing number of S. molnari. Alternatively, S. molnari uses an unknown mechanism to induce il-10 as an immunomodulation strategy to deactivate the effector capacities of the immune system. The present data does not exclude any of these scenarios and future studies will investigate whether high concentrations of il-10 are a cause or a consequence of high pathogen burdens.
Using this model of S. molnari infection, we investigated the kinetics of antibody responses in infected fish, including their expression levels, level of specific antibody and number of B lymphocytes, over time. As the present model bypasses the natural route of infection via mucosal surfaces, where the activation of mucosal immunoglobulins could be expected, we focused solely on the response of systemic IgM [46, 92, 93]. Previously, myxozoan infections were reported to show a delayed development of antibody responses, with specific immunoglobulins detectable after 6–8 weeks . However, we observed an early onset of the reaction to S. molnari, with the first significant increase in IgM levels detectable at 14 dpi. Such timing is in agreement with the known time span for antibody responses in common carp towards model antigens such as LPS or KLH, which at 21 °C takes 10–28 days . Expression data were further supported by a 9-fold increase in the absolute numbers of IgM+ B cells in the blood. Such a prominent reaction of B lymphocytes with excessive expression of IgM points to a possible hypergammaglobulinemia, resulting in a polyclonal expansion and production of immunoglobulins with only a minor proportion of specific antibodies, a phenomena that has been reported in two other myxozoan infections [46, 95,96,97], in other parasites of fish (T. borreli ) and mammals (Trypanosoma cruzi ; Plasmodium chabaudi ).
To elucidate the nature of the B cell responses, we investigated the presence of specific antibodies to S. molnari by SDS PAGE of parasite antigens and subsequent western blotting with sera of infected fish, which identified a single protein of 16 kDa. While the identity of this protein remains to be determined, it is noteworthy that an antigenic protein of the same size (and an additional 165 kDa protein) was characterized in E. leei, using a similar approach . The two proteins could well correspond to common antigens shared among myxozoans. Estensoro et al.  showed the E. leei antigen to be a glycoprotein, and suggested a minicollagen; however, minicollagens would unlikely be expressed in proliferative stages of S. molnari as they are essential components of polar capsules and cnidarian stinging cells [100, 101], expressed during spore formation. Nevertheless, while the western blot clearly demonstrated the acquisition of specific immunity of carp against S. molnari, the low titer of specific antibodies in combination with an extremely high numbers of circulating B cells (up to 9-fold increase in number) indirectly confirms the suspicion that S. molnari, like other myxozoan parasites, manipulates B cells responses.
Our study demonstrates an activation of innate and adaptive humoral immune responses of common carp to the myxozoan S. molnari, during early stage of infection and parasite proliferation. We show that myxozoans use important immunomodulatory and immune evasion strategies such as intracellular disguise, motility, polyclonal activation of B cell responses and skewing of host response to an anti-inflammatory phenotype, in order to be able to successfully proliferate in their fish hosts and produce infective spore stages to continue their life-cycle. We believe that the present comprehensive characterization of host and parasite interactions between S. molnari and common carp represents a solid base for further research investigating newly raised questions in our in vivo model, with a focus on specific immune responses and universal myxozoan antigenic proteins that could in the future be used for targeted therapies and potential vaccine design.
Experimental infection of fish and sampling procedure
In May 2015, fertilized eggs of common carp (Cyprinus carpio) were obtained from the Faculty of Fisheries and Protection of Waters of the University of South Bohemia. They originated from a carp population reared in production ponds in Vodňany (Czech Republic). Eggs were stocked and fish were reared in an experimental recirculation system (UV filtration, ozone) at the animal facility of the Institute of Parasitology, Biology Centre of the Czech Academy of Sciences. Water temperature was kept at 21 ± 1 °C and fish were fed a commercial carp diet (Skretting, Stavanger, Norway) at a daily rate of 2% of biomass. In April 2016, 57 fish (weight 15 ± 4 g), were intraperitoneally injected with motile blood stages of S. molnari, isolated from whole blood of an infected carp population aged 0+. Isolation of the parasites was performed by centrifugation of whole blood from 18 donor fish (0+) in heparinized hematocrit capillaries (4000×g, 5 min). A mixture of white blood cells and parasite stages was collected after breaking of the capillary tubes. The cells were washed in PBS and after centrifugation (800×g, 5 min) they were resuspended, quantified and adjusted to 370 parasites/µl. A volume equaling 14,800 parasites per fish (averaging 1000 parasites per g of body weight, 40 µl of parasite/host cell solution) was intraperitoneally injected into SPF receptor fish. Thirty-six control fish were injected with 40 µl of PBS. This experiment was repeated with 10 specific-pathogen-free carp, to obtain information on B cell numbers by cytometry on a weekly basis. All fish were anesthetized with clove oil before bleeding or injection. Experimental and control carp were maintained in 3 isolated aquaria of 80 liters (2× infected, 1× control), with aeration, filtration, and weekly water changes. Sampling of fish was performed at 1, 2, 4 and 7 days post-infection (dpi) and thereafter weekly until week 9 post-infection, with 5 experimental and 3 control fish sampled per date. No mortalities occurred during the experiment. On each sampling date, blood was drawn with a heparinized syringe and full hematological analyses as well as flow cytometry were performed. Two blood smears were produced and fixed with methanol. Samples of blood, gills, kidney, liver and muscle were stored in TNES urea  for DNA extraction; the head kidney was fixed in RNAlater (Ambion, Europe Ltd, Huntingdon, UK) and frozen. Sera were collected and frozen after spinning the remainder of the blood sample at 4000×g for 10 min.
Full blood was diluted in Hayem’s solution (1:200) and blood cells were counted in a Bürker chamber, with the methodology detailed in Additional file 1: Figure S3. Two replicate blood smears per fish, fixed with methanol, were stained with Kwik–Diff (Richard Allen Scientific, San Diego, USA), according to the manufacturer’s recommendations and permanently mounted with Eukitt (Merck, Darmstadt, Germany). Two times 100 differentially stained cells (leukocytes, thrombocytes, S. molnari blood stages but excluding erythrocytes) were identified morphologically, based on known host and parasite features [13, 17, 103]. The compacted erythrocyte cell volume relative to whole blood volume (hematocrit) was measured in 2 capillaries per fish, drawn from whole blood and centrifuged at 10,000×g for 5 min. Hemoglobin content was estimated by spectrophotometry (540 nm) after lysis of cells and conversion of hemoglobin to cyanmethemoglobin, using Drabkin’s solution (Sigma; 1:250 dilution of blood:Drabkin’s; duplicate measurements). Hemoglobin concentrations were calculated from absorbance levels according to y (g l−1) = (A540–ABLANK) WHb FD (CE d 1000)−1, where A540 is the absorbance at 540 nm, ABLANK is the absorbance of pure Drabkin’s reagent at 540 nm, WHb is the molecular mass (human hemoglobin tetramer; = 64,458), FD is the dilution factor (= 250), CE is the millimolar extinction coefficient of tetrameric cyanomethhemoglobin at 540 nm (= 44), d is the vial light path in cm (= 1) and 1000 converts from mg to g .
To gain insights into the number and proportion of B lymphocytes from 7 dpi onwards and throughout the infection we adapted a protocol for the flow cytometric analysis of the full blood described previously . Briefly, 2 µl of blood of control and infected fish were washed with cold RPMI and stained for 20 min with a monoclonal antibody recognizing the heavy chain of carp IgM (1µg/ml) (Aquatic Diagnostics Ltd, Stirling, UK), followed by staining with goat-anti-mouse IgG Alexa Fluor 488 (2 µg/ml; Thermo Fisher Scientific, Pardubice, Czech Republic). The samples were washed twice and resuspended in 200 µl of RPMI. The proportion and total number of IgM+ B cells were determined using BD FACSCanto II (BD Biosciences, Prague, Czech Republic). Each sample was acquired for 30 s with a flow rate of 60 µl/min.
Parasite localization and quantification
Sphaerospora molnari was quantified in blood, gills, kidney, liver and muscle tissues, using real-time quantitative PCR which was performed using a LightCycler® 480 (Roche, Prague, Czech Republic). DNA was extracted using a modified phenol/chloroform extraction protocol , including RNAse A treatment at 50 µg/ml, for 30 min at 37 °C, before ethanol precipitation. DNA concentration and purity was estimated using a NanoDrop spectrophotometer (Thermo Fisher Scientific) and dilutions of 100 ng/µl were prepared using RNAse/DNAse-free water. Sphaerospora molnari stages were quantified using a TaqMan-based qPCR assay base on the SSU rDNA sequence of the parasite (GenBank accession number JX431511). Briefly, each 25 µl reaction contained 100 ng of DNA, 0.4 µM forward and reverse primers, 0.2 µM 5′ FAM and 3’BHQ1-labelled Taqman probe and 12.5 µl of LightCycler® 480 Probes Master (Roche) as detection chemistry. Samples were run as duplicates and Ct values for each sample were recalculated based on the average reaction efficiencies for each primer set, which were obtained upon comparative quantitation analysis. Ct values were compared after adjustment using an inter-plate control sample. Parasite SSU rDNA copy numbers were calculated relative to the sample with the highest parasite concentration (set to 100) and normalized relative to host β-actin (housekeeping gene). All qPCR primers and probes are detailed in Additional file 1: Table S1. Multicellular blood stages were furthermore quantified on Kwik–Diff-stained blood smears (see above) to be able to relate qPCR values to parasite numbers of varying cell composition in the blood.
In order to elucidate the early parasite migration/uptake from the intraperitoneal injection site, in a second trial we labelled isolated S. molnari blood stages (BS) with carboxyfluorescein succinimidyl ester (CellTrace™ CFSE Cell Proliferation Kit; Thermo Fisher Scientific) at 1:10,000 concentration, washed the cells twice with RPMI medium and injected 100,000 BS per fish into the peritoneal cavity of 15 SPF carp. Peritoneal lavages were performed at 16 h, 40 h, 7 days and 21 days post-injection. Percoll-isolated leukocyte + parasite preparations (57%/34% Percoll interphase collected after centrifugation at 400×g, 20 min) of head kidney, spleen and blood were also prepared and the recovered cells were analyzed by flow cytometry (number of CSFE-labelled S. molnari stages per 100,000 cells). The cells isolated from blood were microscopically screened for motile BS at 21 dpi. Three fish were sampled at each time point.
IgM expression in head kidney
Transcriptomic data obtained from a mixture of white blood cells and S. molnari blood stages was mined for IgM heavy chain sequences of C. carpio. Two isoforms were found, one with a transmembrane domain (on B cell/plasma cell surface; GenBank: MH352353) and one with a secretory tail (secreted form; GenBank: MH352354). We developed TaqMan qPCR assays to differentiate between the expression of the secretory tetramer (IgMsec) and the membrane-bound monomer (IgMmem) (Additional file 1: Table S1). IgM expression was determined in the head kidney of experimental and control fish, after extraction of total RNA from the RNAlater-fixed samples using the NucleoSpin RNA (Macherey-Nagel, Düren, Germany) and cDNA preparation with the Transcriptor High Fidelity cDNA Synthesis Kit (Roche Diagnostics). Both RNA as well as cDNA quality and quantity was determined by NanoDrop measurement. qPCRs were performed on samples adjusted to 100 ng/µl cDNA concentrations and relative to carp β-actin (Additional file 1: Table S1), using the same quality and compatibility measures as above.
Detection of S. molnari specific antibodies by western blot
In order to determine specific antibody responses induced by the infection with S. molnari, we performed western blotting with sera of infected carp. For antigen preparation, blood stages of S. molnari were obtained from infected fish at the peak of infection (28 dpi, see “Results”). Whole blood was diluted in 5 ml of cold RPMI and layered onto Ficoll. Following gradient centrifugation at 800×g for 20 min, a mixture of white blood cells and parasites was collected on the top of the Ficoll layer and washed once with cold RPMI. Parasites were enriched by allowing adhesion of white blood cells to a cell culture dish for 30 min, before being pelleted by centrifugation and resuspended in cold RPMI. The purity of the resulting parasite mixture was assessed under the microscope. The parasite preparation was mixed with 2× Laemmli buffer (Bio-Rad, Prague, Czech Republic) and used as an antigen for the western blotting following the manufacturer’s protocol. Proteins divided by SDS PAGE in a 12% gel were transferred on PVDF membrane using the Trans-Blot Turbo™ Transfer System. Following blocking with 8% of dry milk dissolved in PBS-Tween, the membrane was cut into strips, which were exposed to the pooled sera of infected fish (5 fish per time point) diluted 1:100 in PBS-Tween. After washing with PBS-Tween, mouse anti-carp IgM antibody (see above) was used as a primary antibody to visualize specific carp IgM bound to parasite protein, followed by the incubation with secondary HRP-labeled goat anti-mouse IgG (Thermo Fisher Scientific). Finally, after the addition of ELC substrate, the western blot was visualized on a ChemiDoc Imaging System (Bio-Rad).
Cytokine gene expression
To obtain further insights into the regulation of immune responses the expression of a panel of pro- and anti-inflammatory cytokines was analyzed by qPCR (Additional file 1: Table S1). The head kidney was selected for the analysis as it is considered the principal immune organ in fish combining primary and secondary immune function . Furthermore, the slow blood flow through its stroma presumably allows a direct contact with the parasite once it appears in the blood. The reaction mix contained 0.26 µM forward and reverse primers, 100 ng of cDNA and 7 µl of FastStart Universal SYBR Green Master Mix (Roche) in 16 µl reactions. All samples were run in duplicates, applying the same quality standards as above. The relative expression ratio of each sample was calculated according to Pfaffl , based on the take-off deviation of sample versus controls at each time point and normalized relative to host β-actin (housekeeping gene).
Cytokine and IgM expression differences were evaluated with linear models (LM), using log-transformed response variables. Thereafter, multiple pairwise comparisons evaluating the influence of fixed effect (dpi) were obtained using Tukey’s all-pair comparisons, applying the Bonferroni correction to adjust the P-values (package multcomp, v.1.3-3; ). To investigate differences in cells numbers (lymphocytes, B cells, monocytes, neutrophils and thrombocytes) between days and groups, the data were analyzed using linear models with log-transformed response variables and an interaction of the fixed factors dpi and group (control/experimental). Thereafter, pairwise comparisons evaluating the influence of the fixed effects on cell numbers were obtained as explained above. Parasite numbers were compared between sampling dates using the Kruskal Wallis H test. All analyses were conducted in R (R Core Team, v.3.4.2.). Graphs were prepared in R and colored/assembled in Adobe Photoshop.
blood stage (of S. molnari)
- Ct value:
cycle threshold value
days post infection
hours post infection
linear mixed model
liver stage (of S. molnari)
million years ago
- SDS PAGE:
sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Holzer AS, Bartošová-Sojková P, Born-Torrijos A, Lövy A, Hartigan A, Fiala I. The joint evolution of the Myxozoa and their alternate hosts: a cnidarian recipe for success and vast biodiversity. Mol Biol. 2018;27:1651–66.
Agrawal A. Transposition and evolution of antigen-specific immunity. Science. 2000;290:1715–6.
Bayne CJ. Origins and evolutionary relationships between the innate and adaptive arms of immune systems. Integr Comp Biol. 2003;43:293–9.
Tort L, Balasch JC, Mackenzie S. Fish immune system. A crossroads between innate and adaptive responses. Immunologia. 2003;22:277–86.
Gibson AK, Stoy KS, Gelarden IA, Penley MJ, Lively CM, Morran LT. The evolution of reduced antagonism - a role for host–parasite coevolution. Evolution. 2015;69:2820–30.
Atkinson SD, Bartholomew JL, Lotan T. Myxozoans: ancient metazoan parasites find a home in phylum Cnidaria. Zoology. 2018;129:66–8.
Carraro L, Bertuzzo E, Mari L, Fontes I, Hartikainen H, Strepparava N, et al. Integrated field, laboratory, and theoretical study of PKD spread in a Swiss prealpine river. Proc Natl Acad Sci USA. 2017;114:11992–7.
Okamura B, Hartikainen H, Schmidt-Posthaus H, Wahli T. Life cycle complexity, environmental change and the emerging status of salmonid proliferative kidney disease. Freshw Biol. 2011;56:735–53.
Sitjà-Bobadilla A, Schmidt-Posthaus H, Wahli T, Holland JW, Secombes CJ. Fish immune responses to Myxozoa. In: Okamura B, Gruhl A, Bartholomew JL, editors. Myxozoan evolution, ecology and development. Cham: Springer International Publishing; 2015. p. 253–80.
Eszterbauer E, Atkinson S, Diamant A, Morris D, El-Matbouli Mansour M, Hartikainen H. Myxozoan life cycles: practical approaches and insights. In: Okamura B, Gruhl A, Bartholomew JL, editors. Myxozoan evolution, ecology and development. Cham: Springer International Publishing; 2015. p. 175–98.
Al-Samman A, Molnár K, Székely C, Reiczigel J. Reno-, hepato- and splenomegaly of common carp fingerlings (Cyprinus carpio L.) diseased in swimbladder inflammation caused by Sphaerospora renicola Dyková et Lom, 1982. Acta Vet Hung. 2003;51:321–9.
Dyková I, Lom J. Review of pathogenic myxosporeans in intensive culture of carp (Cyprinus carpio) in Europe. Folia Parasitol. 1988;35:289–307.
Lom J, Dyková I, Pavlásková M, Grupcheva G. Sphaerospora molnari sp.nov. (Myxozoa:Myxosporea), an agent of gill, skin and blood sphaerosporosis of common carp in Europe. Parasitology. 1983;86:529.
Bartošová P, Fiala I, Jirků M, Cinková M, Caffara M, Fioravanti ML, et al. Sphaerospora sensu stricto: taxonomy, diversity and evolution of a unique lineage of myxosporeans (Myxozoa). Mol Phylogenet Evol. 2013;68:93–105.
Patra S, Bartošová-Sojková P, Pecková H, Fiala I, Eszterbauer E, Holzer AS. Biodiversity and host–parasite cophylogeny of Sphaerospora (sensu stricto) (Cnidaria: Myxozoa). Parasit Vectors. 2018;11:347.
Csaba GY. An unidentifiable extracellular sporozoa parasite from the blood of the carp. Parasitol Hungarica. 1976;9:21–4.
Hartigan A, Estensoro I, Vancová M, Bílý T, Patra S, Eszterbauer E, et al. New cell motility model observed in parasitic cnidarian Sphaerospora molnari (Myxozoa:Myxosporea) blood stages in fish. Sci Rep. 2016;6:39093.
Ibarra AM, Gall GAE, Hedrick RP. Susceptibility of two strains of rainbow trout Oncorhynchus mykiss to experimentally induced infections with the myxosporean Ceratomyxa shasta. Dis Aquat Organ. 1991;10:191–4.
Ibarra AM, Gall GAE, Hedrick RP. Experimental infection of rainbow trout, Oncorhynchus mykiss (Walbaum), with cryopreserved developmental and sporogonic stages of the myxosporean Ceratomyxa shasta (Noble). J Fish Dis. 1992;15:353–5.
Palikova M, Papezikova I, Markova Z, Navratil S, Mares J, Mares L, et al. Proliferative kidney disease in rainbow trout (Oncorhynchus mykiss) under intensive breeding conditions: Pathogenesis and hematological and immune parameters. Vet Parasitol. 2017;238:5–16.
Gómez D, Bartholomew J, Sunyer JO. Biology and mucosal immunity to myxozoans. Dev Comp Immunol. 2014;43:243–56.
El-Matbouli M, Hoffmann RW, Mandok C. Light and electron microscopic observations on the route of the triactinomyxon-sporoplasm of Myxobolus cerebralis from epidermis into rainbow trout cartilage. J Fish Biol. 1995;46:919–35.
Ferguson JA, Atkinson SD, Whipps CM, Kent ML. Molecular and morphological analysis of Myxobolus spp. of salmonid fishes with the description of a new Myxobolus species. J Parasitol. 2008;94:1322–34.
Amandi A, Holt RA, Fryer L. Observations on Myxobolus insidiosus (Myxozoa: Myxosporea) a parasite of salmonid fishes. Fish Pathol. 1985;20:287–304.
El-Matbouli M, Hoffmann R, Schoel H, McDowell T, Hedrick R. Whirling disease: host specificity and interaction between the actinosporean stage of Myxobolus cerebralis and rainbow trout Oncorhynchus mykiss. Dis Aquat Organ. 1999;35:1–12.
Bartholomew JL, Smith CE, Rohovec JS, Fryer JL. Characterization of a host response to the myxosporean parasite, Ceratomyxa shasta (Noble), by histology, scanning electron microscopy and immunological techniques. J Fish Dis. 1989;12:509–22.
Alvarez-Pellitero P. Fish immunity and parasite infections: from innate immunity to immunoprophylactic prospects. Vet Immunol Immunopathol. 2008;126:171–98.
Lovy J, Goodwin AE, Speare DJ, Wadowska DW, Wright GM. Histochemical and ultrastructural analysis of pathology and cell responses in gills of channel catfish affected with proliferative gill disease. Dis Aquat Organ. 2011;94:125–34.
Gorgoglione B, Wang T, Secombes CJ, Holland JW. Immune gene expression profiling of proliferative kidney disease in rainbow trout Oncorhynchus mykiss reveals a dominance of anti-inflammatory, antibody and T helper cell-like activities. Vet Res. 2013;44:55.
Sitjà-Bobadilla A. Fish immune responses to myxozoan parasites. Parasite. 2008;15:420–5.
Davey GC, Calduch-Giner JA, Houeix B, Talbot A, Sitjà-Bobadilla A, Prunet P, et al. Molecular profiling of the gilthead sea bream (Sparus aurata L.) response to chronic exposure to the myxosporean parasite Enteromyxum leei. Mol Immunol. 2011;48:2102–12.
Bjork SJ, Zhang Y-A, Hurst CN, Alonso-Naveiro ME, Alexander JD, Sunyer JO, et al. Defenses of susceptible and resistant Chinook salmon (Onchorhynchus tshawytscha) against the myxozoan parasite Ceratomyxa shasta. Fish Shellfish Immunol. 2014;37:87–95.
Braden LM, Rasmussen KJ, Purcell SL, Ellis L, Mahony A, Cho S, et al. Acquired protective immunity in Atlantic salmon Salmo salar against the myxozoan Kudoa thyrsites involves induction of MHIIβ+ CD83+ antigen-presenting cells. Infect Immun. 2018;86:e00556.
Foott JS, Hedrick RP. Seasonal occurrence of the infectious stage of proliferative kidney disease (PKD) and resistance of rainbow trout, Salmo gairdneri Richardson, to reinfection. J Fish Biol. 1987;30:477–83.
Golomazou E, Athanassopoulou F, Karagouni E, Tsagozis P, Tsantilas H, Vagianou S. Experimental transmission of Enteromyxum leei Diamant, Lom and Dykova, 1994 in sharpsnout sea bream, Diplodus puntazzo C. and the effect on some innate immune parameters. Aquaculture. 2006;260:44–53.
Körting W, Kruse P, Steinhagen D. Development of “Csaba cells” in experimentally infected Cyprinus carpio. Angew Parasitol. 1989;30:185–8.
Feist SW, Morris DJ, Alama-Bermejo G, Holzer AS. Cellular process in myxozoans. In: Okamura B, Gruhl A, Bartholomew JL, editors. Myxozoan evolution, ecology and development. Cham: Springer International Publishing; 2015. p. 139–54.
Alama-Bermejo G, Bron JE, Raga JA, Holzer AS. 3D morphology, ultrastructure and development of Ceratomyxa puntazzi stages: first insights into the mechanisms of motility and budding in the Myxozoa. PLoS One. 2012;7:e32679.
Roederer M. Interpretation of cellular proliferation data: Avoid the panglossian. Cytom Part A. 2011;79:95–101.
Kudriashova IV, Naumova AM. Effect of Haemogregarina on the body of carp. Veterinariia. 1978;4:76–8.
Molnár K. Renal sphaerosporosis in the common carp Cypnus carpio L. J Fish Dis. 1980;3:11–9.
Holzer AS, Hartigan A, Patra S, Pecková H, Eszterbauer E. Molecular fingerprinting of the myxozoan community in common carp suffering swim bladder inflammation (SBI) identifies multiple etiological agents. Parasit Vectors. 2014;7:398.
Amezcua Vesely MCMC, Bermejo DA, Montes CL, Acosta-Rodríguez EV, Gruppi A. B-cell response during protozoan parasite infections. J Parasitol Res. 2012;2012:362131.
Bermejo DA, Amezcua Vesely MC, Khan M, Acosta Rodríguez EV, Montes CL, Merino MC, et al. Trypanosoma cruzi infection induces a massive extrafollicular and follicular splenic B-cell response which is a high source of non-parasite-specific antibodies. Immunology. 2011;132:123–33.
Achtman AH, Khan M, MacLennan ICM, Langhorne J. Plasmodium chabaudi chabaudi infection in mice induces strong B cell responses and striking but temporary changes in splenic cell distribution. J Immunol. 2003;171:317–24.
Abos B, Estensoro I, Perdiguero P, Faber M, Hu Y, Díaz Rosales P, et al. Dysregulation of B cell activity during proliferative kidney disease in rainbow trout. Front Immunol. 2018;9:1203.
Joerink M, Groeneveld A, Ducro B, Savelkoul HFJ, Wiegertjes GF. Mixed infection with Trypanoplasma borreli and Trypanosoma carassii induces protection: Involvement of cross-reactive antibodies. Dev Comp Immunol. 2007;31:903–15.
Chistiakov DA, Myasoedova VA, Revin VV, Orekhov AN, Bobryshev YV. The impact of interferon-regulatory factors to macrophage differentiation and polarization into M1 and M2. Immunobiology. 2018;223:101–11.
Du X, Williams DA. Interleukin-11: review of molecular, cell biology, and clinical use. Blood. 1997;89:3897–908.
Zou J, Secombes CJ. The function of fish cytokines. Biology. 2016;5:23.
Holzer AS, Sommerville C, Wootten R. Tracing the route of Sphaerospora truttae from the entry locus to the target organ of the host, Salmo salar L., using an optimized and specific in situ hybridization technique. J Fish Dis. 2003;26:647–55.
Bjork SJ, Bartholomew JL. Invasion of Ceratomyxa shasta (Myxozoa) and comparison of migration to the intestine between susceptible and resistant fish hosts. Int J Parasitol. 2010;40:1087–95.
Nagy JA. Lymphatic and nonlymphatic pathways of peritoneal absorption in mice: physiology versus pathology. Blood Purif. 1992;10:148–62.
Dundas K, Shears MJ, Sun Y, Hopp CS, Crosnier C, Metcalf T, et al. Alpha-v–containing integrins are host receptors for the Plasmodium falciparum sporozoite surface protein, TRAP. Proc Natl Acad Sci USA. 2018;115:4477–82.
Morris DJ. Cell formation by myxozoan species is not explained by dogma. Proc Biol Sci. 2010;277:2565–70.
Abdallah DSA, Denkers EY. Neutrophils cast extracellular traps in response to protozoan parasites. Front Immunol. 2012;3:382.
Denkers EY, Schneider AG, Cohen SB, Butcher BA. Phagocyte responses to protozoan infection and how Toxoplasma gondii meets the challenge. PLoS Pathog. 2012;8:e1002794.
Sheel M, Engwerda CR. The diverse roles of monocytes in inflammation caused by protozoan parasitic diseases. Trends Parasitol. 2012;28:408–16.
Estensoro I, Mulero I, Redondo MJ, Alvarez-Pellitero P, Mulero V, Sitjà-Bobadilla A. Modulation of leukocytic populations of gilthead sea bream (Sparus aurata) by the intestinal parasite Enteromyxum leei (Myxozoa: Myxosporea). Parasitology. 2014;141:425–40.
Cuesta A, Salinas I, Rodriguez A, Munoz P, Sitjà-Bobadilla A, Alvarez-Pellitero P, et al. Cell-mediated cytotoxicity is the main innate immune mechanism involved in the cellular defence of gilthead seabream (Teleostei: Sparidae) against Enteromyxum leei (Myxozoa). Parasite Immunol. 2006;28:657–65.
Forthal DN. Functions of antibodies. Microbiol Spectr. 2014;2:1–17.
Molnár K. Myxobolus pavlovskii (Achmerov, 1954) (Myxosporidia) infection in the silver carp and bighead. Acta Vet Acad Sci Hung. 1979;27:207–16.
Kaup FJ, Kuhn EM, Körting W. Light and electron microscopic studies on the sporogenesis of Sphaerospora molnari in the gill lamellas of carp (Cyprinus carpio). Berl Munch Tierarztl Wochenschr. 1995;108:206–14.
Baska F, Molnár K. Blood stages of Sphaerospora spp (Myxosporea) in cyprinid fishes. Dis Aquat Organ. 1988;5:23–8.
Dyková I, Lom J, Körting W. Light and electron microscopic observations on the swimbladder stages of Sphaerospora renicola, a parasite of carp (Cyprinus carpio). Parasitol Res. 1990;76:228–37.
Supamattaya K, Fischer-Scherl T, Hoffmann RW. Boonyaratpalin S. Sphaerospora epinepheli n. sp. (Myxosporea: Sphaerosporidae) observed in grouper (Epinephelus malabaricus). J Protozool. 1991;38:448–54.
Lom J, Pavlásková M, Dyková I. Notes on kidney-infecting species of the genus Sphaerospora Thelohan (Myxosporea), including a new species S. gobionis sp. nov., and on myxosporean life cycle stages in the blood of some freshwater fish. J Fish Dis. 1985;8:221–32.
McGeorge J, Sommerville C, Wootten R. Light and electron microscope observations on extrasporogonic and sporogonic stages of a myxosporean parasite of the genus Sphaerospora Thelohan, 1892 from Atlantic salmon, Salmo solar L., in Scotland. J Fish Dis. 1994;17:227–38.
Azevedo C, Clemente SCDS, Casal G, Matos P, Oliveira E, Al-Quraishy S, et al. Ultrastructure of the plasmodial development of Myxobolus insignis (Myxozoa), infecting the Amazonian fish Semaprochilodus insignis (Prochilodontidae). Acta Protozool. 2013;52:91–7.
Current WL, Janovy JJ, Knight SA. Myxosoma funduli Kudo (Myxosporida) in Fundulus kansae: ultrastructure of the plasmodium wall and of sporogenesis. J Protozool. 1979;26:574–83.
Stehr C, Whitaker DJ. Host–parasite interaction of the myxosporeans Kudoa paniformis Kabata & Whitaker, 1981 and Kudoa thyrsites (Gilchrist, 1924) in the muscle of Pacific whiting, Merluccius productus (Ayres): an ultrastructural study. J Fish Dis. 1986;9:505–17.
Itoe MA, Sampaio JL, Cabal GG, Real E, Zuzarte-Luis V, March S, et al. Host cell phosphatidylcholine is a key mediator of malaria parasite survival during liver stage infection. Cell Host Microbe. 2014;16:778–86.
Posfai D, Sylvester K, Reddy A, Ganley JG, Wirth J, Cullen QE, et al. Plasmodium parasite exploits host aquaporin-3 during liver stage malaria infection. PLoS Pathog. 2018;14:e1007057.
Agop-Nersesian C, Niklaus L, Wacker R, Theo Heussler V. Host cell cytosolic immune response during Plasmodium liver stage development. FEMS Microbiol Rev. 2018;42:324–34.
Nagasawa T, Nakayasu C, Rieger AM, Barreda DR, Somamoto T, Nakao M. Phagocytosis by thrombocytes is a conserved innate immune mechanism in lower vertebrates. Front Immunol. 2014;5:445.
Linke B, Schreiber Y, Picard-Willems B, Slattery P, Nüsing RM, Harder S, et al. Activated platelets induce an anti-inflammatory response of monocytes/macrophages through cross-regulation of PGE2 and cytokines. Mediators Inflamm. 2017;2017:1463216.
Findly RC, Zhao X, Noe J, Camus AC, Dickerson HW. B cell memory following infection and challenge of channel catfish with Ichthyophthirius multifiliis. Dev Comp Immunol. 2013;39:302–11.
Bailey C, Strepparava N, Wahli T, Segner H. Exploring the immune response, tolerance and resistance in proliferative kidney disease of salmonids. Dev Comp Immunol. 2019;90:165–75.
Pérez-Cordón G, Estensoro I, Benedito-Palos L, Calduch-Giner JA, Sitjà-Bobadilla A, Pérez-Sánchez J. Interleukin gene expression is strongly modulated at the local level in a fish-parasite model. Fish Shellfish Immunol. 2014;37:201–8.
Bailey C, Schmidt-Posthaus H, Segner H, Wahli T, Strepparava N. Are brown trout Salmo trutta fario and rainbow trout Oncorhynchus mykiss two of a kind? A comparative study of salmonids to temperature-influenced Tetracapsuloides bryosalmonae infection. J Fish Dis. 2018;41:191–8.
Ronza P, Robledo D, Bermúdez R, Losada AP, Pardo BG, Sitjà-Bobadilla A, et al. RNA-seq analysis of early enteromyxosis in turbot (Scophthalmus maximus): new insights into parasite invasion and immune evasion strategies. Int J Parasitol. 2016;46:507–17.
Piazzon MC, Estensoro I, Calduch-Giner JA, del Pozo R, Picard-Sánchez A, Pérez-Sánchez J, Sitjà-Bobadilla A. Hints on T cell responses in a fish-parasite model: Enteromyxum leei induces differential expression of T cell signature molecules depending on the organ and the infection status. Parasit Vectors. 2018;11:443.
Piazzon MC, Wentzel AS, Wiegertjes GF, Forlenza M. Carp Il10a and Il10b exert identical biological activities in vitro, but are differentially regulated in vivo. Dev Comp Immunol. 2017;67:350–60.
Costa MM, Maehr T, Diaz-Rosales P, Secombes CJ, Wang T. Bioactivity studies of rainbow trout (Oncorhynchus mykiss) interleukin-6: effects on macrophage growth and antimicrobial peptide gene expression. Mol Immunol. 2011;48:1903–16.
Marinho FV, Alves CC, de Souza SC, da Silva CMG, Cassali GD, Oliveira SC, et al. Schistosoma mansoni tegument (Smteg) induces IL-10 and modulates experimental airway Inflammation. PLoS ONE. 2016;11:e0160118.
Forlenza M, Magez S, Scharsack JP, Westphal A, Savelkoul HFJ, Wiegertjes GF, et al. Receptor-mediated and lectin-like activities of carp (Cyprinus carpio) TNF-alpha. J Immunol. 2009;183:5319–32.
Ronza P, Bermúdez R, Losada AP, Sitjà-Bobadilla A, Pardo BG, Quiroga MI. Immunohistochemical detection and gene expression of TNFα in turbot (Scophthalmus maximus) enteromyxosis. Fish Shellfish Immunol. 2015;47:368–76.
Kaneda M, Odaka T, Suetake H, Tahara D, Miyadai T. Teleost IL-6 promotes antibody production through STAT3 signaling via IL-6R and gp130. Dev Comp Immunol. 2012;38:224–31.
Abós B, Wang T, Castro R, Granja AG, Leal E, Havixbeck J, et al. Distinct differentiation programs triggered by IL-6 and LPS in teleost IgM+ B tells in the absence of germinal centers. Sci Rep. 2016;6:30004.
Piazzon MC, Savelkoul HFJ, Pietretti D, Wiegertjes GF, Forlenza M. Carp Il10 has anti-inflammatory activities on phagocytes, promotes proliferation of memory T cells, and regulates B cell differentiation and antibody secretion. J Immunol. 2015;194:187–99.
Mege J-L, Meghari S, Honstettre A, Capo C, Raoult D. The two faces of interleukin 10 in human infectious diseases. Lancet Infect Dis. 2006;6:557–69.
Mashoof S, Criscitiello MF. Fish immunoglobulins. Biology (Basel). 2016;5:45.
Zhang Y, Salinas I, Li J, Parra D, Bjork S, Xu Z, et al. IgT, a primitive immunoglobulin class specialized in mucosal immunity. Nat Immunol. 2010;11:827–35.
Engelsma MY, Hougee S, Nap D, Hofenk M, Rombout JHWM, Van Muiswinkel WB, et al. Multiple acute temperature stress affects leucocyte populations and antibody responses in common carp, Cyprinus carpio L. Fish Shellfish Immunol. 2003;15:397–410.
Olesen N, Jørgensen P. Quantification of serum immunoglobulin in rainbow trout Salmo gairdneri under various environmental conditions. Dis Aquat Organ. 1985;1:183–9.
Piazzon MC, Galindo-Villegas J, Pereiro P, Estensoro I, Calduch-Giner JA, Gómez-Casado E, et al. Differential modulation of IgT and IgM upon parasitic, bacterial, viral, and dietary challenges in a perciform fish. Front Immunol. 2016;7:637.
Estensoro I, Calduch-Giner JA, Kaushik S, Pérez-Sánchez J, Sitjà-Bobadilla A. Modulation of the IgM gene expression and IgM immunoreactive cell distribution by the nutritional background in gilthead sea bream (Sparus aurata) challenged with Enteromyxum leei (Myxozoa). Fish Shellfish Immunol. 2012;33:401–10.
Saeij JPJ, De Vries BJ, Wiegertjes GF. The immune response of carp to Trypanoplasma borreli: Kinetics of immune gene expression and polyclonal lymphocyte activation. Dev Comp Immunol. 2003;27:859–74.
Estensoro I, Álvarez-Pellitero P, Sitjà-Bobadilla A. Antigenic characterization of Enteromyxum leei (Myxozoa: Myxosporea). Dis Aquat Organ. 2013;106:149–62.
Engel U, Ozbek S, Streitwolf-Engel R, Petri B, Lottspeich F, Holstein TW, et al. Nowa, a novel protein with minicollagen Cys-rich domains, is involved in nematocyst formation in Hydra. J Cell Sci. 2002;115:3923–34.
Adamczyk P, Zenkert C, Balasubramanian PG, Yamada S, Murakoshi S, Sugahara K, et al. A non-sulfated chondroitin stabilizes membrane tubulation in cnidarian organelles. J Biol Chem. 2010;285:25613–23.
Asahida T, Kobayashi T, Saitoh K, Nakayama I. Tissue preservation and total DNA extraction form fish stored at ambient temperature using buffers containing high concentration of urea. Fish Sci. 1996;62:727–30.
Ferguson HW, Bjerkas E, Evensen O. Systemic pathology of fish: a text and atlas of normal tissues in teleosts and their responses in disease. 2nd ed. London: Scotian Press; 2006.
Clark TD, Eliason EJ, Sandblom E, Hinch SG, Farrell AP. Calibration of a hand-held haemoglobin analyser for use on fish blood. J Fish Biol. 2008;73:2587–95.
Korytář T, Dang Thi H, Takizawa F, Köllner B. A multicolour flow cytometry identifying defined leukocyte subsets of rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 2013;35:2017–9.
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.
Zapata A, Diez B, Cejalvo T, Gutiérrez-De Frías C, Cortés A. Ontogeny of the immune system of fish. Fish Shellfish Immunol. 2006;20:126–36.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biometrical J. 2008;50:346–63.
We would like to thank David Gela for providing fertilized carp eggs as well as Pavel Lepič, Tomas Douda and Joana Pimentel for fish keeping and recirculation system maintenance.
This work was supported by the European Commission under the H2020 Programme (project reference 634429, ParaFishControl), the Czech Science Foundation (project reference 19-28399X for ASH, 19-28389Y for TK) and the Technology Agency of the Czech Republic (TG02010016, salary support for AB-T). Fish rearing at the partner organization Faculty of Fisheries and Protection of Waters was supported by the Ministry of Education, Youth and Sports of the Czech Republic (project CENAKVA no. CZ.1.05/2.1.00/01.0024) and the Ministry of Agriculture of the Czech Republic (project NAZV QK1820354) for CY-R and TP.
Availability of data and materials
All data supporting the findings of this study are available within the article and its additional file.
TK performed the flow cytometry analysis, western-blotting experiments and co-wrote the manuscript. GFW and ASW provided training in cell isolation methods, antibodies, qPCR primers and assistance with writing the manuscript. EZ trained and assisted with hematological analysis. AT carried out the experiments, collected all samples and performed hematological analyses and cell counts. RŠ assisted in qPCR assay design, while qPCR analyses were performed by ML (S. molnari), SP and SB-M (cytokines), and VS (IgM). AB-T did the statistical analysis of the data. CY-R and TP were responsible for fish rearing and helped with the aquaculture system setup. ASH planned and designed the study, supervised the findings and analyses of the data, and co-wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Animal procedures were performed in accordance with Czech legislation (section 29 of Act no. 246/1992 Coll., on Protection of animals against cruelty, as amended by Act no. 77/2004 Coll.) and animal handling complied with the relevant European guidelines on animal welfare (Directive 2010/63/EU on the protection of animals used for scientific purposes) and the recommendations of the Federation of Laboratory Animal Science Associations.
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.